title: Modulation theory and resonant regimes for dispersive shock waves in nematic liquid crystals

abstract: A full analysis of all regimes for optical dispersive shock wave (DSW) propagation in nematic liquid crystals is undertaken. These dispersive shock waves are generated from step initial conditions for the optical field and are resonant in that linear diffractive waves are in resonance with the DSW, resulting in a resonant linear wavetrain propagating ahead of it. It is found that there are six regimes, which are distinct and require different solution methods. In previous studies, the same solution method was used for all regimes, which does not yield solutions in full agreement with numerical solutions. Indeed, the standard DSW structure disappears for sufficiently large initial jumps. Asymptotic theory, approximate methods or Whitham modulation theory are used to find solutions for these resonant dispersive shock waves in a given regime. The solutions are found to be in excellent agreement with numerical solutions of the nematic equations in all regimes. It is found that for small initial jumps, the resonant wavetrain is unstable, but that it stabilises above a critical jump height. It is additionally found that the DSW is unstable, except for small jump heights for which there is no resonance and large jump heights for which there is no standard DSW structure.


title: When black box algorithms are (not) appropriate: a principled prediction-problem ontology

abstract: In the 1980s a new, extraordinarily productive way of reasoning about algorithms emerged. In this paper, we introduce the term "outcome reasoning" to refer to this form of reasoning. Though outcome reasoning has come to dominate areas of data science, it has been under-discussed and its impact under-appreciated. For example, outcome reasoning is the primary way we reason about whether ``black box'' algorithms are performing well. In this paper we analyze outcome reasoning's most common form (i.e., as "the common task framework") and its limitations. We discuss why a large class of prediction-problems are inappropriate for outcome reasoning. As an example, we find the common task framework does not provide a foundation for the deployment of an algorithm in a real world situation. Building off of its core features, we identify a class of problems where this new form of reasoning can be used in deployment. We purposefully develop a novel framework so both technical and non-technical people can discuss and identify key features of their prediction problem and whether or not it is suitable for outcome reasoning.


title: Integrating data science ethics into an undergraduate major: A case study

abstract: We present a programmatic approach to incorporating ethics into an undergraduate major in statistical and data sciences. We discuss departmental-level initiatives designed to meet the National Academy of Sciences recommendation for integrating ethics into the curriculum from top-to-bottom as our majors progress from our introductory courses to our senior capstone course, as well as from side-to-side through co-curricular programming. We also provide six examples of data science ethics modules used in five different courses at our liberal arts college, each focusing on a different ethical consideration. The modules are designed to be portable such that they can be flexibly incorporated into existing courses at different levels of instruction with minimal disruption to syllabi. We connect our efforts to a growing body of literature on the teaching of data science ethics, present assessments of our effectiveness, and conclude with next steps and final thoughts.


title: Dynamics of a stochastic excitable system with slowly adapting feedback

abstract: We study an excitable active rotator with slowly adapting nonlinear feedback and noise. Depending on the adaptation and the noise level, this system may display noise-induced spiking, noise-perturbed oscillations, or stochastic busting. We show how the system exhibits transitions between these dynamical regimes, as well as how one can enhance or suppress the coherence resonance, or effectively control the features of the stochastic bursting. The setup can be considered as a paradigmatic model for a neuron with a slow recovery variable or, more generally, as an excitable system under the influence of a nonlinear control mechanism. We employ a multiple timescale approach that combines the classical adiabatic elimination with averaging of rapid oscillations and stochastic averaging of noise-induced fluctuations by a corresponding stationary Fokker-Planck equation. This allows us to perform a numerical bifurcation analysis of a reduced slow system and to determine the parameter regions associated with different types of dynamics. In particular, we demonstrate the existence of a region of bistability, where the noise-induced switching between a stationary and an oscillatory regime gives rise to stochastic bursting.


title: Progenitors of early-time interacting supernovae

abstract: We compute an extensive set of early-time spectra of supernovae interacting with circumstellar material using the radiative transfer code CMFGEN. Our models are applicable to events observed from 1 to a few days after explosion. Using these models, we constrain the progenitor and explosion properties of a sample of 17 observed interacting supernovae at early-times. Because massive stars have strong mass loss, these spectra provide valuable information about supernova progenitors, such as mass-loss rates, wind velocities, and surface abundances. We show that these events span a wide range of explosion and progenitor properties, exhibiting supernova luminosities in the 1e8 to 1e12 Lsun range, temperatures from 10000 to 60000 K, progenitor mass-loss rates from a few 1e-4 up to 1 Msun/yr, wind velocities from 100 to 800 km/s, and surface abundances from solar-like to H-depleted. Our results suggest that many progenitors of supernovae interacting with circumstellar material have significantly increased mass-loss before explosion compared to what massive stars show during the rest of their lifetimes. We also infer a lack of correlation between surface abundances and mass-loss rates. This may point to the pre-explosion mass-loss mechanism being independent of stellar mass. We find that the majority of these events have CNO-processed surface abundances. In the single star scenario this points to a preference towards high-mass RSGs as progenitors of interacting SNe, while binary evolution could impact this conclusion. Our models are publicly available and readily applicable to analyze results from ongoing and future large scale surveys such as the Zwicky Transient Factory.


title: A new method for constructing squeezed states for the isotropic 2D harmonic oscillator

abstract: We introduce a new method for constructing squeezed states for the 2D isotropic harmonic oscillator. Based on the construction of coherent states in [1], we define a new set of ladder operators for the 2D system as a linear combination of the x and y ladder operators and construct the SU(2) coherent states. The new ladder operators are used for generalizing the squeezing operator to 2D and the SU(2) coherent states play the role of the Fock states in the expansion of the 2D squeezed states. We discuss some properties of the 2D squeezed states.


title: Many-body localization of zero modes

abstract: The celebrated Dyson singularity signals the relative delocalization of single-particle wave functions at the zero-energy symmetry point of disordered systems with a chiral symmetry. Here we show that analogous zero modes in interacting quantum systems can fully localize at sufficiently large disorder, but do so less strongly than nonzero modes, as signifed by their real-space and Fock-space localization characteristics. We demonstrate this effect in a spin-1 Ising chain, which naturally provides a chiral symmetry in an odd-dimensional Hilbert space, thereby guaranteeing the existence of a many-body zero mode at all disorder strengths. In the localized phase, the bipartite entanglement entropy of the zero mode follows an area law, but is enhanced by a system-size-independent factor of order unity when compared to the nonzero modes. Analytically, this feature can be attributed to a specific zero-mode hybridization pattern on neighboring spins. The zero mode also displays a symmetry-induced even-odd and spin-orientation fragmentation of excitations, characterized by real-space spin correlation functions, which generalizes the sublattice polarization of topological zero modes in noninteracting systems, and holds at any disorder strength.


title: A Sub-grid Scale Energy Dissipation Rate Model for Large-eddy Spray Simulations

abstract: In high Reynolds number turbulent flows, energy dissipation refers to the process of energy transfer from kinetic energy to internal energy due to molecular viscosity. In large eddy simulation (LES) with one-equation turbulence models, the energy dissipation process is modeled by a rate term in the transport equation of the subgrid-scale (SGS) kinetic energy. Despite its important role in maintaining a proper energy balance between the resolved and SGS scales, modeling of the energy dissipation rate has received scarce attention. In this paper, a SGS model belonging to the dynamic structure family is developed based on findings from direct numerical simulation (DNS) studies of decaying isotropic turbulence. The model utilizes a Leonard-type term, a SGS viscosity, and a characteristic scaling term to predict the energy dissipation rate in LES. A posteriori tests of the model have been carried out under direct-injection gasoline and diesel engine-like conditions. Spray characteristics such as penetration rates and mixture fractions have been examined. It is found that the current SGS model accurately predicts vapor-phase penetrations across different mesh resolutions under both gasoline and diesel spray conditions, due to its correct scaling of SGS energy dissipation rate with the SGS kinetic energy and LES fitter width. In contrast, the classic model that is widely used in the literature predicts a scaling of energy dissipation rate upon mesh resolution, exhibiting a noticeable mesh dependence.


title: Coherence via Well-Foundedness: Taming Set-Quotients in Homotopy Type Theory

abstract: Suppose we are given a graph and want to show a property for all its cycles (closed chains). Induction on the length of cycles does not work since sub-chains of a cycle are not necessarily closed. This paper derives a principle reminiscent of induction for cycles for the case that the graph is given as the symmetric closure of a locally confluent and (co-)well-founded relation. We show that, assuming the property in question is sufficiently nice, it is enough to prove it for the empty cycle and for cycles given by local confluence. Our motivation and application is in the field of homotopy type theory, which allows us to work with the higher-dimensional structures that appear in homotopy theory and in higher category theory, making coherence a central issue. This is in particular true for quotienting - a natural operation which gives a new type for any binary relation on a type and, in order to be well-behaved, cuts off higher structure (set-truncates). The latter makes it hard to characterise the type of maps from a quotient into a higher type, and several open problems stem from this difficulty. We prove our theorem on cycles in a type-theoretic setting and use it to show coherence conditions necessary to eliminate from set-quotients into 1-types, deriving approximations to open problems on free groups and pushouts. We have formalised the main result in the proof assistant Lean.


title: Proof of the Peres conjecture for contextuality

abstract: A central result in the foundations of quantum mechanics is the Kochen-Specker theorem. In short, it states that quantum mechanics cannot be reconciled with classical models that are noncontextual for ideal measurements. The first explicit derivation by Kochen and Specker was rather complex, but considerable simplifications have been achieved thereafter. We propose a systematic approach to find minimal Hardy-type and Greenberger-Horne-Zeilinger-type (GHZ-type) proofs of the Kochen-Specker theorem, these are characterized by the fact that the predictions of classical models are opposite to the predictions of quantum mechanics. Based on our results, we show that the Kochen-Specker set with 18 vectors from Cabello et al. [A. Cabello et al., Phys. Lett. A 212, 183 (1996)] is the minimal set for any dimension, verifying a longstanding conjecture by Peres. Our results allow to identify minimal contextuality scenarios and to study their usefulness for information processing.


title: A Closed-form Solution for the Kinematics of Asymmetric Miura Vertices

abstract: The Miura vertex is a versatile origami pattern found in a variety of mechanisms. Previous papers have derived and validated a closed-form solution for the kinematics of a symmetric Miura vertex, but the motion of an asymmetric vertex has only been shown numerically. In this paper, we present the trigonometric derivation of a closed-form solution for the folding of an asymmetric Miura vertex.


title: A symmetric protocol to establish service level agreements

abstract: We present a symmetrical protocol to repeatedly negotiate a desired service level between two parties, where the service levels are taken from some totally ordered finite domain. The agreed service level is selected from levels dynamically proposed by both parties and parties can only decrease the desired service level during a negotiation. The correctness of the protocol is stated using modal formulas and its behaviour is explained using behavioural reductions of the external behaviour modulo weak trace equivalence and divergence-preserving branching bisimulation. Our protocol originates from an industrial use case and it turned out to be remarkably tricky to design correctly.


title: Discovery of a Sirius-like binary system with a very strongly magnetic white dwarf. Binarity among magnetic white dwarfs

abstract: Binary systems containing a magnetic white dwarf and a main-sequence star are considered extremely rare, perhaps non-existent. In the course of a search of magnetic fields in high-mass white dwarfs we have discovered a Sirius-like wide binary system composed of a main-sequence G0 star and a $M\sim 1.1\,M_\odot$ white dwarf with a huge (hundreds of MG) magnetic field. This star, WDS J03038+0608B, shows a circular polarisation amplitude of 5% in the continuum, with no evidence of variability in a 1d time-scale, little or no linear polarisation in the blue part of the spectrum, and about 2% linear polarisation in the red part of the optical spectrum. A search in the literature reveals the existence of at least four more binary systems that include a magnetic white dwarf and a non-degenerate companion - three such systems passed unremarked in previous studies. We estimate that up to a few percent of magnetic white dwarfs may be found to occur in wide binary pairs. However, at least four of the five known binary systems with a magnetic white dwarf are too widely separated to be expected to evolve into systems experiencing Roche lobe overflow, and cannot be considered as progenitors of magnetic cataclysmic variable (AM Her and DQ Her) systems.


title: On the Confidence in Bit-Alias Measurement of Physical Unclonable Functions

abstract: Physical Unclonable Functions (PUFs) are modern solutions for cheap and secure key storage. The security level strongly depends on a PUF's unpredictability, which is impaired if certain bits of the PUF response tend towards the same value on all devices. The expectation for the probability of 1 at some position in the response, the Bit-Alias, is a state-of-the-art metric in this regard. However, the confidence interval of the Bit-Alias is never considered, which can lead to an overestimation of a PUF's unpredictability. Moreover, no tool is available to verify if the Bit-Alias is within given limits. This work adapts a method for the calculation of confidence intervals to Bit-Alias. It further proposes a statistical hypothesis test to verify if a PUF design meets given specifications on Bit-Alias or bit-wise entropy. Application to several published PUF designs demonstrates the methods' capabilities. The results prove the need for a high number of samples when the unpredictability of PUFs is tested. The proposed methods are publicly available and should improve the design and evaluation of PUFs in the future.


title: Central limit theorem for mesoscopic eigenvalue statistics of the free sum of matrices

abstract: We consider random matrices of the form $H_N=A_N+U_N B_N U^*_N$, where $A_N$, $B_N$ are two $N$ by $N$ deterministic Hermitian matrices and $U_N$ is a Haar distributed random unitary matrix. We establish a universal Central Limit Theorem for the linear eigenvalue statistics of $H_N$ on all mesoscopic scales inside the regular bulk of the spectrum. The proof is based on studying the characteristic function of the linear eigenvalue statistics, and consists of two main steps: (1) generating Ward identities using the left-translation-invariance of the Haar measure, along with a local law for the resolvent of $H_N$ and analytic subordination properties of the free additive convolution, allow us to derive an explicit formula for the derivative of the characteristic function; (2) a local law for two-point product functions of resolvents is derived using a partial randomness decomposition of the Haar measure. We also prove the corresponding results for orthogonal conjugations.


title: Finite Deformations of Quantum Mechanics

abstract: We investigate modifications of quantum mechanics (QM) that replace the unitary group in a finite dimensional Hilbert space with a finite group and determine the minimal sequence of subgroups necessary to approximate QM arbitrarily closely for general choices of Hamiltonian. This mathematical study reveals novel insights about 't Hooft's Ontological Quantum Mechanics, and the derivation of statistical mechanics from quantum mechanics. We show that Kornyak's proposal to understand QM as classical dynamics on a Hilbert space of one dimension higher than that describing the universe, supplemented by a choice of the value of a naturally conserved quantum operator in that classical evolution can probably be a model of the world we observe.


title: X-ray spectra and polarization from magnetar candidates

abstract: Magnetars are believed to host the strongest magnetic fields in the present universe ($B\gtrsim10^{14}$ G) and the study of their persistent emission in the X-ray band offers an unprecendented opportunity to gain insight into physical processes in the presence of ultra-strong magnetic fields. Up to now, most of our knowledge about magnetar sources came from spectral analysis, which allowed to test the resonant Compton scattering scenario and to probe the structure of the star magnetosphere. On the other hand, radiation emitted from magnetar surface is expected to be strongly polarized and its observed polarization pattern bears the imprint of both scatterings onto magnetospheric charges and QED effects as it propagates in the magnetized vacuum around the star. X-ray polarimeters scheduled to fly in the next years will finally allow to exploit the wealth of information stored in the polarization observables. Here we revisit the problem of assessing the spectro-polarimetric properties of magnetar persistent emission. At variance with previous investigations, proper account for more physical surface emission models is made by considering either a condensed surface or a magnetized atmosphere. Results are used to simulate polarimetric observations with the forthcoming Imaging X-ray Polarimetry Explorer (IXPE). We find that X-ray polarimetry will allow to detect QED vacuum effects for all the emission models we considered and to discriminate among them.


title: Regulation of a single-server queue with customers who dynamically choose their service durations

abstract: In recent years, there is a growing research about queueing models with customers who choose their service durations. In general, the model assumptions in the existing literature imply that every customer knows his service demand when he enters into the service position. Clearly, this property is not consistent with some real-life situations. Thus, motivated by this issue, this work includes a single-server queueing model with customers who dynamically choose their service durations. In this setup, it is shown how to derive a price-regulation (Pigouvian tax) which implies an optimal resource allocation from a social point of view. In particular, this work includes an efficient procedure for computation of the parameters of this price function.


title: Extracting high-level information from gamma-ray burst supernova spectra

abstract: Radiation transport codes are often used in astrophysics to construct spectral models. In this work we demonstrate how producing these models for a time series of data can provide unique information about supernovae (SNe). Unlike previous work, we specifically concentrate on the method for obtaining the best synthetic spectral fits, and the errors associated with the preferred model parameters. We demonstrate how varying the ejecta mass, bolometric luminosity ($L_{bol}$) and photospheric velocity ($v_{ph}$), affects the outcome of the synthetic spectra. As an example we analyze the photospheric phase spectra of the GRB-SN\,2016jca. It is found that for most epochs (where the afterglow subtraction is small) the error on $L_{bol}$ and $v_{ph}$ was $\sim$5\%. The uncertainty on ejecta mass and K.E. was found to be $\sim$20%, although this can be expected to dramatically decrease if models of nebular phase data can be simultaneously produced. We also demonstrate how varying the elemental abundance in the ejecta can produce better synthetic spectral fits. In the case of SN\,2016jca it is found that a decreasing $^{56}$Ni abundance as a function of decreasing velocity produces the best fit models. This could be the case if the $^{56}$Ni was sythesised at the side of the GRB jet, or dredged up from the centre of the explosion. The work presented here can be used as a guideline for future studies on supernovae which use the same or similar radiation transfer code.


title: Magnetic field independent sub-gap states in hybrid Rashba nanowires

abstract: Sub-gap states in semiconducting-superconducting nanowire hybrid devices are controversially discussed as potential topologically non-trivial quantum states. One source of ambiguity is the lack of an energetically and spatially well defined tunnel spectrometer. Here, we use quantum dots directly integrated into the nanowire during the growth process to perform tunnel spectroscopy of discrete sub-gap states in a long nanowire segment. In addition to sub-gap states with a standard magnetic field dependence, we find topologically trivial sub-gap states that are independent of the external magnetic field, i.e. that are pinned to a constant energy as a function of field. We explain this effect qualitatively and quantitatively by taking into account the strong spin-orbit interaction in the nanowire, which can lead to a decoupling of Andreev bound states from the field due to a spatial spin texture of the confined eigenstates.


title: Confirmation of WASP-107b's extended Helium atmosphere with Keck II/NIRSPEC

abstract: We present the detection of helium in the extended atmosphere of the sub-Saturn WASP-107b using high resolution ($R \approx 25000$) near-infrared spectra from Keck II/NIRSPEC. We find peak excess absorption of $7.26 \pm 0.24\%$ (30$\sigma$) centered on the HeI triplet at 10833A. The amplitude and shape of the helium absorption profile is in excellent agreement with previous observations of escaping helium from this planet made by CARMENES and HST. This suggests there is no significant temporal variation in the signature of escaping helium from the planet over a two year baseline. This result demonstrates Keck II/NIRSPEC's ability to detect atmospheric escape in exoplanets, making it a useful instrument to further our understanding of the evaporation of exoplanetary atmospheres via ground-based observations of HeI.


title: Data-Driven Approach for Uncertainty Propagation and Reachability Analysis in Dynamical Systems

abstract: In this paper, we propose a data-driven approach for uncertainty propagation and reachability analysis in a dynamical system. The proposed approach relies on the linear lifting of a nonlinear system using linear Perron-Frobenius (P-F) and Koopman operators. The uncertainty can be characterized in terms of the moments of a probability density function. We demonstrate how the P-F and Koopman operators are used for propagating the moments. Time-series data is used for the finite-dimensional approximation of the linear operators, thereby enabling data-driven approach for moment propagation. Simulation results are presented to demonstrate the effectiveness of the proposed method.


title: Globalization for geometric partial comodules

abstract: We discuss globalization for geometric partial comodules in a monoidal category with pushouts and we provide a concrete procedure to construct it, whenever it exists. The mild assumptions required by our approach make it possible to apply it in a number of contexts of interests, recovering and extending numerous ad hoc globalization constructions from the literature in some cases and providing obstruction for globalization in some other cases.


title: LOcAl DEcisions on Replicated States (LOADER) in programmable data planes: programming abstraction and experimental evaluation

abstract: Programmable data planes recently emerged as a prominent innovation in Software Defined Networking (SDN), by permitting support of stateful flow processing functions over hardware network switches specifically designed for network processing. Unlike early SDN solutions such as OpenFlow, modern stateful data planes permit to keep (and dynamically update) local per-flow states inside network switches, thus dramatically improving reactiveness of network applications to state changes. Still, also in stateful data planes, the control and update of non-local states is assumed to be completely delegated to a centralized controller and thus accessed only at the price of extra delay. Our LOADER proposal aims at contrasting the apparent dichotomy between local states and global states. We do so by introducing a new possibility: permit to take localized (in-switch) decisions not only on local states but also on replicated global states, thus providing support for network-wide applications without incurring the drawbacks of classical approaches. To this purpose, i) we provide high-level programming abstractions devised to define the states and the update logic of a generic network-wide application, and ii) we detail the underlying low level state management and replication mechanisms. We then show LOADER's independence of the stateful data plane technology employed, by implementing it over two distinct stateful data planes (P4 switches and OPP - Open Packet Processor - switches), and by experimentally validating both implementations in an emulated testbed using a simple distributed Deny-of-Service (DoS) detection application.


title: Efficient orthogonal control of tunnel couplings in a quantum dot array

abstract: Electrostatically-defined semiconductor quantum dot arrays offer a promising platform for quantum computation and quantum simulation. However, crosstalk of gate voltages to dot potentials and inter-dot tunnel couplings complicates the tuning of the device parameters. To date, crosstalk to the dot potentials is routinely and efficiently compensated using so-called virtual gates, which are specific linear combinations of physical gate voltages. However, due to exponential dependence of tunnel couplings on gate voltages, crosstalk to the tunnel barriers is currently compensated through a slow iterative process. In this work, we show that the crosstalk on tunnel barriers can be efficiently characterized and compensated for, using the fact that the same exponential dependence applies to all gates. We demonstrate efficient calibration of crosstalk in a quadruple quantum dot array and define a set of virtual barrier gates, with which we show orthogonal control of all inter-dot tunnel couplings. Our method marks a key step forward in the scalability of the tuning process of large-scale quantum dot arrays.


title: Streaming Complexity of Spanning Tree Computation

abstract: The semi-streaming model is a variant of the streaming model frequently used for the computation of graph problems. It allows the edges of an $n$-node input graph to be read sequentially in $p$ passes using $\tilde{O}(n)$ space. In this model, some graph problems, such as spanning trees and $k$-connectivity, can be exactly solved in a single pass; while other graph problems, such as triangle detection and unweighted all-pairs shortest paths, are known to require $\tilde{\Omega}(n)$ passes to compute. For many fundamental graph problems, the tractability in these models is open. In this paper, we study the tractability of computing some standard spanning trees. Our results are: (1) Maximum-Leaf Spanning Trees. This problem is known to be APX-complete with inapproximability constant $\rho\in[245/244,2)$. By constructing an $\varepsilon$-MLST sparsifier, we show that for every constant $\varepsilon > 0$, MLST can be approximated in a single pass to within a factor of $1+\varepsilon$ w.h.p. (albeit in super-polynomial time for $\varepsilon \le \rho-1$ assuming $\mathrm{P} \ne \mathrm{NP}$). (2) BFS Trees. It is known that BFS trees require $\omega(1)$ passes to compute, but the na\"{i}ve approach needs $O(n)$ passes. We devise a new randomized algorithm that reduces the pass complexity to $O(\sqrt{n})$, and it offers a smooth tradeoff between pass complexity and space usage. (3) DFS Trees. The current best algorithm by Khan and Mehta {[}STACS 2019{]} takes $\tilde{O}(h)$ passes, where $h$ is the height of computed DFS trees. Our contribution is twofold. First, we provide a simple alternative proof of this result, via a new connection to sparse certificates for $k$-node-connectivity. Second, we present a randomized algorithm that reduces the pass complexity to $O(\sqrt{n})$, and it also offers a smooth tradeoff between pass complexity and space usage.


title: An inverse problem for Moore-Gibson-Thompson equation arising in high intensity ultrasound

abstract: In this article we study the inverse problem of recovering a space-dependent coefficient of the Moore-Gibson-Thompson (MGT) equation, from knowledge of the trace of the solution on some open subset of the boundary. We obtain the Lipschitz stability for this inverse problem, and we design a convergent algorithm for the reconstruction of the unknown coefficient. The techniques used are based on Carleman inequalities for wave equations and properties of the MGT equation.


title: Huge linear magnetoresistance due to open orbits in $\gamma$-PtBi$_2$

abstract: Some single-crystalline materials present an electrical resistivity which decreases between room temperature and low temperatures at zero magnetic field as in a good metal and switches to a nearly semiconductinglike behavior at low temperatures with the application of a magnetic field. Often, this is accompanied by a huge and nonsaturating linear magnetoresistance which remains difficult to explain. Here we present a systematic study of the magnetoresistance in single-crystal $\gamma$-PtBi$_2$. We observe that the angle between the magnetic field and the crystalline $c$ axis fundamentally changes the magnetoresistance, going from a saturating to a nonsaturating magnetic field dependence. In between, there is one specific angle where the magnetoresistance is perfectly linear with the magnetic field. We show that the linear dependence of the nonsaturating magnetoresistance is due to the formation of open orbits in the Fermi surface of $\gamma$-PtBi$_2$.


title: Active learning phase boundaries of a quantum many-body system

abstract: We describe how to use techniques from the field of Machine Learning to direct a variational energy minimization scheme to search for phase boundaries of a quantum many-body system. The modeled physical system presents states of finite momentum condensate, also known as FFLO states, as well as a uniform superfluid phase-all of which is interesting in its own right; however, a full description of the multitude of phase boundaries is expensive from a computational standpoint. In this work, we treat the output of the energy minimization as a labeled sythetic data set to train a support vector classifier to separate states of finite momentum condensate from superfluid and normal states. We can then use the trained support vector classifier to refocus the minimizer to intensify its calculations near the boundary separating each of the three regions. Doing so will preclude using the minimizer to perform expensive calculations deep within the normal or superfluid regions, resulting in more efficient use of compute time. The application of the procedure we describe is straightforward and should be applicable in any computational search of phase boundaries.


title: Exploiting Cloze Questions for Few Shot Text Classification and Natural Language Inference

abstract: Some NLP tasks can be solved in a fully unsupervised fashion by providing a pretrained language model with "task descriptions" in natural language (e.g., Radford et al., 2019). While this approach underperforms its supervised counterpart, we show in this work that the two ideas can be combined: We introduce Pattern-Exploiting Training (PET), a semi-supervised training procedure that reformulates input examples as cloze-style phrases to help language models understand a given task. These phrases are then used to assign soft labels to a large set of unlabeled examples. Finally, standard supervised training is performed on the resulting training set. For several tasks and languages, PET outperforms supervised training and strong semi-supervised approaches in low-resource settings by a large margin.


title: Flat extensions of abstract polytopes

abstract: We consider the problem of constructing an abstract $(n+1)$-polytope $Q$ with $k$ facets isomorphic to a given $n$-polytope $P$, where $k \geq 3$. In particular, we consider the case where we want $Q$ to be $(n-2,n)$-flat, meaning that every $(n-2)$-face is incident to every $n$-face (facet). We show that if $P$ admits such a flat extension for a given $k$, then the facet graph of $P$ is $(k-1)$-colorable. Conversely, we show that if the facet graph is $(k-1)$-colorable and $k-1$ is prime, then $P$ admits a flat extension for that $k$. We also show that if $P$ is facet-bipartite, then for every even $k$, there is a flat extension $P|k$ such that every automorphism of $P$ extends to an automorphism of $P|k$. Finally, if $P$ is a facet-bipartite $n$-polytope and $Q$ is a vertex-bipartite $m$-polytope, we describe a flat amalgamation of $P$ and $Q$, an $(m+n-1)$-polytope that is $(n-2,n)$-flat, with $n$-faces isomorphic to $P$ and co-$(n-2)$-faces isomorphic to $Q$.


title: Iterator-Based Temporal Logic Task Planning

abstract: Temporal logic task planning for robotic systems suffers from state explosion when specifications involve large numbers of discrete locations. We provide a novel approach, particularly suited for tasks specifications with universally quantified locations, that has constant time with respect to the number of locations, enabling synthesis of plans for an arbitrary number of them. We propose a hybrid control framework that uses an iterator to manage the discretised workspace hiding it from a plan enacted by a discrete event controller. A downside of our approach is that it incurs in increased overhead when executing a synthesised plan. We demonstrate that the overhead is reasonable for missions of a fixed-wing Unmanned Aerial Vehicle in simulated and real scenarios for up to 700000 locations.


title: Stochastic Finite State Control of POMDPs with LTL Specifications

abstract: Partially observable Markov decision processes (POMDPs) provide a modeling framework for autonomous decision making under uncertainty and imperfect sensing, e.g. robot manipulation and self-driving cars. However, optimal control of POMDPs is notoriously intractable. This paper considers the quantitative problem of synthesizing sub-optimal stochastic finite state controllers (sFSCs) for POMDPs such that the probability of satisfying a set of high-level specifications in terms of linear temporal logic (LTL) formulae is maximized. We begin by casting the latter problem into an optimization and use relaxations based on the Poisson equation and McCormick envelopes. Then, we propose an stochastic bounded policy iteration algorithm, leading to a controlled growth in sFSC size and an any time algorithm, where the performance of the controller improves with successive iterations, but can be stopped by the user based on time or memory considerations. We illustrate the proposed method by a robot navigation case study.


title: Quantum Cosmology of Fab Four John Theory with Conformable Fractional Derivative

abstract: We study a quantization via fractional derivative of a nonminimal derivative coupling cosmological theory, namely, the Fab Four John theory. Its Hamiltonian version presents the issue of fractional powers in the momenta. That problem is solved here by the application of the so-called conformable fractional derivative. This leads to a Wheeler-DeWitt equation of second order, showing that a Bohm-de Broglie interpretation can be constructed. That combination of fractional quantization and Bohmian interpretation provides us a new quantization method, in which the quantum potential is the criterion to say if a quantum solution is acceptable or not to be further studied. We show that a wide range of solutions for the scale factor is possible. Among all of those, a bouncing solution analogous to the perfect fluid cosmology seems to deserve special attention.


title: Non-loose negative torus knots

abstract: We study Legendrian and transverse realizations of the negative torus knots $T_{(p,-q)}$ in all contact structures on the $3$-sphere. We give a complete classification of the strongly non-loose transverse realizations and the strongly non-loose Legendrian realizations with the Thurston-Bennequin invariant smaller than $-pq$. Additionally, we show that the strongly non-loose transverse realizations $T$ are classified by their non-zero invariants $\mathfrak T(T)$ in the minus version of the knot Floer homology. However, not all the elements of $HFK^-(T_{(p,q)})$ can be realized. Along the way, we relate our Legendrian realizations to the tight contact structures on the Legendrian surgeries along them. Specifically, we realize all tight structures on the lens spaces $L(pq+1,p^2)$ as a single Legendrian surgery on a Legendrian $T_{(p,-q)}$, and we relate transverse realizations in overtwisted structures to the non-fillable tight structures on the large negative surgeries along the underlying knots.


title: The finite $N$ origin of the Bardeen-Moshe-Bander phenomenon and its extension at $N=\infty$ by singular fixed points

abstract: We study the $O(N)$ model in dimension three (3$d$) at large and infinite $N$ and show that the line of fixed points found at $N=\infty$ --the Bardeen-Moshe-Bander (BMB) line-- has an intriguing origin at finite $N$. The large $N$ limit that allows us to find the BMB line must be taken on particular trajectories in the $(d,N)$-plane: $d=3-\alpha/N$ and not at fixed dimension $d=3$. Our study also reveals that the known BMB line is only half of the true line of fixed points, the second half being made of singular fixed points. The potentials of these singular fixed points show a cusp for a finite value of the field and their finite $N$ counterparts a boundary layer.


title: Dynamical Detection of Level Repulsion in the One-Particle Aubry-Andr\'e Model

abstract: The analysis of level statistics provides a primary method to detect signatures of chaos in the quantum domain. However, for experiments with ion traps and cold atoms, the energy levels are not as easily accessible as the dynamics. In this work, we discuss how properties of the spectrum that are usually associated with chaos can be directly detected from the evolution of the number operator in the one-dimensional, noninteracting Aubry-Andr\'e model. Both the quantity and the model are studied in experiments with cold atoms. We consider a single-particle and system sizes experimentally reachable. By varying the disorder strength within values below the critical point of the model, level statistics similar to those found in random matrix theory are obtained. Dynamically, these properties of the spectrum are manifested in the form of a dip below the equilibration point of the number operator. This feature emerges at times that are experimentally accessible. This work is a contribution to a special issue dedicated to Shmuel Fishman.


title: Stellar metallicities from SkyMapper photometry I: A study of the Tucana II ultra-faint dwarf galaxy

abstract: We present a study of the ultra-faint Milky Way dwarf satellite galaxy Tucana II using deep photometry from the 1.3m SkyMapper telescope at Siding Spring Observatory, Australia. The SkyMapper filter-set contains a metallicity-sensitive intermediate-band $v$ filter covering the prominent Ca II K feature at 3933.7A. When combined with photometry from the SkyMapper $u, g$, and $i$ filters, we demonstrate that $v$ band photometry can be used to obtain stellar metallicities with a precision of $\sim0.20$dex when [Fe/H] $> -2.5$, and $\sim0.34$dex when [Fe/H] $< -2.5$. Since the $u$ and $v$ filters bracket the Balmer Jump at 3646A, we also find that the filter-set can be used to derive surface gravities. We thus derive photometric metallicities and surface gravities for all stars down to a magnitude of $g\sim20$ within $\sim$75 arcminutes of Tucana II. Photometric metallicity and surface gravity cuts remove nearly all foreground contamination. By incorporating Gaia proper motions, we derive quantitative membership probabilities which recover all known members on the red giant branch of Tucana II. Additionally, we identify multiple likely new members in the center of the system and candidate members several half-light radii from the center of the system. Finally, we present a metallicity distribution function derived from the photometric metallicities of likely Tucana II members. This result demonstrates the utility of wide-field imaging with the SkyMapper filter-set in studying UFDs, and in general, low surface brightness populations of metal-poor stars. Upcoming work will clarify the membership status of several distant stars identified as candidate members of Tucana II.


title: FixMatch: Simplifying Semi-Supervised Learning with Consistency and Confidence

abstract: Semi-supervised learning (SSL) provides an effective means of leveraging unlabeled data to improve a model's performance. In this paper, we demonstrate the power of a simple combination of two common SSL methods: consistency regularization and pseudo-labeling. Our algorithm, FixMatch, first generates pseudo-labels using the model's predictions on weakly-augmented unlabeled images. For a given image, the pseudo-label is only retained if the model produces a high-confidence prediction. The model is then trained to predict the pseudo-label when fed a strongly-augmented version of the same image. Despite its simplicity, we show that FixMatch achieves state-of-the-art performance across a variety of standard semi-supervised learning benchmarks, including 94.93% accuracy on CIFAR-10 with 250 labels and 88.61% accuracy with 40 -- just 4 labels per class. Since FixMatch bears many similarities to existing SSL methods that achieve worse performance, we carry out an extensive ablation study to tease apart the experimental factors that are most important to FixMatch's success. We make our code available at https://github.com/google-research/fixmatch.


title: BLE Beacons for Indoor Positioning at an Interactive IoT-Based Smart Museum

abstract: The Internet of Things (IoT) can enable smart infrastructures to provide advanced services to the users. New technological advancement can improve our everyday life, even simple tasks as a visit to the museum. In this paper, an indoor localization system is presented, to enhance the user experience in a museum. In particular, the proposed system relies on Bluetooth Low Energy (BLE) beacons proximity and localization capabilities to automatically provide the users with cultural contents related to the observed artworks. At the same time, an RSS-based technique is used to estimate the location of the visitor in the museum. An Android application is developed to estimate the distance from the exhibits and collect useful analytics regarding each visit and provide a recommendation to the users. Moreover, the application implements a simple Kalman filter in the smartphone, without the need of the Cloud, to improve localization precision and accuracy. Experimental results on distance estimation, location, and detection accuracy show that BLE beacon is a promising solution for an interactive smart museum. The proposed system has been designed to be easily extensible to the IoT technologies and its effectiveness has been evaluated through experimentation.


title: A conservation law for testing methods of prediction of the seismic wave response of a protuberance emerging from flat ground

abstract: We establish the equations which translate a conservation law for the problem of the seismic response of an above-ground structure (e.g., building, hill or mountain) of arbitrary shape and inquire whether both the implicit (formal) and explicit (numerical) solutions for the response obey this law for the case of a cylindrical, rectangular protuberance. Both the low-order (poor approximations of the response) as well as higher-order (supposedly better approximations) turn out to satisfy the conservation of flux relation, which means that the satisfaction of this relation is a necessary, but not sufficient, means for determining whether a solution to the scattering problem is valid.


title: Estimating international trade status of countries from global liner shipping networks

abstract: Maritime shipping is a backbone of international trade and, thus, the world economy. Cargo-loaded vessels travel from one country's port to another via an underlying port-to-port transport network, contributing to international trade values of countries en route. We hypothesize that ports that involve trans-shipment activities serve as a third-party broker to mediate trade between two foreign countries and contribute to the corresponding country's status in international trade. We test this hypothesis using a port-level dataset of global liner shipping services. We propose two indices that quantify the importance of countries in the global liner shipping network and show that they explain a large amount of variation in individual countries' international trade values and related measures. These results support a long-standing view in maritime economics, which has yet to be directly tested, that countries that are strongly integrated into the global maritime transportation network have enhanced access to global markets and trade opportunities.


title: On the Existence of Minimal Hypersurfaces with Arbitrarily Large Area and Morse Index

abstract: We show that a bumpy closed Riemannian manifold $(M^{n+1}, g)$ $(3 \leq n+1 \leq 7)$ admits a sequence of connected closed embedded two-sided minimal hypersurfaces whose areas and Morse indices both tend to infinity. This improves a previous result by O. Chodosh and C. Mantoulidis on connected minimal hypersurfaces with arbitrarily large area.


title: Understanding the role of phenotypic switching in cancer drug resistance

abstract: The emergence of acquired drug resistance in cancer represents a major barrier to treatment success. While research has traditionally focused on genetic sources of resistance, recent findings suggest that cancer cells can acquire transient resistant phenotypes via epigenetic modifications and other non-genetic mechanisms. Although these resistant phenotypes are eventually relinquished by individual cells, they can temporarily 'save' the tumor from extinction and enable the emergence of more permanent resistance mechanisms. These observations have generated interest in the potential of epigenetic therapies for long-term tumor control or eradication. In this work, we develop a mathematical model to study how phenotypic switching at the single-cell level affects resistance evolution in cancer. We highlight unique features of non-genetic resistance, probe the evolutionary consequences of epigenetic drugs and explore potential therapeutic strategies. We find that even short-term epigenetic modifications and stochastic fluctuations in gene expression can drive long-term drug resistance in the absence of any bona fide resistance mechanisms. We also find that an epigenetic drug that slightly perturbs the average retention of the resistant phenotype can turn guaranteed treatment failure into guaranteed success. Lastly, we find that combining an epigenetic drug with an anti-cancer agent can significantly outperform monotherapy, and that treatment outcome is heavily affected by drug sequencing.


title: Identifying New X-ray Binary Candidates in M31 using Random Forest Classification

abstract: Identifying X-ray binary (XRB) candidates in nearby galaxies requires distinguishing them from possible contaminants including foreground stars and background active galactic nuclei. This work investigates the use of supervised machine learning algorithms to identify high-probability X-ray binary candidates. Using a catalogue of 943 Chandra X-ray sources in the Andromeda galaxy, we trained and tested several classification algorithms using the X-ray properties of 163 sources with previously known types. Amongst the algorithms tested, we find that random forest classifiers give the best performance and work better in a binary classification (XRB/non-XRB) context compared to the use of multiple classes. Evaluating our method by comparing with classifications from visible-light and hard X-ray observations as part of the Panchromatic Hubble Andromeda Treasury, we find compatibility at the 90% level, although we caution that the number of sources in common is rather small. The estimated probability that an object is an X-ray binary agrees well between the random forest binary and multiclass approaches and we find that the classifications with the highest confidence are in the X-ray binary class. The most discriminating X-ray bands for classification are the 1.7-2.8, 0.5-1.0, 2.0-4.0, and 2.0-7.0 keV photon flux ratios. Of the 780 unclassified sources in the Andromeda catalogue, we identify 16 new high-probability X-ray binary candidates and tabulate their properties for follow-up.


title: Estimating locomotor demands during team play from broadcast-derived tracking data

abstract: The introduction of optical tracking data across sports has given rise to the ability to dissect athletic performance at a level unfathomable a decade ago. One specific area that has seen substantial benefit is sports science, as high resolution coordinate data permits sports scientists to have to-the-second estimates of external load metrics, such as acceleration load and high speed running distance, traditionally used to understand the physical toll a game takes on an athlete. Unfortunately, collecting this data requires installation of expensive hardware and paying costly licensing fees to data providers, restricting its availability. Algorithms have been developed that allow a traditional broadcast feed to be converted to x-y coordinate data, making tracking data easier to acquire, but coordinates are available for an athlete only when that player is within the camera frame. Obviously, this leads to inaccuracies in player load estimates, limiting the usefulness of this data for sports scientists. In this research, we develop models that predict offscreen load metrics and demonstrate the viability of broadcast-derived tracking data for understanding external load in soccer.


title: Effective equidistribution of expanding horospheres in the locally symmetric space SO(d) \ SL(d,R) / SL(d,Z)

abstract: We use a dictionary between lattice point counting inside dilated d-dimensional ellipsoids (Euclidean counting) and counting of lifts of a closed horosphere that intersect a ball of increasing radius, to obtain two types of results. Firstly, via an $L^2$-integral error estimate for Euclidean counting, we prove effective equidistribution results for a family of expanding horospheres in the locally symmetric space SO(d) \ SL(d,R) / SL(d,Z). Secondly, we derive from uniform error estimates in Euclidean counting, error terms for counting SL(d,Z)-orbit points in a certain increasing family of subsets in SO(d) \ SL(d,R) (which we call truncated chimneys), and for counting the number of lifts of a closed horosphere that intersect a ball with large radius.


title: The Echo Index and multistability in input-driven recurrent neural networks

abstract: A recurrent neural network (RNN) possesses the echo state property (ESP) if, for a given input sequence, it ``forgets'' any internal states of the driven (nonautonomous) system and asymptotically follows a unique, possibly complex trajectory. The lack of ESP is conventionally understood as a lack of reliable behaviour in RNNs. Here, we show that RNNs can reliably perform computations under a more general principle that accounts only for their local behaviour in phase space. To this end, we formulate a generalisation of the ESP and introduce an echo index to characterise the number of simultaneously stable responses of a driven RNN. We show that it is possible for the echo index to change with inputs, highlighting a potential source of computational errors in RNNs due to characteristics of the inputs driving the dynamics.


title: Minimal $2$-Spheres and Optimal Foliations in $3$-Spheres with Arbitrary Metric

abstract: In this paper, we prove that the $3$-sphere endowed with an arbitrary Riemannian metric either contains at least two embedded minimal $2$-spheres or admits an optimal foliation by $2$-spheres. This generalizes recent results by Haslhofer-Ketover (Duke Math. J. 2019), where the existence of optimal foliations and minimal $2$-spheres has been established under the additional assumption that the metric is generic. In light of recent examples by Wang-Zhou, where min-max for some non-bumpy metrics on the 3-sphere produces higher multiplicities, our results are in a certain sense sharp.


title: Charging assisted by thermalization

abstract: A system in thermal equilibrium with a bath will generally be in an athermal state, if the system-bath coupling is strong. In some cases, it will be possible to extract work from that athermal state, after disconnecting the system from the bath. We use this observation to devise a battery charging and storing unit, simply consisting of a system, acting as the battery, and a bath. The charging cycle---connect, let thermalize, disconnect, extract work---requires very little external control and the charged state of the battery, being a part of global thermal equilibrium, can be maintained indefinitely and for free. The efficiency, defined as the ratio of the extractable work stored in the battery and the total work spent on connecting and disconnecting, is always $\leq 1$, which is a manifestation of the second law of thermodynamics. Moreover, coupling, being a resource for the device, is also a source of dissipation: the entropy production per charging cycle is always significant, strongly limiting the efficiency in all coupling strength regimes. We show that our general results also hold for generic microcanonical baths. We illustrate our theory on the Caldeira-Leggett model with a harmonic oscillator (the battery) coupled to a harmonic bath, for which we derive general asymptotic formulas in both weak and ultrastrong coupling regimes, for arbitrary Ohmic spectral densities. We show that the efficiency can be increased by connecting several copies of the battery to the bath. Finally, as a side result, we derive a general formula for Gaussian ergotropy, that is, the maximal work extractable by Gaussian unitary operations from Gaussian states of multipartite continuous-variable systems.


title: Stochastic Approximation versus Sample Average Approximation for population Wasserstein barycenters

abstract: In the machine learning and optimization community, there are two main approaches for the convex risk minimization problem, namely, the Stochastic Approximation (SA) and the Sample Average Approximation (SAA). In terms of oracle complexity (required number of stochastic gradient evaluations), both approaches are considered equivalent on average (up to a logarithmic factor). The total complexity depends on the specific problem, however, starting from work \cite{nemirovski2009robust} it was generally accepted that the SA is better than the SAA. % Nevertheless, in case of large-scale problems SA may run out of memory as storing all data on one machine and organizing online access to it can be impossible without communications with other machines. SAA in contradistinction to SA allows parallel/distributed calculations. We show that for the Wasserstein barycenter problem this superiority can be inverted. We provide a detailed comparison by stating the complexity bounds for the SA and the SAA implementations calculating barycenters defined with respect to optimal transport distances and entropy-regularized optimal transport distances. As a byproduct, we also construct confidence intervals for the barycenter defined with respect to entropy-regularized optimal transport distances in the $\ell_2$-norm. The preliminary results are derived for a general convex optimization problem given by the expectation in order to have other applications besides the Wasserstein barycenter problem.


title: Intelligent Bandwidth Allocation for Latency Management in NG-EPON using Reinforcement Learning Methods

abstract: A novel intelligent bandwidth allocation scheme in NG-EPON using reinforcement learning is proposed and demonstrated for latency management. We verify the capability of the proposed scheme under both fixed and dynamic traffic loads scenarios to achieve <1ms average latency. The RL agent demonstrates an efficient intelligent mechanism to manage the latency, which provides a promising IBA solution for the next-generation access network.


title: Field-Driven Gapless Spin Liquid in the Spin-1 Kitaev Honeycomb Model

abstract: Recent proposals for spin-1 Kitaev materials, such as honeycomb Ni oxides with heavy elements of Bi and Sb, have shown that these compounds naturally give rise to antiferromagnetic (AFM) Kitaev couplings. Conceptual interest in such AFM Kitaev systems has been sparked by the observation of a transition to a gapless $U(1)$ spin liquid at intermediate field strengths in the AFM spin-1/2 Kitaev model. However, all hitherto known spin-1/2 Kitaev materials exhibit ferromagnetic bond-directional exchanges. Here we discuss the physics of the spin-1 Kitaev model in a magnetic field and show, by extensive numerical analysis, that for AFM couplings it exhibits an extended gapless quantum spin liquid at intermediate field strengths. The close analogy to its spin-1/2 counterpart suggests that this gapless spin liquid is a $U(1)$ spin liquid with a neutral Fermi surface, that gives rise to enhanced thermal transport signatures.


title: Fundamental Plane of BOSS galaxies: Correlations with galaxy properties, density field and impact on RSD measurements

abstract: Fundamental plane of elliptical galaxies can be used to predict the intrinsic size of galaxies and has a number of plausible application to study cosmology and galaxy physics. We present a detailed analysis of the fundamental plane of the SDSS-III BOSS LOWZ and CMASS galaxies. For the standard fundamental plane, we find a strong redshift evolution for the mean residual and show that it is primarily driven by the redshift evolution of the surface brightness of the galaxies. After correcting for the redshift evolution, the FP residuals are strongly correlated with the galaxy properties and some observational systematics. We show that the variations in the FP between the central and satellite galaxies, that have been observed in the literature, can primarily be explained by the correlation of the FP with the galaxy luminosity. We also measure the cross correlations of the FP residuals with the galaxy density field. The amplitude of the cross correlations depends on the galaxy properties and environment with brighter and redder galaxies showing stronger correlation. In general, galaxies in denser environments (higher galaxy bias ) show stronger correlations. We also compare FP amplitude with the amplitudes of intrinsic alignments of galaxy shapes (IA), finding the two to be correlated. Finally, using the FP residuals we also study the impact of intrinsic alignments on the constraint of growth rate using redshift space distortions. We do not observe any significant trends in measurements of the growth rate $f$ as function of the amplitude of FP-density correlations, resulting in null detection of the effects of IA on the RSD measurements.


title: A Re-Examination of the Evidence used by Hooge et al (2018) "Is human classification by experienced untrained observers a gold standard in fixation detection?"

abstract: Hooge et al. asked the question: "Is human classification by experienced untrained observers a gold standard in fixation detection?" They conclude the answer is no. If they had entitled their paper: "Is human classification by experienced untrained observers a gold standard in fixation detection when data quality is very poor, data are error-filled, data presentation was not optimal, and the analysis was seriously flawed?", I would have no case to make. In the present report, I will present evidence to support my view that this latter title is justified. The low quality data assessment is based on using a relatively imprecise eye-tracker, the absence of head restraint for any subjects, and the use of infants as the majority of subjects (60 of 70 subjects). Allowing subjects with more than 50% missing data (as much as 95%) is also evidence of low quality data. The error-filled assessment is based on evidence that a number of the "fixations" classified by "experts" have obvious saccades within them, and that, apparently, a number of fixations were classified on the basis of no signal at all. The evidence for non-optimal data presentation stems from the fact that, in a number of cases, perfectly good data was not presented to the coders. The flaws in the analysis are evidenced by the fact that entire stretches of missing data were considered classified, and that the measurement of saccade amplitude was based on many cases in which there was no saccade at all. Without general evidence to the contrary, it is correct to assume that some human classifiers under some conditions may meet the criteria for a gold standard, and classifiers under other conditions may not. This conditionality is not recognized by Hooge et al. A fair assessment would conclude that whether or not humans can be considered a gold standard is still very much an open question.


title: Structural Properties of an S-system Model of Mycobacterium Tuberculosis Gene Regulation

abstract: Magombedze and Mulder in 2013 studied the gene regulatory system of Mycobacterium Tuberculosis (Mtb) by partitioning this into three subsystems based on putative gene function and role in dormancy/latency development. Each subsystem, in the form of S-system, is represented by an embedded chemical reaction network (CRN), defined by a species subset and a reaction subset induced by the set of digraph vertices of the subsystem. For the embedded networks of S-system, we showed interesting structural properties and proved that all S-system CRNs (with at least two species) are discordant. Analyzing the subsystems as subnetworks, where arcs between vertices belonging to different subsystems are retained, we formed a digraph homomorphism from the corresponding subnetworks to the embedded networks. Lastly, we explored the modularity concept of CRN in the context of digraph.


title: Generalized Dobrushin Ergodicity Coefficient and Ergodicities of Non-homogeneous Markov Chains

abstract: In our earlier paper, a generalized Dobrushin ergodicity coefficient of Markov operators (acting on abstract state spaces) with respect to a projection $P$, has been introduced and studied. It turned out that the introduced coefficient was more effective than the usual ergodicity coefficient. In the present work, by means of a left consistent Markov projections and the generalized Dobrushin's ergodicity coefficient, we investigate uniform and weak $P$-ergodicities of non-homogeneous discrete Markov chains (NDMC) on abstract state spaces. It is easy to show that uniform $P$-ergodicity implies a weak one, but in general the reverse is not true. Therefore, some conditions are provided together with weak $P$-ergodicity of NDMC which imply its uniform $P$-ergodicity. Furthermore, necessary and sufficient conditions are found by means of the Doeblin's condition for the weak $P$-ergodicity of NDMC. The weak $P$-ergodicity is also investigated in terms of perturbations. Several perturbative results are obtained which allow us to produce nontrivial examples of uniform and weak $P$-ergodic NDMC. Moreover, some category results are also obtained. We stress that all obtained results have potential applications in the classical and non-commutative probabilities.


title: Asynchronous Consensus Algorithm

abstract: This document describes a new consensus algorithm which is asynchronous and uses gossip based message dissemination between nodes. The current version of the algorithm does not cover the case of a node failure or significantly delayed response. This is the subject of further research of the algorithm. An outline of a new design for trust-less payment system is given in appendices.


title: Superexponential Interactions and the Dynamical Unfolding of Confined Degrees of Freedom

abstract: We explore a two-body system with superexponential interactions that serves as a fundamental building block for a route to complexity. While being of striking simplicity this highly nonlinear interaction yields a plethora of intriguing properties and a rich dynamics. It exhibits a spatial region where the dynamics occurs in a channel characterised by a transversally confined and longitudinally unbounded motion and additionally two distinct regions where the dynamics is asymptotically free. A deconfinement transition via two saddle points connects the dynamics in the channel with the asymptotically free motion. The scattering functions show plateau and peak structures that can be interpreted in terms of corresponding correlation diagrams. These are intimately related to the varying anharmonicity of the transverse motion while moving along the longitudinal dimension of the channel. We perform a comprehensive analysis of the scattering transition for energies below and above the saddle points. Possible variants and extensions of the superexponential interaction to many-body system are briefly discussed.


title: Towards Digital Twins for the Description of Automotive Software Systems

abstract: We present models for automotive software that capture quantitative and qualitative aspects of software systems and the underlying hardware architecture. In particular, we consider different levels of computing power. These range from controllers up to the cloud. We present a modeling approach for software deployment taking different automotive requirements such as criticality, latency, memory, computational resources, and communication into account. Our models capture automotive software and hardware system configurations and can serve as digital twins that are digital counterparts of (usually) physical entities. Furthermore, we highlight connected research areas and challenges.


title: Tomographic and entropic analysis of modulated signals

abstract: We study an application of the quantum tomography framework for the time-frequency analysis of modulated signals. In particular, we calculate optical tomographic representations and Wigner-Ville distributions for signals with amplitude and frequency modulations. We also consider time-frequency entropic relations for modulated signals, which are naturally associated with the Fourier analysis. A numerical toolbox for calculating optical time-frequency tomograms based on pseudo Wigner-Ville distributions for modulated signals is provided.


title: Towards Comparability in Non-Intrusive Load Monitoring: On Data and Performance Evaluation

abstract: Non-Intrusive Load Monitoring (NILM) comprises of a set of techniques that provide insights into the energy consumption of households and industrial facilities. Latest contributions show significant improvements in terms of accuracy and generalisation abilities. Despite all progress made concerning disaggregation techniques, performance evaluation and comparability remains an open research question. The lack of standardisation and consensus on evaluation procedures makes reproducibility and comparability extremely difficult. In this paper, we draw attention to comparability in NILM with a focus on highlighting the considerable differences amongst common energy datasets used to test the performance of algorithms. We divide discussion on comparability into data aspects, performance metrics, and give a close view on evaluation processes. Detailed information on pre-processing as well as data cleaning methods, the importance of unified performance reporting, and the need for complexity measures in load disaggregation are found to be the most urgent issues in NILM-related research. In addition, our evaluation suggests that datasets should be chosen carefully. We conclude by formulating suggestions for future work to enhance comparability.


title: Adaptive Large Neighborhood Search for Circle Bin Packing Problem

abstract: We address a new variant of packing problem called the circle bin packing problem (CBPP), which is to find a dense packing of circle items to multiple square bins so as to minimize the number of used bins. To this end, we propose an adaptive large neighborhood search (ALNS) algorithm, which uses our Greedy Algorithm with Corner Occupying Action (GACOA) to construct an initial layout. The greedy solution is usually in a local optimum trap, and ALNS enables multiple neighborhood search that depends on the stochastic annealing schedule to avoid getting stuck in local minimum traps. Specifically, ALNS perturbs the current layout to jump out of a local optimum by iteratively reassigns some circles and accepts the new layout with some probability during the search. The acceptance probability is adjusted adaptively using simulated annealing that fine-tunes the search direction in order to reach the global optimum. We benchmark computational results against GACOA in heterogeneous instances. ALNS always outperforms GACOA in improving the objective function, and in several cases, there is a significant reduction on the number of bins used in the packing.


title: An Image Enhancing Pattern-based Sparsity for Real-time Inference on Mobile Devices

abstract: Weight pruning has been widely acknowledged as a straightforward and effective method to eliminate redundancy in Deep Neural Networks (DNN), thereby achieving acceleration on various platforms. However, most of the pruning techniques are essentially trade-offs between model accuracy and regularity which lead to impaired inference accuracy and limited on-device acceleration performance. To solve the problem, we introduce a new sparsity dimension, namely pattern-based sparsity that comprises pattern and connectivity sparsity, and becoming both highly accurate and hardware friendly. With carefully designed patterns, the proposed pruning unprecedentedly and consistently achieves accuracy enhancement and better feature extraction ability on different DNN structures and datasets, and our pattern-aware pruning framework also achieves pattern library extraction, pattern selection, pattern and connectivity pruning and weight training simultaneously. Our approach on the new pattern-based sparsity naturally fits into compiler optimization for highly efficient DNN execution on mobile platforms. To the best of our knowledge, it is the first time that mobile devices achieve real-time inference for the large-scale DNN models thanks to the unique spatial property of pattern-based sparsity and the help of the code generation capability of compilers.


title: Diffraction of magnetic field in a device of circular rings

abstract: A device consisted of a set of circular rings, the centers of which lie on an axis, behaves like a solenoid when the ratio of its radius and distance between two successive rings is greater than one. As this ratio decreases, the device deviates a lot from the solenoid behavior.In the same way, a diffraction phenomenon for magnetic field appears when currents of random direction flow through the rings. This phenomenon demonstrates a critical behavior. Thus an extension of diffraction phenomenon can be done beyond the classical wave diffraction. A possible application of this device could be the diffraction of electronic beams.


title: SMAPGAN: Generative Adversarial Network Based Semi-Supervised Styled Map Tiles Generating Method

abstract: Traditional online map tiles, widely used on the Internet such as Google Map and Baidu Map, are rendered from vector data. Timely updating online map tiles from vector data, of which the generating is time-consuming, is a difficult mission. It is a shortcut to generate map tiles in time from remote sensing images, which can be acquired timely without vector data. However, this mission used to be challenging or even impossible. Inspired by image-to-image translation (img2img) techniques based on generative adversarial networks (GAN), we proposed a semi-supervised Generation of styled map Tiles based on Generative Adversarial Network (SMAPGAN) model to generate styled map tiles directly from remote sensing images. In this model, we designed a semi-supervised learning strategy to pre-train SMAPGAN on rich unpaired samples and fine-tune it on limited paired samples in reality. We also designed image gradient L1 loss and image gradient structure loss to generate a styled map tile with global topological relationships and detailed edge curves of objects, which are important in cartography. Moreover, we proposed edge structural similarity index (ESSI) as a metric to evaluate the quality of topological consistency between generated map tiles and ground truths. Experimental results present that SMAPGAN outperforms state-of-the-art (SOTA) works according to mean squared error, structural similarity index, and ESSI. Also, SMAPGAN won more approval than SOTA in the human perceptual test on the visual realism of cartography. Our work shows that SMAPGAN is potentially a new paradigm to produce styled map tiles. Our implementation of the SMAPGAN is available at https://github.com/imcsq/SMAPGAN.


title: Revealing the Two-Electron Cusp in the Ground States of He and H2 via Quasifree Double Photoionization

abstract: We report on kinematically complete measurements and ab initio non-perturbative calculations of double ionization of He and H2 by a single 800 eV circularly polarized photon. We confirm the quasifree mechanism of photoionization for H2 and show how it originates from the two-electron cusp in the ground state of a two-electron target. Our approach establishes a new method for mapping electrons relative to each other and provides valuable insight into photoionization beyond the electric-dipole approximation.


title: Annealing and replica-symmetry in Deep Boltzmann Machines

abstract: In this paper we study the properties of the quenched pressure of a multi-layer spin-glass model (a deep Boltzmann Machine in artificial intelligence jargon) whose pairwise interactions are allowed between spins lying in adjacent layers and not inside the same layer nor among layers at distance larger than one. We prove a theorem that bounds the quenched pressure of such a K-layer machine in terms of K Sherrington-Kirkpatrick spin glasses and use it to investigate its annealed region. The replica-symmetric approximation of the quenched pressure is identified and its relation to the annealed one is considered. The paper also presents some observation on the model's architectural structure related to machine learning. Since escaping the annealed region is mandatory for a meaningful training, by squeezing such region we obtain thermodynamical constraints on the form factors. Remarkably, its optimal escape is achieved by requiring the last layer to scale sub-linearly in the network size.

categories: cond-mat.stat-mech cond-mat.dis-nn math-ph math.MP


title: TEASER: Fast and Certifiable Point Cloud Registration

abstract: We propose the first fast and certifiable algorithm for the registration of two sets of 3D points in the presence of large amounts of outlier correspondences. We first reformulate the registration problem using a Truncated Least Squares (TLS) cost that is insensitive to a large fraction of spurious correspondences. Then, we provide a general graph-theoretic framework to decouple scale, rotation, and translation estimation, which allows solving in cascade for the three transformations. Despite the fact that each subproblem is still non-convex and combinatorial in nature, we show that (i) TLS scale and (component-wise) translation estimation can be solved in polynomial time via adaptive voting, (ii) TLS rotation estimation can be relaxed to a semidefinite program (SDP) and the relaxation is tight, even in the presence of extreme outlier rates, and (iii) the graph-theoretic framework allows drastic pruning of outliers by finding the maximum clique. We name the resulting algorithm TEASER (Truncated least squares Estimation And SEmidefinite Relaxation). While solving large SDP relaxations is typically slow, we develop a second fast and certifiable algorithm, named TEASER++, that uses graduated non-convexity to solve the rotation subproblem and leverages Douglas-Rachford Splitting to efficiently certify global optimality. For both algorithms, we provide theoretical bounds on the estimation errors, which are the first of their kind for robust registration problems. Moreover, we test their performance on standard, object detection, and the 3DMatch benchmarks, and show that (i) both algorithms dominate the state of the art and are robust to more than 99% outliers, (ii) TEASER++ can run in milliseconds, and (iii) TEASER++ is so robust it can also solve problems without correspondences, where it largely outperforms ICP and it is more accurate than Go-ICP while being orders of magnitude faster.


title: Microscopic pairing mechanism, order parameter, and disorder sensitivity in moir\'e superlattices: Applications to twisted double-bilayer graphene

abstract: Starting from a continuum-model description, we develop a microscopic weak-coupling theory for superconductivity in twisted double-bilayer graphene. We study both electron-phonon and entirely electronic pairing mechanisms. In each case, the leading superconducting instability transforms under the trivial representation, $A$, of the point group $C_3$ of the system, while the subleading pairing phases belong to the $E$ channel. We explicitly compute the momentum dependence of the associated order parameters and find that the leading state has no nodal points for electron-phonon pairing but exhibits six sign changes on the Fermi surface if the Coulomb interaction dominates. On top of these system-specific considerations, we also present general results relevant to other correlated graphene-based moir\'e superlattice systems. We show that, irrespective of microscopic details, triplet pairing will be stabilized if the collective electronic fluctuations breaking the enhanced $\text{SU}(2)_+ \times \text{SU}(2)_-$ spin symmetry of these systems are odd under time reversal, even when the main $\text{SU}(2)_+ \times \text{SU}(2)_-$-symmetric part of the pairing glue is provided by phonons. Furthermore, we discuss the disorder sensitivity of the candidate pairing states and demonstrate that the triplet phase is protected against disorder on the moir\'e scale.


title: The Drell-Yan cross section to third order in the strong coupling constant

abstract: We present phenomenological results for the inclusive cross section for the production of a lepton-pair via virtual photon exchange at next-to-next-to-next-to-leading order (N$^3$LO) in perturbative QCD. In line with the case of Higgs production, we find that the hadronic cross section receives corrections at the percent level, and the residual dependence on the perturbative scales is reduced. However, unlike in the Higgs case, we observe that the uncertainty band derived from scale variation is no longer contained in the band of the previous order.


title: Second order Fermi reacceleration mechanisms and large scale synchrotron radio emission in intra-cluster bridges

abstract: Radio observations at low frequencies with the Low Frequency Array (LOFAR) start discovering gigantic radio bridges connecting pairs of massive galaxy clusters. These observations probe unexplored mechanisms of in situ particle acceleration that operate on volumes of several Mpc$^3$. Numerical simulations suggest that such bridges are dynamically complex and that weak shocks and super-Alfv\'{e}nic turbulence can be driven across the entire volume of these regions. In this Letter we explore, for the first time, the role of second order Fermi mechanisms for the reacceleration of relativistic electrons interacting with turbulence in these peculiar regions. We assume the turbulent energy flux measured in simulations and adopt a scenario in which relativistic particles scatter with magnetic field lines diffusing in super-Alfv\'enic turbulence and magnetic fields are amplified by the same turbulence. We show that steep spectrum and volume filling synchrotron emission can be generated in the entire intra-cluster bridge region thus providing a natural explanation for radio bridges. Consequently, radio observations have the potential to probe the dissipation of energy on scales larger than galaxy clusters and second order Fermi mechanisms operating in physical regimes that are still poorly explored. This has a potential impact on several branches of astrophysics and cosmology.


title: A versatile tool for cluster lensing source reconstruction. I. methodology and illustration on sources in the Hubble Frontier Field Cluster MACS J0717.5+3745

abstract: We describe a general purpose method to reconstruct the intrinsic properties of sources lensed by the gravitational potential of foreground clusters of galaxies. The tool Lenstruction is implemented in the publicly available multi-purpose gravitational lensing software Lenstronomy, in order to provide an easy and fast solution to this common astrophysical problem. The tool is based on forward modeling the appearance of the source in the image plane, taking into account the distortion by lensing and the instrumental point spread function (PSF). For singly-imaged sources a global lens model in the format of the Hubble Frontier Fields (HFF) lensing maps is required as a starting point. For multiply-imaged sources, the tool can also fit and apply first (deflection), second (shear, convergence), and third order (flexion) corrections to the local gravitational potential to improve the reconstruction, depending on the quality of the data. We illustrate the performance and features of the code with two examples of multiply-imaged systems taken from the Hubble Frontier Fields, starting from five different publicly available cluster models. We find that, after our correction, the relative magnification - and other lensing properties - between the multiple images become robustly constrained. Furthermore, we find that scatter between models of the reconstructed source size and magnitude is reduced. The code and jupyter notebooks are publicly available.


title: Nine localised deviations from Keplerian rotation in the DSHARP circumstellar disks: Kinematic evidence for protoplanets carving the gaps

abstract: We present evidence for localised deviations from Keplerian rotation, i.e., velocity "kinks", in 8 of 18 circumstellar disks observed by the DSHARP program: DoAr 25, Elias 2-27, GW Lup, HD 143006, HD 163296, IM Lup, Sz 129 and WaOph 6. Most of the kinks are detected over a small range in both radial extent and velocity, suggesting a planetary origin, but for some of them foreground contamination prevents us from measuring their spatial and velocity extent. Because of the DSHARP limited spectral resolution and signal-to-noise in the 12CO J=2-1 line, as well as cloud contamination, the kinks are usually detected in only one spectral channel, and will require confirmation. The strongest circumstantial evidence for protoplanets in the absence of higher spectral resolution data and additional tracers is that, upon deprojection, we find that all of the candidate planets lie within a gap and/or at the end of a spiral detected in dust continuum emission. This suggests that a significant fraction of the dust gaps and spirals observed by ALMA in disks are caused by embedded protoplanets.


title: Turbulence in a self-gravitating molecular cloud core

abstract: Externally driven interstellar turbulence plays an important role in shaping the density structure in molecular clouds. Here we study the dynamical role of internally driven turbulence in a self-gravitating molecular cloud core. Depending on the initial conditions and evolutionary stages, we find that a self-gravitating core in the presence of gravity-driven turbulence can undergo constant, decelerated, and accelerated infall, and thus has various radial velocity profiles. In the gravity-dominated central region, a higher level of turbulence results in a lower infall velocity, a higher density, and a lower mass accretion rate. As an important implication of this study, efficient reconnection diffusion of magnetic fields against the gravitational drag naturally occurs due to the gravity-driven turbulence, without invoking externally driven turbulence.


title: Parton Distributions in Nucleons and Nuclei

abstract: We review the current status of Parton Distribution Function (PDF) determinations for unpolarized and longitudinally polarized protons and for unpolarized nuclei, which are probed by high-energy hadronic scattering in perturbative Quantum Chromodynamics (QCD). We present the established theoretical framework, the experimental information, and the methodological aspects inherent to any modern PDF extraction. Furthermore, we summarize the present knowledge of PDFs and discuss their limitations in both accuracy and precision relevant to advance our understanding of QCD proton substructure and pursue our quest for precision in the Standard Model and beyond. In this respect, we highlight various achievements, discuss contemporary issues in PDF analyses, and outline future directions of progress.


title: Field Tensor Network States

abstract: We define a class of tensor network states for spin systems where the individual tensors are functionals of fields. The construction is based on the path integral representation of correlators of operators in quantum field theory. These tensor network states are infinite dimensional versions of matrix product states and projected entangled pair states. We find the field-tensor that generates the Haldane-Shastry wave function and extend it to two dimensions. We give evidence that the latter underlies the topological chiral state described by the Kalmeyer-Laughlin wave function.


title: Higher-Order Entanglement and Many-Body Invariants for Higher-Order Topological Phases

abstract: We discuss how strongly interacting higher-order symmetry protected topological (HOSPT) phases can be characterized from the entanglement perspective: First, we introduce a topological many-body invariant which reveals the non-commutative algebra between flux operator and $C_n$ rotations. We argue that this invariant denotes the angular momentum carried by the instanton which is closely related to the discrete Wen-Zee response and fractional corner charge. Second, we define a new entanglement property, dubbed `higher-order entanglement', to scrutinize and differentiate various higher-order topological phases from a hierarchical sequence of the entanglement structure. We support our claims by numerically studying a super-lattice Bose-Hubbard model that exhibits different HOSPT phases.

categories: cond-mat.str-el cond-mat.mes-hall cond-mat.quant-gas hep-th


title: A stellar census in globular clusters with MUSE -- Extending the CaT-metallicity relation below the horizontal branch and applying it to multiple populations

abstract: Aims. We use the spectra of more than 30,000 red giant branch (RGB) stars in 25 globular clusters (GC), obtained within the MUSE survey of Galactic globular clusters, to calibrate the Ca II triplet (CaT) metallicity relation and derive metallicities for all individual stars. We investigate the overall metallicity distributions as well as those of the different populations within each cluster. Methods. The Ca II triplet in the near-infrared at 8498, 8542, and 8662 AA is visible in stars with spectral types between F and M and can be used to determine their metallicities. In this work, we calibrate the relation using average cluster metallicities from literature and MUSE spectra, and extend it below the horizontal branch - a cutoff that has traditionally been made to avoid a non-linear relation - using a quadratic function. In addition to the classic relation based on V-VHB we also present calibrations based on absolute magnitude and luminosity. The obtained relations are then used to calculate metallicities for all the stars in the sample and to derive metallicity distributions for different populations within a cluster, which have been separated using so-called "chromosome maps" based on HST photometry. Results. We show that, despite the relatively low spectral resolution of MUSE (R=1900-3700) we can derive single star metallicities with a mean statistical intra-cluster uncertainty of ~0.12 dex. We present metallicity distributions for the RGB stars in 25 GCs, and investigate the different metallicities of the populations P3 (and higher) in so-called metal-complex or Type II clusters, finding metallicity variations in all of them. We also detected unexpected metallicity variations in the Type I cluster NGC 2808 and confirm the Type II status of NGC 7078.


title: Crescent states in charge-imbalanced polariton condensates

abstract: We study two-dimensional charge-imbalanced electron-hole systems embedded in an optical microcavity. We find that strong coupling to photons favors states with pairing at zero or small center of mass momentum, leading to a condensed state with spontaneously broken time-reversal and rotational symmetry, and unpaired carriers that occupy an anisotropic crescent-shaped sliver of momentum space. The crescent state is favoured at moderate charge imbalance, while a Fulde--Ferrel--Larkin--Ovchinnikov-like state --- with pairing at large center of mass momentum --- occurs instead at strong imbalance. The crescent state stability results from long-range Coulomb interactions in combination with extremely long-range photon-mediated interactions.


title: Measurement of the B-band Galaxy Luminosity Function with Approximate Bayesian Computation

abstract: The galaxy Luminosity Function (LF) is a key observable for galaxy formation, evolution studies and for cosmology. In this work, we propose a novel technique to forward model wide-field broad-band galaxy surveys using the fast image simulator UFig and measure the LF of galaxies in the B-band. We use Approximate Bayesian Computation (ABC) to constrain the galaxy population model parameters of the simulations and match data from CFHTLS. We define a number of distance metrics between the simulated and the survey data. By exploring the parameter space of the galaxy population model through ABC to find the set of parameters that minimize these distance metrics, we obtain constraints on the LFs of blue and red galaxies as a function of redshift. We find that $\mathrm{M^*}$ fades by $\Delta \mathrm{M}^*_{\mathrm{0.1-1.0,b}} = 0.68 \pm 0.52$ and $\Delta \mathrm{M}^*_{\mathrm{0.1-1.0,r}} = 0.54 \pm 0.48$ magnitudes between redshift $\mathrm{z = 1}$ and $\mathrm{z = 0.1}$ for blue and red galaxies, respectively. We also find that $\phi^*$ for blue galaxies stays roughly constant between redshift $\mathrm{z = 0.1}$ and $\mathrm{z=1}$, while for red galaxies it decreases by $\sim 35\%$. We compare our results to other measurements, finding good agreement at all redshifts, for both blue and red galaxies. To further test our results, we compare the redshift distributions for survey and simulated data. We use the spectroscopic redshift distribution from the VIMOS Public Extragalactic Redshift Survey (VIPERS) and we apply the same selection in colours and magnitudes on our simulations. We find a good agreement between the survey and the simulated redshift distributions. We provide best-fit values and uncertainties for the parameters of the LF. This work offers excellent prospects for measuring other galaxy population properties as a function of redshift using ABC.


title: Linear relativistic thermoelastic rod

abstract: We derive and analyze the linearized hyperbolic equations describing a relativistic heat-conducting elastic rod. We construct a decreasing energy integral for these equations, compute the associated characteristic propagation speeds, and prove that the solutions decay in time by using a Fourier decomposition. For comparison purposes, we obtain analogous results for the classical system with heat waves, in which the finite propagation speed of heat is kept but the other relativistic terms are neglected, and also for the usual classical system.


title: Study of the variable broadband emission of Markarian 501 during the most extreme Swift X-ray activity

abstract: Mrk501 is a very high-energy (VHE) gamma-ray blazar located at z=0.034. During a period of two weeks in July 2014, the highest X-ray activity of Mrk501 was observed in ~14 years of operation of the Neil Gehrels Swift Gamma-ray Burst Observatory. We characterize the broadband variability of Mrk501 from radio to VHE gamma rays, and evaluate whether it can be interpreted within theoretical scenarios widely used to explain the broadband emission from blazars. The temporal evolution of the most prominent and variable segments of the SED is described with a one-zone synchrotron self-Compton model with variations in the break energy of the electron energy distribution (EED), and with some adjustments in the magnetic field strength and spectral shape of the EED. A narrow feature at ~3 TeV was observed in the VHE spectrum measured on 2014 July 19 (MJD 56857.98), which is the day with the highest X-ray flux ($>0.3$ keV) measured during the entire Swift mission. This feature is inconsistent with the classical analytic functions to describe the measured VHE spectra (power law, log-parabola, and log-parabola with exponential cutoff) at more than 3$\sigma$. A fit with a log-parabola plus a narrow component is preferred over the fit with a single log-parabola at more than 4$\sigma$, and a dedicated Monte Carlo simulation estimated the significance of this extra component to be larger than 3$\sigma$. Under the assumption that this VHE spectral feature is real, we show that it can be reproduced with three distinct theoretical scenarios: a) a pileup in the EED due to stochastic acceleration; b) a structured jet with two-SSC emitting regions, with one region dominated by an extremely narrow EED; and c) an emission from an IC pair cascade.


title: The Clustering of X-ray Luminous Quasars

abstract: The clustering of active galactic nuclei (AGN) sheds light on their typical large (Mpc-scale) environments, which can constrain the growth and evolution of supermassive black holes. Here we measure the clustering of luminous X-ray-selected AGN in the Stripe 82X and XMM-XXL-North surveys around the peak epoch of black hole growth, in order to investigate the dependence of luminosity on large-scale AGN environment. We compute the auto-correlation function of AGN in two luminosity bins, $10^{43}\leq L_X<10^{44.5}$ erg s$^{-1}$ at $z\sim 0.8$ and $L_X\geq 10^{44.5}$ erg s$^{-1}$ at $z\sim 1.8$, and calculate the AGN bias taking into account the redshift distribution of the sources using three different methods. Our results show that while the less luminous sample has an inferred typical halo mass that is smaller than for the more luminous AGN, the host halo mass may be less dependent on luminosity than suggested in previous work. Focusing on the luminous sample, we calculate a typical host halo mass of $\sim 10^{13}$ M$_{\odot}~h^{-1}$, which is similar to previous measurements of moderate-luminosity X-ray AGN and significantly larger than the values found for optical quasars of similar luminosities and redshifts. We suggest that the clustering differences between different AGN selection techniques are dominated by selection biases, and not due to a dependence on AGN luminosity. We discuss the limitations of inferring AGN triggering mechanisms from halo masses derived by large-scale bias.