title: Fermionic criticality is shaped by Fermi surface topology: a case study of the Tomonaga-Luttinger liquid

abstract: We perform a unitary renormalization group (URG) study of the 1D fermionic Hubbard model. The formalism generates a family of effective Hamiltonians and many-body eigenstates arranged holographically across the tensor network from UV to IR. The URG is realized as a quantum circuit, leading to the entanglement holographic mapping (EHM) tensor network description. A topological $\Theta$-term of the projected Hilbert space of the degrees of freedom at the Fermi surface are shown to govern the nature of RG flow towards either the gapless Tomonaga-Luttinger liquid or gapped quantum liquid phases. This results in a nonperturbative version of the Berezenskii-Kosterlitz-Thouless (BKT) RG phase diagram, revealing a line of intermediate coupling stable fixed points, while the nature of RG flow around the critical point is identical to that obtained from the weak-coupling RG analysis. This coincides with a phase transition in the many-particle entanglement, as the entanglement entropy RG flow shows distinct features for the critical and gapped phases depending on the value of the topological $\Theta$-term. We demonstrate the Ryu-Takyanagi entropy bound for the many-body eigenstates comprising the EHM network, concretizing the relation to the holographic duality principle. The scaling of the entropy bound also distinguishes the gapped and gapless phases, implying the generation of very different holographic spacetimes across the critical point. Finally, we treat the Fermi surface as a quantum impurity coupled to the high energy electronic states. A thought-experiment is devised in order to study entanglement entropy generated by isolating the impurity, and propose ways by which to measure it by studying the quantum noise and higher order cumulants of the full counting statistics.


title: Analytic gradients for spin multiplets in natural orbital functional theory

abstract: Analytic energy gradients with respect to nuclear motion are derived for non-singlet compounds in the natural orbital functional theory. We exploit the formulation for multiplets in order to obtain a simple formula valid for any many-electron system in its ground mixed state with a total spin S and all possible spin projection Sz values. We demonstrate that the analytic gradients can be obtained without resorting to linear response theory or involving iterative procedures. It is required a single evaluation, so integral derivatives can be computed on-the-fly along the calculation and thus improve the effectiveness of screening by the Schwarz inequality. Results for small and medium size molecules with many spin multiplicities are shown. Our results are compared with experimental data and accurate theoretical equilibrium geometries.


title: Rare decays of the top quark mediated by Z' gauge bosons and flavor violation

abstract: The rare top quark decays mediated by a new neutral massive gauge boson that is predicted in models with extended gauge symmetries are studied. We focus on the processes $t\to cV, uV$ induced at the one loop level, where $V =\gamma, g$, by considering different extended models. It is found that, within a broad range of mass of the new neutral gauge boson, the models predict branching ratios for the decays in study that are competitive with respect to the corresponding branching ratios in the standard model. In order to establish bound on our branching ratios, we consider the recent experimental bounds as $m_{Z^\prime}\geq$ 3.8-4.5 TeV, depending on the model, which also impose restrictions on our calculation. Even in this case, the resulting branching ratios are of the same order of magnitude as that predicted by the standard model. It should be noted that for the case of two models studied here, since no experimental bound exists to compare with, the results could be important, as they are, in the best of cases, two orders of magnitude larger than the predicted by the standard model.


title: Don't Explain without Verifying Veracity: An Evaluation of Explainable AI with Video Activity Recognition

abstract: Explainable machine learning and artificial intelligence models have been used to justify a model's decision-making process. This added transparency aims to help improve user performance and understanding of the underlying model. However, in practice, explainable systems face many open questions and challenges. Specifically, designers might reduce the complexity of deep learning models in order to provide interpretability. The explanations generated by these simplified models, however, might not accurately justify and be truthful to the model. This can further add confusion to the users as they might not find the explanations meaningful with respect to the model predictions. Understanding how these explanations affect user behavior is an ongoing challenge. In this paper, we explore how explanation veracity affects user performance and agreement in intelligent systems. Through a controlled user study with an explainable activity recognition system, we compare variations in explanation veracity for a video review and querying task. The results suggest that low veracity explanations significantly decrease user performance and agreement compared to both accurate explanations and a system without explanations. These findings demonstrate the importance of accurate and understandable explanations and caution that poor explanations can sometimes be worse than no explanations with respect to their effect on user performance and reliance on an AI system.


title: Inverse scattering transform for the focusing nonlinear Schrodinger equation with counterpropagating flows

abstract: The inverse scattering transform for the focusing nonlinear Schrodinger equation is presented for a general class of initial conditions whose asymptotic behavior at infinity consists of counterpropagating waves. The formulation takes into account the branched nature of the two asymptotic eigenvalues of the associated scattering problem. The Jost eigenfunctions and scattering coefficients are defined explicitly as single-valued functions on the complex plane with jump discontinuities along certain branch cuts. The analyticity properties, symmetries, discrete spectrum, asymptotics and behavior at the branch points are discussed explicitly. The inverse problem is formulated as a matrix Riemann-Hilbert problem with poles. Reductions to all cases previously discussed in the literature are explicitly discussed. The scattering data associated to a few special cases consisting of physically relevant Riemann problems are explicitly computed.


title: Heuristics in experiments with infinitely large strategy spaces

abstract: We introduce a new methodology that enables detection of the onset of convergence towards Nash equilibria in simple repeated games with infinitely large strategy spaces, thereby revealing the heuristics used in decision-making. The method works by constraining on a special finite subset of strategies, called decoupled strategies. We show how the technique can be applied to understand price formation in financial market experiments by introducing a predictive measure {\Delta}D: the different between positive decoupled strategies (recommending to buy) and negative decoupled strategies (recommending to sell). Using {\Delta}D we illustrate how the method can predict (at certain special times) participants' actions with a high success rate in a series of experiments


title: Lim Ulrich sequences and Lech's conjecture

abstract: The long standing Lech's conjecture in commutative algebra states that for a flat local extension $(R,\mathfrak{m})\to (S,\mathfrak{n})$ of Noetherian local rings, we have an inequality on the Hilbert--Samuel multiplicities: $e(R)\leq e(S)$. In general the conjecture is wide open when $\dim R>3$, even in equal characteristic. In this paper, we prove Lech's conjecture in all dimensions, provided $(R,\mathfrak{m})$ is a standard graded ring over a perfect field localized at the homogeneous maximal ideal. We introduce the notions of lim Ulrich and weakly lim Ulrich sequences. Roughly speaking these are sequences of finitely generated modules that are not necessarily Cohen--Macaulay, but asymptotically behave like Ulrich modules. We prove that the existence of these sequences imply Lech's conjecture. Though the existence of Ulrich modules is known in very limited cases, we construct weakly lim Ulrich sequences for all standard graded domains over perfect fields of positive characteristic.


title: A Hierarchical Approach to Multi-Energy Demand Response: From Electricity to Multi-Energy Applications

abstract: Due to proliferation of energy efficiency measures and availability of the renewable energy resources, traditional energy infrastructure systems (electricity, heat, gas) can no longer be operated in a centralized manner under the assumption that consumer behavior is inflexible, i.e. cannot be adjusted in return for an adequate incentive. To allow for a less centralized operating paradigm, consumer-end perspective and abilities should be integrated in current dispatch practices and accounted for in switching between different energy sources not only at the system but also at the individual consumer level. Since consumers are confined within different built environments, this paper looks into an opportunity to control energy consumption of an aggregation of many residential, commercial and industrial consumers, into an ensemble. This ensemble control becomes a modern demand response contributor to the set of modeling tools for multi-energy infrastructure systems.


title: Automorphic Schwarzian equations and integrals of weight 2 forms

abstract: In this paper, we investigate the non-modular solutions to the Schwarz differential equation $\{f,\tau \}=sE_4(\tau)$ where $E_4(\tau)$ is the weight 4 Eisenstein series and $s$ is a complex parameter. In particular, we provide explicit solutions for each $s=2\pi^2(n/6)^2$ with $n\equiv 1\mod 12$. These solutions are obtained as integrals of meromorphic weight 2 modular forms. As a consequence, we find explicit solutions to the differential equation $\displaystyle y''+\frac{\pi^2n^2}{36}\,E_4\,y=0$ for each $n\equiv 1\mod 12$ generalizing the work of Hurwitz and Klein on the case $n=1$. Our investigation relies on the theory of equivariant functions on the complex upper half-plane. This paper supplements a previous work where we determine all the parameters $s$ for which the above Schwarzian equation has a modular solution.


title: Coverage and Deployment Analysis of Narrowband Internet of Things in the Wild

abstract: Narrowband Internet of Things (NB-IoT) is gaining momentum as a promising technology for massive Machine Type Communication (mMTC). Given that its deployment is rapidly progressing worldwide, measurement campaigns and performance analyses are needed to better understand the system and move toward its enhancement. With this aim, this paper presents a large scale measurement campaign and empirical analysis of NB-IoT on operational networks, and discloses valuable insights in terms of deployment strategies and radio coverage performance. The reported results also serve as examples showing the potential usage of the collected dataset, which we make open-source along with a lightweight data visualization platform.


title: Heuristic-Based Weak Learning for Automated Decision-Making

abstract: Machine learning systems impact many stakeholders and groups of users, often disparately. Prior studies have reconciled conflicting user preferences by aggregating a high volume of manually labeled pairwise comparisons, but this technique may be costly or impractical. How can we lower the barrier to participation in algorithm design? Instead of creating a simplified labeling task for a crowd, we suggest collecting ranked decision-making heuristics from a focused sample of affected users. With empirical data from two use cases, we show that our weak learning approach, which requires little to no manual labeling, agrees with participants' pairwise choices nearly as often as fully supervised approaches.


title: Observing Left-Right Symmetry in the Cosmic Microwave Background

abstract: We consider the possibility of probing left-right symmetric model (LRSM) via cosmic microwave background (CMB). We adopt the minimal LRSM with Higgs doublets, also known as the doublet left-right model (DLRM), where all fermions including the neutrinos acquire masses only via their couplings to the Higgs bidoublet. Due to the Dirac nature of light neutrinos, there exist additional relativistic degrees of freedom which can thermalise in the early universe by virtue of their gauge interactions corresponding to the right sector. We constrain the model from Planck 2018 bound on the effective relativistic degrees of freedom and also estimate the prospects for planned CMB Stage IV experiments to constrain the model further. We find that $W_R$ boson mass below 4.06 TeV can be ruled out from Planck 2018 bound at $2\sigma$ CL in the exact left-right symmetric limit which is equally competitive as the LHC bounds from dijet resonance searches. On the other hand Planck 2018 bound at $1\sigma$ CL can rule out a much larger parameter space out of reach of present direct search experiments, even in the presence of additional relativistic degrees of freedom around the TeV corner. We also study the consequence of these constraints on dark matter in DLRM by considering a right handed real fermion quintuplet to be the dominant dark matter component in the universe.


title: Cubic forms, anomaly cancellation and modularity

abstract: Motivated by the cubic forms and anomaly cancellation formulas of Witten-Freed-Hopkins, we give some new cubic forms on spin, spin$^c$, spin$^{w_2}$ and orientable 12-manifolds respectively. We relate them to $\eta$-invariants when the manifolds are with boundary, and mod 2 indices on 10 dimensional characteristic submanifolds when the manifolds are spin$^c$ or spin$^{w_2}$. Our method of producing these cubic forms is a combination of (generalized) Witten classes and the character of the basic representation of affine $E_8$.


title: Stability of dynamical topology against dynamical noise in quantum quenches

abstract: Equilibrium topological phases are robust against weak static disorder but may break down in the strong disorder regime. Here we explore the stability of the quench-induced emergent dynamical topology in the presence of dynamical noise. We develop an analytic theory and show that for weak noise, the quantum dynamics induced by quenching an initial trivial phase to Chern insulating regime exhibits robust emergent topology on certain momentum subspaces called band inversion surfaces (BISs). The dynamical topology is protected by the minimal oscillation frequency over the BISs, mimicking a bulk gap of the dynamical phase. Singularities emerge in the quench dynamics, with the minimal oscillation frequency vanishing on the BISs if increasing noise to critical strength, manifesting a dynamical topological transition, beyond which the emergent topology breaks down. Two types of dynamical transitions are predicted. Interestingly, we predict a sweet spot in the critical transition when noise couples to all three spin components in the same strength, in which case the dynamical topology survives at arbitrarily strong noise regime. This work unveils novel features of the dynamical topology under dynamical noise, which can be probed with control in experiment.


title: Maximal subgroups of non-torsion Grigorchuk-Gupta-Sidki groups

abstract: A Grigorchuk-Gupta-Sidki (GGS-)group is a subgroup of the automorphism group of the $p$-adic tree for an odd prime $p$, generated by one rooted automorphism and one directed automorphism. Pervova proved that all torsion GGS-groups do not have maximal subgroups of infinite index. Here we extend the result to non-torsion GGS-groups, which include the weakly regular branch, but not branch, GGS-group.


title: Differential Machine Learning

abstract: Differential machine learning combines automatic adjoint differentiation (AAD) with modern machine learning (ML) in the context of risk management of financial Derivatives. We introduce novel algorithms for training fast, accurate pricing and risk approximations, online, in real-time, with convergence guarantees. Our machinery is applicable to arbitrary Derivatives instruments or trading books, under arbitrary stochastic models of the underlying market variables. It effectively resolves computational bottlenecks of Derivatives risk reports and capital calculations. Differential ML is a general extension of supervised learning, where ML models are trained on examples of not only inputs and labels but also differentials of labels wrt inputs. It is also applicable in many situations outside finance, where high quality first-order derivatives wrt training inputs are available. Applications in Physics, for example, may leverage differentials known from first principles to learn function approximations more effectively. In finance, AAD computes pathwise differentials with remarkable efficacy so differential ML algorithms provide extremely effective pricing and risk approximations. We can produce fast analytics in models too complex for closed form solutions, extract the risk factors of complex transactions and trading books, and effectively compute risk management metrics like reports across a large number of scenarios, backtesting and simulation of hedge strategies, or regulations like XVA, CCR, FRTB or SIMM-MVA. TensorFlow implementation is available on https://github.com/differential-machine-learning


title: A Redundancy-Guided Approach for the Hazard Analysis of Digital Instrumentation and Control Systems in Advanced Nuclear Power Plants

abstract: Digital instrumentation and control (I&C) upgrades are a vital research area for nuclear industry. Despite their performance benefits, deployment of digital I&C in nuclear power plants (NPPs) has been limited. Digital I&C systems exhibit complex failure modes including common cause failures (CCFs) which can be difficult to identify. This paper describes the development of a redundancy-guided application of the Systems-Theoretic Process Analysis (STPA) and Fault Tree Analysis (FTA) for the hazard analysis of digital I&C in advanced NPPs. The resulting Redundancy-guided System-theoretic Hazard Analysis (RESHA) is applied for the case study of a representative state-of-the-art digital reactor trip system. The analysis qualitatively and systematically identifies the most critical CCFs and other hazards of digital I&C systems. Ultimately, RESHA can help researchers make informed decisions for how, and to what degree, defensive measures such as redundancy, diversity, and defense-in-depth can be used to mitigate or eliminate the potential hazards of digital I&C systems.


title: $(1+\varepsilon)$-moments suffice to characterise the GFF

abstract: We show that there is "no stable free field of index $\alpha\in (1,2)$", in the following sense. It was proved in a previous work by the authors, that subject to a \emph{fourth moment assumption}, any random generalised function on a domain $D$ of the plane, satisfying conformal invariance and a natural domain Markov property, must be a constant multiple of the Gaussian free field. In this article we show that the existence of $(1+\varepsilon)$-moments is sufficient for the same conclusion. A key idea is a new way of exploring the field, where (instead of looking at the more standard circle averages) we start from the boundary and discover averages of the field with respect to a certain "hitting density" of It\^o excursions.


title: Quantum mean field games

abstract: Quantum games represent the really 21st century branch of game theory, tightly linked to the modern development of quantum computing and quantum technologies. The main accent in these developments so far was made on stationary or repeated games. In the previous paper of the author the truly dynamic quantum game theory was initiated with strategies chosen by players in real time. Since direct continuous observations are known to destroy quantum evolutions (so-called quantum Zeno paradox) the necessary new ingredient for quantum dynamic games represented the theory of non-direct observations and the corresponding quantum filtering. Another remarkable 21st century branch of game theory represent the so-called mean-field games (MFG), with impressive and ever growing development. In this paper we are merging these two exciting new branches of game theory. Building a quantum analog of MFGs requires the full reconstruction of its foundations and methodology, because in $N$-particle quantum evolution particles are not separated in individual dynamics and the key concept of the classical MFG theory, the empirical measure defined as the sum of Dirac masses of the positions of the players, is not applicable in quantum setting. As a preliminary result we derive the new nonlinear stochastic Schr\"odinger equation, as the limit of continuously observed and controlled system of large number of interacting quantum particles, the result that may have an independent value. We then show that to a control quantum system of interacting particles there corresponds a special system of classical interacting particles with the identical limiting MFG system, defined on an appropriate Riemanian manifold. Solutions of this system are shown to specify approximate Nash equilibria for $N$-agent quantum games.


title: Defining an intrinsic stickiness parameter of stock price returns

abstract: We introduce a non linear pricing model of individual stock returns that defines a stickiness parameter of the returns. The pricing model resembles the capital asset pricing model used in finance but has a non linear component inspired from models of earth quake tectonic plate movements. The link to tectonic plate movements happens, since price movements of a given stock index is seen adding stress to its components of individual stock returns, in order to follow the index. How closely individual stocks follow the indexs price movements, can then be used to define their stickiness


title: Surface Pressure of Liquid Interfaces Laden with Micron-Sized Particles

abstract: We consider the surface pressure of a colloid-laden liquid interface. As micron-sized particles of suitable wettability can be irreversibly bound to the liquid interface on experimental timescales, we use the canonical ensemble to derive an expression for the surface pressure of a colloid-laden interface. We use this expression to show that adsorption of particles with only hard-core interactions has a negligible effect on surface pressures from typical Langmuir-trough measurements. Moreover, we show that Langmuir-trough measurements cannot be used to extract typical interparticle potentials. Finally, we argue that the dependence of measured surface pressure on surface fraction can be explained by particle coordination number at low to intermediate particle surface fractions. At high surface fractions, where the particles are jammed and cannot easily rearrange, contact-line sliding and/or deformations of the liquid interface at the length scale of the particles play a pivotal role.


title: Analysis and Prediction of COVID-19 Pandemic in Pakistan using Time-dependent SIR Model

abstract: The current outbreak is known as Coronavirus Disease or COVID-19 caused by the virus SAR-COV-2 which continues to wreak havoc across the globe. The World Health Organization (WHO) has declared the outbreak a Public Health Emergency of International Concern. In Pakistan, the spread of the virus is on the rise with the number of infected people and causalities rapidly increasing. In the absence of proper vaccination and treatment, to reduce the number of infections and casualties, the only option so far is to educate people regarding preventive measures and to enforce countrywide lock-down. Any strategy about the preventive measures needs to be based upon detailed analysis of the COVID-19 outbreak and accurate scientific predictions. In this paper, we conduct mathematical and numerical analysis to come up with reliable and accurate predictions of the outbreak in Pakistan. The time-dependent Susceptible-Infected-Recovered (SIR) model is used to fit the data and provide future predictions. The turning point of the peak of the pandemic is defined as the day when the transmission rate becomes less than the recovering rate. We have predicted that the outbreak will reach its maximum peak occurring from late May to 9 June with unrecovered number of Infectives in the range 20000-47000 and the cumulative number of infected cases in the range of 57500-153100. The number of Infectives will remain at the lower end in the lock-down scenario but can rapidly double or triple if the spread of the epidemic is not curtailed and localized. The uncertainty on single day projection in our analysis after April 15 is found to be within 5\%.


title: It's Easier to Translate out of English than into it: Measuring Neural Translation Difficulty by Cross-Mutual Information

abstract: The performance of neural machine translation systems is commonly evaluated in terms of BLEU. However, due to its reliance on target language properties and generation, the BLEU metric does not allow an assessment of which translation directions are more difficult to model. In this paper, we propose cross-mutual information (XMI): an asymmetric information-theoretic metric of machine translation difficulty that exploits the probabilistic nature of most neural machine translation models. XMI allows us to better evaluate the difficulty of translating text into the target language while controlling for the difficulty of the target-side generation component independent of the translation task. We then present the first systematic and controlled study of cross-lingual translation difficulties using modern neural translation systems. Code for replicating our experiments is available online at https://github.com/e-bug/nmt-difficulty.


title: Atomic Hydrogen Clues to the Formation of Counterrotating Stellar Discs

abstract: We present interferometric HI observations of six double-disc stellar counterrotator ("2$\sigma$") galaxies from the Atlas3D sample. Three are detected in HI emission; two of these are new detections. NGC 7710 shows a modestly asymmetric HI disc, and the atomic gas in PGC 056772 is centrally peaked but too poorly resolved to identify the direction of rotation. IC 0719, the most instructive system in this study, shows an extended, strongly warped disc of 43 kpc diameter, with a faint tail extending towards its neighbor IC 0718. The gas has likely been accreted from this external source during an encounter whose geometry directed the gas into misaligned retrograde orbits (with respect to the primary stellar body of IC 0719). In the interior, where dynamical time-scales are shorter, the HI has settled into the equatorial plane forming the retrograde secondary stellar disc. This is the first direct evidence that a double-disc stellar counterrotator could be formed through the accretion of retrograde gas. However, the dominant formation pathway for the formation of $2\sigma$ galaxies is still unclear. The Atlas3D sample shows some cases of the retrograde accretion scenario and also some cases in which a scenario based on an unusually well-aligned merger is more likely.


title: Manifold Proximal Point Algorithms for Dual Principal Component Pursuit and Orthogonal Dictionary Learning

abstract: We consider the problem of maximizing the $\ell_1$ norm of a linear map over the sphere, which arises in various machine learning applications such as orthogonal dictionary learning (ODL) and robust subspace recovery (RSR). The problem is numerically challenging due to its nonsmooth objective and nonconvex constraint, and its algorithmic aspects have not been well explored. In this paper, we show how the manifold structure of the sphere can be exploited to design fast algorithms for tackling this problem. Specifically, our contribution is threefold. First, we present a manifold proximal point algorithm (ManPPA) for the problem and show that it converges at a sublinear rate. Furthermore, we show that ManPPA can achieve a quadratic convergence rate when applied to the ODL and RSR problems. Second, we propose a stochastic variant of ManPPA called StManPPA, which is well suited for large-scale computation, and establish its sublinear convergence rate. Both ManPPA and StManPPA have provably faster convergence rates than existing subgradient-type methods. Third, using ManPPA as a building block, we propose a new approach to solving a matrix analog of the problem, in which the sphere is replaced by the Stiefel manifold. The results from our extensive numerical experiments on the ODL and RSR problems demonstrate the efficiency and efficacy of our proposed methods.


title: Sub-Image Anomaly Detection with Deep Pyramid Correspondences

abstract: Nearest neighbor (kNN) methods utilizing deep pre-trained features exhibit very strong anomaly detection performance when applied to entire images. A limitation of kNN methods is the lack of segmentation map describing where the anomaly lies inside the image. In this work we present a novel anomaly segmentation approach based on alignment between an anomalous image and a constant number of the similar normal images. Our method, Semantic Pyramid Anomaly Detection (SPADE) uses correspondences based on a multi-resolution feature pyramid. SPADE is shown to achieve state-of-the-art performance on unsupervised anomaly detection and localization while requiring virtually no training time.


title: Coulomb-blockade-controlled single-electron point source

abstract: Coulomb blockade is a fundamental phenomenon in physics enabling transfer of individual electrons one by one into electrically isolated nanostructures such as nanowires or quantum dots and thereby creation of sources of single electrons. Nowadays, solid-state single-electron sources are key elements of the emerging new technologies of quantum information processing and single-electron electronics. Moreover, advanced research in free-electron quantum optics and developments in electron microscopy require the point sources of free electrons in vacuum, which can be controlled on a single electron level. However, up to now, single-electron vacuum guns were not realized in practice. The problems to be solved include the formation of stable tip-shaped heterostructured emitters and control of liberated electrons in time and energy domain. Here we overcome these challenges by creating a field emission (FE) electron source based on a carbon nanowire coupled to an ultra-sharp diamond tip by a tunnel junction. Using energy spectroscopy, we directly observe Coulomb oscillations of the electron Fermi level in the nanowire at room temperature and FE currents up to 1 $\mu$A. We reveal that the oscillations are suppressed at high FE current either by Joule heating or due to the high tunneling resistance, depending on the nanowire charging time, ranging from about 25 fs to 1 ps. We anticipate that the combination of introduced carbon single-electron sources with laser-induced gating is highly promising for the creation of coherent ultrashort free-electron bunches of interest for low-energy electron holography and ultrafast electron or X-ray imaging and spectroscopy.


title: Classification-Based Anomaly Detection for General Data

abstract: Anomaly detection, finding patterns that substantially deviate from those seen previously, is one of the fundamental problems of artificial intelligence. Recently, classification-based methods were shown to achieve superior results on this task. In this work, we present a unifying view and propose an open-set method, GOAD, to relax current generalization assumptions. Furthermore, we extend the applicability of transformation-based methods to non-image data using random affine transformations. Our method is shown to obtain state-of-the-art accuracy and is applicable to broad data types. The strong performance of our method is extensively validated on multiple datasets from different domains.


title: Combinatorial and accessible weak model categories

abstract: In a previous work, we have introduced a weakening of Quillen model categories called weak model categories. They still allow all the usual constructions of model category theory, but are easier to construct and are in some sense better behaved. In this paper we continue to develop their general theory by introducing combinatorial and accessible weak model categories. We give simple necessary and sufficient conditions under which such a weak model category can be extended into a left and/or right semi-model category. As an application, we recover Cisinski-Olschok theory and generalize it to weak and semi-model categories. We also provide general existence theorems for both left and right Bousfield localization of combinatorial and accessible weak model structures, which combined with the results above gives existence results for left and right Bousfield localization of combinatorial and accessible left and right semi-model categories, generalizing previous results of Barwick. Surprisingly,we show that any left or right Bousfield localization of an accessible or combinatorial Quillen model category always exists, without properness assumptions, and is simultaneously both a left and a right semi-model category, without necessarily being a Quillen model category itself.


title: An obscured AGN population hidden in the VIPERS galaxies: identification through spectral energy distribution decomposition

abstract: The detection of X-ray emission constitutes a reliable and efficient tool for the selection of Active Galactic Nuclei (AGNs), although it may be biased against the most heavily absorbed AGNs. Simple mid-IR broad-band selection criteria identify a large number of luminous and absorbed AGNs, yet again host contamination could lead to non-uniform and incomplete samples. Spectral Energy Distribution (SED) decomposition is able to decouple the emission from the AGN versus that from star-forming regions, revealing weaker AGN components. We aim to identify the obscured AGN population in the VIPERS survey in the CFHTLS W1 field through SED modelling. We construct SEDs for 6,860 sources and identify 160 AGNs at a high confidence level using a Bayesian approach. Using optical spectroscopy, we confirm the nature of ~85% of the AGNs. Our AGN sample is highly complete (~92%) compared to mid-IR colour selected AGNs, including a significant number of galaxy-dominated systems with lower luminosities. In addition to the lack of X-ray emission (80%), the SED fitting results suggest that the majority of the sources are obscured. We use a number of diagnostic criteria in the optical, infrared and X-ray regime to verify these results. Interestingly, only 35% of the most luminous mid-IR selected AGNs have X-ray counterparts suggesting strong absorption. Our work emphasizes the importance of using SED decomposition techniques to select a population of type II AGNs, which may remain undetected by either X-ray or IR colour surveys.


title: The effectiveness of backward contact tracing in networks

abstract: Discovering and isolating infected individuals is a cornerstone of epidemic control. Because many infectious diseases spread through close contacts, contact tracing is a key tool for case discovery and control. However, although contact tracing has been performed widely, the mathematical understanding of contact tracing has not been fully established and it has not been clearly understood what determines the efficacy of contact tracing. Here, we reveal that, compared with "forward" tracing---tracing to whom disease spreads, "backward" tracing---tracing from whom disease spreads---is profoundly more effective. The effectiveness of backward tracing is due to simple but overlooked biases arising from the heterogeneity in contacts. Using simulations on both synthetic and high-resolution empirical contact datasets, we show that even at a small probability of detecting infected individuals, strategically executed contact tracing can prevent a significant fraction of further transmissions. We also show that---in terms of the number of prevented transmissions per isolation---case isolation combined with a small amount of contact tracing is more efficient than case isolation alone. By demonstrating that backward contact tracing is highly effective at discovering super-spreading events, we argue that the potential effectiveness of contact tracing has been underestimated. Therefore, there is a critical need for revisiting current contact tracing strategies so that they leverage all forms of biases. Our results also have important consequences for digital contact tracing because it will be crucial to incorporate the capability for backward and deep tracing while adhering to the privacy-preserving requirements of these new platforms.


title: Stably embedded submodels of Henselian valued fields

abstract: We show a transfer principle for the property that all types realised in a given elementary extension are definable. It can be written as follows: a Henselian valued fields is stably embedded in an elementary extension if and only if its value group is stably embedded in its corresponding extension, its residue field is stably embedded in its corresponding extension, and the extension of valued fields satisfies a certain algebraic condition. We show for instance that all types over the Hahn field $\mathbb{R}((\mathbb{Z}))$ are definable. Similarly, all types over the quotient field of the Witt ring $W(\mathbb{F}_p^{\text{alg}})$ are definable. This extends a work of Cubides and Delon and of Cubides and Ye.


title: Dimensional crossover and phase transitions in coupled chains: Density matrix renormalization group results

abstract: Quasi-one-dimensional (Q1D) systems, i.e., three- and two-dimensional (3D/2D) arrays composed of weakly coupled one-dimensional lattices of interacting quantum particles, exhibit rich and fascinating physics. They are studied across various areas of condensed matter and ultracold atomic lattice-gas physics, and are often marked by dimensional crossover as the coupling between one-dimensional systems is increased or temperature decreased, i.e., the Q1D system goes from appearing largely 1D to largely 3D. Phase transitions occurring along the crossover can strongly enhance this effect. Understanding these crossovers and associated phase transitions can be challenging due to the very different elementary excitations of 1D systems compared to higher-dimensional ones. In the present work, we combine numerical matrix product state (MPS) methods with mean-field (MF) theory to study paradigmatic cases of dimensional crossovers and the associated phase transitions in systems of both hard-core and soft-core lattice bosons, with relevance to both condensed matter physics and ultracold atomic gases. We show that the superfluid-to-insulator transition is a first order one, as opposed to the isotropic cases and calculate transition temperatures for the superfluid states, finding excellent agreement with analytical theory. At the same time, our MPS+MF approach keeps functioning well where the current analytical framework cannot be applied. We further confirm the qualitative and quantitative reliability of our approach by comparison to exact quantum Monte Carlo calculations for the full 3D arrays.


title: SLEDGE: A Simple Yet Effective Baseline for COVID-19 Scientific Knowledge Search

abstract: With worldwide concerns surrounding the Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2), there is a rapidly growing body of literature on the virus. Clinicians, researchers, and policy-makers need a way to effectively search these articles. In this work, we present a search system called SLEDGE, which utilizes SciBERT to effectively re-rank articles. We train the model on a general-domain answer ranking dataset, and transfer the relevance signals to SARS-CoV-2 for evaluation. We observe SLEDGE's effectiveness as a strong baseline on the TREC-COVID challenge (topping the learderboard with an nDCG@10 of 0.6844). Insights provided by a detailed analysis provide some potential future directions to explore, including the importance of filtering by date and the potential of neural methods that rely more heavily on count signals. We release the code to facilitate future work on this critical task at https://github.com/Georgetown-IR-Lab/covid-neural-ir


title: Inflation in Motion: Unitarity Constraints in Effective Field Theories with Broken Lorentz Symmetry

abstract: During inflation, there is a preferred reference frame in which the expansion of the background spacetime is spatially isotropic. In contrast to Minkowski spacetime, observables can depend on the velocity of the system with respect to this cosmic rest frame. We derive new constraints from radiative stability and unitarity on effective field theories with such spontaneously broken Lorentz symmetry. In addition to a maximum energy scale, there is now also a critical velocity at which the theory breaks down. The theory therefore has different resolving power in time and in space, and we show that these can only coincide if cubic Lorentz-violating interactions are absent. Applying these bounds to the Effective Field Theory of Inflation, we identify the region of parameter space in which inflation can be both single-field and weakly coupled on subhorizon scales. This can be implemented as a theoretical prior, and we illustrate this explicitly using Planck observational constraints on the primordial bispectrum.


title: CODA-19: Using a Non-Expert Crowd to Annotate Research Aspects on 10,000+ Abstracts in the COVID-19 Open Research Dataset

abstract: This paper introduces CODA-19, a human-annotated dataset that codes the Background, Purpose, Method, Finding/Contribution, and Other sections of 10,966 English abstracts in the COVID-19 Open Research Dataset. CODA-19 was created by 248 crowd workers from Amazon Mechanical Turk within 10 days, and achieved labeling quality comparable to that of experts. Each abstract was annotated by nine different workers, and the final labels were acquired by majority vote. The inter-annotator agreement (Cohen's kappa) between the crowd and the biomedical expert (0.741) is comparable to inter-expert agreement (0.788). CODA-19's labels have an accuracy of 82.2% when compared to the biomedical expert's labels, while the accuracy between experts was 85.0%. Reliable human annotations help scientists access and integrate the rapidly accelerating coronavirus literature, and also serve as the battery of AI/NLP research, but obtaining expert annotations can be slow. We demonstrated that a non-expert crowd can be rapidly employed at scale to join the fight against COVID-19.


title: Fast Dynamic Cuts, Distances and Effective Resistances via Vertex Sparsifiers

abstract: We present a general framework of designing efficient dynamic approximate algorithms for optimization on undirected graphs. In particular, we develop a technique that, given any problem that admits a certain notion of vertex sparsifiers, gives data structures that maintain approximate solutions in sub-linear update and query time. We illustrate the applicability of our paradigm to the following problems. (1) A fully-dynamic algorithm that approximates all-pair maximum-flows/minimum-cuts up to a nearly logarithmic factor in $\tilde{O}(n^{2/3})$ amortized time against an oblivious adversary, and $\tilde{O}(m^{3/4})$ time against an adaptive adversary. (2) An incremental data structure that maintains $O(1)$-approximate shortest path in $n^{o(1)}$ time per operation, as well as fully dynamic approximate all-pair shortest path and transshipment in $\tilde{O}(n^{2/3+o(1)})$ amortized time per operation. (3) A fully-dynamic algorithm that approximates all-pair effective resistance up to an $(1+\epsilon)$ factor in $\tilde{O}(n^{2/3+o(1)} \epsilon^{-O(1)})$ amortized update time per operation. The key tool behind result (1) is the dynamic maintenance of an algorithmic construction due to Madry [FOCS' 10], which partitions a graph into a collection of simpler graph structures (known as j-trees) and approximately captures the cut-flow and metric structure of the graph. The $O(1)$-approximation guarantee of (2) is by adapting the distance oracles by [Thorup-Zwick JACM `05]. Result (3) is obtained by invoking the random-walk based spectral vertex sparsifier by [Durfee et al. STOC `19] in a hierarchical manner, while carefully keeping track of the recourse among levels in the hierarchy.


title: The Expander Hierarchy and its Applications to Dynamic Graph Algorithms

abstract: We introduce a notion for hierarchical graph clustering which we call the expander hierarchy and show a fully dynamic algorithm for maintaining such a hierarchy on a graph with $n$ vertices undergoing edge insertions and deletions using $n^{o(1)}$ update time. An expander hierarchy is a tree representation of graphs that faithfully captures the cut-flow structure and consequently our dynamic algorithm almost immediately implies several results including: (1) The first fully dynamic algorithm with $n^{o(1)}$ worst-case update time that allows querying $n^{o(1)}$-approximate conductance, $s$-$t$ maximum flows, and $s$-$t$ minimum cuts for any given $(s,t)$ in $O(\log^{1/6} n)$ time. Our results are deterministic and extend to multi-commodity cuts and flows. The key idea behind these results is a fully dynamic algorithm for maintaining a tree flow sparsifier, a notion introduced by R\"acke [FOCS'02] for constructing competitive oblivious routing schemes. (2) A deterministic fully dynamic connectivity algorithm with $n^{o(1)}$ worst-case update time. This significantly simplifies the recent algorithm by Chuzhoy et al.~that uses the framework of Nanongkai et al. [FOCS'17]. (3) The first non-trivial deterministic fully dynamic treewidth decomposition algorithm on constant-degree graphs with $n^{o(1)}$ worst-case update time that maintains a treewidth decomposition of width $\text{tw}(G)\cdot n^{o(1)}$ where $\text{tw}(G)$ denotes the treewidth of the current graph. Our technique is based on a new stronger notion of the expander decomposition, called the boundary-linked expander decomposition. This decomposition is more robust against updates and better captures the clustering structure of graphs. Given that the expander decomposition has proved extremely useful in many fields, we expect that our new notion will find more future applications.


title: Nonlinear response properties in the inferior colliculus: Common sound-frequency dependence for linear and nonlinear responses

abstract: Neurons in the main center of convergence in the auditory midbrain, the central nucleus of the inferior colliculus (ICC) have been shown to display either linear significant receptive fields, or both, linear and nonlinear significant receptive fields. In this study, we used reverse correlation to probe linear and nonlinear response properties of single neurons in the cat ICC. The receptive fields display areas of stimulus parameters leading to enhanced or inhibited spiking activity, and thus allow investigating the interplay to process complex sounds. Spiking responses were obtained from neural recordings of anesthetized cats in response to dynamic moving ripple (DMR) stimuli. The DMR sound contains amplitude and frequency modulations and allows systematically mapping neural preferences. Correlations of the stimulus envelope that preferably excite neurons can be mapped with the spike-triggered covariance. The spike-triggered average and -covariance were computed for the envelope of the DMR, separately for each frequency carrier (spanning a range of 0-5.5 octaves). This enables studying processing of the sound envelope, and to investigate whether nonlinearities are more pronounced at the neurons' preferred frequencies rather than at other frequencies. We find that more than half of the neurons (n=120) display significant nonlinear response properties at least at one frequency carrier. Nonlinearities are dominant at the neuron's best frequency. The nonlinear preferences can have either the same or opposite temporal receptive field pattern (e.g. on-off) as the linear preferences. No relationship to other neural properties such as feature-selectivity, phase-locking, or the like has been found. Thus, these nonlinearities do not seem to be linked to a specific type of neuron but to be inherent to ICC neurons indicating a diverse range of filtering characteristics.


title: A comparative study of host galaxy properties between Fast Radio Bursts and stellar transients

abstract: Recent arcsecond localizations of Fast Radio Bursts and identifications of their host galaxies confirmed their extragalactic origin.While FRB 121102 resides in the bright region of a dwarf star forming galaxy, other FRBs reside in more massive galaxies and are related to older stellar populations. We compare the host galaxy properties of {nine} FRBs with those of several types of stellar transients: from young to old population, long duration gamma ray bursts (LGRBs), superluminous supernovae (SLSNe), Type Ib/Ic supernovae (SN Ibc), Type II supernovae (SN II), type Ia supernovae (SN Ia), and short duration gamma ray bursts (SGRBs). We find that as a whole sample, the stellar mass and star formation rate of the FRB host galaxies prefer a medium to old population, and are against a young population similar to LGRBs and SLSNe by a null probability 0.02. Individually, the FRB 121102 host is consistent with that of young population objects; the FRB 180924 environment is similar to that of SGRBs; and the FRB 190523 environment is similar to those of SN Ia. These results are consistent with the magnetar engine model for FRBs, if both magnetars produced from extreme explosions (GRBs/SLSNe) and from regular channels (e.g. those producing Galactic magnetars) can produce FRBs.


title: Community Detection Clustering via Gumbel Softmax

abstract: Recently, in many systems such as speech recognition and visual processing, deep learning has been widely implemented. In this research, we are exploring the possibility of using deep learning in community detection among the graph datasets. Graphs have gained growing traction in different fields, including social networks, information graphs, the recommender system, and also life sciences. In this paper, we propose a method of community detection clustering the nodes of various graph datasets. We cluster different category datasets that belong to Affiliation networks, Animal networks, Human contact networks, Human social networks, Miscellaneous networks. The deep learning role in modeling the interaction between nodes in a network allows a revolution in the field of science relevant to graph network analysis. In this paper, we extend the gumbel softmax approach to graph network clustering. The experimental findings on specific graph datasets reveal that the new approach outperforms traditional clustering significantly, which strongly shows the efficacy of deep learning in graph community detection clustering. We do a series of experiments on our graph clustering algorithm, using various datasets: Zachary karate club, Highland Tribe, Train bombing, American Revolution, Dolphins, Zebra, Windsurfers, Les Mis\'erables, Political books.


title: Context-Oriented Behavioral Programming

abstract: Modern systems require programmers to develop code that dynamically adapts to different contexts, leading to the evolution of new context-oriented programming languages. These languages introduce new software-engineering challenges, such as: how to maintain and keep the separation of concerns of the codebase? how to model the changing behaviors? how to verify the system behavior? and more. This paper introduces Context-Oriented Behavioral Programming(COBP) - a novel paradigm for developing context-aware systems, centered on natural and incremental specification of context-dependent behaviors. As the name suggests, we combine behavioral-programming(BP) - a scenario-based modeling paradigm - with context idioms that explicitly specify when scenarios are relevant and what information they need. The core idea is to connect the behavioral model with a data model that represents the context, allowing an intuitive connection between the models via update and select queries. Combining behavioral-programming with context-oriented programming brings the best of the two worlds, solving issues that arise when using each of the approaches in separation. We begin with providing abstract semantics for COBP, laying the foundations for applying reasoning algorithms to context-aware behavioral programs. We then exemplify the semantics with formal specifications of systems, including a variant of Conway's Game of Life. Finally, we present a JavaScript-based implementation of the paradigm and provide two case studies of real-life context-aware systems (one in robotics and another in IoT) that were developed using this tool. Throughout the examples and case studies, we provide design patterns and a methodology for coping with the above challenges.


title: ARES I: WASP-76 b, A Tale of Two HST Spectra

abstract: We analyse the transmission and emission spectra of the ultra-hot Jupiter WASP-76b, observed with the G141 grism of the Hubble Space Telescope's Wide Field Camera 3 (WFC3). We reduce and fit the raw data for each observation using the open-source software Iraclis before performing a fully Bayesian retrieval using the publicly available analysis suite TauRex 3. Previous studies of the WFC3 transmission spectra of WASP-76 b found hints of titanium oxide (TiO) and vanadium oxide (VO) or non-grey clouds. Accounting for a fainter stellar companion to WASP-76, we reanalyse this data and show that removing the effects of this background star changes the slope of the spectrum, resulting in these visible absorbers no longer being detected, eliminating the need for a non-grey cloud model to adequately fit the data but maintaining the strong water feature previously seen. However, our analysis of the emission spectrum suggests the presence of TiO and an atmospheric thermal inversion, along with a significant amount of water. Given the brightness of the host star and the size of the atmospheric features, WASP-76 b is an excellent target for further characterisation with HST, or with future facilities, to better understand the nature of its atmosphere, to confirm the presence of TiO and to search for other optical absorbers.


title: Jamming and force distribution in growing epithelial tissue

abstract: We investigate morphologies of proliferating cellular tissues using a newly developed numerical simulation model for mechanical cell division and migration in 2D. The model is applied to a bimodal mixture consisting of stiff cells with a low growth potential and soft cells with a high growth potential; cancer cells are typically considered to be softer than healthy cells. In an even mixture, the soft cells develop into a tissue matrix and the stiff cells into a dendrite-like network structure. When soft cells are placed inside a tissue consisting of stiff cells (to model cancer growth), the soft cells develop to a fast growing tumor-like structure that gradually evacuates the stiff cell matrix. The model also demonstrates 1) how soft cells orient themselves in the direction of the largest effective stiffness as predicted by the theory of Bischofs and Schwarz (Proc. Natl. Acad. Sci U.S.A., 100, 9274--9279 (2003) and 2) that the orientation and force generation continue a few cell rows behind the soft-stiff interface. With increasing inter-cell friction, tumor growth slows down and cell death occurs. The contact force distribution between cells is demonstrated to be highly sensitive to cell type mixtures and cell-cell interactions, which indicates that local mechanical forces can be useful as a regulator of tissue formation. The results shed new light on established experimental data.


title: A generalized Gauss curvature flow related to the Orlicz-Minkowski problem

abstract: In this paper a generalized Gauss curvature flow about a convex hypersurface in the Euclidean $n$-space is studied. This flow is closely related to the Orlicz-Minkowski problem, which involves Gauss curvature and a function of support function. Under some appropriate assumptions, we prove the long-time existence and convergence of this flow. As a byproduct, two existence results of solutions to the even Orlicz-Minkowski problem are obtained, one of which improves the known result.


title: Complex a priori bounds for multicritical circle maps with bounded type rotation number

abstract: In this paper we study homeomorphisms of the circle with several critical points and bounded type rotation number. We prove complex a priori bounds for these maps. As an application, we get that bi-cubic circle maps with same bounded type rotation number are $C^{1+\alpha}$ rigid.


title: Single-Copy Certification of Two-Qubit Gates without Entanglement

abstract: A quantum state transformation can be generally approximated by single- and two-qubit gates. This, however, does not hold with noisy intermediate-scale quantum technologies due to the errors appearing in the gate operations, where errors of two-qubit gates such as controlled-NOT and SWAP operations are dominated. In this work, we present a cost efficient single-copy certification for a realization of a two-qubit gate in the presence of depolarization noise, where it is aimed to identify if the realization is noise-free, or not. It is shown that entangled resources such as entangled states and a joint measurement are not necessary for the purpose, i.e., a noise-free two-qubit gate is not needed to certify an implementation of a two-qubit gate. A proof-of-principle demonstration is presented with photonic qubits.


title: A Theory of the Saving Rate of the Rich

abstract: Empirical evidence suggests that the rich have higher propensity to save than do the poor. While this observation may appear to contradict the homotheticity of preferences, we theoretically show that that is not the case. Specifically, we consider an income fluctuation problem with homothetic preferences and general shocks and prove that consumption functions are asymptotically linear, with an exact analytical characterization of asymptotic marginal propensities to consume (MPC). We provide necessary and sufficient conditions for the asymptotic MPCs to be zero. We calibrate a model with standard constant relative risk aversion utility and show that zero asymptotic MPCs are empirically plausible, implying that our mechanism has the potential to accommodate a large saving rate of the rich and high wealth inequality (small Pareto exponent) as observed in the data.


title: Bit-Interleaved Coded Multiple Beamforming with Perfect Coding in Massive MIMO Systems

abstract: This letter investigates bit-interleaved coded multiple beamforming (BICMB) with perfect coding in millimeter-wave (mm-wave) massive multiple-input multiple-output (MIMO) systems to achieve both maximum diversity gain and multiplexing gain. Using perfect coding with BICMB enables us to do this. We show that by using BICMB and perfect coding, the diversity gain becomes independent from the number of transmitted data streams and the number of antennas in each remote antenna unit (RAU) at the transmitter and the receiver. The assumption is that the perfect channel state information (CSI) is known at both the transmitter and the receiver and the number of antennas goes to infinity. This latter assumption can be relaxed by a large number of antennas in each RAU, similar to the case for all massive MIMO research. Simulation results show that when the perfect channel state information assumption is satisfied, the use of BICMB with perfect coding results in the diversity gain values predicted by the analysis.


title: Time-lapse Study of Neural Networks Using Phase Imaging with Computational Specificity (PICS)

abstract: In life sciences, fluorescent labeling techniques are used to study molecular structures and interactions of cells. However, this type of cell imaging has its own limitations, one of which is that the process of staining the cells could be toxic to the cells and possibly damage them. We are specifically interested in time-lapse imaging of live neurons to study their growth and proliferation. Neurodegenerative diseases are characterized by phenotypic differences in neuron growth and arborization. This thesis proposes a label-free digital staining method using the deep convolutional neural network to address the issues with the previous cell imaging method. Our results show that a deep neural network, when trained on phase images with correct fluorescent labels, can correctly learn the necessary morphological information to successfully predict MAP2 and Tau labels. This inference, in turn, allows us to classify axons from dendrites in live, unlabeled neurons.


title: Equilibrium mean-field-like statistical models with KPZ scaling

abstract: We have considered three different "one-body" statistical systems involving Brownian excursions, which possess for fluctuations Kardar-Parisi-Zhang scaling with the critical exponent $\nu=\frac{1}{3}$. In all models imposed external constraints push the underlying stochastic process to a large deviation regime. Specifically, we have considered fluctuations for: (i) Brownian excursions on non-uniform finite trees with linearly growing branching originating from the mean-field approximation of the Dumitriu-Edelman representation of matrix models, (ii) (1+1)D "magnetic" Dyck paths within the strip of finite width, (iii) inflated ideal polymer ring with fixed gyration radius. In the latter problem cutting off the long-ranged spatial fluctuations and leaving only the "typical" modes for stretched paths, we ensure the KPZ-like scaling for bond fluctuations. To the contrary, summing up all normal modes, we get the Gaussian behavior. In all considered models, KPZ fluctuations emerge in presence of two complementary conditions: (i) the trajectories are pushed to a large deviation region of a phase space, and (ii) the trajectories are leaning on an impenetrable boundary.


title: On the fourth order Cattaneo equation of heat conduction with memory

abstract: The well known heat equation with finite speed of propagation proposed by Cattaneo is obtained by a more general fourth order PDE when a certain (small) parameters is put equal to zero. It seems that this fourth order equation has been essentially overlooked in the literature, in particular when it is subject to non homogeneous boundary conditions. In this paper we examine its well posedness and the asymptotic behavior when certain coefficients tend to singular values.


title: Compositionality of the MSO+U Logic

abstract: We prove that the MSO+U logic is compositional in the following sense: whether an MSO+U formula holds in a tree T depends only on MSO+U-definable properties of the root of T and of subtrees of T starting directly below the root. Another kind of compositionality follows: every MSO+U formula whose all free variables range only over finite sets of nodes (in particular, whose all free variables are first-order) can be rewritten into an MSO formula having access to properties of subtrees definable by MSO+U sentences (without free variables).


title: Infinitely many hyperelliptic curves with exactly two rational points: Part II

abstract: In the previous paper, Hirakawa and the author determined the set of rational points of a certain infinite family of hyperelliptic curves $C^{(p;i,j)}$ parametrized by a prime number $p$ and integers $i$, $j$. In the proof, we used the standard $2$-descent argument and a Lutz-Nagell theorem that was proven by Grant. In this paper, we extend the above work. By using the descent theorem, the proof for $j=2$ is reduced to elliptic curves of rank $0$ that are independent of $p$. On the other hand, for odd $j$, we consider another hyperelliptic curve $C'^{(p;i,j)}$ whose Jacobian variety is isogenous to that of $C^{(p;i,j)}$, and prove that the Mordell-Weil rank of the Jacobian variety of $C'^{(p;i,j)}$ is $0$ by $2$-descent. Then, we determine the set of rational points of $C^{(p;i,j)}$ by using the Lutz-Nagell type theorem.


title: The universal DAHA of type $(C_1^\vee,C_1)$ and Leonard triples

abstract: Assume that $\mathbb F$ is an algebraically closed field and $q$ is a nonzero scalar in $\mathbb F$ that is not a root of unity. The universal Askey--Wilson algebra $\triangle_q$ is a unital associative $\mathbb F$-algebra generated by $A,B, C$ and the relations state that each of $$ A+\frac{q BC-q^{-1} CB}{q^2-q^{-2}}, \qquad B+\frac{q CA-q^{-1} AC}{q^2-q^{-2}}, \qquad C+\frac{q AB-q^{-1} BA}{q^2-q^{-2}} $$ is central in $\triangle_q$. The universal DAHA $\mathfrak H_q$ of type $(C_1^\vee,C_1)$ is a unital associative $\mathbb F$-algebra generated by $\{t_i^{\pm 1}\}_{i=0}^3$ and the relations state that \begin{gather*} t_it_i^{-1}=t_i^{-1} t_i=1 \quad \hbox{for all $i=0,1,2,3$}; \\ \hbox{$t_i+t_i^{-1}$ is central} \quad \hbox{for all $i=0,1,2,3$}; \\ t_0t_1t_2t_3=q^{-1}. \end{gather*} It was given an $\mathbb F$-algebra homomorphism $\triangle_q\to \mathfrak H_q$ that sends \begin{eqnarray*} A &\mapsto & t_1 t_0+(t_1 t_0)^{-1}, \\ B &\mapsto & t_3 t_0+(t_3 t_0)^{-1}, \\ C &\mapsto & t_2 t_0+(t_2 t_0)^{-1}. \end{eqnarray*} Therefore any $\mathfrak H_q$-module can be considered as a $\triangle_q$-module. Let $V$ denote a finite-dimensional irreducible $\mathfrak H_q$-module. In this paper we show that $A,B,C$ are diagonalizable on $V$ if and only if $A,B,C$ act as Leonard triples on all composition factors of the $\triangle_q$-module $V$.


title: SurvLIME-Inf: A simplified modification of SurvLIME for explanation of machine learning survival models

abstract: A new modification of the explanation method SurvLIME called SurvLIME-Inf for explaining machine learning survival models is proposed. The basic idea behind SurvLIME as well as SurvLIME-Inf is to apply the Cox proportional hazards model to approximate the black-box survival model at the local area around a test example. The Cox model is used due to the linear relationship of covariates. In contrast to SurvLIME, the proposed modification uses $L_{\infty }$-norm for defining distances between approximating and approximated cumulative hazard functions. This leads to a simple linear programming problem for determining important features and for explaining the black-box model prediction. Moreover, SurvLIME-Inf outperforms SurvLIME when the training set is very small. Numerical experiments with synthetic and real datasets demonstrate the SurvLIME-Inf efficiency.


title: Application-oriented mathematical algorithms for group testing

abstract: We have a large number of samples and we want to find the infected ones using as few number of tests as possible. We can use group testing which tells about a small group of people whether at least one of them is infected. Group testing is particularly efficient if the infection rate is low. The goal of this article is to summarize and extend the mathematical knowledge about the most efficient group testing algorithms, focusing on real-life applications instead of pure mathematical motivations and approaches.


title: Jointly Sparse Support Recovery via Deep Auto-encoder with Applications in MIMO-based Grant-Free Random Access for mMTC

abstract: In this paper, a data-driven approach is proposed to jointly design the common sensing (measurement) matrix and jointly support recovery method for complex signals, using a standard deep auto-encoder for real numbers. The auto-encoder in the proposed approach includes an encoder that mimics the noisy linear measurement process for jointly sparse signals with a common sensing matrix, and a decoder that approximately performs jointly sparse support recovery based on the empirical covariance matrix of noisy linear measurements. The proposed approach can effectively utilize the feature of common support and properties of sparsity patterns to achieve high recovery accuracy, and has significantly shorter computation time than existing methods. We also study an application example, i.e., device activity detection in Multiple-Input Multiple-Output (MIMO)-based grant-free random access for massive machine type communications (mMTC). The numerical results show that the proposed approach can provide pilot sequences and device activity detection with better detection accuracy and substantially shorter computation time than well-known recovery methods.


title: Mechanism Design and Blockchains

abstract: Game theory is often used as a tool to analyze decentralized systems and their properties, in particular, blockchains. In this note, we take the opposite view. We argue that blockchains can and should be used to implement economic mechanisms because they can help to overcome problems that occur if trust in the mechanism designer cannot be assumed. Mechanism design deals with the allocation of resources to agents, often by extracting private information from them. Some mechanisms are immune to early information disclosure, while others may heavily depend on it. Some mechanisms have to randomize to achieve fairness and efficiency. Both issues, information disclosure, and randomness require trust in the mechanism designer. If there is no trust, mechanisms can be manipulated. We claim that mechanisms that use randomness or sequential information disclosure are much harder, if not impossible, to audit. Therefore, centralized implementation is often not a good solution. We consider some of the most frequently used mechanisms in practice and identify circumstances under which manipulation is possible. We propose a decentralized implementation of such mechanisms, that can be, in practical terms, realized by blockchain technology. Moreover, we argue in which environments a decentralized implementation of a mechanism brings a significant advantage.


title: Recurrence Relations for Values of the Riemann Zeta Function in Odd Integers

abstract: It is commonly known that $\zeta(2k) = q_{k}\frac{\zeta(2k + 2)}{\pi^2}$ with known rational numbers $q_{k}$. In this work we construct recurrence relations of the form $\sum_{k = 1}^{\infty}r_{k}\frac{\zeta(2k + 1)}{\pi^{2k}} = 0$ and show that series representations for the coefficients $r_{k} \in \mathbb{R}$ can be computed explicitly.


title: Deep Constraint-based Propagation in Graph Neural Networks

abstract: The popularity of deep learning techniques renewed the interest in neural architectures able to process complex structures that can be represented using graphs, inspired by Graph Neural Networks (GNNs). We focus our attention on the originally proposed GNN model of Scarselli et al. 2009, which encodes the state of the nodes of the graph by means of an iterative diffusion procedure that, during the learning stage, must be computed at every epoch, until the fixed point of a learnable state transition function is reached, propagating the information among the neighbouring nodes. We propose a novel approach to learning in GNNs, based on constrained optimization in the Lagrangian framework. Learning both the transition function and the node states is the outcome of a joint process, in which the state convergence procedure is implicitly expressed by a constraint satisfaction mechanism, avoiding iterative epoch-wise procedures and the network unfolding. Our computational structure searches for saddle points of the Lagrangian in the adjoint space composed of weights, nodes state variables and Lagrange multipliers. This process is further enhanced by multiple layers of constraints that accelerate the diffusion process. An experimental analysis shows that the proposed approach compares favourably with popular models on several benchmarks.


title: Primes in arithmetic progressions and semidefinite programming

abstract: Assuming the generalized Riemann hypothesis, we give asymptotic bounds on the size of intervals that contain primes from a given arithmetic progression using the approach developed by Carneiro, Milinovich and Soundararajan [Comment. Math. Helv. 94, no. 3 (2019)]. For this we extend the Guinand-Weil explicit formula over all Dirichlet characters modulo $q \geq 3$, and we reduce the associated extremal problems to convex optimization problems that can be solved numerically via semidefinite programming.


title: Lessons from the Ramond sector

abstract: We revisit the consistency of torus partition functions in (1+1)$d$ fermionic conformal field theories, combining traditional ingredients of modular invariance/covariance with a modernized understanding of bosonization/fermionization dualities. Various lessons can be learned by simply examining the oft-ignored Ramond sector. For several extremal/kinky modular functions in the bootstrap literature, we can either rule out or identify the underlying theory. We also revisit the ${\cal N} = 1$ Maloney-Witten partition function by calculating the spectrum in the Ramond sector, and further extending it to include the modular sum of seed Ramond characters. Finally, we perform the full ${\cal N} = 1$ RNS modular bootstrap and obtain new universal results on the existence of relevant deformations preserving different amounts of supersymmetry.


title: Bayesian constraints on the astrophysical neutrino source population from IceCube data

abstract: We present constraints on an astrophysical population of neutrino sources imposed by recent data from the IceCube neutrino observatory. By using the IceCube point source search method to model the detection of sources, our detection criterion is more sensitive than using the observation of high-energy neutrino multiplets for source identification. We frame the problem as a Bayesian hierarchical model to connect the high-level population parameters to the IceCube data, allowing us to consistently account for all relevant sources of uncertainty in our model assumptions. Our results show that sources with a local density of $n_0 \gtrsim 10^{-7}$ $\rm{Mpc}^{-3}$ and luminosity $L \lesssim 10^{43}$ erg/s are the most likely candidates, but that populations of rare sources with $n_0 \simeq 10^{-9}$ $\rm{Mpc}^{-3}$ and $L \simeq 10^{45}$ erg/s can still be consistent with the IceCube observations. We demonstrate that these conclusions are strongly dependent on the source evolution considered, for which we consider a wide range of models. In doing so, we present realistic, model-independent constraints on the population parameters that reflect our current state of knowledge from astrophysical neutrino observations. We also use our framework to investigate constraints in the case of possible source detections and future instrument upgrades. Our approach is flexible and can be used to model specific source cases and extended to include multi-messenger information.


title: Unquantized thermal Hall effect in quantum spin liquids with spinon Fermi surfaces

abstract: Recent theoretical studies have found quantum spin liquid states with spinon Fermi surfaces upon the application of a magnetic field on a gapped state with topological order. We investigate the thermal Hall conductivity across this transition, describing how the quantized thermal Hall conductivity of the gapped state changes to an unquantized thermal Hall conductivity in the gapless spinon Fermi surface state. We consider two cases, both of potential experimental interest: the state with non-Abelian Ising topological order on the honeycomb lattice, and the state with Abelian chiral spin liquid topological order on the triangular lattice.


title: Dark Higgs Dark Matter

abstract: A new $U(1)$ dark gauge group coupled to the Standard Model (SM) via the kinetic mixing portal provides a natural dark matter candidate in the form of the Higgs field, $h_d$, responsible for generating the mass of the dark photon, $\gamma_d$. We show that the condition $m_{h_d}\leq m_{\gamma_d}$, together with smallness of the kinetic mixing parameter, $\epsilon$, and/or dark gauge coupling, $g_d$, leads the dark Higgs to be sufficiently metastable to constitute dark matter. We analyze the Universe's thermal history and show that both freeze-in, ${\rm SM}\to \{\gamma_d, h_d\}$, and freeze-out, $ \{\gamma_d, h_d\} \to {\rm SM}$, processes can lead to viable dark Higgs dark matter with a sub-GeV mass and a kinetic mixing parameter in the range $10^{-13}\lesssim\epsilon\lesssim10^{-6}$. Observable signals in astrophysics and cosmology include modifications to primordial elemental abundances, altered energetics of supernovae explosions, dark Higgs decays in the late Universe, and dark matter self-interactions.


title: Polarized jet fragmentation functions

abstract: We develop the theoretical framework needed to study the distribution of hadrons with general polarization inside jets, with and without transverse momentum measured with respect to the standard jet axis. The key development in this paper, referred to as "polarized jet fragmentation functions", opens up new opportunities to study both collinear and transverse momentum dependent (TMD) fragmentation functions. As two examples of the developed framework, we study longitudinally polarized collinear $\Lambda$ and transversely polarized TMD $\Lambda$ production inside jets in both $pp$ and $ep$ collisions. We find that both observables have high potential in constraining spin-dependent fragmentation functions with sizeable asymmetries predicted, in particular, at the future Electron-Ion Collider.


title: Creating and manipulating a Laughlin-type $\nu=1/3$ fractional quantum Hall state on a quantum computer with linear depth circuits

abstract: Here we present an efficient quantum algorithm to generate an equivalent many-body state to Laughlin's $\nu=1/3$ fractional quantum Hall state on a digitized quantum computer. Our algorithm only uses quantum gates acting on neighboring qubits in a quasi-one-dimensional setting, and its circuit depth is linear in the number of qubits, i.e., the number of Landau orbitals in the second quantized picture. We identify correlation functions that serve as signatures of the Laughlin state and discuss how to obtain them on a quantum computer. We also discuss a generalization of the algorithm for creating quasiparticles in the Laughlin state. This paves the way for several important studies, including quantum simulation of nonequilibrium dynamics and braiding of quasiparticles in quantum Hall states.


title: Signatures of a topological Weyl loop in Co$_3$Sn$_2$S$_2$

abstract: The search for novel topological phases of matter in quantum magnets has emerged as a frontier of condensed matter physics. Here we use state-of-the-art angle-resolved photoemission spectroscopy (ARPES) to investigate single crystals of Co$_3$Sn$_2$S$_2$ in its ferromagnetic phase. We report for the first time signatures of a topological Weyl loop. From fundamental symmetry considerations, this magnetic Weyl loop is expected to be gapless if spin-orbit coupling (SOC) is strictly zero but gapped, with possible Weyl points, under finite SOC. We point out that high-resolution ARPES results to date cannot unambiguously resolve the SOC gap anywhere along the Weyl loop, leaving open the possibility that Co$_3$Sn$_2$S$_2$ hosts zero Weyl points or some non-zero number of Weyl points. On the surface of our samples, we further observe a possible Fermi arc, but we are unable to clearly verify its topological nature using the established counting criteria. As a result, we argue that from the point of view of photoemission spectroscopy the presence of Weyl points and Fermi arcs in Co$_3$Sn$_2$S$_2$ remains ambiguous. Our results have implications for ongoing investigations of Co$_3$Sn$_2$S$_2$ and other topological magnets.


title: The kinematics of globular cluster populations in the E-MOSAICS simulations and their implications for the assembly history of the Milky Way

abstract: We present a detailed comparison of the Milky Way (MW) globular cluster (GC) kinematics with the 25 Milky Way-mass cosmological simulations from the E-MOSAICS project. While the MW falls within the kinematic distribution of GCs spanned by the simulations, the relative kinematics of its metal-rich ($[\rm{Fe/H}]>-1.2$) versus metal-poor ($[\rm{Fe/H}]<-1.2$), and inner ($r<8\rm{kpc}$) versus outer ($r>8\rm{kpc}$) populations are atypical for its mass. To understand the origins of these features, we perform a comprehensive statistical analysis of the simulations, and find 18 correlations describing the assembly of $L^*$ galaxies and their dark matter haloes based on their GC population kinematics. The correlations arise because the orbital distributions of accreted and in-situ GCs depend on the masses and accretion redshifts of accreted satellites, driven by the combined effects of dynamical fraction, tidal stripping, and dynamical heating. Because the kinematics of in-situ/accreted GCs are broadly traced by the metal-rich/metal-poor and inner/outer populations, the observed GC kinematics are a sensitive probe of galaxy assembly. We predict that relative to the population of $L^*$ galaxies, the MW assembled its dark matter and stellar mass rapidly through a combination of in-situ star formation, more than a dozen low-mass mergers, and $1.4\pm1.2$ early ($z=3.1\pm1.3$) major merger. The rapid assembly period ended early, limiting the fraction of accreted stars. We conclude by providing detailed quantitative predictions for the assembly history of the MW.


title: The [OIII]$+$H$\beta$ Equivalent Width Distribution at z$\simeq$7: Implications for the Contribution of Galaxies to Reionization

abstract: We quantify the distribution of [OIII]+H$\beta$ line strengths at z$\simeq$7 using a sample of 20 bright (M$_{\mathrm{UV}}$ $\lesssim$ $-$21) galaxies. We select these systems over wide-area fields (2.3 deg$^2$ total) using a new colour-selection which precisely selects galaxies at z$\simeq$6.63$-$6.83, a redshift range where blue Spitzer/IRAC [3.6]$-$[4.5] colours unambiguously indicate strong [OIII]$+$H$\beta$ emission. These 20 galaxies suggest a log-normal [OIII]$+$H$\beta$ EW distribution with median EW = 759$^{+112}_{-113}$ $\mathrm{\mathring{A}}$ and standard deviation = 0.26$^{+0.06}_{-0.05}$ dex. We find no evidence for strong variation in this EW distribution with UV luminosity. The typical [OIII]+H$\beta$ EW at z$\simeq$7 implied by our sample is considerably larger than that in massive star forming galaxies at z$\simeq$2, consistent with a shift toward larger average sSFR (4.4 Gyr$^{-1}$) and lower metallicities (0.16 Z$_\odot$). We also find evidence for the emergence of a population with yet more extreme nebular emission ([OIII]+H$\beta$ EW$>$1200 $\mathrm{\mathring{A}}$) that is rarely seen at lower redshifts. These objects have extremely large sSFR ($>$30 Gyr$^{-1}$), as would be expected for systems undergoing a burst or upturn in star formation. While this may be a short-lived phase, our results suggest that 20% of the z$\simeq$7 population has such extreme nebular emission, implying that galaxies likely undergo intense star formation episodes regularly at z$>$6. We argue that this population may be among the most effective ionizing agents in the reionization era, both in terms of photon production efficiency and escape fraction. We furthermore suggest that galaxies passing through this large sSFR phase are likely to be very efficient in forming bound star clusters.


title: Quantum advantage in simulating stochastic processes

abstract: We investigate the problem of simulating classical stochastic processes through quantum dynamics, and present three scenarios where memory or time quantum advantages arise. First, by introducing and analysing a quantum version of the embeddability problem for stochastic matrices, we show that quantum memoryless dynamics can simulate classical processes that necessarily require memory. Second, by extending the notion of space-time cost of a stochastic process $P$ to the quantum domain, we prove an advantage of the quantum cost of simulating $P$ over the classical cost. Third, we demonstrate that the set of classical states accessible via Markovian master equations with quantum controls is larger than the set of those accessible with classical controls, leading, e.g., to a potential advantage in cooling protocols.


title: Non-equilibrium readiness and accuracy of Gaussian Quantum Thermometers

abstract: The dimensionality of a thermometer is key in the design of quantum thermometry schemes. In general, the phenomenology that is typical of finite-dimensional quantum thermometry does not apply to infinite dimensional ones. We analyse the dynamical and metrological features of non-equilibrium Gaussian Quantum Thermometers: on one hand, we highlight how quantum entanglement can enhance the readiness of composite Gaussian thermometers; on the other hand, we show that non-equilibrium conditions do not guarantee the best sensitivities in temperature estimation, thus suggesting the reassessment of the working principles of quantum thermometry.


title: Merger of Dark Matter Axion Clumps and Resonant Photon Emission

abstract: A portion of light scalar dark matter, especially axions, may organize into gravitationally bound clumps (stars) and be present in large number in the galaxy today. It is therefore of utmost interest to determine if there are novel observational signatures of this scenario. Work has shown that for moderately large axion-photon couplings, such clumps can undergo parametric resonance into photons, for clumps above a critical mass $M^{\star}_c$ determined precisely by some of us in Ref. [1]. In order to obtain a clump above the critical mass in the galaxy today would require mergers. In this work we perform full 3-dimensional simulations of pairs of axion clumps and determine the conditions under which mergers take place through the emission of scalar waves, including analyzing head-on and non-head-on collisions, phase dependence, and relative velocities. Consistent with other work in the literature, we find that the final mass from the merger $M^{\star}_{\text{final}}\approx 0.7(M^{\star}_1+M^{\star}_2)$ is larger than each of the original clump masses (for $M^{\star}_1\sim M^{\star}_2$). Hence, it is possible for sub-critical mass clumps to merge and become super-critical and therefore undergo parametric resonance into photons. We find that mergers are expected to be kinematically allowed in the galaxy today for high Peccei-Quinn scales, which is strongly suggested by unification ideas, although the collision rate is small. While mergers can happen for axions with lower Peccei-Quinn scales due to statistical fluctuations in relative velocities, as they have a high collision rate. We estimate the collision and merger rates within the Milky Way galaxy today. We find that a merger leads to a flux of energy on earth that can be appreciable and we mention observational search strategies.

categories: hep-ph astro-ph.CO astro-ph.HE gr-qc hep-th


title: Transport across twist angle domains in moir\'e graphene

abstract: Many of the experiments in twisted bilayer graphene (TBG) differ from each other in terms of the details of their phase diagrams. Few controllable aspects aside, this discrepancy is largely believed to be arising from the presence of a varying degree of twist angle inhomogeneity across different samples. Real space maps indeed reveal TBG devices splitting into several large domains of different twist angles. Motivated by these observations, we study the quantum mechanical tunneling across a domain wall (DW) that separates two such regions. We show that the tunneling of the moir\'e particles can be understood by the formation of an effective step potential at the DW. The height of this step potential is simply a measure of the difference in twist angles. These computations lead us to identify the global transport signatures for detecting and quantifying the local twist angle variations. In particular, Using Landauer-B\"uttiker formalism we compute single-channel conductance ($dI/dV$) and Fano factor for shot noise (ratio of noise power and mean current). A zero-bias, sub-meV transport gap is observed in the conductance which scales with the height of the step potential. One of the key findings of our work is that transport in presence of twist angle inhomogeneity is "noisy", though sub-Poissonian. In particular, the differential Fano factor peaks near the van Hove energies corresponding to the domains in the sample. The location and the strength of the peak is simply a measure of the degree of twist angle inhomogeneity.


title: The e-MERLIN Galaxy Evolution Survey (e-MERGE): Overview and Survey Description

abstract: We present an overview and description of the eMERLIN Galaxy Evolution survey (eMERGE) Data Release 1 (DR1), a large program of high-resolution 1.5 GHz radio observations of the GOODS-N field comprising $\sim140$ hours of observations with eMERLIN and $\sim40$ hours with the Very Large Array (VLA). We combine the long baselines of eMERLIN (providing high angular resolution) with the relatively closely-packed antennas of the VLA (providing excellent surface brightness sensitivity) to produce a deep 1.5 GHz radio survey with the sensitivity ($\sim 1.5\mu$Jy beam$^{-1}$), angular resolution ($0.2"$--$0.7"$) and field-of-view ($\sim15' \times 15'$) to detect and spatially resolve star-forming galaxies and AGN at $z\gtrsim 1$. The goal of eMERGE is to provide new constraints on the deep, sub-arcsecond radio sky which will be surveyed by SKA1-mid. In this initial publication, we discuss our data analysis techniques, including steps taken to model in-beam source variability over a $\sim20$ year baseline and the development of new point spread function/primary beam models to seamlessly merge eMERLIN and VLA data in the $uv$ plane. We present early science results, including measurements of the luminosities and/or linear sizes of $\sim500$ galaxes selected at 1.5 GHz. In combination with deep Hubble Space Telescope observations, we measure a mean radio-to-optical size ratio of $r_{\rm eMERGE}/r_{\rm HST}\sim1.02\pm0.03$, suggesting that in most high-redshift galaxies, the $\sim$GHz continuum emission traces the stellar light seen in optical imaging. This is the first in a series of papers which will explore the $\sim$kpc-scale radio properties of star-forming galaxies and AGN in the GOODS-N field observed by eMERGE DR1.


title: The Cosmic Ultraviolet Baryon Survey (CUBS) I. Overview and the diverse environments of Lyman limit systems at z<1

abstract: We present initial results from the Cosmic Ultraviolet Baryon Survey (CUBS). CUBS is designed to map diffuse baryonic structures at redshift z<~1 using absorption-line spectroscopy of 15 UV-bright QSOs with matching deep galaxy survey data. CUBS QSOs are selected based on their NUV brightness to avoid biases against the presence of intervening Lyman Limit Systems (LLSs) at zabs<1. We report five new LLSs of log N(HI)/cm^-2 >~ 17.2 over a total redshift survey pathlength of dz=9.3, and a number density of n(z)=0.43 (-0.18, +0.26). Considering all absorbers with log N(HI)/cm^-2 > 16.5 leads to n(z)=1.08 (-0.25, +0.31) at z<1. All LLSs exhibit a multi-component structure and associated metal transitions from multiple ionization states such as CII, CIII, MgII, SiII, SiIII, and OVI absorption. Differential chemical enrichment levels as well as ionization states are directly observed across individual components in three LLSs. We present deep galaxy survey data obtained using the VLT-MUSE integral field spectrograph and the Magellan Telescopes, reaching sensitivities necessary for detecting galaxies fainter than 0.1L* at d<~300 physical kpc (pkpc) in all five fields. A diverse range of galaxy properties is seen around these LLSs, from a low-mass dwarf galaxy pair, a co-rotating gaseous halo/disk, a star-forming galaxy, a massive quiescent galaxy, to a galaxy group. The closest galaxies have projected distances ranging from d=15 to 72 pkpc and intrinsic luminosities from ~0.01L* to ~3L*. Our study shows that LLSs originate in a variety of galaxy environments and trace gaseous structures with a broad range of metallicities.


title: Global torques and stochasticity as the drivers of massive black hole pairing in the young Universe

abstract: The forthcoming Laser Interferometer Space Antenna (LISA) will probe the population of coalescing massive black hole (MBH) binaries up to the onset of structure formation. Here we simulate the galactic-scale pairing of $\sim10^6 M_\odot$ MBHs in a typical, non-clumpy main-sequence galaxy embedded in a cosmological environment at $z = 7-6$. In order to increase our statistical sample, we adopt a strategy that allows us to follow the evolution of six secondary MBHs concomitantly. We find that the magnitude of the dynamical-friction induced torques is significantly smaller than that of the large-scale, stochastic gravitational torques arising from the perturbed and morphologically evolving galactic disc, suggesting that the standard dynamical friction treatment is inadequate for realistic galaxies at high redshift. The dynamical evolution of MBHs is very stochastic, and a variation in the initial orbital phase can lead to a drastically different time-scale for the inspiral. Most remarkably, the development of a galactic bar in the host system either significantly accelerates the inspiral by dragging a secondary MBH into the centre, or ultimately hinders the orbital decay by scattering the MBH in the galaxy outskirts. The latter occurs more rarely, suggesting that galactic bars overall promote MBH inspiral and binary coalescence. The orbital decay time can be an order of magnitude shorter than what would be predicted relying on dynamical friction alone. The stochasticity, and the important role of global torques, have crucial implications for the rates of MBH coalescences in the early Universe: both have to be accounted for when making predictions for the upcoming LISA observatory.


title: Powering Central Compact Objects with a Tangled Crustal Magnetic Field

abstract: Central Compact Objects (CCOs) are X-ray sources with luminosity ranging between $10^{32}$-$10^{34}$ erg~s$^{-1}$, located at the centres of supernova remnants. Some of them have been confirmed to be neutron stars. Timing observations have allowed the estimation of their dipole magnetic field, placing them in the range $\sim10^{10}$-$10^{11}$ G. The decay of their weak dipole fields, mediated by the Hall effect and Ohmic dissipation, cannot provide sufficient thermal energy to power their X-ray luminosity, as opposed to magnetars whose X-ray luminosities are comparable. Motivated by the question of producing high X-ray power through magnetic field decay while maintaining a weak dipole field, we explore the evolution of a crustal magnetic field that does not consist of an ordered axisymmetric structure, but rather comprises a tangled configuration. This can be the outcome of a non-self-excited dynamo, buried inside the crust by fallback material following the supernova explosion. We find that such initial conditions lead to the emergence of the magnetic field from the surface of the star and the formation of a dipolar magnetic field component. An internal tangled magnetic field of the order of $10^{14}$ G can provide sufficient Ohmic heating to the crust and power CCOs, while the dipole field it forms is approximately $10^{10}$ G, as observed in CCOs.


title: Detecting axisymmetric magnetic fields using gravity modes in intermediate-mass stars

abstract: Context: Angular momentum (AM) transport models of stellar interiors require improvements to explain the strong extraction of AM from stellar cores that is observed with asteroseismology. One of the often invoked mediators of AM transport are internal magnetic fields, even though their properties, observational signatures and influence on stellar evolution are largely unknown. Aims: We study how a fossil, axisymmetric internal magnetic field affects period spacing patterns of dipolar gravity mode oscillations in main-sequence stars with masses of 1.3, 2.0 and 3.0 M$_\odot$. We assess the influence of fundamental stellar parameters on the magnitude of pulsation mode frequency shifts. Methods: We compute dipolar gravity mode frequency shifts due to a fossil, axisymmetric poloidal-toroidal internal magnetic field for a grid of stellar evolution models, varying stellar fundamental parameters. Rigid rotation is taken into account using the traditional approximation of rotation and the influence of the magnetic field is computed using a perturbative approach. Results: We find magnetic signatures for dipolar gravity mode oscillations in terminal-age main-sequence stars that are measurable for a near-core field strength larger than $10^{5}$ G. The predicted signatures differ appreciably from those due to rotation. Conclusions: Our formalism demonstrates the potential for the future detection and characterization of strong fossil, axisymmetric internal magnetic fields in gravity-mode pulsators near the end of core-hydrogen burning from Kepler photometry, if such fields exist.


title: A Blast from the Pasture: Studying the environment of AT 2018cow with MUSE

abstract: AT 2018cow was the nearest and best studied example of a new breed of extra-galactic, luminous and rapidly-evolving transient. Both the progenitor systems and explosion mechanisms of these rapid transients remain a mystery - the energetics, spectral signatures, and timescales make them challenging to interpret in established classes of supernovae and tidal disruption events. The rich, multi-wavelength data-set of AT 2018cow has still left several interpretations viable to explain the nature of this event. In this paper we analyse integral-field spectroscopic data of the host galaxy, CGCG 137-068, to compare environmental constraints with leading progenitor models. We find the explosion site of AT 2018cow to be very typical of core-collapse supernovae (known to form from stars with MZAMS ~8-25M), and infer a young stellar population age at the explosion site of few times 10Myr, at slightly sub-solar metallicity. When comparing to expectations for exotic intermediate-mass black hole (IMBH) tidal disruption events,we find no evidence for a potential host system of the IMBH. In particular, there are no abrupt changes in metallicity or kinematics in the vicinity of the explosion site, arguing against the presence of a distinct host system. The proximity of AT 2018cow to strong star-formation in the host galaxy makes us favour a massive stellar progenitor for this event.