title: An equilibrated a posteriori error estimator for arbitrary-order N\'ed\'elec elements for magnetostatic problems

abstract: We present a novel \textit{a posteriori} error estimator for N\'ed\'elec elements for magnetostatic problems that is constant-free, i.e. it provides an upper bound on the error that does not involve a generic constant. The estimator is based on equilibration of the magnetic field and only involves small local problems that can be solved in parallel. Such an error estimator is already available for the lowest-degree N\'ed\'elec element [D. Braess, J. Sch\"oberl, \textit{Equilibrated residual error estimator for edge elements}, Math. Comp. 77 (2008)] and requires solving local problems on vertex patches. The novelty of our estimator is that it can be applied to N\'ed\'elec elements of arbitrary degree. Furthermore, our estimator does not require solving problems on vertex patches, but instead requires solving problems on only single elements, single faces, and very small sets of nodes. We prove reliability and efficiency of the estimator and present several numerical examples that confirm this.


title: Towards a Unified Approach to Electromagnetic Analysis by Objects Embedded in Multilayers

abstract: In this paper, an efficient and accurate unified approach is proposed to solve transverse magnetic scattering problems by multilayer embedded objects. In the proposed approach, an equivalent current density is derived when the equivalent theorem is recursively applied on each boundary from inner to exterior regions. Then, the objects are replaced by the background medium and the equivalent electric current density enforced on the outermost boundary remains fields in the exterior region unchanged. Therefore, the scattering problems by the multilayer embedded objects can be solved with the electric field integral equation (EFIE). Compared with other two-region formulations (TRFs), like the PMCHWT, CTF, the proposed approach shows significant benefits: only the surface electric current density instead of both the electric and magnetic current densities is required to model the complex objects. Furthermore, the single electric current density is only enforced on the outermost boundary of objects. Therefore, the overall count of unknowns can be significantly reduced, especially when the number of boundaries between different homogeneous media is large. At last, several numerical experiments are performed to validate the accuracy and efficiency of the proposed approach.


title: Dynamic Boundary Guarding Against Radially Incoming Targets

abstract: We introduce a dynamic vehicle routing problem in which a single vehicle seeks to guard a circular perimeter against radially inward moving targets. Targets are generated uniformly as per a Poisson process in time with a fixed arrival rate on the boundary of a circle with a larger radius and concentric with the perimeter. Upon generation, each target moves radially inward toward the perimeter with a fixed speed. The aim of the vehicle is to maximize the capture fraction, i.e., the fraction of targets intercepted before they enter the perimeter. We first obtain a fundamental upper bound on the capture fraction which is independent of any policy followed by the vehicle. We analyze several policies in the low and high arrival rates of target generation. For low arrival, we propose and analyze a First-Come-First-Served and a Look-Ahead policy based on repeated computation of the path that passes through maximum number of unintercepted targets. For high arrival, we design and analyze a policy based on repeated computation of Euclidean Minimum Hamiltonian path through a fraction of existing targets and show that it is within a constant factor of the optimal. Finally, we provide a numerical study of the performance of the policies in parameter regimes beyond the scope of the analysis.


title: Antibacterial and antifungal activity of various extracts of Bacopa monnieri

abstract: Bacopa monnieri was also known as Brahmi. Theses plant has been considered as medicinal plants in Ayurvedic systems for centuries. This present study was carried out to evaluate the antibacterial and antifungal properties of ethanolic, diethyl ether, ethyl acetate, and aqueous extracts of B. monnieri. Agar disc diffusion tests were carried out to determine the antimicrobial effects of ethanolic, diethyl ether, ethyl acetate and aqueous extracts of B. monnieri against Gram-positive (Staphylococcus aureus ATCC 25923), Gram-negative (Escherichia coli ATCC 25922) bacterial strains and antifungal strains (Aspergillus flavus, and Candida albicans). Among the various extracts, diethyl ether extracts of B. monnieri has an antibacterial potency against Staphylococcus aureus (gram positive), and ethyl acetate extract showed effects on E. coli (gram negative) at higher concentrations of 300 microgam mL-1. The ethanolic extract has potent antifungal activity against the fungus (Aspergillus flavus, and Candida albicans) compared to diethyl ether and ethyl acetate-ether. Both extracts (diethyl ether and ethyl acetate) has a minimum antifungal effect while these extracts showed more inhibitory effects on tested bacteria. Inhibitory effects of the aqueous extract were not observed at all concentration (100 microgam mL-1, 200 microgam mL-1and 300 microgam mL-1) against the entire examined bacterial and fungal species. The finding of this study indicates the antibacterial and antifungal activity exhibited by the ethanolic, diethyl ether, and ethyl acetate extract of B. monnieri can be validated to the management of medicinal plant throughout the world.


title: Azimuthal asymmetries of back-to-back $\pi^\pm-(\pi^0,\eta,\pi^\pm)$ pairs in $e^+e^-$ annihilation

abstract: This work reports the first observation of azimuthal asymmetries around the thrust axis in $e^+e^-$ annihilation of pairs of back-to-back charged pions in one hemisphere, and $\pi^0$ and $\eta$ mesons in the opposite hemisphere. These results are complemented by a new analysis of pairs of back-to-back charged pions. The $\pi^0$ and $\eta$ asymmetries rise with the relative momentum $z$ of the detected hadrons as well as with the transverse momentum with respect to the thrust axis. These asymmetries are sensitive to the Collins fragmentation function $H_1^{\perp}$ and provide complementary information to previous measurements with charged pions and kaons in the final state. In particular, the $\eta$ final states will provide additional information on the flavor structure of $H_1^{\perp}$. This is the first measurement of the explicit transverse-momentum dependence of the Collins fragmentation function from Belle data. It uses a dataset of 980.4~fb$^{-1}$ collected by the Belle experiment at or near a center-of-mass energy of 10.58 GeV.


title: Anomalous scaling of dynamical large deviations of stationary Gaussian processes

abstract: Employing the optimal fluctuation method (OFM), we study the large deviation function of long-time averages $(1/T)\int_{-T/2}^{T/2} x^n(t) dt$, $n=1,2, \dots$, of centered stationary Gaussian processes. These processes are correlated and, in general, non-Markovian. We show that the anomalous scaling with time of the large-deviation function, recently observed for $n>2$ for the particular case of the Ornstein-Uhlenbeck process, holds for a whole class of stationary Gaussian processes.


title: A multi-fidelity neural network surrogate sampling method for uncertainty quantification

abstract: We propose a multi-fidelity neural network surrogate sampling method for the uncertainty quantification of physical/biological systems described by ordinary or partial differential equations. We first generate a set of low/high-fidelity data by low/high-fidelity computational models, e.g. using coarser/finer discretizations of the governing differential equations. We then construct a two-level neural network, where a large set of low-fidelity data are utilized in order to accelerate the construction of a high-fidelity surrogate model with a small set of high-fidelity data. We then embed the constructed high-fidelity surrogate model in the framework of Monte Carlo sampling. The proposed algorithm combines the approximation power of neural networks with the advantages of Monte Carlo sampling within a multi-fidelity framework. We present two numerical examples to demonstrate the accuracy and efficiency of the proposed method. We show that dramatic savings in computational cost may be achieved when the output predictions are desired to be accurate within small tolerances.


title: Visual Question Answering using Deep Learning: A Survey and Performance Analysis

abstract: The Visual Question Answering (VQA) task combines challenges for processing data with both Visual and Linguistic processing, to answer basic `common sense' questions about given images. Given an image and a question in natural language, the VQA system tries to find the correct answer to it using visual elements of the image and inference gathered from textual questions. In this survey, we cover and discuss the recent datasets released in the VQA domain dealing with various types of question-formats and robustness of the machine-learning models. Next, we discuss about new deep learning models that have shown promising results over the VQA datasets. At the end, we present and discuss some of the results computed by us over the vanilla VQA model, Stacked Attention Network and the VQA Challenge 2017 winner model. We also provide the detailed analysis along with the challenges and future research directions.


title: Rethinking the Number of Channels for the Convolutional Neural Network

abstract: Latest algorithms for automatic neural architecture search perform remarkable but few of them can effectively design the number of channels for convolutional neural networks and consume less computational efforts. In this paper, we propose a method for efficient automatic architecture search which is special to the widths of networks instead of the connections of neural architecture. Our method, functionally incremental search based on function-preserving, will explore the number of channels rapidly while controlling the number of parameters of the target network. On CIFAR-10 and CIFAR-100 classification, our method using minimal computational resources (0.4~1.3 GPU-days) can discover more efficient rules of the widths of networks to improve the accuracy by about 0.5% on CIFAR-10 and a~2.33% on CIFAR-100 with fewer number of parameters. In particular, our method is suitable for exploring the number of channels of almost any convolutional neural network rapidly.


title: Privacy with Surgical Robotics: Challenges in Applying Contextual Privacy Theory

abstract: The use of connected surgical robotics to automate medical procedures presents new privacy challenges. We argue that conventional patient consent protocols no longer work. Indeed robots that replace human surgeons take on an extraordinary level of responsibility. Surgeons undergo years of training and peer review in a strongly regulated environment, and derive trust via a patient's faith in the hospital system. Robots on the other hand derive trust differently, via the integrity of the software that governs their operation. From a privacy perspective, there are two fundamental shifts. First, the threat model has shifted from one where the humans involved were untrusted to one where the robotic software is untrusted. Second, the basic unit of privacy control is no longer a medical record, but is replaced by four new basic units: the subject on which the robot is taking action; the tools used by the robot; the sensors (i.e data) the robot can access; and, finally access to monitoring and calibration services which afford correct operation of the robot. We suggest that contextual privacy provides useful theoretical tools to solve the privacy problems posed by surgical robots. However, it also poses some challenges: not least that the complexity of the contextual-privacy policies, if rigorously specified to achieve verification and enforceability, will be exceedingly high to directly expose to humans that review contextual privacy policies. A medical robot works with both information and physical material. While informational norms allow for judgements about contextual integrity and the transmission principle governs the constraints applied on information transfer, nothing is said about material property. Certainly, contextual privacy provides an anchor for useful notions of privacy in this scenario and thus should be considered to be extended to cover both information and material flows.


title: Empirical Study of Diachronic Word Embeddings for Scarce Data

abstract: Word meaning change can be inferred from drifts of time-varying word embeddings. However, temporal data may be too sparse to build robust word embeddings and to discriminate significant drifts from noise. In this paper, we compare three models to learn diachronic word embeddings on scarce data: incremental updating of a Skip-Gram from Kim et al. (2014), dynamic filtering from Bamler and Mandt (2017), and dynamic Bernoulli embeddings from Rudolph and Blei (2018). In particular, we study the performance of different initialisation schemes and emphasise what characteristics of each model are more suitable to data scarcity, relying on the distribution of detected drifts. Finally, we regularise the loss of these models to better adapt to scarce data.


title: Holographic fundamental matter in multilayered media

abstract: We describe a strongly coupled layered system in 3+1 dimensions by means of a top-down D-brane construction. Adjoint matter is encoded in a large-$N_c$ stack of D3-branes, while fundamental matter is confined to $(2+1)$-dimensional defects introduced by a large-$N_f$ stack of smeared D5-branes. To the anisotropic Lifshitz-like background geometry, we add a single flavor D7-brane treated in the probe limit. Such bulk setup corresponds to a partially quenched approximation for the dual field theory. The holographic model sheds light on the anisotropic physics induced by the layered structure, allowing one to disentangle flavor physics along and orthogonal to the layers as well as identifying distinct scaling laws for various dynamical quantities. We study the thermodynamics and the fluctuation spectrum with varying valence quark mass or baryon chemical potential. We also focus on the density wave propagation in both the hydrodynamic and collisionless regimes where analytic methods complement the numerics, while the latter provides the only resource to address the intermediate transition regime.


title: Optimal Wireless Resource Allocation with Random Edge Graph Neural Networks

abstract: We consider the problem of optimally allocating resources across a set of transmitters and receivers in a wireless network. The resulting optimization problem takes the form of constrained statistical learning, in which solutions can be found in a model-free manner by parameterizing the resource allocation policy. Convolutional neural networks architectures are an attractive option for parameterization, as their dimensionality is small and does not scale with network size. We introduce the random edge graph neural network (REGNN), which performs convolutions over random graphs formed by the fading interference patterns in the wireless network. The REGNN-based allocation policies are shown to retain an important permutation equivariance property that makes them amenable to transference to different networks. We further present an unsupervised model-free primal-dual learning algorithm to train the weights of the REGNN. Through numerical simulations, we demonstrate the strong performance REGNNs obtain relative to heuristic benchmarks and their transference capabilities.


title: Understanding Bias in Machine Learning

abstract: Bias is known to be an impediment to fair decisions in many domains such as human resources, the public sector, health care etc. Recently, hope has been expressed that the use of machine learning methods for taking such decisions would diminish or even resolve the problem. At the same time, machine learning experts warn that machine learning models can be biased as well. In this article, our goal is to explain the issue of bias in machine learning from a technical perspective and to illustrate the impact that biased data can have on a machine learning model. To reach such a goal, we develop interactive plots to visualizing the bias learned from synthetic data.


title: 3D Bounding Box Estimation for Autonomous Vehicles by Cascaded Geometric Constraints and Depurated 2D Detections Using 3D Results

abstract: 3D object detection is one of the most important tasks in 3D vision perceptual system of autonomous vehicles. In this paper, we propose a novel two stage 3D object detection method aimed at get the optimal solution of object location in 3D space based on regressing two additional 3D object properties by a deep convolutional neural network and combined with cascaded geometric constraints between the 2D and 3D boxes. First, we modify the existing 3D properties regressing network by adding two additional components, viewpoints classification and the center projection of the 3D bounding box s bottom face. Second, we use the predicted center projection combined with similar triangle constraint to acquire an initial 3D bounding box by a closed-form solution. Then, the location predicted by previous step is used as the initial value of the over-determined equations constructed by 2D and 3D boxes fitting constraint with the configuration determined with the classified viewpoint. Finally, we use the recovered physical world information by the 3D detections to filter out the false detection and false alarm in 2D detections. We compare our method with the state-of-the-arts on the KITTI dataset show that although conceptually simple, our method outperforms more complex and computational expensive methods not only by improving the overall precision of 3D detections, but also increasing the orientation estimation precision. Furthermore our method can deal with the truncated objects to some extent and remove the false alarm and false detections in both 2D and 3D detections.


title: Deep learning networks for selection of persistent scatterer pixels in multi-temporal SAR interferometric processing

abstract: In multi-temporal SAR interferometry (MT-InSAR), persistent scatterer (PS) pixels are used to estimate geophysical parameters, essentially deformation. Conventionally, PS pixels are selected on the basis of the estimated noise present in the spatially uncorrelated phase component along with look-angle error in a temporal interferometric stack. In this study, two deep learning architectures, namely convolutional neural network for interferometric semantic segmentation (CNN-ISS) and convolutional long short term memory network for interferometric semantic segmentation (CLSTM-ISS), based on learning spatial and spatio-temporal behaviour respectively, were proposed for selection of PS pixels. These networks were trained to relate the interferometric phase history to its classification into phase stable (PS) and phase unstable (non-PS) measurement pixels using ~10,000 real world interferometric images of different study sites containing man-made objects, forests, vegetation, uncropped land, water bodies, and areas affected by lengthening, foreshortening, layover and shadowing. The networks were trained using training labels obtained from the Stanford method for Persistent Scatterer Interferometry (StaMPS) algorithm. However, pixel selection results, when compared to a combination of R-index and a classified image of the test dataset, reveal that CLSTM-ISS estimates improved the classification of PS and non-PS pixels compared to those of StaMPS and CNN-ISS. The predicted results show that CLSTM-ISS reached an accuracy of 93.50%, higher than that of CNN-ISS (89.21%). CLSTM-ISS also improved the density of reliable PS pixels compared to StaMPS and CNN-ISS and outperformed StaMPS and other conventional MT-InSAR methods in terms of computational efficiency.


title: Generalized Integrated Gradients: A practical method for explaining diverse ensembles

abstract: We introduce Generalized Integrated Gradients (GIG), a formal extension of the Integrated Gradients (IG) (Sundararajan et al., 2017) method for attributing credit to the input variables of a predictive model. GIG improves IG by explaining a broader variety of functions that arise from practical applications of ML in domains like financial services. GIG is constructed to overcome limitations of Shapley (1953) and Aumann-Shapley (1974), and has desirable properties when compared to other approaches. We prove GIG is the only correct method, under a small set of reasonable axioms, for providing explanations for mixed-type models or games. We describe the implementation, and present results of experiments on several datasets and systems of models.


title: Valley-selective Klein tunneling through superlattice barrier in graphene

abstract: Graphene electrons feature a pair of massless Dirac cones of opposite pseudospin chirality at two valleys. Klein tunneling refers to the intriguing capability of these chiral electrons to penetrate through high and wide potential barrier. The two valleys have been treated independently in the literature, where time reversal symmetry dictates that neither the normal incidence transmission nor the angle-averaged one can have any valley polarization. Here we show that, when intervalley scattering by barrier is accounted, graphene electrons normally incident at a superlattice barrier can experience a fully valley-selective Klein tunneling, i.e. perfect transmission in one valley, and perfect reflection in the other. Intervalley backscattering creates staggered pseudospin gaps in the superlattice barrier, which, combined with the valley contrast in pseudospin chirality, determines the valley polarity of Klein tunneling. The angle averaged transmission can have a net valley polarization of 20% for a 5-period barrier, and exceed 75% for a 20-period barrier. Our finding points to an unexpected opportunity to realize valley functionalities in graphene electronics.


title: Help, Anna! Visual Navigation with Natural Multimodal Assistance via Retrospective Curiosity-Encouraging Imitation Learning

abstract: Mobile agents that can leverage help from humans can potentially accomplish more complex tasks than they could entirely on their own. We develop "Help, Anna!" (HANNA), an interactive photo-realistic simulator in which an agent fulfills object-finding tasks by requesting and interpreting natural language-and-vision assistance. An agent solving tasks in a HANNA environment can leverage simulated human assistants, called ANNA (Automatic Natural Navigation Assistants), which, upon request, provide natural language and visual instructions to direct the agent towards the goals. To address the HANNA problem, we develop a memory-augmented neural agent that hierarchically models multiple levels of decision-making, and an imitation learning algorithm that teaches the agent to avoid repeating past mistakes while simultaneously predicting its own chances of making future progress. Empirically, our approach is able to ask for help more effectively than competitive baselines and, thus, attains higher task success rate on both previously seen and previously unseen environments. We publicly release code and data at https://github.com/khanhptnk/hanna . A video demo is available at https://youtu.be/18P94aaaLKg .


title: Tethered single-legged molecular spiders on independent 1D tracks

abstract: We study the motion of random walkers with residence time bias between first and subsequent visits to a site, as a model for synthetic molecular walkers composed of coupled DNAzyme legs known as molecular spiders. The mechanism of the transient superdiffusion has been explained via the emergence of a boundary between the new and the previously visited sites, and the tendency of the multi-legged spider to cling to this boundary, provided residence time for a first visit to a site is longer than for subsequent visits. Using both kinetic Monte Carlo simulation and an analytical approach, we model a system that consists of single-legged walkers, each on its own one-dimensional track, connected by a "leash", i.e., a kinematic constraint that no two spiders can be more than a certain distance apart. Even though a single one-legged walker does not at all exhibit directional, superdiffusive motion, we find that a team of one-legged walkers on parallel tracks, connected by a flexible tether, does enjoy a superdiffusive transient. Furthermore, the one-legged walker teams exhibit a greater expected number of steps per boundary period and are able to diffuse more quickly through the product sea than two-legged walkers, which leads to longer periods of superdiffusion.


title: Sharp pointwise estimates for the gradients of solutions to linear parabolic second order equation in the layer

abstract: We deal with solutions of the Cauchy problem to linear both homogeneous and nonhomogeneous parabolic second order equations with real constant coefficients in the layer ${\mathbb R}^{n+1}_T={\mathbb R}^n\times (0, T)$, where $n\geq 1$ and $T<\infty$. The homogeneous equation is considered with initial data in $L^p({\mathbb R}^n)$, $1\leq p \leq \infty $. For the nonhomogeneous equation we suppose that initial function is equal to zero and the function in the right-hand side belongs to $f\in L^p({\mathbb R}^{n+1}_T)\cap C^\alpha \big (\bar{{\mathbb R}^{n+1}_T} \big ) $ , $p>n+2$ and $\alpha \in (0, 1)$. Explicit formulas for the sharp coefficients in pointwise estimates for the length of the gradient to solutions to these problems are obtained.


title: Optimal Planning and Operation of Multi-Frequency HVac Transmission Systems

abstract: Low-frequency high-voltage alternating-current (LF-HVac) transmission scheme has been recently proposed as an alternative solution to conventional 50/60-Hz HVac and high-voltage direct-current (HVdc) schemes for bulk power transfer. This paper proposes an optimal planning and operation for loss minimization in a multi-frequency HVac transmission system. In such a system, conventional HVac and LF-HVac grids are interconnected using back-to-back (BTB) converters. The dependence of system MW losses on converter dispatch as well as the operating voltage and frequency in the LF-HVac is discussed and compared with that of HVdc transmission. Based on the results of the loss analysis, multi-objective optimization formulations for both planning and operation stages are proposed. The planning phase decides a suitable voltage level for the LF-HVac grid, while the operation phase determines the optimal operating frequency and power dispatch of BTB converters, generators, and shunt capacitors. A solution approach that effectively handles the variations of transmission line parameters with the rated voltage and operating frequency in the LF-HVac grid is proposed. The proposed solutions of the planning and operation stages are evaluated using a multi-frequency HVac system. The results show a significant loss reduction and improved voltage regulation during a 24-hour simulation.


title: Learning Elastic Constitutive Material and Damping Models

abstract: Commonly used linear and nonlinear constitutive material models in deformation simulation contain many simplifications and only cover a tiny part of possible material behavior. In this work we propose a framework for learning customized models of deformable materials from example surface trajectories. The key idea is to iteratively improve a correction to a nominal model of the elastic and damping properties of the object, which allows new forward simulations with the learned correction to more accurately predict the behavior of a given soft object. Space-time optimization is employed to identify gentle control forces with which we extract necessary data for model inference and to finally encapsulate the material correction into a compact parametric form. Furthermore, a patch based position constraint is proposed to tackle the challenge of handling incomplete and noisy observations arising in real-world examples. We demonstrate the effectiveness of our method with a set of synthetic examples, as well with data captured from real world homogeneous elastic objects.


title: Global existence of the solution to Einstein-Yang-Mills-Higgs equations with small initial datum

abstract: The problem involved in this paper is the global existence of the solution to the $\mathfrak{su}(2)$-Einstein-Yang-Mills-Higgs(EYMH) equation. The approach we employ stems from H. Lindblad and I. Rodnianski and is dependent of wave coordinates and Lorentzian gauge conditions. Our main conclusion is that the EYMH system admits global existence provided the initial datum are sufficiently small. To the best of our knowledge, there is no similar result in the area of EYMH equations.


title: On diagram groups over Fibonacci-like semigroup presentations and their generalizations

abstract: We answer the question by Matt Brin on the structure of diagram groups over semigroup presentation ${\mathcal P}=\langle a,b,c\mid a=bc,b=ca,c=ab\rangle$. In the talk on Oberwolfach workshop, Brin conjectured that the diagram group over $\mathcal P$ with base $a$ is isomorphic to the generalized Thompson's group $F_9$. We confirm this conjecture and consider some generalizations of this fact.


title: Local Hole revisited: evidence for bulk motions and self-consistent outflow

abstract: We revisit our mapping of the `Local Hole', a large underdensity in the local galaxy redshift distribution that extends out to redshift, $z\approx0.05$ and a potential source of outflows that may perturb the global expansion rate and thus help mitigate the present `$H_0$ tension'. First, we compare local peculiar velocities measured via the galaxy average redshift-magnitude Hubble diagram, $\bar{z}(m)$, with a simple dynamical outflow model based on the average underdensity in the Local Hole. We find that this outflow model is in good agreement with our peculiar velocity measurements from $\bar{z}(m)$ and not significantly inconsistent with SNIa peculiar velocity measurements from at least the largest previous survey. This outflow could cause an $\approx2-3$\% increase in the local value of Hubble's constant. Second, considering anisotropic motions, we find that the addition of the outflow model may improve the $\bar{z}(m)$ fit of a bulk flow where galaxies are otherwise at rest in the Local Group frame. We conclude that the Local Hole plus neighbouring overdensities such as the Shapley Supercluster may cause outflow and bulk motions out to $\approx150$h$^{-1}$Mpc that are cosmologically significant and that need to be taken into account in estimating Hubble's constant.


title: Polyhedral products, and relations in the commutator subgroup of a right-angled Coxeter group

abstract: We give a criterion of the existence of a presentation with a single relation for the commutator subgroup $RC_{\mathscr{K}}'$ of a right-angled Coxeter group $RC_{\mathscr{K}}$. Namely, we prove that $RC_{\mathscr{K}}'$ is a one-relator group if and only if $\mathscr{K}$ is either a $p$-cycle with $p\geqslant4$ or has the form ($p$-cycle)$*\Delta^q$ with $p\geqslant4$ and $g\geqslant0$, where $\Delta^q$ is a $q$-simplex.


title: Covariant Relativistic Non-Equilibrium Thermodynamics of Multi-Component Systems

abstract: Non-equilibrium and equilibrium thermodynamics of an interacting component in a relativistic multi-component system is discussed covariantly by exploiting an entropy identity. The special case of the corresponding free component is considered. Equilibrium conditions and especially the multi-component Killing relation of the 4-temperature are discussed. Two axioms characterize the mixture: additivity of the energy momentum tensors and additivity of the 4-entropies of the components generating those of the mixture. The resulting quantities of a single component and of the mixture as a whole, energy, energy flux, momentum flux, stress tensor, entropy, entropy flux, supply and production are derived. Finally, a general relativistic 2-component mixture is discussed with respect to their gravitation generating energy-momentum tensors.


title: Scattering of solutions to NLW by Inward Energy Decay

abstract: The topic of this paper is a semi-linear, energy sub-critical, defocusing wave equation $\partial_t^2 u - \Delta u = - |u|^{p -1} u$ in the 3-dimensional space ($3\leq p<5$) whose initial data are radial and come with a finite energy. In this work we prove scattering in the positive time direction by only assuming the inward part of the energy decays at a certain rate, as long as the total energy is finite, regardless of the decay rate or size of the outward energy. More precisely, we assume the initial data comes with a finite energy and \[ \int_{{\mathbb R}^3} \max\{1,|x|^\kappa\}\ (\ |\nabla u_0(x)\cdot \frac{x}{|x|} + \frac{u_0(x)}{|x|} + u_1(x)\ |^2 + \frac{2}{p+1}|u_0(x)|^{p+1}\ ) dx < \infty. \] Here $\kappa\geq \kappa_0(p) = \frac{5-p}{p+1}$ is a constant. If $\kappa>\kappa_0(p)$, we can also prove $\|u\|_{L^p L^{2p}}({\mathbb R}^+ \times {\mathbb R}^3)< +\infty$ and give an explicit rate of $u$'s convergence to a free wave.


title: On the density of Cayley graphs of R.Thompson's group $F$ in symmetric generators

abstract: By the density of a finite graph we mean its average vertex degree. For an $m$-generated group, the density of its Cayley graph in a given set of generators, is the supremum of densities taken over all its finite subgraphs. It is known that a group with $m$ generators is amenable iff the density of the corresponding Cayley graph equals $2m$. A famous problem on the amenability of R.\,Thompson's group $F$ is still open. What is known due to the result by Belk and Brown, is that the density of its Cayley graph in the standard set of group generators $\{x_0,x_1\}$, is at least $3.5$. This estimate has not been exceeded so far. For the set of symmetric generators $S=\{x_1,\bar{x}_1\}$, where $\bar{x}_1=x_1x_0^{-1}$, the same example gave the estimate only $3$. There was a conjecture that for this generating set the equality holds. If so, $F$ would be non-amenable, and the symmetric generating set had doubling property. This means that for any finite set $X\subset F$, the inequality $|S^{\pm1}X|\ge2|X|$ holds. In this paper we disprove this conjecture showing that the density of the Cayley graph of $F$ in symmetric generators $S$ strictly exceeds $3$. Moreover, we show that even larger generating set $S_0=\{x_0,x_1,\bar{x}_1\}$ does not have doubling property.


title: Projectively self-concordant barriers

abstract: Self-concordance is the most important property required for barriers in convex programming. It is intrinsically linked to the affine structure of the underlying space. Here we introduce an alternative notion of self-concordance which is linked to the projective structure. A function on a set $X \subset A^n$ in an $n$-dimensional affine space is projectively self-concordant if and only if it can be extended to an affinely self-concordant logarithmically homogeneous function on the conic extension $K \subset V^{n+1}$ of $X$ in the $(n+1)$-dimensional vector space obtained by homogenization of $A^n$. The feasible sets in conic programs, notably linear and semi-definite programs, are naturally equipped with projectively self-concordant barriers. However, the interior-point methods used to solve these programs employ only affine self-concordance. We show that estimates used in the analysis of interior-point methods are tighter for projective self-concordance, in particular inner and outer approximations of the set. This opens the way to a better tuning of parameters in interior-points algorithms to allow larger steps and hence faster convergence. Projective self-concordance is also a useful tool in the theoretical analysis of logarithmically homogeneous barriers on cones.


title: Does the ratio of Laplace transforms of powers of a function identify the function?

abstract: We study the following question: if $f$ is a nonzero measurable function on $[0,\infty)$ and $m$ and $n$ distinct nonnegative integers, does the ratio $\widehat{f^n}/\widehat{f^m}$ of the Laplace transforms of the powers $f^n$ and $f^m$ of $f$ uniquely determine $f$? The answer is yes if one of $m, n$ is zero, by the inverse Laplace transform. Under some assumptions on the smoothness of $f$ we show that the answer in the general case is also affirmative. The question arose from a problem in economics, specifically in auction theory where $f$ is the cumulative distribution function of a certain random variable. This is also discussed in the paper.


title: Agent-based model for tumour-analysis using Python+Mesa

abstract: The potential power provided and possibilities presented by computation graphs has steered most of the available modeling techniques to re-implementing, utilization and including the complex nature of System Biology (SB). To model the dynamics of cellular population, we need to study a plethora of scenarios ranging from cell differentiation to tumor growth and etcetera. Test and verification of a model in research means running the model multiple times with different or in some cases identical parameters, to see how the model interacts and if some of the outputs would change regarding different parameters. In this paper, we will describe the development and implementation of a new agent-based model using Python. The model can be executed using a development environment (based on Mesa, and extremely simplified for convenience) with different parameters. The result is collecting large sets of data, which will allow an in-depth analysis in the microenvironment of the tumor by the means of network analysis.


title: Magnon polaron formed by selectively coupled coherent magnon and phonon modes of a surface patterned ferromagnet

abstract: Strong coupling between two quanta of different excitations leads to the formation of a hybridized state which paves a way for exploiting new degrees of freedom to control phenomena with high efficiency and precision. A magnon polaron is the hybridized state of a phonon and a magnon, the elementary quanta of lattice vibrations and spin waves in a magnetically-ordered material. A magnon polaron can be formed at the intersection of the magnon and phonon dispersions, where their frequencies coincide. The observation of magnon polarons in the time domain has remained extremely challenging because the weak interaction of magnons and phonons and their short lifetimes jeopardize the strong coupling required for the formation of a hybridized state. Here, we overcome these limitations by spatial matching of magnons and phonons in a metallic ferromagnet with a nanoscale periodic surface pattern. The spatial overlap of the selected phonon and magnon modes formed in the periodic ferromagnetic structure results in a high coupling strength which, in combination with their long lifetimes allows us to find clear evidence of an optically excited magnon polaron. We show that the symmetries of the localized magnon and phonon states play a crucial role in the magnon polaron formation and its manifestation in the optically excited magnetic transients.


title: Optimal translational-rotational invariant dictionaries for images

abstract: We provide the construction of a set of square matrices whose translates and rotates provide a Parseval frame that is optimal for approximating a given dataset of images. Our approach is based on abstract harmonic analysis techniques. Optimality is considered with respect to the quadratic error of approximation of the images in the dataset with their projection onto a linear subspace that is invariant under translations and rotations. In addition, we provide an elementary and fully self-contained proof of optimality, and the numerical results from datasets of natural images.

categories: eess.IV cs.CV cs.NA math.FA math.NA stat.ML


title: Scoring Strategic Agents

abstract: I introduce a model of predictive scoring. A receiver wants to predict a sender's quality. An intermediary observes multiple features of the sender and aggregates them into a score. Based on the score, the receiver makes a decision. The sender prefers "higher" decisions, and she can distort each feature at a privately known cost. I characterize the scoring rule that maximizes decision accuracy. This rule underweights some features to deter sender distortion, and overweights other features so that the score is correct on average. The receiver prefers this scoring rule to full disclosure because it mitigates his commitment problem.


title: Illiquid Financial Markets and Monetary Policy

abstract: This paper analyzes the role of money in asset markets characterized by search frictions. We develop a dynamic framework that brings together a model for illiquid financial assets `a la Duffie, Garleanu, and Pedersen, and a search-theoretic model of monetary exchange `a la Lagos and Wright. The presence of decentralized financial markets generates an essential role for money, which helps investors re-balance their portfolios. We provide conditions that guarantee the existence of a monetary equilibrium. In this case, asset prices are always above their fundamental value, and this differential represents a liquidity premium. We are able to derive an asset pricing theory that delivers an explicit connection between monetary policy, asset prices, and welfare. We obtain a negative relationship between inflation and equilibrium asset prices. This key result stems from the complementarity between money and assets in our framework.


title: Exact non-symmetric closed line vacuum magnetic fields in a topological torus

abstract: Non-symmetric vacuum magnetic fields with closed magnetic field lines are of interest in the construction of stellarator equilibria. Beyond the result of D.Lortz (ZAMP \textbf{21}, 196 (1970)), few results are available. This work presents a closed-form expression for a class of vacuum magnetic fields in a topological torus with closed field lines. We explicitly obtain the invariants of such a field. We finally show that a three-dimensional low beta magnetohydrodynamic (MHD) equilibrium may be constructed in a topological torus starting with these closed line vacuum magnetic fields.


title: Let's agree to disagree: learning highly debatable multirater labelling

abstract: Classification and differentiation of small pathological objects may greatly vary among human raters due to differences in training, expertise and their consistency over time. In a radiological setting, objects commonly have high within-class appearance variability whilst sharing certain characteristics across different classes, making their distinction even more difficult. As an example, markers of cerebral small vessel disease, such as enlarged perivascular spaces (EPVS) and lacunes, can be very varied in their appearance while exhibiting high inter-class similarity, making this task highly challenging for human raters. In this work, we investigate joint models of individual rater behaviour and multirater consensus in a deep learning setting, and apply it to a brain lesion object-detection task. Results show that jointly modelling both individual and consensus estimates leads to significant improvements in performance when compared to directly predicting consensus labels, while also allowing the characterization of human-rater consistency.


title: Sur la repr\'esentation des entiers par les formes cyclotomiques de grand degr\'e

abstract: For each integer $d\ge 4$, we study the sequence of positive integers which are represented by one at least of the cyclotomic binary forms $\Phi_n(X,Y)$, with $n$ a positive integer satisfying $\varphi(n)\ge d$. The case $d=2$ was studied in our previous work [FLW]. Our demonstration is based on a variant of a statement of [SX] concerning the common values taken by two binary forms of the same degree and non-zero discriminants. All constants are effectively calculable.


title: About Fibonacci trees III: multiple Fibonacci trees

abstract: In this third paper, we revisit the question to which extent the properties of the trees associated to the tilings $\{p,4\}$ of the hyperbolic plane are still true if we consider a finitely generated tree by the same rules but rooted at a black node? What happens if, considering the same distinction between black and white nodes but changing the place of the black son in the rules. What happens if we change the representation of the numbers by another set of digits? We tackle all of these questions in the paper. The present paper is an extension of the previous papers arXiv:1904.12135 and arXiv:1907.04677.


title: Towards the Web of Quantum Chaos Diagnostics

abstract: We study the connections between three quantities that can be used as diagnostics for quantum chaos, i.e., the out-of-time-order correlator (OTOC), Loschmidt echo (LE), and complexity. We generalize the connection between OTOC and LE for infinite dimensions and extend it for higher-order OTOCs and multi-fold LEs. Novel applications of this intrinsic relation are proposed. We also propose a relationship between a specific circuit complexity and LE by using the inverted oscillator model. These relationships signal a deeper connection between these three probes of quantum chaos.


title: Learning a Spatial Field in Minimum Time with a Team of Robots

abstract: We study an informative path-planning problem where the goal is to minimize the time required to learn a spatially varying entity. We use Gaussian Process (GP) regression for learning the underlying field. Our goal is to ensure that the GP posterior variance, which is also the mean square error between the learned and actual fields, is below a predefined value. We study three versions of the problem. In the placement version, the objective is to minimize the number of measurement locations while ensuring that the posterior variance is below a predefined threshold. In the mobile robot version, we seek to minimize the total time required to visit and collect measurements from the measurement locations using a single robot. We also study a multi-robot version where the objective is to minimize the time required by the last robot to return to a common starting location called depot. By exploiting the properties of GP regression, we present constant-factor approximation algorithms. In addition to the theoretical results, we also compare the empirical performance using a real-world dataset, with other baseline strategies.


title: Some lattice models with hyperbolic chaotic attractors

abstract: Examples of one-dimensional lattice systems are considered, in which patterns of different spatial scales arise alternately, so that the spatial phase over a full cycle undergo transformation according to expanding circle map that implies occurrence of Smale-Williams attractors in the multidimensional state space. These models can serve as a basis for design electronic generators of robust chaos within a paradigm of coupled cellular networks. One of the examples is a mechanical pendulum system interesting and demonstrative for research and educational experimental studies.


title: Multi-epoch X-ray spectral analysis of the narrow-line Seyfert 1 galaxy Mrk 478

abstract: A multi-epoch X-ray spectral and variability analysis is conducted for the narrow-line Seyfert 1 (NLS1) active galactic nucleus (AGN) Mrk 478. All available X-ray data from XMM-Newton and Suzaku satellites, spanning from 2001 to 2017, are modelled with a variety of physical models including partial covering, soft-Comptonisation, and blurred reflection, to explain the observed spectral shape and variability over the 16 years. All models are a similar statistical fit to the data sets, though the analysis of the variability between data sets favours the blurred reflection model. In particular, the variability can be attributed to changes in flux of the primary coronal emission. Different reflection models fit the data equally well, but differ in interpretation. The use of reflionx predicts a low disc ionisation and power law dominated spectrum, while relxill predicts a highly ionised and blurred reflection dominated spectrum. A power law dominated spectrum might be more consistent with the normal X-ray-to-UV spectral shape (aox). Both blurred reflection models suggest a rapidly spinning black hole seen at a low inclination angle, and both require a sub-solar (~0.5) abundance of iron. All physical models require a narrow emission feature at 6.7 keV likely attributable to Fe xxv emission, while no evidence for a narrow 6.4 keV line from neutral iron is detected.


title: On the generic part of the cohomology of non-compact unitary Shimura varieties

abstract: We prove that the generic part of the mod l cohomology of Shimura varieties associated to quasi-split unitary groups of even dimension is concentrated above the middle degree, extending our previous work to a non-compact case. The result applies even to Eisenstein cohomology classes coming from the locally symmetric space of the general linear group, and has been used in joint work with Allen, Calegari, Gee, Helm, Le Hung, Newton, Taylor and Thorne to get good control on these classes and deduce potential automorphy theorems without any self-duality hypothesis. Our main geometric result is a computation of the fibers of the Hodge-Tate period map on compactified Shimura varieties, in terms of similarly compactified Igusa varieties.


title: Experimental probe of a complete 3D photonic band gap

abstract: The identification of a complete three-dimensional (3D) photonic band gap in real crystals always employs theoretical or numerical models that invoke idealized crystal structures. Thus, this approach is prone to false positives (gap wrongly assigned) or false negatives (gap missed). Therefore, we propose a purely experimental probe of the 3D photonic band gap that pertains to many different classes of photonic materials. We study position and polarization-resolved reflectivity spectra of 3D inverse woodpile structures that consist of two perpendicular nanopore arrays etched in silicon. We observe intense reflectivity peaks $(R > 90\%)$ typical of high-quality crystals with broad stopbands. We track the stopband width versus pore radius, which agrees much better with the predicted 3D photonic band gap than with a directional stop gap on account of the large numerical aperture used. A parametric plot of s-polarized versus p-polarized stopband width agrees very well with the 3D band gap and is model-free. This practical probe provides fast feedback on the advanced nanofabrication needed for 3D photonic crystals and stimulates practical applications of band gaps in 3D silicon nanophotonics and photonic integrated circuits, photovoltaics, cavity QED, and quantum information processing.


title: Efficient Verification of Pure Quantum States in the Adversarial Scenario

abstract: Efficient verification of pure quantum states in the adversarial scenario is crucial to many applications in quantum information processing, such as blind measurement-based quantum computation and quantum networks. However, little is known about this topic so far. Here we establish a general framework for verifying pure quantum states in the adversarial scenario and clarify the resource cost. Moreover, we propose a simple and general recipe to constructing efficient verification protocols for the adversarial scenario from protocols for the nonadversarial scenario. With this recipe, arbitrary pure states can be verified in the adversarial scenario with almost the same efficiency as in the nonadversarial scenario. Many important quantum states can be verified in the adversarial scenario using local projective measurements with unprecedented high efficiencies.


title: On total Springer representations for the symplectic Lie algebra in characteristic 2 and the exotic case

abstract: Let $W$ be the Weyl group of type $BC_n$. We first provide restriction formulas of the total Springer representations for the symplectic Lie algebra in characteristic 2 and the exotic case to the maximal parabolic subgroup of $W$ which is of type $BC_{n-1}$. Then we show that these two restriction formulas are equivalent, and discuss how the results can be used to examine the existence of affine pavings of Springer fibers corresponding to the symplectic Lie algebra in characteristic 2.


title: Isotropization and change of complexity by gravitational decoupling

abstract: We employ the gravitational decoupling approach for static and spherically symmetric systems to develop a simple and powerful method in order to a) continuously isotropize any anisotropic solution of the Einstein field equations, and b) generate new solutions for self-gravitating distributions with the same or vanishing complexity factor. A few working examples are given for illustrative purposes.


title: Performance of a temporally multiplexed single-photon source with imperfect devices

abstract: Scalable photonic quantum technologies require highly efficient sources of single photons on demand. Although much progress has been done in the field within the last decade, the requirements impose stringent conditions on the efficiency of such devices. One of the most promising approaches is to multiplex a single or several heralded photon sources into temporal modes. In this work we analyze a specific proposal to synchronize photons from a continuous source with an external reference clock using imperfect optical switches, which necessarily degrade the ideal behavior of the devised arrangement. The performance of the source as a sub-poissonian light emitter is studied taking into account losses in the multiplexing arrangement, detector efficiency and dark counts. We estimate a fivefold increase in the single photon probability achieved for 0.5 dB loss switches.


title: VoIPLoc: Passive VoIP call provenance via acoustic side-channels

abstract: We propose VoIPLoc, a novel location fingerprinting technique and apply it to the VoIP call provenance problem. It exploits echo-location information embedded within VoIP audio to support fine-grained location inference. We found consistent statistical features induced by the echo-reflection characteristics of the location into recorded speech. These features are discernible within traces received at the VoIP destination, enabling location inference. We evaluated VoIPLoc by developing a dataset of audio traces received through VoIP channels over the Tor network. We show that recording locations can be fingerprinted and detected remotely with a low false-positive rate, even when a majority of the audio samples are unlabelled. Finally, we note that the technique is fully passive and thus undetectable, unlike prior art. VoIPLoc is robust to the impact of environmental noise and background sounds, as well as the impact of compressive codecs and network jitter. The technique is also highly scalable and offers several degrees of freedom terms of the fingerprintable space.


title: Studying the variability of the X-ray spectral parameters of high-redshift GRBs' afterglows

abstract: The Swift satellite has observed more than a thousand GRBs with X-ray data. Almost a third of them have redshift measurement, too. Here we start to investigate the X-ray spectral fitting of the data considering the low energy part where the N(H) absorption happens. Based on the available more accurate input data we examined the robustness of previous fittings and tested how sensitive the changes of the starting parameters are. We studied the change of the intrinsic hydrogen column density during the outburst for a few events. No significant variability of N(H) column density was identified.


title: Relative entropy for von Neumann subalgebras

abstract: We revisit the connection between index and relative entropy for an inclusion of finite von Neumann algebras. We observe that the Pimsner-Popa index connects to sandwiched Renyi $p$-relative entropy for all $1/2\le p\le \infty$, including Umegaki's relative entropy at $p=1$. Based on that, we introduce a new notation of relative entropy to a subalgebra which generalizes subfactors index. This relative entropy has application in estimating decoherence time of quantum Markov semigroups.


title: Multiple stellar populations in NGC 1866. New clues from Cepheids and Colour-Magnitude Diagram

abstract: We perform a comprehensive study of the stellar populations in the young Large Magellanic Cloud cluster NGC 1866, combining the analysis of its best-studied Cepheids with that of a very accurate colour-magnitude diagram (CMD) obtained from the most recent Hubble Space Telescope photometry. We use a Bayesian method based on new PARSEC stellar evolutionary tracks with overshooting and rotation to obtain ages and initial rotation velocities of five well-studied Cepheids of the cluster. We find that four of the five Cepheids belong to an initially slowly rotating young population (of $ 176 \pm 5$ Myr), while the fifth one is significantly older, either $ 288 \pm 20$ Myr for models with high initial rotational velocity ($\omega_\mathrm{i} \sim 0.9$), or $ 202 \pm 5$ Myr for slowly rotating models. The complementary analysis of the CMD rules out the latter solution while strongly supporting the presence of two distinct populations of $\sim$176 Myr and $\sim$288 Myr, respectively. Moreover, the observed multiple main sequences and the turn-offs indicate that the younger population is mainly made of slowly rotating stars, as is the case of the four younger Cepheids, while the older population is made mainly of initially fast rotating stars, as is the case of the fifth Cepheid. Our study not only reinforces the notion that some young clusters like NGC 1866 harbor multiple populations, but gives also hints that the first population, the older, may inherit the angular momentum from the parent cloud while stars of the second one, the younger, do not.


title: A Microscopic Theory of Intrinsic Timescales in Spiking Neural Networks

abstract: A complex interplay of single-neuron properties and the recurrent network structure shapes the activity of cortical neurons. The single-neuron activity statistics differ in general from the respective population statistics, including spectra and, correspondingly, autocorrelation times. We develop a theory for self-consistent second-order single-neuron statistics in block-structured sparse random networks of spiking neurons. In particular, the theory predicts the neuron-level autocorrelation times, also known as intrinsic timescales, of the neuronal activity. The theory is based on an extension of dynamic mean-field theory from rate networks to spiking networks, which is validated via simulations. It accounts for both static variability, e.g. due to a distributed number of incoming synapses per neuron, and temporal fluctuations of the input. We apply the theory to balanced random networks of generalized linear model neurons, balanced random networks of leaky integrate-and-fire neurons, and a biologically constrained network of leaky integrate-and-fire neurons. For the generalized linear model network with an error function nonlinearity, a novel analytical solution of the colored noise problem allows us to obtain self-consistent firing rate distributions, single-neuron power spectra, and intrinsic timescales. For the leaky integrate-and-fire networks, we derive an approximate analytical solution of the colored noise problem, based on the Stratonovich approximation of the Wiener-Rice series and a novel analytical solution for the free upcrossing statistics. Again closing the system self-consistently, in the fluctuation-driven regime this approximation yields reliable estimates of the mean firing rate and its variance across neurons, the inter-spike interval distribution, the single-neuron power spectra, and intrinsic timescales.


title: On the geometry of K3 surfaces with finite automorphism group and no elliptic fibrations

abstract: Nikulin and Vinberg proved that there are only a finite number of lattices of rank $\geq 3$ that are the N\'eron-Severi group of projective K3 surfaces with a finite automorphism group. The aim of this paper is to provide a more geometric description of such K3 surfaces $X$, when these surfaces have moreover no elliptic fibrations. In that case we show that such K3 surface is either a quartic with special hyperplane sections or a double cover of the plane branched over a smooth sextic curve which has special tangencies properties with some lines, conics or cuspidal cubic curves. We then study the converse i.e. if the geometric description we obtained characterizes these surfaces. In $4$ cases the description is sufficient, in each of the $4$ other cases there is exactly another one possibility which we study. We obtain that at least 5 moduli spaces of K3 surfaces (among the 8 we study) are unirational.


title: ICSrange: A Simulation-based Cyber Range Platform for Industrial Control Systems

abstract: Maintenance staff of Industrial Control Systems (ICS) is generally not aware about information technologies, and even less about cyber security problems. The scary impact of cyber attacks in the industrial world calls for tools to train defensive skills and test effective security measures. Cyber range offers this opportunity, but current research is lacking cost-effective solutions verticalized for the industrial domain. This work proposes ICSrange, a simulation-based cyber range platform for Industrial Control Systems. ICSrange adopts Commercial-Off-The-Shelf (COTS) technologies to virtualize an enterprise network connected to Industrial Control Systems. ICSrange is the outcome of a preliminary study intended to investigate challenges and opportunities to build a configurable and extensible cyber range with simulated industrial processes. Literature shows that testbeds based on realistic mock-ups are effectively employed to develop complex exploits like Advanced Persistent Threats (APTs), hence motivating their usage to train and test security in ICS. We prove the effectiveness of ICSrange through the execution of a multi-staged attack that breaches an enterprise network and progressively intrudes a simulated ICS with water tanks. The attack mimics lateral movements as observed in APTs.


title: Formation of Power-law Electron Energy Spectra in Three-dimensional Low-$\beta$ Magnetic Reconnection

abstract: While observations have suggested that power-law electron energy spectra are a common outcome of strong energy release during magnetic reconnection, e.g., in solar flares, kinetic simulations have not been able to provide definite evidence of power-laws in energy spectra of non-relativistic reconnection. By means of 3D large-scale fully kinetic simulations, we study the formation of power-law electron energy spectra in non-relativistic low-$\beta$ reconnection. We find that both the global spectrum integrated over the entire domain and local spectra within individual regions of the reconnection layer have power-law tails with a spectral index $p \sim 4$ in the 3D simulation, which persist throughout the non-linear reconnection phase until saturation. In contrast, the spectrum in the 2D simulation rapidly evolves and quickly becomes soft. We show that 3D effects such as self-generated turbulence and chaotic magnetic field lines enable the transport of high-energy electrons across the reconnection layer and allow them to access several main acceleration regions. This leads to a sustained and nearly constant acceleration rate for electrons at different energies. We construct a model that explains the observed power-law spectral index in terms of the dynamical balance between particle acceleration and escape from main acceleration regions, which are defined based upon a threshold for the curvature drift acceleration term. This result could be important for explaining the formation of power-law energy spectrum in solar flares.

categories: astro-ph.SR astro-ph.HE physics.plasm-ph physics.space-ph


title: Subradiance-protected excitation spreading in the generation of collimated photon emission from an atomic array

abstract: We show how an initial localized radiative excitation in a two-dimensional array of cold atoms can be converted into highly-directional coherent emission of light by protecting the spreading of the excitation across the array in a subradiant collective eigenmode with a lifetime orders of magnitude longer than that of an isolated atom. We demonstrate how to reach two such strongly subradiant modes, a uniform one where all the dipoles are oscillating in phase normal to the plane and an antiferromagnetic mode where each dipole is $\pi$ out of phase with its nearest neighbor. The excitation, which can consist of a single photon, is then released from the protected subradiant eigenmode by controlling the Zeeman level shifts of the atoms. Hence, an original localized excitation which emits in all directions is transferred to a delocalized subradiance-protected excitation, with a probabilistic emission of a photon only along the axis perpendicular to the plane of the atoms. This protected spreading and directional emission could potentially be used to link stages in a quantum information or quantum computing architecture.

categories: cond-mat.quant-gas physics.atom-ph physics.optics quant-ph


title: The Effect of the Approach to Gas Disk Gravitational Instability on the Rapid Formation of Gas Giant Planets

abstract: Observational evidence suggests that gas disk instability may be responsible for the formation of at least some gas giant exoplanets, particularly massive or distant gas giants. With regard to close-in gas giants, Boss (2017) used the $\beta$ cooling approximation to calculate hydrodynamical models of inner gas disk instability, finding that provided disks with low values of the initial minimum Toomre stability parameter (i.e., $Q_i < 2$ inside 20 au) form, fragmentation into self-gravitating clumps could occur even for $\beta$ as high as 100 (i.e., extremely slow cooling). Those results implied that the evolution of disks toward low $Q_i$ must be taken into account. This paper presents such models: initial disk masses of 0.091 $M_\odot$ extending from 4 to 20 au around a 1 $M_\odot$ protostar, with a range (1 to 100) of $\beta$ cooling parameters, the same as in Boss (2017), but with all the disks starting with $Q_i = 2.7$, i.e., gravitationally stable, and allowed to cool from their initial outer disk temperature of 180 K to as low as 40 K. All the disks eventually fragment into at least one dense clump. The clumps were again replaced by virtual protoplanets (VPs) and the masses and orbits of the resulting ensemble of VPs compare favorably with those of Boss (2017), supporting the claim that disk instability can form gas giants rapidly inside 20 au, provided that sufficiently massive protoplanetary disks exist.


title: Rayleigh-Taylor unstable flames at higher Reynolds number

abstract: Rayleigh-Taylor (RT) unstable flames are a key component of Type Ia and Iax supernovae explosions, but their complex hydrodynamics is still not well understood. These flames are affected not only by the RT instability, but also by the turbulence it generates. Both processes can increase the flame speed by stretching and wrinkling the flame. This makes it hard to choose a subgrid model for the flame speed in full star Type Ia or Iax simulations. Commonly used subgrid models get around this difficulty by assuming that either the RT instability or turbulence is dominant and sets the flame speed. In previous work, we evaluated the physical assumptions and predictive abilities of these two types of models by analysing a large parameter study of 3D direct numerical simulations of RT unstable flames. Surprisingly, we found that the flame dynamics is dominated by the RT instability and that RT unstable flames are very different from turbulent flames. In particular, RT unstable flames are thinner rather than thicker when turbulence is strong. In addition, none of the turbulent flame speed models adequately predicted the flame speed. We also showed that the RT flame speed model failed when the RT instability was strong, suggesting that geometrical burning effects also influence the flame speed. However, these results depended on simulations with $Re\lesssim720$. In this paper, we extend the parameter study to higher Reynolds number and show that the basic conclusions of our previous study still hold when the RT-generated turbulence is stronger.


title: The influence of diversity on the measurement of functional impairment: An international validation of the Amsterdam IADL Questionnaire in 8 countries

abstract: INTRODUCTION: To understand the potential influence of diversity on the measurement of functional impairment in dementia, we aimed to investigate possible bias caused by age, gender, education, and cultural differences. METHODS: 3,571 individuals (67.1 {\pm} 9.5 years old, 44.7% female) from the Netherlands, Spain, France, United States, United Kingdom, Greece, Serbia and Finland were included. Functional impairment was measured using the Amsterdam IADL Questionnaire. Item bias was assessed using differential item functioning (DIF) analysis. RESULTS: There were some differences in activity endorsement. A few items showed statistically significant DIF. However, there was no evidence of meaningful item bias: effect sizes were low ({\Delta}R2 range 0-0.03). Impact on total scores was minimal. DISCUSSION: The results imply a limited bias for age, gender, education and culture in the measurement of functional impairment. This study provides an important step in recognizing the potential influence of diversity on primary outcomes in dementia research.


title: Space- and time-crystallization effects in multicomponent superfluids

abstract: We observe that space- and time-crystallization effects in multicomponent superfluids---while having the same physical origin and mathematical description as in the single-component case---are conceptually much more straightforward. Specifically, the values of the temporal and spatial periods are absolute rather than relative, and the broken translation symmetry in space and/or time can be revealed with experiments involving only one equilibrium sample. We discuss two realistic setups---one with cold atoms and another one with bilayer superconductors---for observation of space and time crystallization in two-component counterflow superfluids.


title: Differentially Private SQL with Bounded User Contribution

abstract: Differential privacy (DP) provides formal guarantees that the output of a database query does not reveal too much information about any individual present in the database. While many differentially private algorithms have been proposed in the scientific literature, there are only a few end-to-end implementations of differentially private query engines. Crucially, existing systems assume that each individual is associated with at most one database record, which is unrealistic in practice. We propose a generic and scalable method to perform differentially private aggregations on databases, even when individuals can each be associated with arbitrarily many rows. We express this method as an operator in relational algebra, and implement it in an SQL engine. To validate this system, we test the utility of typical queries on industry benchmarks, and verify its correctness with a stochastic test framework we developed. We highlight the promises and pitfalls learned when deploying such a system in practice, and we publish its core components as open-source software.


title: On Orthogonal Vector Edge Coloring

abstract: Given a graph $G$ and a positive integer $d$, an orthogonal vector $d$-coloring of $G$ is an assignment $f$ of vectors of $\mathbb{R}^d$ to $V(G)$ in such a way that adjacent vertices receive orthogonal vectors. The orthogonal chromatic number of $G$, denoted by $\chi_v(G)$, is the minimum $d$ for which $G$ admits an orthogonal vector $d$-coloring. This notion has close ties with the notions of Lov\'asz Theta Function, quantum chromatic number, and many other problems, and even though this and related metrics have been extensively studied over the years, we have found that there is a gap in the knowledge concerning the edge version of the problem. In this article, we discuss this version and its relation with other insteresting known facts, and pose a question about the orthogonal chromatic index of cubic graphs.


title: Mape_Maker: A Scenario Creator

abstract: We describe algorithms for creating probabilistic scenarios for the situation when the underlying forecast methodology is modeled as being more (or less) accurate than it has been historically. Such scenarios can be used in studies that extend into the future and may need to consider the possibility that forecast technology will improve. Our approach can also be used to generate alternative realizations of renewable energy production that are consistent with historical forecast accuracy, in effect serving as a method for creating families of realistic alternatives -- which are often critical in simulation-based analysis methodologies


title: Ramsey numbers of path-matchings, covering designs and 1-cores

abstract: A path-matching of order $p$ is a vertex disjoint union of nontrivial paths spanning $p$ vertices. Burr and Roberts, and Faudree and Schelp determined the 2-color Ramsey number of path-matchings. In this paper we study the multicolor Ramsey number of path-matchings. Given positive integers $r, p_1, \dots, p_r$, define $R^{PM}(p_1, \dots, p_r)$ to be the smallest integer $n$ such that in any $r$-coloring of the edges of $K_n$ there exists a path-matching of color $i$ and order at least $p_i$ for some $i\in [r]$. Our main result is that for $r\geq 2$ and $p_1\geq \dots\geq p_r\geq 2$, if $p_1\geq 2r-2$, then \[R^{PM}(p_1, \dots, p_r)= p_1- (r-1) + \sum_{i=2}^{r}\left\lceil\frac{p_i}{3}\right\rceil.\] Perhaps surprisingly, we show that when $p_1<2r-2$, it is possible that $R^{PM}(p_1, \dots, p_r)$ is larger than $p_1- (r-1) + \sum_{i=2}^{r}\left\lceil\frac{p_i}{3}\right\rceil$, but in any case we determine the correct value to within a constant (depending on $r$); i.e. \[p_1- (r-1) + \sum_{i=2}^{r}\left\lceil\frac{p_i}{3}\right\rceil \leq R^{PM}(p_1, \dots, p_r)\leq \left\lceil p_1-\frac{r}{3}+\sum_{i=2}^r\frac{p_i}{3}\right\rceil.\] As a corollary we get that in every $r$-coloring of $K_n$ there is a monochromatic path-matching of order at least $3\left\lfloor\frac{n}{r+2}\right\rfloor$, which is essentially best possible. We also determine $R^{PM}(p_1, \dots, p_r)$ in all cases when the number of colors is at most 4. The proof of the main result uses a minimax theorem for path-matchings derived from a result of Las Vergnas (extending Tutte's 1-factor theorem) to show that the value of $R^{PM}(p_1, \dots, p_r)$ depends on the block sizes in covering designs (which can be also formulated in terms of monochromatic 1-cores in colored complete graphs). Then we obtain the result above by giving estimates on the block sizes in covering designs in the arbitrary (non-uniform) case.


title: Structural and electronic phenomena at oxyfluoride KTaO$_3$/K$M$F$_3$ ($M$ = Zn and Ni) superlattices: Rashba splitting and 2DEG

abstract: Here, we present the theoretical analysis of the structural and electronic degrees of freedom of two different oxide/fluoride perovskite superlattices, KTaO$_3$/KZnF$_3$ and KTaO$_3$/KNiF$_3$. Using first-principles calculations, we found the appearance of a two-dimensional electron, 2DEG, and hole, 2DHG, gases as a function of the number of layers of the different pristine materials. We demonstrate that the phonon-dynamics at the KTO/K$M$F superlattices play a crucial role in the appearance of these effects. Additionally, our results reveal a rather sizeable Rashba-type spin-splitting at these interfaces in comparison to another oxide/oxide counterparts.


title: Spatially resolved X-ray study of supernova remnants that host magnetars: Implication of their fossil field origin

abstract: Magnetars are regarded as the most magnetized neutron stars in the Universe. Aiming to unveil what kinds of stars and supernovae can create magnetars, we have performed a state-of-the-art spatially resolved spectroscopic X-ray study of the supernova remnants (SNRs) Kes 73, RCW 103, and N49, which host magnetars 1E 1841-045, 1E 161348-5055, and SGR 0526-66, respectively. The three SNRs are O- and Ne-enhanced and are evolving in the interstellar medium with densities of >1--2 cm$^{-3}$. The metal composition and dense environment indicate that the progenitor stars are not very massive. The progenitor masses of the three magnetars are constrained to be < 20 Msun (11--15 Msun for Kes 73, < 13 Msun for RCW 103, and ~13 --17 Msun for N49). Our study suggests that magnetars are not necessarily made from very massive stars, but originate from stars that span a large mass range. The explosion energies of the three SNRs range from $10^{50}$ erg to ~2$\times 10^{51}$ erg, further refuting that the SNRs are energized by rapidly rotating (millisecond) pulsars. We report that RCW 103 is produced by a weak supernova explosion with significant fallback, as such an explosion explains the low explosion energy (~$10^{50}$ erg), small observed metal masses ($M_{\rm O}\sim 4\times 10^{-2}$ Msun and $M_{\rm Ne}\sim 6\times 10^{-3}$ Msun), and sub-solar abundances of heavier elements such as Si and S. Our study supports the fossil field origin as an important channel to produce magnetars, given the normal mass range ($M_{\rm ZAMS} < 20$ Msun) of the progenitor stars, the low-to-normal explosion energy of the SNRs, and the fact that the fraction of SNRs hosting magnetars is consistent with the magnetic OB stars with high fields.


title: Progenitors and Explosion Properties of Supernova Remnants Hosting Central Compact Objects: I. RCW 103 Associated with the Peculiar Source 1E 161348-5055

abstract: We present a Chandra and XMM-Newton imaging and spectroscopic study of the supernova remnant (SNR) RCW 103 (G332.4-00.4) containing the Central Compact Object 1E 161348-5055. The high resolution Chandra X-ray images reveal enhanced emission in the south-eastern and north-western regions. Equivalent width line images of Fe L, Mg, Si, and S using XMM-Newton data were used to map the distribution of ejecta. The SNR was sectioned into 56 regions best characterized by two-component thermal models. The harder component (kT~0.6 keV) is adequately fitted by the VPSHOCK non-equilibrium ionization model with an ionization timescale ~ 1E11-1E12 cm^-3 s, and slightly enhanced abundances over solar values. The soft component (kT~0.2 keV), fitted by the APEC model, is well described by plasma in collisional ionization equilibrium with abundances consistent with solar values. Assuming a distance of 3.1 kpc and a Sedov phase of expansion into a uniform medium, we estimate an SNR age of 4.4 kyr, a swept-up mass M_sw=16 f_s^{-1/2} D_{3.1}^{5/2} solar masses, and a low explosion energy E=3.7E49 f_s^{-1/2} D_{3.1}^{5/2} erg. This energy could be an order of magnitude higher if we relax the Sedov assumption, the plasma has a low filling factor, the plasma temperature is under-estimated, or if the SNR is expanding into the progenitor's wind-blown bubble. Standard explosion models did not match the ejecta yields. By comparing the fitted abundances to the most recent core-collapse nucleosynthesis models, our best estimate yields a low-mass progenitor around 12-13 solar masses, lower than previously reported. We discuss degeneracies in the model fitting, particularly the effect of altering the explosion energy on the progenitor mass estimate.


title: A time-dependent Hartree-Fock study of triple-alpha dynamics

abstract: Time-dependent Hartree-Fock calculations have been performed for fusion reactions of He-4 + He-4 -> Be*-8, followed by He-4 + Be*-8 . Depending on the orientation of the initial state, a linear chain vibrational state or a triangular vibration is found in 12C, with transitions between these states observed. The vibrations of the linear chain state and the triangular state occur at ~9 and 4 MeV respectively.


title: Moderate deviations for the range of a transient random walk. II

abstract: We obtain sharp upper and lower bounds for the moderate deviations of the volume of the range of a random walk in dimension five and larger. Our results encompass two regimes: a Gaussian regime for small deviations, and a stretched exponential regime for larger deviations. In the latter regime, we show that conditioned on the moderate deviations event, the walk folds a small part of its range in a ball-like subset. Also, we provide new path properties, in dimension three as well. Besides the key role Newtonian capacity plays in this study, we introduce two original ideas, of general interest, which strengthen the approach developed in \cite{AS}.


title: Plant cells-based biological matrix composites

abstract: The global increase in materials consumption calls for innovative materials, with tailored performance and multi-functionality, that are environmentally sustainable. Composites from renewable resources offer solutions to fulfil these demands but have so far been dominated by hybrid petrochemicals-based matrixes reinforced by natural fillers. Here, we present biological matrix composites with properties comparable to wood and commercial polymers. The biocomposites are obtained from cultured, undifferentiated plant cells, dehydrated and compressed under controlled conditions, forming a lamellar microstructure. Their stiffness and strength surpass that of commercial plastics of similar density, like polystyrene, and low-density polyethylene, while being entirely biodegradable. The properties can be further tuned varying the fabrication process. For example, filler particles can be integrated during fabrication, to vary the mechanical response or introduce new functionalities.


title: The spectral properties of Vandermonde matrices with clustered nodes

abstract: We study rectangular Vandermonde matrices $\mathbf{V}$ with $N+1$ rows and $s$ irregularly spaced nodes on the unit circle, in cases where some of the nodes are "clustered" together -- the elements inside each cluster being separated by at most $h \lesssim {1\over N}$, and the clusters being separated from each other by at least $\theta \gtrsim {1\over N}$. We show that any pair of column subspaces corresponding to two different clusters are nearly orthogonal: the minimal principal angle between them is at most $$\frac{\pi}{2}-\frac{c_1}{N \theta}-c_2 N h,$$ for some constants $c_1,c_2$ depending only on the multiplicities of theclusters. As a result, spectral analysis of $\mathbf{V}_N$ is significantly simplified by reducing the problem to the analysis of each cluster individually. Consequently we derive accurate estimates for 1) all the singular values of $\mathbf{V}$, and 2) componentwise condition numbers for the linear least squares problem. Importantly, these estimates are exponential only in the local cluster multiplicities, while changing at most linearly with $s$.


title: Inhomogeneous Diophantine approximation over fields of formal power series

abstract: We prove a sharp analogue of Minkowski's inhomogeneous approximation theorem over fields of power series $\mathbb{F}_q((T^{-1}))$. Furthermore, we study the approximation to a given point $\underline{y}$ in $\mathbb{F}_q((T^{-1}))^2$ by the $SL_2(\mathbb{F}_q[T])$-orbit of a given point $\underline{x}$ in $\mathbb{F}_q((T^{-1}))^2$.


title: Recognition and classification of the cosmic-ray events in images captured by CMOS/CCD cameras

abstract: Muons and other ionizing radiation produced by cosmic rays and radiative decays affect CMOS/CCD sensor. When particles colliding with sensors atoms cause specific kind of noise on images recorded by cameras. We present a concept and preliminary implementation of method for recognizing those events and algorithms for image processing and their classification by machine learning. Our method consists of analyzing the shape of traces present in images recorded by a camera sensor and metadata related to an image like camera model, GPS location of camera, vertical and horizontal orientation of a camera sensor, timestamp of image acquisition, and other events recognized near-by sensors. The so created feature vectors are classified as either a muon-like event, an electron-like event or the other event, possibly noise. For muon-like events our method estimates azimuth of a muon track. Source of the data is database of CREDO (Cosmic-Ray Extremely Distributed Observatory) project and ESO (European Southern Observatory) archives. The telescope dark frames from ESO are analysed. CREDO project collected so far over 2 millions images of events from many kinds of cameralike: smartphones camera, laptop webcams and Internet of Things cameras localised around the globe.


title: Simultaneous Estimation of Number of Clusters and Feature Sparsity in Clustering High-Dimensional Data

abstract: Estimating the number of clusters (K) is a critical and often difficult task in cluster analysis. Many methods have been proposed to estimate K, including some top performers using resampling approach. When performing cluster analysis in high-dimensional data, simultaneous clustering and feature selection is needed for improved interpretation and performance. To our knowledge, none has investigated simultaneous estimation of K and feature selection in an exploratory cluster analysis. In this paper, we propose a resampling method to meet this gap and evaluate its performance under the sparse K-means clustering framework. The proposed target function balances between sensitivity and specificity of clustering evaluation of pairwise subjects from clustering of full and subsampled data. Through extensive simulations, the method performs among the best over classical methods in estimating K in low-dimensional data. For high-dimensional simulation data, it also shows superior performance to simultaneously estimate K and feature sparsity parameter. Finally, we evaluated the methods in four microarray, two RNA-seq, one SNP and two non-omics datasets. The proposed method achieves better clustering accuracy with fewer selected predictive genes in almost all real applications.


title: Efron-Stein PAC-Bayesian Inequalities

abstract: We prove semi-empirical concentration inequalities for random variables which are given as possibly nonlinear functions of independent random variables. These inequalities describe concentration of random variable in terms of the data/distribution-dependent Efron-Stein (ES) estimate of its variance and they do not require any additional assumptions on the moments. In particular, this allows us to state semi-empirical Bernstein type inequalities for general functions of unbounded random variables, which gives user-friendly concentration bounds for cases where related methods (e.g. bounded differences) might be more challenging to apply. We extend these results to Efron-Stein PAC-Bayesian inequalities which hold for arbitrary probability kernels that define a random, data-dependent choice of the function of interest. Finally, we demonstrate a number of applications, including PAC-Bayesian generalization bounds for unbounded loss functions, empirical Bernstein type generalization bounds, new truncation-free bounds for off-policy evaluation with Weighted Importance Sampling (WIS), and off-policy PAC-Bayesian learning with WIS.


title: Current constraints from cosmogenic neutrinos on the fraction of protons in UHECRs

abstract: Cosmogenic neutrinos are created when ultra-high-energy cosmic rays (UHECRs) interact with extragalactic photon backgrounds. In general, the expected flux of these cosmogenic neutrinos depends on multiple parameters, describing the sources and propagation of UHECRs. In our recent paper (van Vliet et al. 2019), we show that a `sweet spot` occurs at a neutrino energy of $E_{\nu} \sim 1$ EeV. At that energy the flux mainly depends on two parameters, the source evolution and the fraction of protons in UHECRs at Earth for $E_p \gtrsim 30$ EeV. Therefore, with current upper limits on the cosmogenic neutrino flux at $E_{\nu} \sim 1$ EeV and assuming a certain source class, a constraint on the composition of UHECRs can be obtained. This constraint is independent of hadronic interaction models and indicates that the combination of a large proton fraction and a strong source evolution is disfavored. Upcoming neutrino experiments will be able to constrain the fraction of protons in UHECRs even further, and for any realistic model for the evolution of UHECR sources.


title: Revealing the nanoparticle composition of Edvard Munch's The Scream, and implications for paint alteration in iconic early 20th century artworks

abstract: A major motivation for the scientific study of artworks is to understand their states of preservation and ongoing degradation mechanisms. This enables preservation strategies to be developed for irreplaceable works. Intensely-hued cadmium sulphide (CdS) yellow pigments are of particular interest because these are key to the palettes of many important late 19th and early 20th century masters, including Vincent Van Gogh, Pablo Picasso, Henri Matisse, and Edvard Munch. As these paintings age, their cadmium yellow paints are undergoing severe fading, flaking, and discolouration. These effects are associated with photodegradation, the light-facilitated reactions of CdS with oxygen, moisture, and even the paint binding medium. The use of common optical and X-ray methods to characterize the physical state of the pigment is challenging due to the mixing of the various components of the paint at length scales smaller than their resolution. Here, we present an atomic-scale structural and chemical analysis of the CdS pigment in Edvard Munch's The Scream (c. 1910, Munch Museet), enabled by new electron microscope detector technologies. We show that the CdS pigment consists of clusters of defective nanoparticles ~5-10 nm in diameter. It is known from the modern use of such particles in photocatalysis that they are inherently vulnerable to photodegradation. Chlorine doping and a polytype crystal structure further enhance the sensitivity of the CdS pigment to photodegradation. In addition to The Scream, we have also observed this inherently unstable pigment structure in Henri Matisse's Flower Piece (1906, Barnes Foundation). The fundamental understanding of the pigments' nanoscale structures and impurities described here can now be used to predict which paintings are most at risk of photooxidation, and guide the most effective preservation strategies for iconic masterpieces.


title: Constraining the Gravitational-Wave Afterglow From a Binary Neutron Star Coalescence

abstract: Binary neutron star mergers are rich laboratories for physics, accessible with ground-based interferometric gravitational-wave detectors such as the Advanced LIGO and Advanced Virgo. If a neutron star remnant survives the merger, it can emit gravitational waves that might be detectable with the current or next generation detectors. The physics of the long-lived post-merger phase is not well understood and makes modelling difficult. In particular the phase of the gravitational-wave signal is not well modelled. In this paper, we explore methods for using long duration post-merger gravitational-wave signals to constrain the parameters and the properties of the remnant. We develop a phase-agnostic likelihood model that uses only the spectral content for parameter estimation and demonstrate the calculation of a Bayesian upper limit in the absence of a signal. With the millisecond magnetar model, we show that for an event like GW170817, the ellipticity of a long-lived remnant can be constrained to less than about 0.5 in the parameter space used.


title: Gaps of Summands of the Zeckendorf Lattice

abstract: A beautiful theorem of Zeckendorf states that every positive integer has a unique decomposition as a sum of non-adjacent Fibonacci numbers. Such decompositions exist more generally, and much is known about them. First, for any positive linear recurrence {Gn} the number of summands in the legal decompositions for integers in [Gn, Gn+1) converges to a Gaussian distribution. Second, Bower, Insoft, Li, Miller, and Tosteson proved that the probability of a gap between summands in a decomposition which is larger than the recurrence length converges to geometric decay. While most of the literature involves one-dimensional sequences, some recent work by Chen, Guo, Jiang, Miller, Siktar, and Yu have extended these decompositions to d-dimensional lattices, where a legal decomposition is a chain of points such that one moves in all d dimensions to get from one point to the next. They proved that some but not all properties from 1-dimensional sequences still hold. We continue this work and look at the distribution of gaps between terms of legal decompositions, and prove similar to the 1-dimensional cases that when d = 2 the gap vectors converge to a bivariate geometric random variable.


title: State Drug Policy Effectiveness: Comparative Policy Analysis of Drug Overdose Mortality

abstract: Opioid overdose rates have reached an epidemic level and state-level policy innovations have followed suit in an effort to prevent overdose deaths. State-level drug law is a set of policies that may reinforce or undermine each other, and analysts have a limited set of tools for handling the policy collinearity using statistical methods. This paper uses a machine learning method called hierarchical clustering to empirically generate "policy bundles" by grouping states with similar sets of policies in force at a given time together for analysis in a 50-state, 10-year interrupted time series regression with drug overdose deaths as the dependent variable. Policy clusters were generated from 138 binomial variables observed by state and year from the Prescription Drug Abuse Policy System. Clustering reduced the policies to a set of 10 bundles. The approach allows for ranking of the relative effect of different bundles and is a tool to recommend those most likely to succeed. This study shows that a set of policies balancing Medication Assisted Treatment, Naloxone Access, Good Samaritan Laws, Medication Assisted Treatment, Prescription Drug Monitoring Programs and legalization of medical marijuana leads to a reduced number of overdose deaths, but not until its second year in force.


title: Improved BSM Sensitivity in Diboson Processes

abstract: We study $W^+W^-$ and $Zh$ final states at future linear $e^+e^-$ colliders; designing analyses specific to the various final state polarizations allows us to target specific beyond the Standard Model (BSM) effects, parametrized in the form of dimension-6 operators. We find that CLIC can access effects roughly an order of magnitude smaller than HL-LHC or ILC, and two orders of magnitude smaller than LEP. These results are interpreted in the context of well-motivated BSM scenarios---at weak and strong coupling---where we expect correlated effects in Drell-Yann processes. The latter turn out to have better discovery potential, although the diboson processes provide additional discriminating power, potentially furnishing a way to measure the spin and coupling of BSM states.


title: The Fibonacci Quilt Game

abstract: Zeckendorf proved that every positive integer can be expressed as the sum of non-consecutive Fibonacci numbers. This theorem inspired a beautiful game, the Zeckendorf Game. Two players begin with $n \ 1$'s and take turns applying rules inspired by the Fibonacci recurrence, $F_{n+1} = F_n + F_{n-1}$, until a decomposition without consecutive terms is reached; whoever makes the last move wins. We look at a game resulting from a generalization of the Fibonacci numbers, the Fibonacci Quilt sequence. These arise from the two-dimensional geometric property of tiling the plane through the Fibonacci spiral. Beginning with 1 in the center, we place integers in the squares of the spiral such that each square contains the smallest positive integer that does not have a decomposition as the sum of previous terms that do not share a wall. This sequence eventually follows two recurrence relations, allowing us to construct a variation on the Zeckendorf Game, the Fibonacci Quilt Game. While some properties of the Fibonaccis are inherited by this sequence, the nature of its recurrence leads to others, such as Zeckendorf's theorem, no longer holding; it is thus of interest to investigate the generalization of the game in this setting to see which behaviors persist. We prove, similar to the original game, that this game also always terminates in a legal decomposition, give a lower bound on game lengths, show that depending on strategies the length of the game can vary and either player could win, and give a conjecture on the length of a random game.


title: EleAtt-RNN: Adding Attentiveness to Neurons in Recurrent Neural Networks

abstract: Recurrent neural networks (RNNs) are capable of modeling temporal dependencies of complex sequential data. In general, current available structures of RNNs tend to concentrate on controlling the contributions of current and previous information. However, the exploration of different importance levels of different elements within an input vector is always ignored. We propose a simple yet effective Element-wise-Attention Gate (EleAttG), which can be easily added to an RNN block (e.g. all RNN neurons in an RNN layer), to empower the RNN neurons to have attentiveness capability. For an RNN block, an EleAttG is used for adaptively modulating the input by assigning different levels of importance, i.e., attention, to each element/dimension of the input. We refer to an RNN block equipped with an EleAttG as an EleAtt-RNN block. Instead of modulating the input as a whole, the EleAttG modulates the input at fine granularity, i.e., element-wise, and the modulation is content adaptive. The proposed EleAttG, as an additional fundamental unit, is general and can be applied to any RNN structures, e.g., standard RNN, Long Short-Term Memory (LSTM), or Gated Recurrent Unit (GRU). We demonstrate the effectiveness of the proposed EleAtt-RNN by applying it to different tasks including the action recognition, from both skeleton-based data and RGB videos, gesture recognition, and sequential MNIST classification. Experiments show that adding attentiveness through EleAttGs to RNN blocks significantly improves the power of RNNs.