title: The quantum query complexity of composition with a relation

abstract: The negative weight adversary method, $\mathrm{ADV}^\pm(g)$, is known to characterize the bounded-error quantum query complexity of any Boolean function $g$, and also obeys a perfect composition theorem $\mathrm{ADV}^\pm(f \circ g^n) = \mathrm{ADV}^\pm(f) \mathrm{ADV}^\pm(g)$. Belovs gave a modified version of the negative weight adversary method, $\mathrm{ADV}_{rel}^\pm(f)$, that characterizes the bounded-error quantum query complexity of a relation $f \subseteq \{0,1\}^n \times [K]$, provided the relation is efficiently verifiable. A relation is efficiently verifiable if $\mathrm{ADV}^\pm(f_a) = o(\mathrm{ADV}_{rel}^\pm(f))$ for every $a \in [K]$, where $f_a$ is the Boolean function defined as $f_a(x) = 1$ if and only if $(x,a) \in f$. In this note we show a perfect composition theorem for the composition of a relation $f$ with a Boolean function $g$ \[ \mathrm{ADV}_{rel}^\pm(f \circ g^n) = \mathrm{ADV}_{rel}^\pm(f) \mathrm{ADV}^\pm(g) \enspace . \] For an efficiently verifiable relation $f$ this means $Q(f \circ g^n) = \Theta( \mathrm{ADV}_{rel}^\pm(f) \mathrm{ADV}^\pm(g) )$.


title: Analysis of Maxwell-Stefan systems for heat conducting fluid mixtures

abstract: The global-in-time existence of bounded weak solutions to the Maxwell-Stefan-Fourier equations in Fick-Onsager form is proved. The model consists of the mass balance equations for the partial mass densities and and the energy balance equation for the total energy. The diffusion and heat fluxes depend linearly on the gradients of the thermo-chemical potentials and the gradient of the temperature and include the Soret and Dufour effects. The cross-diffusion system exhibits an entropy structure, which originates from the thermodynamic modeling. The lack of positive definiteness of the diffusion matrix is compensated by the fact that the total mass density is constant in time. The entropy estimate yields the a.e. positivity of the partial mass densities and temperature. Also diffusion matrices are considered that degenerate for vanishing partial mass densities.


title: Chemotaxis and Reactions in Biology

abstract: Chemotaxis plays a crucial role in a variety of processes in biology and ecology. Quite often it acts to improve efficiency of biological reactions. One example is the immune system signalling, where infected tissues release chemokines attracting monocytes to fight invading bacteria. Another example is reproduction, where eggs release pheromones that attract sperm. A macro scale example is flower scent appealing to pollinators. In this paper we consider a system of PDE designed to model such processes. Our interest is to quantify the effect of chemotaxis on reaction rates compared to pure reaction-diffusion. We limit consideration to surface chemotaxis, which is well motivated from the point of view of many applications. Our results provide the first insight into situations where chemotaxis can be crucial for reaction success, and where its effect is likely to be limited. The proofs are based on new analytical tools; a significant part of the paper is dedicated to building up the linear machinery that can be useful in more general settings. In particular we establish precise estimates on the rates of convergence to ground state for a class of Fokker-Planck operators with potentials that grow at a logarithmic rate at infinity. These estimates are made possible by a new sharp weak weighted Poincar\'e inequality improving in particular a result of Bobkov and Ledoux.


title: An equivalence criterion for infinite products of Cauchy measures

abstract: We give an equivalence-singularity criterion for infinite products of Cauchy measures under simultaneous shifts of the location and scale parameters. Our result is an extension of Lie and Sullivan's result giving an equivalence-singularity criterion under dilations of scale parameters. Our proof utilizes McCullagh's parameterization of the Cauchy distributions and maximal invariant, and a closed-form formula of the Kullback-Leibler divergence between two Cauchy measures given by Chyzak and Nielsen.


title: Particle-based Energetic Variational Inference

abstract: We introduce a new variational inference (VI) framework, called energetic variational inference (EVI). It minimizes the VI objective function based on a prescribed energy-dissipation law. Using the EVI framework, we can derive many existing Particle-based Variational Inference (ParVI) methods, including the popular Stein Variational Gradient Descent (SVGD) approach. More importantly, many new ParVI schemes can be created under this framework. For illustration, we propose a new particle-based EVI scheme, which performs the particle-based approximation of the density first and then uses the approximated density in the variational procedure, or "Approximation-then-Variation" for short. Thanks to this order of approximation and variation, the new scheme can maintain the variational structure at the particle level, and can significantly decrease the KL-divergence in each iteration. Numerical experiments show the proposed method outperforms some existing ParVI methods in terms of fidelity to the target distribution.


title: m-subharmonic and m-plurisubharmonic functions -- on two problems of Sadullaev

abstract: We show that the spaces of $A$-$m$-subharmonic and $B$-$m$-subharmonic functions differ in sufficently high dimensions. We also prove that the Monge-Amp\`ere type operator $\mathcal M_m$ associated to the space of $m$-plurisubharmonic functions does not allow an integral comparison principle except in the classical cases $m=1$ and $m=n$. These answer in the negative two problems posed by A. Sadullaev.


title: Stochastic modeling of non-linear adsorption with Gaussian kernel density estimators

abstract: Adsorption is a relevant process in many fields, such as product manufacturing or pollution remediation in porous materials. Adsorption takes place at the molecular scale, amenable to be modeled by Lagrangian numerical methods. We have proposed a chemical diffusion-reaction model for the simulation of adsorption, based on the combination of a random walk particle tracking method involving the use of Gaussian Kernel Density Estimators. The main feature of the proposed model is that it can effectively reproduce the nonlinear behavior characteristic of the Langmuir and Freundlich isotherms. In the former, it is enough to add a finite number of sorption sites of homogeneous sorption properties, and to set the process as the combination of the forward and the backward reactions, each one of them with a prespecified reaction rate. To model the Freundlich isotherm instead, typical of low to intermediate range of solute concentrations, there is a need to assign a different equilibrium constant to each specific sorption site, provided they are all drawn from a truncated power-law distribution. Both nonlinear models can be combined in a single framework to obtain a typical observed behavior for a wide range of concentration values.


title: Network positions in active learning environments in physics

abstract: This study uses positional analysis to describe the student interaction networks in four research-based introductory physics curricula. Positional analysis is a technique for simplifying the structure of a network into blocks of actors whose connections are more similar to each other than to the rest of the network. This method describes social structure in a way that is comparable between networks of different sizes and densities and can show large-scale patterns such as hierarchy or brokering among actors. We detail the method and apply it to class sections using Peer Instruction, SCALE-UP, ISLE, and context-rich problems. At the level of detail shown in the blockmodels, most of the curricula are more alike than different, showing a late-term tendency to form coherent subgroups that communicate actively among themselves but have few inter-position links. This pattern may be a network signature of active learning classes, but wider data collection is needed to investigate.


title: Molecular docking and binding mode analysis of selected FDA approved drugs against COVID-19 selected key protein targets: An effort towards drug repurposing to identify the combination therapy to combat COVID-19

abstract: The emergence of COVID-19 has severely compromised the arsenal of antiviral and antibiotic drugs. Drug discovery is a multistep process with a high failure rate, high cost and it takes approximately 10-12 years for the development of new molecules into the clinical candidate. On the other side, drug repurposing also called old drugs for new uses, is an attractive alternative approach for a new application of marketed FDA approved or investigational drugs. In the current pandemic situation raised due to COVID-19, repurposing of existing FDA approved drugs are emerging as the first line of the treatment. The causative viral agent of this highly contagious disease and acute respiratory syndrome coronavirus (SARS-CoV) shares high nucleotide similarity. Therefore, many existing viral targets are structurally expected to be similar to SARS-CoV and likely to be inhibited by the same compounds. Here, we selected three viral key proteins based on their vital role in viral life cycle: ACE2 (helps in entry into the human host), viral nonstructural proteins RNA-dependent RNA polymerase (RdRp) NSP12, and NSP16 which helps in replication, and viral latency (invasion from immunity). The FDA approved drugs chloroquine (CQ), hydroxychloroquine (HCQ), remdesivir (RDV) and arbidol (ABD) are emerging as promising agents to combat COVID-19. Our hypothesis behind the docking studies is to determine the binding affinities of these drugs and identify the key amino acid residues playing a key role in their mechanism of action. The docking studies were carried out through Autodock and online COVID-19 docking server. Further studies on a broad range of FDA approved drugs including few more protein targets, molecular dynamics studies, in-vitro and in-vivo biological evaluation are required to identify the combination therapy targeting various stages of the viral life cycle.


title: Measurement Error in Nutritional Epidemiology: A Survey

abstract: This article reviews bias-correction models for measurement error of exposure variables in the field of nutritional epidemiology. Measurement error usually attenuates estimated slope towards zero. Due to the influence of measurement error, inference of parameter estimate is conservative and confidence interval of the slope parameter is too narrow. Bias-correction in estimators and confidence intervals are of primary interest. We review the following bias-correction models: regression calibration methods, likelihood based models, missing data models, simulation based methods, nonparametric models and sampling based procedures.


title: H$_2$ content of galaxies inside and around intermediate redshift clusters

abstract: Dense environments have an impact on the star formation rate of galaxies. As stars form from molecular gas, looking at the cold molecular gas content of a galaxy gives useful insights on its efficiency in forming stars. However, most galaxies observed in CO (a proxy for the cold molecular gas content) at intermediate redshifts, are field galaxies. Only a handful of studies focused on cluster galaxies. I present new results on the environment of one medium mass cluster from the EDisCS survey at $z\sim0.5$. 27 star-forming galaxies were selected to evenly sample the range of densities encountered inside and around the cluster. We cover a region extending as far as 8 virial radii from the cluster center. Indeed there is ample evidence that star formation quenching starts already beyond 3 cluster virial radius. I discuss our CO(3-2) ALMA observations, which unveil a large fraction of galaxies with low gas-to-stellar mass ratios.


title: Role of a Habitat's Air Humidity in Covid-19 Mortality

abstract: Transient local over-dry environment might be a contributor and an explanation for the observed asynchronous local rises in Covid-19 mortality. We propose that a habitat's air humidity negatively correlate with Covid-19 morbidity and mortality, and support this hypothesis on the example of publicly available data from German federal states.


title: Vacuum polarization contribution to muon $g-2$ as an inverse problem

abstract: We analyze the electromagnetic current correlator at an arbitrary photon invariant mass $q^2$ by exploiting its associated dispersion relation. The dispersion relation is turned into an inverse problem, via which the involved vacuum polarization function $\Pi(q^2)$ at low $q^2$ is solved with the perturbative input of $\Pi(q^2)$ at large $q^2$. It is found that the result for $\Pi(q^2)$, including its first derivative $\Pi^\prime(q^2=0)$, agrees with those from lattice QCD, and its imaginary part accommodates the $e^+e^-$ annihilation data. The corresponding hadronic vacuum polarization contribution $a^{\rm HVP}_\mu= (641^{+65}_{-63})\times 10^{-10}$ to the muon anomalous magnetic moment $g-2$, where the uncertainty arises from the variation of the perturbative input, also agrees with those obtained in other phenomenological and theoretical approaches. We point out that our formalism is equivalent to imposing the analyticity constraint to the phenomenological approach solely relying on experimental data, and can improve the precision of the $a^{\rm HVP}_\mu$ determination in the Standard Model.


title: On the Incidence Rate of Blazhko Stars

abstract: In a recent paper (Kovacs 2018) we examined the incidence rate of the modulated RR Lyrae stars by using the data from Campaigns 01 - 04 of the Kepler K2 mission. We found an observed rate of ~90%, implying near 100% underlying rate, after correcting for detection bias due to observational noise. In this work we extend the sample to Campaign 08 and check the reliability of our earlier estimate. We get the same high rate, stressing the importance of full time series modelling (including systematics) in searching for shallow signals in the presence of large amplitude variabilities.


title: Peer Offloading in Mobile Edge Computing with Worst-Case Response Time Guarantees

abstract: Mobile edge computing (MEC) is a new paradigm that provides cloud computing services at the edge of networks. To achieve better performance with limited computing resources, peer offloading between cooperative edge servers (e.g. MEC- enabled base stations) has been proposed as an effective technique to handle bursty and spatially imbalanced arrival of computation tasks. While various performance metrics of peer offloading policies have been considered in the literatures, the worst-case response time, a common Quality of Service(QoS) requirement in real-time applications, yet receives much less attention. To fill the gap, we formulate the peer offloading problem based on a stochastic arrival model and propose two online algorithms for cases with and without prior knowledge of task arrival rate. Our goal is to maximize the utility function of time-average throughput under constraints of energy consumption and worst-case response time. Both theoretical analysis and numerical results show that our algorithms are able to produce close to optimal performance.


title: On Remote Phonon Scattering

abstract: Polar phonons can induce electric fields in an adjacent layer, whether non-polar or polar, producing remote phonon scattering of electrons. Treatment of remote phonon scattering has been based on the dielectric continuum model which takes only the electrical boundary conditions into account. We show that crystals whose polar modes satisfy both mechanical and electric boundary conditions cannot introduce remote phonon effects in the absence of dispersion. Further, even in the presence of dispersion, remote phonon effects are negligible, as a consequence of the necessity of satisfying mechanical boundary conditions.


title: Tensor Network Rewriting Strategies for Satisfiability and Counting

abstract: We provide a graphical treatment of SAT and #SAT on equal footing. Instances of #SAT can be represented as tensor networks in a standard way. These tensor networks are interpreted by diagrams of the ZH-calculus: a system to reason about tensors over C in terms of diagrams built from simple generators, in which computation may be carried out by transformations of diagrams alone. In general, nodes of ZH diagrams take parameters over C which determine the tensor coefficients; for the standard representation of #SAT instances, the coefficients take the value 0 or 1. Then, by choosing the coefficients of a diagram to range over B, we represent the corresponding instance of SAT. Thus, by interpreting a diagram either over the boolean semiring or the complex numbers, we instantiate either the decision or counting version of the problem. We find that for classes known to be in P, such as 2SAT and #XORSAT, the existence of appropriate rewrite rules allows for efficient simplification of the diagram, producing the solution in polynomial time. In contrast, for classes known to be NP-complete, such as 3SAT, or #P-complete, such as #2SAT, the corresponding rewrite rules introduce hyperedges to the diagrams, in numbers which are not easily bounded above by a polynomial. This diagrammatic approach unifies the diagnosis of the complexity of CSPs and #CSPs and shows promise in aiding tensor network contraction-based algorithms.


title: On conformal field theories based on Takiff superalgebras

abstract: We revisit the construction of conformal field theories based on Takiff algebras and superalgebras that was introduced by Babichenko and Ridout. Takiff superalgebras can be thought of as truncated current superalgebras with Z-grading which arise from taking p copies of a Lie superalgebra g and placing them in the degrees s=0,...,p-1. Using suitably defined non-degenerate invariant forms we show that Takiff superalgebras give rise to families of conformal field theories with central charge c=p sdim(g). The resulting conformal field theories are defined in the standard way, i.e. they lend themselves to a Lagrangian description in terms of a WZW model and their chiral energy momentum tensor is the one obtained naturally from the usual Sugawara construction. In view of their intricate representation theory they provide interesting examples of conformal field theories.


title: A null test of the equivalence principle using relativistic effects in galaxy surveys

abstract: The weak equivalence principle is one of the cornerstone of general relativity. Its validity has been tested with impressive precision in the Solar System, with experiments involving baryonic matter and light. However, on cosmological scales and when dark matter is concerned, the validity of this principle is still unknown. In this paper we construct a null test that probes the validity of the equivalence principle for dark matter. Our test has the strong advantage that it can be applied on data without relying on any modelling of the theory of gravity. It involves a combination of redshift-space distortions and relativistic effects in the galaxy number-count fluctuation, that vanishes if and only if the equivalence principle holds. We show that the null test is very insensitive to typical uncertainties in other cosmological parameters, including the magnification bias parameter, and to non-linear effects, making this a robust null test for modified gravity.


title: General Lieb-Schultz-Mattis type theorems for quantum spin chains

abstract: We develop a general operator algebraic method which focuses on projective representations of symmetry group for proving Lieb-Schultz-Mattis type theorems, i.e., no-go theorems that rule out the existence of a unique gapped ground state (or, more generally, a pure split state), for quantum spin chains with on-site symmetry. We first prove a theorem for translation invariant spin chains that unifies and extends two theorems proved by two of the authors in [OT1]. We then prove a Lieb-Schultz-Mattis type theorem for spin chains that are invariant under the reflection about the origin and not necessarily translation invariant.


title: The R Package stagedtrees for Structural Learning of Stratified Staged Trees

abstract: stagedtrees is an R package which includes several algorithms for learning the structure of staged trees and chain event graphs from data. Score-based and clustering-based algorithms are implemented, as well as various functionalities to plot the models and perform inference. The capabilities of stagedtrees are illustrated using mainly two datasets both included in the package or bundled in R.


title: A generalized one phase Stefan problem as a vanishing viscosity limit

abstract: We study the vanishing viscosity limit of a nonlinear diffusion equation describing chemical reaction interface or the spatial segregation interface of competing species, where the diffusion rate for the negative part of the solution converges to zero. As in the standard one phase Stefan problem, we prove that the positive part of the solution converges uniformly to the solution of a generalized one phase Stefan problem. This information is then employed to determine the limiting equation for the negative part, which is an ordinary differential equation.


title: Small-time asymptotics of hypoelliptic heat kernels near the diagonal, nilpotentization and related results

abstract: We establish small-time asymptotic expansions for heat kernels of hypoelliptic H\"ormander operators in a neighborhood of the diagonal, generalizing former results obtained in particular by M\'etivier and by Ben Arous. The coefficients of our expansions are identified in terms of the nilpotentization of the underlying sub-Riemannian structure. Our approach is purely analytic and relies in particular on local and global subelliptic estimates as well as on the local nature of small-time asymptotics of heat kernels. The fact that our expansions are valid not only along the diagonal but in an asymptotic neighborhood of the diagonal is the main novelty, useful in view of deriving Weyl laws for subelliptic Laplacians. In turn, we establish a number of other results on hypoelliptic heat kernels that are interesting in themselves, such as Kac's principle of not feeling the boundary, asymptotic results for singular perturbations of hypoelliptic operators, global smoothing properties for selfadjoint heat semigroups.


title: On dual transform of fractional Hankel transform

abstract: We deal with a class of one-parameter family of integral transforms of Bargmann type arising as dual transforms of fractional Hankel transform. Their ranges are identified to be special subspaces of the weighted hyperholomorphic left Hilbert spaces, generalizing the slice Bergman space of the second kind. Their reproducing kernel is given by closed expression involving the $\star$-regularization of Gauss hypergeometric function. We also discuss their basic properties such as their boundedness and we determinate their singular values. Moreover, we describe their compactness and membership in $p$-Schatten classes.


title: Anomalous Motion of Particle Levitated by Laguerre-Gaussian beam

abstract: Laguerre-Gaussian (LG) beam has orbital angular momentum (OAM). A particle trapped in an LG beam will rotate about the beam axis, due to the transfer of OAM. The rotation of the particle is usually in the same direction as that of the beam OAM. However, we discovere that when the LG beam is strongly focused, the rotation of the particle and the beam OAM might be in the opposite direction. This anomalous effect is caused by the negative torque on the particle exerted by the focused LG beam, which is similar to the optical pulling force in the linear case. We calculated the scattering force distribution of a micro-particle trapped in an optical tweezers formed by the strongly focused LG beam, and showed that there exist stable trajectories of the particle that controlled by the negative torque. We proposed several necessary conditions for observing the counter-intuitive trajectories. Our work reveals that the strongly trapped micro-particle exhibits diversity of motion patterns.


title: Winning by hiding behind others: An analysis of speed skating data

abstract: In some athletic races, such as cycling and types of speed skating races, athletes have to complete a relatively long distance at a high speed in the presence of direct opponents. To win such a race, athletes are motivated to hide behind others to suppress energy consumption before a final moment of the race. This situation seems to produce a social dilemma: players want to hide behind others, whereas if a group of players attempts to do so, they may all lose to other players that overtake them. To support that speed skaters are involved in such a social dilemma, we analyzed video footage data for 14 mass start skating races to find that skaters that hid behind others to avoid air resistance for a long time before the final lap tended to win. Furthermore, the finish rank of the skaters in mass start races was independent of the record of the same skaters in time-trial races measured in the absence of direct opponents. The results suggest that how to strategically cope with a skater's dilemma may be a key determinant for winning long-distance and high-speed races with direct opponents.


title: Deep Learning Models for Multilingual Hate Speech Detection

abstract: Hate speech detection is a challenging problem with most of the datasets available in only one language: English. In this paper, we conduct a large scale analysis of multilingual hate speech in 9 languages from 16 different sources. We observe that in low resource setting, simple models such as LASER embedding with logistic regression performs the best, while in high resource setting BERT based models perform better. In case of zero-shot classification, languages such as Italian and Portuguese achieve good results. Our proposed framework could be used as an efficient solution for low-resource languages. These models could also act as good baselines for future multilingual hate speech detection tasks. We have made our code and experimental settings public for other researchers at https://github.com/punyajoy/DE-LIMIT.


title: Coupled Multi-Proca Vector Dark Energy

abstract: We study a new class of vector dark energy models where multi-Proca fields $A_\mu^a$ are coupled to cold dark matter by the term $f(X)\tilde{\mathcal{L}}_{m}$ where $f(X)$ is a general function of $X\equiv -\frac{1}{2}A^\mu_ a A^a_\mu$ and $\tilde{\mathcal{L}}_{m}$ is the cold dark matter Lagrangian. From here, we derive the general covariant form of the novel interaction term sourcing the field equations. This result is quite general in the sense that encompasses Abelian and non-Abelian vector fields. In particular, we investigate the effects of this type of coupling in a simple dark energy model based on three copies of canonical Maxwell fields to realize isotropic expansion. The cosmological background dynamics of the model is examined by means of a dynamical system analysis to determine the stability of the emergent cosmological solutions. As an interesting result, we find that the coupling function leads to the existence of a novel scaling solution during the dark matter dominance. Furthermore, the critical points show an early contribution of the vector field in the form of dark radiation and a stable de Sitter-type attractor at late times mimicking dark energy. The cosmological evolution of the system as well as the aforementioned features are verified by numerical computations. Observational constraints are also discussed to put the model in a more phenomenological context in the light of future observations.


title: On Symmetrized Pearson's Type Test for Normality of Autoregression: Power under Local Alternatives

abstract: We consider a stationary linear AR($p$) model with observations subject to gross errors (outliers). The autoregression parameters as well as the distribution function (d.f.) $G$ of innovations are unknown. The distribution of outliers $\Pi$ is unknown and arbitrary, their intensity is $\gamma n^{-1/2}$ with an unknown $\gamma$, $n$ is the sample size. We test the hypothesis for normality of innovations $$\mathbf{H}_\Phi \colon G \in \{\Phi(x/\theta),\,\theta>0\},$$ $\Phi(x)$ is the d.f. $\mathbf{N}(0,1)$. Our test is the special symmetrized Pearson's type test. We find the power of this test under local alternatives $$\mathbf{H}_{1n}(\rho)\colon G(x)=A_n(x):=(1-\rho n^{-1/2})\Phi(x/\theta_0)+\rho n^{-1/2}H(x), $$ $\rho\geq 0,\,\theta_0$ is the unknown (under $\mathbf{H}_\Phi$) variance of innovations. First of all we estimate the autoregression parameters and then using the residuals from the estimated autoregression we construct a kind of empirical distribution function (r.e.d.f.), which is a counterpart of the (inaccessible) e.d.f. of the autoregression innovations. After this we construct the symmetrized variant r.e.d.f. Our test statistic is the functional from symmetrized r.e.d.f. We obtain a stochastic expansion of this symmetrized r.e.d.f. under $\mathbf{H}_{1n}(\rho)$ , which enables us to investigate our test. We establish qualitative robustness of this test in terms of uniform equicontinuity of the limiting power (as functions of $\gamma,\rho$ and $\Pi$) with respect to $\gamma$ in a neighborhood of $\gamma=0$.


title: Self6D: Self-Supervised Monocular 6D Object Pose Estimation

abstract: 6D object pose estimation is a fundamental problem in computer vision. Convolutional Neural Networks (CNNs) have recently proven to be capable of predicting reliable 6D pose estimates even from monocular images. Nonetheless, CNNs are identified as being extremely data-driven, and acquiring adequate annotations is oftentimes very time-consuming and labor intensive. To overcome this shortcoming, we propose the idea of monocular 6D pose estimation by means of self-supervised learning, removing the need for real annotations. After training our proposed network fully supervised with synthetic RGB data, we leverage recent advances in neural rendering to further self-supervise the model on unannotated real RGB-D data, seeking for a visually and geometrically optimal alignment. Extensive evaluations demonstrate that our proposed self-supervision is able to significantly enhance the model's original performance, outperforming all other methods relying on synthetic data or employing elaborate techniques from the domain adaptation realm.


title: Efficient Approximation Algorithms for Adaptive Influence Maximization

abstract: Given a social network $G$ and an integer $k$, the influence maximization (IM) problem asks for a seed set $S$ of $k$ nodes from $G$ to maximize the expected number of nodes influenced via a propagation model. The majority of the existing algorithms for the IM problem are developed only under the non-adaptive setting, i.e., where all $k$ seed nodes are selected in one batch without observing how they influence other users in real world. In this paper, we study the adaptive IM problem where the $k$ seed nodes are selected in batches of equal size $b$, such that the $i$-th batch is identified after the actual influence results of the former $i-1$ batches are observed. In this paper, we propose the first practical algorithm for the adaptive IM problem that could provide the worst-case approximation guarantee of $1-\mathrm{e}^{\rho_b(\varepsilon-1)}$, where $\rho_b=1-(1-1/b)^b$ and $\varepsilon \in (0, 1)$ is a user-specified parameter. In particular, we propose a general framework AdaptGreedy that could be instantiated by any existing non-adaptive IM algorithms with expected approximation guarantee. Our approach is based on a novel randomized policy that is applicable to the general adaptive stochastic maximization problem, which may be of independent interest. In addition, we propose a novel non-adaptive IM algorithm called EPIC which not only provides strong expected approximation guarantee, but also presents superior performance compared with the existing IM algorithms. Meanwhile, we clarify some existing misunderstandings in recent work and shed light on further study of the adaptive IM problem. We conduct experiments on real social networks to evaluate our proposed algorithms comprehensively, and the experimental results strongly corroborate the superiorities and effectiveness of our approach.


title: Habitability is a continuous property of nature

abstract: In their recent comment, Cockell et al. argue that the habitability of an environment is fundamentally a binary property; that is to say, an environment can either support the metabolic processes of a given organism or not. The habitability of different environments, they argue, may have different degrees that could be determined at least in theory by answering the question: 'is this environment habitable to a given organism?' 'More' or 'less' habitable environments could then be related by the number of yes or no answers given to what is fundamentally a series of binary questions and decisions. In my opinion, there are at least three implicit assumptions made for this line of reasoning that are implausible and that sabotage the conclusions. The first is in the genetic diversity of the organisms, which I argue is fundamentally continuous in nature and a binary construction of the sample is not meaningful. The second misconception is in the assumption that an ecosystem can be decomposed into subsets of independent samples. The third problem is in the definition of an environment. The question of the environment is a continuous one in both space and time and thus any concept constructed to be applicable to a sample of environments must be continuous as well. In summary, the question of whether habitability is a binary quantity or not brings us back to the question of whether our concepts of life and of the environments that life thrives in (or not) are binary or non-binary. I argue that the latter is the case, and hence any modern concept of habitability must be continuous too.


title: Acceleration of nonlinear solvers for natural convection problems

abstract: This paper develops an efficient and robust solution technique for the steady Boussinesq model of non-isothermal flow using Anderson acceleration applied to a Picard iteration. After analyzing the fixed point operator associated with the nonlinear iteration to prove that certain stability and regularity properties hold, we apply the authors' recently constructed theory for Anderson acceleration, which yields a convergence result for the Anderson accelerated Picard iteration for the Boussinesq system. The result shows that the leading term in the residual is improved by the gain in the optimization problem, but at the cost of additional higher order terms that can be significant when the residual is large. We perform numerical tests that illustrate the theory, and show that a 2-stage choice of Anderson depth can be advantageous. We also consider Anderson acceleration applied to the Newton iteration for the Boussinesq equations, and observe that the acceleration allows the Newton iteration to converge for significantly higher Rayleigh numbers that it could without acceleration, even with a standard line search.


title: On 2-local diameter-preserving maps between C(X)-spaces

abstract: The 2-locality problem of diameter-preserving maps between C(X)-spaces is addressed in this paper. For any compact Hausdorff space X with at least three points, we give an example of a 2-local diameter-preserving map on C(X) which is not linear. However, we show that for first countable compact Hausdorff spaces X and Y, every 2-local diameter-preserving map from C(X) to C(Y) is linear and surjective up to constants in some sense. This yields the 2-algebraic reflexivity of isometries with respect to the diameter norms on the quotient spaces.


title: Computational study of the cross-link and the entanglement contributions to the elastic properties of model PDMS networks

abstract: We built randomly cross-linked model PDMS networks and used Molecular Dynamics Methods to obtain stress-strain curves. Mooney-Rivlin (MR) analysis was used to estimate the shear moduli. We applied Primitive Path analysis (PPA) and its variation, Phantom Primitive Path analysis (3PA), to estimate the entanglement and the cross-link moduli, respectively. The MR moduli estimates are in good agreement with the sum of the entanglement and the cross-link moduli, and we observe that the stress-strain data collapse to a universal form when reduced with the PPA and 3PA moduli. We studied how the MR parameters $\mathrm{C}_1$, $\mathrm{C}_2$ vary from cross-link to entanglement dominated networks. For the latter, we observed a $40\%$, $60\%$ contribution of $2\,\mathrm{C}_1$, $2\,\mathrm{C}_2$ to the shear modulus, respectively. Finally, we fitted several models to the data. While all fits are good, the estimates for the entanglement and the cross-link moduli vary significantly when compared to our PPA and 3PA benchmarks.


title: Two halves of a meaningful text are statistically different

abstract: Which statistical features distinguish a meaningful text (possibly written in an unknown system) from a meaningless set of symbols? Here we answer this question by comparing features of the first half of a text to its second half. This comparison can uncover hidden effects, because the halves have the same values of many parameters (style, genre {\it etc}). We found that the first half has more different words and more rare words than the second half. Also, words in the first half are distributed less homogeneously over the text in the sense of of the difference between the frequency and the inverse spatial period. These differences hold for the significant majority of several hundred relatively short texts we studied. The statistical significance is confirmed via the Wilcoxon test. Differences disappear after random permutation of words that destroys the linear structure of the text. The differences reveal a temporal asymmetry in meaningful texts, which is confirmed by showing that texts are much better compressible in their natural way (i.e. along the narrative) than in the word-inverted form. We conjecture that these results connect the semantic organization of a text (defined by the flow of its narrative) to its statistical features.


title: The continuum limit of quantum gravity at first order in perturbation theory

abstract: The Wilsonian renormalization group (RG) properties of the conformal factor of the metric are profoundly altered by the fact that it has a wrong-sign kinetic term. The result is a novel perturbative continuum limit for quantum gravity, which is however non-perturbative in $\hbar$. The ultraviolet part of the renormalized trajectory lies outside the diffeomorphism invariant subspace, entering this subspace only in the infrared, below a dynamically generated amplitude suppression scale. Interactions are dressed with coefficient functions of the conformal factor, their form being determined by the RG. In the ultraviolet, the coefficient functions are parametrised by an infinite number of underlying couplings. Choosing these couplings appropriately, the coefficient functions trivialise on entering the diffeomorphism invariant subspace. Here, dynamically generated effective diffeomorphism couplings emerge, including Newton's constant. In terms of the Legendre effective action, we establish the continuum limit to first order, characterising the most general form of such coefficient functions so as to verify universality.


title: A search for new $\gamma$-ray blazars from infrared selected candidates

abstract: We present a systematic study of gamma-ray blazar candidates based on a sample of 40 objects taken from the WIBR catalogue. By using a likelihood analysis, 26 of the 40 sources showed significant gamma-ray signatures $\geq3\sigma$. Using high-energy test statistics (TS) maps, we confirm 8 sources, which are completely new, and show another 15 promising $\gamma$-ray candidates. The results from this analysis show that a multi-frequency approach can help to improve the current description of the gamma-ray sky.


title: Agent-Based Modelling of Malaria Transmission Dynamics

abstract: Recent statistics of malaria shows that over 200 million cases and estimated deaths of nearly half a million occur globally. Africa alone accounts for almost 90% of the cases. Several studies have been conducted to understand the disease transmission dynamics. In particular, mathematical methods have been frequently used to model and understand the disease dynamics and outbreak patterns. Although, mathematical methods have provided good results for homogeneous populations, these methods impose significant limitations for studying malaria dynamics in heterogeneous populations, a result of various factors, e.g. spatial and temporal fluctuations, social networks, human movements pattern etc. This paper proposes an agent-based modelling approach that permits modelling and analysing malaria dynamics for heterogenous populations. Our approach is illustrated using the climate and demographic data for the Tripura, Limpopo and Benin cities. Our agent-based simulation has been validated against the reported cases of malaria collected in the cities mentioned. Furthermore, the efficiency of the proposed model has been compared with the mathematical model used as benchmark. A statistical test confirms the proposed model is robust and has potential for predicting the peak seasons of malaria. This potentially makes our methods a useful tool as an intervention mechanism, which will have impact on hospitals, healthcare providers, health organisations.


title: Offline and online reconstruction for radio interferometric imaging

abstract: Radio astronomy is transitioning to a big-data era due to the emerging generation of radio interferometric (RI) telescopes, such as the Square Kilometre Array (SKA), which will acquire massive volumes of data. In this article we review methods proposed recently to resolve the ill-posed inverse problem of imaging the raw visibilities acquired by RI telescopes in the big-data scenario. We focus on the recently proposed online reconstruction method [4] and the considerable savings in data storage requirements and computational cost that it yields.


title: Faster Stochastic Quasi-Newton Methods

abstract: Stochastic optimization methods have become a class of popular optimization tools in machine learning. Especially, stochastic gradient descent (SGD) has been widely used for machine learning problems such as training neural networks due to low per-iteration computational complexity. In fact, the Newton or quasi-newton methods leveraging second-order information are able to achieve a better solution than the first-order methods. Thus, stochastic quasi-Newton (SQN) methods have been developed to achieve the better solution efficiently than the stochastic first-order methods by utilizing approximate second-order information. However, the existing SQN methods still do not reach the best known stochastic first-order oracle (SFO) complexity. To fill this gap, we propose a novel faster stochastic quasi-Newton method (SpiderSQN) based on the variance reduced technique of SIPDER. We prove that our SpiderSQN method reaches the best known SFO complexity of $\mathcal{O}(n+n^{1/2}\epsilon^{-2})$ in the finite-sum setting to obtain an $\epsilon$-first-order stationary point. To further improve its practical performance, we incorporate SpiderSQN with different momentum schemes. Moreover, the proposed algorithms are generalized to the online setting, and the corresponding SFO complexity of $\mathcal{O}(\epsilon^{-3})$ is developed, which also matches the existing best result. Extensive experiments on benchmark datasets demonstrate that our new algorithms outperform state-of-the-art approaches for nonconvex optimization.


title: Semi-supervised acoustic modelling for five-lingual code-switched ASR using automatically-segmented soap opera speech

abstract: This paper considers the impact of automatic segmentation on the fully-automatic, semi-supervised training of automatic speech recognition (ASR) systems for five-lingual code-switched (CS) speech. Four automatic segmentation techniques were evaluated in terms of the recognition performance of an ASR system trained on the resulting segments in a semi-supervised manner. The system's output was compared with the recognition rates achieved by a semi-supervised system trained on manually assigned segments. Three of the automatic techniques use a newly proposed convolutional neural network (CNN) model for framewise classification, and include a novel form of HMM smoothing of the CNN outputs. Automatic segmentation was applied in combination with automatic speaker diarization. The best-performing segmentation technique was also tested without speaker diarization. An evaluation based on 248 unsegmented soap opera episodes indicated that voice activity detection (VAD) based on a CNN followed by Gaussian mixture modelhidden Markov model smoothing (CNN-GMM-HMM) yields the best ASR performance. The semi-supervised system trained with the resulting segments achieved an overall WER improvement of 1.1% absolute over the system trained with manually created segments. Furthermore, we found that system performance improved even further when the automatic segmentation was used in conjunction with speaker diarization.


title: The covariance matrix of Green's functions and its application to machine learning

abstract: In this paper, a regression algorithm based on Green's function theory is proposed and implemented. We first survey Green's function for the Dirichlet boundary value problem of 2nd order linear ordinary differential equation, which is a reproducing kernel of a suitable Hilbert space. We next consider a covariance matrix composed of the normalized Green's function, which is regarded as aprobability density function. By supporting Bayesian approach, the covariance matrix gives predictive distribution, which has the predictive mean $\mu$ and the confidence interval [$\mu$-2s, $\mu$+2s], where s stands for a standard deviation.


title: A Malaria Control Model Using Mobility Data: An Early Explanation of Kedougou's Case in Senegal

abstract: Studies in malaria control cover many areas such as medicine, sociology, biology, mathematic, physic, computer science and so forth. Researches in the realm of mathematic are conducted to predict the occurrence of the disease and to support the eradication process. Basically, the modeling methodology is predominantly deterministic and differential equation based while selecting clinical and biological features that seem to be important. Yet, if the individual characteristics matter when modeling the disease, the overall estimation of the malaria is not done based on the health status of each individual but in a non-specified percentage of the global population. The goal of this paper is to propose a model that relies on a daily evolution of the individual's state, which depends on their mobility and the characteristics of the area they visit. Thus, the mobility data of a single person moving from one area to another, gathered thanks to mobile networks, is the essential building block to predict the outcome of the disease. We implement our solution and demonstrate its effectiveness through empirical experiments. The results show how promising the model is in providing possible insights into the failure of the disease control in the Kedougou region.


title: Strategic Investment in Energy Markets: A Multiparametric Programming Approach

abstract: An investor has to carefully select the location and size of new generation units it intends to build, since adding capacity in a market affects the profit from units this investor may already own. To capture this closed-loop characteristic, strategic investment (SI) can be posed as a bilevel optimization. By analytically studying a small market, we first show that its objective function can be non-convex and discontinuous. Realizing that existing mixed-integer problem formulations become impractical for larger markets and increasing number of scenarios, this work put forth two SI solvers: a grid search to handle setups where the candidate investment locations are few, and a stochastic gradient descent approach for otherwise. Both solvers leverage the powerful toolbox of multiparametric programming (MPP), each in a unique way. The grid search entails finding the primal/dual solutions for a large number of optimal power flow (OPF) problems, which nonetheless can be efficiently computed several at once thanks to the properties of MPP. The same properties facilitate the rapid calculation of gradients in a mini-batch fashion, thus accelerating the implementation of a stochastic gradient descent search. Tests on the IEEE 118-bus system using real-world data corroborate the advantages of the novel MPP-aided solvers.


title: Hermite spectral method for Fokker-Planck-Landau equation modeling collisional plasma

abstract: We propose an Hermite spectral method for the Fokker-Planck-Landau (FPL) equation. Both the distribution functions and the collision terms are approximated by series expansions of the Hermite functions. To handle the complexity of the quadratic FPL collision operator, a reduced collision model is built by adopting the quadratic collision operator for the lower-order terms and the diffusive Fokker-Planck operator for the higher-order terms in the Hermite expansion of the reduced collision operator. The numerical scheme is split into three steps according to the Strang splitting, where different expansion centers are employed for different numerical steps to take advantage of the Hermite functions. The standard normalized Hermite basis [36] is adopted during the convection and collision steps to utilize the precalculated coefficients of the quadratic collision terms, while the one constituted by the local macroscopic velocity and temperature is utilized for the acceleration step, by which the effect of the external force can be simplified to an ODE. Projections between different expansion centers are achieved by an algorithm proposed in [29]. Several numerical examples are studied to test and validate our new method


title: In-situ Synthesis of Bismuth/Reduced Graphene Oxide Nanocomposites as High Capacity Anode Materials for Mg-ion Battery

abstract: Herein, in-situ reduction of bismuth and graphene oxide was performed by a solvothermal method under a nitrogen atmosphere, and the resulting Bi on RGO nanocomposites were used as an anode material for Mg-ion batteries. The nanocomposite of 60% Bi:40% RGO is a beneficial anode material, delivering a discharge capacity as high as 413 and 372 mAh per g at the specific current of 39 mA per g in the 1st and 50th cycles, respectively. In addition, it shows high-rate capability with the discharge capacities of 381, 372, 354, 295, and 238 mAh per g at the specific currents of 53, 100, 200, 500, and 700 mA per g, respectively. The better electrochemical performance of the nanocomposite is due to improvement in the electronic conductivity and significant reduction of volume changes during electrochemical cycling. This study demonstrates the bismuth (Bi) on reduced graphene oxide (RGO) nanocomposite as a promising high-capacity anode for magnesium-ion batteries with longer life cycle and high-rate performance.


title: Fast radio bursts from axion stars moving through pulsar magnetospheres

abstract: We study the radio signals generated when an axion star enters the magnetosphere of a neutron star. As the axion star moves through the resonant region where the plasma-induced photon mass becomes equal to the axion mass, the axions can efficiently convert into photons, giving rise to an intense, transient radio signal. We show that a dense axion star with a mass $\sim 10^{-13}M_{\odot}$ composed of $\sim 10\: \mu$eV axions can account for most of the mysterious fast radio bursts.


title: A stochastic reduced-order model for statistical microstructure descriptors evolution

abstract: Integrated Computational Materials Engineering (ICME) models have been a crucial building block for modern materials development, relieving heavy reliance on experiments and significantly accelerating the materials design process. However, ICME models are also computationally expensive, particularly with respect to time integration for dynamics, which hinders the ability to study statistical ensembles and thermodynamic properties of large systems for long time scales. To alleviate the computational bottleneck, we propose to model the evolution of statistical microstructure descriptors as a continuous-time stochastic process using a non-linear Langevin equation, where the probability density function (PDF) of the statistical microstructure descriptors, which are also the quantities of interests (QoIs), are modeled by the Fokker-Planck equation. We discuss how to calibrate the drift and diffusion terms of the Fokker-Planck equation from the theoretical and computational perspectives. The calibrated Fokker-Planck equation can be used as a stochastic reduced-order model (ROM) to simulate the microstructure evolution of statistical microstructure descriptors PDF. Considering statistical microstructure descriptors in the microstructure evolution as QoIs, we demonstrate our proposed methodology in three integrated computational materials engineering (ICME) models: kinetic Monte Carlo, phase field, and molecular dynamics simulations.


title: RTLola Cleared for Take-Off: Monitoring Autonomous Aircraft

abstract: The autonomous control of unmanned aircraft is a highly safety-critical domain with great economic potential in a wide range of application areas, including logistics, agriculture, civil engineering, and disaster recovery. We report on the development of a dynamic monitoring framework for the DLR ARTIS (Autonomous Rotorcraft Testbed for Intelligent Systems) family of unmanned aircraft based on the formal specification language RTLola. RTLola is a stream-based specification language for real-time properties. An RTLola specification of hazardous situations and system failures is statically analyzed in terms of consistency and resource usage and then automatically translated into an FPGA-based monitor. Our approach leads to highly efficient, parallelized monitors with formal guarantees on the noninterference of the monitor with the normal operation of the autonomous system.


title: Peak-structure in self-energy of cuprate superconductors

abstract: The recently deduced normal and anomalous self-energies from photoemission spectra of cuprate superconductors via the machine learning technique are calling for an explanation. Here the normal and anomalous self-energies in cuprate superconductors are analyzed within the framework of the kinetic-energy-driven superconductivity. It is shown that the exchanged spin excitations give rise to the well-pronounced low-energy peak-structures in both the normal and anomalous self-energies, however, they do not cancel in the total self-energy. In particular, the peak-structure in the normal self-energy is mainly responsible for the peak-dip-hump structure in the single-particle excitation spectrum, and can persist into the normal-state, while the sharp peak in the anomalous self-energy gives rise to a crucial contribution to the superconducting gap, and vanishes in the normal-state. Moreover, the evolution of the peak-structure with doping and momentum are also analyzed.


title: PFNN: A Penalty-Free Neural Network Method for Solving a Class of Second-Order Boundary-Value Problems on Complex Geometries

abstract: We present PFNN, a penalty-free neural network method, to efficiently solve a class of second-order boundary-value problems on complex geometries. To reduce the smoothness requirement, the original problem is reformulated to a weak form so that the evaluations of high-order derivatives are avoided. Two neural networks, rather than just one, are employed to construct the approximate solution, with one network satisfying the essential boundary conditions and the other handling the rest part of the domain. In this way, an unconstrained optimization problem, instead of a constrained one, is solved without adding any penalty terms. The entanglement of the two networks is eliminated with the help of a length factor function that is scale invariant and can adapt with complex geometries. We prove the convergence of the PFNN method and conduct numerical experiments on a series of linear and nonlinear second-order boundary-value problems to demonstrate that PFNN is superior to several existing approaches in terms of accuracy, flexibility and robustness.


title: Thermodynamics of structure-forming systems

abstract: Structure-forming systems are ubiquitous in nature, ranging from atoms building molecules to self-assembly of colloidal amphibolic particles. The understanding of the underlying thermodynamics of such systems remains an important problem. Here we derive the entropy for structure-forming systems that differs from Boltzmann-Gibbs entropy by a term that explicitly captures clustered states. For large systems and low concentrations, the approach is equivalent to the grand-canonical ensemble; for small systems, we find significant deviations. We derive the detailed fluctuation theorem and Crooks' work fluctuation theorem for structure-forming systems. The connection to the theory of particle self-assembly is discussed. We apply the results to several physical systems. We present the phase diagram for patchy particles described by the Kern-Frenkel potential. We show that the Curie-Weiss model with molecule structures exhibits a first-order phase transition.


title: Global well-posedness of the half space problem of the Navier-Stokes equations in critical function spaces of limiting case

abstract: In this paper, we study the initial-boundary value problem of the Navier-Stokes equations in half-space. Let a solenoidal initial velocity be given in the function space $ \dot{B}_{p\infty,0}^{ -1 + n/p}({\mathbb R}^n_+)$ for $ \frac{n}3< p < n$. We prove the global in time existence of weak solution $u\in L^\infty(0,\infty; \dot B^{-1 +n/p}_{p\infty}({\mathbb R}^n_+))$, when the given initial velocity has small norm in function space $ \dot{B}_{p\infty,0}^{-1 + n/p} ({\mathbb R}^n_+)$, where $ \frac{n}3< p< n$.


title: Solving Newton's Equations of Motion with Large Timesteps using Recurrent Neural Networks based Operators

abstract: Classical molecular dynamics simulations are based on solving Newton's equations of motion. Using a small timestep, numerical integrators such as Verlet generate trajectories of particles as solutions to Newton's equations. We introduce operators derived using recurrent neural networks that accurately solve Newton's equations utilizing sequences of past trajectory data, and produce energy-conserving dynamics of particles using timesteps up to 4000 times larger compared to the Verlet timestep. We demonstrate significant speedup in many example problems including 3D systems of up to 16 particles.

categories: physics.comp-ph cond-mat.soft cs.LG stat.ML


title: On the Ultraviolet Limit of the Pauli-Fierz Hamiltonian in the Lieb-Loss Model

abstract: Two decades ago, Lieb and Loss proposed to approximate the ground state energy of a free, nonrelativistic electron coupled to the quantized radiation field by the infimum $E_{\alpha, \Lambda}$ of all expectation values $\langle \phi_{el} \otimes \psi_{ph} | H_{\alpha, \Lambda} (\phi_{el} \otimes \psi_{ph}) \rangle$, where $H_{\alpha, \Lambda}$ is the corresponding Hamiltonian with fine structure constant $\alpha >0$ and ultraviolet cutoff $\Lambda < \infty$, and $\phi_{el}$ and $\psi_{ph}$ are normalized electron and photon wave functions, respectively. Lieb and Loss showed that $c \alpha^{1/2} \Lambda^{3/2} \leq E_{\alpha, \Lambda} \leq c^{-1} \alpha^{2/7} \Lambda^{12/7}$ for some constant $c >0$. In the present paper we prove the existence of a constant $C < \infty$, such that \begin{align*} \bigg| \frac{E_{\alpha, \Lambda}}{F[1] \, \alpha^{2/7} \, \Lambda^{12/7}} - 1 \bigg| \ \leq \ C \, \alpha^{4/105} \, \Lambda^{-4/105} \end{align*} holds true, where $F[1] >0$ is an explicit universal number. This result shows that Lieb and Loss' upper bound is actually sharp and gives the asymptotics of $E_{\alpha, \Lambda}$ uniformly in the limit $\alpha \to 0$ and in the ultraviolet limit $\Lambda \to \infty$.


title: Cosmological Constraints on Scalar Field Dark Matter

abstract: This paper aims to put constraints on the parameters of the Scalar Field Dark Matter (SFDM) model, when dark matter is described by a free real scalar field filling the whole Universe, plus a cosmological constant term. By using a compilation of 51 $H(z)$ data and 1048 Supernovae data from Panteon, a lower limit for the mass of the scalar field was obtained, $m \geq 5.1\times 10^{-34} $eV and $H_0=69.5^{+2.0}_{-2.1}\text{ km s}^{-1}\text{Mpc}^{-1}$. Also, the present dark matter density parameter was obtained as $\Omega_\phi = 0.230^{+0.033}_{-0.031}$ at $2\sigma$ confidence level. The results are in good agreement to standard model of cosmology, showing that SFDM model is viable in describing the dark matter content of the universe.


title: Certifiable Robustness to Adversarial State Uncertainty in Deep Reinforcement Learning

abstract: Deep Neural Network-based systems are now the state-of-the-art in many robotics tasks, but their application in safety-critical domains remains dangerous without formal guarantees on network robustness. Small perturbations to sensor inputs (from noise or adversarial examples) are often enough to change network-based decisions, which was recently shown to cause an autonomous vehicle to swerve into another lane. In light of these dangers, numerous algorithms have been developed as defensive mechanisms from these adversarial inputs, some of which provide formal robustness guarantees or certificates. This work leverages research on certified adversarial robustness to develop an online certifiably robust for deep reinforcement learning algorithms. The proposed defense computes guaranteed lower bounds on state-action values during execution to identify and choose a robust action under a worst-case deviation in input space due to possible adversaries or noise. Moreover, the resulting policy comes with a certificate of solution quality, even though the true state and optimal action are unknown to the certifier due to the perturbations. The approach is demonstrated on a Deep Q-Network policy and is shown to increase robustness to noise and adversaries in pedestrian collision avoidance scenarios and a classic control task. This work extends one of our prior works with new performance guarantees, extensions to other RL algorithms, expanded results aggregated across more scenarios, an extension into scenarios with adversarial behavior, comparisons with a more computationally expensive method, and visualizations that provide intuition about the robustness algorithm.


title: Polar tangential angles and free elasticae

abstract: In this note we investigate the behavior of the polar tangential angle of a general plane curve, and in particular prove its monotonicity for certain curves of monotone curvature. As an application we give (non)existence results for an obstacle problem involving free elasticae.


title: Geometry of AdS black hole thermodynamics in extended phase space

abstract: We consider the geometry of anti-de-Sitter (AdS) black hole thermodynamics in four dimensions, where the equation of state in the extended phase space formalism allows explicit comparison with normal fluid systems. We show that for the two-dimensional parameter manifolds considered here, the scalar curvature is proportional to the thermodynamic volume. This allows us to critically examine the applicability of geometric methods in black hole thermodynamics in extended phase space. We show how several standard features that are expected to hold in normal fluid systems impose severe restrictions on the black hole parameters, leading to the fact that several results in the current literature on the geometry of thermodynamics in extended phase space may be physically invalid. These are true for both charged and rotating AdS black holes. As a by-product of our analysis, we examine a conjecture regarding the equality of the correlation lengths of co-existing phases near criticality, in charged AdS black hole backgrounds, and find reasonable validity.


title: What's so special about BERT's layers? A closer look at the NLP pipeline in monolingual and multilingual models

abstract: Peeking into the inner workings of BERT has shown that its layers resemble the classical NLP pipeline, with progressively more complex tasks being concentrated in later layers. To investigate to what extent these results also hold for a language other than English, we probe a Dutch BERT-based model and the multilingual BERT model for Dutch NLP tasks. In addition, through a deeper analysis of part-of-speech tagging, we show that also within a given task, information is spread over different parts of the network and the pipeline might not be as neat as it seems. Each layer has different specialisations, so that it may be more useful to combine information from different layers, instead of selecting a single one based on the best overall performance.


title: SpaceNet 6: Multi-Sensor All Weather Mapping Dataset

abstract: Within the remote sensing domain, a diverse set of acquisition modalities exist, each with their own unique strengths and weaknesses. Yet, most of the current literature and open datasets only deal with electro-optical (optical) data for different detection and segmentation tasks at high spatial resolutions. optical data is often the preferred choice for geospatial applications, but requires clear skies and little cloud cover to work well. Conversely, Synthetic Aperture Radar (SAR) sensors have the unique capability to penetrate clouds and collect during all weather, day and night conditions. Consequently, SAR data are particularly valuable in the quest to aid disaster response, when weather and cloud cover can obstruct traditional optical sensors. Despite all of these advantages, there is little open data available to researchers to explore the effectiveness of SAR for such applications, particularly at very-high spatial resolutions, i.e. <1m Ground Sample Distance (GSD). To address this problem, we present an open Multi-Sensor All Weather Mapping (MSAW) dataset and challenge, which features two collection modalities (both SAR and optical). The dataset and challenge focus on mapping and building footprint extraction using a combination of these data sources. MSAW covers 120 km^2 over multiple overlapping collects and is annotated with over 48,000 unique building footprints labels, enabling the creation and evaluation of mapping algorithms for multi-modal data. We present a baseline and benchmark for building footprint extraction with SAR data and find that state-of-the-art segmentation models pre-trained on optical data, and then trained on SAR (F1 score of 0.21) outperform those trained on SAR data alone (F1 score of 0.135).


title: Dynamical cloud formation traced by atomic and molecular gas

abstract: Context: Atomic and molecular cloud formation is a dynamical process. However, kinematic signatures of these processes are still observationally poorly constrained. Methods: Targeting the cloud-scale environment of the prototypical infrared dark cloud G28.3, we employ spectral line imaging observations of the two atomic lines HI and [CI] as well as molecular lines observations in 13CO in the 1--0 and 3--2 transitions. The analysis comprises investigations of the kinematic properties of the different tracers, estimates of the mass flow rates, velocity structure functions, a Histogram of Oriented Gradients (HOG) study as well as comparisons to simulations. Results: The central IRDC is embedded in a more diffuse envelope of cold neutral medium (CNM) traced by HI self-absorption (HISA) and molecular gas. The spectral line data as well as the HOG and structure function analysis indicate a possible kinematic decoupling of the HI from the other gas compounds. Spectral analysis and position-velocity diagrams reveal two velocity components that converge at the position of the IRDC. Estimated mass flow rates appear rather constant from the cloud edge toward the center. The velocity structure function analysis is consistent with gas flows being dominated by the formation of hierarchical structures. Conclusions: The observations and analysis are consistent with a picture where the IRDC G28 is formed at the center of two converging gas flows. While the approximately constant mass flow rates are consistent with a self-similar, gravitationally driven collapse of the cloud, external compression by, e.g., spiral arm shocks or supernovae explosions cannot be excluded yet. Future investigations should aim at differentiating the origin of such converging gas flows.


title: Unsupervised Multimodal Video-to-Video Translation via Self-Supervised Learning

abstract: Existing unsupervised video-to-video translation methods fail to produce translated videos which are frame-wise realistic, semantic information preserving and video-level consistent. In this work, we propose UVIT, a novel unsupervised video-to-video translation model. Our model decomposes the style and the content, uses the specialized encoder-decoder structure and propagates the inter-frame information through bidirectional recurrent neural network (RNN) units. The style-content decomposition mechanism enables us to achieve style consistent video translation results as well as provides us with a good interface for modality flexible translation. In addition, by changing the input frames and style codes incorporated in our translation, we propose a video interpolation loss, which captures temporal information within the sequence to train our building blocks in a self-supervised manner. Our model can produce photo-realistic, spatio-temporal consistent translated videos in a multimodal way. Subjective and objective experimental results validate the superiority of our model over existing methods. More details can be found on our project website: https://uvit.netlify.com


title: Computationally efficient models for the dominant and sub-dominant harmonic modes of precessing binary black holes

abstract: We present IMRPhenomXPHM, a phenomenological frequency-domain model for the gravitational-wave signal emitted by quasi-circular precessing binary black holes, which incorporates multipoles beyond the dominant quadrupole in the precessing frame. The model is a precessing extension of IMRPhenomXHM (Garc\'ia-Quir\'os 2020), based on approximate maps between aligned-spin waveform modes in the co-precessing frame and precessing waveform modes in the inertial frame, which is commonly referred to as "twisting up" the non-precessing waveforms. IMRPhenomXPHM includes IMRPhenomXP as a special case, the restriction to the dominant quadrupole contribution in the co-precessing frame. We implement two alternative mappings, one based on a single-spin PN approximation, as used in IMRPhenomPv2 (Hannam 2013), and one based on the double-spin MSA approach (Chatziioannou 2017). We include a detailed discussion of conventions used in the description of precessing binaries and of all choices made in constructing the model. The computational cost of \phXPHM is further reduced by extending the interpolation technique of (C. Garc\'ia-Quir\'os 2020) to the Euler angles. The accuracy, speed, robustness and modularity of the IMRPhenomX family will make these models productive tools for gravitational wave astronomy in the current era of greatly increased number and diversity of detected events.


title: Effect of Autonomous Driving on Traffic Breakdown in Mixed Traffic Flow: A Critical Mini-Review

abstract: A critical analysis of the effect of autonomous driving vehicles on traffic breakdown in mixed traffic flow is made. Classical adaptive cruise control (ACC) and ACC based on three-phase traffic theory (TPACC) are considered. We show that within a wide range of dynamic parameters of classical ACC, the ACC-vehicles can initiate traffic breakdown and reduce highway capacity. Contrarily, in the same range of parameters of TPACC, the TPACC-vehicles do not initiate traffic breakdown. To understand physical reasons for the effect of classical ACC- and TPACC-vehicles on traffic breakdown, we introduce a model of ACC that can be considered a combination of dynamic features of classical ACC and TPACC. With the use of this model, we find how the amplitude of a local speed disturbance caused by the ACC in a vicinity of a bottleneck and the probability of traffic breakdown depend on the dynamic parameters of the ACC.


title: Weak KAM approach to first-order Mean Field Games with state constraints

abstract: We study the asymptotic behavior of solutions to the constrained MFG system as the time horizon $T$ goes to infinity. For this purpose, we analyze first Hamilton-Jacobi equations with state constraints from the viewpoint of weak KAM theory, constructing a Mather measure for the associated variational problem. Using these results, we show that a solution to the constrained ergodic mean field games system exists and the ergodic constant is unique. Finally, we prove that any solution of the first-order constrained MFG problem on $[0,T]$ converges to the solution of the ergodic system as $T \to +\infty$.


title: Nilpotent groups of automorphisms of families of Riemann surfaces

abstract: In this article we extend results of Zomorrodian to determine upper bounds for the order of a nilpotent group of automorphisms of a complex $d$-dimensional family of compact Riemann surfaces, where $d \geqslant 1.$ We provide conditions under which these bounds are sharp. In addition, for the one-dimensional case we construct and describe an explicit family attaining the bound for infinitely many genera. We obtain similar results for the case of $p$-groups of automorphisms.


title: C++ Modules in ROOT and Beyond

abstract: C++ Modules come in C++20 to fix the long-standing build scalability problems in the language. They provide an io-efficient, on-disk representation capable to reduce build times and peak memory usage. ROOT employs the C++ modules technology further in the ROOT dictionary system to improve its performance and reduce the memory footprint. ROOT with C++ Modules was released as a technology preview in fall 2018, after intensive development during the last few years. The current state is ready for production, however, there is still room for performance optimizations. In this talk, we show the roadmap for making the technology default in ROOT. We demonstrate a global module indexing optimization which allows reducing the memory footprint dramatically for many workflows. We will report user feedback on the migration to ROOT with C++ Modules.


title: The growth rate over trees of any family of set defined by a monadic second order formula is semi-computable

abstract: Monadic second order logic can be used to express many classical notions of sets of vertices of a graph as for instance: dominating sets, induced matchings, perfect codes, independent sets or irredundant sets. Bounds on the number of sets of any such family of sets are interesting from a combinatorial point of view and have algorithmic applications. Many such bounds on different families of sets over different classes of graphs are already provided in the literature. In particular, Rote recently showed that the number of minimal dominating sets in trees of order $n$ is at most $95^{\frac{n}{13}}$ and that this bound is asymptotically sharp up to a multiplicative constant. We build on his work to show that what he did for minimal dominating sets can be done for any family of sets definable by a monadic second order formula. We first show that, for any monadic second order formula over graphs that characterizes a given kind of subset of its vertices, the maximal number of such sets in a tree can be expressed as the \textit{growth rate of a bilinear system}. This mostly relies on well known links between monadic second order logic over trees and tree automata and basic tree automata manipulations. Then we show that this "growth rate" of a bilinear system can be approximated from above.We then use our implementation of this result to provide bounds on the number of independent dominating sets, total perfect dominating sets, induced matchings, maximal induced matchings, minimal perfect dominating sets, perfect codes and maximal irredundant sets on trees. We also solve a question from D. Y. Kang et al. regarding $r$-matchings and improve a bound from G\'orska and Skupie\'n on the number of maximal matchings on trees. Remark that this approach is easily generalizable to graphs of bounded tree width or clique width (or any similar class of graphs where tree automata are meaningful).


title: Directional square functions

abstract: Quantitative formulations of Fefferman's counterexample for the ball multiplier are naturally linked to square function estimates for conical and directional multipliers. In this article we develop a novel framework for these square function estimates, based on a directional embedding theorem for Carleson sequences and multi-parameter time-frequency analysis techniques. As applications we prove sharp or quantified bounds for Rubio de Francia type square functions of conical multipliers and of multipliers adapted to rectangles pointing along $N$ directions. A suitable combination of these estimates yields a new and currently best-known logarithmic bound for the Fourier restriction to an $N$-gon, improving on previous results of A. Cordoba. Our directional Carleson embedding extends to the weighted setting, yielding previously unknown weighted estimates for directional maximal functions and singular integrals.


title: Hi Sigma, do I have the Coronavirus?: Call for a New Artificial Intelligence Approach to Support Health Care Professionals Dealing With The COVID-19 Pandemic

abstract: Just like your phone can detect what song is playing in crowded spaces, we show that Artificial Intelligence transfer learning algorithms trained on cough phone recordings results in diagnostic tests for COVID-19. To gain adoption by the health care community, we plan to validate our results in a clinical trial and three other venues in Mexico, Spain and the USA . However, if we had data from other on-going clinical trials and volunteers, we may do much more. For example, for confirmed stay-at-home COVID-19 patients, a longitudinal audio test could be developed to determine contact-with-hospital recommendations, and for the most critical COVID-19 patients a success ratio forecast test, including patient clinical data, to prioritize ICU allocation. As a challenge to the engineering community and in the context of our clinical trial, the authors suggest distributing cough recordings daily, hoping other trials and crowdsourcing users will contribute more data. Previous approaches to complex AI tasks have either used a static dataset or were private efforts led by large corporations. All existing COVID-19 trials published also follow this paradigm. Instead, we suggest a novel open collective approach to large-scale real-time health care AI. We will be posting updates at https://opensigma.mit.edu. Our personal view is that our approach is the right one for large scale pandemics, and therefore is here to stay - will you join?


title: Percolation with plasticity for neuromorphic systems

abstract: We develop a theory of percolation with plasticity systems (PWPs) rendering properties of interest for neuromorphic computing. Unlike the standard percolation between two large electrodes, they have multiple ($N\gg 1$) interfaces and exponentially large number ($N!$) of conductive pathways between them. These pathways consist of non-ohmic random resistors that can undergo bias induced nonvolatile modifications (plasticity). The neuromorphic properties of PWPs include: multi-valued memory, high dimensionality and nonlinearity capable of transforming input data into spatiotemporal patterns, tunably fading memory ensuring outputs that depend more on recent inputs, and no need for massive interconnects. A few conceptual examples of functionality here are random number generation, matrix-vector multiplication, and associative memory. Understanding PWP topology, statistics, and operations opens a field of its own calling upon further theoretical and experimental insights.


title: Exact analytical solution of a novel modified nonlinear Schr\"odinger equation: solitary quantum waves on a lattice

abstract: A novel modified nonlinear Schr\"odinger equation is presented. Through a travelling wave ansatz, the equation is transformed into a nonlinear ODE which is then solved exactly and analytically. The soliton solution is characterised in terms of waveform and wave speed, and the dependence of these properties upon parameters in the equation is detailed. It is discovered that some parameter settings yield unique waveforms while others yield degeneracy, with two distinct waveforms per set of parameter values. The uni-waveform and bi-waveform regions of parameter space are identified. It is also found that each waveform has two modes of propagation with shared directionality but distinct speeds. Finally, the equation is shown to be a model for the propagation of a quantum mechanical exciton, such as an electron, through a collectively-oscillating plane lattice with which the exciton interacts. The physical implications of the soliton solution are discussed.


title: Reaction-diffusion equation on thin porous media

abstract: We consider a reaction-diffusion equation on a 3D thin porous media of thickness $\varepsilon$ which is perforated by periodically distributed cylinders of size $\varepsilon$. On the boundary of the cylinders we prescribe a dynamical boundary condition of pure-reactive type. As $\varepsilon\to 0$, in the 2D limit the resulting reaction-diffusion equation has a source term coming from the dynamical-type boundary conditions imposed on boundaries of the original 3D domain.


title: Estimating state occupation and transition probabilities in non-Markov multi-state models subject to both random left-truncation and right-censoring

abstract: The Aalen-Johansen estimator generalizes the Kaplan-Meier estimator for independently left-truncated and right-censored survival data to estimating the transition probability matrix of a time-inhomogeneous Markov model with finite state space. Such multi-state models have a wide range of applications for modelling complex courses of a disease over the course of time, but the Markov assumption may often be in doubt. If censoring is entirely unrelated to the multi-state data, it has been noted that the Aalen-Johansen estimator, standardized by the initial empirical distribution of the multi-state model, still consistently estimates the state occupation probabilities. Recently, this result has been extended to transition probabilities using landmarking, which is, inter alia, useful for dynamic prediction. We complement these results in three ways. Firstly, delayed study entry is a common phenomenon in observational studies, and we extend the earlier results to multi-state data also subject to left-truncation. Secondly, we present a rigorous proof of consistency of the Aalen-Johansen estimator for state occupation probabilities, on which also correctness of the landmarking approach hinges, correcting, simplifying and extending the earlier result. Thirdly, our rigorous proof motivates wild bootstrap resampling. Our developments for left-truncation are motivated by a prospective observational study on the occurrence and the impact of a multi-resistant infectious organism in patients undergoing surgery. Both the real data example and simulation studies are presented. Studying wild bootstrap is motivated by the fact that, unlike drawing with replacement from the data, it is desirable to have a technique that works both with non-Markov models subject to random left-truncation and right-censoring and with Markov models where left-truncation and right-censoring need not be entirely random.


title: Emergence of complex structures from nonlinear interactions and noise in coevolving networks

abstract: We study the joint effect of the non-linearity of interactions and noise on coevolutionary dynamics. We choose the coevolving voter model as a prototype framework for this problem. By numerical simulations and analytical approximations we find three main phases that differ in the absolute magnetization and the size of the largest component: a consensus phase, a coexistence phase, and a dynamical fragmentation phase. More detailed analysis reveals inner differences in these phases, allowing us to divide two of them further. In the consensus phase we can distinguish between a weak or alternating consensus (switching between two opposite consensus states), and a strong consensus, in which the system remains in the same state for the whole realization of the stochastic dynamics. Additionally, weak and strong consensus phases scale differently with the system size. The strong consensus phase exists for superlinear interactions and it is the only consensus phase that survives in the thermodynamic limit. In the coexistence phase we distinguish a fully-mixing phase (both states well mixed in the network) and a structured coexistence phase, where the number of links connecting nodes in different states (active links) drops significantly due to the formation of two homogeneous communities of opposite states connected by a few links. The structured coexistence phase is an example of emergence of community structure from not exclusively topological dynamics, but coevolution. Our numerical observations are supported by an analytical description using a pair approximation approach and an ad-hoc calculation for the transition between the coexistence and dynamical fragmentation phases. Our work shows how simple interaction rules including the joint effect of non-linearity, noise, and coevolution lead to complex structures relevant in the description of social systems.


title: Quantum vs. classical algorithms for solving the heat equation

abstract: Quantum computers are predicted to outperform classical ones for solving partial differential equations, perhaps exponentially. Here we consider a prototypical PDE - the heat equation in a rectangular region - and compare in detail the complexities of ten classical and quantum algorithms for solving it, in the sense of approximately computing the amount of heat in a given region. We find that, for spatial dimension $d \ge 2$, there is an at most quadratic quantum speedup using an approach based on applying amplitude estimation to an accelerated classical random walk. However, an alternative approach based on a quantum algorithm for linear equations is never faster than the best classical algorithms.


title: Learning a low dimensional manifold of real cancer tissue with PathologyGAN

abstract: Application of deep learning in digital pathology shows promise on improving disease diagnosis and understanding. We present a deep generative model that learns to simulate high-fidelity cancer tissue images while mapping the real images onto an interpretable low dimensional latent space. The key to the model is an encoder trained by a previously developed generative adversarial network, PathologyGAN. We study the latent space using 249K images from two breast cancer cohorts. We find that the latent space encodes morphological characteristics of tissues (e.g. patterns of cancer, lymphocytes, and stromal cells). In addition, the latent space reveals distinctly enriched clusters of tissue architectures in the high-risk patient group.


title: Gender Detection on Social Networks using Ensemble Deep Learning

abstract: Analyzing the ever-increasing volume of posts on social media sites such as Facebook and Twitter requires improved information processing methods for profiling authorship. Document classification is central to this task, but the performance of traditional supervised classifiers has degraded as the volume of social media has increased. This paper addresses this problem in the context of gender detection through ensemble classification that employs multi-model deep learning architectures to generate specialized understanding from different feature spaces.


title: Compass-aligned Distributional Embeddings for Studying Semantic Differences across Corpora

abstract: Word2vec is one of the most used algorithms to generate word embeddings because of a good mix of efficiency, quality of the generated representations and cognitive grounding. However, word meaning is not static and depends on the context in which words are used. Differences in word meaning that depends on time, location, topic, and other factors, can be studied by analyzing embeddings generated from different corpora in collections that are representative of these factors. For example, language evolution can be studied using a collection of news articles published in different time periods. In this paper, we present a general framework to support cross-corpora language studies with word embeddings, where embeddings generated from different corpora can be compared to find correspondences and differences in meaning across the corpora. CADE is the core component of our framework and solves the key problem of aligning the embeddings generated from different corpora. In particular, we focus on providing solid evidence about the effectiveness, generality, and robustness of CADE. To this end, we conduct quantitative and qualitative experiments in different domains, from temporal word embeddings to language localization and topical analysis. The results of our experiments suggest that CADE achieves state-of-the-art or superior performance on tasks where several competing approaches are available, yet providing a general method that can be used in a variety of domains. Finally, our experiments shed light on the conditions under which the alignment is reliable, which substantially depends on the degree of cross-corpora vocabulary overlap.


title: Dark defect charge dynamics in bulk chemical-vapor-deposition-grown diamonds probed via nitrogen vacancy centers

abstract: Although chemical vapor deposition (CVD) is one of the preferred routes to synthetic diamond crystals, a full knowledge of the point defects produced during growth is still incomplete. Here we exploit the charge and spin properties of nitrogen-vacancy (NV) centers in type-1b CVD diamond to expose an optically and magnetically dark point defect, so far virtually unnoticed despite an abundance comparable to (if not greater than) that of substitutional nitrogen. Indirectly-detected photo-luminescence spectroscopy indicates a donor state 1.6 eV above the valence band, although the defect's microscopic structure and composition remain elusive. Our results may prove relevant to the growing set of applications that rely on CVD-grown single crystal diamond.


title: Quantum speedups of some general-purpose numerical optimisation algorithms

abstract: We give quantum speedups of several general-purpose numerical optimisation methods for minimising a function $f:\mathbb{R}^n \to \mathbb{R}$. First, we show that many techniques for global optimisation under a Lipschitz constraint can be accelerated near-quadratically. Second, we show that backtracking line search, an ingredient in quasi-Newton optimisation algorithms, can be accelerated up to quadratically. Third, we show that a component of the Nelder-Mead algorithm can be accelerated by up to a multiplicative factor of $O(\sqrt{n})$. Fourth, we show that a quantum gradient computation algorithm of Gily\'en et al. can be used to approximately compute gradients in the framework of stochastic gradient descent. In each case, our results are based on applying existing quantum algorithms to accelerate specific components of the classical algorithms, rather than developing new quantum techniques.


title: Prospects and limits of SIR-type Mathematical Models to Capture the COVID-19 Pandemic

abstract: For the description of a pandemic mathematical models could be interesting. Both for physicians and politicians as a base for decisions to treat the disease. The responsible estimation of parameters is a main issue of mathematical pandemic models. Especially a good choice of $\beta$ as the number of others that one infected person encounters per unit time (per day) influences the adequateness of the results of the model. For the actual COVID-19 pandemic some aspects of the parameter choice will be discussed. Because of the incompatibility of the data of the Johns-Hopkins-University to the data of the German Robert-Koch-Institut we use the COVID-19 data of the European Centre for Disease Prevention and Control (ECDC) as a base for the parameter estimation. Two different mathematical methods for the data analysis will be discussed in this paper and possible sources of trouble will be shown. As example of the parameter choice serve the data of the USA and the UK. The resulting parameters will be used estimated and used in W.\,O. Kermack and A.\,G. McKendrick's SIR model. Strategies for the commencing and ending of social and economic shutdown measures are discussed. The numerical solution of the ordinary differential equation system of the modified SIR model is being done with a Runge-Kutta integration method of fourth order. At the end the applicability of the SIR model could be shown essentially. Suggestions about appropriate points in time at which to commence with lockdown measures based on the acceleration rate of infections conclude the paper.


title: Mass spectrum of pseudo-scalar glueballs from a Bethe-Salpeter approach with the rainbow-ladder truncation

abstract: We suggest a framework based on the rainbow approximation to the Dyson-Schwinger and Bethe-Salpeter equations with effective parameters adjusted to lattice QCD data to calculate the masses of the ground and excited states of pseudo-scalar glueballs. The structure of the truncated Bethe-Salpeter equation with the gluon and ghost propagators as solutions of the truncated Dyson-Schwinger equations is analysed in Landau gauge. Both, the Bete-Salpeter and Dyson-Schwinger equations, are solved numerically within the same rainbow-ladder truncation with the same effective parameters which ensure consistency of the approach. We found that with a set of parameters, which provides a good description of the lattice data within the Dyson-Schwinger approach, the solutions of the Bethe Salpeter equation for the pseudo-scalar glueballs exhibit a rich mass spectrum which also includes the ground and excited states predicted by lattice calculations. The obtained mass spectrum contains also several intermediate excitations beyond the lattice approaches. The partial Bethe-Salpeter amplitudes of the pseudo-scalar glueballs are presented as well.


title: Contrastive Examples for Addressing the Tyranny of the Majority

abstract: Computer vision algorithms, e.g. for face recognition, favour groups of individuals that are better represented in the training data. This happens because of the generalization that classifiers have to make. It is simpler to fit the majority groups as this fit is more important to overall error. We propose to create a balanced training dataset, consisting of the original dataset plus new data points in which the group memberships are intervened, minorities become majorities and vice versa. We show that current generative adversarial networks are a powerful tool for learning these data points, called contrastive examples. We experiment with the equalized odds bias measure on tabular data as well as image data (CelebA and Diversity in Faces datasets). Contrastive examples allow us to expose correlations between group membership and other seemingly neutral features. Whenever a causal graph is available, we can put those contrastive examples in the perspective of counterfactuals.


title: Shadow of a Naked Singularity without Photon Sphere

abstract: It is generally believed that the shadows of either a black hole or naked singularity arise due to photon spheres developing in these spacetimes. Here we propose a new spherically symmetric naked singularity spacetime solution of Einstein equations which has no photon sphere, and we show that the singularity casts a shadow in the absence of the photon sphere. We discuss some novel features of this shadow and the lightlike geodesics in this spacetime. We compare the shadow of the naked singularity here with shadows cast by Schwarzschild black hole and the first type of Joshi-Malafarina-Narayan (JMN1) naked singularity, where for the last two spacetimes the shadow is formed due to the presence of a photon sphere. It is seen, in particular, that the size of shadow of the singularity is considerably smaller than that of a black hole. Our analysis shows that the shadow of this naked singularity is distinguishable from the shadow of a Schwarzschild black hole and the JMN1 naked singularity. These results are useful and important in the context of recent observations of shadow of the M87 galactic center.


title: Adaptive degenerate space method for source term estimation using a backward Lagrangian stochastic model

abstract: The general problem of characterizing gas source parameters based on concentration measurements is known to be a difficult task. As many inverse problems, one of the main obstacles for accurate estimation is the non-uniqueness of solution, induced by the lack of sufficient information. As the number of detectors is lowered, which is more than a plausible scenario for many practical situations, the number of possible solutions that can characterize the source increases dramatically, leading to severe errors. In this paper, a Lagrangian stochastic based method for identifying these suspected points, which will be referred to as 'degenerate space', is formulated and analysed. Then, a new procedure for quantitative prediction of the effect of deploying a new detector in space is used to design an adaptive scheme for source term estimation. This scheme has been tested for several scenarios, differing by the location of the initial detectors, and is shown to reduce dramatically the degeneracy formed by insufficient information. The combined formulation of degenerate space with the new adaptive scheme is shown to give improved accuracy, and in particular for a relatively small number of detectors.

categories: physics.comp-ph astro-ph.IM physics.ao-ph physics.flu-dyn physics.geo-ph