title: Exchange interactions in two-state systems: rare earth pyrochlores

abstract: The general form of the nearest neighbour exchange interaction for rare earth pyrochlores is derived based on symmetry. Generally, the rare earth angular momentum degeneracy is lifted by the crystal electric field (CEF) into singlets and doublets. When the CEF ground state is a doublet that is well-separated from the first excited state the CEF ground state doublet can be treated as a pseudo-spin of some kind. The general form of the nearest neighbour exchange interaction for pseudo-spins on the pyrochlore lattice is derived for three different types of pseudo-spins. The methodology presented in this paper can be applied to other two-state spin systems with a high space group symmetry.


title: Systolic inequalities for K3 surfaces via stability conditions

abstract: We introduce the notions of categorical systoles and categorical volumes of Bridgeland stability conditions on triangulated categories. We prove that for any projective K3 surface, there exists a constant C depending only on the rank and discriminant of its Picard group, such that $$\mathrm{sys}(\sigma)^2\leq C\cdot\mathrm{vol}(\sigma)$$ holds for any stability condition on the derived category of coherent sheaves on the K3 surface. This is an algebro-geometric generalization of a classical systolic inequality on two-tori. We also discuss applications of this inequality in symplectic geometry.


title: Heterogeneous Turbulent Jets of the Mutually Immiscible Liquids: the Mixing and Heat Transfer Processes

abstract: Many natural and technical processes deal with the turbulent mixing and heat transfer in the jet of mutually immiscible liquids, which represent an important class of the modern multiphase systems dynamics. The differential equations for axisymmetrical two-dimensional stationary flow and the integral correlations in a cylindrical coordinate system are considered for the free and confined jets. The parameters of the turbulent mixing in a free jet and in two-phase flow in a chamber are modelled and analyzed. One example is also done for high-temperature flow of the liquid metal melt cooled by water in a simulation of the hypothetic severe accidents at the nuclear power plant is presented and compared to the experimental data. The algorithm and the results obtained may be of interest for some research and industrial tasks, where calculation of the parameters of multiphase turbulent mixing and heat transfer are important.


title: Strict monotonicity of percolation thresholds under covering maps

abstract: We answer a question of Benjamini and Schramm by proving that under reasonable conditions, quotienting a graph strictly increases the value of its percolation critical parameter $p_c$. More precisely, let $\mathcal{G}=(V,E)$ be a quasi-transitive graph with $p_c(\mathcal{G})<1$, and let $G$ be a nontrivial group that acts freely on $V$ by graph automorphisms. Assume that $\mathcal{H}:=\mathcal{G}/G$ is quasi-transitive. Then one has $p_c(\mathcal{G})<p_c(\mathcal{H})$. We provide results beyond this setting: we treat the case of general covering maps and provide a similar result for the uniqueness parameter $p_u$, under an additional assumption of boundedness of the fibres. The proof makes use of a coupling built by lifting the exploration of the cluster, and an exploratory counterpart of Aizenman-Grimmett's essential enhancements.


title: New formulas for the Laplacian of distance functions and applications

abstract: The goal of the paper is to prove an exact representation formula for the Laplacian of the distance (and more generally for an arbitrary 1-Lipschitz function) in the framework of metric measure spaces satisfying Ricci curvature lower bounds in a synthetic sense (more precisely in essentially non-branching MCP(K,N)-spaces). Such a representation formula makes apparent the classical upper bounds and also some new lower bounds, together with a precise description of the singular part. The exact representation formula for the Laplacian of 1-Lipschitz functions (in particular for distance functions) holds also (and seems new) in a general complete Riemannian manifold. We apply these results to prove the equivalence of CD(K,N) and a dimensional Bochner inequality on signed distance functions. Moreover we obtain a measure-theoretic Splitting Theorem for infinitesimally Hilbertian essentially non-branching spaces verifying MCP(0,N).


title: Optimisation-based representations for branching processes

abstract: It is shown that a certain functional of a branching process has representations in terms of both a maximisation problem and a minimisation problem. A consequence of these representation is that upper and lower bounds on the functional can be found easily, yielding a non-asymptotic Trotter product formula. As an application, the speed of the right-most particle of a branching Levy process is calculated.


title: Flow From Motion: A Deep Learning Approach

abstract: Wearable devices have the potential to enhance sports performance, yet they are not fulfilling this promise. Our previous studies with 6 professional tennis coaches and 20 players indicate that this could be due the lack of psychological or mental state feedback, which the coaches claim to provide. Towards this end, we propose to detect the flow state, mental state of optimal performance, using wearables data to be later used in training. We performed a study with a professional tennis coach and two players. The coach provided labels about the players' flow state while each player had a wearable device on their racket holding wrist. We trained multiple models using the wearables data and the coach labels. Our deep neural network models achieved around 98% testing accuracy for a variety of conditions. This suggests that the flow state or what coaches recognize as flow, can be detected using wearables data in tennis which is a novel result. The implication for the HCI community is that having access to such information would allow for design of novel hardware and interaction paradigms that would be helpful in professional athlete training.


title: A simple, space constrained NIRIM type reactor for chemical vapour deposition of diamond

abstract: In this paper the design of a simple, space constrained chemical vapour deposition reactor for diamond growth is detailed. Based on the design by NIRIM, the reactor is composed of a quartz discharge tube placed within a 2.45 GHz waveguide to create the conditions required for metastable growth of diamond. Utilising largely off-the-shelf components and a modular design, the reactor allows for easy modification, repair, and cleaning between growth runs. The elements of the reactor design are laid out with the CAD files, parts list, and control files made easily available to enable replication. Finally, the quality of nanocrystalline diamond films produced are studied with SEM and Raman spectroscopy, with the observation of clear faceting and a large diamond fraction suggesting the design offers deposition of diamond with minimal complexity.


title: Design optimisation and post-trial analysis in group sequential stepped-wedge cluster randomised trials

abstract: Recently, methodology was presented to facilitate the incorporation of interim analyses in stepped-wedge (SW) cluster randomised trials (CRTs). Here, we extend this previous discussion. We detail how the stopping boundaries, allocation sequences, and per-cluster per-period sample size of a group sequential SW-CRT can be optimised. We then describe methods by which point estimates, p-values, and confidence intervals, which account for the sequential nature of the design, can be calculated. We demonstrate that optimal sequential designs can reduce the expected required number of measurements under the null hypothesis, compared to the classical design, by up to 30%, with no cost to the maximal possible required number of measurements. Furthermore, the adjusted analysis procedure almost universally reduces the average bias in the point estimate, and consistently provides a confidence interval with coverage close to the nominal level. In contrast, the coverage of a naive 95% confidence interval is observed to range between 92 and 98%. Methodology is now readily available for the efficient design and analysis of group sequential SW-CRTs. In scenarios in which there are substantial ethical or financial reasons to terminate a SW-CRT as soon as possible, trialists should strongly consider a group sequential approach.


title: Large nematic susceptibility in the double-Q $C_{4}$ magnetic phase of Ba$_{1-x}$Na$_{x}$Fe$_{2}$As$_{2}$

abstract: The nematic susceptibility of Ba$_{1-x}$Na$_{x}$Fe$_{2}$As$_{2}$ single crystals is studied by measuring the Young's modulus using a three-point-bending setup in a capacitance dilatometer over a wide doping range. Particular emphasis is placed on the behavior within the double-Q antiferromagnetic $C_{4}$ re-entrant phase. Here, we surprisingly still observe a sizable nematic susceptibility in spite of the well-developed magnetic order, consistent with recent theoretical calculations. Outside the $C_{4}$ re-entrant region, the behavior is very similar to that of K-doped BaFe$_{2}$As$_{2}$. A significant coupling of superconductivity to the shear modulus is observed and is discussed in terms of strong competition between the superconducting and magnetic phases.


title: Efficient Interactive Annotation of Segmentation Datasets with Polygon-RNN++

abstract: Manually labeling datasets with object masks is extremely time consuming. In this work, we follow the idea of Polygon-RNN to produce polygonal annotations of objects interactively using humans-in-the-loop. We introduce several important improvements to the model: 1) we design a new CNN encoder architecture, 2) show how to effectively train the model with Reinforcement Learning, and 3) significantly increase the output resolution using a Graph Neural Network, allowing the model to accurately annotate high-resolution objects in images. Extensive evaluation on the Cityscapes dataset shows that our model, which we refer to as Polygon-RNN++, significantly outperforms the original model in both automatic (10% absolute and 16% relative improvement in mean IoU) and interactive modes (requiring 50% fewer clicks by annotators). We further analyze the cross-domain scenario in which our model is trained on one dataset, and used out of the box on datasets from varying domains. The results show that Polygon-RNN++ exhibits powerful generalization capabilities, achieving significant improvements over existing pixel-wise methods. Using simple online fine-tuning we further achieve a high reduction in annotation time for new datasets, moving a step closer towards an interactive annotation tool to be used in practice.


title: The continuum self-similar tree

abstract: We introduce the continuum self-similar tree (CSST) and characterize it topologically. We apply this to answer a question of Curien about the topology of the continuum random tree (CRT). We also give a topological characterization of other trees with branch points of finite or infinite valences.


title: A sufficient condition for the Kolmogorov 4/5 law for stationary martingale solutions to the 3D Navier-Stokes equations

abstract: We prove that statistically stationary martingale solutions of the 3D Navier-Stokes equations on $\mathbb{T}^3$ subjected to white-in-time (colored-in-space) forcing satisfy the Kolmogorov 4/5 law (in an averaged sense and over a suitable inertial range) using only the assumption that the kinetic energy is $o(\nu^{-1})$ as $\nu \rightarrow 0$ (where $\nu$ is the inverse Reynolds number). This plays the role of a weak anomalous dissipation. No energy balance or additional regularity is assumed (aside from that satisfied by all martingale solutions from the energy inequality). If the force is statistically homogeneous, then any homogeneous martingale solution satisfies the spherically averaged 4/5 law pointwise in space. An additional hypothesis of approximate isotropy in the inertial range gives the traditional version of the Kolmogorov law. We demonstrate a necessary condition by proving that energy balance and an additional quantitative regularity estimate as $\nu \rightarrow 0$ imply that the 4/5 law (or any similar scaling law) cannot hold.


title: Are small radii of compact stars ruled out by GW170817/AT2017gfo?

abstract: The detection of GW170817 and its electromagnetic counterparts allows to constrain the equation of state of dense matter in new and complementary ways. Very stiff equations of state are ruled out by the upper limit on the average tidal deformability, $\tilde\Lambda \lesssim 800$, imposed by the detected gravitational wave signal. A lower limit, $\tilde\Lambda \gtrsim 400$, can also be extracted by considering the large amount of ejected matter which powers the kilonova AT2017gfo. By using several microscopic nucleonic equations of state, we first confirm the existence of a monotonic relation between $R_{1.5}$ (the radius of the $1.5M_{\odot}$ configuration) and $\tilde\Lambda$. This translates the limits on $\tilde\Lambda$ into limits on the radius: $11.8\,\text{km} \lesssim R_{1.5} \lesssim 13.1\,$km. We then show that the monotonic relation is violated, if a second branch of compact stars composed of quark matter exists, as in the two-families or the twin-stars scenarios. In particular, it is possible to fulfill the limits on $\tilde\Lambda$ while having $R_{1.5}$ significantly smaller than $12\,$km. In both those scenarios the event GW170817/AT2017gfo originates from the merger of a hadronic star and a star containing quark matter.


title: Constraining Primordial Black Holes with the EDGES 21-cm Absorption Signal

abstract: The EDGES experiment has recently measured an anomalous global 21-cm spectrum due to hydrogen absorptions at redshifts of about $z\sim 17$. Model independently, the unusually low temperature of baryons probed by this observable sets strong constraints on any physical process that transfers energy into the baryonic environment at such redshifts. Here we make use of the 21-cm spectrum to derive bounds on the energy injection due to a possible population of ${\cal O}(1-100) M_\odot$ primordial black holes, which induce a wide spectrum of radiation during the accretion of the surrounding gas. After calculating the total radiative intensity of a primordial black hole population, we estimate the amount of heat and ionisations produced in the baryonic gas and compute the resulting thermal history of the Universe with a modified version of RECFAST code. Finally, by imposing that the temperature of the gas at $z\sim 17$ does not exceed the indications of EDGES, we constrain the possible abundance of primordial black holes. Depending on uncertainties related to the accretion model, we find that ${\cal O}(10) M_\odot$ primordial black holes can only contribute to a fraction $f_{\rm PBH}<(1-10^{-3})$ of the total dark matter abundance.


title: Proactive Received Power Prediction Using Machine Learning and Depth Images for mmWave Networks

abstract: This study demonstrates the feasibility of the proactive received power prediction by leveraging spatiotemporal visual sensing information toward the reliable millimeter-wave (mmWave) networks. Since the received power on a mmWave link can attenuate aperiodically due to a human blockage, the long-term series of the future received power cannot be predicted by analyzing the received signals before the blockage occurs. We propose a novel mechanism that predicts a time series of the received power from the next moment to even several hundred milliseconds ahead. The key idea is to leverage the camera imagery and machine learning (ML). The time-sequential images can involve the spatial geometry and the mobility of obstacles representing the mmWave signal propagation. ML is used to build the prediction model from the dataset of sequential images labeled with the received power in several hundred milliseconds ahead of when each image is obtained. The simulation and experimental evaluations using IEEE 802.11ad devices and a depth camera show that the proposed mechanism employing convolutional LSTM predicted a time series of the received power in up to 500 ms ahead at an inference time of less than 3 ms with a root-mean-square error of 3.5 dB.


title: Model for Metal-Insulator Transition in Graphene Superlattices and Beyond

abstract: We propose a two-orbital Hubbard model on an emergent honeycomb lattice to describe the low-energy physics of twisted bilayer graphene. Our model provides a theoretical basis for studying metal-insulator transition, Landau level degeneracy lifting and unconventional superconductivity that are recently observed.


title: Electron correlations and pre-collision in the re-collision picture of high harmonic generation

abstract: We discuss the seminal three-step model and the recollision picture in the context of high harmonic generation in molecules. In particular, we stress the importance of multi-electron correlation during the first and the third of the three steps of the process: (1) the strong field ionization and (3) the recombination. We point out how accurate account for multi-electron correlations during the third, recombination, step allows one to gauge the importance of pre-collision: the term coined by J. H. Eberly to describe unusual pathways during the first, ionization, step.


title: Comparing Published Scientific Journal Articles to Their Pre-Print Versions -- Extended Version

abstract: Academic publishers claim that they add value to scholarly communications by coordinating reviews and contributing and enhancing text during publication. These contributions come at a considerable cost: U.S. academic libraries paid $1.7 billion for serial subscriptions in 2008 alone. Library budgets, in contrast, are flat and not able to keep pace with serial price inflation. We have investigated the publishers' value proposition by conducting a comparative study of pre-print papers from two distinct science, technology, and medicine (STM) corpora and their final published counterparts. This comparison had two working assumptions: 1) if the publishers' argument is valid, the text of a pre-print paper should vary measurably from its corresponding final published version, and 2) by applying standard similarity measures, we should be able to detect and quantify such differences. Our analysis revealed that the text contents of the scientific papers generally changed very little from their pre-print to final published versions. These findings contribute empirical indicators to discussions of the added value of commercial publishers and therefore should influence libraries' economic decisions regarding access to scholarly publications.


title: HAMLET: Interpretable Human And Machine co-LEarning Technique

abstract: Efficient label acquisition processes are key to obtaining robust classifiers. However, data labeling is often challenging and subject to high levels of label noise. This can arise even when classification targets are well defined, if instances to be labeled are more difficult than the prototypes used to define the class, leading to disagreements among the expert community. Here, we enable efficient training of deep neural networks. From low-confidence labels, we iteratively improve their quality by simultaneous learning of machines and experts. We call it Human And Machine co-LEarning Technique (HAMLET). Throughout the process, experts become more consistent, while the algorithm provides them with explainable feedback for confirmation. HAMLET uses a neural embedding function and a memory module filled with diverse reference embeddings from different classes. Its output includes classification labels and highly relevant reference embeddings as explanation. We took the study of brain monitoring at intensive care unit (ICU) as an application of HAMLET on continuous electroencephalography (cEEG) data. Although cEEG monitoring yields large volumes of data, labeling costs and difficulty make it hard to build a classifier. Additionally, while experts agree on the labels of clear-cut examples of cEEG patterns, labeling many real-world cEEG data can be extremely challenging. Thus, a large minority of sequences might be mislabeled. HAMLET has shown significant performance gain against deep learning and other baselines, increasing accuracy from 7.03% to 68.75% on challenging inputs. Besides improved performance, clinical experts confirmed the interpretability of those reference embeddings in helping explaining the classification results by HAMLET.


title: Bounded Context Switching for Valence Systems

abstract: We study valence systems, finite-control programs over infinite-state memories modeled in terms of graph monoids. Our contribution is a notion of bounded context switching (BCS). Valence systems generalize pushdowns, concurrent pushdowns, and Petri nets. In these settings, our definition conservatively generalizes existing notions. The main finding is that reachability within a bounded number of context switches is in NP, independent of the memory (the graph monoid). Our proof is genuinely algebraic, and therefore contributes a new way to think about BCS. In addition, we exhibit a class of storage mechanisms for which BCS reachability belongs to P.


title: MOrdReD: Memory-based Ordinal Regression Deep Neural Networks for Time Series Forecasting

abstract: Time series forecasting is ubiquitous in the modern world. Applications range from health care to astronomy, and include climate modelling, financial trading and monitoring of critical engineering equipment. To offer value over this range of activities, models must not only provide accurate forecasts, but also quantify and adjust their uncertainty over time. In this work, we directly tackle this task with a novel, fully end-to-end deep learning method for time series forecasting. By recasting time series forecasting as an ordinal regression task, we develop a principled methodology to assess long-term predictive uncertainty and describe rich multimodal, non-Gaussian behaviour, which arises regularly in applied settings. Notably, our framework is a wholly general-purpose approach that requires little to no user intervention to be used. We showcase this key feature in a large-scale benchmark test with 45 datasets drawn from both, a wide range of real-world application domains, as well as a comprehensive list of synthetic maps. This wide comparison encompasses state-of-the-art methods in both the Machine Learning and Statistics modelling literature, such as the Gaussian Process. We find that our approach does not only provide excellent predictive forecasts, shadowing true future values, but also allows us to infer valuable information, such as the predictive distribution of the occurrence of critical events of interest, accurately and reliably even over long time horizons.


title: The $\{1,s\}$-weighted Davenport constant in $C_n^k$

abstract: Let $G$ be a finite abelian group and let $\varnothing \neq A \subset \mathbb Z$. The $A$-weighted Davenport constant of $G$ is the smallest positive integer ${\sf D}_A(G)$ such that every sequence $x_1 \boldsymbol{\cdot} {\dots} \boldsymbol{\cdot} x_{{\sf D}_A(G)}$ over $G$ has a non-empty subsequence $(x_{j_i})_i$ such that ${\varepsilon_1} x_{j_1} + {\varepsilon_2} x_{j_2} + {\dots} + {\varepsilon_t} x_{j_t} = 0$ for some $\varepsilon_1, \varepsilon_2, {\dots}, \varepsilon_t \in A$. In this paper, we obtain both upper and lower bounds for ${\sf D}_{\{1,s\}}(C_n^k)$, where $C_n$ denotes the cyclic group of order $n$, $s^2 \equiv 1 \pmod n$ and $s \not\equiv \pm1 \pmod n$. These bounds become sharp in some "small" cases.


title: A note on Stokes' problem in dense granular media using the $\mu(I)$--rheology

abstract: The classical Stokes' problem describing the fluid motion due to a steadily moving infinite wall is revisited in the context of dense granular flows of mono-dispersed beads using the recently proposed $\mu(I)$--rheology. In Newtonian fluids, molecular diffusion brings about a self-similar velocity profile and the boundary layer in which the fluid motion takes place increases indefinitely with time $t$ as $\sqrt{\nu t}$, where $\nu$ is the kinematic viscosity. For a dense granular visco-plastic liquid, it is shown that the local shear stress, when properly rescaled, exhibits self-similar behaviour at short-time scales and it then rapidly evolves towards a steady-state solution. The resulting shear layer increases in thickness as $\sqrt{\nu_g t}$ analogous to a Newtonian fluid where $\nu_g$ is an equivalent granular kinematic viscosity depending not only on the intrinsic properties of the granular media such as grain diameter $d$, density $\rho$ and friction coefficients but also on the applied pressure $p_w$ at the moving wall and the solid fraction $\phi$ (constant). In addition, the $\mu(I)$--rheology indicates that this growth continues until reaching the steady-state boundary layer thickness $\delta_s = \beta_w (p_w/\phi \rho g )$, independent of the grain size, at about a finite time proportional to $\beta_w^2 (p_w/\rho g d)^{3/2} \sqrt{d/g}$, where $g$ is the acceleration due to gravity and $\beta_w = (\tau_w - \tau_s)/\tau_s$ is the relative surplus of the steady-state wall shear-stress $\tau_w$ over the critical wall shear stress $\tau_s$ (yield stress) that is needed to bring the granular media into motion... (see article for a complete abstract).


title: A library of second-order models for synchronous machines

abstract: This paper presents a library of second-order models for synchronous machines that can be utilized in power system dynamic performance analysis and control design tasks. The models have a similar structure to the classical model in that they consist of two dynamic states, the power angle and the angular speed. However, unlike the classical model, the models find applications beyond first swing stability analysis; for example, they can also be utilized in transient stability studies. The models are developed through a systematic reduction of a nineteenth-order model, using singular perturbation techniques, and they are validated by comparing their voltage, frequency, and phase profiles with that of the high-order model and that of the classical model.


title: Optical conductivity of the Weyl semimetal NbP

abstract: The optical properties of (001)-oriented NbP single crystals have been studied in a wide spectral range from 6 meV to 3 eV from room temperature down to 10 K. The itinerant carriers lead to a Drude-like contribution to the optical response; we can further identify two pronounced phonon modes and interband transitions starting already at rather low frequencies. By comparing our experimental findings to the calculated interband optical conductivity, we can assign the features observed in the measured conductivity to certain interband transitions. In particular, we find that transitions between the electronic bands spilt by spin-orbit coupling dominate the interband conductivity of NbP below 100 meV. At low temperatures, the momentum-relaxing scattering rate of the itinerant carriers in NbP is very small, leading to macroscopic characteristic length scales of the momentum relaxation of approximately 0.5 $\mu$m.


title: A many-sorted polyadic modal logic

abstract: This paper presents a many-sorted polyadic modal logic that generalizes some of the existing approaches. The algebraic semantics has led us to a many-sorted generalization of boolean algebras with operators, for which we prove the analogue of the J\'onsson-Tarski theorem. While the transition from the mono-sorted logic to many-sorted one is a smooth process, we see our system as a step towards deepening the connection between modal logic and program verification, since our system can be seen as the propositional fragment of Matching logic, a first-order logic for specifying and reasoning about programs.


title: Secure and Reliable Biometric Access Control for Resource-Constrained Systems and IoT

abstract: With the emergence of the Internet-of-Things (IoT), there is a growing need for access control and data protection on low-power, pervasive devices. Biometric-based authentication is promising for IoT due to its convenient nature and lower susceptibility to attacks. However, the costs associated with biometric processing and template protection are nontrivial for smart cards, key fobs, and so forth. In this paper, we discuss the security, cost, and utility of biometric systems and develop two major frameworks for improving them. First, we introduce a new framework for implementing biometric systems based on physical unclonable functions (PUFs) and hardware obfuscation that, unlike traditional software approaches, does not require nonvolatile storage of a biometric template/key. Aside from reducing the risk of compromising the biometric, the nature of obfuscation also provides protection against access control circumvention via malware and fault injection. The PUF provides non-invertibility and non-linkability. Second, a major requirement of the proposed PUF/obfuscation approach is that a reliable (robust) key be generated from the users input biometric. We propose a noiseaware biometric quantization framework capable of generating unique, reliable keys with reduced enrollment time and denoising costs. Finally, we conduct several case studies. In the first, the proposed noise-aware approach is compared to our previous approach for multiple biometric modalities, including popular ones (fingerprint and iris) and emerging cardiovascular ones (ECG and PPG). The results show that ECG provides the best tradeoff between reliability, key length, entropy, and cost. In the second and third case studies, we demonstrate how reliability, denoising costs, and enrollment times can be simultaneously improved by modeling subject intra-variations for ECG.


title: Cost optimal scenarios of a future highly renewable European electricity system: Exploring the influence of weather data, cost parameters and policy constraints

abstract: Cost optimal scenarios derived from models of a highly renewable electricity system depend on the specific input data, cost assumptions and system constraints. Here this influence is studied using a techno-economic optimisation model for a networked system of 30 European countries, taking into account the capacity investment and operation of wind, solar, hydroelectricity, natural gas power generation, transmission, and different storage options. A considerable robustness of total system costs to the input weather data and to moderate changes in the cost assumptions is observed. Flat directions in the optimisation landscape around cost-optimal configurations often allow system planners to choose between different technology options without a significant increase in total costs, for instance by replacing onshore with offshore wind power capacity in case of public acceptance issues. Exploring a range of carbon dioxide emission limits shows that for scenarios with moderate transmission expansion, a reduction of around 57% compared to 1990 levels is already cost optimal. For stricter carbon dioxide limits, power generated from gas turbines is at first replaced by generation from increasing renewable capacities. Non-hydro storage capacities are only built for low-emission scenarios, in order to provide the necessary flexibility to meet peaks in the residual load.


title: On the Diffusion of Sticky Particles in 1D

abstract: The 1D Ising model is the simplest Hamiltonian-based model in statistical mechanics. The sim- plest interacting particle process is the Symmetric Exclusion Process (SEP), a 1D lattice gas of particles that hop symmetrically and cannot overlap. Combining the two gives a model for sticky particle diffusion, SPM, which is described here. SPM dynamics are based on SEP with short-range interaction, allowing flow due to non-equilibrium boundary conditions. We prove that SPM is also a detailed-balance respecting, particle-conserving, Monte Carlo (MC) description of the Ising model. Neither the Ising model nor SEP have a phase transition in 1D, but the SPM exhibits a non-equilibrium transition from a diffusing to a blocked state as stickiness increases. This transition manifests in peaks in the MC density fluctuation, a change in the dependency of flow-rate on stickiness, and odd structure in the eigenspectrum of the transition rate matrix. We derive and solve a fully non-linear, analytic, mean-field solution, which has a crossover from a positive to a negative diffusion constant where the MC shows the transition. The negative diffusion constant in fact indicates a breakdown of the mean-field approximation, with close to zero flow and breaking into a two-phase mixture, and thus the mean field theory successfully predicts its own demise. We also present an analytic solution for the flow via direct analysis of the transition rate matrix. The simplicity of the model suggest a wide range of possible applications.


title: Robust principal components for irregularly spaced longitudinal data

abstract: Consider longitudinal data $x_{ij},$ with $i=1,...,n$ and $j=1,...,p_{i},$ where $x_{ij}$ is the $j-$th observation of the random function $X_{i}\left( .\right) $ observed at time $t_{j}.$ The goal of this paper is to develop a parsimonious representation of the data by a linear combination of a set of $q$ smooth functions $H_{k}\left( .\right) $ ($k=1,..,q)$ in the sense that $x_{ij}\approx\mu_{j}+\sum_{k=1}^{q}\beta_{ki}H_{k}\left( t_{j}\right) ,$ such that it fulfills three goals: it is resistant to atypical $X_{i}$'s ('case contamination'), it is resistant to isolated gross errors at some $t_{ij}$ ('cell contamination'), and it can be applied when some of the $x_{ij}$ are missing ('irregularly spaced' ---or 'incomplete'-- data). Two approaches will be proposed for this problem. One deals with the three goals stated above, and is based on ideas similar to MM-estimation (Yohai 1987). The other is a simple and fast estimator which can be applied to complete data with case- and cellwise contamination, and is based on applying a standard robust principal components estimate and smoothing the principal directions. Experiments with real and simulated data suggest that with complete data the simple estimator outperforms its competitors, while the MM estimator is competitive for incomplete data. Keywords: Principal components, MM-estimator, longitudinal .data, B-splines, incomplete data.


title: Observation of $\psi(3686) \to \eta' e^+ e^- $

abstract: Using a data sample of $448.1 \times 10^6$ $\psi(3686)$ events collected with the BESIII detector at the BEPCII collider, we report the first observation of the electromagnetic Dalitz decay $\psi(3686) \to \eta' e^+ e^-$, with significances of 7.0$\sigma$ and 6.3$\sigma$ when reconstructing the $\eta'$ meson via its decay modes $\eta'\to\gamma \pi^+ \pi^-$ and $\eta'\to\pi^+\pi^-\eta$ ($\eta \to \gamma\gamma$), respectively. The weighted average branching fraction is determined to be $\mathcal{B}(\psi(3686) \to \eta' e^+ e^-)= (1.90 \pm 0.25 \pm 0.11) \times 10^{-6}$, where the first uncertainty is statistical and the second systematic.


title: On the Runtime Analysis of the Clearing Diversity-Preserving Mechanism

abstract: Clearing is a niching method inspired by the principle of assigning the available resources among a niche to a single individual. The clearing procedure supplies these resources only to the best individual of each niche: the winner. So far, its analysis has been focused on experimental approaches that have shown that clearing is a powerful diversity-preserving mechanism. Using rigorous runtime analysis to explain how and why it is a powerful method, we prove that a mutation-based evolutionary algorithm with a large enough population size, and a phenotypic distance function always succeeds in optimising all functions of unitation for small niches in polynomial time, while a genotypic distance function requires exponential time. Finally, we prove that with phenotypic and genotypic distances clearing is able to find both optima for Twomax and several general classes of bimodal functions in polynomial expected time. We use empirical analysis to highlight some of the characteristics that makes it a useful mechanism and to support the theoretical results.


title: Random Nodal Lengths and Wiener Chaos

abstract: In this survey we collect some of the recent results on the "nodal geometry" of random eigenfunctions on Riemannian surfaces. We focus on the asymptotic behavior, for high energy levels, of the nodal length of Gaussian Laplace eigenfunctions on the torus (arithmetic random waves) and on the sphere (random spherical harmonics). We give some insight on both Berry's cancellation phenomenon and the nature of nodal length second order fluctuations (non-Gaussian on the torus and Gaussian on the sphere) in terms of chaotic components. Finally we consider the general case of monochromatic random waves, i.e. Gaussian random linear combination of eigenfunctions of the Laplacian on a compact Riemannian surface with frequencies from a short interval, whose scaling limit is Berry's Random Wave Model. For the latter we present some recent results on the asymptotic distribution of its nodal length in the high energy limit (equivalently, for growing domains).


title: Parameterized Intractability of Even Set and Shortest Vector Problem from Gap-ETH

abstract: The $k$-Even Set problem is a parameterized variant of the Minimum Distance Problem of linear codes over $\mathbb F_2$, which can be stated as follows: given a generator matrix $\mathbf A$ and an integer $k$, determine whether the code generated by $\mathbf A$ has distance at most $k$. Here, $k$ is the parameter of the problem. The question of whether $k$-Even Set is fixed parameter tractable (FPT) has been repeatedly raised in literature and has earned its place in Downey and Fellows' book (2013) as one of the "most infamous" open problems in the field of Parameterized Complexity. In this work, we show that $k$-Even Set does not admit FPT algorithms under the (randomized) Gap Exponential Time Hypothesis (Gap-ETH) [Dinur'16, Manurangsi-Raghavendra'16]. In fact, our result rules out not only exact FPT algorithms, but also any constant factor FPT approximation algorithms for the problem. Furthermore, our result holds even under the following weaker assumption, which is also known as the Parameterized Inapproximability Hypothesis (PIH) [Lokshtanov et al.'17]: no (randomized) FPT algorithm can distinguish a satisfiable 2CSP instance from one which is only $0.99$-satisfiable (where the parameter is the number of variables). We also consider the parameterized $k$-Shortest Vector Problem (SVP), in which we are given a lattice whose basis vectors are integral and an integer $k$, and the goal is to determine whether the norm of the shortest vector (in the $\ell_p$ norm for some fixed $p$) is at most $k$. Similar to $k$-Even Set, this problem is also a long-standing open problem in the field of Parameterized Complexity. We show that, for any $p > 1$, $k$-SVP is hard to approximate (in FPT time) to some constant factor, assuming PIH. Furthermore, for the case of $p = 2$, the inapproximability factor can be amplified to any constant.


title: The Krzy\.z Conjecture and an Entropy Conjecture

abstract: We show that if the minimum entropy for a polynomial with roots on the unit circle is attained by polynomials with equally spaced roots, then, under a generic hypothesis about the nature of the extremum, the Krzyz conjecture on the maximum modulus of the Taylor coefficients of a holomorphic function that maps the disk to the punctured disk is true.


title: On the Importance of Stereo for Accurate Depth Estimation: An Efficient Semi-Supervised Deep Neural Network Approach

abstract: We revisit the problem of visual depth estimation in the context of autonomous vehicles. Despite the progress on monocular depth estimation in recent years, we show that the gap between monocular and stereo depth accuracy remains large$-$a particularly relevant result due to the prevalent reliance upon monocular cameras by vehicles that are expected to be self-driving. We argue that the challenges of removing this gap are significant, owing to fundamental limitations of monocular vision. As a result, we focus our efforts on depth estimation by stereo. We propose a novel semi-supervised learning approach to training a deep stereo neural network, along with a novel architecture containing a machine-learned argmax layer and a custom runtime (that will be shared publicly) that enables a smaller version of our stereo DNN to run on an embedded GPU. Competitive results are shown on the KITTI 2015 stereo dataset. We also evaluate the recent progress of stereo algorithms by measuring the impact upon accuracy of various design criteria.


title: CliCR: A Dataset of Clinical Case Reports for Machine Reading Comprehension

abstract: We present a new dataset for machine comprehension in the medical domain. Our dataset uses clinical case reports with around 100,000 gap-filling queries about these cases. We apply several baselines and state-of-the-art neural readers to the dataset, and observe a considerable gap in performance (20% F1) between the best human and machine readers. We analyze the skills required for successful answering and show how reader performance varies depending on the applicable skills. We find that inferences using domain knowledge and object tracking are the most frequently required skills, and that recognizing omitted information and spatio-temporal reasoning are the most difficult for the machines.


title: Efficient algorithms to discover alterations with complementary functional association in cancer

abstract: Recent large cancer studies have measured somatic alterations in an unprecedented number of tumours. These large datasets allow the identification of cancer-related sets of genetic alterations by identifying relevant combinatorial patterns. Among such patterns, mutual exclusivity has been employed by several recent methods that have shown its effectivenes in characterizing gene sets associated to cancer. Mutual exclusivity arises because of the complementarity, at the functional level, of alterations in genes which are part of a group (e.g., a pathway) performing a given function. The availability of quantitative target profiles, from genetic perturbations or from clinical phenotypes, provides additional information that can be leveraged to improve the identification of cancer related gene sets by discovering groups with complementary functional associations with such targets. In this work we study the problem of finding groups of mutually exclusive alterations associated with a quantitative (functional) target. We propose a combinatorial formulation for the problem, and prove that the associated computation problem is computationally hard. We design two algorithms to solve the problem and implement them in our tool UNCOVER. We provide analytic evidence of the effectiveness of UNCOVER in finding high-quality solutions and show experimentally that UNCOVER finds sets of alterations significantly associated with functional targets in a variety of scenarios. In addition, our algorithms are much faster than the state-of-the-art, allowing the analysis of large datasets of thousands of target profiles from cancer cell lines. We show that on one such dataset from project Achilles our methods identify several significant gene sets with complementary functional associations with targets.


title: 3D Human Pose Estimation in the Wild by Adversarial Learning

abstract: Recently, remarkable advances have been achieved in 3D human pose estimation from monocular images because of the powerful Deep Convolutional Neural Networks (DCNNs). Despite their success on large-scale datasets collected in the constrained lab environment, it is difficult to obtain the 3D pose annotations for in-the-wild images. Therefore, 3D human pose estimation in the wild is still a challenge. In this paper, we propose an adversarial learning framework, which distills the 3D human pose structures learned from the fully annotated dataset to in-the-wild images with only 2D pose annotations. Instead of defining hard-coded rules to constrain the pose estimation results, we design a novel multi-source discriminator to distinguish the predicted 3D poses from the ground-truth, which helps to enforce the pose estimator to generate anthropometrically valid poses even with images in the wild. We also observe that a carefully designed information source for the discriminator is essential to boost the performance. Thus, we design a geometric descriptor, which computes the pairwise relative locations and distances between body joints, as a new information source for the discriminator. The efficacy of our adversarial learning framework with the new geometric descriptor has been demonstrated through extensive experiments on widely used public benchmarks. Our approach significantly improves the performance compared with previous state-of-the-art approaches.


title: Proton Colliders at the Energy Frontier

abstract: Since the CERN ISR, hadron colliders have defined the energy frontier. Noteworthy are the conversion of the Super Proton Synchrotron (SPS) into a proton-antiproton collider, the Tevatron collider, as well as the abandoned SSC in the United States. Hadron colliders are likely to determine the pace of particle-physics progress also during the next hundred years. Discoveries at past hadron colliders were essential for establishing the so-called Standard Model of particle physics. The world's present flagship collider, the LHC, including its high-luminosity upgrade HL-LHC, is set to operate through the second half of the 2030's. Further increases of the energy reach during the 21st century require another, still more powerful hadron collider. Three options for a next hadron collider are presently under investigation. The Future Circular Collider (FCC) study, hosted by CERN, is designing a 100 TeV collider, to be installed inside a new 100 km tunnel in the Lake Geneva basin. A similar 100 km collider, called SppC, is being pursued by CAS-IHEP in China. In either machine, for the first time in hadron storage rings, synchrotron radiation damping will be significant. In parallel, the synchrotron-radiation power emitted inside the cold magnets becomes an important design constraint. One important difference between FCC and SppC is the magnet technology. FCC uses 16 Tesla magnets based on Nb3Sn superconductor, while SppC magnets shall be realized with cables made from iron-based high-temperature superconductor. Initially the SppC magnets are assumed to provide a more moderate dipole field of 12 T, but they can later be pushed to a final ultimate field of 24 T. A third collider presently under study is the High-Energy LHC (HE-LHC), which is a higher energy collider in the existing LHC tunnel, exploiting the FCC magnet technology in order to double the LHC energy at significantly higher luminosity.


title: ExoMol molecular line lists - XXVI: spectra of SH and NS

abstract: Line lists for the sulphur-containing molecules SH (the mercapto radical) and NS are computed as part of the ExoMol project. These line lists consider transitions within the $X$ ${}^2\Pi$ ground state for $^{32}$SH, $^{33}$SH, $^{34}$SH and $^{\text{32}}$SD, and $^{14}$N$^{32}$S, $^{14}$N$^{33}$S, $^{14}$N$^{34}$S, $^{14}$N$^{36}$S and $^{15}$N$^{32}$S. Ab initio potential energy (PEC) and spin-orbit coupling (SOC) curves are computed and then improved by fitting to experimentally observed transitions. Fully ab initio dipole moment curves (DMCs) computed at high level of theory are used to produce the final line lists. For SH, our fit gives a root-mean-square (rms) error of 0.03 cm$^{-1}$ between the observed ($v_{\rm max}=4$, $J_{\rm max} = 34.5$) and calculated transitions wavenumbers; this is extrapolated such that all $X$ $^2\Pi$ rotational-vibrational-electronic (rovibronic) bound states are considered. For $^{\text{32}}$SH the resulting line list contains about 81000 transitions and 2300 rovibronic states, considering levels up to $v_{\rm max} = 14$ and $J_{\rm max} = 60.5$. For NS the refinement used a combination of experimentally determined frequencies and energy levels and led to an rms fitting error of 0.002 cm$^{-1}$. Each NS calculated line list includes around 2.8 million transitions and 31000 rovibronic states with a vibrational range up to $v=53$ and rotational range to $J=235.5$, which covers up to 23000 cm$^{-1}$. Both line lists should be complete for temperatures up to 5000 K. Example spectra simulated using this line list are shown and comparisons made to the existing data in the CDMS database. The line lists are available from the CDS (http://cdsarc.u-strasbg.fr) and ExoMol (www.exomol.com) data bases.

categories: astro-ph.SR astro-ph.EP astro-ph.GA physics.chem-ph


title: Correcting the Mistaken Identification of Nonequilibrium Microscopic Work

abstract: The energy change dE_k for the kth microstate is erroneously equated with the external work done on the microstate. It ignores the ubiquitous internal energy change d_iW_k due to force imbalance between the internal and external forces. We show that this contribution is present even in a reversible process, which is a surprise. We show that the correct identification is dE_k=-dW_k, where dW_k is the generalized work done by the microstate. We prove that the thermodynamic average of the internal work gives dissipation and is not captured by the external work. The latter effectively sets d_iW_k =0 and results in no dissipation. Using dW_k to account for irreversibility, we obtain a new work relation that works even for free expansion, where the Jarzynski equality fails. In the new work relation, dW_k depends only on the energies of the initial and final states and not on the actual process. This makes the new relation very different from the Jarzynski equality. The correction has far-reaching consequences and requires reassessment of current applications of external work in theoretical physics.


title: Neutron star merger GW170817 strongly constrains doubly coupled bigravity

abstract: We study the implications of the recent detection of gravitational waves emitted by a pair of merging neutron stars and their electromagnetic counterpart, events GW170817 and GRB170817A, on the viability of the doubly coupled bimetric models of cosmic evolution, where the two metrics couple directly to matter through a composite, effective metric. We demonstrate that the bounds on the speed of gravitational waves place strong constraints on the doubly coupled models, forcing either the two metrics to be proportional at the background level or the models to become singly coupled. Proportional backgrounds are particularly interesting as they provide stable cosmological solutions with phenomenologies equivalent to that of $\Lambda$CDM at the background level as well as for linear perturbations, while nonlinearities are expected to show deviations from the standard model.


title: Alignment-Free Sequence Analysis and Applications

abstract: Genome and metagenome comparisons based on large amounts of next-generation sequencing (NGS) data pose significant challenges for alignment-based approaches due to the huge data size and the relatively short length of the reads. Alignment-free approaches based on the counts of word patterns in NGS data do not depend on the complete genome and are generally computationally efficient. Thus, they contribute significantly to genome and metagenome comparison. Recently, novel statistical approaches have been developed for the comparison of both long and shotgun sequences. These approaches have been applied to many problems including the comparison of gene regulatory regions, genome sequences, metagenomes, binning contigs in metagenomic data, identification of virus-host interactions, and detection of horizontal gene transfers. We provide an updated review of these applications and other related developments of word-count based approaches for alignment-free sequence analysis.


title: On the search for Galactic supernova remnant PeVatrons with current TeV instruments

abstract: The supernova remnant hypothesis for the origin of Galactic cosmic rays has passed several tests, but the firm identification of a supernova remnant pevatron, considered to be a decisive step to prove the hypothesis, is still missing. While a lot of hope has been placed in next-generation instruments operating in the multi-TeV range, it is possible that current gamma-ray instruments, operating in the TeV range, could pinpoint these objects or, most likely, identify a number of promising targets for instruments of next generation. Starting from the assumption that supernova remnants are indeed the sources of Galactic cosmic rays, and therefore must be pevatrons for some fraction of their lifetime, we investigate the ability of current instruments to detect such objects, or to identify the most promising candidates.


title: From nonlocal to local Cahn-Hilliard equation

abstract: In this paper we prove the convergence of a nonlocal version of the Cahn-Hilliard equation to its local counterpart as the nonlocal convolution kernel is scaled using suitable approximations of a Dirac delta in a periodic boundary conditions setting. This convergence result strongly relies on the dynamics of the problem. More precisely, the $H^{-1}$-gradient flow structure of the equation allows to deduce uniform $H^1$ estimates for solutions of the nonlocal Cahn-Hilliard equation and, together with a Poincar\'e type inequality by Ponce, provides the compactness argument that allows to prove the convergence result.


title: Resilient Active Information Gathering with Mobile Robots

abstract: Applications of safety, security, and rescue in robotics, such as multi-robot target tracking, involve the execution of information acquisition tasks by teams of mobile robots. However, in failure-prone or adversarial environments, robots get attacked, their communication channels get jammed, and their sensors may fail, resulting in the withdrawal of robots from the collective task, and consequently the inability of the remaining active robots to coordinate with each other. As a result, traditional design paradigms become insufficient and, in contrast, resilient designs against system-wide failures and attacks become important. In general, resilient design problems are hard, and even though they often involve objective functions that are monotone or submodular, scalable approximation algorithms for their solution have been hitherto unknown. In this paper, we provide the first algorithm, enabling the following capabilities: minimal communication, i.e., the algorithm is executed by the robots based only on minimal communication between them; system-wide resiliency, i.e., the algorithm is valid for any number of denial-of-service attacks and failures; and provable approximation performance, i.e., the algorithm ensures for all monotone (and not necessarily submodular) objective functions a solution that is finitely close to the optimal. We quantify our algorithm's approximation performance using a notion of curvature for monotone set functions. We support our theoretical analyses with simulated and real-world experiments, by considering an active information gathering scenario, namely, multi-robot target tracking.


title: Lessons from HAWC PWNe observations: the diffusion constant is not a constant; Pulsars remain the likeliest sources of the anomalous positron fraction; Cosmic rays are trapped for long periods of time in pockets of inefficient diffusion

abstract: Recent TeV observations of nearby pulsars with the HAWC telescope have been interpreted as evidence that the diffusion of high-energy electrons and positrons within pulsar wind nebulae is highly inefficient compared to the rest of the interstellar medium. If the diffusion coefficient well outside the nebula is close to the value inferred for the region inside the nebula, high-energy electrons and positrons produced by the two observed pulsars could not contribute significantly to the local measured cosmic-ray flux. The HAWC Collaboration thus concluded that, under the assumption of isotropic and homogeneous diffusion, the two pulsars are ruled out as sources of the anomalous high-energy positron flux. Here, we argue that since the diffusion coefficient is likely not spatially homogeneous, the assumption leading to this conclusion is flawed. We solve the diffusion equation with a radially dependent diffusion coefficient, and show that the pulsars observed by HAWC produce potentially perfect matches to the observed high-energy positron fluxes. We also study the implications of inefficient diffusion within pulsar wind nebulae on Galactic scales, and show that cosmic rays are likely to have very long residence times in regions of inefficient diffusion. We describe how this prediction can be tested with studies of the diffuse Galactic emission.


title: Ground-state properties of dilute Bose systems with synthetic dispersion laws

abstract: Experimental advances in synthesizing spin-orbit couplings in cold atomic Bose gases promise to create single-particle dispersion laws featuring energy minima that are degenerate on a ring or a sphere in momentum space. We show that for arbitrary space dimensionality the ground-state properties of a dilute system of spin-orbit coupled Bose particles with such dispersion and short-range repulsive interactions are universal: the chemical potential exhibits a quadratic dependence on the particle density as found in a one-dimensional free Fermi gas.


title: Domain transfer convolutional attribute embedding

abstract: In this paper, we study the problem of transfer learning with the attribute data. In the transfer learning problem, we want to leverage the data of the auxiliary and the target domains to build an effective model for the classification problem in the target domain. Meanwhile, the attributes are naturally stable cross different domains. This strongly motives us to learn effective domain transfer attribute representations. To this end, we proposed to embed the attributes of the data to a common space by using the powerful convolutional neural network (CNN) model. The convolutional representations of the data points are mapped to the corresponding attributes so that they can be effective embedding of the attributes. We also represent the data of different domains by a domain-independent CNN, ant a domain-specific CNN, and combine their outputs with the attribute embedding to build the classification model. An joint learning framework is constructed to minimize the classification errors, the attribute mapping error, the mismatching of the domain-independent representations cross different domains, and to encourage the the neighborhood smoothness of representations in the target domain. The minimization problem is solved by an iterative algorithm based on gradient descent. Experiments over benchmark data sets of person re-identification, bankruptcy prediction, and spam email detection, show the effectiveness of the proposed method.


title: Early-Universe Constraints on Dark Matter-Baryon Scattering and their Implications for a Global 21cm Signal

abstract: We present and compare several cosmological constraints on the cross section for elastic scattering between dark matter (DM) and baryons, for cross sections with a range of power-law dependences on the DM-baryon relative velocity $v$, especially focusing on the case of $\sigma \propto v^{-4}$. We study constraints spanning a wide range of epochs in cosmological history, from pre-recombination distortions to the blackbody spectrum and anisotropies of the cosmic microwave background (CMB), to modifications to the intergalactic medium temperature and the resulting 21cm signal, and discuss the allowed signals in the latter channels given the constraints from the former. We improve previous constraints on DM-baryon scattering from the CMB anisotropies, demonstrate via principal component analysis that the effect on the CMB can be written as a simple function of DM mass, and map out the redshifts dominating this signal. We show that given high-redshift constraints on DM-baryon scattering, a $v^{-4}$ scaling of the cross section for light DM would be sufficient to explain the deep 21cm absorption trough recently claimed by the EDGES experiment, if 100% of the DM scatters with baryons. For millicharged DM models proposed to explain the observation, where only a small fraction of the DM interacts, we estimate that a PIXIE-like future experiment measuring CMB spectral distortion could test the relevant parameter space.


title: A Scalable Empirical Bayes Approach to Variable Selection in Generalized Linear Models

abstract: A new empirical Bayes approach to variable selection in the context of generalized linear models is developed. The proposed algorithm scales to situations in which the number of putative explanatory variables is very large, possibly much larger than the number of responses. The coefficients in the linear predictor are modeled as a three-component mixture allowing the explanatory variables to have a random positive effect on the response, a random negative effect, or no effect. A key assumption is that only a small (but unknown) fraction of the candidate variables have a non-zero effect. This assumption, in addition to treating the coefficients as random effects facilitates an approach that is computationally efficient. In particular, the number of parameters that have to be estimated is small, and remains constant regardless of the number of explanatory variables. The model parameters are estimated using a Generalized Alternating Maximization algorithm which is scalable, and leads to significantly faster convergence compared with simulation-based fully Bayesian methods.


title: Antigravity from a spacetime defect

abstract: We argue that there may exist spacetime defects embedded in Minkowski spacetime, which have negative active gravitational mass. One such spacetime defect then repels a test particle, corresponding to what may be called "antigravity."


title: DJAM: distributed Jacobi asynchronous method for learning personal models

abstract: Processing data collected by a network of agents often boils down to solving an optimization problem. The distributed nature of these problems calls for methods that are, themselves, distributed. While most collaborative learning problems require agents to reach a common (or consensus) model, there are situations in which the consensus solution may not be optimal. For instance, agents may want to reach a compromise between agreeing with their neighbors and minimizing a personal loss function. We present DJAM, a Jacobi-like distributed algorithm for learning personalized models. This method is implementation-friendly: it has no hyperparameters that need tuning, it is asynchronous, and its updates only require single-neighbor interactions. We prove that DJAM converges with probability one to the solution, provided that the personal loss functions are strongly convex and have Lipschitz gradient. We then give evidence that DJAM is on par with state-of-the-art methods: our method reaches a solution with error similar to the error of a carefully tuned ADMM in about the same number of single-neighbor interactions.


title: Note on the truncated generalizations of Gauss' square exponent theorem

abstract: In this note, we investigate J.-C. Liu's work on truncated Gauss' square exponent theorem and obtain more truncations. We also discuss some possible multiple summation extensions of Liu's results.


title: Too Hot, Too Cold or Just Right? Implications of a 21-cm Signal for Dark Matter Annihilation and Decay

abstract: Measurements of the temperature of the baryons at the end of the cosmic dark ages can potentially set very precise constraints on energy injection from exotic sources, such as annihilation or decay of the dark matter. However, additional effects that lower the gas temperature can substantially weaken the expected constraints on exotic energy injection, whereas additional radiation backgrounds can conceal the effect of an increased gas temperature in measurements of the 21-cm hyperfine transition of neutral hydrogen. Motivated in part by recent claims of a detection of 21-cm absorption from a redshift of 17 by the EDGES experiment, we derive the constraints on dark matter annihilation and decay that can be placed in the presence of extra radiation backgrounds or effects that modify the gas temperature, such as dark matter-baryon scattering and early baryon-photon decoupling. We find that if the EDGES observation is confirmed, then constraints on light dark matter decaying or annihilating to electrons will in most scenarios be stronger than existing state-of-the-art limits from the cosmic microwave background, potentially by several orders of magnitude. More generally, our results allow mapping any future measurement of the global 21-cm signal into constraints on dark matter annihilation and decay, within the broad range of scenarios we consider.


title: Comparison of local density functionals based on electron gas and finite systems

abstract: A widely used approximation to the exchange-correlation functional in density functional theory is the local density approximation (LDA), typically derived from the properties of the homogeneous electron gas (HEG). We previously introduced a set of alternative LDAs constructed from one-dimensional systems of one, two, and three electrons that resemble the HEG within a finite region. We now construct a HEG-based LDA appropriate for spinless electrons in one dimension and find that it is remarkably similar to the finite LDAs. As expected, all LDAs are inadequate in low-density systems where correlation is strong. However, exploring the small but significant differences between the functionals, we find that the finite LDAs give better densities and energies in high-density exchange-dominated systems, arising partly from a better description of the self-interaction correction.


title: The Poisson bracket invariant on surfaces

abstract: We study the Poisson bracket invariant, which measures the level of Poisson noncommutativity of a smooth partition of unity, on closed symplectic surfaces. Motivated by a general conjecture of Polterovich and building on preliminary work of Buhovsky--Tanny, we prove that for any smooth partition of unity subordinate to an open cover by discs of area at most $c$, and under some localization condition on the cover when the surface is a sphere, then the product of the Poisson bracket invariant with $c$ is bounded from below by a universal constant. Similar results were obtained recently by Buhovsky--Logunov--Tanny for open covers consisting of displaceable sets on all closed surfaces, and their approach was extended by Shi--Lu to open covers by nondisplaceable discs. We investigate the sharpness of all these results.


title: Origin of Mott insulating behavior and superconductivity in twisted bilayer graphene

abstract: A remarkable recent experiment has observed Mott insulator and proximate superconductor phases in twisted bilayer graphene when electrons partly fill a nearly flat mini-band that arises a `magic' twist angle. However, the nature of the Mott insulator, origin of superconductivity and an effective low energy model remain to be determined. We propose a Mott insulator with intervalley coherence that spontaneously breaks U(1) valley symmetry, and describe a mechanism that selects this order over the competing magnetically ordered states favored by the Hunds coupling. We also identify symmetry related features of the nearly flat band that are key to understanding the strong correlation physics and constrain any tight binding description. First, although the charge density is concentrated on the triangular lattice sites of the moir$\text{\'e }$ pattern, the Wannier states of the tight-binding model must be centered on different sites which form a honeycomb lattice. Next, spatially localizing electrons derived from the nearly flat band necessarily breaks valley and other symmetries within any mean-field treatment, which is suggestive of a valley-ordered Mott state, and also dictates that additional symmetry breaking is present to remove symmetry-enforced band contacts. Tight-binding models describing the nearly flat mini-band are derived, which highlight the importance of further neighbor hopping and interactions. We discuss consequences of this picture for superconducting states obtained on doping the valley ordered Mott insulator. We show how important features of the experimental phenomenology may be explained and suggest a number of further experiments for the future. We also describe a model for correlated states in trilayer graphene heterostructures and contrast it with the bilayer case.


title: Non-minimally Coupled Pseudoscalar Inflaton

abstract: We consider a scenario in which the inflaton $\phi$ is a pseudoscalar field non-minimally coupled to gravity through a term of the form ${\cal X} R \phi^2$. The pseudoscalar is also coupled to a $U(1)$ gauge field (or an ensemble of ${\cal N}$ gauge fields) through an axial coupling of the form $\phi F \tilde{F}$. After M. M. Anber and L. Sorbo, Phys. Rev. D 81, 043534 (2010), Ref. [1], it is well known that this axial coupling leads to a production of gauge particles which acts as a friction term in the dynamics of the inflaton, producing a slow-roll regime even in presence of a steep potential. A remarkable result in this scenario, is that the spectrum of the chiral gravitational waves sourced by the scalar-gauge field interplay can be enhanced due to the non-minimal coupling with gravity, leading to measurable signatures, while maintaining agreement with current observational constraints on $n_s$ and $r$. The inclusion of non-minimal coupling could be helpful to alleviate tensions with non-Gaussianity bounds in models including axial couplings.


title: On the Formation and Properties of Fluid Shocks and Collisionless Shock Waves in Astrophysical Plasmas

abstract: When two plasmas collide, their interaction can be mediated by collisionless plasma instabilities or binary collisions between particles of each shell. By comparing the maximum growth rate of the collisionless instabilities with the collision frequency between particles of the shells, we determine the critical density separating the collisionless formation from the collisional formation of the resulting shock waves. This critical density is also the density beyond which the shock downstream is field free, as plasma instabilities do not have time to develop electromagnetic patterns. We further determine the conditions on the shells initial density and velocity for the downstream to be collisional. If these quantities fulfill the determined conditions, the collisionality of the downstream also prevents the shock from accelerating particles or generating strong magnetic fields. We compare the speed of sound with the relative speed of collision between the two shells, thus determining the portion of the parameters space where strong shock formation is possible for both classical and degenerate plasmas. Finally, we discuss the observational consequences in several astrophysical settings.

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


title: English verb regularization in books and tweets

abstract: The English language has evolved dramatically throughout its lifespan, to the extent that a modern speaker of Old English would be incomprehensible without translation. One concrete indicator of this process is the movement from irregular to regular (-ed) forms for the past tense of verbs. In this study we quantify the extent of verb regularization using two vastly disparate datasets: (1) Six years of published books scanned by Google (2003--2008), and (2) A decade of social media messages posted to Twitter (2008--2017). We find that the extent of verb regularization is greater on Twitter, taken as a whole, than in English Fiction books. Regularization is also greater for tweets geotagged in the United States relative to American English books, but the opposite is true for tweets geotagged in the United Kingdom relative to British English books. We also find interesting regional variations in regularization across counties in the United States. However, once differences in population are accounted for, we do not identify strong correlations with socio-demographic variables such as education or income.


title: Lenstronomy: multi-purpose gravitational lens modelling software package

abstract: We present Lenstronomy, a multi-purpose open-source gravitational lens modeling python package. Lenstronomy is able to reconstruct the lens mass and surface brightness distributions of strong lensing systems using forward modelling. Lenstronomy supports a wide range of analytic lens and light models in arbitrary combination. The software is also able to reconstruct complex extended sources (Birrer et. al 2015) as well as being able to model point sources. We designed Lenstronomy to be stable, flexible and numerically accurate, with a clear user interface that could be deployed across different platforms. Throughout its development, we have actively used Lenstronomy to make several measurements including deriving constraints on dark matter properties in strong lenses, measuring the expansion history of the universe with time-delay cosmography, measuring cosmic shear with Einstein rings and decomposing quasar and host galaxy light. The software is distributed under the MIT license. The documentation, starter guide, example notebooks, source code and installation guidelines can be found at https://lenstronomy.readthedocs.io.


title: The Extreme Faint End of the UV Luminosity Function at $z\sim6$ Through Gravitational Telescopes: a comprehensive assessment of strong lensing uncertainties

abstract: With the Hubble Frontier Fields program, gravitational lensing has provided a powerful way to extend the study of the ultraviolet luminosity function (LF) of galaxies at $z \sim 6$ down to unprecedented magnitude limits. At the same time, significant discrepancies between different studies were found at the very faint end of the LF. In an attempt to understand such disagreements, we present a comprehensive assessment of the uncertainties associated with the lensing models and the size distribution of galaxies. We use end-to-end simulations from the source plane to the final LF that account for all lensing effects and systematic uncertainties by comparing several mass models. In addition to the size distribution, the choice of lens model leads to large differences at magnitudes fainter than $M_{UV} = -15~$ AB mag, where the magnification factor becomes highly uncertain. We perform MCMC simulations that include all these uncertainties at the individual galaxy level to compute the final LF, allowing, in particular, a crossover between magnitude bins. The best LF fit, using a modified Schechter function that allows for a turnover at faint magnitudes, gives a faint-end slope of $\alpha = -2.01_{-0.14}^{+0.12}$, a curvature parameter of $\beta = 0.48_{-0.25}^{+0.49}$, and a turnover magnitude of $M_{T} = -14.93_{-0.52}^{+0.61}$. Most importantly our procedure shows that robust constraints on the LF at magnitudes fainter than $M_{UV} = -15~$ AB remain unrealistic. More accurate lens modeling and future observations of lensing clusters with the James Webb Space Telescope can reliably extend the UV LF to fainter magnitudes.


title: Accurate nucleon electromagnetic form factors from dispersively improved chiral effective field theory

abstract: We present a theoretical parametrization of the nucleon electromagnetic form factors (FFs) based on a combination of chiral effective field theory and dispersion analysis. The isovector spectral functions on the two-pion cut are computed using elastic unitarity, chiral pion-nucleon amplitudes, and timelike pion FF data. Higher-mass isovector and isoscalar t-channel states are described by effective poles, whose strength is fixed by sum rules (charges, radii). Excellent agreement with the spacelike proton and neutron FF data is achieved up to Q^2 \sim 1 GeV^2. Our parametrization provides proper analyticity and theoretical uncertainty estimates and can be used for low-Q^2 FF studies and proton radius extraction.


title: Shedding light on the Type Ia supernova extinction puzzle: dust location found

abstract: The colour evolution of reddened Type Ia supernovae can place strong constraints on the location of dust and help address the question of whether the observed extinction stems from the interstellar medium or from circumstellar material surrounding the progenitor. Here we analyse BV photometry of 48 reddened Type Ia supernovae from the literature and estimate the dust location from their $B-V$ colour evolution. We find a time-variable colour excess $E(B-V)$ for 15 supernovae in our sample and constrain dust to distances between 0.013 and 45 pc ($4\times10^{16}-10^{20}$ cm). For the remaining supernovae, we obtain a constant $E(B-V)$ evolution and place lower limits on the dust distance from the explosion. In all the 48 supernovae, the inferred dust location is compatible with an interstellar origin for the extinction. This is corroborated by the observation that supernovae with relatively nearby dust ($\lesssim$ 1 pc) are located close to the center of their host galaxy, in high-density dusty regions where interactions between the supernova radiation and interstellar clouds close by are likely to occur. For supernovae showing time-variable $E(B-V)$, we identify a potential preference for low $R_\mathrm{V}$ values, unusually strong sodium absorption and blue-shifted and time-variable absorption features. Within the interstellar framework, this brings evidence to a proposed scenario where cloud-cloud collisions induced by the supernova radiation pressure can shift the grain size distribution to smaller values and enhance the abundance of sodium in the gaseous phase.


title: Upper limits on the presence of central massive black holes in two ultra-compact dwarf galaxies in Centaurus A

abstract: The recent discovery of massive black holes (BHs) in the centers of high-mass ultra compact dwarf galaxies (UCDs) suggests that at least some are the stripped nuclear star clusters of dwarf galaxies. We present the first study that investigates whether such massive BHs, and therefore stripped nuclei, also exist in low-mass ($M<10^{7}M_{\odot}$) UCDs. We constrain the BH masses of two UCDs located in Centaurus A (UCD320 and UCD330) using Jeans modeling of the resolved stellar kinematics from adaptive optics VLT/SINFONI data. No massive BHs are found in either UCD. We find a $3\,\sigma$ upper limit on the central BH mass in UCD\,330 of $M_{\bullet}<1.0\times10^{5}M_{\odot}$, which corresponds to 1.7\% of the total mass. This excludes a high mass fraction BH and would only allow a low-mass BHs similar to those claimed to be detected in Local Group GCs. For UCD320, poorer data quality results in a less constraining $3\,\sigma$ upper limit of $M_{\bullet}<1\times10^{6}M_{\odot}$, which is equal to 37.7\% of the total mass. The dynamical $M/L$ of UCD320 and UCD330 are not inflated compared to predictions from stellar population models. The non-detection of BHs in these low-mass UCDs is consistent with the idea that elevated dynamical $M/L$s do indicate the presence of a substantial BH. Despite not detecting massive BHs, these systems could still be stripped nuclei. The strong rotation ($v/\sigma$ of 0.3 to 0.4) in both UCDs and the two-component light profile in UCD330 support the idea that these UCDs may be stripped nuclei of low-mass galaxies where the BH occupation fraction is not yet known.


title: Higgs Couplings at High Scales

abstract: We study the off-shell production of the Higgs boson at the LHC to probe Higgs physics at higher energy scales utilizing the process $g g \rightarrow h^{*} \rightarrow ZZ$. We focus on the energy scale dependence of the off-shell Higgs propagation, and of the top quark Yukawa coupling, $y_t (Q^2)$. Extending our recent study in arXiv:1710.02149, we first discuss threshold effects in the Higgs propagator due to the existence of new states, such as a gauge singlet scalar portal, and a possible continuum of states in a conformal limit, both of which would be difficult to discover in other traditional searches. We then examine the modification of $y_t (Q^2)$ from its Standard Model (SM) prediction in terms of the renormalization group running of the top Yukawa, which could be significant in the presence of large flat extra-dimensions. Finally, we explore possible strongly coupled new physics in the top-Higgs sector that can lead to the appearance of a non-local $Q^2$-dependent form factor in the effective top-Higgs vertex. We find that considerable deviations compared to the SM prediction in the invariant mass distribution of the $Z$-boson pair can be conceivable, and may be probed at a $2\sigma$-level at the high-luminosity 14 TeV HL-LHC for a new physics scale up to $\mathcal{O}(1 {~\rm TeV})$, and at the upgraded 27 TeV HE-LHC for a scale up to $\mathcal{O}(3 {~\rm TeV})$. For a few favorable scenarios, $5\sigma$-level observation may be possible at the HE-LHC for a scale of about $\mathcal{O}(1 {~\rm TeV})$.


title: Geometric mean of bimetric spacetimes

abstract: We use the geometric mean to parametrize metrics in the Hassan-Rosen ghost-free bimetric theory and pose the initial-value problem. The geometric mean of two positive definite symmetric matrices is a well-established mathematical notion which can be, under certain conditions, extended to quadratic forms having the Lorentzian signature, say metrics $g$ and $f$. In such a case, the null cone of the geometric mean metric $h$ is in the middle of the null cones of $g$ and $f$ appearing as a geometric average of a bimetric spacetime. The parametrization based on $h$ ensures the reality of the square root in the ghost-free bimetric interaction potential. Subsequently, we derive the standard $n+1$ decomposition in a frame adapted to the geometric mean and state the initial-value problem, that is, the evolution equations, the constraints, and the preservation of the constraints equation.


title: The holographic interpretation of $J \bar T$-deformed CFTs

abstract: Recently, a non-local yet possibly UV-complete quantum field theory has been constructed by deforming a two-dimensional CFT by the composite operator $J \bar T$, where $J$ is a chiral $U(1)$ current and $\bar T$ is a component of the stress tensor. Assuming the original CFT was a holographic CFT, we work out the holographic dual of its $J \bar T$ deformation. We find that the dual spacetime is still AdS$_3$, but with modified boundary conditions that mix the metric and the Chern-Simons gauge field dual to the $U(1)$ current. We show that when the coefficient of the chiral anomaly for $J$ vanishes, the energy and thermodynamics of black holes obeying these modified boundary conditions precisely reproduce the previously derived field theory spectrum and thermodynamics. Our proposed holographic dictionary can also reproduce the field-theoretical spectrum in presence of the chiral anomaly, upon a certain assumption that we justify. The asymptotic symmetry group associated to these boundary conditions consists of two copies of the Virasoro and one copy of the $U(1)$ Ka\v{c}-Moody algebra, just as before the deformation; the only effect of the latter is to modify the spacetime dependence of the right-moving Virasoro generators, whose action becomes state-dependent and effectively non-local.


title: Properties of thermal quantum states: locality of temperature, decay of correlations, and more

abstract: We review several properties of thermal states of spin Hamiltonians with short range interactions. In particular, we focus on those aspects in which the application of tools coming from quantum information theory has been specially successful in the recent years. This comprises the study of the correlations at finite and zero temperature, the stability against distant and/or weak perturbations, the locality of temperature and their classical simulatability. For the case of states with a finite correlation length, we overview the results on their energy distribution and the equivalence of the canonical and microcanonical ensemble.

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


title: Full Counting Statistics in the Transverse Field Ising Chain

abstract: We consider the full probability distribution for the transverse magnetization of a finite subsystem in the transverse field Ising chain. We derive a determinant representation of the corresponding characteristic function for general Gaussian states. We consider applications to the full counting statistics in the ground state, finite temperature equilibrium states, non-equilibrium steady states and time evolution after global quantum quenches. We derive an analytical expression for the time and subsystem size dependence of the characteristic function at sufficiently late times after a quantum quench. This expression features an interesting multiple light-cone structure.


title: Non-power-law universality in one-dimensional quasicrystals

abstract: We have investigated scaling properties of the Aubry-Andr\'e model and related one-dimensional quasiperiodic Hamiltonians near their localisation transitions. We find numerically that the scaling of characteristic energies near the ground state, usually captured by a single dynamical exponent, does not obey a power law relation. Instead, the scaling behaviour depends strongly on the correlation length in a manner governed by the continued fraction expansion of the irrational number $\beta$ describing incommensurability in the system. This dependence is, however, found to be universal between a range of models sharing the same value of $\beta$. For the Aubry-Andr\'e model, we explain this behaviour in terms of a discrete renormalisation group protocol which predicts rich critical behaviour. This result is complemented by studies of the expansion dynamics of a wave packet under the Aubry-Andr\'e model at the critical point. Anomalous diffusion exponents are derived in terms of multifractal (R\'enyi) dimensions of the critical spectrum; non-power-law universality similar to that found in ground state dynamics is observed between a range of critical tight-binding Hamiltonians.

categories: cond-mat.dis-nn cond-mat.quant-gas cond-mat.stat-mech physics.atom-ph quant-ph


title: Measurement of ambient neutrons in an underground laboratory at Kamioka Observatory

abstract: Ambient neutrons are one of the most serious backgrounds for underground experiments searching for rare events. The ambient neutron flux in an underground laboratory of Kamioka Observatory was measured using a $\mathrm{^3He}$ proportional counter with various moderator setups. Since the detector response largely depends on the spectral shape, the energy spectra of the neutrons transported from the rock to the laboratory are estimated by Monte-Carlo simulations. The ratio of the thermal neutron flux to the total neutron flux was found to depend on the thermalizing efficiency of the rock. Therefore, the ratio of the count rate without a moderator to that with a moderator was used to determine this parameter. Consequently, the most-likely neutron spectrum predicted by the simulations for the parameters determined by the experimental results was obtained. The result suggests an interesting spectral shape, which has not been indicated in previous studies. The total ambient neutron flux is $(23.52 \pm 0.68 \ \mathrm{_{stat.}} ^{+1.87}_{-2.13} \ \mathrm{_{sys.}}) \times 10^{-6}$ cm$^{-2}$ s$^{-1}$. This result, especially the energy spectrum information, could be a new and important input for estimating the background in current and future experiments in the underground laboratory at Kamioka Observatory.


title: Flag fault-tolerant error correction, measurement, and quantum computation for cyclic CSS codes

abstract: Flag qubits have recently been proposed in syndrome extraction circuits to detect high-weight errors arising from fewer faults. The use of flag qubits allows the construction of fault-tolerant protocols with the fewest number of ancillas known to-date. In this work, we prove some critical properties of CSS codes constructed from classical cyclic codes that enable the construction of a flag fault-tolerant error correction scheme. We then develop fault-tolerant protocols as well as a family of circuits for flag fault-tolerant error correction and operator measurement, requiring only four ancilla qubits and applicable to cyclic CSS codes of distance 3. The measurement protocol can be further used for logical Clifford gate implementation via quantum gate teleportation. We also provide examples of cyclic CSS codes with large encoding rates.


title: The chemical composition of NGC5824, a globular cluster without iron spread but with an extreme Mg-Al anticorrelation

abstract: NGC5824 is a massive Galactic globular cluster suspected to have an intrinsic spread in its iron content, according to the strength of the calcium triplet lines. We present chemical abundances of 117 cluster giant stars using high-resolution spectra acquired with the multi-object spectrograph FLAMES. The metallicity distribution of 87 red giant branch stars is peaked at [Fe/H]=-2.11+-0.01 dex, while that derived from 30 asymptotic giant branch stars is peaked at [Fe/H]=-2.20+-0.01 dex. Both the distributions are compatible with a null spread, pointing out that this cluster did not retain the ejecta of supernovae. The small iron abundance offset between the two groups of stars is similar to those already observed among red and asymptotic giant branch stars in other clusters. The lack of intrinsic iron spread rules out the possibility that NGC5824 is the remnant of a disrupted dwarf galaxy, as previously suggested. We also find evidence of the chemical anomalies usually observed in globular clusters, namely the Na-O and the Mg-Al anticorrelations. In particular, NGC5824 exhibits a huge range of [Mg/Fe] abundance, observed only in a few metal-poor and/or massive clusters. We conclude that NGC5824 is a normal globular cluster, without spread in [Fe/H] but with an unusually large spread in [Mg/Fe], possibly due to an efficient self-enrichment driven by massive asymptotic giant branch stars.


title: Predicting the Future with Transformational States

abstract: An intelligent observer looks at the world and sees not only what is, but what is moving and what can be moved. In other words, the observer sees how the present state of the world can transform in the future. We propose a model that predicts future images by learning to represent the present state and its transformation given only a sequence of images. To do so, we introduce an architecture with a latent state composed of two components designed to capture (i) the present image state and (ii) the transformation between present and future states, respectively. We couple this latent state with a recurrent neural network (RNN) core that predicts future frames by transforming past states into future states by applying the accumulated state transformation with a learned operator. We describe how this model can be integrated into an encoder-decoder convolutional neural network (CNN) architecture that uses weighted residual connections to integrate representations of the past with representations of the future. Qualitatively, our approach generates image sequences that are stable and capture realistic motion over multiple predicted frames, without requiring adversarial training. Quantitatively, our method achieves prediction results comparable to state-of-the-art results on standard image prediction benchmarks (Moving MNIST, KTH, and UCF101).


title: Semi-analytic forecasts for JWST - I. UV luminosity functions at z = 4 - 10

abstract: In anticipation of the upcoming deployment of the James Webb Space Telescope (JWST), we present high-redshift predictions by the well-established Santa Cruz semi-analytic model. We update the models by re-calibrating them after adopting cosmological parameters consistent with recent constraints from Planck. We provide predictions for rest-frame UV luminosity functions for galaxy populations over a wide range of $M_\text{UV}$ from $\sim-6$ to $\sim-24$ between $z = 4 - 10$. In addition, we present the corresponding predictions for observed-frame galaxy number counts in different redshift bins in the full set of NIRCam filters. We provide predictions of the quantitative effect on these observables of varying the physical recipes implemented in the models, such as the molecular gas depletion time (star formation efficiency) scalings or the scalings of outflow rates driven by stars and supernovae with galaxy circular velocity. Based on these results, we discuss what may be learned about the physical processes that shape galaxy formation from JWST observations of galaxy number densities at different intrinsic luminosities. All data tables for the results presented in this work are available at https://www.simonsfoundation.org/semi-analytic-forecasts-for-jwst/.


title: Core formation from self-heating dark matter

abstract: Cosmological simulations of the $\Lambda$CDM model suggest that the dark matter halos of dwarf galaxies are denser in their center than what observational data of such galaxies imply. In this letter, we propose a novel solution to this problem by invoking a certain class of dark matter self-heating processes. As we will argue, such processes lead to the formation of dark matter cores at late times by considerably reducing the inner mass density of dwarf-sized halos. For deriving concrete results, we focus on semi-annihilating dark matter scenarios and model the inner region of dark matter halos as a gravothermal fluid. An important aspect of this new solution is that the semi-annihilation effects are much more prominent in dwarf-sized halos than in the more massive halos that host galaxies and clusters, even if the corresponding cross sections are the same. Furthermore, the preferred parameter space for solving the small-scale problem suggests a thermal dark matter candidate with a mass below the GeV scale, which can be probed in dark matter direct and indirect detection experiments.


title: A New Test of Copper and Zinc Abundances in Late-Type Stars Using Ultraviolet Cu II and Zn II Lines

abstract: We present new abundances derived from Cu I, Cu II, Zn I, and Zn II lines in six warm (5766 < Teff < 6427 K), metal-poor (-2.50 < [Fe/H] < -0.95) dwarf and subgiant (3.64 < log g < 4.44) stars. These abundances are derived from archival high-resolution ultraviolet spectra from the Space Telescope Imaging Spectrograph on board the Hubble Space Telescope and ground-based optical spectra from several observatories. Ionized Cu and Zn are the majority species, and abundances derived from Cu II and Zn II lines should be largely insensitive to departures from local thermodynamic equilibrium (LTE). We find good agreement between the [Zn/H] ratios derived separately from Zn I and Zn II lines, suggesting that departures from LTE are, at most, minimal (< 0.1 dex). We find that the [Cu/H] ratios derived from Cu II lines are 0.36 +/- 0.06 dex larger than those derived from Cu I lines in the most metal-poor stars ([Fe/H] < -1.8), suggesting that LTE underestimates the Cu abundance derived from Cu I lines. The deviations decrease in more metal-rich stars. Our results validate previous theoretical non-LTE calculations for both Cu and Zn, supporting earlier conclusions that the enhancement of [Zn/Fe] in metal-poor stars is legitimate, and the deficiency of [Cu/Fe] in metal-poor stars may not be as large as previously thought.


title: Dimensional regularization vs methods in fixed dimension with and without $\gamma_5$

abstract: We study the Lorentz and Dirac algebra, including antisymmetric $\epsilon$ tensors and the $\gamma_5$ matrix, in implicit gauge-invariant regularization/renormalization methods defined in fixed integer dimensions. They include constrained differential, implicit and four-dimensional renormalization. We find that these fixed-dimension methods face the same difficulties as the different versions of dimensional regularization. We propose a consistent procedure in these methods, similar to the consistent version of regularization by dimensional reduction.


title: Discovery of a diffuse optical line emitting halo in the core of the Centaurus cluster of galaxies: Line emission outside the protection of the filaments

abstract: We present the discovery of diffuse optical line emission in the Centaurus cluster seen with the MUSE IFU. The unparalleled sensitivity of MUSE allows us to detect the faint emission from these structures which extend well beyond the bounds of the previously known filaments. Diffuse structures (emission surrounding the filaments, a northern shell and an extended Halo) are detected in many lines typical of the nebulae in cluster cores ([NII]$_{\lambda 6548\&6583}$ ,[SII]$_{\lambda 6716\&6731}$, [OI]$_{\lambda 6300}$, [OIII]$_{\lambda 4959\&5007}$ etc.) but are more than an order of magnitude fainter than the filaments, with the faint halo only detected through the brightest line in the spectrum ([NII]$_{\lambda 6583}$). These structures are shown to be kinematically distinct from the stars in the central galaxy and have different physical and excitation states to the filaments. Possible origins are discussed for each structure in turn and we conclude that shocks and/or pressure imbalances are resulting in gas dispersed throughout the cluster core, formed from either disrupted filaments or direct cooling, which is not confined to the bright filaments.