title: Ratios of the hadronic contributions to the lepton $g-2$ from Lattice QCD+QED simulations

abstract: The ratios among the leading-order (LO) hadronic vacuum polarization (HVP) contributions to the anomalous magnetic moments of electron, muon and tau-lepton, $a_{\ell=e,\mu \tau}^{HVP,LO}$, are computed using lattice QCD+QED simulations. The results include the effects at order $O(\alpha_{em}^2)$ as well as the electromagnetic and strong isospin-breaking corrections at orders $O(\alpha_{em}^3)$ and $O(\alpha_{em}^2(m_u-m_d))$, respectively, where $(m_u-m_d)$ is the $u$- and $d$-quark mass difference. We employ the gauge configurations generated by the Extended Twisted Mass Collaboration with $N_f=2+1+1$ dynamical quarks at three values of the lattice spacing ($a \simeq 0.062, 0.082, 0.089$ fm) with pion masses in the range 210 - 450 MeV. We show that in the case of the electron-muon ratio the hadronic uncertainties in the numerator and in the denominator largely cancel out, while in the cases of the electron-tau and muon-tau ratios such a cancellation does not occur. For the electron-muon ratio we get $R_{e/\mu } \equiv (m_\mu/m_e)^2 (a_e^{HVP,LO} / a_\mu^{HVP,LO}) = 1.1456~(83)$ with an uncertainty of $\simeq 0.7 \%$. Our result, which represents an accurate Standard Model (SM) prediction, agrees very well with the estimate obtained using the results of dispersive analyses of the experimental $e^+ e^- \to$ hadrons data. Instead, it differs by $\simeq 2.7$ standard deviations from the value expected from present electron and muon (g - 2) experiments after subtraction of the current estimates of the QED, electro-weak, hadronic light-by-light and higher-order HVP contributions, namely $R_{e/\mu} = 0.575~(213)$. An improvement of the precision of both the experiment and the QED contribution to the electron (g - 2) by a factor of $\simeq 2$ could be sufficient to reach a tension with our SM value of the ratio $R_{e/\mu }$ at a significance level of $\simeq 5$ standard deviations.


title: Parallel Quantum Simulation of Large Systems on Small Quantum Computers

abstract: Tensor networks permit computational and entanglement resources to be concentrated in interesting regions of Hilbert space. Implemented on NISQ machines they allow simulation of quantum systems that are much larger than the computational machine itself. This is achieved by parallelising the quantum simulation. Here, we demonstrate this in the simplest case; an infinite, translationally invariant quantum spin chain. We provide Cirq and Qiskit code that translate infinite, translationally invariant matrix product state (iMPS) algorithms to finite-depth quantum circuit machines, allowing the representation, optimisation and evolution arbitrary one-dimensional systems. Illustrative simulated output of these codes for achievable circuit sizes is given.


title: An Efficient Machine Learning Approach for Accurate Short Term Solar Power Prediction

abstract: Solar based electricity generations have experienced strong and impactful growth in recent years. The regulation, scheduling, dispatching, and unit commitment of intermittent solar power is dependent on the accuracy of the forecasting methods. In this paper, a robust Expanded Extreme Learning Machine (EELM) is proposed to accurately predict solar power for different time horizons and weather conditions. The proposed EELM technique has no randomness due to the absence of random input layer weights and takes very less time to predict the solar power efficiently. The performance of the proposed EELM is validated through historical data collected from the National Renewable Energy Laboratory (NREL) through various performance metrics. The efficacy of the proposed EELM method is evaluated against basic ELM and Functional Link Neural Network (FLNN) for 5 minutes and 1 hour ahead time horizon.


title: Population dynamics of Varroa mite and honeybee: Effects of parasitism with age structure and seasonality

abstract: Honeybees play an important role in the production of many agricultural crops and in sustaining plant diversity in undisturbed ecosystems. The rapid decline of honeybee populations have sparked great concern worldwide. Previous studies have shown that the parasitic Varroa mite could be the main reason for colony losses. In order to understand how mites affect population dynamics of honeybees and a colony health, we propose a brood-adult bee-mite model in which the time lag from brood to adult is taken into account. Noting that the dynamics of a honeybee colony varies with respect to season, we validate the model and perform parameter estimations under both constant and fluctuating seasonality scenarios. Our analytical and numerical studies reveal the following: (a) In the presence of parasite mites, the large time lag from brood to adult could destabilize population dynamics and drive the colony to collapse; but the small natural mortality of the adult population can promote a mite-free colony when time lag is small or at an intermediate level; (b) Small brood' infestation rates could stabilize all populations at the unique interior equilibrium under constant seasonality while may drive the mite population to die out when seasonality is taken into account; (c) High brood' infestation rates can destabilize the colony dynamics leading to population collapse depending on initial population size under constant and seasonal conditions; (d) Results from sensitivity analysis indicate the queen's egg-laying may have the greatest effect on colony population size. The brood death rate and the colony size at which brood survivability is the half maximal were also shown to be highly sensitive with an inverse correlation to the colony population size. Our results provide insights on the effects of seasonality on the dynamics.


title: Properties of the LWR model with time delay

abstract: In this article, we investigate theoretical and numerical properties of the first-order Lighthill-Whitham-Richards (LWR) traffic flow model with time delay. Since standard results from the literature are not directly applicable to the delayed model, we mainly focus on the numerical analysis of the proposed finite difference discretization. The simulation results also show that the delay model is able to capture Stop & Go waves.


title: Online and Real-time Object Tracking Algorithm with Extremely Small Matrices

abstract: Online and Real-time Object Tracking is an interesting workload that can be used to track objects (e.g., car, human, animal) in a series of video sequences in real-time. For simple object tracking on edge devices, the output of object tracking could be as simple as drawing a bounding box around a detected object and in some cases, the input matrices used in such computation are quite small (e.g., 4x7, 3x3, 5x5, etc). As a result, the amount of actual work is low. Therefore, a typical multi-threading based parallelization technique can not accelerate the tracking application; instead, a throughput based parallelization technique where each thread operates on independent video sequences is more rewarding. In this paper, we share our experience in parallelizing a Simple Online and Real-time Tracking (SORT) application on shared-memory multicores.


title: Effective Medium Transformation: the Case of Stratified Magnetic Structures

abstract: Effective medium theory replaces a given fine-scale heterostructure with a homogeneous one in such a way that the physically measurable quantities, e.g. reaction fields and losses, remain approximately the same. This Letter shows that the very nature of the physical problem may change upon homogenization. A specific example is a stratified nonlinear magnetic and conducting medium, where a low frequency excitation induces eddy currents. It is shown that the appropriate coarse-scale (homogeneous) model is, counter intuitively, magnetostatic, with an effective complex valued BH curve whose real and imaginary parts represent active and reactive losses in the sample. Similar situations may arise in other physical and engineering applications, notably, in diffusion problems with boundary layers.


title: To Tweet or Not to Tweet: Covertly Manipulating a Twitter Debate on Vaccines Using Malware-Induced Misperceptions

abstract: Trolling and social bots have been proven as powerful tactics for manipulating the public opinion and sowing discord among Twitter users. This effort requires substantial content fabrication and account coordination to evade Twitter's detection of nefarious platform use. In this paper we explore an alternative tactic for covert social media interference by inducing misperceptions about genuine, non-trolling content from verified users. This tactic uses a malware that covertly manipulates targeted words, hashtags, and Twitter metrics before the genuine content is presented to a targeted user in a covert man-in-the-middle fashion. Early tests of the malware found that it is capable of achieving a similar goal as trolls and social bots, that is, silencing or provoking social media users to express their opinion in polarized debates on social media. Following this, we conducted experimental tests in controlled settings (N=315) where the malware covertly manipulated the perception in a Twitter debate on the risk of vaccines causing autism. The empirical results demonstrate that inducing misperception is an effective tactic to silence users on Twitter when debating polarizing issues like vaccines. We used the findings to propose a solution for countering the effect of the malware-induced misperception that could also be used against trolls and social bots on Twitter.


title: Tactile sensing and computing on a random network of conducting fluid channels

abstract: Liquid electronic sensors are typically based on regular arrays of channels filled with a conductive liquid. We propose that a random planar network of conducting liquid allows us for a wider spectrum of electrical responses and localisation of tactile stimuli. We also speculate that a computation protocol can be implemented on such a network, featuring mechanical inputs and electrical outputs. Our results pave a way towards future developments on sensing and computing wearables with disordered sensing networks structure.


title: Concerning Quantum Identification Without Entanglement

abstract: Identification schemes are interactive protocols typically involving two parties, a prover, who wants to provide evidence of his or her identity and a verifier, who checks the provided evidence and decide whether it comes or not from the intended prover. In this paper, we comment on a recent proposal for quantum identity authentication from Zawadzki, and give a concrete attack upholding theoretical impossibility results from Lo and Buhrman et al. More precisely, we show that using a simple strategyan adversary may indeed obtain non-negligible information on the shared identification secret. While the security of a quantum identity authentication scheme is not formally defined in [1], it is clear that such a definition should somehow imply that an external entity may gain no information on the shared identification scheme (even if he actively participates injecting messages in a protocol execution, which is not assumed in our attack strategy).


title: Engineering Fast High-Fidelity Quantum Operations With Constrained Interactions

abstract: Understanding how to tailor quantum dynamics to achieve a desired evolution is a crucial problem in almost all quantum technologies. We present a very general method for designing high-efficiency control sequences that are always fully compatible with experimental constraints on available interactions and their tunability. Our approach reduces in the end to finding control fields by solving a set of time-independent linear equations. We illustrate our method by applying it to a number of physically-relevant problems: the strong-driving limit of a two-level system, fast squeezing in a parametrically driven cavity, the leakage problem in transmon qubit gates, and the acceleration of SNAP gates in a qubit-cavity system.


title: Tuned Hybrid Non-Uniform Subdivision Surfaces with Optimal Convergence Rates

abstract: This paper presents an enhanced version of our previous work, hybrid non-uniform subdivision surfaces [19], to achieve optimal convergence rates in isogeometric analysis. We introduce a parameter $\lambda$ ($\frac{1}{4}<\lambda<1$) to control the rate of shrinkage of irregular regions, so the method is called tuned hybrid non-uniform subdivision (tHNUS). Our previous work corresponds to the case when $\lambda=\frac{1}{2}$. While introducing $\lambda$ in hybrid subdivision significantly complicates the theoretical proof of $G^1$ continuity around extraordinary vertices, reducing $\lambda$ can recover the optimal convergence rates when tuned hybrid subdivision functions are used as a basis in isogeometric analysis. From the geometric point of view, the tHNUS retains comparable shape quality as [19] under non-uniform parameterization. Its basis functions are refinable and the geometric mapping stays invariant during refinement. Moreover, we prove that a tuned hybrid subdivision surface is globally $G^1$-continuous. From the analysis point of view, tHNUS basis functions form a non-negative partition of unity, are globally linearly independent, and their spline spaces are nested. We numerically demonstrate that tHNUS basis functions can achieve optimal convergence rates for the Poisson's problem with non-uniform parameterization around extraordinary vertices.


title: Non-equilibrium dynamics of initially spherical vesicles in general flow

abstract: Many vesicles have a spherical resting shape and exposure to fluid flows induces an exchange between sub-optical area and visible (systematic) deformation, while the total area is conserved. The dynamics which controls the exchange between sub-optical and visible area depends on membrane properties such as bending rigidity and initial tension. Conversely, observation of these dynamics can be used to determine the membrane properties. The goal of this work is to create a general numerical model which accounts for the exchange between sub-optical and visible area. Unlike prior modeling efforts, the model does not pre-assume a shape type, such as nearly-spherical, or applied flow field, allowing the model to capture a wider variety of flow conditions. Based on implicit interface tracking and using a volume-preserving multiphase Navier-Stokes solver, the model is compared to several experimental results, showing excellent agreement. It is used to explore regimes not possible with previously published models, such as the variable viscosity case, and how these system properties influence experimentally measurable parameters such as deformation parameter. By creating a more generalized framework for the modeling of vesicles with sub-optical area, it will now be possible to make predictions on vesicle material properties from a wider variety of experimental results.


title: Magnetoelastic waves in thin films

abstract: This paper discusses the physics of magnetoelasticity and magnetoelastic waves as well as their mathematical description. Magnetoelastic waves occur as a result of strong coupling between spin waves and elastic waves in magnetostrictive ferromagnetic media. In a first part, the basic behavior of spin waves is reviewed in both bulk ferromagnets as well as in thin films. Next, elastic waves are discussed with a focus on thin films. Finally, the interactions between the elastic and magnetic domains are described and it is shown how this results in the formation of magnetoelastic waves. Based on the description of bulk magnetoelastic waves, a theory for magnetoelastic waves in thin films is developed and their dispersion relations are derived. It is shown that the behavior strongly depends on the geometry of the system, especially on the polarization of spin and elastic waves and the direction of the magnetization of the magnetostrictive ferromagnet.


title: Study of galaxy morphology and merging time of two interacting galaxies under different initial rotation and orientation configurations

abstract: Using the GADGET-2 N-body code, we make a study of the galaxy morphology and merging time due to two interacting galaxies (for the same types and different sizes and masses, $1:1$ and $1:10$ ratio masses) merging due to gravity interaction. This is done for different initial relative orientation and rotation of these galaxies (modes of interaction) but with the same relative bulge separation and the same relative initial velocities. It was found that the resulting galaxy morphology resemble many of the observed galaxies in our Universe, and that, in general, a binary galaxy system with 1:10 mass ratio has larger merging time than a binary galaxy system with 1:1 mass ratio. This difference is due to the different evolution of the masses during the interaction in both cases. For the case with a 1:10 mass ratio, the global mass maximum is located at the end evolution, meaning that the second galaxy increases its mass constantly. For the case with mass ratio 1:1, the global maximum is located around $t=0.35$ Gy, causing a reduction of the merging time.


title: Enabling Efficient and Flexible FPGA Virtualization for Deep Learning in the Cloud

abstract: FPGAs have shown great potential in providing low-latency and energy-efficient solutions for deep neural network (DNN) inference applications. Currently, the majority of FPGA-based DNN accelerators in the cloud run in a time-division multiplexing way for multiple users sharing a single FPGA, and require re-compilation with $\sim$100 s overhead. Such designs lead to poor isolation and heavy performance loss for multiple users, which are far away from providing efficient and flexible FPGA virtualization for neither public nor private cloud scenarios. To solve these problems, we introduce a novel virtualization framework for instruction architecture set (ISA) based on DNN accelerators by sharing a single FPGA. We enable the isolation by introducing a two-level instruction dispatch module and a multi-core based hardware resources pool. Such designs provide isolated and runtime-programmable hardware resources, further leading to performance isolation for multiple users. On the other hand, to overcome the heavy re-compilation overheads, we propose a tiling-based instruction frame package design and two-stage static-dynamic compilation. Only the light-weight runtime information is re-compiled with $\sim$1 ms overhead, thus the performance is guaranteed for the private cloud. Our extensive experimental results show that the proposed virtualization design achieves 1.07-1.69x and 1.88-3.12x throughput improvement over previous static designs using the single-core and the multi-core architectures, respectively.


title: Translationally cold trapped CCl+ reactions with acetylene (C2H2)

abstract: Ion-neutral chemical reactions are important in several areas of chemistry, including in some regions of the interstellar medium, planetary atmospheres, and comets. Reactions of CCl+ with C2H2 are measured and the main products include c-C3H2+ and l-C3H+, both relevant in extraterrestrial environments. Accurate branching ratios are obtained, which favor formation of c-C3H2+ over l-C3H+ by a factor of four. Measured rate constants are on the order of Langevin and complementary electronic structure calculations are used to aid in the interpretation of experimental results.


title: Real-time information retrieval from Identity cards

abstract: Information is frequently retrieved from valid personal ID cards by the authorised organisation to address different purposes. The successful information retrieval (IR) depends on the accuracy and timing process. A process which necessitates a long time to respond is frustrating for both sides in the exchange of data. This paper aims to propose a series of state-of-the-art methods for the journey of an Identification card (ID) from the scanning or capture phase to the point before Optical character recognition (OCR). The key factors for this proposal are the accuracy and speed of the process during the journey. The experimental results of this research prove that utilising the methods based on deep learning, such as Efficient and Accurate Scene Text (EAST) detector and Deep Neural Network (DNN) for face detection, instead of traditional methods increase the efficiency considerably.


title: Proximity-Induced Odd-Frequency Superconductivity in a Topological Insulator

abstract: At an interface between a topological insulator (TI) and a conventional superconductor (SC), superconductivity has been predicted to change dramatically and exhibit novel correlations. In particular, the induced superconductivity by an $s$-wave SC in a TI can develop an order parameter with a $p$-wave component. Here we present experimental evidence for an unexpected proximity-induced novel superconducting state in a thin layer of the prototypical TI, Bi$_2$Se$_3$, proximity coupled to Nb. From depth-resolved magnetic field measurements below the superconducting transition temperature of Nb, we observe a local enhancement of the magnetic field in Bi$_2$Se$_3$ that exceeds the externally applied field, thus supporting the existence of an intrinsic paramagnetic Meissner effect arising from an odd-frequency superconducting state. Our experimental results are complemented by theoretical calculations supporting the appearance of such a component at the interface which extends into the TI. This state is topologically distinct from the conventional Bardeen-Cooper-Schrieffer state it originates from. To the best of our knowledge, these findings present a first observation of bulk odd-frequency superconductivity in a TI. We thus reaffirm the potential of the TI-SC interface as a versatile platform to produce novel superconducting states.


title: Arbitrarily many independent observers can share the nonlocality of a single maximally entangled qubit pair

abstract: Alice and Bob each have half of a pair of entangled qubits. Bob measures his half and then passes his qubit to a second Bob who measures again and so on. The goal is to maximize the number of Bobs that can have an expected violation of the Clauser-Horne-Shimony-Holt (CHSH) Bell inequality with the single Alice. This scenario was introduced in [Phys. Rev. Lett. 114, 250401 (2015)] where the authors mentioned evidence that when the Bobs act independently and with unbiased inputs then at most two of them can expect to violate the CHSH inequality with Alice. Here we show that, contrary to this evidence, arbitrarily many independent Bobs can have an expected CHSH violation with the single Alice. Our proof is constructive and our measurement strategies can be generalized to work with a larger class of two-qubit states that includes all pure entangled two-qubit states. Since violation of a Bell inequality is necessary for device-independent tasks, our work represents a step towards an eventual understanding of the limitations on how much device-independent randomness can be robustly generated from a single pair of qubits.


title: Data-Driven Inference of Representation Invariants

abstract: A representation invariant is a property that holds of all values of abstract type produced by a module. Representation invariants play important roles in software engineering and program verification. In this paper, we develop a counterexample-driven algorithm for inferring a representation invariant that is sufficient to imply a desired specification for a module. The key novelty is a type-directed notion of visible inductiveness, which ensures that the algorithm makes progress toward its goal as it alternates between weakening and strengthening candidate invariants. The algorithm is parameterized by an example-based synthesis engine and a verifier, and we prove that it is sound and complete for first-order modules over finite types, assuming that the synthesizer and verifier are as well. We implement these ideas in a tool called Hanoi, which synthesizes representation invariants for recursive data types. Hanoi not only handles invariants for first-order code, but higher-order code as well. In its back end, Hanoi uses an enumerative synthesizer called Myth and an enumerative testing tool as a verifier. Because Hanoi uses testing for verification, it is not sound, though our empirical evaluation shows that it is successful on the benchmarks we investigated.


title: Reverse Faber-Krahn inequality for a truncated laplacian operator

abstract: In this paper we prove a reverse Faber-Krahn inequality for the principal eigenvalue $\mu_1(\Omega)$ of the fully nonlinear eigenvalue problem \[ \label{eq} \left\{\begin{array}{r c l l} -\lambda_N(D^2 u) & = & \mu u & \text{in }\Omega, \\ u & = & 0 & \text{on }\partial \Omega. \end{array}\right. \] Here $ \lambda_N(D^2 u)$ stands for the largest eigenvalue of the Hessian matrix of $u$. More precisely, we prove that, for an open, bounded, convex domain $\Omega \subset \mathbb{R}^N$, the inequality \[ \mu_1(\Omega) \leq \frac{\pi^2}{[\text{diam}(\Omega)]^2} = \mu_1(B_{\text{diam}(\Omega)/2}),\] where $\text{diam}(\Omega)$ is the diameter of $\Omega$, holds true. The inequality actually implies a stronger result, namely, the maximality of the ball under a diameter constraint. Furthermore, we discuss the minimization of $\mu_1(\Omega)$ under different kinds of constraints.


title: A Review of Multi-Objective Deep Learning Speech Denoising Methods

abstract: This paper presents a review of multi-objective deep learning methods that have been introduced in the literature for speech denoising. After stating an overview of conventional, single objective deep learning, and hybrid or combined conventional and deep learning methods, a review of the mathematical framework of the multi-objective deep learning methods for speech denoising is provided. A representative method from each speech denoising category, whose codes are publicly available, is selected and a comparison is carried out by considering the same public domain dataset and four widely used objective metrics. The comparison results indicate the effectiveness of the multi-objective method compared with the other methods, in particular when the signal-to-noise ratio is low. Possible future improvements that can be achieved are also mentioned.


title: Naked singularities in quadratic f(R) gravity

abstract: We find new static, spherically symmetric, and asymptotically flat vacuum solutions without horizon in Starobinsky's quadratic f(R) gravity. We systematically classify these solutions by an asymptotic analysis around the origin and find seven different integer Frobenius families. We numerically solve the exact equations of motion by a double-shooting method and specify boundary conditions by matching the numerical solution to the analytic solution of the linearised field equations in the weak field regime. We find that all integer Frobenius families can be connected to asymptotically flat solutions and trace out lines in the parameter space, allowing to ultimately relate all free parameters to the total mass at infinity.


title: Advanced Flow-Based Multilevel Hypergraph Partitioning

abstract: The balanced hypergraph partitioning problem is to partition a hypergraph into $k$ disjoint blocks of bounded size such that the sum of the number of blocks connected by each hyperedge is minimized. We present an improvement to the flow-based refinement framework of KaHyPar-MF, the current state-of-the-art multilevel $k$-way hypergraph partitioning algorithm for high-quality solutions. Our improvement is based on the recently proposed HyperFlowCutter algorithm for computing bipartitions of unweighted hypergraphs by solving a sequence of incremental maximum flow problems. Since vertices and hyperedges are aggregated during the coarsening phase, refinement algorithms employed in the multilevel setting must be able to handle both weighted hyperedges and weighted vertices -- even if the initial input hypergraph is unweighted. We therefore enhance HyperFlowCutter to handle weighted instances and propose a technique for computing maximum flows directly on weighted hypergraphs. We compare the performance of two configurations of our new algorithm with KaHyPar-MF and seven other partitioning algorithms on a comprehensive benchmark set with instances from application areas such as VLSI design, scientific computing, and SAT solving. Our first configuration, KaHyPar-HFC, computes slightly better solutions than KaHyPar-MF using significantly less running time. The second configuration, KaHyPar-HFC*, computes solutions of significantly better quality and is still slightly faster than KaHyPar-MF. Furthermore, in terms of solution quality, both configurations also outperform all other competing partitioners.


title: FFR V1.0: Fon-French Neural Machine Translation

abstract: Africa has the highest linguistic diversity in the world. On account of the importance of language to communication, and the importance of reliable, powerful and accurate machine translation models in modern inter-cultural communication, there have been (and still are) efforts to create state-of-the-art translation models for the many African languages. However, the low-resources, diacritical and tonal complexities of African languages are major issues facing African NLP today. The FFR is a major step towards creating a robust translation model from Fon, a very low-resource and tonal language, to French, for research and public use. In this paper, we describe our pilot project: the creation of a large growing corpora for Fon-to-French translations and our FFR v1.0 model, trained on this dataset. The dataset and model are made publicly available.


title: Network Structure and Collective Intelligence in the Diffusion of Innovation

abstract: When multiple innovations compete for adoption, historical chance leading to early advantage can generate lock-in effects that allow suboptimal innovations to succeed at the expense of superior alternatives. Research on the diffusion of innovafacetion has identified many possible sources of early advantage, but these mechanisms can benefit both optimal and suboptimal innovations. This paper moves beyond chance-as-explanation to identify structural principles that systematically impact the likelihood that the optimal strategy will spread. A formal model of innovation diffusion shows that the network structure of organizational relationships can systematically impact the likelihood that widely adopted innovations will be payoff optimal. Building on prior diffusion research, this paper focuses on the role of central actors i.e. well-connected people or firms. While contagion models of diffusion highlight the benefits of central actors for spreading innovations further and faster, the present analysis reveals a dark side to this influence: the mere presence of central actors in a network increases rates of adoption but also increases the likelihood of suboptimal outcomes. This effect, however, does not represent a speed-optimality tradeoff, as dense networks are both fast and optimal. This finding is consistent with related research showing that network centralization undermines collective intelligence.


title: Automorphism Groups of Endomorphisms of $\mathbb{P}^1 (\bar{\mathbb{F}}_p)$

abstract: For any algebraically closed field $K$ and any endomorphism $f$ of $\mathbb{P}^1(K)$ of degree at least 2, the automorphisms of $f$ are the M\"obius transformations that commute with $f$, and these form a finite subgroup of $\operatorname{PGL}_2(K)$. In the moduli space of complex dynamical systems, the locus of maps with nontrivial automorphisms has been studied in detail and there are techniques for constructing maps with prescribed automorphism groups that date back to Klein. We study the corresponding questions when $K$ is the algebraic closure $\bar{\mathbb{F}}_p$ of a finite field. We use the classification of finite subgroups of $\operatorname{PGL}_2(\bar{\mathbb{F}}_p)$ to show that every finite subgroup is realizable as an automorphism group. To construct examples, we use methods from modular invariant theory. Then, we calculate the locus of maps over $\bar{\mathbb{F}}_p$ of degree $2$ with nontrivial automorphisms, showing how the geometry and possible automorphism groups depend on the prime $p$.


title: Interference as an information-theoretic game

abstract: The double slit experiment provides a clear demarcation between classical and quantum theory, while multi-slit experiments demarcate quantum and higher-order interference theories. In this work we show that these experiments pertain to a broader class of processes, which can be formulated as information-processing tasks, providing a clear cut between classical, quantum and higher-order theories. The tasks involve two parties and communication between them with the goal of winning certain parity games. We show that the order of interference is in one-to-one correspondence with the parity order of these games. Furthermore, we prove the order of interference to be additive under composition of systems both in classical and quantum theory. The latter result can be used as a (semi)device-independent witness of the number of particles in the quantum setting. Finally, we extend our game formulation within the generalized probabilistic framework and prove that tomographic locality implies the additivity of the order of interference under composition. These results shed light on the operational meaning of the order of interference and can be important for the identification of the information-theoretic principles behind second-order interference in quantum theory.


title: A Decomposition Theorem for Aronszajn Lines

abstract: We show that under the proper forcing axiom the class of all Aronszajn lines behave like $\sigma$-scattered orders under the embeddability relation. In particular, we are able to show that the class of better quasi order labeled fragmented Aronszajn lines is itself a better quasi order. Moreover, we show that every better quasi order labeled Aronszajn line can be expressed as a finite sum of labeled types which are algebraically indecomposable. By encoding lines with finite labeled trees, we are also able to deduce a decomposition result, that for every Aronszajn line $L$ there is integer n such that for any finite colouring of $L$ there is subset $L^\prime$ of $L$ isomorphic to $L$ which uses no more than n colours.


title: The Mira-Titan Universe. III. Emulation of the Halo Mass Function

abstract: We construct an emulator for the halo mass function over group and cluster mass scales for a range of cosmologies, including the effects of dynamical dark energy and massive neutrinos. The emulator is based on the recently completed Mira-Titan Universe suite of cosmological $N$-body simulations. The main set of simulations spans 111 cosmological models with 2.1 Gpc boxes. We extract halo catalogs in the redshift range $z=[0.0, 2.0]$ and for masses $M_{200\mathrm{c}}\geq 10^{13}M_\odot/h$. The emulator covers an 8-dimensional hypercube spanned by {$\Omega_\mathrm{m}h^2$, $\Omega_\mathrm{b}h^2$, $\Omega_\nu h^2$, $\sigma_8$, $h$, $n_s$, $w_0$, $w_a$}; spatial flatness is assumed. We obtain smooth halo mass functions by fitting piecewise second-order polynomials to the halo catalogs and employ Gaussian process regression to construct the emulator while keeping track of the statistical noise in the input halo catalogs and uncertainties in the regression process. For redshifts $z\lesssim1$, the typical emulator precision is better than $2\%$ for $10^{13}-10^{14} M_\odot/h$ and $<10\%$ for $M\simeq 10^{15}M_\odot/h$. For comparison, fitting functions using the traditional universal form for the halo mass function can be biased at up to 30\% at $M\simeq 10^{14}M_\odot/h$ for $z=0$. Our emulator is publicly available at \url{https://github.com/SebastianBocquet/MiraTitanHMFemulator}.


title: STRIDES: Spectroscopic and photometric characterization of the environment and effects of mass along the line of sight to the gravitational lenses DES J0408-5354 and WGD 2038-4008

abstract: In time-delay cosmography, three of the key ingredients are 1) determining the velocity dispersion of the lensing galaxy, 2) identifying galaxies and groups along the line of sight with sufficient proximity and mass to be included in the mass model, and 3) estimating the external convergence $\kappa_\mathrm{ext}$ from less massive structures that are not included in the mass model. We present results on all three of these ingredients for two time-delay lensed quasar systems, DES J0408-5354 and WGD 2038-4008. We use the Gemini, Magellan and VLT telescopes to obtain spectra to both measure the stellar velocity dispersions of the main lensing galaxies and to identify the line-of-sight galaxies in these systems. Next, we identify 10 groups in DES J0408-5354 and 2 groups in WGD 2038-4008using a group-finding algorithm. We then identify the most significant galaxy and galaxy-group perturbers using the "flexion shift" criterion. We determine the probability distribution function of the external convergence $\kappa_\mathrm{ext}$ for both of these systems based on our spectroscopy and on the DES-only multiband wide-field observations. Using weighted galaxy counts, calibrated based on the Millennium Simulation, we find that DES J0408-5354 is located in a significantly underdense environment, leading to a tight (width $\sim3\%$), negative-value $\kappa_\mathrm{ext}$ distribution. On the other hand, WGD 2038-4008 is located in an environment of close to unit density, and its low source redshift results in a much tighter $\kappa_\mathrm{ext}$ of $\sim1\%$, as long as no external shear constraints are imposed.


title: Correlation between UV resilience and wavelength of narrow DIBs

abstract: Carriers of diffuse interstellar bands (DIBs) still need to be identified. In a recent paper, we reported a correlation between the DIB wavelength and the apparent UV resilience (or boost) of their carriers. We proposed that this might be an indication of the important role of conjugated elongated molecules among the DIB carriers. The aim of this paper is to further understand the origin of this correlation. The analysis of 509 optical DIBs on the lines of sight of HD 183143 and/or HD 204827 reported in the literature shows that this correlation mainly implies the 386 narrow DIBs with a band width < 1.1 A, which include most of the identified DIBs of the C2 and zeta families, while the majority of the 123 broader DIBs, including the identified sigma DIBs, do not display such a correlation. We present a possible origin of this correlation from very strong bands of large conjugated elongated molecules, such as carbon chains, polyacenes, or other catacondensed polycyclic aromatic hydrocarbons. The total amount of carbon contained in all the carriers of these narrow DIBs is a very small fraction of the interstellar carbon if their oscillator strengths are <~1. The amount of carbon locked in the carriers of the broader DIBs is higher, especially if their oscillator strengths are significantly weaker.


title: Floquet spectrum for anisotropic and tilted Dirac materials under linearly polarized light at all field intensities

abstract: The Floquet spectrum in an anisotropic tilted Dirac semimetal modulated by linearly polarized light is addressed through the solution of the time-dependent Schr\"odinger equation for the two-dimensional Dirac Hamiltonian via the Floquet theorem. The time-dependent wave functions and the quasienergy spectrum of the two-dimensional Dirac Hamiltonian under the normal incidence of linearly polarized waves are obtained for an arbitrarily intense electromagnetic radiation. We applied a set of unitary transformations to reduce the Schr\"odinger equation to an ordinary second-order differential Hill equation with complex coefficients. Through the stability analysis of this differential equation, the weak and strong field regimes are clearly distinguished in the quasi-spectrum. In the weak electric field regime, above a certain threshold given by the field parameters, the spectrum mostly resembles that of free electrons in graphene. Below this threshold, in the strong electric field regime, the spectrum abruptly becomes highly anisotropic and a gap opens up. As an example, we apply the results to the particular case of borophene.


title: Gaussian-Dirichlet Random Fields for Inference over High Dimensional Categorical Observations

abstract: We propose a generative model for the spatio-temporal distribution of high dimensional categorical observations. These are commonly produced by robots equipped with an imaging sensor such as a camera, paired with an image classifier, potentially producing observations over thousands of categories. The proposed approach combines the use of Dirichlet distributions to model sparse co-occurrence relations between the observed categories using a latent variable, and Gaussian processes to model the latent variable's spatio-temporal distribution. Experiments in this paper show that the resulting model is able to efficiently and accurately approximate the temporal distribution of high dimensional categorical measurements such as taxonomic observations of microscopic organisms in the ocean, even in unobserved (held out) locations, far from other samples. This work's primary motivation is to enable deployment of informative path planning techniques over high dimensional categorical fields, which until now have been limited to scalar or low dimensional vector observations.


title: Energy storage in magnetic textures driven by vorticity flow

abstract: An experimentally feasible energy-storage concept is formulated based on vorticity (hydro)dynamics within an easy-plane insulating magnet. The free energy, associated with the magnetic winding texture, is built up in a circular easy-plane magnetic structure by injecting a vorticity flow in the radial direction. The latter is accomplished by electrically induced spin-transfer torque, which pumps energy into the magnetic system in proportion to the vortex flux. The resultant magnetic metastable state with a finite winding number can be maintained for a long time because the process of its relaxation via phase slips is exponentially suppressed when the temperature is well below the Curie temperature. We propose to characterize the vorticity-current interaction underlying the energy-loading mechanism through its contribution to the effective electric inductance in the rf response. Our proposal may open an avenue for naturally powering spintronic circuits and nontraditional magnet-based neuromorphic networks.


title: Pedestrian Detection with Wearable Cameras for the Blind: A Two-way Perspective

abstract: Blind people have limited access to information about their surroundings, which is important for ensuring one's safety, managing social interactions, and identifying approaching pedestrians. With advances in computer vision, wearable cameras can provide equitable access to such information. However, the always-on nature of these assistive technologies poses privacy concerns for parties that may get recorded. We explore this tension from both perspectives, those of sighted passersby and blind users, taking into account camera visibility, in-person versus remote experience, and extracted visual information. We conduct two studies: an online survey with MTurkers (N=206) and an in-person experience study between pairs of blind (N=10) and sighted (N=40) participants, where blind participants wear a working prototype for pedestrian detection and pass by sighted participants. Our results suggest that both of the perspectives of users and bystanders and the several factors mentioned above need to be carefully considered to mitigate potential social tensions.


title: Robinson-Schensted correspondence for unit interval orders

abstract: The Stanley-Stembridge conjecture associates a symmetric function to each natural unit interval order $\mathcal P$. In this paper, we define relations \`a la Knuth on the symmetric group for each $\mathcal P$ and conjecture that the associated $\mathcal P$-Knuth equivalence classes are Schur-positive, refining theorems of Gasharov, Brosnan-Chow, and Guay-Paquet. The resulting equivalence graphs fit into the framework of D graphs studied by Assaf. Furthermore, we conjecture that the Schur expansion is given by column-readings of $\mathcal P$-tableaux that occur in the equivalence class. We prove these conjectures for $\mathcal P$ avoiding two specific suborders by introducing $\mathcal P$-analog of Robinson-Schensted insertion, giving an answer to a long standing question of Chow.


title: Quasiclassical theory of $C_4$-symmetric magnetic order in disordered multiband metals

abstract: Recent experimental studies performed in the normal state of iron-based superconductors have discovered the existence of the $C_4$-symmetric (tetragonal) itinerant magnetic state. This state can be described as a spin density wave with two distinct magnetic vectors ${\vec Q}_1$ and ${\vec Q}_2$. Given an itinerant nature of magnetism in iron-pnictides, we develop a quasiclassical theory of tetragonal magnetic order in disordered three-band metal with anisotropic band structure. Within our model we find that the $C_4$-symmetric magnetism competes with the $C_2$-symmetric state with a single ${\vec Q}$ magnetic structure vector. Our main results is that disorder promotes tetragonal magnetic state which is in agreement with earlier theoretical studies.


title: SaccadeNet: A Fast and Accurate Object Detector

abstract: Object detection is an essential step towards holistic scene understanding. Most existing object detection algorithms attend to certain object areas once and then predict the object locations. However, neuroscientists have revealed that humans do not look at the scene in fixed steadiness. Instead, human eyes move around, locating informative parts to understand the object location. This active perceiving movement process is called \textit{saccade}. %In this paper, Inspired by such mechanism, we propose a fast and accurate object detector called \textit{SaccadeNet}. It contains four main modules, the \cenam, the \coram, the \atm, and the \aggatt, which allows it to attend to different informative object keypoints, and predict object locations from coarse to fine. The \coram~is used only during training to extract more informative corner features which brings free-lunch performance boost. On the MS COCO dataset, we achieve the performance of 40.4\% mAP at 28 FPS and 30.5\% mAP at 118 FPS. Among all the real-time object detectors, %that can run faster than 25 FPS, our SaccadeNet achieves the best detection performance, which demonstrates the effectiveness of the proposed detection mechanism.


title: Quantifying deviations from separability in space-time functional processes

abstract: The estimation of covariance operators of spatio-temporal data is in many applications only computationally feasible under simplifying assumptions, such as separability of the covariance into strictly temporal and spatial factors.Powerful tests for this assumption have been proposed in the literature. However, as real world systems, such as climate data are notoriously inseparable, validating this assumption by statistical tests, seems inherently questionable. In this paper we present an alternative approach: By virtue of separability measures, we quantify how strongly the data's covariance operator diverges from a separable approximation. Confidence intervals localize these measures with statistical guarantees. This method provides users with a flexible tool, to weigh the computational gains of a separable model against the associated increase in bias. As separable approximations we consider the established methods of partial traces and partial products, and develop weak convergence principles for the corresponding estimators. Moreover, we also prove such results for estimators of optimal, separable approximations, which are arguably of most interest in applications. In particular we present for the first time statistical inference for this object, which has been confined to estimation previously. Besides confidence intervals, our results encompass tests for approximate separability. All methods proposed in this paper are free of nuisance parameters and do neither require computationally expensive resampling procedures nor the estimation of nuisance parameters. A simulation study underlines the advantages of our approach and its applicability is demonstrated by the investigation of German annual temperature data.


title: Gryffin: An algorithm for Bayesian optimization of categorical variables informed by expert knowledge

abstract: Designing functional molecules and advanced materials requires complex design choices: tuning continuous process parameters such as temperatures or flow rates, while simultaneously selecting catalysts or solvents. To date, the development of data-driven experiment planning strategies for autonomous experimentation has largely focused on continuous process parameters despite the urge to devise efficient strategies for the selection of categorical variables. Here, we introduce Gryffin, a general purpose optimization framework for the autonomous selection of categorical variables driven by expert knowledge. Gryffin augments Bayesian optimization based on kernel density estimation with smooth approximations to categorical distributions. Leveraging domain knowledge in the form of physicochemical descriptors, Gryffin can significantly accelerate the search for promising molecules and materials. Gryffin can further highlight relevant correlations between the provided descriptors to inspire physical insights and foster scientific intuition. In addition to comprehensive benchmarks, we demonstrate the capabilities and performance of Gryffin on three examples in materials science and chemistry: (i) the discovery of non-fullerene acceptors for organic solar cells, (ii) the design of hybrid organic-inorganic perovskites for light harvesting, and (iii) the identification of ligands and process parameters for Suzuki-Miyaura reactions. Our results suggest that Gryffin, in its simplest form, is competitive with state-of-the-art categorical optimization algorithms. However, when leveraging domain knowledge provided via descriptors, Gryffin outperforms other approaches while simultaneously refining this domain knowledge to promote scientific understanding.


title: Going in circles is the way forward: the role of recurrence in visual inference

abstract: Biological visual systems exhibit abundant recurrent connectivity. State-of-the-art neural network models for visual recognition, by contrast, rely heavily or exclusively on feedforward computation. Any finite-time recurrent neural network (RNN) can be unrolled along time to yield an equivalent feedforward neural network (FNN). This important insight suggests that computational neuroscientists may not need to engage recurrent computation, and that computer-vision engineers may be limiting themselves to a special case of FNN if they build recurrent models. Here we argue, to the contrary, that FNNs are a special case of RNNs and that computational neuroscientists and engineers should engage recurrence to understand how brains and machines can (1) achieve greater and more flexible computational depth, (2) compress complex computations into limited hardware, (3) integrate priors and priorities into visual inference through expectation and attention, (4) exploit sequential dependencies in their data for better inference and prediction, and (5) leverage the power of iterative computation.


title: A mathematical discussion of Pons Viver's implementation of L\"owdin's spin projection operator

abstract: Recently the molecular electronic structure theories for efficiently treating static (or strong) correlation in a black-box manner have attracted much attention. In these theories, a spin projection operator is used to recover the spin symmetry of a broken-symmetry Slater determinant. Very recently, Pons Viver proposed the practical and exact implementation of L\"owdin's spin projection operator [Int. J. Quantum Chem. 119, e25770 (2019)]. In the present study, we attempt to supply mathematical proofs to Pons Viver's proposals and show a condition for establishing Pons Viver's implementation. Moreover, we explicitly derive the (spin projected) extended Hartree-Fock (EHF) equations on the basis of the model of common orbitals (i.e., closed-shell orbitals used in the restricted open-shell Hartree-Fock (ROHF) method), which was combined by Pons Viver with the EHF method.


title: Charged multi-hadron systems in lattice QCD+QED

abstract: Systems with the quantum numbers of up to twelve charged and neutral pseudoscalar mesons, as well as one-, two-, and three-nucleon systems, are studied using dynamical lattice quantum chromodynamics and quantum electrodynamics (QCD+QED) calculations and effective field theory. QED effects on hadronic interactions are determined by comparing systems of charged and neutral hadrons after tuning the quark masses to remove strong isospin breaking effects. A non-relativistic effective field theory, which perturbatively includes finite-volume Coulomb effects, is analyzed for systems of multiple charged hadrons and found to accurately reproduce the lattice QCD+QED results. QED effects on charged multi-hadron systems beyond Coulomb photon exchange are determined by comparing the two- and three-body interaction parameters extracted from the lattice QCD+QED results for charged and neutral multi-hadron systems.


title: Tiling large areas with luminescent solar concentrators: geometry effects

abstract: Luminescence solar concentrators act as semi-transparent photovoltaic cells and can be applied to large surface areas in modern urban environment. In this paper their optical efficiencies were analytically derived for different unit shapes as simple, integral-free expressions. Perimeter length increase for the same enclosed area requires more photodetectors, yet reduces optical losses per unit. An exact interplay between these two parameters is presented, which can be directly applied to the analysis of tiling for large areas. An explicit expression for the critical size of an LSC unit, above which its inner part becomes inactive, has been obtained. Fabrication of devices larger than specified by this limit should be discouraged from efficiency considerations.


title: On the Emerging Area of Biocybersecurity and Relevant Considerations

abstract: Biocybersecurity is a novel space for the 21st century that meets our innovations in biotechnology and computing head on. Within this space, many considerations are open for and demand consideration as groups endeavor to develop products and policies that adequately ensure asset management and protection. Herein, simplified and brief exploration is given followed by some surface discussion of impacts. These impacts concern the end user, ethical and legal considerations, international proceedings, business, and limitations. It is hoped that this will be helpful in future considerations towards biocybersecurity policy developments and implementations. Notice: This article has been queued for publication in the Proceedings of the 2020 Future of Information and Communication Conference (FICC)


title: Machine learning as a model for cultural learning: Teaching an algorithm what it means to be fat

abstract: As we navigate our cultural environment, we learn cultural biases, like those around gender, social class, health, and body weight. It is unclear, however, exactly how public culture becomes private culture. In this paper, we provide a theoretical account of such cultural learning. We propose that neural word embeddings provide a parsimonious and cognitively plausible model of the representations learned from natural language. Using neural word embeddings, we extract cultural schemata about body weight from New York Times articles. We identify several cultural schemata that link obesity to gender, immorality, poor health, and low socioeconomic class. Such schemata may be subtly but pervasively activated in public culture; thus, language can chronically reproduce biases. Our findings reinforce ongoing concerns that machine learning can also encode, and reproduce, harmful human biases.


title: Miniature Robot Path Planning for Bridge Inspection: Min-Max Cycle Cover-Based Approach

abstract: We study the problem of planning the deployments of a group of mobile robots. While the problem and formulation can be used for many different problems, here we use a bridge inspection as the motivating application for the purpose of exposition. The robots are initially stationed at a set of depots placed throughout the bridge. Each robot is then assigned a set of sites on the bridge to inspect and, upon completion, must return to the same depot where it is stored. The problem of robot planning is formulated as a rooted min-max cycle cover problem, in which the vertex set consists of the sites to be inspected and robot depots, and the weight of an edge captures either (i) the amount of time needed to travel from one end vertex to the other vertex or (ii) the necessary energy expenditure for the travel. In the first case, the objective function is the total inspection time, whereas in the latter case, it is the maximum energy expenditure among all deployed robots. We propose a novel algorithm with approximation ratio of $5 + \epsilon$, where $0<\epsilon<1$. In addition, the computational complexity of the proposed algorithm is shown to be $O\big( n^2+2^{m-1} n \log(n+k) \big)$, where $n$ is the number of vertices, and $m$ is the number of depots.


title: Using ray-traversal for 3D particle matching in the context of particle tracking velocimetry in fluid mechanics

abstract: An innovative method based on the traversal of rays, originating from detected particles, through a three-dimensional grid of voxels is presented. The methodology has as main advantage that the outcome of the method is independent of the order of the input; the order of the cameras and the order of the rays presented as input to the algorithm does not influence the outcome. The algorithm finds matches in decreasing value of match quality, ensuring that globally best matches are matched before worse matches. The time complexity of the algorithm is found to scale efficiently with the number of cameras and particles. A variety of show-cases are given to exemplify the algorithm for different geometries and different number of cameras. The method is designed for the tracking of tracer or inertial particles in fluid mechanics, for which the particle size generally ranges from O({\mu}m)--O(cm). The method, however, does not impose a size limit on the particles.


title: Results on Total and Elastic Cross Sections in Proton-Proton Collisions at $\sqrt{s} = 200$ GeV

abstract: We report results on the total and elastic cross sections in proton-proton collisions at $\sqrt{s}=200$ GeV obtained with the Roman Pot setup of the STAR experiment at the Relativistic Heavy Ion Collider (RHIC). The elastic differential cross section was measured in the squared four-momentum transfer range $0.045 \leq -t \leq 0.135$ GeV$^2$. The value of the exponential slope parameter $B$ of the elastic differential cross section $d\sigma/dt \sim e^{-Bt}$ in the measured $-t$ range was found to be $B = 14.32 \pm 0.09 (stat.)^{\scriptstyle +0.13}_{\scriptstyle -0.28} (syst.)$ GeV$^{-2}$. The total cross section $\sigma_{tot}$, obtained from extrapolation of the $d\sigma/dt$ to the optical point at $-t = 0$, is $\sigma_{tot} = 54.67 \pm 0.21 (stat.) ^{\scriptstyle +1.28}_{\scriptstyle -1.38} (syst.)$ mb. We also present the values of the elastic cross section $\sigma_{el} = 10.85 \pm 0.03 (stat.) ^{\scriptstyle +0.49}_{\scriptstyle -0.41}(syst.)$ mb, the elastic cross section integrated within the STAR $t$-range $\sigma^{det}_{el} = 4.05 \pm 0.01 (stat.) ^{\scriptstyle+0.18}_{\scriptstyle -0.17}(syst.)$ mb, and the inelastic cross section $\sigma_{inel} = 43.82 \pm 0.21 (stat.) ^{\scriptstyle +1.37}_{\scriptstyle -1.44} (syst.)$ mb. The results are compared with the world data.


title: Cycle Text-To-Image GAN with BERT

abstract: We explore novel approaches to the task of image generation from their respective captions, building on state-of-the-art GAN architectures. Particularly, we baseline our models with the Attention-based GANs that learn attention mappings from words to image features. To better capture the features of the descriptions, we then built a novel cyclic design that learns an inverse function to maps the image back to original caption. Additionally, we incorporated recently developed BERT pretrained word embeddings as our initial text featurizer and observe a noticeable improvement in qualitative and quantitative performance compared to the Attention GAN baseline.


title: Sixteen overlooked open clusters in the fourth Galactic quadrant. A combined analysis of UBVI photometry and Gaia DR2 with ASteCA

abstract: Aims: This paper has two main objectives: (1) To determine the intrinsic properties of 16 faint and mostly unstudied open clusters in the poorly known sector of the Galaxy at 270$^\circ-$300$^\circ$, to probe the Milky Way structure in future investigations. (2) To address previously reported systematics in Gaia DR2 parallaxes by comparing the cluster distances derived from photometry with those derived from parallaxes. Methods: Deep UBVI photometry of 16 open clusters was carried out. Observations were reduced and analyzed in an automaticway using the ASteCA package to get individual distances, reddening, masses, ages and metallicities. Photometric distances were compared to those obtained from a Bayesian analysis of Gaia DR2 parallaxes. Results: Ten out of the 16 clusters are true or highly probable open clusters. Two of them are quite young and follow the trace of the Carina Arm and the already detected warp. The rest of the clusters are placed in the interarm zone between the Perseus and Carina Arms as expected for older objects. We found that the cluster van den Berg-Hagen 85 is 7.5$\times$10$^9$ yrs old becoming then one of the oldest open cluster detected in our Galaxy so far. The relationship of these ten clusters with the Galaxy structure in the solar neighborhood is discussed. The comparison of distances from photometry and parallaxes data, in turn, reveals a variable level of disagreement. Conclusions: Various zero point corrections for Gaia DR2 parallax data recently reported were considered for a comparison between photometric and parallax based distances. The results tend to improve with some of these corrections. Photometric distance analysis suggest an average correction of $\sim$+0.026 mas (to be added to the parallaxes). The correction may have a more intricate distance dependency, but addressing that level of detail will require a larger cluster sample.


title: Integrating Crowdsourcing and Active Learning for Classification of Work-Life Events from Tweets

abstract: Social media, especially Twitter, is being increasingly used for research with predictive analytics. In social media studies, natural language processing (NLP) techniques are used in conjunction with expert-based, manual and qualitative analyses. However, social media data are unstructured and must undergo complex manipulation for research use. The manual annotation is the most resource and time-consuming process that multiple expert raters have to reach consensus on every item, but is essential to create gold-standard datasets for training NLP-based machine learning classifiers. To reduce the burden of the manual annotation, yet maintaining its reliability, we devised a crowdsourcing pipeline combined with active learning strategies. We demonstrated its effectiveness through a case study that identifies job loss events from individual tweets. We used Amazon Mechanical Turk platform to recruit annotators from the Internet and designed a number of quality control measures to assure annotation accuracy. We evaluated 4 different active learning strategies (i.e., least confident, entropy, vote entropy, and Kullback-Leibler divergence). The active learning strategies aim at reducing the number of tweets needed to reach a desired performance of automated classification. Results show that crowdsourcing is useful to create high-quality annotations and active learning helps in reducing the number of required tweets, although there was no substantial difference among the strategies tested.


title: MetNet: A Neural Weather Model for Precipitation Forecasting

abstract: Weather forecasting is a long standing scientific challenge with direct social and economic impact. The task is suitable for deep neural networks due to vast amounts of continuously collected data and a rich spatial and temporal structure that presents long range dependencies. We introduce MetNet, a neural network that forecasts precipitation up to 8 hours into the future at the high spatial resolution of 1 km$^2$ and at the temporal resolution of 2 minutes with a latency in the order of seconds. MetNet takes as input radar and satellite data and forecast lead time and produces a probabilistic precipitation map. The architecture uses axial self-attention to aggregate the global context from a large input patch corresponding to a million square kilometers. We evaluate the performance of MetNet at various precipitation thresholds and find that MetNet outperforms Numerical Weather Prediction at forecasts of up to 7 to 8 hours on the scale of the continental United States.


title: Scalable Deployment of AI Time-series Models for IoT

abstract: IBM Research Castor, a cloud-native system for managing and deploying large numbers of AI time-series models in IoT applications, is described. Modelling code templates, in Python and R, following a typical machine-learning workflow are supported. A knowledge-based approach to managing model and time-series data allows the use of general semantic concepts for expressing feature engineering tasks. Model templates can be programmatically deployed against specific instances of semantic concepts, thus supporting model reuse and automated replication as the IoT application grows. Deployed models are automatically executed in parallel leveraging a serverless cloud computing framework. The complete history of trained model versions and rolling-horizon predictions is persisted, thus enabling full model lineage and traceability. Results from deployments in real-world smart-grid live forecasting applications are reported. Scalability of executing up to tens of thousands of AI modelling tasks is also evaluated.


title: Graphene-ionic liquid interfacial potential drop from DFT-MD simulations

abstract: Ionic liquids (IL) are promising electrolytes for electrochemical applications due to their remarkable stability and high charge density. Molecular dynamics simulations are essential for better understanding the complex phenomena occurring at the electrode-IL interface. In this work, we have studied the interface between graphene and 1-ethyl-3-methyl-imidazolium tetrafluoroborate IL, using density functional theory-based molecular dynamics simulations at variable surface charge densities. We have disassembled the electrical double layer potential drop into two main components: one involving atomic charges and the other - dipoles. The latter component arises due to the electronic polarisation of the surface and is related to concepts hotly debated in the literature, such as the Thomas-Fermi screening length, effective surface charge plane, and quantum capacitance.


title: Coronavirus Geographic Dissemination at Chicago and its Potential Proximity to Public Commuter Rail

abstract: The community spread of coronavirus at great Chicago area has severely threatened the residents health, family and normal activities. CDC daily updates on infected cases on County level are not satisfying to address publics concern on virus spread. On March 20th, NBC5 published case information of 435 coronavirus infections. The data is relative comprehensive and of high value for understanding on the virus spread patterns at Chicago. Data engineering, natural language processing and Google map technology are applied to organize the data and retrieve geographic information of the virus. The analysis shows community spread in Chicago areas has a potential proximity relation with public commuter rail. Residents nearby major public commuter rails need limit outdoor activities during the outbreak and even the post-peak time.


title: Metrics and Optimization of Internal Poses for Highly Redundant Truss-Like Serialized Parallel Manipulators

abstract: This paper presents a kinematic definition of a serialized Stewart platform designed for autonomous in-space assembly called an Assembler. The Assemblers architecture describes problems inherent to the inverse kinematics of over-actuated mixed kinematic systems. This paper also presents a methodology for optimizing poses. In order to accomplish this with the Assembler system, an algorithm for finding a feasible solution to its inverse kinematics was developed with a wrapper for a nonlinear optimization algorithm designed to minimize the magnitude of forces incurred by each actuator. A simulated version of an Assembler was placed into a number of representative poses, and the positions were optimized. The results of these optimizations are discussed in terms of actuator forces, reachability of the platform, and applicability to high-payload structure assembly capabilities.


title: Knowledge Graph Alignment using String Edit Distance

abstract: In this work, we propose a novel knowledge graph alignment technique based upon string edit distance that exploits the type information between entities and can find similarity between relations of any arity


title: Scaling Confirmation of the Thermodynamic Dislocation Theory

abstract: We show that the thermodynamic dislocation theory (TDT) predicts a scaling relation between stresses, strain rates, and temperatures for steady-state deformations of crystalline solids, and that this relation is accurately obeyed by a wide range of experimental data for both aluminum and copper. Unlike conventional phenomenological dislocation theories, the TDT is based on the second law of thermodynamics. Its success implies that descriptions of solid deformation that are not based on the statistical mechanics of nonequilibrium processes cannot be relied upon to be predictive. Thus there is an urgent need -- and a new opportunity -- to revitalize this central part of materials physics.


title: Deep Convolutional Recurrent Autoencoders for Flow Field Prediction

abstract: In this paper, an end-to-end nonlinear model reduction methodology is presented based on the convolutional recurrent autoencoder networks. The methodology is developed in the context of the overall data-driven reduced-order model framework proposed in the paper. The basic idea behind the methodology is to obtain the low dimensional representations via convolutional neural networks and evolve these low dimensional features via recurrent neural networks in the time domain. The high dimensional representations are constructed from the evolved low dimensional features via transpose convolutional neural networks. With an unsupervised training strategy, the model serves as an end to end tool which can evolve the flow state of the nonlinear dynamical system. The convolutional recurrent autoencoder network model is applied to the problem of flow past bluff bodies for the first time. To demonstrate the effectiveness of the methodology, two canonical problems namely the flow past a plain cylinder and the flow past side-by-side cylinders are explored in this paper. Pressure and velocity fields of the unsteady flow are predicted in future via the convolutional recurrent autoencoder model. The performance of the model is satisfactory for both the problems. Specifically, the multiscale nature and the gap flow dynamics of the side-by-side cylinders are captured by the proposed data-driven model reduction methodology. The error metrics, the normalized squared error, and the normalized reconstruction error are considered for the assessment of the data-driven framework.


title: A Flexible Job Shop Scheduling Representation of the Autonomous In-Space Assembly Task Assignment Problem

abstract: As in-space exploration increases, autonomous systems will play a vital role in building the necessary facilities to support exploration. To this end, an autonomous system must be able to assign tasks in a scheme that efficiently completes all of the jobs in the desired project. This research proposes a flexible job shop problem (FJSP) representation to characterize an autonomous assembly project and then proposes both a mixed integer programming (MIP) solution formulation and a reinforcement learning (RL) solution formulation. The MIP formulation encodes all of the constraints and interjob dynamics a priori and was able to solve for the optimal solution to minimize the makespan. The RL formulation did not converge to an optimal solution but did successfully learn implicitly interjob dynamics through interaction with the reward function. Future work will include developing a solution formulation that utilizes the strengths of both proposed solution methods to handle scaling in size and complexity.


title: Dissipative generation of pure steady states and a gambler's ruin problem

abstract: We consider an open quantum system, with dissipation applied only to a part of its degrees of freedom, evolving via a quantum Markov dynamics. We demonstrate that, in the Zeno regime of large dissipation, the relaxation of the quantum system towards a pure quantum state is linked to the evolution of a classical Markov process towards a single absorbing state. The rates of the associated classical Markov process are determined by the original quantum dynamics. Extension of this correspondence to absorbing states with internal structure allows us to establish a general criterion for having a Zeno-limit nonequilibrium stationary state of arbitrary finite rank. An application of this criterion is illustrated in the case of an open XXZ spin-1/2 chain dissipatively coupled at its edges to baths with fixed and different polarizations. For this system, we find exact nonequilibrium steady-state solutions of ranks 1 and 2.

categories: quant-ph cond-mat.stat-mech math-ph math.MP


title: Modeling the early evolution of the COVID-19 in Brazil: results from a Susceptible-Infectious-Quarantined-Recovered (SIQR) model

abstract: The world evolution of the Severe acute respiratory syndrome coronavirus 2 (SARS-Cov2 or simply COVID-19) led the World Health Organization to declare it a pandemic. The disease appeared in China in December 2019, and it has spread fast around the world, specially in european countries like Italy and Spain. The first reported case in Brazil was recorded in February 26, and after that the number of cases growed fast. In order to slow down the initial growth of the disease through the country, confirmed positive cases were isolated to not transmit the disease. To better understand the early evolution of COVID-19 in Brazil, we apply a Susceptible-Infectious-Quarantined-Recovered (SIQR) model to the analysis of data from the Brazilian Department of Health, obtained from February 26, 2020 through March 25, 2020. Based on analyical and numerical results, as well on the data, the basic reproduction number is estimated to $R_{0}=5.25$. In addition, we estimate that the ratio unidentified infectious individuals and confirmed cases at the beginning of the epidemic is about $10$, in agreement with previous studies. We also estimated the epidemic doubling time to be $2.72$ days.


title: Q-Learning in Regularized Mean-field Games

abstract: In this paper, we introduce a regularized mean-field game and study learning of this game under an infinite-horizon discounted reward function. Regularization is introduced by adding a strongly concave regularization function to the one-stage reward function in the classical mean-field game model. We establish a value iteration based learning algorithm to this regularized mean-field game using fitted Q-learning. The regularization term in general makes reinforcement learning algorithm more robust to the system components. Moreover, it enables us to establish error analysis of the learning algorithm without imposing restrictive convexity assumptions on the system components, which are needed in the absence of a regularization term.


title: Shouts and Murmurs: Combining Individual Gravitational-Wave Sources with the Stochastic Background to Measure the History of Binary Black Hole Mergers

abstract: One of the goals of gravitational-wave astronomy is to quantify the evolution of the compact binary merger rate with redshift. The redshift distribution of black hole mergers would offer considerable information about their evolutionary history, including their progenitor formation rate, the dependence of black hole formation on stellar metallicity, and the time delay distribution between formation and merger. Efforts to measure the binary redshift distribution are currently limited, however, by the detection range of existing instruments, which can individually resolve compact binary merger events only out to $z\lesssim1$. We present a novel strategy with which to measure the redshift distribution of binary black hole mergers well beyond the detection range of current instruments. By synthesizing direct detections of individually resolved mergers with \textit{indirect} searches for the stochastic gravitational-wave background due to unresolved distant sources, we can glean information about the peak redshift, $z_p$, at which the binary black hole merger rate attains its maximum, even when this redshift is beyond the detection horizon. Using data from Advanced LIGO and Virgo's first and second observing runs, we employ this strategy to place joint constraints on $z_p$ and the slope $\alpha$ with which the binary merger rate increases at low redshifts, ruling out merger rates that grow faster than $\alpha \gtrsim 7$ and peak beyond $z_p \gtrsim 1.5$. Looking ahead, we project that approximately one year of observation with design-sensitivity Advanced LIGO will further break remaining degeneracies, enabling a direct measurement of the peak redshift of the binary black hole merger history.


title: Predicting the benefit of wake steering on the annual energy production of a wind farm using large eddy simulations and Gaussian process regression

abstract: In recent years, wake steering has been established as a promising method to increase the energy yield of a wind farm. Current practice in estimating the benefit of wake steering on the annual energy production (AEP) consists of evaluating the wind farm with simplified surrogate models, casting a large uncertainty on the estimated benefit. This paper presents a framework for determining the benefit of wake steering on the AEP, incorporating simulation results from a surrogate model and large eddy simulations in order to reduce the uncertainty. Furthermore, a time-varying wind direction is considered for a better representation of the ambient conditions at the real wind farm site. Gaussian process regression is used to combine the two data sets into a single improved model of the energy gain. This model estimates a 0.60% gain in AEP for the considered wind farm, which is a 76% increase compared to the estimate of the surrogate model.


title: Not All Features Are Equal: Discovering Essential Features for Preserving Prediction Privacy

abstract: When receiving machine learning services from the cloud, the provider does not need to receive all features; in fact, only a subset of the features are necessary for the target prediction task. Discerning this subset is the key problem of this work. We formulate this problem as a gradient-based perturbation maximization method that discovers this subset in the input feature space with respect to the functionality of the prediction model used by the provider. After identifying the subset, our framework, Cloak, suppresses the rest of the features using utility-preserving constant values that are discovered through a separate gradient-based optimization process. We show that Cloak does not necessarily require collaboration from the service provider beyond its normal service, and can be applied in scenarios where we only have black-box access to the service provider's model. We theoretically guarantee that Cloak's optimizations reduce the upper bound of the Mutual Information (MI) between the data and the sifted representations that are sent out. Experimental results show that Cloak reduces the mutual information between the input and the sifted representations by 85.01% with only a negligible reduction in utility (1.42%). In addition, we show that Cloak greatly diminishes adversaries' ability to learn and infer non-conducive features.


title: Discovery of new members of the nearby young stellar association in Cepheus

abstract: Young field stars are hardly distinguishable from older ones because their space motion rapidly mixes them with the stellar population of the Galactic plane. Nevertheless, a careful target selection allows for young stars to be spotted throughout the sky. We aim to identify additional sources associated with the four young comoving stars that we discovered towards the CO Cepheus void and to provide a comprehensive view of the Cepheus association. Based on multivariate analysis methods, we have built an extended sample of 193 young star candidates, which are the optical and infrared counterparts of ROSAT All-Sky Survey and XMM-Newton X-ray sources. From optical spectroscopic observations, we measured their radial velocity with the cross-correlation technique. We derived their atmospheric parameters and projected rotational velocity with the code ROTFIT. We applied the subtraction of inactive templates to measure the lithium equivalent width, from which we infer their lithium abundance and age. Finally, we studied their kinematics using the second Gaia data release. Our sample is mainly composed of young or active stars and multiple systems. We identify two distinct populations of young stars that are spatially and kinematically separated. Those with an age between 100 and 300 Myr are mostly projected towards the Galactic plane. In contrast, 23 of the 37 sources younger than 30 Myr are located in the CO Cepheus void, and 21 of them belong to the stellar kinematic group that we previously reported in this sky area. We report a total of 32 bona fide members and nine candidates for this nearby (distance = 157$\pm$10 pc) young (age = 10-20 Myr) stellar association. According to the spatial distribution of its members, the original cluster is already dispersed and partially mixed with the local population of the Galactic plane.


title: Data Integration by combining big data and survey sample data for finite population inference

abstract: The statistical challenges in using big data for making valid statistical inference in the finite population have been well documented in literature. These challenges are due primarily to statistical bias arising from under-coverage in the big data source to represent the population of interest and measurement errors in the variables available in the data set. By stratifying the population into a big data stratum and a missing data stratum, we can estimate the missing data stratum by using a fully responding probability sample, and hence the population as a whole by using a data integration estimator. By expressing the data integration estimator as a regression estimator, we can handle measurement errors in the variables in big data and also in the probability sample. We also propose a fully nonparametric classification method for identifying the overlapping units and develop a bias-corrected data integration estimator under misclassification errors. Finally, we develop a two-step regression data integration estimator to deal with measurement errors in the probability sample. An advantage of the approach advocated in this paper is that we do not have to make unrealistic missing-at-random assumptions for the methods to work. The proposed method is applied to the real data example using 2015-16 Australian Agricultural Census data.


title: Sobolev inequalities with jointly concave weights on convex cones

abstract: Using optimal mass transport arguments, we prove weighted Sobolev inequalities of the form \[\left(\int_E |u(x)|^q\,\omega(x) \,dx\right)^{1/q}\leq K_0\,\left(\int_E |\nabla u(x)|^p\,\sigma(x)\,dx\right)^{1/p},\ \ u\in C_0^\infty(\mathbb R^n),\ \ \ \ \ \ {\rm (WSI)}\] where $p\geq 1$ and $q>0$ is the corresponding Sobolev critical exponent. Here $E\subseteq \mathbb R^n$ is an open convex cone, and $\omega,\sigma:E\to (0,\infty)$ are two homogeneous weights verifying a general concavity-type structural condition. The constant $K_0= K_0(n, p, q, \omega, \sigma) >0$ is given by an explicit formula. Under mild regularity assumptions on the weights, we also prove that $K_0$ is optimal in (WSI) if and only if $\omega$ and $\sigma$ are equal up to a multiplicative factor. Several previously known results, including the cases for monomials and radial weights, are covered by our statement. Further examples and applications to PDEs are also provided.


title: Molecular Simulation of Covalent Bond Dynamics in Liquid Silicon

abstract: Many atomic liquids can form transient covalent bonds reminiscent of those in the corresponding solid states. These directional interactions dictate many important properties of the liquid state, necessitating a quantitative, atomic-scale understanding of bonding in these complex systems. A prototypical example is liquid silicon, wherein transient covalent bonds give rise to local tetrahedral order and consequent non-trivial effects on liquid state thermodynamics and dynamics. To further understand covalent bonding in liquid silicon, and similar liquids, we present an ab initio simulation-based approach for quantifying the structure and dynamics of covalent bonds in condensed phases. Through the examination of structural correlations among silicon nuclei and maximally localized Wannier function centers, we develop a geometric criterion for covalent bonds in liquid Si. We use this to monitor the dynamics of transient covalent bonding in the liquid state and estimate a covalent bond lifetime. We compare covalent bond dynamics to other processes in liquid Si and similar liquids and suggest experiments to measure the covalent bond lifetime.


title: Learning To Solve Differential Equations Across Initial Conditions

abstract: Recently, there has been a lot of interest in using neural networks for solving partial differential equations. A number of neural network-based partial differential equation solvers have been formulated which provide performances equivalent, and in some cases even superior, to classical solvers. However, these neural solvers, in general, need to be retrained each time the initial conditions or the domain of the partial differential equation changes. In this work, we posit the problem of approximating the solution of a fixed partial differential equation for any arbitrary initial conditions as learning a conditional probability distribution. We demonstrate the utility of our method on Burger's Equation.


title: Traffic assignment models. Numerical aspects

abstract: In this book we describe BMW traffic assignment model and Nesterov-dePalma model. We consider Entropy model for demand matrix. Based on this models we build multi-stage traffic assignment models. The equilibrium in such models can be found from convex-concave saddle-point problem. We show how to solve this problem by using special combination of universal gradient method and Sinkhorn's algorithm.


title: $\ell$-torsion bounds for the class group of number fields with an $\ell$-group as Galois group

abstract: We describe the relations among the $\ell$-torsion conjecture, a conjecture of Malle giving an upper bound for the number of extensions, and the discriminant multiplicity conjecture. We prove that the latter two conjectures are equivalent in some sense. Altogether, the three conjectures are equivalent for the class of solvable groups. We then prove the $\ell$-torsion conjecture for $\ell$-groups and the other two conjectures for nilpotent groups.


title: Zero-shot and few-shot time series forecasting with ordinal regression recurrent neural networks

abstract: Recurrent neural networks (RNNs) are state-of-the-art in several sequential learning tasks, but they often require considerable amounts of data to generalise well. For many time series forecasting (TSF) tasks, only a few dozens of observations may be available at training time, which restricts use of this class of models. We propose a novel RNN-based model that directly addresses this problem by learning a shared feature embedding over the space of many quantised time series. We show how this enables our RNN framework to accurately and reliably forecast unseen time series, even when there is little to no training data available.


title: Using constraint structure and an improved object detection network to detect the 12^{th} Vertebra from CT images with a limited field of view for image-guided radiotherapy

abstract: Image guidance has been widely used in radiation therapy. Correctly identifying the bounding box of the anatomical landmarks from limited field of views is the key to success. In image-guided radiation therapy (IGRT), the detection of those landmarks like the 12th vertebra (T12) still requires tedious manual inspections and annotations; and superior-inferior misalignment to the wrong vertebral body is still relatively common. It is necessary to develop an automated approach to detect those landmarks from images. The challenges of training a model to identify T12 vertebrae automatically mainly are high shape similarity between T12 and neighboring vertebrae, limited annotated data, and class imbalance. This study proposed a novel 3D detection network, requiring only a small amount of training data. Our approach has the following innovations, including 1) the introduction of an auxiliary network to build constraint feature map for improving the model's generalization, especially when the constraint structure is easier to be detected than the main one; 2) an improved detection head and target functions for accurate bounding box detection; and 3) an improved loss functions to address the high class imbalance. Our proposed network was trained, validated and tested on anotated CT images from 55 patients and demonstrated accurate distinguish T12 vertebra from its neighboring vertebrae of high shape similarity. Our proposed algorithm yielded the bounding box center and size errors of 3.98\pm2.04mm and 16.83\pm8.34mm, respectively. Our approach significantly outperformed state-of-the-arts Retina-Net3D in average precision (AP) at IoU thresholds of 0.35 and 0.5, with AP increasing from 0 to 95.4 and 0 to 64.7, respectively. In summary, our approach has a great potential to be integrated into the clinical workflow to improve the safety of IGRT.


title: Representing Multi-Robot Structure through Multimodal Graph Embedding for the Selection of Robot Teams

abstract: Multi-robot systems of increasing size and complexity are used to solve large-scale problems, such as area exploration and search and rescue. A key decision in human-robot teaming is dividing a multi-robot system into teams to address separate issues or to accomplish a task over a large area. In order to address the problem of selecting teams in a multi-robot system, we propose a new multimodal graph embedding method to construct a unified representation that fuses multiple information modalities to describe and divide a multi-robot system. The relationship modalities are encoded as directed graphs that can encode asymmetrical relationships, which are embedded into a unified representation for each robot. Then, the constructed multimodal representation is used to determine teams based upon unsupervised learning. We perform experiments to evaluate our approach on expert-defined team formations, large-scale simulated multi-robot systems, and a system of physical robots. Experimental results show that our method successfully decides correct teams based on the multifaceted internal structures describing multi-robot systems, and outperforms baseline methods based upon only one mode of information, as well as other graph embedding-based division methods.


title: Dynamic Mode Decomposition for Construction of Reduced-Order Models of Hyperbolic Problems with Shocks

abstract: Construction of reduced-order models (ROMs) for hyperbolic conservation laws is notoriously challenging mainly due to the translational property and nonlinearity of the governing equations. While the Lagrangian framework for ROM construction resolves the translational issue, it is valid only before a shock forms. Once that occurs, characteristic lines cross each other and projection from a high-fidelity model space onto a ROM space distorts a moving grid, resulting in numerical instabilities. We address this grid distortion issue by developing a physics-aware dynamic mode decomposition (DMD) method based on hodograph transformation. The latter provides a map between the original nonlinear system and its linear counterpart, which coincides with the Koopman operator. This strategy is consistent with the spirit of physics-aware DMDs in that it retains information about shock dynamics. Several numerical examples are presented to validate the proposed physics-aware DMD approach to constructing accurate ROMs.


title: A generalization of primitive sets and a conjecture of Erd\H{o}s

abstract: A set of integers greater than 1 is primitive if no element divides another. Erd\H{o}s proved in 1935 that the sum of $1/(n \log n)$ for $n$ running over a primitive set $A$ is universally bounded over all choices for $A$. In 1988 he asked if this universal bound is attained by the set of prime numbers. We answer the Erd\H{o}s question in the affirmative for 2-primitive sets. Here a set is 2-primitive if no element divides the product of 2 other elements.


title: Liouville-type theorems for sign-changing solutions to nonlocal elliptic inequalities and systems with variable-exponent nonlinearities

abstract: We consider the fractional elliptic inequality with variable-exponent nonlinearity $$ (-\Delta)^{\frac{\alpha}{2}} u+\lambda\, \Delta u \geq |u|^{p(x)}, \quad x\in\mathbb{R}^N, $$ where $N\geq 1$, $\alpha\in (0,2)$, $\lambda\in\mathbb{R}$ is a constant, $p: \mathbb{R}^N\to (1,\infty)$ is a measurable function, and $(-\Delta)^{\frac{\alpha}{2}}$ is the fractional Laplacian operator of order $\frac{\alpha}{2}$. A Liouville-type theorem is established for the considered problem. Namely, we obtain sufficient conditions under which the only weak solution is the trivial one. Next, we extend our study to systems of fractional elliptic inequalities with variable-exponent nonlinearities. Besides the consideration of variable-exponent nonlinearities, the novelty of this work consists in investigating sign-changing solutions to the considered problems. Namely, to the best of our knowledge, only nonexistence results of positive solutions to fractional elliptic problems were invetigated previously. Our approach is based on the nonlinear capacity method combined with a pointwise estimate of the fractional Laplacian of some test functions, which was derived by Fujiwara (2018) (see also Dao and Reissig (2019)). Note that the standard nonlinear capacity method cannot be applied to the considered problems due to the change of sign of solutions.


title: Adversarial System Variant Approximation to Quantify Process Model Generalization

abstract: In process mining, process models are extracted from event logs using process discovery algorithms and are commonly assessed using multiple quality dimensions. While the metrics that measure the relationship of an extracted process model to its event log are well-studied, quantifying the level by which a process model can describe the unobserved behavior of its underlying system falls short in the literature. In this paper, a novel deep learning-based methodology called Adversarial System Variant Approximation (AVATAR) is proposed to overcome this issue. Sequence Generative Adversarial Networks are trained on the variants contained in an event log with the intention to approximate the underlying variant distribution of the system behavior. Unobserved realistic variants are sampled either directly from the Sequence Generative Adversarial Network or by leveraging the Metropolis-Hastings algorithm. The degree by which a process model relates to its underlying unknown system behavior is then quantified based on the realistic observed and estimated unobserved variants using established process model quality metrics. Significant performance improvements in revealing realistic unobserved variants are demonstrated in a controlled experiment on 15 ground truth systems. Additionally, the proposed methodology is experimentally tested and evaluated to quantify the generalization of 60 discovered process models with respect to their systems.


title: A Collective Learning Framework to Boost GNN Expressiveness

abstract: Graph Neural Networks (GNNs) have recently been used for node and graph classification tasks with great success, but GNNs model dependencies among the attributes of nearby neighboring nodes rather than dependencies among observed node labels. In this work, we consider the task of inductive node classification using GNNs in supervised and semi-supervised settings, with the goal of incorporating label dependencies. Because current GNNs are not universal (i.e., most-expressive) graph representations, we propose a general collective learning approach to increase the representation power of any existing GNN. Our framework combines ideas from collective classification with self-supervised learning, and uses a Monte Carlo approach to sampling embeddings for inductive learning across graphs. We evaluate performance on five real-world network datasets and demonstrate consistent, significant improvement in node classification accuracy, for a variety of state-of-the-art GNNs.