title: Automorphism groups and new constructions of maximum additive rank metric codes with restrictions

abstract: Let $d, n \in \mathbb{Z}^+$ such that $1\leq d \leq n$. A $d$-code $\mathcal{C} \subset \mathbb{F}_q^{n \times n}$ is a subset of order $n$ square matrices with the property that for all pairs of distinct elements in $\mathcal{C}$, the rank of their difference is greater than or equal to $d$. A $d$-code with as many as possible elements is called a maximum $d$-code. The integer $d$ is also called the minimum distance of the code. When $d<n$, a classical example of such an object is the so-called generalized Gabidulin code. There exist several classes of maximum $d$-codes made up respectively of symmetric, alternating and hermitian matrices. In this article we focus on such examples. Precisely, we determine their automorphism groups and solve the equivalence issue for them. Finally, we exhibit a maximum symmetric $2$-code which is not equivalent to the one with same parameters known so far.


title: Abnormality Detection in Musculoskeletal Radiographs with Convolutional Neural Networks(Ensembles) and Performance Optimization

abstract: Musculoskeletal conditions affect more than 1.7 billion people worldwide based on a study by Global Burden Disease, and they are the second greatest cause of disability[1,2]. The diagnosis of these conditions vary but mostly physical exams carried out and image tests. There are few imaging and diagnostic experts while there is a huge workload of radiograph examinations which might affect diagnostic accuracy. We built machine learning models to perform abnormality detection using the available musculoskeletal public dataset [3]. Convolutional Neural Networks (CNN) were used as are the most successful models in performing various tasks such as classification and object detection [4]. The development of the models involved theoretical study, iterative prototyping, and empirical evaluation of the results. The current model, 169 layer DenseNet, by Pranav et al.(2018) on the abnormality detection task, the performance was lower than the worst radiologist in 5 out of the 7 studies, and the overall model performance was lower than the best radiologist. We developed the ensemble200 model which scored 0.66 Cohen Kappa which was lower than the DenseNet model (Pranav et al, 2018) but the model performance with the F1 score outperforms the DenseNet model and its Cohen Kappa score variability with the different studies is lower as the best cohen kappa score on the upper extremity studies is 0.7408 (Wrist) and the lowest is (0.5844) hand. The ensemble200 model outperformed DenseNet model on the finger studies with a Cohen Kappa score of 0.653 showing reduced performance variability on the model performance.


title: Numerical Simulations of Cosmic-Ray Acceleration at Core-Collapse Supernovae

abstract: Core-collapse supernovae exploding in dense winds are favorable sites for cosmic-ray (CR) acceleration to very high energies. We present our CR-radiation-hydrodynamics simulations of the explosion of a red supergiant. We study the evolution of the shock wave during the first day following core collapse, and estimate the time at which CR acceleration can start. We then calculate the maximum CR energy at the forward shock as a function of time, and show that it may already exceed 100 TeV only a few hours after shock breakout from the surface of the star.


title: Bayesian Network Based Label Correlation Analysis For Multi-label Classifier Chain

abstract: Classifier chain (CC) is a multi-label learning approach that constructs a sequence of binary classifiers according to a label order. Each classifier in the sequence is responsible for predicting the relevance of one label. When training the classifier for a label, proceeding labels will be taken as extended features. If the extended features are highly correlated to the label, the performance will be improved, otherwise, the performance will not be influenced or even degraded. How to discover label correlation and determine the label order is critical for CC approach. This paper employs Bayesian network (BN) to model the label correlations and proposes a new BN-based CC method (BNCC). First, conditional entropy is used to describe the dependency relations among labels. Then, a BN is built up by taking nodes as labels and weights of edges as their dependency relations. A new scoring function is proposed to evaluate a BN structure, and a heuristic algorithm is introduced to optimize the BN. At last, by applying topological sorting on the nodes of the optimized BN, the label order for constructing CC model is derived. Experimental comparisons demonstrate the feasibility and effectiveness of the proposed method.


title: A note on Euler number of locally conformally K\"{a}hler manifolds

abstract: Let $M^{2n}$ be a compact Riemannian manifold of non-positive (resp. negative) sectional curvature. We call $(M,J,\theta)$ a $d$(bounded) locally conformally K\"{a}hler manifold if the lifted Lee form $\tilde{\theta}$ on the universal covering space of $M$ is $d$(bounded). We shown that if $M^{2n}$ is homeomorphic to a $d$(bounded) LCK manifold, then its Euler number satisfies the inequality $(-1)^{n}\chi(M^{2n})\geq$ (resp. $>$) $0$.


title: On the maximum number of minimal connected dominating sets in convex bipartite graphs

abstract: The enumeration of minimal connected dominating sets is known to be notoriously hard for general graphs. Currently, it is only known that the sets can be enumerated slightly faster than $\mathcal{O}^{*}(2^n)$ and the algorithm is highly nontrivial. Moreover, it seems that it is hard to use bipartiteness as a structural aide when constructing enumeration algorithms. Hence, to the best of our knowledge, there is no known input-sensitive algorithm for enumerating minimal dominating sets, or one of their related sets, in bipartite graphs better than that of general graphs. In this paper, we provide the first input-sensitive enumeration algorithm for some non trivial subclass of bipartite graphs, namely the convex graphs. We present an algorithm to enumerate all minimal connected dominating sets of convex bipartite graphs in time $\mathcal{O}(1.7254^{n})$ where $n$ is the number of vertices of the input graph. Our algorithm implies a corresponding upper bound for the number of minimal connected dominating sets for this graph class. We complement the result by providing a convex bipartite graph, which have at least $3^{(n-2)/3}$ minimal connected dominating sets.


title: Vulnerable Window of Yield Strength for Swelling-Driven Fracture of Phase-Transforming Battery Materials

abstract: Despite numerous experimental and theoretical investigations of the mechanical behavior of high-capacity Si and Ge Li-ion battery anodes, our basic understanding of swelling-driven fracture in these materials remains limited. Existing theoretical studies have provided insights into elasto-plastic deformations caused by large volume change phase transformations, but have not modeled fracture explicitly beyond Griffith's criterion. Here we use a multi-physics phase-field approach to model self-consistently anisotropic phase transformation, elasto-plastic deformation, and crack initiation and propagation during lithiation of Si nanopillars. Our results reveal the existence of a vulnerable window of yield strength inside which pillars fracture during lithiation. They identify two different modes of fracture inside that window with and without surface localization of plastic deformation prior to fracture for lower and higher yield strength, respectively, and highlight the importance of taking into account this localization to accurately predict the onset of fracture within Griffith theory. The results further demonstrate how the increased robustness of hollow nanopillars can be understood as a direct effect of anode geometry on the size of this vulnerable window. Those insights provide an improved theoretical basis for designing mechanically stable phase-transforming battery materials undergoing large volume changes.


title: Review of the AC Loss Computation for HTS using the H-formulation

abstract: This article presents a review of the finite element method (FEM) model based on the $H$ formulation of Maxwell's equations used to calculate AC losses in high temperature superconductor (HTS) tapes, cables and windings for different applications. This model, which uses the components of the magnetic field as state variables, has been gaining a great popularity and has been in use in tens of research groups around the world. This contribution first reviews the equations on which the model is based and their implementation in finite element method programs for different cases, such 2D longitudinal and axis-symmetric geometries, 3D geometries. Modeling strategies to tackle large number of HTS tapes, such as multi-scale and homogenization methods, are also introduced. Then, the second part of the article reviews the applications for which the $H$ formulations has been used to calculate AC losses, ranging from individual tapes, to complex cables and large magnet windings. Afterwards, a section is dedicated to the discussion of the $H$ formulation in terms of computational efficiency and easiness of implementation. Its pros and cons are listed. Finally, the last section draws the main conclusions.


title: Topological Nearly Entropy on Nearly Compact Spaces

abstract: In our previous paper [9], we have introduced topological nearly entropy, Ent_N (f) by restricting X into a class of nearly compact spaces. In the present paper, some additional properties of this notion are studied. Furthermore, we introduce another new notion of topological nearly entropy of f denoted by Ent_n (f) when the whole space X itself is nearly compact. We show the relationship between these two notions for the class of nearly compact subspaces. We also propose new space, namely, R-space in studying the topological nearly entropy on nearly compact and Hausdorff space. As a consequence, the topological nearly entropy of f and it restriction f|K coincides. Finally, some fundamental properties of topological nearly entropy for product space are obtained.


title: Long-term Operation of the Multi-Wire-Proportional-Chambers of the LHCb Muon System

abstract: The muon detector of LHCb, which comprises 1368 multi-wire-proportional-chambers (MWPC) for a total area of 435 m2, is the largest instrument of its kind exposed to such a high-radiation environment. In nine years of operation, from 2010 until 2018, we did not observe appreciable signs of ageing of the detector in terms of reduced performance. However, during such a long period, many chamber gas gaps suffered from HV trips. Most of the trips were due to Malter-like effects, characterised by the appearance of local self-sustained high currents, presumably originating from impurities induced during chamber production. Very effective, though long, recovery procedures were implemented with a HV training of the gaps in situ while taking data. The training allowed most of the affected chambers to be returned to their full functionality and the muon detector efficiency to be kept close to 100%. The possibility of making the recovery faster and even more effective by adding a small percentage of oxygen in the gas mixture has been studied and successfully tested.


title: Invariant-mass distribution of top-quark pairs and top-quark mass determination

abstract: We investigate the invariant-mass distribution of top-quark pairs near the $2m_t$ threshold, which has strong impact on the determination of the top-quark mass $m_t$. We show that higher-order non-relativistic corrections lead to large contributions which are not included in the state-of-the-art theoretical predictions. We derive a factorization formula to resum such corrections to all orders in the strong-coupling, and calculate necessary ingredients to perform the resummation at next-to-leading power. We combine the resummation with fixed-order results and present phenomenologically relevant numeric results. We find that the resummation effect significantly enhances the differential cross section in the threshold region, and makes the theoretical prediction more compatible with experimental data. We estimate that using our prediction in the determination of $m_t$ will lead to a value closer to the result of direct measurement.


title: A Condition Monitoring Concept Studied at the MST Prototype for the Cherenkov Telescope Array

abstract: The Cherenkov Telescope Array (CTA) is a future ground-based gamma-ray observatory that will provide unprecedented sensitivity and angular resolution for the detection of gamma rays with energies above a few tens of GeV. In comparison to existing instruments (like H.E.S.S., MAGIC, and VERITAS) the sensitivity will be improved by installing two extended arrays of telescopes in the northern and southern hemisphere, respectively. A large number of planned telescopes (>100 in total) motivates the application of predictive maintenance techniques to the individual telescopes. A constant and automatic condition monitoring of the mechanical telescope structure and of the drive system (motors, gears) is considered for this purpose. The condition monitoring system aims at detecting degradations well before critical errors occur; it should help to ensure long-term operation and to reduce the maintenance efforts of the observatory. We present approaches for the condition monitoring of the structure and the drive system of Medium-Sized Telescopes (MSTs), respectively. The overall concept has been developed and tested at the MST prototype for CTA in Berlin. The sensors used, the joint data acquisition system, possible analysis methods (like Operational Modal Analysis, OMA, and Experimental Modal Analysis, EMA) and first performance results are discussed.


title: Quantum Mechanics of Proteins in Explicit Water: The Role of Plasmon-Like Solute-Solvent Interactions

abstract: Quantum-mechanical van der Waals dispersion interactions play an essential role for both intra-protein and protein-water interactions -- the two main driving forces for the structure and dynamics of proteins in aqueous solution. Typically, these interactions are only treated phenomenologically via pairwise potential terms in classical force fields. Here, we use an explicit quantum-mechanical approach based on density-functional tight-binding with the many-body dispersion formalism, which allows us to demonstrate the unexpected relevance of the many-body character of dispersion interactions for protein energetics and the protein-water interaction. In contrast to commonly employed pairwise approaches, many-body effects significantly decrease the relative stability of the native state in the absence of water. In an aqueous environment, the collective character of the protein-water van der Waals interaction counteracts this effect and stabilizes native conformations and transition states. This stabilization arises due to a high degree of delocalization and collectivity of protein-water dispersion interactions, suggesting a remarkable persistence of long-range electron correlation through aqueous environments. Our findings are exemplified on prototypical showcases of proteins forming $\beta$-sheets, hairpins, and helices, emphasizing the crucial role of plasmon-like solute-solvent interactions in biomolecular systems.


title: An attempt at beating the 3D U-Net

abstract: The U-Net is arguably the most successful segmentation architecture in the medical domain. Here we apply a 3D U-Net to the 2019 Kidney and Kidney Tumor Segmentation Challenge and attempt to improve upon it by augmenting it with residual and pre-activation residual blocks. Cross-validation results on the training cases suggest only very minor, barely measurable improvements. Due to marginally higher dice scores, the residual 3D U-Net is chosen for test set prediction. With a Composite Dice score of 91.23 on the test set, our method outperformed all 105 competing teams and won the KiTS2019 challenge by a small margin.


title: Anisotropic electron populations in BL Lac jets: consequences for the observed emission

abstract: We investigate the impact on the properties of high-energy emitting BL Lac objects of a population of electrons with an anisotropic momentum distribution. We adopt a simple phenomenological description of the anisotropy, in which the most energetic electrons have a small pitch angle and the least energetic electrons are isotropic, as proposed by Sobacchi \& Lyubarsky (2019). We explore (i) a simple model that assumes a phenomenological shape for the electron energy distribution, and (ii) a self-consistent scheme in which the electrons follow a distribution which is the result of the balance between injection and radiative losses (we include the effects of the anisotropy on the synchrotron cooling rate). Considering the BL Lac object Mkn 421 as representative of the entire class, we show that in both cases the emission can be satisfactorily reproduced under equipartition between the magnetic field and the relativistic electrons. This is in better agreement with the idea that jets are launched as Poynting dominated flows with respect to the standard isotropic scenario, which requires both a low magnetization and a low radiative efficiency to reproduce the observed emission. The hard spectrum predicted for the inverse Compton continuum at TeV energies could be used as a potential test of the anisotropic model.


title: Growth temperature dependence of partially Fe filled MWCNT using chemical vapor deposition

abstract: This investigation deals with the effect of growth temperature on the growth behavior of Fe filled multiwalled carbon nanotubes (MWCNTs). Carbon nanotube (CNT) synthesis was carried out in a thermal chemical vapor deposition (CVD) reactor in the temperature range 650 to 950 1C using propane as the carbon source, Fe as the catalyst material, and Si as the catalyst support. Atomic force microscopy (AFM) analysis of the catalyst exhibits that at elevated temperature clusters of catalyst coalesce and form macroscopic islands. Field emission scanning electron microscopy (FESEM) results show that with increased growth temperature the average diameter of the nanotubes increases but their density decreases. High-resolution transmission electron microscopy (HRTEM) studies suggest that the nanotubes have multi-walled structure with partial Fe filling for all growth temperatures. The X-ray diffraction (XRD) pattern of the grown materials indicates that they are graphitic in nature. The characterization of nanotubes by Raman spectroscopy reveals that the optimized growth temperature for Fe filled CNTs is 850 1C, in terms of quality. A simple model for the growth of Fe filled carbon nanotubes is proposed.


title: On the initial geometry of a vacuum cosmological spacetime

abstract: In the first part of this paper we consider expanding vacuum cosmological spacetimes with a free $T^N$-action. Among them, we give evidence that Gowdy spacetimes have AVTD (asymptotically velocity term dominated) behavior for their initial geometry, in any dimension. We then give sufficient conditions to reach a similar conclusion about a $T^2$-invariant four dimensional nonGowdy spacetime. In the second part of the paper we consider vacuum cosmological spacetimes with crushing singularities. We introduce a monotonic quantity to characterize Kasner spacetimes. Assuming scale-invariant curvature bounds and local volume bounds, we give results about causal pasts.


title: $\mathbb{Z}_n$ solitons in intertwined topological phases

abstract: Topological phases of matter can support fractionalized quasi-particles localized at topological defects. The current understanding of these exotic excitations, based on the celebrated bulk-defect correspondence, typically relies on crude approximations where such defects are replaced by a static classical background coupled to the matter sector. In this work, we explore the strongly-correlated nature of symmetry-protected topological defects by focusing on situations where such defects arise spontaneously as dynamical solitons in intertwined topological phases, where symmetry breaking coexists with topological symmetry protection. In particular, we focus on the $\mathbb{Z}_2$ Bose-Hubbard model, a one-dimensional chain of interacting bosons coupled to $\mathbb{Z}_2$ fields, and show how solitons with $\mathbb{Z}_n$ topological charges appear for particle/hole dopings about certain commensurate fillings, extending the results of [1] beyond half filling. We show that these defects host fractionalized bosonic quasi-particles, forming bound states that travel through the system unless externally pinned, and repel each other giving rise to a fractional soliton lattice for sufficiently high densities. Moreover, we uncover the topological origin of these fractional bound excitations through a pumping mechanism, where the quantization of the inter-soliton transport allows us to establish a generalized bulk-defect correspondence. This in-depth analysis of dynamical topological defects bound to fractionalized quasi-particles, together with the possibility of implementing our model in cold-atomic experiments, paves the way for further exploration of exotic topological phenomena in strongly-correlated systems.

categories: cond-mat.quant-gas cond-mat.mes-hall cond-mat.str-el quant-ph


title: First-principles prediction into robust high-performance photovoltaic double perovskites A$_{2}$SiI$_{6}$ (A = K, Rb, Cs)

abstract: Despite the exceeding 23\% photovoltaic efficiency achieved in organic-inorganic hybrid perovskite solar cells obtaining, the stable materials with desirable band gap are rare and are highly desired. With the aid of first-principles calculations, we predict a new promising family of nontoxic inorganic double perovskites (DPs), namely, silicon (Si)-based halides A$_{2}$SiI$_{6}$ (A = K, Rb, Cs; X = Cl, Br, I). This family containing the earth-abundant Si could be applied for perovskite solar cells (PSCs). Particularly A$_{2}$SiI$_{6}$ exhibits superb physical traits, including suitable band gaps of 0.84-1.15 eV, dispersive lower conduction bands, small carrier effective masses, wide photon absorption in the visible range. Importantly, the good stability at high temperature renders them as promising optical absorbers for solar cells.


title: Comments on "Scattering Cancellation-Based Cloaking for the Maxwell--Cattaneo Heat Waves"

abstract: A number of errors, both mathematical and conceptual, are identified, in a recent article by Farhat \textit{et al.}\ [Phys.\ Rev.\ Appl.\ \textbf{11}, 044089 (2019)] on cloaking of thermal waves in solids, and corrected. The differences between the two thermal flux laws considered in the latter article are also critically discussed, specifically showing that the chosen model does not, in fact, correspond to the Maxwell--Cattaneo hyperbolic (wave) theory of heat transfer.


title: Superstability, noetherian rings and pure-semisimple rings

abstract: We uncover a connection between the model-theoretic notion of superstability and that of noetherian rings and pure-semisimple rings. We characterize noetherian rings via superstability of the class of left modules with embeddings. $\mathbf{Theorem.}$ For a ring $R$ the following are equivalent. - $R$ is left noetherian. - The class of left $R$-modules with embeddings is superstable. - For every $\lambda \geq |R| + \aleph_0$, there is $\chi \geq \lambda$ such that the class of left $R$-modules with embeddings has uniqueness of limit models of cardinality $\chi$. - Every limit model in the class of left $R$-modules with embeddings is $\Sigma$-injective. We characterize left pure-semisimple rings via superstability of the class of left modules with pure embeddings. $\mathbf{Theorem.}$ For a ring $R$ the following are equivalent. - $R$ is left pure-semisimple. - The class of left $R$-modules with pure embeddings is superstable. - There exists $\lambda \geq (|R| + \aleph_0)^+$ such that the class of left $R$-modules with pure embeddings has uniqueness of limit models of cardinality $\lambda$. - Every limit model in the class of left $R$-modules with pure embeddings is $\Sigma$-pure-injective. We think that both equivalences provide evidence that that the notion of superstability could shed light in the understanding of algebraic concepts. As this paper is aimed at model theorists and algebraists an effort was made to provide the background for both.


title: Simulation of temperature profile for the electron- and the lattice-systems in laterally structured layered conductors

abstract: Electrons in operating microelectronic semiconductor devices are accelerated by locally varying strong electric field to acquire effective electron temperatures nonuniformly distributing in nanoscales and largely exceeding the temperature of host crystal lattice. The thermal dynamics of electrons and the lattice are hence nontrivial and its understanding at nanoscales is decisively important for gaining higher device performance. Here, we propose and demonstrate that in layered conductors nonequilibrium nature between the electrons and the lattice can be explicitly pursued by simulating the conducting layer by separating it into two physical sheets representing, respectively, the electron- and the lattice-subsystems. We take, as an example of simulating GaAs devices, a 35nm thick 1um wide U-shaped conducting channel with 15nm radius of curvature at the inner corner of the U-shaped bend, and find a remarkable hot spot to develop due to hot electron generation at the inner corner. The hot spot in terms of the electron temperature achieves a significantly higher temperature and is of far sharper spatial distribution when compared to the hot spot in terms of the lattice temperature. Similar simulation calculation made on a metal (NiCr) narrow lead of the similar geometry shows that a hot spot shows up as well at the inner corner, but its strength and the spatial profiles are largely different from those in semiconductor devices. The remarkable difference between the semiconductor and the metal is interpreted to be due to the large difference in the electron specific heat, rather than the difference in the electron phonon interaction. This work will provide useful hints to deeper understanding of the nonequilibrium properties of electrical conductors, through a simple and convenient method for modeling nonequilibrium layered conductors.


title: On the three-dimensional stability of Poiseuille flow in a finite-length duct

abstract: The stability of a three-dimensional, incompressible, viscous flow through a finite-length duct is studied. A divergence-free basis technique is used to formulate the weak form of the problem. A SUPG (streamingline upwind Petrov-Galerkin) based scheme for eigenvalue problems is proposed to stabilize the solution. With proper boundary condtions, the least-stable eigenmodes and decay rates are computed. It is again found that the flows are asymptotically stable for all $Re$ up to $2500$. It is discovered that the least-stable eigenmodes have a boundary-layer-structure at high $Re$, although the Poiseuille base flow does not exhibits such structure. At these Reynolds numbers, the eigenmodes are dominant in the vicinity of the duct wall and are convected downstream. The boundary-layer-structure brings singularity to the modes at high $Re$ with unbounded perturbation gradient. It is shown that due to the singular structure of the least-stable eigenmodes, the linear Navier-Stoker operator tends to have pseudospectrua and the nonlinear mechanism kicks in when the perturbation energy is still small at high $Re$. The decreasing stable region as $Re$ increases is a result of both the decreasing decay rate and the singular structure of the least-stable modes. The result demonstrated that at very high $Re$, linearization of Navier-Stokes equation for duct flow may not be a good model problem with physical disturbances.


title: Special Toeplitz operators on a class of bounded Hartogs domains

abstract: We introduce a wider class of bounded Hartogs domains, which contains some generalizations of the classical Hartogs triangle. A sharp criteria for the $L^p-L^q$ boundedness of the Toeplitz operator with symbol $K^{-t}$ is obtained on these domains, where $K$ is the Bergman kernel on diagonal and $t\geq 0$. It generalizes the results by Chen and Beberok in the case $1<p<\infty$.


title: Bound on FWER for correlated normal distribution

abstract: In this paper,our main focus is to obtain an asymptotic bound on the family wise error rate (FWER) for Bonferroni-type procedure in the simultaneous hypotheses testing problem when the observations corresponding to individual hypothesis are correlated. In particular, we have considered the sequence of null hypotheses H_{0i} : X_i follows N(0,1) , (i=1,2,....,n) and equicorrelated structure of the sequence (X_1,....,X_n). Distribution free bound on FWER under equicorrelated setup can be found in Tong(2014). But the upper bound provided in Tong(2014) is not a bounded quantity as the no. of hypotheses(n) gets larger and larger and as a result,FWER is highly overestimated for the choice of a particular distribution (e.g.- normal). In the equicorrelated normal setup, we have shown that FWER asymptotically is a convex function (as a function of correlation (rho)) and hence an upper bound on the FWER of Bonferroni-(alpha) procedure is alpha(1-\rho).This implies,Bonferroni's method actually controls the FWER at a much smaller level than the desired level of significance under the positively correlated case and necessitates a correlation correction.


title: Fixed Points of Augmented Generalized Happy Functions II: Oases and Mirages

abstract: An augmented generalized happy function $S_{[c,b]}$ maps a positive integer to the sum of the squares of its base $b$ digits plus $c$. For $b\geq 2$ and $k \in \mathbb{Z}^+$, a $k$-desert base $b$ is a set of $k$ consecutive non-negative integers $c$ for each of which $S_{[c,b]}$ has no fixed points. In this paper, we examine a complementary notion, a $k$-oasis base $b$, which we define to be a set of $k$ consecutive non-negative integers $c$ for each of which $S_{[c,b]}$ has a fixed point. In particular, after proving some basic properties of oases base $b$, we compute bounds on the lengths of oases base $b$ and compute the minimal examples of maximal length oases base $b$ for small values of $b$.


title: Two invariant surface-tensors determine CSL of massive body wave function

abstract: Decoherence of massive body wave function under Continuous Spontaneous Localization is reconsidered. It is shown for homogeneous probes with wave functions narrow in position and angle that decoherence is a surface effect. Corresponding new surface integrals are derived as the main result. Probe's constant density and two completely geometric surface-dependent invariant tensors encode full dependence of positional and angular decoherence of masses, irrespective of their microscopic structure. The two surface-tensors offer a new insight into CSL and a flexible approach to design laboratory test masses.


title: Multi-Strain Age-Structured Dengue Transmission Model: Analysis and Optimal Control

abstract: Dengue is a serious health problem in the Philippines. In 2016, the Department of Health of the country launched a dengue vaccination campaign using Dengvaxia. However, the campaign was mired with controversy and the use of Dengvaxia was banned in the country. This study proposes a mathematical model that represents the dynamics of the transmission of dengue with its four strains. Considering that the Dengvaxia vaccine was intended to be given only to people aging from 9 to 45 years old, the human population is divided into two age groups: from 9-45 years old and the rest of the population. Using this model and optimal control theory we simulate what could have been the effect of Dengvaxia in the number of dengue cases and then compare this with the other usual dengue intervention strategies. Results show that the best implementation of the usual strategies is better than that of Dengvaxia.


title: Neural Blind Deconvolution Using Deep Priors

abstract: Blind deconvolution is a classical yet challenging low-level vision problem with many real-world applications. Traditional maximum a posterior (MAP) based methods rely heavily on fixed and handcrafted priors that certainly are insufficient in characterizing clean images and blur kernels, and usually adopt specially designed alternating minimization to avoid trivial solution. In contrast, existing deep motion deblurring networks learn from massive training images the mapping to clean image or blur kernel, but are limited in handling various complex and large size blur kernels. To connect MAP and deep models, we in this paper present two generative networks for respectively modeling the deep priors of clean image and blur kernel, and propose an unconstrained neural optimization solution to blind deconvolution. In particular, we adopt an asymmetric Autoencoder with skip connections for generating latent clean image, and a fully-connected network (FCN) for generating blur kernel. Moreover, the SoftMax nonlinearity is applied to the output layer of FCN to meet the non-negative and equality constraints. The process of neural optimization can be explained as a kind of "zero-shot" self-supervised learning of the generative networks, and thus our proposed method is dubbed SelfDeblur. Experimental results show that our SelfDeblur can achieve notable quantitative gains as well as more visually plausible deblurring results in comparison to state-of-the-art blind deconvolution methods on benchmark datasets and real-world blurry images. The source code is available at https://github.com/csdwren/SelfDeblur


title: HI observations of IC 10 with the DRAO synthesis telescope

abstract: HI observations of the nearby blue compact dwarf galaxy IC 10 obtained with the Dominion Radio Astrophysical Observatory synthesis telescope (DRAO), for a total integration of $\sim$1000 hours, are presented. We confirm the NW faint 21 cm HI emission feature discovered in GBT observations. The HI feature has an HI mass of 4.7 $\times 10^{5}$ M$_{\odot}$, which is only $\sim$ 0.6$\%$ of the total HI mass of the galaxy (7.8 $\times 10^{7}$ M$_{\odot}$). In the inner disk, the rotation curve of IC 10 rises steeply, then flattens until the last point where it rises again, with a maximum velocity of 30 km s$^{-1}$. Based on our mass models, the kinematics of the inner disk of IC 10 can be described without the need of a dark matter halo. However, this does not exclude the possible presence of dark matter on a larger scale. It is unlikely that the disturbed features seen in the outer HI disk of IC 10 are caused by an interaction with M 31. Features seen from our simulations are larger and at lower surface density than can be reached by current observations. The higher velocity dispersions seen in regions where several distinct HI features meet with the main core of IC 10 suggests that there is ongoing accretion.


title: MetaAdvDet: Towards Robust Detection of Evolving Adversarial Attacks

abstract: Deep neural networks (DNNs) are vulnerable to adversarial attack which is maliciously implemented by adding human-imperceptible perturbation to images and thus leads to incorrect prediction. Existing studies have proposed various methods to detect the new adversarial attacks. However, new attack methods keep evolving constantly and yield new adversarial examples to bypass the existing detectors. It needs to collect tens of thousands samples to train detectors, while the new attacks evolve much more frequently than the high-cost data collection. Thus, this situation leads the newly evolved attack samples to remain in small scales. To solve such few-shot problem with the evolving attack, we propose a meta-learning based robust detection method to detect new adversarial attacks with limited examples. Specifically, the learning consists of a double-network framework: a task-dedicated network and a master network which alternatively learn the detection capability for either seen attack or a new attack. To validate the effectiveness of our approach, we construct the benchmarks with few-shot-fashion protocols based on three conventional datasets, i.e. CIFAR-10, MNIST and Fashion-MNIST. Comprehensive experiments are conducted on them to verify the superiority of our approach with respect to the traditional adversarial attack detection methods.


title: Physics-based model of the adaptive-optics corrected point-spread-function

abstract: Context. Adaptive optics (AO) systems greatly increase the resolution of large telescopes, but produce complex point spread function (PSF) shapes, varying in time and across the field of view. This PSF must be accurately known since it provides crucial information about optical systems for design, characterisation, diagnostics and image post processing. Aims. We develop here a model of the AO long exposure PSF, adapted to various seeing conditions and any AO system. This model is made to match accurately both the core of the PSF and its turbulent halo. Methods. The PSF model we develop is based on a parsimonious parameterization of the phase power spectral density with only five parameters to describe circularly symmetric PSFs and seven parameters for asymmetrical ones. Moreover, one of the parameters is directly the Fried parameter r0 of the turbulence s strength. This physical parameter is an asset in the PSF model since it can be correlated with external measurements of the r0, such as phase slopes from the AO real time computer (RTC) or site seeing monitoring. Results. We fit our model against endtoend simulated PSFs using OOMAO tool, and against on sky PSFs from the SPHERE ZIMPOL imager and the MUSE integral field spectrometer working in AO narrowfield mode. Our model matches the shape of the AO PSF both in the core and the halo, with a sub 1 percent relative error for simulated and experimental data. We also show that we retrieve the r0 parameter with subcentimeter precision on simulated data. For ZIMPOL data, we show a correlation of 97 percent between our r0 estimation and the RTC estimation. Finally, MUSE allows us to test the spectral dependency of the fitted r0 parameter. It follows the theoretical $\lambda^{6/5}$ evolution with a standard deviation of 0.3 cm. Evolution of other PSF parameters, such as residual phase variance or aliasing, is also discussed.


title: Clouds of Fluffy Aggregates: How They Form in Exoplanetary Atmospheres and Influence Transmission Spectra

abstract: Transmission spectrum surveys have suggested the ubiquity of high-altitude clouds in exoplanetary atmospheres. Theoretical studies have investigated the formation processes of the high-altitude clouds; however, cloud particles have been commonly approximated as compact spheres, which is not always true for solid mineral particles that likely constitute exoplanetary clouds. Here, we investigate how the porosity of cloud particles evolve in exoplanetary atmospheres and influence the cloud vertical profiles. We first construct a porosity evolution model that takes into account the fractal aggregation and the compression of cloud particle aggregates. Using a cloud microphysical model coupled with the porosity model, we demonstrate that the particle internal density can significantly decrease during the cloud formation. As a result, fluffy-aggregate clouds ascend to altitude much higher than that for compact-sphere clouds assumed so far. We also examine how the fluffy-aggregate clouds affect transmission spectra. We find that the clouds largely obscure the molecular features and produce a spectral slope originated by the scattering properties of aggregates. Finally, we compare the synthetic spectra with the observations of GJ1214 b and find that its flat spectrum could be explained if the atmospheric metallicity is sufficiently high ($\ge100\times$ solar) and the monomer size is sufficiently small ($r_{\rm mon}<1~{\rm {\mu}m}$). The high-metallicity atmosphere may offer the clues to explore the gas accretion processes onto past GJ1214b.


title: Generalized Lens Categories via functors $\mathcal{C}^{\rm op}\to\mathsf{Cat}$

abstract: Lenses have a rich history and have recently received a great deal of attention from applied category theorists. We generalize the notion of lens by defining a category $\mathsf{Lens}_F$ for any category $\mathcal{C}$ and functor $F\colon \mathcal{C}^{\rm op}\to\mathsf{Cat}$, using a variant of the Grothendieck construction. All of the mathematics in this note is straightforward; the purpose is simply to see lenses in a broader context where some closely-related examples, such as ringed spaces and open continuous dynamical systems, can be included.


title: Model Agnostic Defence against Backdoor Attacks in Machine Learning

abstract: Machine Learning (ML) has automated a multitude of our day-to-day decision making domains such as education, employment and driving automation. The continued success of ML largely depends on our ability to trust the model we are using. Recently, a new class of attacks called Backdoor Attacks have been developed. These attacks undermine the user's trust in ML models. In this work, we present NEO, a model agnostic framework to detect and mitigate such backdoor attacks in image classification ML models. For a given image classification model, our approach analyses the inputs it receives and determines if the model is backdoored. In addition to this feature, we also mitigate these attacks by determining the correct predictions of the poisoned images. An appealing feature of NEO is that it can, for the first time, isolate and reconstruct the backdoor trigger. NEO is also the first defence methodology, to the best of our knowledge that is completely blackbox. We have implemented NEO and evaluated it against three state of the art poisoned models. These models include highly critical applications such as traffic sign detection (USTS) and facial detection. In our evaluation, we show that NEO can detect $\approx$88% of the poisoned inputs on average and it is as fast as 4.4 ms per input image. We also reconstruct the poisoned input for the user to effectively test their systems.


title: Cross-Origin State Inference (COSI) Attacks: Leaking Web Site States through XS-Leaks

abstract: In a Cross-Origin State Inference (COSI) attack, an attacker convinces a victim into visiting an attack web page, which leverages the cross-origin interaction features of the victim's web browser to infer the victim's state at a target web site. Multiple instances of COSI attacks have been found in the past under different names such as login detection or access detection attacks. But, those attacks only consider two states (e.g., logged in or not) and focus on a specific browser leak method (or XS-Leak). This work shows that mounting more complex COSI attacks such as deanonymizing the owner of an account, determining if the victim owns sensitive content, and determining the victim's account type often requires considering more than two states. Furthermore, robust attacks require supporting a variety of browsers since the victim's browser cannot be predicted apriori. To address these issues, we present a novel approach to identify and build complex COSI attacks that differentiate more than two states and support multiple browsers by combining multiple attack vectors, possibly using different XS-Leaks. To enable our approach, we introduce the concept of a COSI attack class. We propose two novel techniques to generalize existing COSI attack instances into COSI attack classes and to discover new COSI attack classes. We systematically apply our techniques to existing attacks, identifying 40 COSI attack classes. As part of this process, we discover a novel XS-Leak based on window.postMessage. We implement our approach into Basta-COSI, a tool to find COSI attacks in a target web site. We apply Basta-COSI to test four stand-alone web applications and 58 popular web sites, finding COSI attacks against each of them.


title: After You, Please: Browser Extensions Order Attacks and Countermeasures

abstract: Browser extensions are small applications executed in the browser context that provide additional capabilities and enrich the user experience while surfing the web. The acceptance of extensions in current browsers is unquestionable. For instance, Chrome's official extension repository has more than 63,000 extensions, with some of them having more than 10M users. When installed, extensions are pushed into an internal queue within the browser. The order in which each extension executes depends on a number of factors, including their relative installation times. In this paper, we demonstrate how this order can be exploited by an unprivileged malicious extension (i.e., one with no more permissions than those already assigned when accessing web content) to get access to any private information that other extensions have previously introduced. Our solution does not require modifying the core browser engine as it is implemented as another browser extension. We prove that our approach effectively protects the user against usual attackers (i.e., any other installed extension) as well as against strong attackers having access to the effects of all installed extensions (i.e., knowing who did what). We also prove soundness and robustness of our approach under reasonable assumptions.


title: A Case for Electron-Astrophysics

abstract: A grand-challenge problem at the forefront of physics is to understand how energy is transported and transformed in plasmas. This fundamental research priority encapsulates the conversion of plasma-flow and electromagnetic energies into particle energy, either as heat or some other form of energisation. The smallest characteristic scales, at which electron dynamics determines the plasma behaviour, are the next frontier in space and astrophysical plasma research. The analysis of astrophysical processes at these scales lies at the heart of the field of electron-astrophysics. Electron scales are the ultimate bottleneck for dissipation of plasma turbulence, which is a fundamental process not understood in the electron-kinetic regime. Since electrons are the most numerous and most mobile plasma species in fully ionised plasmas and are strongly guided by the magnetic field, their thermal properties couple very efficiently to global plasma dynamics and thermodynamics.

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


title: On quantum operations of photon subtraction and photon addition

abstract: The conventional photon subtraction and photon addition transformations, $\varrho \rightarrow t a \varrho a^{\dag}$ and $\varrho \rightarrow t a^{\dag} \varrho a$, are not valid quantum operations for any constant $t>0$ since these transformations are not trace nonincreasing. For a fixed density operator $\varrho$ there exist fair quantum operations, ${\cal N}_{-}$ and ${\cal N}_{+}$, whose conditional output states approximate the normalized outputs of former transformations with an arbitrary accuracy. However, the uniform convergence for some classes of density operators $\varrho$ has remained essentially unknown. Here we show that, in the case of photon addition operation, the uniform convergence takes place for the energy-second-moment-constrained states such that ${\rm tr}[\varrho H^2] \leq E_2 < \infty$, $H = a^{\dag}a$. In the case of photon subtraction, the uniform convergence takes place for the energy-second-moment-constrained states with nonvanishing energy, i.e., the states $\varrho$ such that ${\rm tr}[\varrho H] \geq E_1 >0$ and ${\rm tr}[\varrho H^2] \leq E_2 < \infty$. We prove that these conditions cannot be relaxed and generalize the results to the cases of multiple photon subtraction and addition.


title: Transversality of sections on elliptic surfaces with applications to elliptic divisibility sequences and geography of surfaces

abstract: We consider elliptic surfaces $\mathcal{E}$ over a field $k$ equipped with zero section $O$ and another section $P$ of infinite order. If $k$ has characteristic zero, we show there are only finitely many points where $O$ is tangent to a multiple of $P$. Equivalently, there is a finite list of integers such that if $n$ is not divisible by any of them, then $nP$ is not tangent to $O$. Such tangencies can be interpreted as unlikely intersections. If $k$ has characteristic zero or $p>3$ and $\mathcal{E}$ is very general, then we show there are no tangencies between $O$ and $nP$. We apply these results to square-freeness of elliptic divisibility sequences and to geography of surfaces. In particular, we construct mildly singular surfaces of arbitrary fixed geometric genus with $K$ ample and $K^2$ unbounded.


title: Higher-order and mixed QCD-QED corrections for Drell-Yan

abstract: In this article, we discuss about the Drell-Yan process focusing on the computation of radiative corrections for vector boson production. First, we describe the $q_T$-resummation formalism and its application to obtain higher-order QCD corrections. We briefly summarize the state-of-the-art in these calculations, and motivate the need of more refined theoretical predictions. After that, we center into the inclusion of mixed QCD-QED higher-order contributions, both for fixed order and resummed calculations. In particular, we present a framework based on the $q_T$-resummation formalism to simultaneously handle soft gluon and photon radiation. Finally, we discuss future extensions of this approach and its phenomenological relevance.


title: Knapsack Problem variants of QAOA for battery revenue optimisation

abstract: We implement two Quantum Approximate Optimisation Algorithm (QAOA) variants for a battery revenue optimisation problem, equivalent to the weakly NP-hard Knapsack Problem. Both approaches investigate how to tackle constrained problems with QAOA. A first 'constrained' approach introduces a quadratic penalty to enforce the constraint to be respected strictly and reformulates the problem into an Ising Problem. However, simulations on IBM's simulator highlight non-convergent results for intermediate depth ($ p\leq 50$). A second 'relaxed' approach applies the QAOA with a non-Ising target function to compute a linear penalty, running in time $O(p(\log_2 n)^3)$ and needing $O(n \log n)$ qubits. Simulations reveal an exponential improvement over the number of depth levels and obtain approximations about $0.95$ of the optimum with shallow depth ($p \leq 10$).


title: An Index for Sequencing Reads Based on The Colored de Bruijn Graph

abstract: In this article, we show how to transform a colored de Bruijn graph (dBG) into a practical index for processing massive sets of sequencing reads. Similar to previous works, we encode an instance of a colored dBG of the set using BOSS and a color matrix C. To reduce the space requirements, we devise an algorithm that produces a smaller and more sparse version of C. The novelties in this algorithm are (i) an incomplete coloring of the graph and (ii) a greedy coloring approach that tries to reuse the same colors for different strings when possible. We also propose two algorithms that work on top of the index; one is for reconstructing reads, and the other is for contig assembly. Experimental results show that our data structure uses about half the space of the plain representation of the set (1 Byte per DNA symbol) and that more than 99% of the reads can be reconstructed just from the index.


title: A simple decomposition of European temperature variability capturing the variance from days to a decade

abstract: We analyze European temperature variability from station data with the method of detrended fluctuation analysis. This method is known to give a scaling exponent indicating long range correlations in time for temperature anomalies. However, by a more careful look at the fluctuation function we are able to explain the emergent scaling behaviour by short time relaxation, the yearly cycle and one additional process. It turns out that for many stations this interannual variability is an oscillatory mode with a period length of approximately 7-8 years, which is consistent with results of other methods. We discuss the spatial patterns in all parameters and validate the finding of the 7-8 year period by comparing stations with and without this mode.


title: A Universality Theorem for Nested Polytopes

abstract: In a nutshell, we show that polynomials and nested polytopes are topological, algebraic and algorithmically equivalent. Given two polytops $A\subseteq B$ and a number $k$, the Nested Polytope Problem (NPP) asks, if there exists a polytope $X$ on $k$ vertices such that $A\subseteq X \subseteq B$. The polytope $A$ is given by a set of vertices and the polytope $B$ is given by the defining hyperplanes. We show a universality theorem for NPP. Given an instance $I$ of the NPP, we define the solutions set of $I$ as $$ V'(I) = \{(x_1,\ldots,x_k)\in \mathbb{R}^{k\cdot n} : A\subseteq \text{conv}(x_1,\ldots,x_k) \subseteq B\}.$$ As there are many symmetries, induced by permutations of the vertices, we will consider the \emph{normalized} solution space $V(I)$. Let $F$ be a finite set of polynomials, with bounded solution space. Then there is an instance $I$ of the NPP, which has a rationally-equivalent normalized solution space $V(I)$. Two sets $V$ and $W$ are rationally equivalent if there exists a homeomorphism $f : V \rightarrow W$ such that both $f$ and $f^{-1}$ are given by rational functions. A function $f:V\rightarrow W$ is a homeomorphism, if it is continuous, invertible and its inverse is continuous as well. As a corollary, we show that NPP is $\exists \mathbb{R}$-complete. This implies that unless $\exists \mathbb{R} =$ NP, the NPP is not contained in the complexity class NP. Note that those results already follow from a recent paper by Shitov. Our proof is geometric and arguably easier.


title: A new functional RG flow: regulator-sourced 2PI versus average 1PI

abstract: We derive the renormalization group evolution of the quartic scalar theory with spontaneous symmetry breaking from an alternative flow equation, obtained within the externally sourced two-particle irreducible framework due to Garbrecht and Millington. In order to make a straightforward comparison with the evolution from the standard Wetterich-Morris-Ellwanger equation, we employ the Litim regulator, work to lowest order in the derivative expansion and neglect anomalous scaling. By this means, we illustrate the leading differences between analytic expressions for the resulting threshold and (non-perturbative) beta functions. In four dimensions, we find that the positions of the potential minima and the cosmological constant evolve more rapidly with scale compared to the standard approach, whereas the quartic coupling evolves more slowly, albeit by a small amount. These differences may have implications for the asymptotic safety programme, as well as our understanding of the non-perturbative scale evolution of the Standard Model Higgs sector.


title: On cylindrical regression in three-dimensional Euclidean space

abstract: The three-dimensional cylindrical regression problem is a problem of finding a cylinder best fitting a group of points in three-dimensional Euclidean space. The words best fitting are usually understood in the sense of the minimum root mean square deflection of the given points from a cylinder to be found. In this form the problem has no analytic solution. If one replaces the root mean square averaging by a certain biquadratic averaging, the resulting problem has an almost analytic solution. This solution is reproduced in the present paper in a coordinate-free form.


title: Global well-posedness for the 2-D inhomogeneous incompressible Navier-Stokes system with large initial data in critical spaces

abstract: Without any smallness assumption, we prove the global unique solvability of the 2-D incompressible inhomogeneous Navier-Stokes equations with initial data in the critical Besov space, which is almost the energy space in the sense that they have the same scaling in terms of this 2-D system.


title: Biases in retrieving planetary signals in the presence of quasi-periodic stellar activity

abstract: Gaussian process regression is a widespread tool used to mitigate stellar correlated noise in radial velocity time series. It is particularly useful to search for and determine the properties of signals induced by small-size, low-mass planets ($R_p<4R_{\rm \oplus}$, $m_p<10M_{\rm \oplus}$). By using extensive simulations based on a quasi-periodic representation of the stellar activity component, we investigate the ability in retrieving the planetary parameters in 16 different realistic scenarios. We analyse systems composed by one planet and host stars having different levels of activity, focusing on the challenging case represented by low-mass planets, with Doppler semi-amplitudes in the range 1-3 $m s^{-1}$. We consider many different configurations for the quasi-periodic stellar activity component, as well as different combinations of the observing epochs. We use commonly-employed analysis tools to search for and characterize the planetary signals in the datasets. The goal of our injection-recovery statistical analysis is twofold. First, we focus on the problem of planet mass determination. Then, we analyse in a statistical way periodograms obtained with three different algorithms, in order to explore some of their general properties, as the completeness and reliability in retrieving the injected planetary and stellar activity signals with low false alarm probabilities. This work is intended to provide some understanding of the biases introduced in the planet parameters inferred from the analysis of radial velocity time series that contain correlated signals due to stellar activity. It also aims to motivate the use and encourage the improvement of extensive simulations for planning spectroscopic follow-up observations.


title: Should we test the model assumptions before running a model-based test?

abstract: Statistical methods are based on model assumptions, and it is statistical folklore that a method's model assumptions should be checked before applying it. This can be formally done by running one or more misspecification tests of model assumptions before running a method that requires these assumptions; here we focus on model-based tests. A combined test procedure can be defined by specifying a protocol in which first model assumptions are tested and then, conditionally on the outcome, a test is run that requires or does not require the tested assumptions. Although such an approach is often taken in practice, much of the literature that investigated this is surprisingly critical of it. Our aim is to explore conditions under which model checking is advisable or not advisable. For this, we review results regarding such "combined procedures" in the literature, we review and discuss controversial views on the role of model checking in statistics, and we present a general setup in which we can show that preliminary model checking is advantageous, which implies conditions for making model checking worthwhile.


title: Gyrokinetic Simulations on Many- and Multi-core Architectures with the Global Electromagnetic Particle-In-Cell Code ORB5

abstract: Gyrokinetic codes in plasma physics need outstanding computational resources to solve increasingly complex problems, requiring the effective exploitation of cutting-edge HPC architectures. This paper focuses on the enabling of ORB5, a state-of-the-art, first-principles-based gyrokinetic code, on modern parallel hybrid multi-core, multi-GPU systems. ORB5 is a Lagrangian, Particle-In-Cell (PIC), finite element, global, electromagnetic code, originally implementing distributed parallelism through MPI-based on domain decomposition and domain cloning. In order to support multi/many cores devices, the code has been completely refactored. Data structures have been re-designed to ensure efficient memory access, enhancing data locality. Multi-threading has been introduced through OpenMP on the CPU and adopting OpenACC to support GPU acceleration. MPI can still be used in combination with the two approaches. The performance results obtained using the full production ORB5 code on the Summit system at ORNL, on Piz Daint at CSCS and on the Marconi system at CINECA are presented, showing the effectiveness and performance portability of the adopted solutions: the same source code version was used to produce all results on all architectures.


title: Constructing cospectral signed graphs

abstract: A well--known fact in Spectral Graph Theory is the existence of pairs of isospectral nonisomorphic graphs (known as PINGS). The work of A.J. Schwenk (in 1973) and of C. Godsil and B. McKay (in 1982) shed some light on the explanation of the presence of isospectral graphs, and they gave routines to construct PINGS. Here, we consider the Godsil-McKay--type routines developed for graphs, whose adjacency matrices are $(0,1)$-matrices, to the level of signed graphs, whose adjacency matrices allow the presence of $-1$'s. We show that, with suitable adaption, such routines can be successfully ported to signed graphs, and we can build pairs of cospectral switching nonisomorphic signed graphs.


title: Preliminary design of a device to assist handwriting in children with movement disorders

abstract: This paper presents the development of a new passive assistive handwriting device, which aims to stabilize the motion of people living with movement disorders. Many people living with conditions such as cerebral palsy, stroke, muscular dystrophy or dystonia experience upper limbs impairments (muscle spasticity, unselective motor control, muscle weakness or tremors) and are unable to write or draw on their own. The proposed device is designed to be fixed on a table. A pen is attached to the device using a pen holder, which maintains the pen in a fixed orientation. The user interacts with the device using a handle while mechanical dampers and inertia contribute to the stabilization of the user's movements. The overall mechanical design of the device is first presented, followed by the design of the pen holder mechanism.


title: Lowest-degree triple Massey products in moment-angle complexes

abstract: We give a combinatorial classification of non-trivial triple Massey products of three dimensional classes in the cohomology of a moment-angle complex. This work improves on a result by Denham and Suciu (2007) by considering triple Massey products with non-trivial indeterminacy.


title: Neutron valence structure from nuclear deep inelastic scattering

abstract: Mechanisms of spin-flavor SU(6) symmetry breaking in Quantum Chromodynamics (QCD) are studied via an extraction of the free neutron structure function from a global analysis of deep inelastic scattering (DIS) data on the proton and on nuclei from $A = 2$ (deuterium) to 208 (lead). Modification of the structure function of nucleons bound in atomic nuclei (known as the EMC effect) are consistently accounted for within the framework of a universal modification of nucleons in short-range correlated (SRC) pairs. Our extracted neutron-to-proton structure function ratio $F_2^n/F_2^p$ becomes constant for $x_B \ge 0.6$, equalling $0.47 \pm 0.04$ as $x_B \rightarrow 1$, in agreement with theoretical predictions of perturbative QCD and the Dyson Schwinger equation, and in disagreement with predictions of the Scalar Diquark dominance model. We also predict $F_2^{^3\mathrm{He}}/F_2^{^3\mathrm{H}}$, recently measured, yet unpublished, by the MARATHON collaboration, the nuclear correction function that is needed to extract $F_2^n/F_2^p$ from $F_2^{^3\mathrm{He}}/F_2^{^3\mathrm{H}}$, and the theoretical uncertainty associated with this extraction.


title: Brachistochrone on a Velodrome

abstract: The Brachistochrone problem, which describes the curve that carries a particle under gravity in a vertical plane from one height to another in the shortest time, is one of the most famous studies in classical physics. There is a similar problem in track cycling, where a cyclist aims to find the trajectory on the curved sloping surface of a velodrome that results in the minimum lap time. In this paper we extend the classical Brachistochrone problem to find the optimum cycling trajectory in a velodrome, treating the cyclist as an active particle. Starting with two canonical cases of cycling on a sloping plane and a cone, where analytical solutions are found, we then solve the problem numerically on the reconstructed surface of the velodrome in Montigny le Bretonneux, France. Finally, we discuss the parameters of the problem and the effects of fatigue.


title: Designing an All-Carbon Membrane for Water Desalination

abstract: We design an all-carbon membrane for the filtration and desalination of water. A unique layered assembly of carbon nanostructures including graphite oxide (GO), buckypaper consisting of carbon nanotubes, and a strong carbon fabric provides high mechanical strength and thermal stability, resilience to harsh chemical cleaning agents and electrical conductivity, thus addressing major shortcomings of commercial reverse osmosis membranes. We use ab initio density functional theory calculations to obtain atomic-level insight into the permeation of water molecules in-between GO layers and across in-layer vacancy defects. Our calculations elucidate the reason for selective rejection of solvated Na$^+$ ions in an optimized GO membrane that is structurally stabilized in a sandwich arrangement in-between layers of buckypaper, which are protected on both sides by strong carbon fabric layers.


title: Differences in the solar cycle variability of simple and complex active regions during 1996-2018

abstract: Aims. Our aim is to examine the solar cycle variability of magnetically simple and complex active region. Methods. We studied simple ($\alpha$ and $\beta$) and complex ($\beta\gamma$ and $\beta\gamma\delta$) active regions based on the Mount Wilson magnetic classification by applying our newly developed daily approach. We analyzed the daily number of the simple active regions (SARs) and compared that to the abundance of the complex active regions (CARs) over the entire solar cycle 23 and cycle 24 until December 2018. Results. We show that CARs evolve differently over the solar cycle from SARs. The time evolution of SARs and CARs on different hemispheres also shows differences, even though on average their latitudinal distributions are shown to be similar. The time evolution of SARs closely follows that of the sunspot number, and their maximum abundance was observed to occur during the early maximum phase, while that of the CARs was seen roughly two years later. We furthermore found that the peak of CARs was reached before the latitudinal width of the activity band starts to decease. Conclusions. Our results suggest that the active region formation process is a competition between the large-scale dynamo (LSD) and the small-scale dynamo (SSD) near the surface, the former varying cyclically and the latter being independent of the solar cycle. During solar maximum, LSD is dominant, giving a preference to SARs, while during the declining phase the relative role of SSD increases. Therefore, a preference for CARs is seen due to the influence of the SSD on the emerging flux.


title: Conservative Link Adaptation for Ultra Reliable Low Latency Communications

abstract: Ultra reliable low latency communications (URLLC) is one of the most promising and demanding services in 5G systems. This service requires very low latency of less than $1-10$ ms and very high transmission reliability: the acceptable packet loss ratio is about $10^{-5}$. To satisfy such strict requirements, many issues shall be solved. This paper focuses on the link adaptation problem, i.e., the selection of a modulation and coding scheme (MCS) for transmission based on the received channel quality indicator (CQI) reports. On the one hand, link adaptation should select a robust MCS to provide high reliability. On the other hand, it should select the highest possible MCS to reduce channel resource consumption. The paper shows that even for one URLLC user, link adaptation is still a challenging problem, especially in highly-variant channels. To solve this problem, a conservative link adaptation algorithm is designed. The algorithm estimates the strongest channel degradation at the time moment of the actual packet transmission and selects an MCS taking into account the worst degradation. The obtained results show that the proposed algorithm is efficient in terms of both the packet loss ratio and the channel resource consumption.


title: Risk-Control Strategies

abstract: In this paper, we consider the pricing of derivative products that involve dynamic hedging strategies and payments within the planning horizon. Equity-indexed annuities (EIAs), Guaranteed investment certificate (GIC), American and Barrier options are typical examples of these products. Our exploration involves evaluation under different assumptions related to the way the risk is tailored by the issuer. The unified constrained discrete stochastic dynamic programming framework presented in this paper makes use of sequential local minimizing strategies related to stochastic transitions. This sequential minimizations takes into account all intermediate requirements and involves several dynamic risk measures modelling. To demonstrate the flexibility of this framework we present numerical examples featuring GICs and point-to-point EIAs.


title: Detection of multiple stellar populations in extragalactic massive clusters with JWST

abstract: The discovery both through spectroscopy and photometry of multiple stellar populations in Galactic globular clusters, and in Magellanic Clouds' massive intermediate-age and old clusters, has led to a major change in our views about the formation of these objects. To date, none of the proposed scenarios is able to explain quantitatively all chemical patterns observed in individual clusters, and an extension of the study of multiple populations to resolved extragalactic massive clusters beyond the Magellanic Clouds would be welcome, for it would enable to investigate and characterize the presence of multiple populations in different environments and age ranges. To this purpose, the James Webb Space Telescope can potentially play a major role. On the one hand, the JWST promises direct observations of proto-globular cluster candidates at high redshift; on the other hand, it can potentially push to larger distances the sample of resolved clusters with detected multiple populations. In this paper we have addressed this second goal. Using theoretical stellar spectra and stellar evolution models, we have investigated the effect of multiple population chemical patterns on synthetic magnitudes in the JWST infrared NIRCam filters. We have identified the colours (F150W-F460M), (F115W-F460M) and pseudocolours C_{F150W,F460M,F115W}=(F150W-F460M)-(F460M-F115W), C_{F150W,F277W,F115W}=(F150W-F277W)-(F277W-F115W), as diagnostics able to reveal the presence of multiple populations along the red giant branches of old and intermediate age clusters. Using the available on-line simulator for the NIRCam detector, we have estimated that multiple populations can be potentially detected --depending on the exposure times, exact filter combination used, plus the extent of the abundance variations and the cluster [Fe/H]--out to a distance of about 5Mpc (approximately the distance to the M83 group).


title: Outer linear measure of connected sets via Steiner trees

abstract: We resurrect an old definition of the linear measure of a metric continuum in terms of Steiner trees, independently due to Menger (1930) and Choquet (1938). We generalise it to any metric space and provide a proof of a little-known theorem of Choquet that it coincides with the outer linear measure for any connected metric space. As corollaries we obtain simple proofs of Go{\l}\k{a}b's theorem (1928) on the lower semicontinuity of linear measure of continua and a theorem of Bogn\'ar (1989) on the linear measure of the closure of a set. We do not use any measure theory apart from the definition of outer linear measure.


title: Learning Aberrance Repressed Correlation Filters for Real-Time UAV Tracking

abstract: Traditional framework of discriminative correlation filters (DCF) is often subject to undesired boundary effects. Several approaches to enlarge search regions have been already proposed in the past years to make up for this shortcoming. However, with excessive background information, more background noises are also introduced and the discriminative filter is prone to learn from the ambiance rather than the object. This situation, along with appearance changes of objects caused by full/partial occlusion, illumination variation, and other reasons has made it more likely to have aberrances in the detection process, which could substantially degrade the credibility of its result. Therefore, in this work, a novel approach to repress the aberrances happening during the detection process is proposed, i.e., aberrance repressed correlation filter (ARCF). By enforcing restriction to the rate of alteration in response maps generated in the detection phase, the ARCF tracker can evidently suppress aberrances and is thus more robust and accurate to track objects. Considerable experiments are conducted on different UAV datasets to perform object tracking from an aerial view, i.e., UAV123, UAVDT, and DTB70, with 243 challenging image sequences containing over 90K frames to verify the performance of the ARCF tracker and it has proven itself to have outperformed other 20 state-of-the-art trackers based on DCF and deep-based frameworks with sufficient speed for real-time applications.


title: Spectral properties of spin-orbital polarons as a fingerprint of orbital order

abstract: Transition metal oxides are a rich group of materials with very interesting physical properties that arise from the interplay of the charge, spin, orbital, and lattice degrees of freedom. One interesting consequence of this, encountered in systems with orbital degeneracy, is the coexistence of long range magnetic and orbital order, and the coupling between them. In this paper we develop and study an effective spin-orbital superexchange model for $e_g^3$ systems and use it to investigate the spectral properties of a charge (hole) injected into the system, which is relevant for photoemission spectroscopy. Using an accurate, semi-analytical, magnon expansion method, we gain insight into various physical aspects of these systems and demonstrate a number of subtle effects, such as orbital to magnetic polaron crossover, the coupling between orbital and magnetic order, as well as the orbital order driving the system towards one-dimensional quantum spin liquid behavior. Our calculations also suggest a potentially simple experimental verification of the character of the orbital order in the system, something that is not easily accessible through most experimental techniques.


title: Koopman Representations of Dynamic Systems with Control

abstract: The design and analysis of optimal control policies for dynamical systems can be complicated by nonlinear dependence in the state variables. Koopman operators have been used to simplify the analysis of dynamical systems by mapping the flow of the system onto a space of observables where the dynamics are linear (and possibly infinte). This paper focuses on the development of consistent Koopman representations for controlled dynamical system. We introduce the concept of dynamical consistency for Koopman representations and analyze several existing and proposed representations deriving necessary constraints on the dynamical system, observables, and Koopman operators. Our main result is a hybrid formulation which independently and jointly observes the state and control inputs. This formulation admits a relatively large space of dynamical systems compared to earlier formulations while keeping the Koopman operator independent of the state and control inputs. More generally, this work provides an analysis framework to evaluate and rank proposed simplifications to the general Koopman representation for controlled dynamical systems.


title: The Variance of the Number of Zeros for Complex Random Polynomials Spanned by OPUC

abstract: Let $\{\varphi_k\}_{k=0}^\infty $ be a sequence of orthonormal polynomials on the unit circle (OPUC) with respect to a probability measure $ \mu $. We study the variance of the number of zeros of random linear combinations of the form $$ P_n(z)=\sum_{k=0}^{n}\eta_k\varphi_k(z), $$ where $\{\eta_k\}_{k=0}^n $ are complex-valued random variables. Under the assumption that the distribution for each $\eta_k$ satisfies certain uniform bounds for the fractional and logarithmic moments, for the cases when $\{\varphi_k\}$ are regular in the sense of Ullman-Stahl-Totik or are such that the measure of orthogonality $\mu$ satisfies $d\mu(\theta)=w(\theta)d\theta$ where $w(\theta)=v(\theta)\prod_{j=1}^J|\theta - \theta_j|^{\alpha_j}$, with $v(\theta)\geq c>0$, $\theta,\theta_j\in [0,2\pi)$, and $\alpha_j>0$, we give a quantitative estimate on the the variance of the number of zeros of $P_n$ in sectors that intersect the unit circle. When $\{\varphi_k\}$ are real-valued on the real-line from the Nevai class and $\{\eta_k\}$ are i.i.d.~complex-valued standard Gaussian, we prove a formula for the limiting value of variance of the number of zeros of $P_n$ in annuli that do not contain the unit circle.


title: Real Clifford algebras and their spinors for relativistic fermions

abstract: Real Clifford algebras for arbitrary number of space and time dimensions as well as their representations in terms of spinors are reviewed and discussed. The Clifford algebras are classified in terms of isomorphic matrix algebras of real, complex or quaternionic type. Spinors are defined as elements of minimal or quasi-minimal left ideals within the Clifford algebra and as representations of the pin and spin groups. Two types of Dirac adjoint spinors are introduced carefully. The relation between mathematical structures and applications to describe relativistic fermions is emphasized throughout.


title: A New Model for Electron-Capture Supernovae in Galactic Chemical Evolution

abstract: We examine the contribution of electron-capture supernovae (ECSNe), low-mass SNe from collapsing Fe cores (FeCCSNe), and rotating massive stars to the chemical composition of the Galaxy. Our model includes contributions to chemical evolution from both thermonuclear ECSNe (tECSNe) and gravitational collapse ECSNe (cECSNe). We show that if ECSNe are predominantly gravitational collapse SNe but about 15% are partial thermonuclear explosions, the model is able to reproduce the solar abundances of several important and problematic isotopes including $^{48}$Ca, $^{50}$Ti and $^{54}$Cr together with $^{58}$Fe, $^{64}$Ni, $^{82}$Se and $^{86}$Kr and several of the Zn--Zr isotopes. A model in which no cECSNe occur, only tECSNe with low-mass FeCCSNe or rotating massive stars, proves also very successful at reproducing the solar abundances for these isotopes. Despite the small mass range for the progenitors of ECSNe and low-mass FeCCSNe, the large production factors suffice for the solar inventory of the above isotopes. Our model is compelling because it introduces no new tensions with the solar abundance distribution for a Milky Way model -- only tending to improve the model predictions for several isotopes. The proposed astrophysical production model thus provides a natural and elegant way to explain one of the last uncharted territories on the periodic table of astrophysical element production.


title: Dynamical Nucleon-Pion System via Basis Light-Front Quantization

abstract: We present the first application of the Basis Light-Front Quantization method to study a simple chiral model of the nucleon-pion system via an ab initio, non-perturbative, Hamiltonian approach. As a test problem, we consider the physical proton as the relativistic bound state of the nucleon-pion system. Based on the chiral model of the nucleon-pion system, we construct the mass-squared matrix of the system within our light-front basis representation. We obtain the proton's mass and the corresponding light-front wave function by solving the eigenvalue problem of the mass-squared matrix. With the resulting boost-invariant light-front wave function, we also compute the proton's parton distribution function.


title: A Survey of Phase Classification Techniques for Characterizing Variable Application Behavior

abstract: Adaptable computing is an increasingly important paradigm that specializes system resources to variable application requirements, environmental conditions, or user requirements. Adapting computing resources to variable application requirements (or application phases) is otherwise known as phase-based optimization. Phase-based optimization takes advantage of application phases, or execution intervals of an application, that behave similarly, to enable effective and beneficial adaptability. In order for phase-based optimization to be effective, the phases must first be classified to determine when application phases begin and end, and ensure that system resources are accurately specialized. In this paper, we present a survey of phase classification techniques that have been proposed to exploit the advantages of adaptable computing through phase-based optimization. We focus on recent techniques and classify these techniques with respect to several factors in order to highlight their similarities and differences. We divide the techniques by their major defining characteristics---online/offline and serial/parallel. In addition, we discuss other characteristics such as prediction and detection techniques, the characteristics used for prediction, interval type, etc. We also identify gaps in the state-of-the-art and discuss future research directions to enable and fully exploit the benefits of adaptable computing.


title: Tuning Algorithms and Generators for Efficient Edge Inference

abstract: A surge in artificial intelligence and autonomous technologies have increased the demand toward enhanced edge-processing capabilities. Computational complexity and size of state-of-the-art Deep Neural Networks (DNNs) are rising exponentially with diverse network models and larger datasets. This growth limits the performance scaling and energy-efficiency of both distributed and embedded inference platforms. Embedded designs at the edge are constrained by energy and speed limitations of available processor substrates and processor to memory communication required to fetch the model coefficients. While many hardware accelerator and network deployment frameworks have been in development, a framework is needed to allow the variety of existing architectures, and those in development, to be expressed in critical parts of the flow that perform various optimization steps. Moreover, premature architecture-blind network selection and optimization diminish the effectiveness of schedule optimizations and hardware-specific mappings. In this paper, we address these issues by creating a cross-layer software-hardware design framework that encompasses network training and model compression that is aware of and tuned to the underlying hardware architecture. This approach leverages the available degrees of DNN structure and sparsity to create a converged network that can be partitioned and efficiently scheduled on the target hardware platform, minimizing data movement, and improving the overall throughput and energy. To further streamline the design, we leverage the high-level, flexible SoC generator platform based on RISC-V ROCC framework. This integration allows seamless extensions of the RISC-V instruction set and Chisel-based rapid generator design. Utilizing this approach, we implemented a silicon prototype in a 16 nm TSMC process node achieving record processing efficiency of up to 18 TOPS/W.


title: Biologically inspired sleep algorithm for artificial neural networks

abstract: Sleep plays an important role in incremental learning and consolidation of memories in biological systems. Motivated by the processes that are known to be involved in sleep generation in biological networks, we developed an algorithm that implements a sleep-like phase in artificial neural networks (ANNs). After initial training phase, we convert the ANN to a spiking neural network (SNN) and simulate an offline sleep-like phase using spike-timing dependent plasticity rules to modify synaptic weights. The SNN is then converted back to the ANN and evaluated or trained on new inputs. We demonstrate several performance improvements after applying this processing to ANNs trained on MNIST, CUB200 and a motivating toy dataset. First, in an incremental learning framework, sleep is able to recover older tasks that were otherwise forgotten in the ANN without sleep phase due to catastrophic forgetting. Second, sleep results in forward transfer learning of unseen tasks. Finally, sleep improves generalization ability of the ANNs to classify images with various types of noise. We provide a theoretical basis for the beneficial role of the brain-inspired sleep-like phase for the ANNs and present an algorithmic way for future implementations of the various features of sleep in deep learning ANNs. Overall, these results suggest that biological sleep can help mitigate a number of problems ANNs suffer from, such as poor generalization and catastrophic forgetting for incremental learning.


title: Bessel SPDEs with general Dirichlet boundary conditions

abstract: We generalise the integration by parts formulae obtained in arXiv:1811.00518v5 [math.PR] to Bessel bridges on $[0,1]$ with arbitrary boundary values, as well as Bessel processes with arbitrary initial conditions. This allows us to write, formally, the corresponding dynamics using renormalised local times, thus extending the Bessel SPDEs of arXiv:1811.00518v5 [math.PR] to general Dirichlet boundary conditions. We also prove a dynamical result for the case of dimension $2$, by providing a weak construction of the gradient dynamics corresponding to a $2$-dimensional Bessel bridge.


title: Toward quantitative fractography using convolutional neural networks

abstract: The science of fractography revolves around the correlation between topographic characteristics of the fracture surface and the mechanisms and external conditions leading to their creation. While being a topic of investigation for centuries, it has remained mostly qualitative to date. A quantitative analysis of fracture surfaces is of prime interest for both the scientific community and the industrial sector, bearing the potential for improved understanding on the mechanisms controlling the fracture process and at the same time assessing the reliability of computational models currently being used for material design. With new advances in the field of image analysis, and specifically with machine learning tools becoming more accessible and reliable, it is now feasible to automate the process of extracting meaningful information from fracture surface images. Here, we propose a method of identifying and quantifying the relative appearance of intergranular and transgranular fracture events from scanning electron microscope images. The newly proposed method is based on a convolutional neural network algorithm for semantic segmentation. The proposed method is extensively tested and evaluated against two ceramic material systems ($Al_2O_3$,$MgAl_2O_4$) and shows high prediction accuracy, despite being trained on only one material system ($MgAl_2O_4$). While here attention is focused on brittle fracture characteristics, the method can be easily extended to account for other fracture morphologies, such as dimples, fatigue striations, etc.


title: Design, Modeling, and Control of Norma: a Slider & Pendulum-Driven Spherical Robot

abstract: This paper discusses the design, modeling, and control of Norma, a novel 2 DOF mobile spherical robot (SR). The propelling mechanism of this robot consists of two actuators: a slider, and a rotational pendulum located on the SR's diagonal shaft. The slider can translate along the shaft and shift the robot's center of gravity towards the robot's sides. The pendulum rotates around the shaft to propel the SR to roll forward and backward. These two actuators enable the SR to perform both rolling and turning maneuvers as a nonholonomic robot. The advantage of the proposed mechanical design lies in its convenience of physical implementation, agility, and accurate mathematical model. The Euler Lagrange approach is utilized to derive the dynamics of the proposed mechanical structure using minimum simplifications possible. Further, a path tracking control scheme is introduced for a smooth trajectory. Finally, simulations are carried out in MATLAB to verify the accuracy of the mathematical model and the effectiveness of the controller against experimental results.


title: Parallel mathematical models of dynamic objects

abstract: The paper deals with the developing of the methodological backgrounds for the modeling and simulation of complex dynamical objects. Such backgrounds allow us to perform coordinate transformation and formulate the algorithm of its usage for transforming the serial mathematical model into parallel ones. This algorithm is based on partial fraction decomposition of the transfer function of a dynamic object. Usage of proposed algorithms is one of the ways to decrease calculation time and improve PC usage while a simulation is being performed. We prove our approach by considering the example of modeling and simulating of fourth order dynamical object with various eigenvalues. This example shows that developed parallel model is stable, well-convergent, and high-accuracy model. There is no defined any calculation errors between well-known serial model and proposed parallel one. Nevertheless, the proposed approach's usage allows us to reduce calculation time by more than 20% by using several CPU's cores while calculations are being performed.


title: A construction of support $\tau$-tilting modules over $\tau$-tilting finite algebras

abstract: The notion of (semi)bricks, regarded as a generalization of (semi)simple modules, appeared in a paper of Ringel in 1976. In recent years, there has been several new developments motivated by links to {\tau}-tilting theory studied by Demonet-Iyama-Jasso and Asai. In this paper, we will discuss how to glue semibricks along a recollement with the intermediate extension functor similar to gluing simple modules by Beilinson-Bernstein-Deligne. As an application, we investigate the behavior of {\tau}-tilting finite under recollements of module categories of algebras. Moreover, we give some examples to show the construction of support {\tau}-tilting modules over {\tau}-tilting finite algebras by gluing semibricks via recollements.


title: On Convergence of Distributed Approximate Newton Methods: Globalization, Sharper Bounds and Beyond

abstract: The DANE algorithm is an approximate Newton method popularly used for communication-efficient distributed machine learning. Reasons for the interest in DANE include scalability and versatility. Convergence of DANE, however, can be tricky; its appealing convergence rate is only rigorous for quadratic objective, and for more general convex functions the known results are no stronger than those of the classic first-order methods. To remedy these drawbacks, we propose in this paper some new alternatives of DANE which are more suitable for analysis. We first introduce a simple variant of DANE equipped with backtracking line search, for which global asymptotic convergence and sharper local non-asymptotic convergence rate guarantees can be proved for both quadratic and non-quadratic strongly convex functions. Then we propose a heavy-ball method to accelerate the convergence of DANE, showing that nearly tight local rate of convergence can be established for strongly convex functions, and with proper modification of algorithm the same result applies globally to linear prediction models. Numerical evidence is provided to confirm the theoretical and practical advantages of our methods.


title: Constraints on the binary black hole nature of GW151226 and GW170608 from the measurement of spin-induced quadrupole moments

abstract: According to the "no-hair" conjecture, a Kerr black hole (BH) is completely described by its mass and spin. In particular, the spin-induced quadrupole moment of a Kerr BH with mass $m$ and dimensionless spin $\chi$ can be written as $Q=-\kappa\,m^3\chi^2$, where $\kappa_{\rm BH}=1$. Thus by measuring the spin-induced quadrupole parameter $\kappa$, we can test the binary black hole nature of compact binaries and distinguish them from binaries comprised of other exotic compact objects, as proposed in [N. V. Krishnendu et al., PRL 119, 091101 (2017)]. Here, we present a Bayesian framework to carry out this test where we measure the symmetric combination of individual spin-induced quadrupole moment parameters fixing the anti-symmetric combination to be zero. The analysis is restricted to the inspiral part of the signal as the spin-induced deformations are not modeled in the post-inspiral regime. We perform detailed simulations to investigate the applicability of this method for compact binaries of different masses and spins and also explore various degeneracies in the parameter space which can affect this test. We then apply this method to the gravitational wave events, GW151226 and GW170608 detected during the first and second observing runs of Advanced LIGO and Advanced Virgo detectors. We find the two events to be consistent with binary black hole mergers in general relativity. By combining information from several more of such events in future, this method can be used to set constraints on the black hole nature of the population of compact binaries that are detected by the Advanced LIGO and Advanced Virgo detectors.


title: Provable bounds for the Korteweg-de Vries reduction in multi-component Nonlinear Schrodinger Equation

abstract: We study the dynamics of multi-component Bose gas described by the Vector Nonlinear Schr\"{o}dinger Equation (VNLS), aka the Vector Gross--Pitaevskii Equation (VGPE) . Through a Madelung transformation, the VNLS can be reduced to coupled hydrodynamic equations in terms of multiple density and velocity fields. Using a multi-scaling and a perturbation method along with the Fredholm alternative, we reduce the problem to a Korteweg de-Vries (KdV) system. This is of great importance to study more transparently, the obscure features hidden in VNLS. This ensures that hydrodynamic effects such as dispersion and nonlinearity are captured at an equal footing. Importantly, before studying the KdV connection, we provide a rigorous analysis of the linear problem. We write down a set of theorems along with proofs and associated corollaries that shine light on the conditions of existence and nature of eigenvalues and eigenvectors of the linear problem. This rigorous analysis is paramount for understanding the nonlinear problem and the KdV connection. We provide strong evidence of agreement between VNLS systems and KdV equations by using soliton solutions as a platform for comparison. Our results are expected to be relevant not only for cold atomic gases, but also for nonlinear optics and other branches where VNLS equations play a defining role.

categories: math-ph cond-mat.quant-gas math.MP nlin.SI physics.optics


title: Reduced-order control using low-rank Dynamic Mode Decomposition

abstract: In this work we perform full-state LQR feedback control of fluid flows using non-intrusive data-driven reduced-order models. We propose a model reduction method called low-rank Dynamic Mode Decomposition (lrDMD) that solves for a rank-constrained linear representation of the dynamical system. lrDMD is shown to have lower data reconstruction error compared to standard Optimal Mode Decomposition (OMD) and Dynamic Mode Decomposition (DMD), but with an increased computational cost arising from solving a non-convex matrix optimization problem. We demonstrate model order reduction on the complex linearized Ginzburg-Landau equation in the globally unstable regime and on the unsteady flow over a flat plate at a high angle of attack. In both cases, low-dimensional full-state feedback controller is constructed using reduced-order models constructed using DMD, OMD and lrDMD. It is shown that lrDMD stabilizes the Ginzburg-Landau system with a lower order controller and is able to suppress vortex shedding from an inclined flat plate at a cost lower than either DMD or OMD. It is further shown that lrDMD yields an improved estimate of the adjoint system, for a given rank, relative to DMD and OMD.


title: Properties of the cumulated deficient binary digit sum

abstract: The sequence A268289 from the On-Line Encyclopedia of Integer Sequences, namely the cumulated differences between the number of digits 1 and the number of digits 0 in the binary expansion of consecutive integers, is studied here. This sequence happens to match a sequence of cardinalities of some specific sets. Furthermore, it can also be expressed by using the Takagi function. The three different definitions have their own properties and combining them together lead to some new identities.


title: Grid dissections of tangential quadrilaterals

abstract: For any integer $n\ge 2$, a square can be partitioned into $n^2$ smaller squares via a checkerboard-type dissection. Does there such a class-preserving grid dissection exist for some other types of quadrilaterals? For instance, is it true that a tangential quadrilateral can be partitioned into $n^2$ smaller tangential quadrilaterals using an $n\times n$ grid dissection? We prove that the answer is affirmative for every integer $n\ge 2$.