title: Learning relevant features for statistical inference

abstract: Given two views of data, we consider the problem of finding the features of one view which can be most faithfully inferred from the other. We find that these are also the most correlated variables in the sense of deep canonical correlation analysis (DCCA). Moreover, we show that these variables can be used to construct a non-parametric representation of the implied joint probability distribution, which can be thought of as a classical version of the Schmidt decomposition of quantum states. This representation can be used to compute the expectations of functions over one view of data conditioned on the other, such as Bayesian estimators and their standard deviations. We test the approach using inference on occluded MNIST images, and show that our representation contains multiple modes. Surprisingly, when applied to supervised learning (one dataset consists of labels), this approach automatically provides regularization and faster convergence compared to the cross-entropy objective. We also explore using this approach to discover salient independent variables of a single dataset.


title: WASP-169, WASP-171, WASP-175 and WASP-182: three hot Jupiters and one bloated sub-Saturn mass planet discovered by WASP-south

abstract: We present the discovery of four new giant planets from \WASP, three hot Jupiters and one bloated sub-Saturn mass planet; WASP-169b, WASP-171b, WASP-175b and WASP-182b. Besides the discovery photometry from \wasp\ we use radial velocity measurements from CORALIE and HARPS as well as follow-up photometry from EulerCam, TRAPPIST-North and -South and SPECULOOS. WASP-169b is a low density Jupiter ($M=0.561 \pm 0.061~\mathrm{M_{Jup}}, R=1.304^{+0.150}_{-0.073} ~\mathrm{R_{Jup}}$) orbiting a V=12.17 F8 sub-giant in a 5.611~day orbit. WASP-171b is a typical hot Jupiter ($M=1.084 \pm 0.094~\mathrm{M_{Jup}}, R=0.98^{+0.07}_{-0.04} ~\mathrm{R_{Jup}}$, $P=3.82~\mathrm{days}$) around a V=13.05 G0 star. We find a linear drift in the radial velocities of WASP-171 spanning 3.5 years, indicating the possibility of an additional outer planet or stellar companion. WASP-175b is an inflated hot Jupiter ($M=0.99 \pm 0.13~\mathrm{M_{Jup}}, R=1.208 \pm 0.081 ~\mathrm{R_{Jup}}$, $P=3.07~\mathrm{days}$) around a V=12.04 F7 star, which possibly is part of a binary system with a star 7.9\arcsec\ away. WASP-182b is a bloated sub-Saturn mass planet ($M=0.148 \pm 0.011~\mathrm{M_{Jup}}, R=0.850\pm 0.030 ~\mathrm{R_{Jup}}$) around a metal rich V=11.98 G5 star ([Fe/H]$=0.27 \pm 0.11$). With a orbital period of $P=3.377~\mathrm{days}$, it sits right in the apex of the sub-Jovian desert, bordering the upper- and lower edge of the desert in both the mass-period and radius-period plane. WASP-169b, WASP-175b and WASP-182b are promising targets for atmospheric characterisation through transmission spectroscopy, with expected transmission signals of 121, 150 and 264 ppm respectively.


title: Mean Field Approach for Configuring Population Dynamics on a Biohybrid Neuromorphic System

abstract: Real-time coupling of cell cultures to neuromorphic circuits necessitates a neuromorphic network that replicates biological behaviour both on a per-neuron and on a population basis, with a network size comparable to the culture. We present a large neuromorphic system composed of 9 chips, with overall 2880 neurons and 144M conductance-based synapses. As they are realized in a robust switched-capacitor fashion, individual neurons and synapses can be configured to replicate with high fidelity a wide range of biologically realistic behaviour. In contrast to other exploration/heuristics-based approaches, we employ a theory-guided mesoscopic approach to configure the overall network to a range of bursting behaviours, thus replicating the statistics of our targeted in-vitro network. The mesoscopic approach has implications beyond our proposed biohybrid, as it allows a targeted exploration of the behavioural space, which is a non-trivial task especially in large, recurrent networks.


title: Minimizing Perceived Image Quality Loss Through Adversarial Attack Scoping

abstract: Neural networks are now actively being used for computer vision tasks in security critical areas such as robotics, face recognition, autonomous vehicles yet their safety is under question after the discovery of adversarial attacks. In this paper we develop simplified adversarial attack algorithms based on a scoping idea, which enables execution of fast adversarial attacks that minimize structural image quality (SSIM) loss, allows performing efficient transfer attacks with low target inference network call count and opens a possibility of an attack using pen-only drawings on a paper for the MNIST handwritten digit dataset. The presented adversarial attack analysis and the idea of attack scoping can be easily expanded to different datasets, thus making the paper's results applicable to a wide range of practical tasks.


title: Momentum-resolved lattice dynamics of parent and electron-doped Sr$_{2}$IrO$_{4}$

abstract: The mixing of orbital and spin character in the wave functions of the $5d$ iridates has led to predictions of strong couplings among their lattice, electronic and magnetic degrees of freedom. As well as realizing a novel spin-orbit assisted Mott-insulating ground state, the perovskite iridate Sr$_{2}$IrO$_{4}$ has strong similarities with the cuprate La$_{2}$CuO$_{4}$, which on doping hosts a charge-density wave that appears intimately connected to high-temperature superconductivity. These phenomena can be sensitively probed through momentum-resolved measurements of the lattice dynamics, made possible by meV-resolution inelastic x-ray scattering. Here we report the first such measurements for both parent and electron-doped Sr$_{2}$IrO$_{4}$. We find that the low-energy phonon dispersions and intensities in both compounds are well described by the same nonmagnetic density functional theory calculation. In the parent compound, no changes of the phonons on magnetic ordering are discernible within the experimental resolution, and in the doped compound no anomalies are apparent due to charge-density waves. These measurements extend our knowledge of the lattice properties of (Sr$_{1-x}$La$_{x}$)$_{2}$IrO$_{4}$ and constrain the couplings of the phonons to magnetic and charge order.


title: Calibration of quantum sensors by neural networks

abstract: Introducing quantum sensors as solution to real-world problem demands reliability and controllability outside laboratory conditions. Producers and operators ought to be assumed to have limited resources ready available for calibration, and yet, they should be able to trust the devices. Neural networks are almost ubiquitous for similar tasks for classical sensors: here we show the applications of this technique to calibrating a quantum photonic sensor. This is based on a set of training data, collected only relying on the available probe states, hence reducing overheads. We found that covering finely the parameter space is key to achieve uncertainties close to their ultimate level. This technique has potential to become the standard approach to calibrate quantum sensors.


title: Estimations of Integrated Information Based on Algorithmic Complexity and Dynamic Querying

abstract: The concept of information has emerged as a language in its own right, bridging several disciplines that analyze natural phenomena and man-made systems. Integrated information has been introduced as a metric to quantify the amount of information generated by a system beyond the information generated by its elements. Yet, this intriguing notion comes with the price of being prohibitively expensive to calculate, since the calculations require an exponential number of sub-divisions of a system. Here we introduce a novel framework to connect algorithmic randomness and integrated information and a numerical method for estimating integrated information using a perturbation test rooted in algorithmic information dynamics. This method quantifies the change in program size of a system when subjected to a perturbation. The intuition behind is that if an object is random then random perturbations have little to no effect to what happens when a shorter program but when an object has the ability to move in both directions (towards or away from randomness) it will be shown to be better integrated as a measure of sophistication telling apart randomness and simplicity from structure. We show that an object with a high integrated information value is also more compressible, and is, therefore, more sensitive to perturbations. We find that such a perturbation test quantifying compression sensitivity provides a system with a means to extract explanations--causal accounts--of its own behaviour. Our technique can reduce the number of calculations to arrive at some bounds or estimations, as the algorithmic perturbation test guides an efficient search for estimating integrated information. Our work sets the stage for a systematic exploration of connections between algorithmic complexity and integrated information at the level of both theory and practice.


title: Cell damage and mitigation in Swiss albino mice: experiment and modelling

abstract: Chronic exposure to inorganic arsenic is a potential cause of carcinogenesis. It elicits its potential by generation of ROS, leading to DNA, protein and lipid damage. Therefore, the deleterious effect of arsenic can be mitigated by quenching ROS using antioxidants. There is a homology between the protein coding regions of mice and human. Effect of these alterations in human can be mimicked in mice. Therefore to understand the underlying mechanism of arsenic toxicity and its amelioration by black tea, studies have been conducted in mice model. Long term exposure to iAs leads to tumour growth, which has been found to be alleviated by black tea. Observations reveal that black tea has two salutary effects on the growth of tumour: the rate of growth of damaged cells was appreciably reduced and an early saturation of the level of damage is achieved. To take the experimental findings further, the experimental data have been modelled with simple dynamical equations. The curves obtained from \textit{in vivo} studies have been fitted with the data obtained from the model. The corresponding steady states and their stabilities are analyzed.


title: Absence of magnetic thermal conductivity in the quantum spin liquid candidate EtMe3Sb[Pd(dmit)2]2 -- revisited

abstract: We present the ultralow-temperature specific heat and thermal conductivity measurements on single crystals of triangular-lattice organic compound EtMe$_3$Sb[Pd(dmit)$_2$]$_2$, which has long been considered as a gapless quantum spin liquid candidate. In specific heat measurements, a finite linear term is observed, consistent with the previous work [S. Yamashita $et$ $al.$, Nat. Commun. {\bf 2}, 275 (2011)]. However, we do not observe a finite residual linear term in the thermal conductivity measurements, and the thermal conductivity does not change in a magnetic field of 6 Tesla. These results are in sharp contrast to previous thermal conductivity measurements on EtMe$_3$Sb[Pd(dmit)$_2$]$_2$ [M. Yamashita $et$ $al.$ Science {\bf 328}, 1246 (2010)], in which a huge residual linear term was observed and attributed to highly mobile gapless excitations, likely the spinons of a quantum spin liquid. In this context, the true ground state of EtMe$_3$Sb[Pd(dmit)$_2$]$_2$ has to be reconsidered.


title: Is coding a relevant metaphor for building AI? A commentary on "Is coding a relevant metaphor for the brain?", by Romain Brette

abstract: Brette contends that the neural coding metaphor is an invalid basis for theories of what the brain does. Here, we argue that it is an insufficient guide for building an artificial intelligence that learns to accomplish short- and long-term goals in a complex, changing environment.


title: Background Model of the CUPID-0 Experiment

abstract: CUPID-0 is the first large mass array of enriched Zn$^{82}$Se scintillating low temperature calorimeters, operated at LNGS since 2017. During its first scientific runs, CUPID-0 collected an exposure of 9.95 kg yr. Thanks to the excellent rejection of $\alpha$ particles, we attained the lowest background ever measured with thermal detectors in the energy region where we search for the signature of $^{82}$Se neutrinoless double beta decay. In this work we develop a model to reconstruct the CUPID-0 background over the whole energy range of experimental data. We identify the background sources exploiting their distinctive signatures and we assess their extremely low contribution (down to $\sim10^{-4}$ counts/(keV kg yr)) in the region of interest for $^{82}$Se neutrinoless double beta decay search. This result represents a crucial step towards the comprehension of the background in experiments based on scintillating calorimeters and in next generation projects such as CUPID.


title: Fallopian tube anatomy predicts pregnancy and pregnancy outcomes after tubal reversal surgery

abstract: We conducted this study to determine whether fallopian tube anatomy can predict the likelihood of pregnancy and pregnancy outcomes after tubal sterilization reversal. We built a flexible, non-parametric, multivariate model via generalized additive models to assess the effects of the following tubal parameters observed during tubal reparative surgery: tubal lengths; differences in tubal segment location, and diameters at the anastomosis sites; and, fibrosis of the tubal muscularis. In this study population, age and tubal length - in that order - were the primary factors predicting the likelihood of pregnancy. For pregnancy outcomes, tubal length was the most influential predictor of birth and ectopic pregnancy, while age was the primary predictor of miscarriage. Segment location and diameters contributed slightly to the odds of miscarriage and ectopic pregnancy. Tubal muscularis fibrosis had a little apparent effect. This study is the first to show that a statistical learning predictive model based on fallopian tube anatomy can predict pregnancy and pregnancy outcome probabilities after tubal reversal surgery.


title: Spike-Based Winner-Take-All Computation: Fundamental Limits and Order-Optimal Circuits

abstract: Winner-Take-All (WTA) refers to the neural operation that selects a (typically small) group of neurons from a large neuron pool. It is conjectured to underlie many of the brain's fundamental computational abilities. However, not much is known about the robustness of a spike-based WTA network to the inherent randomness of the input spike trains. In this work, we consider a spike-based $k$--WTA model wherein $n$ randomly generated input spike trains compete with each other based on their underlying statistics, and $k$ winners are supposed to be selected. We slot the time evenly with each time slot of length $1\, ms$, and model the $n$ input spike trains as $n$ independent Bernoulli processes. The Bernoulli process is a good approximation of the popular Poisson process but is more biologically relevant as it takes the refractory periods into account. Due to the randomness in the input spike trains, no circuits can guarantee to successfully select the correct winners in finite time. We focus on analytically characterizing the minimal amount of time needed so that a target minimax decision accuracy (success probability) can be reached. We first derive an information-theoretic lower bound on the decision time. We show that to have a (minimax) decision error $\le \delta$ (where $\delta \in (0,1)$), the computation time of any WTA circuit is at least \[ ((1-\delta) \log(k(n -k)+1) -1)T_{\mathcal{R}}, \] where $T_{\mathcal{R}}$ is a difficulty parameter of a WTA task that is independent of $\delta$, $n$, and $k$. We then design a simple WTA circuit whose decision time is \[ O( \log\frac{1}{\delta}+\log k(n-k))T_{\mathcal{R}}). \] It turns out that for any fixed $\delta \in (0,1)$, this decision time is order-optimal in terms of its scaling in $n$, $k$, and $T_{\mathcal{R}}$.


title: Efficient single input-output layer spiking neural classifier with time-varying weight model

abstract: This paper presents a supervised learning algorithm, namely, the Synaptic Efficacy Function with Meta-neuron based learning algorithm (SEF-M) for a spiking neural network with a time-varying weight model. For a given pattern, SEF-M uses the learning algorithm derived from meta-neuron based learning algorithm to determine the change in weights corresponding to each presynaptic spike times. The changes in weights modulate the amplitude of a Gaussian function centred at the same presynaptic spike times. The sum of amplitude modulated Gaussian functions represents the synaptic efficacy functions (or time-varying weight models). The performance of SEF-M is evaluated against state-of-the-art spiking neural network learning algorithms on 10 benchmark datasets from UCI machine learning repository. Performance studies show superior generalization ability of SEF-M. An ablation study on time-varying weight model is conducted using JAFFE dataset. The results of the ablation study indicate that using a time-varying weight model instead of single weight model improves the classification accuracy by 14%. Thus, it can be inferred that a single input-output layer spiking neural network with time-varying weight model is computationally more efficient than a multi-layer spiking neural network with long-term or short-term weight model.


title: Cowen-Douglas operators and the third of Halmos' ten problems

abstract: Let $T$ be a bounded linear operator on a complex separable infinite dimensional Hilbert space $\mathcal{H}$. $T$ is called intransitive if it leaves invariant spaces other than 0 or the whole space $\mathcal{H}$; otherwise it is transitive. In 1970, P. R. Halmos raised ten open problems on operator theory. In the past more than 50 years, nine of Halmos' ten problems were answered, but only the third one has made little progress. The third problem of Halmos is the following: if an intransitive operator has an inverse, is its inverse also intransitive? In this paper, we establish a set of theoretical systems with the help of Cowen-Douglas operators and spectral analysis. We give an affirmative answer to this problem under certain spectral conditions, which make essential progress in the research of Halmos' third problem. As the first application, we show that for an invertible hyponormal operator $T$, if $T^{-1}$ is intransitive and int$\sigma(T^{-1})^{\land}$ is not connected, then $T$ is also intransitive. As the second application, we show that if $T^{-1}$ has a proper strictly cyclic invariant subspace and there exists a bounded open set $\Omega$ which is a connected component of $\rho(T^{-1})$ such that $\Omega\cap \mathcal{U}_0=\emptyset$, where $\mathcal{U}_0$ is the connected component of $int(\sigma(T^{-1})^\land)$ containing zero point, then $T$ is intransitive.


title: Thermal conductivity of the quantum spin liquid candidate EtMe3Sb[Pd(dmit)2]2: No evidence of mobile gapless excitations

abstract: The thermal conductivity $\kappa$ of the quasi-2D organic spin-liquid candidate EtMe$_3$Sb[Pd(dmit)$_2$]$_2$ (dmit-131) was measured at low temperatures, down to 0.07 K. We observe a vanishingly small residual linear term $\kappa_0/T$, in $\kappa/T$ vs $T$ as $T \to 0$. This shows that the low-energy excitations responsible for the sizeable residual linear term $\gamma$ in the specific heat $C$, seen in $C/T$ vs $T$ as $T \to 0,$ are localized. We conclude that there are no mobile gapless excitations in this spin liquid candidate, in contrast with a prior study of dmit-131 that reported a large $\kappa_0/T$ value [Yamashita et al., Science 328, 1246 (2010)]. Our study shows that dmit-131 is in fact similar to $\kappa$-(BEDT-TTF)$_2$Cu$_2$(CN)$_3$, another quasi-2D organic spin-liquid candidate where a vanishingly small $\kappa_0/T$ and a sizeable $\gamma$ are seen. We attribute heat conduction in these organic insulators without magnetic order to phonons undergoing strong spin-phonon scattering, as observed in several other spin-liquid materials.


title: Optimizing Search API Queries for Twitter Topic Classifiers Using a Maximum Set Coverage Approach

abstract: Twitter has grown to become an important platform to access immediate information about major events and dynamic topics. As one example, recent work has shown that classifiers trained to detect topical content on Twitter can generalize well beyond the training data. Since access to Twitter data is hidden behind a limited search API, it is impossible (for most users) to apply these classifiers directly to the Twitter unfiltered data streams ("firehose"). Rather, applications must first decide what content to retrieve through the search API before filtering that content with topical classifiers. Thus, it is critically important to query the Twitter API relative to the intended topical classifier in a way that minimizes the amount of negatively classified data retrieved. In this paper, we propose a sequence of query optimization methods that generalize notions of the maximum coverage problem to find the subset of query terms within the API limits that cover most of the topically relevant tweets without sacrificing precision. We evaluate the proposed methods on a large dataset of Twitter data collected during 2013 and 2014 labeled using manually curated hashtags for eight topics. Among many insights, our analysis shows that the best of the proposed methods can significantly outperform the firehose on precision and F1-score while achieving high recall within strict API limitations.


title: From coupled-wire construction of quantum Hall states to wave functions and hydrodynamics

abstract: In this paper we use a close connection between the coupled wire construction (CWC) of Abelian quantum Hall states and the theory of composite bosons to extract the Laughlin wave function and the hydrodynamic effective theory in the bulk, including the Wen-Zee topological action, directly from the CWC. We show how rotational invariance can be recovered by fine-tuning the interactions. A simple recipe is also given to construct general Abelian quantum Hall states desceibed by the multi-component Wen-Zee action.


title: Big Math and the One-Brain Barrier A Position Paper and Architecture Proposal

abstract: Over the last decades, a class of important mathematical results have required an ever increasing amount of human effort to carry out. For some, the help of computers is now indispensable. We analyze the implications of this trend towards "big mathematics", its relation to human cognition, and how machine support for big math can be organized. The central contribution of this position paper is an information model for "doing mathematics", which posits that humans very efficiently integrate four aspects: inference, computation, tabulation, and narration around a well-organized core of mathematical knowledge. The challenge for mathematical software systems is that these four aspects need to be integrated as well. We briefly survey the state of the art.


title: Exponential Random Graph models for Little Networks

abstract: Statistical models for social networks have enabled researchers to study complex social phenomena that give rise to observed patterns of relationships among social actors and to gain a rich understanding of the interdependent nature of social ties and actors. Much of this research has focused on social networks within medium to large social groups. To date, these advances in statistical models for social networks, and in particular, of Exponential-Family Random Graph Models (ERGMS), have rarely been applied to the study of small networks, despite small network data in teams, families, and personal networks being common in many fields. In this paper, we revisit the estimation of ERGMs for small networks and propose using exhaustive enumeration when possible. We developed an R package that implements the estimation of pooled ERGMs for small networks using Maximum Likelihood Estimation (MLE), called "ergmito". Based on the results of an extensive simulation study to assess the properties of the MLE estimator, we conclude that there are several benefits of direct MLE estimation compared to approximate methods and that this creates opportunities for valuable methodological innovations that can be applied to modeling social networks with ERGMs.


title: Note on bimetric causal diagrams

abstract: A new kind of diagrams is presented, showing the causal structure of bimetric interactions.


title: Towards joint sound scene and polyphonic sound event recognition

abstract: Acoustic Scene Classification (ASC) and Sound Event Detection (SED) are two separate tasks in the field of computational sound scene analysis. In this work, we present a new dataset with both sound scene and sound event labels and use this to demonstrate a novel method for jointly classifying sound scenes and recognizing sound events. We show that by taking a joint approach, learning is more efficient and whilst improvements are still needed for sound event detection, SED results are robust in a dataset where the sample distribution is skewed towards sound scenes.


title: Genuine infinitesimal bendings of submanifolds

abstract: A basic question in submanifold theory is whether a given isometric immersion $f\colon M^n\to\R^{n+p}$ of a Riemannian manifold of dimension $n\geq 3$ into Euclidean space with low codimension $p$ admits, locally or globally, a genuine infinitesimal bending. That is, if there exists a genuine smooth variation of $f$ by immersions that are isometric up to the first order. Until now only the hypersurface case $p=1$ was well understood. We show that a strong necessary local condition to admit such a bending is the submanifold to be ruled and give a lower bound for the dimension of the rulings. In the global case, we describe the situation of compact submanifolds of dimension $n\geq 5$ in codimension $p=2$.


title: A Refined Scaling Law for Spatially Coupled LDPC Codes Over the Binary Erasure Channel

abstract: We propose a refined scaling law to predict the finite-length performance in the waterfall region of spatially coupled low-density parity-check codes over the binary erasure channel. In particular, we introduce some improvements to the scaling law proposed by Olmos and Urbanke that result in a better agreement between the predicted and simulated frame error rate. We also show how the scaling law can be extended to predict the bit error rate performance.


title: BIT: Biologically Inspired Tracker

abstract: Visual tracking is challenging due to image variations caused by various factors, such as object deformation, scale change, illumination change and occlusion. Given the superior tracking performance of human visual system (HVS), an ideal design of biologically inspired model is expected to improve computer visual tracking. This is however a difficult task due to the incomplete understanding of neurons' working mechanism in HVS. This paper aims to address this challenge based on the analysis of visual cognitive mechanism of the ventral stream in the visual cortex, which simulates shallow neurons (S1 units and C1 units) to extract low-level biologically inspired features for the target appearance and imitates an advanced learning mechanism (S2 units and C2 units) to combine generative and discriminative models for target location. In addition, fast Gabor approximation (FGA) and fast Fourier transform (FFT) are adopted for real-time learning and detection in this framework. Extensive experiments on large-scale benchmark datasets show that the proposed biologically inspired tracker performs favorably against state-of-the-art methods in terms of efficiency, accuracy, and robustness. The acceleration technique in particular ensures that BIT maintains a speed of approximately 45 frames per second.


title: Expanded Combinatorial Designs as Tool to Model Network Slicing in 5G

abstract: The network slice management function (NSMF) in 5G has a task to configure the network slice instances and to combine network slice subnet instances from the new-generation radio access network and the core network into an end-to-end network slice instance. In this paper, we propose a mathematical model for network slicing based on combinatorial designs such as Latin squares and rectangles and their conjugate forms. We extend those designs with attributes that offer different levels of abstraction. For one set of attributes we prove a stability Lemma for the necessary conditions to reach a stationary ergodic stage. We also introduce a definition of utilization ratio function and offer an algorithm for its maximization. Moreover, we provide algorithms that simulate the work of NSMF with randomized or optimized strategies, and we report the results of our implementation, experiments and simulations for one set of attributes.


title: Metallic state in a strongly interacting spinless two-valley electron system in two dimensions

abstract: We have studied the strongly interacting, two-valley two-dimensional (2D) electron system in ultrahigh mobility SiGe/Si/SiGe quantum wells in parallel magnetic fields strong enough to completely polarize the electron spins thus making the electron system "spinless". It occurs that the metallic temperature dependence of the resistivity, although weaker than that in the absence of magnetic field, still remains strong even when the spin degree of freedom is removed. Several independent methods have been used to establish the existence of the genuine MIT in the spinless two-valley 2D system. This is in contrast to the previous results obtained on more disordered silicon samples, where the polarizing magnetic field causes a complete quench of the metallic temperature behavior.


title: The Theorem Prover Museum -- Conserving the System Heritage of Automated Reasoning

abstract: We present the Theorem Prover Museum, and initiative to conserve -- and make publicly available -- the sources and source-related artefacts of automated reasoning systems. Theorem provers have been at the forefront of Artificial Intelligence, stretching the limits of computation, and incubating many innovations we take for granted today. Without the systems themselves as preserved cultural artefacts, future historians will have difficulties to study the history of science and engineering in our discipline.


title: Disentangle contributions to small-system collectivity via scans of light nucleus-nucleus collisions

abstract: The observation of multi-particle azimuthal correlations in high-energy small-system collisions has led to intense debate on its origin and the possible coexistence from two competing theoretical scenarios: one based on initial-state intrinsic momentum anisotropy (ISM), and the other based on final-state collective response to the collision geometry (FSM). To complement the previous scan of asymmetric collision systems ($p$+Au, $d$+Au and He+Au), we propose a scan of small symmetric collision systems at RHIC, such as C+C, O+O, Al+Al and Ar+Ar $\sqrt{s_{\mathrm{NN}}}=0.2$ TeV, to provide further insights in disentangling contributions from these two scenarios. These symmetric small systems have the advantage of providing a better controlled initial geometry dominated by the average shape of the overlap region, as opposed to fluctuation-driven geometries in asymmetric systems. A transport model is employed to investigate the expected geometry response in the FSM scenario. Different trends of elliptic flow with increasing charge particle multiplicity are observed between symmetric and asymmetric systems, while triangular flow appears to show a similar behavior. Furthermore, a comparison of O+O collisions at $\sqrt{s_{\mathrm{NN}}}=0.2$ TeV and at $\sqrt{s_{\mathrm{NN}}}=2.76-7$ TeV, as proposed at the LHC, provides a unique opportunity to disentangle the collision geometry effects at nucleon level from those arising from subnucleon fluctuations.


title: Regression-Enhanced Random Forests

abstract: Random forest (RF) methodology is one of the most popular machine learning techniques for prediction problems. In this article, we discuss some cases where random forests may suffer and propose a novel generalized RF method, namely regression-enhanced random forests (RERFs), that can improve on RFs by borrowing the strength of penalized parametric regression. The algorithm for constructing RERFs and selecting its tuning parameters is described. Both simulation study and real data examples show that RERFs have better predictive performance than RFs in important situations often encountered in practice. Moreover, RERFs may incorporate known relationships between the response and the predictors, and may give reliable predictions in extrapolation problems where predictions are required at points out of the domain of the training dataset. Strategies analogous to those described here can be used to improve other machine learning methods via combination with penalized parametric regression techniques.


title: Network Investment Game with Wardrop Followers

abstract: We study a two-sided network investment game consisting of two sets of players, called providers and users. The game is set in two stages. In the first stage, providers aim to maximize their profit by investing in bandwidth of cloud computing services. The investments of the providers yield a set of usable services for the users. In the second stage, each user wants to process a task and therefore selects a bundle of services so as to minimize the total processing time. We assume the total processing time to be separable over the chosen services and the processing time of each service to depend on the utilization of the service and the installed bandwidth. We provide insights on how competition between providers affects the total costs of the users and show that every game on a series-parallel graph can be reduced to an equivalent single edge game when analyzing the set of subgame perfect Nash equilibria.


title: A model independent parametrization of the late time cosmic acceleration: constraints on the parameters from recent observations

abstract: In this work, we have considered a model independent approach to study the nature of the late time cosmic acceleration. We have used the Pade approximation to parametrize the comoving distance. Consequently, from this comoving distance, we derive a parameterization for the Hubble parameter. Our parametrization is completely analytic and valid for late-time and matter dominated eras only. This parametrization possesses sub-percentage accuracy compared to any arbitrary cosmological model or parametrization up to matter dominated era. Using this parametrization, we put constraints on the parameters from recent low redshift cosmological observations including Planck 2018 distance priors. Our results show that the $\Lambda$CDM model is 1$\sigma$ to 2$\sigma$ away for lower redshifts. We find that the phantom crossing is allowed by all the combinations of dataset considered. We also find that the dynamical dark energy models are preferable at lower redshifts. Our study also shows that, at lower redshifts ($z<0.5$), phantom models are allowed at almost 1$\sigma$ confidence level.


title: GumDrop at the DISRPT2019 Shared Task: A Model Stacking Approach to Discourse Unit Segmentation and Connective Detection

abstract: In this paper we present GumDrop, Georgetown University's entry at the DISRPT 2019 Shared Task on automatic discourse unit segmentation and connective detection. Our approach relies on model stacking, creating a heterogeneous ensemble of classifiers, which feed into a metalearner for each final task. The system encompasses three trainable component stacks: one for sentence splitting, one for discourse unit segmentation and one for connective detection. The flexibility of each ensemble allows the system to generalize well to datasets of different sizes and with varying levels of homogeneity.


title: Projection bands and atoms in pervasive pre-Riesz spaces

abstract: In vector lattices, the concept of a projection band is a basic tool. We deal with projection bands in the more general setting of an Archimedean pre-Riesz space $X$. We relate them to projection bands in a vector lattice cover $Y$ of $X$. If $X$ is pervasive, then a projection band in $X$ extends to a projection band in $Y$, whereas the restriction of a projection band $B$ in $Y$ is not a projection band in $X$, in general. We give conditions under which the restriction of $B$ is a projection band in $X$. We introduce atoms and discrete elements in $X$ and show that every atom is discrete. The converse implication is true, provided $X$ is pervasive. In this setting, we link atoms in $X$ to atoms in $Y$. If $X$ contains an atom $a>0$, we show that the principal band generated by $a$ is a projection band. Using atoms in a finite dimensional Archimedean pre-Riesz space $X$, we establish that $X$ is pervasive if and only if it is a vector lattice.


title: EvryFlare I: Long-term Evryscope Monitoring of Flares from the Cool Stars Across Half the Southern Sky

abstract: We search for superflares from 4,068 cool stars in 2+ years of Evryscope photometry, focusing on those with high-cadence data from both Evryscope and TESS. The Evryscope array of small telescopes observed 575 flares from 284 stars, with a median energy of 10^34.0 erg. Since 2016, Evryscope has enabled the detection of rare events from all stars observed by TESS through multi-year, high-cadence continuous observing. We report ~2X the previous largest number of 10^34 erg high-cadence flares from nearby cool stars. We find 8 flares with amplitudes of 3+ g' magnitudes, with the largest reaching 5.6 magnitudes and releasing 10^36.2 erg. We observe a 10^34 erg superflare from TOI-455 (LTT 1445), a mid-M with a rocky planet candidate. We measure the superflare rate per flare-star and quantify the average flaring of active stars as a function of spectral type, including superflare rates, FFDs, and typical flare amplitudes in g'. We confirm superflare morphology is broadly consistent with magnetic re-connection. We estimate starspot coverage necessary to produce superflares, and hypothesize maximum-allowed superflare energies and waiting-times between flares corresponding to 100% coverage of the stellar hemisphere. We observe decreased flaring at high galactic latitudes. We explore the effects of superflares on ozone loss to planetary atmospheres: we observe 1 superflare with sufficient energy to photo-dissociate all ozone in an Earth-like atmosphere in one event. We find 17 stars that may deplete an Earth-like atmosphere via repeated flaring. Of the 1822 stars around which TESS may discover temperate rocky planets, we observe 14.6% +/- 2% emit large flares.


title: Unlocking CO Depletion in Protoplanetary Disks II. Primordial C/H Predictions Inside the CO Snowline

abstract: CO is thought to be the main reservoir of volatile carbon in protoplanetary disks, and thus the primary initial source of carbon in the atmospheres of forming giant planets. However, recent observations of protoplanetary disks point towards low volatile carbon abundances in many systems, including at radii interior to the CO snowline. One potential explanation is that gas phase carbon is chemically reprocessed into less volatile species, which are frozen on dust grain surfaces as ice. This mechanism has the potential to change the primordial C/H ratio in the gas. However, current observations primarily probe the upper layers of the disk. It is not clear if the low volatile carbon abundances extend to the midplane, where planets form. We have run a grid of 198 chemical models, exploring how the chemical reprocessing of CO depends on disk mass, dust grain size distribution, temperature, cosmic ray and X-ray ionization rate, and initial water abundance. Building on our previous work focusing on the warm molecular layer, here we analyze the results for our grid of models in the disk midplane at 12 au. We find that either an ISM level cosmic ray ionization rate or the presence of UV photons due to a low dust surface density are needed to chemically reduce the midplane CO gas abundance by at least an order of magnitude within 1 Myr. In the majority of our models CO does not undergo substantial reprocessing by in situ chemistry and there is little change in the gas phase C/H and C/O ratios over the lifetime of the typical disk. However, in the small sub-set of disks where the disk midplane is subject to a source of ionization or photolysis, the gas phase C/O ratio increases by up to nearly 9 orders of magnitude due to conversion of CO into volatile hydrocarbons.


title: A Data Ecosystem to Support Machine Learning in Materials Science

abstract: Facilitating the application of machine learning to materials science problems will require enhancing the data ecosystem to enable discovery and collection of data from many sources, automated dissemination of new data across the ecosystem, and the connecting of data with materials-specific machine learning models. Here, we present two projects, the Materials Data Facility (MDF) and the Data and Learning Hub for Science (DLHub), that address these needs. We use examples to show how MDF and DLHub capabilities can be leveraged to link data with machine learning models and how users can access those capabilities through web and programmatic interfaces.


title: Interpretable and Generalizable Person Re-Identification with Query-Adaptive Convolution and Temporal Lifting

abstract: For person re-identification, existing deep networks often focus on representation learning. However, without transfer learning, the learned model is fixed as is, which is not adaptable for handling various unseen scenarios. In this paper, beyond representation learning, we consider how to formulate person image matching directly in deep feature maps. We treat image matching as finding local correspondences in feature maps, and construct query-adaptive convolution kernels on the fly to achieve local matching. In this way, the matching process and results are interpretable, and this explicit matching is more generalizable than representation features to unseen scenarios, such as unknown misalignments, pose or viewpoint changes. To facilitate end-to-end training of this architecture, we further build a class memory module to cache feature maps of the most recent samples of each class, so as to compute image matching losses for metric learning. Through direct cross-dataset evaluation, the proposed Query-Adaptive Convolution (QAConv) method gains large improvements over popular learning methods (about 10%+ mAP), and achieves comparable results to many transfer learning methods. Besides, a model-free temporal cooccurrence based score weighting method called TLift is proposed, which improves the performance to a further extent, achieving state-of-the-art results in cross-dataset person re-identification. Code is available at https://github.com/ShengcaiLiao/QAConv.


title: On the counting problem in inverse Littlewood--Offord theory

abstract: Let $\epsilon_1, \dotsc, \epsilon_n$ be i.i.d. Rademacher random variables taking values $\pm 1$ with probability $1/2$ each. Given an integer vector $\boldsymbol{a} = (a_1, \dotsc, a_n)$, its concentration probability is the quantity $\rho(\boldsymbol{a}):=\sup_{x\in \mathbb{Z}}\Pr(\epsilon_1 a_1+\dots+\epsilon_n a_n = x)$. The Littlewood-Offord problem asks for bounds on $\rho(\boldsymbol{a})$ under various hypotheses on $\boldsymbol{a}$, whereas the inverse Littlewood-Offord problem, posed by Tao and Vu, asks for a characterization of all vectors $\boldsymbol{a}$ for which $\rho(\boldsymbol{a})$ is large. In this paper, we study the associated counting problem: How many integer vectors $\boldsymbol{a}$ belonging to a specified set have large $\rho(\boldsymbol{a})$? The motivation for our study is that in typical applications, the inverse Littlewood-Offord theorems are only used to obtain such counting estimates. Using a more direct approach, we obtain significantly better bounds for this problem than those obtained using the inverse Littlewood--Offord theorems of Tao and Vu and of Nguyen and Vu. Moreover, we develop a framework for deriving upper bounds on the probability of singularity of random discrete matrices that utilizes our counting result. To illustrate the methods, we present the first `exponential-type' (i.e., $\exp(-n^c)$ for some positive constant $c$) upper bounds on the singularity probability for the following two models: (i) adjacency matrices of dense signed random regular digraphs, for which the previous best known bound is $O(n^{-1/4})$ due to Cook; and (ii) dense row-regular $\{0,1\}$-matrices, for which the previous best known bound is $O_{C}(n^{-C})$ for any constant $C>0$ due to Nguyen.


title: Power Consumption and Energy-Efficiency for In-Band Full-Duplex Wireless Systems

abstract: This paper presents an analytical model of power consumption for In-Band Full-Duplex (IBFD) Wireless Local-Area Networks (WLANs). Energy-efficiency is compared for both Half-Duplex (HD) and IBFD networks. The presented analytical model closely matches the results generated by simulation. For a given traffic scenario, IBFD systems exhibit higher power consumption, however at improved energy efficiency when compared to equivalent HD WLANs.


title: A Convex Body Associated to the Busemann Random Simplex Inequality and the Petty conjecture

abstract: Given $L$ a convex body, the $L_p$-Busemann Random Simplex Inequality is closely related to the centroid body $\Gamma_p L$ for $p=1$ and $2$, and only in these cases it can be proved using the $L_p$-Busemann-Petty centroid inequality. We define a convex body $N_p L$ and prove an isoperimetric inequality for $(N_p L)^\circ$ that is equivalent to the $L_p$-Busemann Random Simplex Inequality. As applications, we give a simple proof of a general functional version of the Busemann Random Simplex Inequality and study a dual theory related to Petty's conjectured inequality. More precisely, we prove dual versions of the $L_p$-Busemann Random Simplex Inequality for sets and functions by means of the $p$-affine surface area measure, and we prove that the Petty conjecture is equivalent to an $L_1$-Sharp Affine Sobolev-type inequality that is stronger than (and directly implies) the Sobolev-Zhang inequality.


title: On the Kullback-Leibler divergence between location-scale densities

abstract: We show that the $f$-divergence between any two densities of potentially different location-scale families can be reduced to the calculation of the $f$-divergence between one standard density with another location-scale density. It follows that the $f$-divergence between two scale densities depends only on the scale ratio. We then report conditions on the standard distribution to get symmetric $f$-divergences: First, we prove that all $f$-divergences between densities of a location family are symmetric whenever the standard density is even, and second, we illustrate a generic symmetric property with the calculation of the Kullback-Leibler divergence between scale Cauchy distributions. Finally, we show that the minimum $f$-divergence of any query density of a location-scale family to another location-scale family is independent of the query location-scale parameters.


title: DenseNet Models for Tiny ImageNet Classification

abstract: In this paper, we present two image classification models on the Tiny ImageNet dataset. We built two very different networks from scratch based on the idea of Densely Connected Convolution Networks. The architecture of the networks is designed based on the image resolution of this specific dataset and by calculating the Receptive Field of the convolution layers. We also used some non-conventional techniques related to image augmentation and Cyclical Learning Rate to improve the accuracy of our models. The networks are trained under high constraints and low computation resources. We aimed to achieve top-1 validation accuracy of 60%; the results and error analysis are also presented.


title: Hopfological algebra for infinite dimensional Hopf algebras

abstract: We consider "Hopfological" techniques as in \cite{Ko} but for infinite dimensional Hopf algebras, under the assumption of being co-Frobenius. In particular, $H=k[{\mathbb Z}]\#k[x]/x^2$ is the first example, whose corepresentations category is d.g. vector spaces. Motivated by this example we define the "Homology functor" (we prove it is homological) for any co-Frobenius algebra, with coefficients in $H$-comodules, that recover usual homology of a complex when $H=k[{\mathbb Z}]\#k[x]/x^2$. Another easy example of co-Frobenius Hopf algebra gives the category of "mixed complexes" and we see (by computing an example) that this homology theory differs from cyclic homology, although there exists a long exact sequence analogous to the SBI-sequence. Finally, because we have a tensor triangulated category, its $K_0$ is a ring, and we prove a "last part of a localization exact sequence" for $K_0$ that allows us to compute -or describe- $K_0$ of some family of examples, giving light of what kind of rings can be categorified using this techniques.


title: A bound for the conductor of an open subgroup of GL2 associated to an elliptic curve

abstract: Given an elliptic curve $E$ without complex multiplication defined over a number field $K$, consider the image of the Galois representation defined by letting Galois act on the torsion of $E$. Serre's open image theorem implies that there is a positive integer $m$ for which the Galois image is completely determined by its reduction modulo $m$. In this note, we prove a bound on the smallest such $m$ in terms of standard invariants associated with $E$. The bound is sharp and improves upon previous results.


title: B-meson charged current anomalies: the post-Moriond status

abstract: In this note, we discuss the impact of the recent Belle result on the various theoretical explanations of the $R_D$ and $R_{D^*}$ anomalies. The pure tensor explanation, which was strongly disfavoured by the measurements of $F_L^{D^*}$ and high-$p_T$ $p \, p \to \tau \, \nu$ searches before Moriond, is now completely allowed because of reduction of the experimental world-average. Moreover, the pure right-chiral vector solution (involving right-chiral neutrinos) has now moved into the $2\sigma$ allowed range of the LHC $p \, p \to \tau \, \nu$ searches. We also critically re-examine the bound on $\mathcal{B}(B_c^- \to \tau^- \bar{\nu}_\tau)$ from LEP data and show that the bound is considerably weaker than the number $10\%$ often used in the recent literature.


title: Pulsating in unison at optical and X-ray energies: simultaneous high-time resolution observations of the transitional millisecond pulsar PSR J1023+0038

abstract: PSR J1023+0038 is the first millisecond pulsar discovered to pulsate in the visible band; such a detection took place when the pulsar was surrounded by an accretion disk and also showed X-ray pulsations. We report on the first high time resolution observational campaign of this transitional pulsar in the disk state, using simultaneous observations in the optical (TNG, NOT, TJO), X-ray (XMM-Newton, NuSTAR, NICER), infrared (GTC) and UV (Swift) bands. Optical and X-ray pulsations were detected simultaneously in the X-ray high intensity mode in which the source spends $\sim$ 70% of the time, and both disappeared in the low mode, indicating a common underlying physical mechanism. In addition, optical and X-ray pulses were emitted within a few km, had similar pulse shape and distribution of the pulsed flux density compatible with a power-law relation $F_{\nu} \propto \nu^{-0.7}$ connecting the optical and the 0.3-45 keV X-ray band. Optical pulses were detected also during flares with a pulsed flux reduced by one third with respect to the high mode; the lack of a simultaneous detection of X-ray pulses is compatible with the lower photon statistics. We show that magnetically channeled accretion of plasma onto the surface of the neutron star cannot account for the optical pulsed luminosity ($\sim 10^{31}$ erg/s). On the other hand, magnetospheric rotation-powered pulsar emission would require an extremely efficient conversion of spin-down power into pulsed optical and X-ray emission. We then propose that optical and X-ray pulses are instead produced by synchrotron emission from the intrabinary shock that forms where a striped pulsar wind meets the accretion disk, within a few light cylinder radii away, $\sim$ 100 km, from the pulsar.


title: Data-driven Computing in Elasticity via Chebyshev Approximation

abstract: This paper proposes a data-driven approach for computing elasticity by means of a non-parametric regression approach rather than an optimization approach. The Chebyshev approximation is utilized for tackling the material data-sets non-linearity of the elasticity. Also, additional efforts have been taken to compare the results with several other state-of-the-art methodologies.


title: Guaranteed and robust $L_2$-norm a posteriori error estimates for 1D linear advection problems

abstract: We propose a reconstruction-based a posteriori error estimate for linear advection problems in one space dimension. In our framework, a stable variational ultra-weak formulation is adopted, and the equivalence of the $L_2$-norm of the error with the dual graph norm of the residual is established. This dual norm is showed to be localizable over vertex-based patch subdomains of the computational domain under the condition of the orthogonality of the residual to the piecewise affine hat functions. We show that this condition is valid for some well-known numerical methods including continuous/discontinuous Petrov--Galerkin and discontinuous Galerkin methods. Consequently, a well-posed local problem on each patch is identified, which leads to a global conforming reconstruction of the discrete solution. We prove that this reconstruction provides a guaranteed upper bound on the $L_2$ error. Moreover, up to a constant, it also gives local lower bounds on the $L_2$ error, where the generic constant is proven to be independent of mesh-refinement, polynomial degree of the approximation, and the advective velocity. This leads to robustness of our estimates with respect to the advection as well as the polynomial degree. All the above properties are verified in a series of numerical experiments, additionally leading to asymptotic exactness. Motivated by these results, we finally propose a heuristic extension of our methodology to any space dimension, achieved by solving local least-squares problems on vertex-based patches. Though not anymore guaranteed, the resulting error indicator is numerically robust with respect to both advection velocity and polynomial degree, for a collection of two-dimensional test cases including discontinuous solutions.


title: Observational constraints of a new unified dark fluid and the $H_0$ tension

abstract: Unified cosmological models have received a lot of attention in astrophysics community for explaining both the dark matter and dark energy evolution. The Chaplygin cosmologies, a well known name in this group have been investigated matched with observations from different sources. Obviously, Chaplygin cosmologies have to obey restrictions in order to be consistent with the observational data. As a consequence, alternative unified models, differing from Chaplygin model, are of special interest. In the present work we consider a specific example of such a unified cosmological model, that is quantified by only a single parameter $\mu$, that can be considered as a minimal extension of the $\Lambda$-cold dark matter cosmology. We investigate its observational boundaries together with an analysis of the universe at large scale. Our study shows that at early time the model behaves like a dust, and as time evolves, it mimics a dark energy fluid depicting a clear transition from the early decelerating phase to the late cosmic accelerating phase. Finally, the model approaches the cosmological constant boundary in an asymptotic manner. We remark that for the present unified model, the estimations of $H_0$ are slightly higher than its local estimation and thus alleviating the $H_0$ tension.


title: $\alpha$-Logarithmic negativity

abstract: The logarithmic negativity of a bipartite quantum state is a widely employed entanglement measure in quantum information theory, due to the fact that it is easy to compute and serves as an upper bound on distillable entanglement. More recently, the $\kappa$-entanglement of a bipartite state was shown to be the first entanglement measure that is both easily computable and has a precise information-theoretic meaning, being equal to the exact entanglement cost of a bipartite quantum state when the free operations are those that completely preserve the positivity of the partial transpose [Wang and Wilde, Phys. Rev. Lett. 125(4):040502, July 2020]. In this paper, we provide a non-trivial link between these two entanglement measures, by showing that they are the extremes of an ordered family of $\alpha$-logarithmic negativity entanglement measures, each of which is identified by a parameter $\alpha\in[ 1,\infty] $. In this family, the original logarithmic negativity is recovered as the smallest with $\alpha=1$, and the $\kappa$-entanglement is recovered as the largest with $\alpha=\infty$. We prove that the $\alpha $-logarithmic negativity satisfies the following properties: entanglement monotone, normalization, faithfulness, and subadditivity. We also prove that it is neither convex nor monogamous. Finally, we define the $\alpha$-logarithmic negativity of a quantum channel as a generalization of the notion for quantum states, and we show how to generalize many of the concepts to arbitrary resource theories.

categories: quant-ph cond-mat.stat-mech cs.IT hep-th math-ph math.IT math.MP


title: Mini-survey of the northern sky to Dec <+30

abstract: We propose an extension of the LSST survey to cover the northern sky to DEC < +30 (accessible at airmass <1.8). This survey will increase the LSST sky coverage by ~9,600 square degrees from 18,900 to 28,500 square degrees (a 50% increase) but use only 0.6-2.5% of the time depending on the synergies with other surveys. This increased area addresses a wide range of science cases that enhance all of the primary LSST science goals by significant amounts. The science enabled includes: increasing the area of the sky accessible for follow-up of multi-messenger transients including gravitational waves, mapping the milky way halo and halo dwarfs including discovery of RR Lyrae stars in the outer galactic halo, discovery of z>7 quasars in combination Euclid, enabling a second generation DESI and other spectroscopic surveys, and enhancing all areas of science by improving synergies with Euclid, WFIRST, and unique northern survey facilities. This white paper is the result of the Tri-Agency Working Group (TAG) appointed to develop synergies between missions and presents a unified plan for northern coverage. The range of time estimates reflects synergies with other surveys. If the modified DESC WFD survey, the ecliptic plane mini survey, and the north galactic spur mini survey are executed this plan would only need 0.6% of the LSST time, however if none of these are included the overall request is 2.5% of the 10 year survey life. In other words, the majority of these observations are already suggested as part of these other surveys and the intent of this white paper is to propose a unified baseline plan to carry out a broad range of objectives to facilitate a combination of multiple science objectives. A companion white paper gives Euclid specific science goals, and we support the white papers for southern extensions of the LSST survey.


title: Enhancing LSST Science with Euclid Synergy

abstract: This white paper is the result of the Tri-Agency Working Group (TAG) appointed to develop synergies between missions and is intended to clarify what LSST observations are needed in order to maximally enhance the combined science output of LSST and Euclid. To facilitate LSST planning we provide a range of possible LSST surveys with clear metrics based on the improvement in the Dark Energy figure of merit (FOM). To provide a quantifiable metric we present five survey options using only between 0.3 and 3.8% of the LSST 10 year survey. We also provide information so that the LSST DDF cadence can possibly be matched to those of \emph{Euclid} in common deep fields, SXDS, COSMOS, CDFS, and a proposed new LSST deep field (near the Akari Deep Field South). Co-coordination of observations from the Large Synoptic Survey Telescope (LSST) and Euclid will lead to a significant number of synergies. The combination of optical multi-band imaging from LSST with high resolution optical and near-infrared photometry and spectroscopy from \emph{Euclid} will not only improve constraints on Dark Energy, but provide a wealth of science on the Milky Way, local group, local large scale structure, and even on first galaxies during the epoch of reionization. A detailed paper has been published on the Dark Energy science case (Rhodes et al.) by a joint LSST/Euclid working group as well as a white paper describing LSST/Euclid/WFIRST synergies (Jain et al.), and we will briefly describe other science cases here. A companion white paper argues the general science case for an extension of the LSST footprint to the north at airmass < 1.8, and we support the white papers for southern extensions of the LSST survey.


title: Minimally-Entangled State Preparation of Localized Wavefunctions on Quantum Computers

abstract: Initializing a single site of a lattice scalar field theory into an arbitrary state with support throughout the quantum register requires ${\cal O}(2^{n_Q})$ entangling gates on a quantum computer with $n_Q$ qubits per site. It is conceivable that, instead, initializing to functions that are good approximations to states may have utility in reducing the number of required entangling gates. In the case of a single site of a non-interacting scalar field theory, initializing to a symmetric exponential wavefunction requires $n_Q-1$ entangling gates, compared with the $2^{n_Q-1} + n_Q-3 + \delta_{n_Q,1}$ required for a symmetric Gaussian wavefunction. In this work, we explore the initialization of 1-site ($n_Q=4$), 2-site ($n_Q=3$) and 3-site ($n_Q=3$) non-interacting scalar field theories with symmetric exponential wavefunctions using IBM's quantum simulators and quantum devices (Poughkeepsie and Tokyo). With the digitizations attainable with $n_Q = 3,4$, these tensor-product wavefunctions are found to have large overlap with a Gaussian wavefunction, and provide a suitable low-noise initialization for improvement and \emph{Somma Inflation}. In performing these simulations, we have employed a workflow that interleaves calibrations to mitigate systematic errors in production. The calibrations allow tolerance cuts on gate performance including the fidelity of the symmetrizing Hadamard gate, both in vacuum ($|{\bf 0}\rangle^{\otimes n_Q}$) and in medium ($n_Q-1$ qubits initialized to an exponential function). The results obtained in this work are relevant to systems beyond scalar field theories, such as the deuteron radial wavefunction, 2- and 3-dimensional cartesian-space wavefunctions, and non-relativistic multi-nucleon systems built on a localized eigenbasis.


title: Characterizing the nonlinear structure of shared variability in cortical neuron populations using latent variable models

abstract: Sensory neurons often have variable responses to repeated presentations of the same stimulus, which can significantly degrade the stimulus information contained in those responses. This information can in principle be preserved if variability is shared across many neurons, but depends on the structure of the shared variability and its relationship to sensory encoding at the population level. The structure of this shared variability in neural activity can be characterized by latent variable models, although they have thus far typically been used under restrictive mathematical assumptions. Here we introduce two nonlinear latent variable models for analyzing large-scale neural recordings. We first present a general nonlinear latent variable model that is agnostic to the stimulus tuning properties of the individual neurons, and is hence well suited for exploring neural populations whose tuning properties are not well characterized. This motivates a second class of model, the Generalized Affine Model, which simultaneously determines each neuron's stimulus selectivity and a set of latent variables that modulate these stimulus-driven responses both additively and multiplicatively. While these approaches can detect very general nonlinear relationships in shared neural variability, we find that neural activity recorded in anesthetized primary visual cortex (V1) is best described by a single additive and single multiplicative latent variable, i.e. an `affine model'. In contrast, application of the same models to recordings in awake macaque prefrontal cortex discover more general nonlinearities to compactly describe the population response variability. These results thus demonstrate how nonlinear latent variable models can be used to describe population response variability, and suggest that a range of methods is necessary to study different brain regions under different experimental conditions.


title: Saddlepoint Approximations for Short-Packet Wireless Communications

abstract: In recent years, the derivation of nonasymptotic converse and achievability bounds on the maximum coding rate as a function of the error probability and blocklength has gained attention in the information theory literature. While these bounds are accurate for many scenarios of interest, they need to be evaluated numerically for most wireless channels of practical interest, and their evaluation is computationally demanding. This paper presents saddlepoint approximations of state-of-the-art converse and achievability bounds for noncoherent, single-antenna, Rayleigh block-fading channels. These approximations can be calculated efficiently and are shown to be accurate for SNR values as small as 0 dB and blocklengths of 168 channel uses or more.


title: A toxicology-informed, safer by design approach for the fabrication of transparent electrodes based on silver nanowires

abstract: Fabrication of silver nanowires (AgNWs) with fine and independent control of both the diameter (from 30 to 120 nm) and length (from 5 to 120 $\mu$m) by concomitant addition of co-nucleants and temperature control is demonstrated, and used for the preparation of size standards. Percolating random networks were fabricated using these standards and their optoelectronic properties were measured and compared with regard to the nanowire dimensions. The transparent electrodes appear suitable for various applications and exhibit excellent performances (e.g. 16 ohm sq --1 at 93% transparency), with haze values varying from 1.6 to 26.2%. Besides, in vitro toxicological studies carried out on murine macrophages with the same size standards revealed that AgNWs are weakly toxic (no toxicity observed below 50 $\mu$g mL --1 Ag), in particular compared to other silver nanoparticles. Short AgNWs (4 $\mu$m) appeared to be slightly more toxic than longer AgNWs (10 and 20 $\mu$m). Conversely, long AgNWs (20 $\mu$m) induced a more prolonged pro-inflammatory response in murine macrophages. These results contribute, in a safer by design approach, to promoting the use of short AgNWs. The global knowledge dealing with the combination of nanowire dimensions associated with optoelectronic performances and related toxicity should encourage the rational use of AgNWs, and guide the choice of the most adequate AgNW dimensions in an integrated approach.


title: Age of Information for Updates with Distortion

abstract: We consider an information update system where an information receiver requests updates from an information provider in order to minimize its age of information. The updates are generated at the transmitter as a result of completing a set of tasks such as collecting data and performing computations. We refer to this as the update generation process. We model the $quality$ (i.e., $distortion$) of an update as an increasing (resp. decreasing) function of the processing time spent while generating the update at the transmitter. While processing longer at the transmitter results in a better quality (lower distortion) update, it causes the update to age. We determine the age-optimal policies for the update request times at the receiver and update processing times at the transmitter subject to a minimum required quality (maximum allowed distortion) constraint on the updates.


title: Ultrasensitive all-optical thermometry using nanodiamonds with high concentration of silicon-vacancy centres and multiparametric data analysis

abstract: Nanoscale thermometry is paramount to study primary processes of heat transfer in solids and is a subject of hot debate in cell biology. Here we report ultrafast temperature sensing using all-optical thermometry exploiting synthetic nanodiamonds with silicon-vacancy (SiV) centres embedded at a high concentration. Using multi-parametric analysis of photoluminescence (PL) of these centres, we have achieved an intrinsic noise floor of about 10 mK Hz$^{-1/2}$, which is a thousand-fold increase in the readout speed in comparison to the current record values demonstrated with all-optical methods of comparable spatial-resolution and precision. Our thermometers are smaller than 250-nm across but can detect a 0.4$^\circ$C change of temperature in a measurement taking only 0.001 second. The exceptional sensitivity and simplicity of these thermometers enable a wide range of applications such as temperature monitoring and mapping within intracellular regions and in state-of-the-art solid-state electronic nanodevices.


title: Generated Loss, Augmented Training, and Multiscale VAE

abstract: The variational autoencoder (VAE) framework remains a popular option for training unsupervised generative models, especially for discrete data where generative adversarial networks (GANs) require workaround to create gradient for the generator. In our work modeling US postal addresses, we show that our discrete VAE with tree recursive architecture demonstrates limited capability of capturing field correlations within structured data, even after overcoming the challenge of posterior collapse with scheduled sampling and tuning of the KL-divergence weight $\beta$. Worse, VAE seems to have difficulty mapping its generated samples to the latent space, as their VAE loss lags behind or even increases during the training process. Motivated by this observation, we show that augmenting training data with generated variants (augmented training) and training a VAE with multiple values of $\beta$ simultaneously (multiscale VAE) both improve the generation quality of VAE. Despite their differences in motivation and emphasis, we show that augmented training and multiscale VAE are actually connected and have similar effects on the model.


title: Phenomenology of GeV-scale scalar portal

abstract: We review and revise the phenomenology of the scalar portal -- a new scalar particle with the mass in GeV range that mixes with the Higgs boson. In particular, we consider production channels $B\to S K_1(1270)$ and $B\to S K_0^*(700)$ and show that their contribution is significant. We extend the previous analysis by comparing the production of scalars from decays of mesons, of the Higgs bosons and direct production via proton bremsstrahlung, deep inelastic scattering and coherent scattering on nuclei. Relative efficiency of the production channels depends on the energy of the beam and we consider the energies of DUNE, SHiP and LHC-based experiments. We present our results in the form directly suitable for calculations of experimental sensitivities.


title: Supercritical percolation on nonamenable graphs: Isoperimetry, analyticity, and exponential decay of the cluster size distribution

abstract: Let $G$ be a connected, locally finite, transitive graph, and consider Bernoulli bond percolation on $G$. We prove that if $G$ is nonamenable and $p > p_c(G)$ then there exists a positive constant $c_p$ such that \[\mathbf{P}_p(n \leq |K| < \infty) \leq e^{-c_p n}\] for every $n\geq 1$, where $K$ is the cluster of the origin. We deduce the following two corollaries: 1. Every infinite cluster in supercritical percolation on a transitive nonamenable graph has anchored expansion almost surely. This answers positively a question of Benjamini, Lyons, and Schramm (1997). 2. For transitive nonamenable graphs, various observables including the percolation probability, the truncated susceptibility, and the truncated two-point function are analytic functions of $p$ throughout the supercritical phase.


title: Trend-Based Networking Driven by Big Data Telemetry for SDN and Traditional Networks

abstract: Organizations face a challenge of accurately analyzing network data and providing automated action based on the observed trend. This trend-based analytics is beneficial to minimize the downtime and improve the performance of the network services, but organizations use different network management tools to understand and visualize the network traffic with limited abilities to dynamically optimize the network. This research focuses on the development of an intelligent system that leverages big data telemetry analysis in Platform for Network Data Analytics (PNDA) to enable comprehensive trend-based networking decisions. The results include a graphical user interface (GUI) done via a web application for effortless management of all subsystems, and the system and application developed in this research demonstrate the true potential for a scalable system capable of effectively benchmarking the network to set the expected behavior for comparison and trend analysis. Moreover, this research provides a proof of concept of how trend analysis results are actioned in both a traditional network and a software-defined network (SDN) to achieve dynamic, automated load balancing.


title: Latent Variable Algorithms for Multimodal Learning and Sensor Fusion

abstract: Multimodal learning has been lacking principled ways of combining information from different modalities and learning a low-dimensional manifold of meaningful representations. We study multimodal learning and sensor fusion from a latent variable perspective. We first present a regularized recurrent attention filter for sensor fusion. This algorithm can dynamically combine information from different types of sensors in a sequential decision making task. Each sensor is bonded with a modular neural network to maximize utility of its own information. A gating modular neural network dynamically generates a set of mixing weights for outputs from sensor networks by balancing utility of all sensors' information. We design a co-learning mechanism to encourage co-adaption and independent learning of each sensor at the same time, and propose a regularization based co-learning method. In the second part, we focus on recovering the manifold of latent representation. We propose a co-learning approach using probabilistic graphical model which imposes a structural prior on the generative model: multimodal variational RNN (MVRNN) model, and derive a variational lower bound for its objective functions. In the third part, we extend the siamese structure to sensor fusion for robust acoustic event detection. We perform experiments to investigate the latent representations that are extracted; works will be done in the following months. Our experiments show that the recurrent attention filter can dynamically combine different sensor inputs according to the information carried in the inputs. We consider MVRNN can identify latent representations that are useful for many downstream tasks such as speech synthesis, activity recognition, and control and planning. Both algorithms are general frameworks which can be applied to other tasks where different types of sensors are jointly used for decision making.


title: Motzkin Intervals and Valid Hook Configurations

abstract: We define a new natural partial order on Motzkin paths that serves as an intermediate step between two previously-studied partial orders. We provide a bijection between valid hook configurations of $312$-avoiding permutations and intervals in these new posets. We also show that valid hook configurations of permutations avoiding $132$ (or equivalently, $231$) are counted by the same numbers that count intervals in the Motzkin-Tamari posets that Fang recently introduced, and we give an asymptotic formula for these numbers. We then proceed to enumerate valid hook configurations of permutations avoiding other collections of patterns. We also provide enumerative conjectures, one of which links valid hook configurations of $312$-avoiding permutations, intervals in the new posets we have defined, and certain closed lattice walks with small steps that are confined to a quarter plane.


title: Landau Level Degeneracy in Twisted Bilayer Graphene: Role of Symmetry Breaking

abstract: The degeneracy of Landau levels flanking charge neutrality in twisted bilayer graphene is known to change from eight-fold to four-fold when the twist angle is reduced to values near the magic angle of $\approx 1.05^\circ$. This degeneracy lifting has been reproduced in experiments by multiple groups, and is known to occur even in devices which do not harbor the correlated insulators and superconductors. We propose $C_3$ symmetry breaking as an explanation of such robust degeneracy lifting, and support our proposal by numerical results on the Landau level spectrum in near-magic-angle twisted bilayer graphene. Motivated by recent experiments, we further consider the effect of $C_2$ symmetry breaking on the Landau levels.


title: Metaheuristics for Min-Power Bounded-Hops Symmetric Connectivity Problem

abstract: We consider a Min-Power Bounded-Hops Symmetric Connectivity problem that consists of the construction of communication spanning tree on a given graph, where the total energy consumption spent for the data transmission is minimized and the maximum number of edges between two nodes is bounded by some predefined constant. We focus on the planar Euclidian case of this problem where the nodes are placed at the random uniformly spread points on a square and the power cost necessary for the communication between two network elements is proportional to the squared distance between them. Since this is an NP-hard problem, we propose different heuristics based on the following metaheuristics: genetic local search, variable neighborhood search, and ant colony optimization. We perform a posteriori comparative analysis of the proposed algorithms and present the obtained results in this paper.


title: Normalizers and permutational isomorphisms in simply-exponential time

abstract: We show that normalizers and permutational isomorphisms of permutation groups given by generating sets can be computed in time simply exponential in the degree of the groups. The result is obtained by exploiting canonical forms for permutation groups (up to permutational isomorphism).


title: Low Energy Consequences of Loop Quantum Gravity

abstract: As the loop quantum gravity is based on polymer quantization, we will argue that the polymer length (like string length) can be several orders larger than the Planck length, and this can have low energy consequences. We will demonstrate that a short distance modification of a quantum system by polymer quantization and by string theoretical considerations can produce similar behavior. Moreover, it will be demonstrated that a family of different deformed Heisenberg algebras can produce similar low energy effects. We will analyze such polymer correction to a degenerate Fermi gases in a harmonic trap, and its polymer corrected thermodynamics.


title: Real-time observation of electronic, vibrational, and rotational dynamics in nitric oxide with attosecond soft X-ray pulses at 400 eV

abstract: Photoinduced quantum dynamics in molecules have hierarchical temporal structures with different energy scales that are associated with electron and nuclear motions. Femtosecond-to-attosecond transient absorption spectroscopy (TAS) using high-harmonic generation (HHG) with a photon energy below 300 eV has been a powerful tool to observe such electron and nuclear dynamics in a table-top manner. However, comprehensive measurements of the electronic, vibrational, and rotational molecular dynamics have not yet been achieved. Here, we demonstrate HHG-based TAS at the nitrogen K-edge (400 eV) for the first time, and observe all the electronic, vibrational, and rotational degrees of freedom in a nitric oxide molecule at attosecond to sub-picosecond time scales. This method of employing core-to-valence transitions offers an all-optical approach to reveal complete molecular dynamics in photochemical reactions with element and electronic state specificity.


title: Discrete convolution operators and Riesz systems generated by actions of abelian groups

abstract: We study the bounded endomorphisms of $\ell_{N}^2(G)=\ell^2(G)\times \dots \times\ell^2(G)$ that commute with translations, where $G$ is a discrete abelian group. It is shown that they form a C*-algebra isomorphic to the C*-algebra of $N\times N$ matrices with entries in $L^\infty(\widehat{G})$, where $\widehat{G}$ is the dual space of $G$. Characterizations of when these endomorphisms are invertible, and expressions for their norms and for the norms of their inverses, are given. These results allow us to study Riesz systems that arise from the action of $ G $ on a finite set of elements of a Hilbert space.


title: Internal Structure and CO$_2$ Reservoirs of Habitable Water-Worlds

abstract: Water-worlds are water-rich (>1 wt% H$_2$O) exoplanets. The classical models of water-worlds considered layered structures determined by the phase boundaries of pure water. However, water-worlds are likely to possess comet-like compositions, with between ~3 mol% to 30 mol% CO$_2$ relative to water. In this study, we build an interior structure model of habitable (i.e. surface-liquid-ocean-bearing) water-worlds using the latest results from experimental data on the CO$_2$-H$_2$O system, to explore the CO$_2$ budget and to localize the main CO$_2$ reservoirs inside of these planets. We show that CO$_2$ dissolved in the ocean and trapped inside of a clathrate layer cannot accommodate a cometary amount of CO$_2$ if the planet accretes more than 11 wt% of volatiles (CO$_2$ + H$_2$O) during its formation. We propose a new, potentially dominant, CO$_2$ reservoir for water-worlds: CO$_2$ buried inside of the high-pressure water ice mantle as CO$_2$ ices or (H$_2$CO$_3$ . H$_2$O), monohydrate of carbonic acid. If insufficient amounts of CO$_2$ are sequestered either in this reservoir or the planet's iron core, habitable zone water-worlds could generically be stalled in their cooling before liquid oceans have a chance to condense.


title: Ionization bias and the ghost proximity effect near $z\gtrsim6$ quasars in the shadow of proximate absorption systems

abstract: The larger-than-expected scatter in the opacity of the Ly$\alpha$ forest suggests that the metagalactic ionizing background is strongly fluctuating at $z > 5.5$. Models for ionizing background fluctuations predict a strong positive bias on large scales, so the environments of massive $>10^{12}\,{\rm M}_\odot$ dark matter halos, e.g. $z\sim6$ quasar hosts, would be ideal laboratories to constrain the sources of ionizing photons. While the quasars themselves should overwhelm any plausible ionizing photon contribution from neighboring galaxies, proximate damped Ly$\alpha$ absorbers (DLAs) have recently been discovered in the foreground of $z\sim6$ quasars, and the Ly$\alpha$ forest in the shadow of these DLAs could probe the local ionization environment. Using Gpc$^3$ simulations of $z=6$ ionizing background fluctuations, we show that while the Ly$\alpha$ forest signal from ionization bias around a quasar host halo should be strong, it is likely suppressed by the associated intergalactic matter overdensity. We also show that the quasar itself may still overwhelm the clustering signal via a "ghost" of the proximity effect from the quasar radiation causing a large-scale bias in the ionizing photon mean free path. This ghost proximity effect is sensitive to the lifetime and geometry of quasar emission, potentially unlocking a new avenue for constraining these fundamental quasar properties. Finally, we present observations of a $z\sim6$ quasar with a proximate DLA which shows a strong excess in Ly$\alpha$ forest transmission at the predicted location of the ghost proximity effect.


title: Dimensional crossover in a layered ferromagnet detected by spin correlation driven distortions

abstract: Magneto-elastic distortions are commonly detected across magnetic long-range ordering (LRO) transitions. In principle, they are also induced by the magnetic short-range ordering (SRO) that precedes a LRO transition, which contains information about short-range correlations and energetics that are essential for understanding how LRO is established. However these distortions are difficult to resolve because the associated atomic displacements are exceedingly small and do not break symmetry. Here we demonstrate high-multipole nonlinear optical polarimetry as a sensitive and mode selective probe of SRO induced distortions using CrSiTe$_3$ as a testbed. This compound is composed of weakly bonded sheets of nearly isotropic ferromagnetically interacting spins that, in the Heisenberg limit, would individually be impeded from LRO by the Mermin-Wagner theorem. Our results show that CrSiTe$_3$ evades this law via a two-step crossover from two- to three-dimensional magnetic SRO, manifested through two successive and previously undetected totally symmetric distortions above its Curie temperature.


title: The effect of tides on the Sculptor dwarf spheroidal galaxy

abstract: Dwarf spheroidal galaxies (dSphs) appear to be some of the most dark matter dominated objects in the Universe. Their dynamical masses are commonly derived using the kinematics of stars under the assumption of equilibrium. However, these objects are satellites of massive galaxies (e.g.\ the Milky Way) and thus can be influenced by their tidal fields. We investigate the implication of the assumption of equilibrium focusing on the Sculptor dSph by means of ad-hoc $N$-body simulations tuned to reproduce the observed properties of Sculptor following the evolution along some observationally motivated orbits in the Milky Way gravitational field. For this purpose, we used state-of-the-art spectroscopic and photometric samples of Sculptor's stars. We found that the stellar component of the simulated object is not directly influenced by the tidal field, while $\approx 30\%-60\%$ the mass of the more diffuse DM halo is stripped. We conclude that, considering the most recent estimate of the Sculptor proper motion, the system is not affected by the tides and the stellar kinematics represents a robust tracer of the internal dynamics. In the simulations that match the observed properties of Sculptor, the present-day dark-to-luminous mass ratio is $\approx 6$ within the stellar half-light radius ($\approx0.3$ kpc) and $>50$ within the maximum radius of the analysed dataset ($\approx1.5^\circ\approx2$ kpc).


title: The stellar host in star-forming low-mass galaxies: Evidence for two classes

abstract: The morphological evolution of star-forming galaxies provides important clues to understand their physical properties, as well as the triggering and quenching mechanisms of star formation. We aim at connecting morphology and star-formation properties of low-mass galaxies (median stellar mass $\sim$ 10$^{8.5}$ M$_{\odot}$) at low redshift ($z<0.36$). We use a sample of medium-band selected star-forming galaxies from the GOODS-North field. H$\alpha$ images for the sample are created combining both spectral energy distribution fits and HST data. Using them, we mask the star forming regions to obtain an unbiased two-dimensional model of the light distribution of the host galaxies. For this purpose we use $\texttt{PHI}$, a new Bayesian photometric decomposition code. We apply it independently to 7 HST bands assuming a S\'ersic surface brightness model. Star-forming galaxy hosts show low S\'ersic index (with median $n$ $\sim$ 0.9), as well as small sizes (median $R_e$ $\sim$ 1.6 kpc), and negligible change of the parameters with wavelength (except for the axis ratio, which grows with wavelength). Using a clustering algorithm, we find two different classes of star-forming galaxies: A more compact, redder, and high-$n$ (class A) and a more extended, bluer and lower-$n$ one (class B). We also find evidence that the first class is more spheroidal-like. In addition, we find that 48% of the analyzed galaxies present negative color gradients (only 5% are positive). The host component of low-mass star-forming galaxies at $z<0.36$ separates into two different classes, similar to what has been found for their higher mass counterparts. The results are consistent with an evolution from class B to class A. Several mechanisms from the literature, like minor and major mergers, and violent disk instability, can explain the physical process behind the likely transition between the classes. [abridged]


title: Variational Quantum Unsampling on a Quantum Photonic Processor

abstract: Quantum algorithms for Noisy Intermediate-Scale Quantum (NISQ) machines have recently emerged as new promising routes towards demonstrating near-term quantum advantage (or supremacy) over classical systems. In these systems samples are typically drawn from probability distributions which --- under plausible complexity-theoretic conjectures --- cannot be efficiently generated classically. Rather than first define a physical system and then determine computational features of the output state, we ask the converse question: given direct access to the quantum state, what features of the generating system can we efficiently learn? In this work we introduce the Variational Quantum Unsampling (VQU) protocol, a nonlinear quantum neural network approach for verification and inference of near-term quantum circuits outputs. In our approach one can variationally train a quantum operation to unravel the action of an unknown unitary on a known input state; essentially learning the inverse of the black-box quantum dynamics. While the principle of our approach is platform independent, its implementation will depend on the unique architecture of a specific quantum processor. Here, we experimentally demonstrate the VQU protocol on a quantum photonic processor. Alongside quantum verification, our protocol has broad applications; including optimal quantum measurement and tomography, quantum sensing and imaging, and ansatz validation.


title: $\mathtt{bimEX}$: A Mathematica package for exact computations in $3+1$ bimetric relativity

abstract: We present $\mathtt{bimEX}$, a Mathematica package for exact computations in 3$+$1 bimetric relativity. It is based on the $\mathtt{xAct}$ bundle, which can handle computations involving both abstract tensors and their components. In this communication, we refer to the latter case as concrete computations. The package consists of two main parts. The first part involves the abstract tensors, and focuses on how to deal with multiple metrics in $\mathtt{xAct}$. The second part takes an ansatz for the primary variables in a chart as the input, and returns the covariant BSSN bimetric equations in components in that chart. Several functions are implemented to make this process as fast and user-friendly as possible. The package has been used and tested extensively in spherical symmetry and was the workhorse in obtaining the bimetric covariant BSSN equations and reproducing the bimetric 3$+$1 equations in the spherical polar chart.


title: Survey of gravitationally-lensed objects in HSC imaging (SuGOHI). III. Statistical strong lensing constraints on the stellar IMF of CMASS galaxies

abstract: Context: The determination of the stellar initial mass function (IMF) of massive galaxies is one of the open problems in cosmology. Strong gravitational lensing is one of the few methods that allow us to constrain the IMF outside of the Local Group. Aims: The goal of this study is to statistically constrain the distribution in the IMF mismatch parameter, defined as the ratio between the true stellar mass of a galaxy and that inferred assuming a reference IMF, of massive galaxies from the BOSS CMASS sample. Methods: We take 23 strong lenses drawn from the CMASS sample, measure their Einstein radii and stellar masses using multi-band photometry from the Hyper Suprime-Cam survey, then fit a model distribution for the IMF mismatch parameter and dark matter halo mass to the whole sample. We use a prior on halo mass from weak lensing measurements and account for strong lensing selection effects in our model. Results: Assuming an NFW density profile for the dark matter distribution, we infer a value $\mu_{\mathrm{IMF}} = -0.04\pm0.11$ for the average base-10 logarithm of the IMF mismatch parameter, defined with respect to a Chabrier IMF. A Salpeter IMF is in tension with our measurements. Conclusions: Our results are consistent with a scenario in which the region of massive galaxies where the IMF normalization is significantly heavier than that of the Milky Way is much smaller than the scales $5\sim10$~kpc probed by the Einstein radius of the lenses in our sample, as recent spatially resolved studies of the IMF in massive galaxies suggest. The MCMC chains describing the posterior probability distribution of the model are available online, together with the code used to obtain them.


title: JINGLE, a JCMT legacy survey of dust and gas for galaxy evolution studies: II. SCUBA-2 850 {\mu}m data reduction and dust flux density catalogues

abstract: We present the SCUBA-2 850 ${\mu}m$ component of JINGLE, the new JCMT large survey for dust and gas in nearby galaxies, which with 193 galaxies is the largest targeted survey of nearby galaxies at 850 ${\mu}m$. We provide details of our SCUBA-2 data reduction pipeline, optimised for slightly extended sources, and including a calibration model adjusted to match conventions used in other far-infrared data. We measure total integrated fluxes for the entire JINGLE sample in 10 infrared/submillimetre bands, including all WISE, Herschel-PACS, Herschel-SPIRE and SCUBA-2 850 ${\mu}m$ maps, statistically accounting for the contamination by CO(J=3-2) in the 850 ${\mu}m$ band. Of our initial sample of 193 galaxies, 191 are detected at 250 ${\mu}m$ with a $\geq$ 5${\sigma}$ significance. In the SCUBA-2 850 ${\mu}m$ band we detect 126 galaxies with $\geq$ 3${\sigma}$ significance. The distribution of the JINGLE galaxies in far-infrared/sub-millimetre colour-colour plots reveals that the sample is not well fit by single modified-blackbody models that assume a single dust-emissivity index $(\beta)$. Instead, our new 850 ${\mu}m$ data suggest either that a large fraction of our objects require $\beta < 1.5$, or that a model allowing for an excess of sub-mm emission (e.g., a broken dust emissivity law, or a very cold dust component 10 K) is required. We provide relations to convert far-infrared colours to dust temperature and $\beta$ for JINGLE-like galaxies. For JINGLE the FIR colours correlate more strongly with star-formation rate surface-density rather than the stellar surface-density, suggesting heating of dust is greater due to younger rather than older stellar-populations, consistent with the low proportion of early-type galaxies in the sample.


title: Two-dimensional Fermi liquid dynamics with density and quadrupolar interactions

abstract: We consider a Fermi liquid model with density-density as well as quadrupolar forward scattering interactions parametrized by the Landau parameters $F_0$ and $F_2$. Using bosonization and a decimation technique, we compute collective modes and spectral functions for a huge range of interactions, ranging from strong repulsion to strong attraction in either angular momentum channels. We present a dynamical phase diagram showing a region of parameters where the collective modes structure changes abruptly, possibly signaling a dynamical phase transition.


title: Breaking Mirror Twin Hypercharge

abstract: The Twin Higgs scenario stabilizes the Higgs mass through an approximate global symmetry and has remained natural in the face of increasingly stringent LHC bounds on colored top partners. Two basic structural questions in this framework concern the nature of the twin hypercharge gauge symmetry and the origin of the $\mathbb{Z}_2$ symmetry breaking needed to achieve the correct vacuum alignment. Both questions are addressed in a simple extension of the Mirror Twin Higgs model with an exact $\mathbb{Z}_2$ symmetry and a scalar field that spontaneously breaks both twin hypercharge and $\mathbb{Z}_2$. Due to the $\mathbb{Z}_2$ symmetry and an approximate $U(2)$ symmetry in the potential, a new hypercharge scalar appears in the visible sector and, like the Higgs, is a pseudo-Nambu-Goldstone boson with a weak-scale mass. Couplings between the hypercharge scalar and matter provide a new dynamical source of twin sector fermion masses. Depending on the nature and size of these couplings, a variety of experimental signatures may arise, including quark and lepton flavor violation, neutrino masses and mixings as well as direct collider probes of the hypercharged scalar. These signals are correlated with the twin matter spectrum, which can differ dramatically from the visible one, including dynamical realizations of fraternal-like scenarios.


title: Young frames for quantum chemistry

abstract: Quantum chemistry often considers atoms and molecules with non-zero spin. In such cases, the need for proper spin functions results in the theory of configuration state functions. Here, we consider the construction of such wavefunctions using the symmetric group and more specifically Young projectors. We discuss the formalism and detail an example to illustrate the theory. Additionally, we consider the pros and cons of specific implementations of spin symmetry in quantum simulation.