title: Percolation framework reveals limits of privacy in Conspiracy, Dark Web, and Blockchain networks

abstract: We consider the privacy of interactions between individuals in a network. For many networks, while nodes are anonymous to outside observers, the existence of a link between individuals implies the possibility of one node revealing identifying information about its neighbor. Moreover, while the identities of the accounts are likely hidden to an observer, the network of interaction between two anonymous accounts is often available. For example, in blockchain cryptocurrencies, transactions between two anonymous accounts are published openly. Here we consider what happens if one (or more) parties in such a network are deanonymized by an outside identity. These compromised individuals could leak information about others with whom they interacted, which could then cascade to more and more nodes' information being revealed. We use a percolation framework to analyze the scenario outlined above and show for different likelihoods of individuals possessing information on their counter-parties, the fraction of accounts that can be identified and the idealized minimum number of steps from a deanonymized node to an anonymous node (a measure of the effort required to deanonymize that individual). We further develop a greedy algorithm to estimate the \emph{actual} number of steps that will be needed to identify a particular node based on the noisy information available to the attacker. We apply our framework to three real-world networks: (1) a blockchain transaction network, (2) a network of interactions on the dark web, and (3) a political conspiracy network. We find that in all three networks, beginning from one compromised individual, it is possible to deanonymize a significant fraction of the network ($>50$%) within less than 5 steps. Overall these results provide guidelines for investigators seeking to identify actors in anonymous networks, as well as for users seeking to maintain their privacy.


title: Minmax Hierarchies, Minimal Fibrations and a PDE based Proof of the Willmore Conjecture

abstract: We introduce a general scheme that permits to generate successive min-max problems for producing critical points of higher and higher indices to Palais-Smale Functionals in normal Banach manifolds equipped with complete Finsler structures. We call the resulting tree of minmax problems a minmax hierarchy. We give several examples and in particular we explain how to implement this scheme in the framework of the viscosity method introduced by the author some years ago in order to give a new proof of the Willmore conjecture after the famous result by Marques and Neves.


title: The $\gamma$-ray sky seen at X-ray energies II: the Swift hunt of Fermi BL Lac objects among unidentified gamma-ray sources

abstract: Nearly 50% of all sources detected by the Fermi Large Area Telescope are classified as blazars or blazar candidates, one of the most elusive classes of active galaxies. Additional blazars can also be hidden within the sample of unidentified or unassociated $\gamma$-ray sources (UGSs) that constitute about one-third of all gamma-ray sources detected to date. We recently confirmed that the large majority of Fermi blazars of the BL Lac subclass have an X-ray counterpart. Using the X-ray properties of a BL Lac training set and combining these with archival multifrequency information, we aim to search for UGSs that could have a BL Lac source within their $\gamma$-ray positional uncertainty regions. We reduced and analyzed the Swift X-ray observations of a selected sample of 327 UGSs. We then compared the X-ray fluxes and hardness ratios of all sources detected in the pointed fields with those of known Fermi BL Lacs. We find at least one X-ray source, lying within the $\gamma$-ray positional uncertainty at 95% confidence level, for 223 UGSs and a total of 464 X-ray sources in all fields analyzed. The X-ray properties of a large fraction of them, eventually combined with radio, infrared, and optical information, exhibit BL Lac multi-frequency behavior, thus allowing us to select high-confidence BL Lac candidates; some of them were recently observed during our optical spectroscopic campaign which confirmed their nature. We find that out of 50 X-ray sources that were confirmed as BL Lacs through optical spectroscopy, 12 do not show canonical mid-infrared or radio BL Lac properties. This indicates that the selection of X-ray BL Lac candidates is a strong method to find new counterparts within Fermi UGSs. Finally, we pinpoint a sample of 32 Swift/XRT candidate counterparts to Fermi UGSs that are most likely BL Lac objects.


title: Efficient ancilla-free reversible and quantum circuits for the Hidden Weighted Bit function

abstract: The Hidden Weighted Bit function plays an important role in the study of classical models of computation. A common belief is that this function is exponentially hard for the implementation by reversible ancilla-free circuits, even though introducing a small number of ancillae allows a very efficient implementation. In this paper, we refute the exponential hardness conjecture by developing a polynomial-size reversible ancilla-free circuit computing the Hidden Weighted Bit function. Our circuit has size $O(n^{6.42})$, where $n$ is the number of input bits. We also show that the Hidden Weighted Bit function can be computed by a quantum ancilla-free circuit of size $O(n^2)$. The technical tools employed come from a combination of Theoretical Computer Science (Barrington's theorem) and Physics (simulation of fermionic Hamiltonians) techniques.


title: Predicting Illegal Fishing on the Patagonia Shelf from Oceanographic Seascapes

abstract: Many of the world's most important fisheries are experiencing increases in illegal fishing, undermining efforts to sustainably conserve and manage fish stocks. A major challenge to ending illegal, unreported, and unregulated (IUU) fishing is improving our ability to identify whether a vessel is fishing illegally and where illegal fishing is likely to occur in the ocean. However, monitoring the oceans is costly, time-consuming, and logistically challenging for maritime authorities to patrol. To address this problem, we use vessel tracking data and machine learning to predict illegal fishing on the Patagonian Shelf, one of the world's most productive regions for fisheries. Specifically, we focus on Chinese fishing vessels, which have consistently fished illegally in this region. We combine vessel location data with oceanographic seascapes -- classes of oceanic areas based on oceanographic variables -- as well as other remotely sensed oceanographic variables to train a series of machine learning models of varying levels of complexity. These models are able to predict whether a Chinese vessel is operating illegally with 69-96% confidence, depending on the year and predictor variables used. These results offer a promising step towards preempting illegal activities, rather than reacting to them forensically.


title: Geometric Style Transfer

abstract: Neural style transfer (NST), where an input image is rendered in the style of another image, has been a topic of considerable progress in recent years. Research over that time has been dominated by transferring aspects of color and texture, yet these factors are only one component of style. Other factors of style include composition, the projection system used, and the way in which artists warp and bend objects. Our contribution is to introduce a neural architecture that supports transfer of geometric style. Unlike recent work in this area, we are unique in being general in that we are not restricted by semantic content. This new architecture runs prior to a network that transfers texture style, enabling us to transfer texture to a warped image. This form of network supports a second novelty: we extend the NST input paradigm. Users can input content/style pair as is common, or they can chose to input a content/texture-style/geometry-style triple. This three image input paradigm divides style into two parts and so provides significantly greater versatility to the output we can produce. We provide user studies that show the quality of our output, and quantify the importance of geometric style transfer to style recognition by humans.


title: Block Basis for Coinvariants of Modular Pseudo-reflection Groups

abstract: As a sequel of \cite{Ou}, in this shot note, we investigate block basis for coinvariants of finite modular pseudo-reflection groups. We are particularly interested in the case where $ G $ is a subgroup of the parabolic subgroups of $ GL_n(q) $ which generalizes the Weyl groups of restricted Cartan type Lie algebra.


title: Invariant Brauer group of an abelian variety

abstract: We study a new object that can be attached to an abelian variety or a complex torus: the invariant Brauer group, as recently defined by Yang Cao. Over the field of complex numbers this is an elementary abelian 2-group with an explicit upper bound on the rank. We exhibit many cases in which the invariant Brauer group is zero, and construct complex abelian varieties in every dimension starting with 2, both simple and non-simple, with invariant Brauer group of order 2. We also address the situation in finite characteristic and over non-closed fields.


title: A Reinforcement Learning Approach for Fast Frequency Control in Low-Inertia Power Systems

abstract: The electric grid is undergoing a major transition from fossil fuel-based power generation to renewable energy sources, typically interfaced to the grid via power electronics. The future power systems are thus expected to face increased control complexity and challenges pertaining to frequency stability due to lower levels of inertia and damping. As a result, the frequency control and development of novel ancillary services is becoming imperative. This paper proposes a data-driven control scheme, based on Reinforcement Learning (RL), for grid-forming Voltage Source Converters (VSCs), with the goal of exploiting their fast response capabilities to provide fast frequency control to the system. A centralized RL-based controller collects generator frequencies and adjusts the VSC power output, in response to a disturbance, to prevent frequency threshold violations. The proposed control scheme is analyzed and its performance evaluated through detailed time-domain simulations of the IEEE 14-bus test system.


title: Dirac-Harper Theory for One Dimensional Moir\'e Superlattices

abstract: We study a Dirac Harper model for moir\'e bilayer superlattices where layer antisymmetric strain periodically modulates the interlayer coupling between two honeycomb lattices in one spatial dimension. Discrete and continuum formulations of this model are analyzed. For sufficiently long moir\'e period the we find low energy spectra that host a manifold of weakly dispersive bands arising from a hierarchy of momentum and position dependent mass inversions. We analyze their charge distributions, mode count and valley-coherence using exact symmetries of the lattice model and approximate symmetries of a four-flavor version of the Jackiw-Rebbi one dimensional solution.


title: Basis Decompositions and a Mathematica Package for Modular Graph Forms

abstract: Modular graph forms (MGFs) are a class of non-holomorphic modular forms which naturally appear in the low-energy expansion of closed-string genus-one amplitudes and have generated considerable interest from pure mathematicians. MGFs satisfy numerous non-trivial algebraic- and differential relations which have been studied extensively in the literature and lead to significant simplifications. In this paper, we systematically combine these relations to obtain basis decompositions of all two- and three-point MGFs of total modular weight $w+\bar{w}\leq12$, starting from just two well-known identities for banana graphs. Furthermore, we study previously known relations in the integral representation of MGFs, leading to a new understanding of holomorphic subgraph reduction as Fay identities of Kronecker--Eisenstein series and opening the door towards decomposing divergent graphs. We provide a computer implementation for the manipulation of MGFs in the form of the $\texttt{Mathematica}$ package $\texttt{ModularGraphForms}$ which includes the basis decompositions obtained.


title: Exponential Convergence of Gradient Methods in Concave Network Zero-sum Games

abstract: Motivated by Generative Adversarial Networks, we study the computation of Nash equilibrium in concave network zero-sum games (NZSGs), a multiplayer generalization of two-player zero-sum games first proposed with linear payoffs. Extending previous results, we show that various game theoretic properties of convex-concave two-player zero-sum games are preserved in this generalization. We then generalize last iterate convergence results obtained previously in two-player zero-sum games. We analyze convergence rates when players update their strategies using Gradient Ascent, and its variant, Optimistic Gradient Ascent, showing last iterate convergence in three settings -- when the payoffs of players are linear, strongly concave and Lipschitz, and strongly concave and smooth. We provide experimental results that support these theoretical findings.


title: ATHENA -- a pre-university study programme at the university of Geneva

abstract: Athena (named after the ancient Greek goddess of wisdom) is a pre-university study programme for mathematics and physics organised by the Faculty of Science at the University of Geneva. It targets pupils enrolled in the final or penultimate year of Secondary II (high school), giving them an opportunity to explore and discover university-level studies in mathematics and physics. The programme aims to enhance pupils interest for the physical and mathematical sciences by introducing them to new topics, all while giving them a taste for student life. It also seeks to promote scientific careers to young pupils, especially to young women, as well as improving the transition between Secondary II and university


title: Impact of Legal Requirements on Explainability in Machine Learning

abstract: The requirements on explainability imposed by European laws and their implications for machine learning (ML) models are not always clear. In that perspective, our research analyzes explanation obligations imposed for private and public decision-making, and how they can be implemented by machine learning techniques.


title: Additive and geometric transversality of fractal sets in the integers

abstract: By juxtaposing ideas from fractal geometry and dynamical systems, Furstenberg proposed a series of conjectures in the late 1960's that explore the relationship between digit expansions with respect to multiplicatively independent bases. In this work, we introduce and study - in the discrete context of the integers - analogues of some of the notions and results surrounding Furstenberg's work. In particular, we define a new class of fractal sets of integers that parallels the notion of $\times r$-invariant sets on the 1-torus and investigate the additive and geometric independence between two such fractal sets when they are structured with respect to multiplicatively independent bases. Our main results in this direction parallel the works of Furstenberg, Hochman-Shmerkin, Shmerkin, Wu, and Lindenstrauss-Meiri-Peres and include: -a classification of all subsets of the positive integers that are simultaneously $\times r$- and $\times s$-invariant; -integer analogues of two of Furstenberg's transversality conjectures pertaining to the dimensions of the intersection $A\cap B$ and the sumset $A+B$ of $\times r$- and $\times s$-invariant sets $A$ and $B$ when $r$ and $s$ are multiplicatively independent; and -a description of the dimension of iterated sumsets $A+A+\cdots+A$ for any $\times r$-invariant set $A$. We achieve these results by combining ideas from fractal geometry and ergodic theory to build a bridge between the continuous and discrete regimes. For the transversality results, we rely heavily on quantitative bounds on the $L^q$-dimensions of projections of restricted digit Cantor measures obtained recently by Shmerkin. We end by outlining a number of open questions and directions regarding fractal subsets of the integers.


title: STaRFlow: A SpatioTemporal Recurrent Cell for Lightweight Multi-Frame Optical Flow Estimation

abstract: We present a new lightweight CNN-based algorithm for multi-frame optical flow estimation. Our solution introduces a double recurrence over spatial scale and time through repeated use of a generic "STaR" (SpatioTemporal Recurrent) cell. It includes (i) a temporal recurrence based on conveying learned features rather than optical flow estimates; (ii) an occlusion detection process which is coupled with optical flow estimation and therefore uses a very limited number of extra parameters. The resulting STaRFlow algorithm gives state-of-the-art performances on MPI Sintel and Kitti2015 and involves significantly less parameters than all other methods with comparable results.


title: Optimal additive quaternary codes of low dimension

abstract: An additive quaternary $[n,k,d]$-code (length $n,$ quaternary dimension $k,$ minimum distance $d$) is a $2k$-dimensional F_2-vector space of $n$-tuples with entries in $Z_2\times Z_2$ (the $2$-dimensional vector space over F_2) with minimum Hamming distance $d.$ We determine the optimal parameters of additive quaternary codes of dimension $k\leq 3.$ The most challenging case is dimension $k=2.5.$ We prove that an additive quaternary $[n,2.5,d]$-code where $d<n-1$ exists if and only if $3(n-d)\geq \lceil d/2\rceil +\lceil d/4\rceil +\lceil d/8\rceil$. In particular we construct new optimal $2.5$-dimensional additive quaternary codes. As a by-product we give a direct proof for the fact that a binary linear $[3m,5,2e]_2$-code for $e<m-1$ exists if and only if the Griesmer bound $3(m-e)\geq \lceil e/2\rceil +\lceil e/4\rceil+\lceil e/8\rceil$ is satisfied.


title: Generic Caldero-Chapoton functions with coefficients and applications to surface cluster algebras

abstract: We realize Derksen-Weyman-Zelevinsky's mutations of representations as densely-defined regular maps on representation spaces, and study the generic values of Caldero-Chapoton functions with coefficients, giving, for instance, a sufficient combinatorial condition for their linear independence. For a quiver with potential $(Q,S)$, we show that if $k$ is a vertex not incident to any oriented 2-cycle, then every generically $\tau$-reduced irreducible component $Z$ of any affine variety of (decorated) representations has a dense open subset $U$ on which Derksen-Weyman-Zelevinsky's mutation of representations $\mu_k$ can be defined consistently as a regular map to an affine variety of (decorated) representations of the Jacobian algebra of the mutated QP $\mu_k(Q,S)$. Our techniques involve only basic linear algebra and elementary algebraic geometry, and do not require to assume Jacobi-finiteness. Thus, the paper yields a new and more general proof of the mutation invariance of generic Caldero-Chapoton functions, generalizing and providing a new natural geometric perspective on results of Derksen-Weyman-Zelevinsky and Plamondon. (For Jacobi-finite non-degenerate quivers with potential, this invariance was shown by Plamondon using the machinery of Ginzburg dg-algebras and Hom-finite generalized cluster categories.) We apply our results, together with results of Mills, Muller and Qin, to prove that for any choice of geometric coefficient systems, not necessarily of full rank, the cluster algebra associated to a possibly punctured surface with at least two marked points on the boundary has the generic Caldero-Chapoton functions as a basis over the Laurent polynomial ring of coefficients.


title: SPIDERS: overview of the X-ray galaxy cluster follow-up and the final spectroscopic data release

abstract: SPIDERS (The SPectroscopic IDentification of eROSITA Sources) is a large spectroscopic programme for X-ray selected galaxy clusters as part of the Sloan Digital Sky Survey-IV (SDSS-IV). We describe the final dataset in the context of SDSS Data Release 16 (DR16): the survey overall characteristics, final targeting strategies, achieved completeness and spectral quality, with special emphasis on its use as a galaxy cluster sample for cosmology applications. SPIDERS now consists of about 27,000 new optical spectra of galaxies selected within 4,000 photometric red sequences, each associated with an X-ray source. The excellent spectrograph efficiency and a robust analysis pipeline yield a spectroscopic redshift measurement success rate exceeding 98%, with a median velocity accuracy of 20 km s$^{-1}$ (at $z=0.2$). Using the catalogue of 2,740 X-ray galaxy clusters confirmed with DR16 spectroscopy, we reveal the three-dimensional map of the galaxy cluster distribution in the observable Universe up to $z\sim0.6$. We highlight the homogeneity of the member galaxy spectra among distinct regions of the galaxy cluster phase space. Aided by accurate spectroscopic redshifts and by a model of the sample selection effects, we compute the galaxy cluster X-ray luminosity function and we present its lack of evolution up to $z=0.6$. Finally we discuss the prospects of forthcoming large multiplexed spectroscopic programmes dedicated to follow up the next generation of all-sky X-ray source catalogues.


title: Partially phase-locked solutions to the Kuramoto model

abstract: The Kuramoto model is a canonical model for understanding phase-locking phenomenon. It is well-understood that, in the usual mean-field scaling, full phase-locking is unlikely and that it is partially phase-locked states that are important in applications. Despite this, while there has been much attention given to the existence and stability of fully phase-locked states in the finite N Kuramoto model, the partially phase-locked states have received much less attention. In this paper, we present two related results. Firstly, we derive an analytical criterion that, for sufficiently strong coupling, guarantees the existence of a partially phase-locked state by proving the existence of an attracting ball around a fixed point of a subset of the oscillators. We also derive a larger invariant ball such that any point in it will asymptotically converge to the attracting ball. Secondly, we consider the large N (thermodynamic) limit for the Kuramoto system with randomly distributed frequencies. Using some results of De Smet and Aeyels on partial entrainment, we derive a deterministic condition giving almost sure existence of a partially entrained state for sufficiently strong coupling when the natural frequencies of the individual oscillators are independent identically distributed random variables, as well as upper and lower bounds on the size of the largest cluster of partially entrained oscillators. Interestingly in a series on numerical experiments we find that the observed size of the largest entrained cluster is predicted extremely well by the upper bound.


title: Traversable wormholes and the Brouwer fixed-point theorem

abstract: The Brouwer fixed-point theorem in topology states that for any continuous mapping $f$ on a compact convex set into itself admits a fixed point, i.e., a point $x_0$ such that $f(x_0)=x_0$. Under certain conditions, this fixed point corresponds to the throat of a traversable wormhole, i.e., $b(r_0)=r_0$ for the shape function $b=b(r)$. The possible existence of wormholes can therefore be deduced from purely mathematical considerations without going beyond the existing physical requirements.


title: Observation of second sound in a rapidly varying temperature field in Ge

abstract: Second sound is known as the thermal transport regime where heat is carried by temperature waves. Its experimental observation was previously restricted to a small number of materials, usually in rather narrow temperature windows. We show that it is possible to overcome these limitations by driving the system with a rapidly varying temperature field. This effect is demonstrated in bulk Ge between 7 kelvin and room temperature, studying the phase lag of the thermal response under a harmonic high frequency external thermal excitation, addressing the relaxation time and the propagation velocity of the heat waves. These results provide a new route to investigate the potential of wave-like heat transport in almost any material, opening opportunities to control heat through its oscillatory nature.


title: The impact of strong recombination on temperature determination in planetary nebulae

abstract: The long-standing difference in chemical abundances determined from optical recombination lines and collisionally excited lines raises questions about our understanding of atomic physics, as well as the assumptions made when determining physical conditions and chemical abundances in astrophysical nebulae. Here, we study the recombination contribution of [O III] 4363 and the validity of the line ratio [O III] 4363/4959 as a temperature diagnostic in planetary nebulae with a high abundance discrepancy. We derive a fit for the recombination coefficient of [O III] 4363 that takes into account the radiative and dielectronic recombinations, for electron temperatures from 200 to 30,000 K. We estimate the recombination contribution of [O III] 4363 for the planetary nebulae Abell 46 and NGC 6778 by subtracting the collisional contribution from the total observed flux. We find that the spatial distribution for the estimated recombination contribution in [O III] 4363 follows that of the O II 4649 recombination line, both peaking in the central regions of the nebula, especially in the case of Abell 46 which has a much higher abundance discrepancy. The estimated recombination contribution reaches up to 70% and 40% of the total [O III] 4363 observed flux, for Abell 46 and NGC 6778, respectively.


title: Class of images of Abel maps on normal surface singularities

abstract: In this paper we investigate Abel maps on normal surface singularities described in \cite{NNI}. We investigate the affine version of the class of the images of Abel maps on normal surface singularities. More precisely we consider the projective clousure of the image of an Abel map, its dual projective variety and we substract from its degree the multiplicity of the infinite hyperplane on the dual variety. In the case of generic singularities we prove explicit combinatorial formulas of this invariant, in the general case we prove an upper bound.


title: Camera-Lidar Integration: Probabilistic sensor fusion for semantic mapping

abstract: An automated vehicle operating in an urban environment must be able to perceive and recognise object/obstacles in a three-dimensional world while navigating in a constantly changing environment. In order to plan and execute accurate sophisticated driving maneuvers, a high-level contextual understanding of the surroundings is essential. Due to the recent progress in image processing, it is now possible to obtain high definition semantic information in 2D from monocular cameras, though cameras cannot reliably provide the highly accurate 3D information provided by lasers. The fusion of these two sensor modalities can overcome the shortcomings of each individual sensor, though there are a number of important challenges that need to be addressed in a probabilistic manner. In this paper, we address the common, yet challenging, lidar/camera/semantic fusion problems which are seldom approached in a wholly probabilistic manner. Our approach is capable of using a multi-sensor platform to build a three-dimensional semantic voxelized map that considers the uncertainty of all of the processes involved. We present a probabilistic pipeline that incorporates uncertainties from the sensor readings (cameras, lidar, IMU and wheel encoders), compensation for the motion of the vehicle, and heuristic label probabilities for the semantic images. We also present a novel and efficient viewpoint validation algorithm to check for occlusions from the camera frames. A probabilistic projection is performed from the camera images to the lidar point cloud. Each labelled lidar scan then feeds into an octree map building algorithm that updates the class probabilities of the map voxels every time a new observation is available. We validate our approach using a set of qualitative and quantitative experimental tests on the USyd Dataset.


title: Habitability Models for Planetary Sciences

abstract: Habitability has been generally defined as the capability of an environment to support life. Ecologists have been using Habitat Suitability Models (HSMs) for more than four decades to study the habitability of Earth from local to global scales. Astrobiologists have been proposing different habitability models for some time, with little integration and consistency between them and different in function to those used by ecologists. In this white paper, we suggest a mass-energy habitability model as an example of how to adapt and expand the models used by ecologists to the astrobiology field. We propose to implement these models into a NASA Habitability Standard (NHS) to standardize the habitability objectives of planetary missions. These standards will help to compare and characterize potentially habitable environments, prioritize target selections, and study correlations between habitability and biosignatures. Habitability models are the foundation of planetary habitability science. The synergy between the methods used by ecologists and astrobiologists will help to integrate and expand our understanding of the habitability of Earth, the Solar System, and exoplanets.


title: XSleepNet: Multi-View Sequential Model for Automatic Sleep Staging

abstract: Automating sleep staging is vital to scale up sleep assessment and diagnosis to serve millions experiencing sleep deprivation and disorders and enable longitudinal sleep monitoring in home environments. Learning from raw polysomnography signals and their derived time-frequency image representations has been prevalent. However, learning from multi-view inputs (e.g., both the raw signals and the time-frequency images) for sleep staging is difficult and not well understood. This work proposes a sequence-to-sequence sleep staging model, XSleepNet, that is capable of learning a joint representation from both raw signals and time-frequency images. Since different views may generalize or overfit at different rates, the proposed network is trained such that the learning pace on each view is adapted based on their generalization/overfitting behavior. In simple terms, the learning on a particular view is speeded up when it is generalizing well and slowed down when it is overfitting. View-specific generalization/overfitting measures are computed on-the-fly during the training course and used to derive weights to blend the gradients from different views. As a result, the network is able to retain the representation power of different views in the joint features which represent the underlying distribution better than those learned by each individual view alone. Furthermore, the XSleepNet architecture is principally designed to gain robustness to the amount of training data and to increase the complementarity between the input views. Experimental results on five databases of different sizes show that XSleepNet consistently outperforms the single-view baselines and the multi-view baseline with a simple fusion strategy. Finally, XSleepNet also outperforms prior sleep staging methods and improves previous state-of-the-art results on the experimental databases.


title: Learning from Data to Optimize Control in Precision Farming

abstract: Precision farming is one way of many to meet a 70 percent increase in global demand for agricultural products on current agricultural land by 2050 at reduced need of fertilizers and efficient use of water resources. The catalyst for the emergence of precision farming has been satellite positioning and navigation followed by Internet-of-Things, generating vast information that can be used to optimize farming processes in real-time. Statistical tools from data mining, predictive modeling, and machine learning analyze pattern in historical data, to make predictions about future events as well as intelligent actions. This special issue presents the latest development in statistical inference, machine learning and optimum control for precision farming.


title: Automatic Detection of COVID-19 Cases on X-ray images Using Convolutional Neural Networks

abstract: In recent months the world has been surprised by the rapid advance of COVID-19. In order to face this disease and minimize its socio-economic impacts, in addition to surveillance and treatment, diagnosis is a crucial procedure. However, the realization of this is hampered by the delay and the limited access to laboratory tests, demanding new strategies to carry out case triage. In this scenario, deep learning models are being proposed as a possible option to assist the diagnostic process based on chest X-ray and computed tomography images. Therefore, this research aims to automate the process of detecting COVID-19 cases from chest images, using convolutional neural networks (CNN) through deep learning techniques. The results can contribute to expand access to other forms of detection of COVID-19 and to speed up the process of identifying this disease. All databases used, the codes built, and the results obtained from the models' training are available for open access. This action facilitates the involvement of other researchers in enhancing these models since this can contribute to the improvement of results and, consequently, the progress in confronting COVID-19.


title: Disease Detectives: Using Mathematics to Forecast the Spread of Infectious Diseases

abstract: The COVID-19 pandemic has led to significant changes in how people are currently living their lives. To determine how to best reduce the effects of the pandemic and start reopening societies, governments have drawn insights from mathematical models of the spread of infectious diseases. In this article, we give an introduction to a family of mathematical models (called "compartmental models") and discuss how the results of analyzing these models influence government policies and human behavior, such as encouraging mask wearing and physical distancing to help slow the spread of the disease.

categories: q-bio.OT math.DS math.HO nlin.AO physics.soc-ph


title: In Response to COVID-19: Configuration Model of the Epidemic Spreading

abstract: A configuration model for epidemic spread, based on a scale-free network, is introduced. We obtain the analytical solutions describing both unstable and stable dynamics of the epidemic spreading, and demonstrate how these regimes can interchange during the epidemic. We apply the model to the COVID-19 case and demonstrate the predictive features of our model.


title: FitzHugh-Nagumo oscillators on complex networks mimic epileptic-seizure-related synchronization phenomena

abstract: We study patterns of partial synchronization in a network of FitzHugh-Nagumo oscillators with empirical structural connectivity measured in human subjects. We report the spontaneous occurrence of synchronization phenomena that closely resemble the ones seen during epileptic seizures in humans. In order to obtain deeper insights into the interplay between dynamics and network topology, we perform long-term simulations of oscillatory dynamics on different paradigmatic network structures: random networks, regular nonlocally coupled ring networks, ring networks with fractal connectivities, and small-world networks with various rewiring probability. Among these networks, a small-world network with intermediate rewiring probability best mimics the findings achieved with the simulations using the empirical structural connectivity. For the other network topologies, either no spontaneously occurring epileptic-seizure-related synchronization phenomena can be observed in the simulated dynamics, or the overall degree of synchronization remains high throughout the simulation. This indicates that a topology with some balance between regularity and randomness favors the self-initiation and self-termination of episodes of seizure-like strong synchronization.


title: Formality of differential graded algebras and complex Lagrangian submanifolds

abstract: Let $i: \mathrm{L} \hookrightarrow \mathrm{X}$ be a compact K\"{a}hler Lagrangian in a holomorphic symplectic variety $\mathrm{X}/\mathbf{C}$. We use deformation quantisation to show that the endomorphism differential graded algebra $\mathrm{RHom}\big(i_*\mathrm{K}_\mathrm{L}^{1/2},i_*\mathrm{K}_\mathrm{L}^{1/2}\big)$ is formal. We prove a generalisation to pairs of Lagrangians, along with auxiliary results on the behaviour of formality in families of $\mathrm{A}_\infty$-modules.


title: Not Your Grandfathers Test Set: Reducing Labeling Effort for Testing

abstract: Building and maintaining high-quality test sets remains a laborious and expensive task. As a result, test sets in the real world are often not properly kept up to date and drift from the production traffic they are supposed to represent. The frequency and severity of this drift raises serious concerns over the value of manually labeled test sets in the QA process. This paper proposes a simple but effective technique that drastically reduces the effort needed to construct and maintain a high-quality test set (reducing labeling effort by 80-100% across a range of practical scenarios). This result encourages a fundamental rethinking of the testing process by both practitioners, who can use these techniques immediately to improve their testing, and researchers who can help address many of the open questions raised by this new approach.


title: Scientific Discovery by Generating Counterfactuals using Image Translation

abstract: Model explanation techniques play a critical role in understanding the source of a model's performance and making its decisions transparent. Here we investigate if explanation techniques can also be used as a mechanism for scientific discovery. We make three contributions: first, we propose a framework to convert predictions from explanation techniques to a mechanism of discovery. Second, we show how generative models in combination with black-box predictors can be used to generate hypotheses (without human priors) that can be critically examined. Third, with these techniques we study classification models for retinal images predicting Diabetic Macular Edema (DME), where recent work showed that a CNN trained on these images is likely learning novel features in the image. We demonstrate that the proposed framework is able to explain the underlying scientific mechanism, thus bridging the gap between the model's performance and human understanding.


title: The determination of the spin and parity of a vector-vector system

abstract: We present a construction of the reaction amplitude for the inclusive production of a resonance decaying to a pair of identical vector particles such as $J/\psi J/\psi$, $\rho\rho$, $\phi\phi$. The method provides the possibility of determining the spin and parity of a resonance in a model-independent way. The methodology is demonstrated using the Standard Model decay of the Higgs boson to four leptons and through angular correlations in the production and decays of $J/\psi$ pairs.


title: Joint Information Theoretic Secrecy and Covert Communication in the Presence of an Untrusted User and Warden

abstract: In this paper, we investigate joint information theoretic secrecy and covert communication in a single-input multi-output (SIMO) system where a transmitter (Alice) is communicating with two legitimate users (Bob and Carol). We consider that an untrusted user and a warden node are also present in the network attempting to attack the secure and covert communications to Bob and Carol, respectively. Specifically, Bob requires secure communications such that his messages from Alice are not decoded by the untrusted user, while Carol requires covert communications such that her messages from Alice are not detected by the warden. To do so, we consider that Alice transmits Carol's messages during selected time slots to hide them from the warden while also transmitting Bob's messages in each time slot contentiously. We formulate an optimization problem with the aim of maximizing the average rate subject to a covert communication requirement and a secure communications constraint. Since the proposed optimization problem is non-convex, we utilize successive convex approximation to obtain a tractable solution. We consider practical assumptions that Alice has imperfect knowledge of the warden's location and imperfect channel state information (CSI) of Bob and Carol. Our numerical examples highlight that the imperfect CSI at Carol has a more detrimental impact on the average rate compared to imperfect CSI at Bob.


title: Predicting Bit Error Rate from Meta Information using Random Forests

abstract: With the increasing power of machine learning-based reasoning, the use of meta-information (e.g., digital signal modulation parameters, channel conditions, etc.) to predict the performance of various signal processing techniques has become feasible. One such problem of practical interest is choosing a proper interference mitigation method based on the meta information of the received signal. Since heuristic table-based methods suffer from limited prediction capability for unseen cases, we propose a recommendation system based on the use of Random Forests (RF). Specifically, RF used to predict the Bit-Error-Rate (BER) of all mitigation approaches so as to determine the approach with the best performance. We found RF can predict BER with high accuracy, and its importance factor demonstrates which input attributes matter most. These BER prediction results can also benefit other functions such as adaptive modulation, channel sensing, beaming selection, etc.


title: Covariance of the matter power spectrum including the survey window function effect: N-body simulations vs. fifth-order perturbation theory on grid

abstract: We present a Next-to-next-to-leading (fifth or NNLO) order calculation for the covariance matrix of the matter power spectrum, taking into account the effect of survey window functions. Using the grid-based calculation scheme for the standard perturbation theory, GridSPT, we quickly generate multiple realizations of the nonlinear density fields to fifth order in perturbation theory, then estimate the power spectrum and the covariance matrix from the sample. To the end, we have obtained the non-Gaussian covariance originated from the one-loop trispectrum without explicitly computing the trispectrum. By comparing the GridSPT calculations with the N-body results, we show that NNLO GridSPT result reproduces the N-body results on quasi-linear scales where SPT accurately models nonlinear matter power spectrum. Incorporating the survey window function effect to GridSPT is rather straightforward, and the resulting NNLO covariance matrix also matches well with the N-body results.


title: Neural Knowledge Extraction From Cloud Service Incidents

abstract: In the last decade, two paradigm shifts have reshaped the software industry - the move from boxed products to services and the widespread adoption of cloud computing. This has had a huge impact on the software development life cycle and the DevOps processes. Particularly, incident management has become critical for developing and operating large-scale services. Incidents are created to ensure timely communication of service issues and, also, their resolution. Prior work on incident management has been heavily focused on the challenges with incident triaging and de-duplication. In this work, we address the fundamental problem of structured knowledge extraction from service incidents. We have built SoftNER, a framework for unsupervised knowledge extraction from service incidents. We frame the knowledge extraction problem as a Named-entity Recognition task for extracting factual information. SoftNER leverages structural patterns like key,value pairs and tables for bootstrapping the training data. Further, we build a novel multi-task learning based BiLSTM-CRF model which leverages not just the semantic context but also the data-types for named-entity extraction. We have deployed SoftNER at Microsoft, a major cloud service provider and have evaluated it on more than 2 months of cloud incidents. We show that the unsupervised machine learning based approach has a high precision of 0.96. Our multi-task learning based deep learning model also outperforms the state of the art NER models. Lastly, using the knowledge extracted by SoftNER we are able to build significantly more accurate models for important downstream tasks like incident triaging.


title: Localization of peripheral reactions and sensitivity to the imaginary potential

abstract: The aim of the present study is to make for the first time in the literature a systematic and quantitative assessment of the evaluation of the imaginary part of the optical potential calculated within the folding model and its consequences on the localization of surface reactions. Comparing theoretical and experimental reaction cross sections, for some light projectiles on a $^9$Be target, it has recently been shown that a single-folded s.f. (light-) nucleus-$^9$Be imaginary optical potential is more accurate than a double-folded d.f. optical potential. Within the eikonal formalism for the cross sections and phase shifts, the single-folded potential was obtained using a n-$^9$Be phenomenological optical potential and microscopic projectile densities. This paper is a follow-up in which we systematically study a series of different light and medium-mass projectile induced reactions on $^9$Be. Our results confirm that the s.f. cross sections are larger than the d.f. cross sections and the effect increases with the projectile mass. Furthermore the strong absorption radius parameter extracted from the $S$ matrices calculated with the s.f. has a stable value $r_s$ =1.3 - 1.4 fm for all projectile masses in the range of incident energies 40-100AMeV. This indicates that a clear geometrical separation can be made between the region of surface reactions, the region of strong absorption into other channels and the region of weak nuclear interaction. The d.f. results are instead much scattered and the separation between surface reactions and other channels does not seem to be consistent. Excellent agreement with recent experimental results confirms the validity of our approach.


title: Experimental Modeling of Cyclists Fatigue and Recovery Dynamics Enabling Optimal Pacing in a Time Trial

abstract: Improving a cyclist performance during a time-trial effort has been a challenge for sport scientists for several decades. There has been a lot of work on understanding the physiological concepts behind it. The concepts of Critical Power (CP) and Anaerobic Work Capacity (AWC) have been discussed often in recent cycling performance related articles. CP is a power that can be maintained by a cyclist for a long time; meaning pedaling at or below this limit, theoretically, can be continued for infinite amount of time. However, there is a limited source of energy for generating power above CP. This limited energy source is AWC. After burning energy from this tank, a cyclist can recover some by pedaling below CP. In this paper we utilize the concepts of CP and AWC to mathematically model muscle fatigue and recovery of a cyclist. Then, the models are used to formulate an optimal control problem for a time trial effort on a 10.3 km course located in Greenville SC. The course is simulated in a laboratory environment using a CompuTrainer. At the end, the optimal simulation results are compared to the performance of one subject on CompuTrainer.


title: Detection of Fe\,{\sc i} Emission in the Day-side Spectrum of WASP-33b

abstract: We analyze the high-resolution emission spectrum of WASP-33b taken using the High Dispersion Spectrograph (R\,$\approx$\,165,000) on the 8.2-m Subaru telescope. The data cover $\lambda$\,$\approx$\,$6170$-$8817$\,\AA, divided over 30 spectral orders. The telluric and stellar lines are removed using a de-trending algorithm, {\sc SysRem}, before cross-correlating with planetary spectral templates. We calculate the templates assuming a 1-D plane-parallel hydrostatic atmosphere including continuum opacity of bound-free H$^{-}$ and Rayleigh scattering by H$_{2}$ with a range of constant abundances of Fe\,{\sc i}. Using a likelihood-mapping analysis, we detect an Fe\,{\sc i} emission signature at 6.4-$\sigma$ located at $K_{\mathrm{p}}$ of 226.0\,$^{+2.1}_{-2.3}$\,km\,s$^{-1}$and $v_{\mathrm{sys}}$ of -3.2\,$^{+2.1}_{-1.8}$\,km\,s$^{-1}$ -- consistent with the planet's expected velocity in the literature. We also confirm the existence of a thermal inversion in the day-side of the planet which is very likely to be caused by the presence of Fe\,{\sc i} and previously-detected TiO in the atmosphere. This makes WASP-33b one of the prime targets to study the relative contributions of both species to the energy budget of an ultra-hot Jupiter.


title: Strong gravitational lensing by a strongly naked null singularity

abstract: We study strong gravitational lensing in a static, spherically symmetric, naked singularity spacetime, without a photon sphere. The nature of the singularity is found to be lightlike. We discuss the characteristic lensing features of this naked singularity in the strong deflection limit. In spite of the absence of a photon sphere in this spacetime, the bending angle of light diverges, as it approaches the singularity. However, unlike black holes, it is found that the nature of this divergence is nonlogarithmic, and we derive an analytic formula for the same. Moreover, the relativistic rings produced due to strong lensing by the singularity are found to be well separated from each other, making them easy to resolve and possibly detect. These features are expected to be important in the study of strong lensing by ultracompact objects, especially ones without event horizons.


title: McKay Matrices for Finite-dimensional Hopf Algebras

abstract: For a finite-dimensional Hopf algebra $A$, the McKay matrix $M_V$ of an $A$-module $V$ encodes the relations for tensoring the simple $A$-modules with $V$. We prove results about the eigenvalues and the right and left (generalized) eigenvectors of $M_V$ by relating them to characters. We show how the projective McKay matrix $Q_V$ obtained by tensoring the projective indecomposable modules of $A$ with $V$ is related to the McKay matrix of the dual module of $V$. We illustrate these results for the Drinfeld double $D_n$ of the Taft algebra by deriving expressions for the eigenvalues and eigenvectors of $M_V$ and $Q_V$ in terms of several kinds of Chebyshev polynomials. For the matrix $N_V$ that encodes the fusion rules for tensoring $V$ with a basis of projective indecomposable $D_n$-modules for the image of the Cartan map, we show that the eigenvalues and eigenvectors also have such Chebyshev expressions.


title: Nonlinear dynamics of emergent traveling waves in a reaction-Cattaneo system

abstract: Standard diffusion equation is based on Brownian motion of the dispersing species without considering persistence in the movement of the individuals. This description allows for the instantaneous spreading of the transported species over an arbitrarily large distances from their original location predicting infinite velocities. This feature is unrealistic particularly while considering biological invasion dynamics and a better description needs the consideration of dispersal with inertia. We here examine the behavior of non-infinitesimal perturbation on the steady state of an one-dimensional reaction-Cattaneo system with a cubic polynomial source term describing population dynamics or flame propagation models. It has been shown analytically that while linear analysis predicts stability of the homogeneous state, consideration of nonlinear contribution leads to a growth of spatiotemporal perturbation as a traveling wave. We show that the presence of a small finite relaxation time of the diffusive flux modifies the speed of the traveling wave. Specifically, we find that the wave speed decays with an increase of a finite relaxation time of flux. Our analytical predictions are well corroborated with the numerical results.


title: Comment on "Casimir effect in a weak gravitational field: Schwinger's approach"

abstract: We show that the statement in F. Sorge [Class. Quant. Grav. 36, no.23, 235006 (2019)] that the Casimir effect receives second order corrections due to gravity is not consistent. We remark especially on the tracing of the proper time Hamiltonian, where the correct procedure is to use the eigenfunctions and eigenvalues of the covariant DAlembertian. After some cancellations we find that the value of the functional W [0] is the same as obtained by Sorge. However, we argue that the proper vacuum energy density carries extra space-time volume terms that cancel over the gravitational correction, returning to the same expression as in Minkowski space-time.


title: An Active-to-Sterile Neutrino Transition Dipole Moment and the XENON1T Excess

abstract: In this short letter, we find that a magnetic transition dipole moment between tau and sterile neutrinos can account for the XENON1T excess events. Unlike the ordinary neutrino dipole moment, the introduction of the new sterile mass scale allows for astrophysical bounds to be suppressed. Interestingly, the best-fit regions that are compatible with the SN1987A imply either boron-8 or CNO neutrinos as the source flux. We find that sterile neutrinos of either $\sim$ 260 keV or in the $\sim$(500 - 800) keV mass range are capable of evading astrophysical constraints while being able to successfully explain the XENON1T event rate. The sterile neutrino in the best fit parameter space may have significant effects on big bang nucleosynthesis (BBN). We show the region in which a low reheating temperature of the Universe may allow the BBN constraints to be alleviated.


title: The EBLM project. VII. Spin-orbit alignment for the circumbinary planet host EBLM J0608-59 A/TOI-1338 A

abstract: A dozen short-period detached binaries are known to host transiting circumbinary planets. In all circumbinary systems so far, the planetary and binary orbits are aligned within a couple of degrees. However, the obliquity of the primary star, which is an important tracer of their formation, evolution, and tidal history, has only been measured in one circumbinary system until now. EBLM J0608-59/TOI-1338 is a low-mass eclipsing binary system with a recently discovered circumbinary planet identified by TESS. Here, we perform high-resolution spectroscopy during primary eclipse to measure the projected stellar obliquity of the primary component. The obliquity is low, and thus the primary star is aligned with the binary and planetary orbits with a projected spin-orbit angle $\beta = 2.8 \pm 17.1$ deg. The rotation period of $18.1 \pm 1.6$ days implied by our measurement of $v\sin{i_\star}$ suggests that the primary has not yet pseudo-synchronized with the binary orbit, but is consistent with gyrochronology and weak tidal interaction with the binary companion. Our result, combined with the known coplanarity of the binary and planet orbits, is suggestive of formation from a single disc. Finally, we considered whether the spectrum of the faint secondary star could affect our measurements. We show through simulations that the effect is negligible for our system, but can lead to strong biases in $v\sin{i_\star}$ and $\beta$ for higher flux ratios. We encourage future studies in eclipse spectroscopy test the assumption of a dark secondary for flux ratios $\gtrsim 1$ ppt.


title: AViD Dataset: Anonymized Videos from Diverse Countries

abstract: We introduce a new public video dataset for action recognition: Anonymized Videos from Diverse countries (AViD). Unlike existing public video datasets, AViD is a collection of action videos from many different countries. The motivation is to create a public dataset that would benefit training and pretraining of action recognition models for everybody, rather than making it useful for limited countries. Further, all the face identities in the AViD videos are properly anonymized to protect their privacy. It also is a static dataset where each video is licensed with the creative commons license. We confirm that most of the existing video datasets are statistically biased to only capture action videos from a limited number of countries. We experimentally illustrate that models trained with such biased datasets do not transfer perfectly to action videos from the other countries, and show that AViD addresses such problem. We also confirm that the new AViD dataset could serve as a good dataset for pretraining the models, performing comparably or better than prior datasets.


title: A Causal Linear Model to Quantify Edge Flow and Edge Unfairness for UnfairEdge Prioritization and Discrimination Removal

abstract: Law enforcement must prioritize sources of unfairness before mitigating their underlying unfairness, considering that they have limited resources. Unlike previous works that only make cautionary claims of discrimination and de-biases data after its generation, this paper attempts to prioritize unfair sources before mitigating their unfairness in the real-world. We assume that a causal bayesian network, representative of the data generation procedure, along with the sensitive nodes, that result in unfairness, are given. We quantify Edge Flow, which is the belief flowing along an edge by attenuating the indirect path influences, and use it to quantify Edge Unfairness. We prove that cumulative unfairness is non-existent in any decision, like judicial bail, towards any sensitive groups, like race, when the edge unfairness is absent, given an error-free linear model of conditional probability. We then measure the potential to mitigate the cumulative unfairness when edge unfairness is decreased. Based on these measures, we propose an unfair edge prioritization algorithm that prioritizes the unfair edges and a discrimination removal procedure that de-biases the generated data distribution. The experimental section validates the specifications used for quantifying the above measures.


title: Global fits of axion-like particles to XENON1T and astrophysical data

abstract: The excess of electron recoil events seen by the XENON1T experiment has been interpreted as a potential signal of axion-like particles (ALPs), either produced in the Sun, or constituting part of the dark matter halo of the Milky Way. It has also been explained as a consequence of trace amounts of tritium in the experiment. We consider the evidence for the solar and dark-matter ALP hypotheses from the combination of XENON1T data and multiple astrophysical probes, including horizontal branch stars, red giants, and white dwarfs. We briefly address the influence of ALP decays and supernova cooling. While the different datasets are in clear tension for the case of solar ALPs, all measurements can be simultaneously accommodated for the case of a sub-dominant fraction of dark-matter ALPs. Nevertheless, this solution requires the tuning of several a priori unknown parameters, such that for our choices of priors a Bayesian analysis shows no strong preference for the ALP interpretation of the XENON1T excess over the background hypothesis.


title: Rheology of mixed motor ensembles

abstract: The rheology of biological cells is not only determined by their cytoskeletal networks, but also by the molecular motors that crosslink and contract them. Recently it has been found that the assemblies of myosin II molecular motors in non-muscle cells are mixtures of fast and slow motor variants. Using computer simulations and an analytical mean field theory of a crossbridge model for myosin II motors, here we show that such motor ensembles effectively behave as active Maxwell elements. We calculate storage and loss moduli as a function of the model parameters and show that the rheological properties cross over from viscous to elastic as one increases the ratio of slow to fast motors. This suggests that cells tune their mechanical properties by regulating the composition of their myosin assemblies.


title: Lyman Continuum Escape Fraction from Low-mass Starbursts at z=1.3

abstract: We present a new constraint on the Lyman Continuum (LyC) escape fraction at z~1.3. We obtain deep, high sensitivity far-UV imaging with the Advanced Camera for Surveys (ACS) Solar Blind Channel (SBC) on the Hubble Space Telescope (HST), targeting 11 star-forming galaxies at 1.2<z<1.4. The galaxies are selected from the 3D-HST survey to have high H$\alpha$ equivalent width (EW) with EW > 190 \AA, low stellar mass (M* < 10^10 M_sun) and U-band magnitude of U<24.2. These criteria identify young, low metallicity star bursting populations similar to the primordial star-forming galaxies believed to have reionized the universe. We do not detect any LyC signal (with S/N >3) in the individual galaxies or in the stack in the far-UV images. We place $3\sigma$ limits on the relative escape fraction of individual galaxies to be f_{esc,rel}<[0.10-0.22] and a stacked $3\sigma$ limit of f_{esc,rel}<0.07. Comparing to the confirmed LyC emitters from the literature, the galaxies in our sample span similar ranges of various galaxy properties including stellar mass, dust attenuation, and star formation rate (SFR). In particular, we compare the distribution of H$\alpha$ and [OIII] EWs of confirmed LyC emitters and non-detections including the galaxies in this study. Finally, we discuss if a dichotomy seen in the distribution of H$\alpha$ EWs can perhaps distinguish the LyC emitters from the non-detections.


title: Representations for Stable Off-Policy Reinforcement Learning

abstract: Reinforcement learning with function approximation can be unstable and even divergent, especially when combined with off-policy learning and Bellman updates. In deep reinforcement learning, these issues have been dealt with empirically by adapting and regularizing the representation, in particular with auxiliary tasks. This suggests that representation learning may provide a means to guarantee stability. In this paper, we formally show that there are indeed nontrivial state representations under which the canonical TD algorithm is stable, even when learning off-policy. We analyze representation learning schemes that are based on the transition matrix of a policy, such as proto-value functions, along three axes: approximation error, stability, and ease of estimation. In the most general case, we show that a Schur basis provides convergence guarantees, but is difficult to estimate from samples. For a fixed reward function, we find that an orthogonal basis of the corresponding Krylov subspace is an even better choice. We conclude by empirically demonstrating that these stable representations can be learned using stochastic gradient descent, opening the door to improved techniques for representation learning with deep networks.


title: Community Network Auto-Regression for High-Dimensional Time Series

abstract: Modeling responses on the nodes of a large-scale network is an important task that arises commonly in practice. This paper proposes a community network vector autoregressive (CNAR) model, which utilizes the network structure to characterize the dependence and intra-community homogeneity of the high dimensional time series. The CNAR model greatly increases the flexibility and generality of the network vector autoregressive (Zhu et al, 2017, NAR) model by allowing heterogeneous network effects across different network communities. In addition, the non-community-related latent factors are included to account for unknown cross-sectional dependence. The number of network communities can diverge as the network expands, which leads to estimating a diverging number of model parameters. We obtain a set of stationary conditions and develop an efficient two-step weighted least-squares estimator. The consistency and asymptotic normality properties of the estimators are established. The theoretical results show that the two-step estimator improves the one-step estimator by an order of magnitude when the error admits a factor structure. The advantages of the CNAR model are further illustrated on a variety of synthetic and real datasets.


title: Quantum Secured Internet Transport

abstract: Quantum computing represents an emerging threat to the public key infrastructure underlying transport layer security (TLS) widely used in the Internet. This paper describes how QKD symmetric keys can be used with TLS to provide quantum computing resistant security for existing Internet applications. We also implement and test a general hybrid key delivery architecture with QKD over long distance fibers between secure sites, and wireless key distribution over short distance within each site Finally we show how this same capability can be extended to a TLS cipher scheme with perfect security.


title: Local Access to Sparse Connected Subgraphs Via Edge Sampling

abstract: We contribute an approach to the problem of locally computing sparse connected subgraphs of dense graphs. In this setting, given an edge in a connected graph $G = (V, E)$, an algorithm locally decides its membership in a sparse connected subgraph $G^* = (V, E^*)$, where $E^* \subseteq E$ and $|E^*| = o(|E|)$. Such an approach to subgraph construction is useful when dealing with massive graphs, where reading in the graph's full network description is impractical. While most prior results in this area require assumptions on $G$ or that $|E'| \le (1+\epsilon)|V|$ for some $\epsilon > 0$, we relax these assumptions. Given a general graph and a parameter $T$, we provide membership queries to a subgraph with $O(|V|T)$ edges using $\widetilde{O}(|E|/T)$ probes. This is the first algorithm to work on general graphs and allow for a tradeoff between its probe complexity and the number of edges in the resulting subgraph. We achieve this result with ideas motivated from edge sparsification techniques that were previously unused in this problem. We believe these techniques will motivate new algorithms for this problem and related ones. Additionally, we describe an efficient method to access any node's neighbor set in a sparsified version of $G$ where each edge is deleted with some i.i.d. probability.


title: Hydrodynamic gradient expansion in linear response theory

abstract: A foundational question in relativistic fluid mechanics concerns the properties of the hydrodynamic gradient expansion at large orders. We establish the precise conditions under which this gradient expansion diverges for a broad class of microscopic theories admitting a relativistic hydrodynamic limit, in the linear regime. Our result does not rely on highly symmetric fluid flows utilized by previous studies of heavy-ion collisions and cosmology. The hydrodynamic gradient expansion diverges whenever energy density or velocity fields have support in momentum space exceeding a critical momentum, and converges otherwise. This critical momentum is an intrinsic property of the microscopic theory and is set by branch point singularities of hydrodynamic dispersion relations.


title: Fermi Surface Geometry

abstract: Motivated by the famous and pioneering mathematical works by Perelman, Hamilton, and Thurston, we introduce the concept of using modern geometrical mathematical classifications of multi-dimensional manifolds to characterize electronic structures and predict non-trivial electron transport phenomena. Here we develop the Fermi Surface Geometry Effect (FSGE), using the concepts of tangent bundles and Gaussian curvature as an invariant. We develop an index, $\mathbb{H}_F$, for describing the the "hyperbolicity" of the Fermi Surface (FS) and show a universal correlation (R$^2$ = 0.97) with the experimentally measured intrinsic anomalous Hall effect of 16 different compounds spanning a wide variety of crystal, chemical, and electronic structure families, including where current methods have struggled. This work lays the foundation for developing a complete theory of geometrical understanding of electronic (and by extension magnonic and phononic) structure manifolds, beginning with Fermi surfaces. In analogy to the broad impact of topological physics, the concepts begun here will have far reaching consequences and lead to a paradigm shift in the understanding of electron transport, moving it to include geometrical properties of the E vs k manifold as well as topological properties.

categories: cond-mat.str-el cond-mat.mtrl-sci physics.comp-ph quant-ph


title: GW190814: Spin and equation of state of a neutron star companion

abstract: The recent discovery by LIGO/Virgo of a merging binary having a $\sim 23 M_\odot$ black hole and a $\sim 2.6 M_\odot$ compact companion has triggered a debate regarding the nature of the secondary, which falls into the so-called mass gap. Here we explore some consequences of the assumption that the secondary was a neutron star (NS). We show with concrete examples of heretofore viable equations of state (EOSs) that rapid uniform rotation may neither be necessary for some EOSs nor sufficient for others to explain the presence of a NS. Absolute upper limits for the maximum mass of a spherical NS derived from GW170817 already suggest that this unknown compact companion might be a slowly or even a nonrotating NS. However several soft NS EOSs favored by GW170817 with maximum spherical masses $\lesssim 2.1 M_\odot$ cannot be invoked to explain this object, even allowing for maximum uniform rotation. By contrast, sufficiently stiff EOSs that yield $2.6 M_\odot$ NSs which are slowly rotating or, in some cases, nonrotating, and are compatible with GW170817 and the results of NICER, can account for the black hole companion.


title: A singularly perturbed system of parabolic equations

abstract: The work is devoted to the construction of the asymptotic behavior of the solution of a singularly perturbed system of equations of parabolic type, in the case when the limit equation has a regular singularity as the small parameter tends to zero. The asymptotics of the solution of such problems contains boundary layer functions.


title: HAT-P-58b -- HAT-P-64b: Seven Planets Transiting Bright Stars

abstract: We report the discovery and characterization of 7 transiting exoplanets from the HATNet survey. The planets, which are hot Jupiters and Saturns transiting bright sun-like stars, include: HAT-P-58b (with mass Mp = 0.37 MJ, radius Rp = 1.33 RJ, and orbital period P = 4.0138 days), HAT-P-59b (Mp = 1.54 MJ, Rp = 1.12 RJ, P = 4.1420 days), HAT-P-60b (Mp = 0.57 MJ, Rp = 1.63 RJ, P = 4.7948 days), HAT-P-61b (Mp = 1.06 MJ, Rp = 0.90 RJ, P = 1.9023 days), HAT-P-62b (Mp = 0.76 MJ, Rp = 1.07 RJ, P = 2.6453 days), HAT-P-63b (Mp = 0.61 MJ, Rp = 1.12 RJ, P = 3.3777 days), and HAT-P-64b (Mp = 0.58 MJ, Rp = 1.70 RJ, P = 4.0072 days). The typical errors on these quantities are 0.06 MJ, 0.03 RJ, and 0.2seconds, respectively. We also provide accurate stellar parameters for each of the hosts stars. With V = 9.710+/-0.050mag, HAT-P-60 is an especially bright transiting planet host, and an excellent target for additional follow-up observations. With Rp = 1.703+/-0.070 RJ, HAT-P-64b is a highly inflated hot Jupiter around a star nearing the end of its main-sequence lifetime, and is among the largest known planets. Five of the seven systems have long-cadence observations by TESS which are included in the analysis. Of particular note is HAT-P-59 (TOI-1826.01) which is within the Northern continuous viewing zone of the TESS mission, and HAT-P-60, which is the TESS candidate TOI-1580.01.


title: Is there a Golden Parachute in Sannikov's principal-agent problem?

abstract: This paper provides a complete review of the continuous-time optimal contracting problem introduced by Sannikov, in the extended context allowing for possibly different discount rates for both parties. The agent's problem is to seek for optimal effort, given the compensation scheme proposed by the principal over a random horizon. Then, given the optimal agent's response, the principal determines the best compensation scheme in terms of running payment, retirement, and lump-sum payment at retirement. A Golden Parachute is a situation where the agent ceases any effort at some positive stopping time, and receives a payment afterwards, possibly under the form of a lump sum payment, or of a continuous stream of payments. We show that a Golden Parachute only exists in certain specific circumstances. This is in contrast with the results claimed by Sannikov, where the only requirement is a positive agent's marginal cost of effort at zero. Namely, we show that there is no Golden Parachute if this parameter is too large. Similarly, in the context of a concave marginal utility, there is no Golden Parachute if the agent's utility function has a too negative curvature at zero. In the general case, we prove that an agent with positive reservation utility is either never retired by the principal, or retired above some given threshold (as in Sannikov's solution). We show that different discount factors induce a face-lifted utility function, which allows to reduce the analysis to a setting similar to the equal discount rates one. Finally, we also confirm that an agent with small reservation utility may have an informational rent, meaning that the principal optimally offers him a contract with strictly higher utility than his participation value.


title: Magnetic effects on fields morphologies and reversals in geodynamo simulations

abstract: The dynamo effect is the most popular candidate to explain the non-primordial magnetic fields of astrophysical objects. Although many systematic studies of parameters have already been made to determine the different dynamical regimes explored by direct numerical geodynamo simulations, it is only recently that the regime corresponding to the outer core of the Earth characterized by a balance of forces between the Coriolis and Lorentz forces is accessible numerically. In most previous studies, the Lorentz force played a relatively minor role. For example, they have shown that a purely hydrodynamic parameter (the local Rossby number $Ro_\ell$ determines the stability domain of dynamos dominated by the axial dipole (dipolar dynamos). In this study, we show that this result cannot hold when the Lorentz force becomes dominant. We model turbulent geodynamo simulations with a strong Lorentz force by varying the important parameters over several orders of magnitude. This method enables us to question previous results and to argue on the applications of numerical dynamos in order to better understand the geodynamo problem. Strong dipolar fields considerably affect the kinetic energy distribution of convective motions which enables the maintenance of this field configuration. The relative importance of each force depends on the spatial length scale, whereas $Ro_\ell$ is a global output parameter which ignores the spatial dependency. We show that inertia does not induce a dipole collapse as long as the Lorentz and the Coriolis forces remain dominant at large length scales.


title: New generalisation of Jacobi's derivative formula

abstract: A stream of new theta relations is obtained. They follow from the general Thomae formula, which is a new result giving expressions for theta derivatives (the zero values of the lowest non-vanishing derivatives of theta functions with singular half-period characteristics) in terms of branch points and the period matrix of a hyperelliptic Riemann surface. The new theta relations contain (i) linear relations on the vector space of first order theta derivatives which are arranged in gradients, (ii) relations between second order theta derivatives and symmetric bilinear forms on the vector space of the gradients, (iii) relations between third order theta derivatives and symmetric trilinear forms on the vector space of the gradients, and (iv) a conjecture regarding higher order theta derivatives. It is shown how the Schottky identity (in the hyperelliptic case) is derived from the obtained relations.


title: Almost automorphically forced flows on $S^1$ or $\mathbb{R}$ in one-dimensional almost periodic semilinear heat equations

abstract: In this paper, we consider the asymptotic dynamics of the skew-product semiflow generated by the following time almost-periodically forced scalar reaction-diffusion equation \begin{equation}\label{eq0} u_{t}=u_{xx}+f(t,u,u_{x}),\,\,t>0,\, 0<x<L \end{equation} with periodic boundary condition \begin{equation} \label{bdc1} u(t,0)=u(t,L),\quad u_x(t,0)=u_x(t,L), \end{equation} where $f$ is uniformly almost periodic in $t$. In particular, we study the topological structure of the limit sets of the skew-product semiflow. It is proved that any compact minimal invariant set (throughout this paper, we refer to it as a minimal set) can be residually embedded into an invariant set of some almost automorphically-forced flow on a circle $S^1=\mathbb{R}/L\mathbb{Z}$. Particularly, if $f(t,u,p)=f(t,u,-p)$, then the flow on a minimal set topologically conjugates to an almost periodically-forced minimal flow on $\mathbb{R}$. Moreover, it is proved that the $\omega$-limit set of any bounded orbit contains at most two minimal sets that cannot be obtained from each other by phase translation. In addition, we further consider the asymptotic dynamics of the skew-product semiflow generated by \eqref{eq0} with Neumann boundary condition \begin{equation*} \label{bcd2} u_x(t,0)=u_x(t,L)=0, \end{equation*} or Dirichlet boundary condition \begin{equation*}\label{bdc3} u(t,0)=u(t,L)=0. \end{equation*} Under certain direct assumptions on $f$, it is proved in this paper that the flow on any minimal set of \eqref{eq0}, with Neumann boundary condition or Dirichlet boundary condition, topologically conjugates to an almost periodically-forced minimal flow on $\mathbb{R}$. Finally, a counterexample is given to show that even for quasi-periodic equations, the results we obtain here cannot be further improved in general.


title: ISINet: An Instance-Based Approach for Surgical Instrument Segmentation

abstract: We study the task of semantic segmentation of surgical instruments in robotic-assisted surgery scenes. We propose the Instance-based Surgical Instrument Segmentation Network (ISINet), a method that addresses this task from an instance-based segmentation perspective. Our method includes a temporal consistency module that takes into account the previously overlooked and inherent temporal information of the problem. We validate our approach on the existing benchmark for the task, the Endoscopic Vision 2017 Robotic Instrument Segmentation Dataset, and on the 2018 version of the dataset, whose annotations we extended for the fine-grained version of instrument segmentation. Our results show that ISINet significantly outperforms state-of-the-art methods, with our baseline version duplicating the Intersection over Union (IoU) of previous methods and our complete model triplicating the IoU.


title: Multi-Domain Image Completion for Random Missing Input Data

abstract: Multi-domain data are widely leveraged in vision applications taking advantage of complementary information from different modalities, e.g., brain tumor segmentation from multi-parametric magnetic resonance imaging (MRI). However, due to possible data corruption and different imaging protocols, the availability of images for each domain could vary amongst multiple data sources in practice, which makes it challenging to build a universal model with a varied set of input data. To tackle this problem, we propose a general approach to complete the random missing domain(s) data in real applications. Specifically, we develop a novel multi-domain image completion method that utilizes a generative adversarial network (GAN) with a representational disentanglement scheme to extract shared skeleton encoding and separate flesh encoding across multiple domains. We further illustrate that the learned representation in multi-domain image completion could be leveraged for high-level tasks, e.g., segmentation, by introducing a unified framework consisting of image completion and segmentation with a shared content encoder. The experiments demonstrate consistent performance improvement on three datasets for brain tumor segmentation, prostate segmentation, and facial expression image completion respectively.


title: Flow-Based Likelihoods for Non-Gaussian Inference

abstract: We investigate the use of data-driven likelihoods to bypass a key assumption made in many scientific analyses, which is that the true likelihood of the data is Gaussian. In particular, we suggest using the optimization targets of flow-based generative models, a class of models that can capture complex distributions by transforming a simple base distribution through layers of nonlinearities. We call these flow-based likelihoods (FBL). We analyze the accuracy and precision of the reconstructed likelihoods on mock Gaussian data, and show that simply gauging the quality of samples drawn from the trained model is not a sufficient indicator that the true likelihood has been learned. We nevertheless demonstrate that the likelihood can be reconstructed to a precision equal to that of sampling error due to a finite sample size. We then apply FBLs to mock weak lensing convergence power spectra, a cosmological observable that is significantly non-Gaussian (NG). We find that the FBL captures the NG signatures in the data extremely well, while other commonly used data-driven likelihoods, such as Gaussian mixture models and independent component analysis, fail to do so. This suggests that works that have found small posterior shifts in NG data with data-driven likelihoods such as these could be underestimating the impact of non-Gaussianity in parameter constraints. By introducing a suite of tests that can capture different levels of NG in the data, we show that the success or failure of traditional data-driven likelihoods can be tied back to the structure of the NG in the data. Unlike other methods, the flexibility of the FBL makes it successful at tackling different types of NG simultaneously. Because of this, and consequently their likely applicability across datasets and domains, we encourage their use for inference when sufficient mock data are available for training.


title: Revealing the intermediate-mass black hole at the heart of the dwarf galaxy NGC 404 with sub-parsec resolution ALMA observations

abstract: We estimate the mass of the intermediate-mass black hole at the heart of the dwarf elliptical galaxy NGC 404 using Atacama Large Millimeter/submillimeter Array (ALMA) observations of the molecular interstellar medium at an unprecedented linear resolution of ~0.5 pc, in combination with existing stellar kinematic information. These ALMA observations reveal a central disc/torus of molecular gas clearly rotating around the black hole. This disc is surrounded by a morphologically and kinematically complex flocculent distribution of molecular clouds, that we resolve in detail. Continuum emission is detected from the central parts of NGC 404, likely arising from the Rayleigh-Jeans tail of emission from dust around the nucleus, and potentially from dusty massive star-forming clumps at discrete locations in the disc. Several dynamical measurements of the black hole mass in this system have been made in the past, but they do not agree. We show here that both the observed molecular gas and stellar kinematics independently require a ~5x10$^5$ Msun black hole once we include the contribution of the molecular gas to the potential. Our best estimate comes from the high-resolution molecular gas kinematics, suggesting the black hole mass of this system is 5.5$^{+4.1}_{-3.8}\times$10$^5$ Msun (at the 99% confidence level), in good agreement with our revised stellar kinematic measurement and broadly consistent with extrapolations from the black hole mass - velocity dispersion and black hole mass - bulge mass relations. This highlights the need to accurately determine the mass and distribution of each dynamically important component around intermediate-mass black holes when attempting to estimate their masses.


title: APOGEE Data and Spectral Analysis from SDSS Data Release 16: Seven Years of Observations Including First Results from APOGEE-South

abstract: The spectral analysis and data products in Data Release 16 (DR16; December 2019) from the high-resolution near-infrared APOGEE-2/SDSS-IV survey are described. Compared to the previous APOGEE data release (DR14; July 2017), APOGEE DR16 includes about 200000 new stellar spectra, of which 100000 are from a new southern APOGEE instrument mounted on the 2.5 m du Pont telescope at Las Campanas Observatory in Chile. DR16 includes all data taken up to August 2018, including data released in previous data releases. All of the data have been re-reduced and re-analyzed using the latest pipelines, resulting in a total of 473307 spectra of 437445 stars. Changes to the analysis methods for this release include, but are not limited to, the use of MARCS model atmospheres for calculation of the entire main grid of synthetic spectra used in the analysis, a new method for filling "holes" in the grids due to unconverged model atmospheres, and a new scheme for continuum normalization. Abundances of the neutron capture element Ce are included for the first time. A new scheme for estimating uncertainties of the derived quantities using stars with multiple observations has been applied, and calibrated values of surface gravities for dwarf stars are now supplied. Compared to DR14, the radial velocities derived for this release more closely match those in the Gaia DR2 data base, and a clear improvement in the spectral analysis of the coolest giants can be seen. The reduced spectra as well as the result of the analysis can be downloaded using links provided in the SDSS DR16 web page.


title: On the Origin of Long-Lived Particles

abstract: MATHUSLA is a proposed large-volume displaced vertex (DV) detector, situated on the surface above CMS and designed to search for long-lived particles (LLPs) produced at the HL-LHC. We show that a discovery of LLPs at MATHUSLA would not only prove the existence of BSM physics, it would also uncover the theoretical origin of the LLPs, despite the fact that MATHUSLA gathers no energy or momentum information on the LLP decay products. Our analysis is simple and robust, making it easily generalizable to include more complex LLP scenarios, and our methods are applicable to LLP decays discovered in ATLAS, CMS, LHCb, or other external detectors. In the event of an LLP detection, MATHUSLA can act as a Level-1 trigger for the main detector, guaranteeing that the LLP production event is read out at CMS. We perform an LLP simplified model analysis to show that combining information from the MATHUSLA and CMS detectors would allow the LLP production mode topology to be determined with as few as $\sim 100$ observed LLP decays. Underlying theory parameters, like the LLP and parent particle masses, can also be measured with $\lesssim 10\%$ precision. Together with information on the LLP decay mode from the geometric properties of the observed DV, it is clear that MATHUSLA and CMS together will be able to characterize any newly discovered physics in great detail.


title: Higher-form Gauge Symmetries in Multipole Topological Phases

abstract: In this article we study field-theoretical aspects of multipolar topological insulators. Previous research has shown that such systems naturally couple to higher-rank tensor gauge fields that arise as a result of gauging dipole or subsystem $U(1)$ symmetries. Here we propose a complementary framework using electric higher-form symmetries. We utilize the fact that gauging 1-form electric symmetries results in a 2-form gauge field which couples naturally to extended line-like objects: Wilson lines. In our context the Wilson lines are electric flux lines associated to the electric polarization of the system. This allows us to define a generalized 2-form Peierls' substitution for dipoles that shows that the off-diagonal components of a rank-2 tensor gauge field $A_{ij}$ can arise as a lattice Peierls factor generated by the background antisymmetric 2-form gauge field. This framework has immediate applications: (i) it allows us to construct a manifestly topological quadrupolar response action given by a Dixmier-Douady invariant -- a generalization of a Chern number for 2-form gauge fields -- which makes plain the quantization of the quadrupole moment in the presence of certain crystal symmetries; (ii) it allows for a clearer interpretation of the rank-2 Berry phase calculation of the quadrupole moment; (iii) it allows for a proof of a generic Lieb-Schultz-Mattis theorem for dipole-conserving systems.


title: Distributed-Memory DMRG via Sparse and Dense Parallel Tensor Contractions

abstract: The Density Matrix Renormalization Group (DMRG) algorithm is a powerful tool for solving eigenvalue problems to model quantum systems. DMRG relies on tensor contractions and dense linear algebra to compute properties of condensed matter physics systems. However, its efficient parallel implementation is challenging due to limited concurrency, large memory footprint, and tensor sparsity. We mitigate these problems by implementing two new parallel approaches that handle block sparsity arising in DMRG, via Cyclops, a distributed memory tensor contraction library. We benchmark their performance on two physical systems using the Blue Waters and Stampede2 supercomputers. Our DMRG performance is improved by up to 5.9X in runtime and 99X in processing rate over ITensor, at roughly comparable computational resource use. This enables higher accuracy calculations via larger tensors for quantum state approximation. We demonstrate that despite having limited concurrency, DMRG is weakly scalable with the use of efficient parallel tensor contraction mechanisms.


title: Is negative kinetic energy meta-stable?

abstract: We explore the possibility that theories with negative kinetic energy (ghosts) can be meta-stable up to cosmologically long times. In classical mechanics, ghosts undergo spontaneous lockdown rather than run-away if weakly-coupled and non-resonant. Physical examples of this phenomenon are shown. In quantum mechanics this leads to meta-stability similar to vacuum decay. In classical field theory, lockdown is broken by resonances and ghosts behave statistically, drifting towards infinite entropy as no thermal equilibrium exists. We analytically and numerically compute the run-away rate finding that it is cosmologically slow in 4-derivative gravity, where ghosts have gravitational interactions only. In quantum field theory the ghost run-away rate is naively infinite in perturbation theory, analogously to what found in early attempts to compute vacuum tunnelling; we do not know the true rate.

categories: hep-th cond-mat.stat-mech gr-qc hep-ph quant-ph


title: Climate & BCG: Effects on COVID-19 Death Growth Rates

abstract: Multiple studies have suggested the spread of COVID-19 is affected by factors such as climate, BCG vaccinations, pollution and blood type. We perform a joint study of these factors using the death growth rates of 40 regions worldwide with both machine learning and Bayesian methods. We find weak, non-significant (< 3$\sigma$) evidence for temperature and relative humidity as factors in the spread of COVID-19 but little or no evidence for BCG vaccination prevalence or $\text{PM}_{2.5}$ pollution. The only variable detected at a statistically significant level (>3$\sigma$) is the rate of positive COVID-19 tests, with higher positive rates correlating with higher daily growth of deaths.


title: Fermi arcs and pseudogap in a lattice model of a doped orthogonal metal

abstract: Since the discovery of the pseudogap and Fermi arc states in underdoped cuprates, the understanding of such non-Fermi-liquid states and the associated violation of Luttinger's theorem have been the central theme in correlated electron systems. However, still lacking is a well-accepted theoretical framework to unambiguously explain these metallic states that are clearly beyond Landau's Fermi liquid and Luttinger's theorem of a Fermi surface and electron filling. Here, we design a lattice model of orthogonal metals with fermion and Ising matter fields coupled to topological order and, by solving the model via unbiased quantum Monte Carlo simulation at generic electron fillings, find that the system gives birth to phenomena of the Fermi arc and pseudogap in the single-particle spectrum that go beyond the Luttinger sum rule with broken Fermi surface but no symmetry breaking. The pseudogap and Fermi arcs coexist with a background of a deconfined Z2 gauge field, and we further find that the confinement transition of the gauge field triggers a superconductivity instability and that the hopping of the gauge-neutral fermions brings the "large" Fermi surface back from the Fermi arc state. Our unbiased numerical results provide a concrete model realization and theoretical framework for the coupling between gauge field and fermions and, in the process, generate the rich phenomena of the pseudogap, the Fermi arc, and superconductivity in generic correlated electron systems.


title: Less than the sum of its parts: the dust-corrected H$\alpha$ luminosity of star-forming galaxies explored at different spatial resolutions with MaNGA and MUSE

abstract: The H$\alpha$ and H$\beta$ emission line luminosities measured in a single integrated spectrum are affected in non-trivial ways by point-to-point variations in dust attenuation in a galaxy. This work investigates the impact of this variation when estimating global H$\alpha$ luminosities corrected for the presence of dust by a global Balmer decrement. Analytical arguments show that the dust-corrected H$\alpha$ luminosity is always underestimated when using the global H$\alpha$/H$\beta$ flux ratio to correct for dust attenuation. We measure this effect on 156 face-on star-forming galaxies from the Mapping Nearby Galaxies at APO (MaNGA) survey. At 1-2 kpc spatial resolution, the effect is small but systematic, with the integrated dust-corrected H$\alpha$ luminosity underestimated by $2$-$4$ per cent (and typically not more than by $10$ per cent), and depends on the specific star formation rate of the galaxy. Given the spatial resolution of MaNGA, these are lower limits for the effect. From Multi Unit Spectroscopic Explorer (MUSE) observations of NGC 628 with a resolution of 36 pc we find the discrepancy between the globally and the point-by-point dust-corrected H$\alpha$ luminosity to be $14 \pm 1$ per cent, which may still underestimate the true effect. We use toy models and simulations to show that the true difference depends strongly on the spatial variance of the H$\alpha$/H$\beta$ flux ratio, and on the slope of the relation between H$\alpha$ luminosity and dust attenuation within a galaxy. Larger samples of higher spatial resolution observations are required to quantify the dependence of this effect as a function of galaxy properties.


title: Dissecting and Modelling Galaxy Assembly Bias

abstract: Understanding the galaxy-halo connection is fundamental for contemporary models of galaxy clustering. The extent to which the haloes' assembly history and environment impact galaxy clustering (a.k.a. galaxy assembly bias; GAB), remains a complex and challenging problem. Using a semi-analytic galaxy formation model, we study the individual contributions of different secondary halo properties to the GAB signal. These are obtained by comparing the clustering of stellar-mass selected samples to that of shuffled samples where the galaxies are randomly reassigned to haloes of fixed mass and a specified secondary halo property. We explore a large range of internal halo properties and environmental measures. We find that commonly used properties like halo age or concentration amount to only 20-30 per cent of the signal, while the smoothed matter density or the tidal anisotropy can account for the full level of GAB (though care should be given to the specific definition). For the "successful" measures, we examine the occupancy variations and the associated changes in the halo occupation function parameters. These are used to create mock catalogues that reproduce the full level of GAB. Finally, we propose a practical modification of the standard halo occupation distribution model, which can be tuned to any level of assembly bias. Fitting the parameters to our semi-analytic model, we demonstrate that the corresponding mock catalogue recovers the target level of GAB as well as the occupancy variations. Our results enable producing realistic mock catalogues and directly inform theoretical modelling of assembly bias and attempts to detect it in the Universe.


title: International Coordination of Multi-Messenger Transient Observations in the 2020s and Beyond: Kavli-IAU White Paper

abstract: This White Paper summarizes the discussions from a five-day workshop, involving 50 people from 18 countries, held in Cape Town, South Africa in February 2020. Convened by the International Astronomical Union's Executive Committee Working Group on Global Coordination of Ground and Space Astrophysics and sponsored by the Kavli Foundation, we discussed existing and potential bottlenecks for transient and multi-messenger astronomy, identifying eight broad areas of concern. Some of these are very similar to the challenges faced by many astronomers engaging in international collaboration, for example, data access policies, funding, theoretical and computational resources and workforce equity. Others, including, alerts, telescope coordination and target-of-opportunity implementation, are strongly linked to the time domain and are particularly challenging as we respond to transients. To address these bottlenecks we offer thirty-five specific recommendations, some of which are simply starting points and require development. These recommendations are not only aimed at collaborative groups and individuals, but also at the various organizations who are essential to making transient collaborations efficient and effective: including the International Astronomical Union, observatories, projects, scientific journals and funding agencies. We hope those involved in transient research will find them constructive and use them to develop collaborations with greater impact and more inclusive teams.

categories: astro-ph.IM astro-ph.GA astro-ph.HE astro-ph.SR