title: Morphisms of Networks of Hybrid Open Systems

abstract: This thesis (defended 10/07/2019) develops a theory of networks of hybrid open systems and morphisms. It builds upon a framework of networks of continuous-time open systems as product and interconnection. We work out categorical notions for hybrid systems, deterministic hybrid systems, hybrid open systems, networks of hybrid open systems, and morphisms of networks of hybrid open systems. We also develop categorical notions for abstract systems, abstract open systems, networks of abstract open systems, and morphisms of networks of abstract open systems. We show that a collection of relations holding among pairs of systems induces a relation between interconnected systems. We use this result for abstract systems to prove a corresponding result for networks of hybrid systems. This result translates as saying that our procedure for building networks preserves morphisms of open systems: a collection of morphisms of (sub)systems is sent to a morphism of networked systems. We thus both justify our formalism and concretize the intuition that a network is a collection of systems pieced together in a certain way.


title: Sharp hypotheses and bispatial inference

abstract: A fundamental class of inferential problems are those characterised by there having been a substantial degree of pre-data (or prior) belief that the value of a model parameter was equal or lay close to a specified value, which may, for example, be the value that indicates the absence of an effect. Standard ways of tackling problems of this type, including the Bayesian method, are often highly inadequate in practice. To address this issue, an inferential framework called bispatial inference is put forward, which can be viewed as both a generalisation and radical reinterpretation of existing approaches to inference that are based on P values. It is shown that to obtain an appropriate post-data density function for a given parameter, it is often convenient to combine a special type of bispatial inference, which is constructed around one-sided P values, with a previously outlined form of fiducial inference. Finally, by using what are called post-data opinion curves, this bispatial-fiducial theory is naturally extended to deal with the general scenario in which any number of parameters may be unknown. The application of the theory is illustrated in various examples, which are especially relevant to the analysis of clinical trial data.


title: A Note on the Swampland Distance Conjecture

abstract: We discuss the Swampland Distance Conjecture in the framework of black hole thermodynamics. In particular, we consider black holes in de Sitter space and we show that the Swampland Distance Conjecture is a consequence of the fact that apparent horizons are always inside cosmic event horizons whenever they exist in the case of fast-roll inflation. In addition, we show that the Bekenstein and the Hubble entropy bounds for the entropy in a region of spacetime lead similarly to the same conjecture.


title: On the Fr\"olicher spectral sequence of the Iwasawa manifold and its small deformations

abstract: We determine the successive pages of the Fr\"olicher spectral sequence of the Iwasawa manifold and some of its small deformations, providing new examples and counterexamples on its properties, including the behaviour under small deformations.


title: Collaborative Machine Learning Markets with Data-Replication-Robust Payments

abstract: We study the problem of collaborative machine learning markets where multiple parties can achieve improved performance on their machine learning tasks by combining their training data. We discuss desired properties for these machine learning markets in terms of fair revenue distribution and potential threats, including data replication. We then instantiate a collaborative market for cases where parties share a common machine learning task and where parties' tasks are different. Our marketplace incentivizes parties to submit high quality training and true validation data. To this end, we introduce a novel payment division function that is robust-to-replication and customized output models that perform well only on requested machine learning tasks. In experiments, we validate the assumptions underlying our theoretical analysis and show that these are approximately satisfied for commonly used machine learning models.


title: Verifiability and Predictability: Interpreting Utilities of Network Architectures for Point Cloud Processing

abstract: In this paper, we diagnose deep neural networks for 3D point cloud processing to explore utilities of different intermediate-layer network architectures. We propose a number of hypotheses on the effects of specific intermediate-layer network architectures on the representation capacity of DNNs. In order to prove the hypotheses, we design five metrics to diagnose various types of DNNs from the following perspectives, information discarding, information concentration, rotation robustness, adversarial robustness, and neighborhood inconsistency. We conduct comparative studies based on such metrics to verify the hypotheses. We further use the verified hypotheses to revise intermediate-layer architectures of existing DNNs and improve their utilities. Experiments demonstrate the effectiveness of our method.


title: Robust Lane Marking Detection Algorithm Using Drivable Area Segmentation and Extended SLT

abstract: In this paper, a robust lane detection algorithm is proposed, where the vertical road profile of the road is estimated using dynamic programming from the v-disparity map and, based on the estimated profile, the road area is segmented. Since the lane markings are on the road area and any feature point above the ground will be a noise source for the lane detection, a mask is created for the road area to remove some of the noise for lane detection. The estimated mask is multiplied by the lane feature map in a bird's eye view (BEV). The lane feature points are extracted by using an extended version of symmetrical local threshold (SLT), which not only considers dark light dark transition (DLD) of the lane markings, like (SLT), but also considers parallelism on the lane marking borders. The segmentation then uses only the feature points that are on the road area. A maximum of two linear lane markings are detected using an efficient 1D Hough transform. Then, the detected linear lane markings are used to create a region of interest (ROI) for parabolic lane detection. Finally, based on the estimated region of interest, parabolic lane models are fitted using robust fitting. Due to the robust lane feature extraction and road area segmentation, the proposed algorithm robustly detects lane markings and achieves lane marking detection with an accuracy of 91% when tested on a sequence from the KITTI dataset.


title: Imaging electric field dynamics with graphene optoelectronics

abstract: The use of electric fields for signalling and control in liquids is widespread, spanning bioelectric activity in cells to electrical manipulation of microstructures in lab-on-a-chip devices. However, an appropriate tool to resolve the spatio-temporal distribution of electric fields over a large dynamic range has yet to be developed. Here we present a label-free method to image local electric fields in real time and under ambient conditions. Our technique combines the unique gate-variable optical transitions of graphene with a critically coupled planar waveguide platform that enables highly sensitive detection of local electric fields with a voltage sensitivity of a few microvolts, a spatial resolution of tens of micrometres and a frequency response over tens of kilohertz. Our imaging platform enables parallel detection of electric fields over a large field of view and can be tailored to broad applications spanning lab-on-a-chip device engineering to analysis of bioelectric phenomena.


title: Fundamental physics with blazar spectra: a critical appraisal

abstract: Very-high-energy (VHE) BL Lac spectra extending above $10 \, \rm TeV$ provide a unique opportunity for testing physics beyond the standard model of elementary particle and alternative blazar emission models. We consider the hadron beam, the photon to axion-like particle (ALP) conversion, and the Lorentz invariance violation (LIV) scenarios by analyzing their consequences and induced modifications to BL Lac spectra. In particular, we consider how different processes can provide similar spectral features (e.g. hard tails) and we discuss the ways they can be disentangled. We use HEGRA data of a high state of Markarian 501 and the HESS spectrum of the extreme BL Lac (EHBL) 1ES 0229+200. In addition, we consider two hypothetical EHBLs similar to 1ES 0229+200 located at redshifts $z=0.3$ and $z=0.5$. We observe that both the hadron beam and the photon-ALP oscillations predict a hard tail extending to energies larger than those possible in the standard scenario. Photon-ALP interaction predicts a peak in the spectra of distant BL Lacs at about $20-30 \, \rm TeV$, while LIV produces a strong peak in all BL Lac spectra around $\sim 100 \, \rm TeV$. The peculiar feature of the photon-ALP conversion model is the production of oscillations in the spectral energy distribution, so that its detection/absence can be exploited to distinguish among the considered models. The above mentioned features coming from the three models may be detected by the upcoming Cherenkov Telescope Array (CTA). Thus, future observations of BL Lac spectra could eventually shed light about new physics and alternative blazar emission models, driving fundamental research towards a specific direction.


title: The perceived effects of group developmental psychology training on agile software development teams

abstract: Research has shown that the maturity of small workgroups from a psychological perspective is intimately connected to team agility. We, therefore, tested if agile team members appreciated group development psychology training. Our results show that the participating teams seem to have a very positive view of group development training and state that they now have a new way of thinking about teamwork and new tools to deal with team-related problems. We, therefore, see huge potential in training agile teams in group development psychology since the positive effects might span over the entire software development organization.


title: Fine-grained Synthesis of Unrestricted Adversarial Examples

abstract: We propose a novel approach for generating unrestricted adversarial examples by manipulating fine-grained aspects of image generation. Unlike existing unrestricted attacks that typically hand-craft geometric transformations, we learn stylistic and stochastic modifications leveraging state-of-the-art generative models. This allows us to manipulate an image in a controlled, fine-grained manner without being bounded by a norm threshold. Our approach can be used for targeted and non-targeted unrestricted attacks on classification, semantic segmentation and object detection models. Our attacks can bypass certified defenses, yet our adversarial images look indistinguishable from natural images as verified by human evaluation. Moreover, we demonstrate that adversarial training with our examples improves performance of the model on clean images without requiring any modifications to the architecture. We perform experiments on LSUN, CelebA-HQ and COCO-Stuff as high resolution datasets to validate efficacy of our proposed approach.


title: There is no maximal decidable expansion of the $\langle \mathbb{N} ,\{ < \} \rangle$ structure

abstract: We are going to prove that if the theory of a structure $\mathcal M=\langle \mathbb{N}, \Sigma \rangle$ is decidable and the standard order $<$ on natural numbers $\mathbb{N}$ is definable in $\mathcal M$, then there is a nontrivial decidable expansion of $\mathcal M$


title: Thermodynamic Equilibrium in General Relativity

abstract: The thermodynamic equilibrium condition for a static self-gravitating fluid in the Einstein theory is defined by the Tolman-Ehrenfest temperature law, $T{\sqrt {g_{00}(x^{i})}} = constant$, according to which the proper temperature depends explicitly on the position within the medium through the metric coefficient $g_{00}(x^{i})$. By assuming the validity of Tolman-Ehrenfest "pocket temperature", Klein also proved a similar relation for the chemical potential, namely, $\mu {\sqrt {g_{00}(x^{i})}} = constant$. In this letter we prove that a more general relation uniting both quantities holds regardless of the equation of state satisfied by the medium, and that the original Tolman-Ehrenfest law form is valid only if the chemical potential vanishes identically. In the general case of equilibrium, the temperature and the chemical potential are intertwined in such a way that only a definite (position dependent) relation uniting both quantities is obeyed. As an illustration of these results, the temperature expressions for an isothermal gas (finite spherical distribution) and a neutron star are also determined.

categories: gr-qc astro-ph.GA cond-mat.stat-mech physics.class-ph


title: Challenges with Extreme Class-Imbalance and Temporal Coherence: A Study on Solar Flare Data

abstract: In analyses of rare-events, regardless of the domain of application, class-imbalance issue is intrinsic. Although the challenges are known to data experts, their explicit impact on the analytic and the decisions made based on the findings are often overlooked. This is in particular prevalent in interdisciplinary research where the theoretical aspects are sometimes overshadowed by the challenges of the application. To show-case these undesirable impacts, we conduct a series of experiments on a recently created benchmark data, named Space Weather ANalytics for Solar Flares (SWAN-SF). This is a multivariate time series dataset of magnetic parameters of active regions. As a remedy for the imbalance issue, we study the impact of data manipulation (undersampling and oversampling) and model manipulation (using class weights). Furthermore, we bring to focus the auto-correlation of time series that is inherited from the use of sliding window for monitoring flares' history. Temporal coherence, as we call this phenomenon, invalidates the randomness assumption, thus impacting all sampling practices including different cross-validation techniques. We illustrate how failing to notice this concept could give an artificial boost in the forecast performance and result in misleading findings. Throughout this study we utilized Support Vector Machine as a classifier, and True Skill Statistics as a verification metric for comparison of experiments. We conclude our work by specifying the correct practice in each case, and we hope that this study could benefit researchers in other domains where time series of rare events are of interest.


title: Quantum Two-Mode Squeezing Radar and Noise Radar: Correlation Coefficients for Target Detection

abstract: Quantum two-mode squeezing (QTMS) radars and noise radars detect targets by correlating the received signal with an internally stored recording. A covariance matrix can be calculated between the two which, in theory, is a function of a single correlation coefficient. This coefficient can be used to decide whether a target is present or absent. We can estimate the correlation coefficient by minimizing the Frobenius norm between the sample covariance matrix and the theoretically expected form of the matrix. Using simulated data, we show that the estimates follow a Rice distribution whose parameters are simple functions of the underlying, "true" correlation coefficient as well as the number of integrated samples. We obtain an explicit expression for the receiver operating characteristic curve that results when the correlation coefficient is used for target detection. This is an important first step toward performance prediction for QTMS radars.


title: Sparse random tensors: Concentration, regularization and applications

abstract: We prove a non-asymptotic concentration inequality for the spectral norm of sparse inhomogeneous random tensors with Bernoulli entries. For an order-$k$ inhomogeneous random tensor $T$ with sparsity $p_{\max}\geq \frac{c\log n}{n }$, we show that $\|T-\mathbb E T\|=O(\sqrt{n p_{\max}}\log^{k-2}(n))$ with high probability. The optimality of this bound up to polylog factors is provided by an information theoretic lower bound. By tensor unfolding, we extend the range of sparsity to $p_{\max}\geq \frac{c\log n}{n^{m}}$ with $1\leq m\leq k-1$ and obtain concentration inequalities for different sparsity regimes. We also provide a simple way to regularize $T$ such that $O(\sqrt{n^{m}p_{\max}})$ concentration still holds down to sparsity $p_{\max}\geq \frac{c}{n^{m}}$ with $k/2\leq m\leq k-1$. We present our concentration and regularization results with two applications: (i) a randomized construction of hypergraphs of bounded degrees with good expander mixing properties, (ii) concentration of sparsified tensors under uniform sampling.


title: Agile Ways of Working: A Team Maturity Perspective

abstract: With the agile approach to managing software development projects comes an increased dependability on well functioning teams, since many of the practices are built on teamwork. The objective of this study was to investigate if, and how, team development from a group psychological perspective is related to some work practices of agile teams. Data were collected from 34 agile teams (200 individuals) from six software development organizations and one university in both Brazil and Sweden using the Group Development Questionnaire (Scale IV) and the Perceptive Agile Measurement (PAM). The result indicates a strong correlation between levels of group maturity and the two agile practices \emph{iterative development} and \emph{retrospectives}. We, therefore, conclude that agile teams at different group development stages adopt parts of team agility differently, thus confirming previous studies but with more data and by investigating concrete and applied agile practices. We thereby add evidence to the hypothesis that an agile implementation and management of agile projects need to be adapted to the group maturity levels of the agile teams.


title: The complexity of promise SAT on non-Boolean domains

abstract: While 3-SAT is NP-hard, 2-SAT is solvable in polynomial time. Austrin, Guruswami, and H\r{a}stad roved a result known as "$(2+\varepsilon)$-SAT is NP-hard" [FOCS'14/SICOMP'17]. They showed that the problem of distinguishing k-CNF formulas that are g-satisfiable (i.e. some assignment satisfies at least g literals in every clause) from those that are not even 1-satisfiable is NP-hard if $\frac{g}{k} < \frac{1}{2}$ and is in P otherwise. We study a generalisation of SAT on arbitrary finite domains, with clauses that are disjunctions of unary constraints, and establish analogous behaviour. Thus we give a dichotomy for a natural fragment of promise constraint satisfaction problems (PCSPs) on arbitrary finite domains. The hardness side is proved using the algebraic approach, via a new general NP-hardness criterion on polymorphisms of the problem, based on a gap version of the Layered Label Cover problem. We show that previously used criteria are insufficient -- the problem hence gives an interesting benchmark of algebraic techniques for proving hardness of approximation problems such as PCSPs.


title: Generalized Ashtekar variables for Palatini f(R) models

abstract: We consider special classes of Palatini f(R) theories, featured by additional Loop Quantum Gravity inspired terms, with the aim of identifying a set of modified Ashtekar canonical variables, which still preserve the SU(2) gauge structure of the standard theory. In particular, we allow for affine connection to be endowed with torsion, which turns out to depend on the additional scalar degree affecting Palatini f(R) gravity, and in this respect we successfully construct a novel Gauss constraint. We analyze the role of the additional scalar field, outlining as it acquires a dynamical character by virtue of a non vanishing Immirzi parameter, and we describe some possible effects on the area operator stemming from such a revised theoretical framework. Finally, we compare our results with earlier studies in literature, discussing differences between metric and Palatini approaches. It is worth noting how the Hamiltonian turns out to be different in the two cases. The results can be reconciled when the analysis is performed in the Einstein frame.


title: Replication-based emulation of the response distribution of stochastic simulators using generalized lambda distributions

abstract: Due to limited computational power, performing uncertainty quantification analyses with complex computational models can be a challenging task. This is exacerbated in the context of stochastic simulators, the response of which to a given set of input parameters, rather than being a deterministic value, is a random variable with unknown probability density function (PDF). Of interest in this paper is the construction of a surrogate that can accurately predict this response PDF for any input parameters. We suggest using a flexible distribution family -- the generalized lambda distribution -- to approximate the response PDF. The associated distribution parameters are cast as functions of input parameters and represented by sparse polynomial chaos expansions. To build such a surrogate model, we propose an approach based on a local inference of the response PDF at each point of the experimental design based on replicated model evaluations. Two versions of this framework are proposed and compared on analytical examples and case studies.


title: Modelagem de Sistemas Audiom\'etricos Usando T\'ecnicas de Computa\c{c}\~ao Flex\'ivel

abstract: In systems identification, the studied phenomena are accompanied by uncertainties, whether arising from measurement data or computational calculations. Interval data provides a valuable way to represent available information on complex problems where uncertainty, inaccuracy, or variability must be taken into account. The present work aims to determine interval parameters for a model considering data measurement uncertainties using the MQ estimator and neural networks. The main objective of this work is to apply this technique in audiometric systems, particularly in the automatic detection of auditory responses taking into consideration the flexible computing concepts that offer solutions tolerant to subjectivity, inaccuracy or uncertainty. It was possible to obtain a model with interval parameters that allow an infinite set of parameters to be evaluated as a limited range. The model was validated using two methods, one with prediction of 2 steps ahead, using the MQ estimator and extended to another representation as neural networks 5 steps delay. It has been shown that interval parameters generate interval results that contain most validation data improving system reliability.


title: A New Characterization of Path Graphs

abstract: Path graphs are intersection graphs of paths in a tree.~In this paper we give a "6\ good characterization" of path graphs, namely, we prove that path graph membership is in $NP\cap CoNP$ without resorting to existing polynomial time algorithms. The characterization is given in terms of the collection of the \emph{attachedness graphs} of a graph, a novel device to deal with the connected components of a graph after the removal of clique separators. On the one hand, the characterization refines and simplifies the characterization of path graphs due to Monma and Wei [C.L.~Monma,~and~V.K.~Wei, Intersection {G}raphs of {P}aths in a {T}ree, J. Combin. Theory Ser. B, 41:2 (1986) 141--181], which we build on, by reducing a constrained vertex coloring problem defined on the \emph{attachedness graphs} to a vertex 2-coloring problem on the same graphs. On the other hand, the characterization allows us to exhibit two exhaustive lists of obstructions to path graph membership in the form of minimal forbidden induced/partial 2-edge colored subgraphs in each of the \emph{attachedness graphs}.


title: EfficientDet: Scalable and Efficient Object Detection

abstract: Model efficiency has become increasingly important in computer vision. In this paper, we systematically study neural network architecture design choices for object detection and propose several key optimizations to improve efficiency. First, we propose a weighted bi-directional feature pyramid network (BiFPN), which allows easy and fast multiscale feature fusion; Second, we propose a compound scaling method that uniformly scales the resolution, depth, and width for all backbone, feature network, and box/class prediction networks at the same time. Based on these optimizations and better backbones, we have developed a new family of object detectors, called EfficientDet, which consistently achieve much better efficiency than prior art across a wide spectrum of resource constraints. In particular, with single model and single-scale, our EfficientDet-D7 achieves state-of-the-art 55.1 AP on COCO test-dev with 77M parameters and 410B FLOPs, being 4x - 9x smaller and using 13x - 42x fewer FLOPs than previous detectors. Code is available at https://github.com/google/automl/tree/master/efficientdet.


title: The Origins and Prevalence of Texture Bias in Convolutional Neural Networks

abstract: Recent work has indicated that, unlike humans, ImageNet-trained CNNs tend to classify images by texture rather than by shape. How pervasive is this bias, and where does it come from? We find that, when trained on datasets of images with conflicting shape and texture, CNNs learn to classify by shape at least as easily as by texture. What factors, then, produce the texture bias in CNNs trained on ImageNet? Different unsupervised training objectives and different architectures have small but significant and largely independent effects on the level of texture bias. However, all objectives and architectures still lead to models that make texture-based classification decisions a majority of the time, even if shape information is decodable from their hidden representations. The effect of data augmentation is much larger. By taking less aggressive random crops at training time and applying simple, naturalistic augmentation (color distortion, noise, and blur), we train models that classify ambiguous images by shape a majority of the time, and outperform baselines on out-of-distribution test sets. Our results indicate that apparent differences in the way humans and ImageNet-trained CNNs process images may arise not primarily from differences in their internal workings, but from differences in the data that they see.


title: Gluing semigroups -- when and how

abstract: Given two semigroups $\langle A\rangle$ and $\langle B\rangle$ in ${\mathbb N}^n$, we wonder when they can be glued, i.e., when there exists a semigroup $\langle C\rangle$ in ${\mathbb N}^n$ such that the defining ideals of the corresponding semigroup rings satisfy that $I_C=I_A+I_B+\langle\rho\rangle$ for some binomial $\rho$. If $n\geq 2$ and $k[A]$ and $k[B]$ are Cohen-Macaulay, we prove that in order to glue them, one of the two semigroups must be degenerate. Then we study the two most degenerate cases: when one of the semigroups is generated by one single element (simple split) and the case where it is generated by at least two elements and all the elements of the semigroup lie on a line. In both cases we characterize the semigroups that can be glued and say how to glue them. Further, in these cases, we conclude that the glued $\langle C\rangle$ is Cohen-Macaulay if and only if both $\langle A\rangle$ and $\langle B\rangle$ are also Cohen-Macaulay. As an application, we characterize precisely the Cohen-Macaulay semigroups that can be glued when $n=2$.


title: Hints, neutrino bounds and WDM constraints from SDSS DR14 Lyman-$\alpha$ and Planck full-survey data

abstract: The Ly-$\alpha$ forest 1D flux power spectrum is a powerful probe of several cosmological parameters. Assuming a $\Lambda$CDM cosmology including massive neutrinos, we find that the latest SDSS DR14 BOSS and eBOSS Ly-$\alpha$ forest data is in very good agreement with current weak lensing constraints on $(\Omega_m, \sigma_8)$ and has the same small level of tension with Planck. We did not identify a systematic effect in the data analysis that could explain this small tension, but we show that it can be reduced in extended cosmological models where the spectral index is not the same on the very different times and scales probed by CMB and Ly-$\alpha$ data. A particular case is that of a $\Lambda$CDM model including a running of the spectral index on top of massive neutrinos. With combined Ly-$\alpha$ and Planck data, we find a slight (3$\sigma$) preference for negative running, $\alpha_s= -0.010 \pm 0.004$ (68% CL). Neutrino mass bounds are found to be robust against different assumptions. In the $\Lambda$CDM model with running, we find $\sum m_\nu <0.11$ eV at the 95% confidence level for combined Ly-$\alpha$ and Planck (temperature and polarisation) data, or $\sum m_\nu < 0.09$ eV when adding CMB lensing and BAO data. We further provide strong and nearly model-independent bounds on the mass of thermal warm dark matter. For a conservative configuration consisting of SDSS data restricted to $z<4.5$ combined with XQ-100 \lya data, we find $m_X > 5.3\;\mathrm{keV}$ (95\%CL).


title: Search to Distill: Pearls are Everywhere but not the Eyes

abstract: Standard Knowledge Distillation (KD) approaches distill the knowledge of a cumbersome teacher model into the parameters of a student model with a pre-defined architecture. However, the knowledge of a neural network, which is represented by the network's output distribution conditioned on its input, depends not only on its parameters but also on its architecture. Hence, a more generalized approach for KD is to distill the teacher's knowledge into both the parameters and architecture of the student. To achieve this, we present a new Architecture-aware Knowledge Distillation (AKD) approach that finds student models (pearls for the teacher) that are best for distilling the given teacher model. In particular, we leverage Neural Architecture Search (NAS), equipped with our KD-guided reward, to search for the best student architectures for a given teacher. Experimental results show our proposed AKD consistently outperforms the conventional NAS plus KD approach, and achieves state-of-the-art results on the ImageNet classification task under various latency settings. Furthermore, the best AKD student architecture for the ImageNet classification task also transfers well to other tasks such as million level face recognition and ensemble learning.


title: Real-Time Emotion Recognition via Attention Gated Hierarchical Memory Network

abstract: Real-time emotion recognition (RTER) in conversations is significant for developing emotionally intelligent chatting machines. Without the future context in RTER, it becomes critical to build the memory bank carefully for capturing historical context and summarize the memories appropriately to retrieve relevant information. We propose an Attention Gated Hierarchical Memory Network (AGHMN) to address the problems of prior work: (1) Commonly used convolutional neural networks (CNNs) for utterance feature extraction are less compatible in the memory modules; (2) Unidirectional gated recurrent units (GRUs) only allow each historical utterance to have context before it, preventing information propagation in the opposite direction; (3) The Soft Attention for summarizing loses the positional and ordering information of memories, regardless of how the memory bank is built. Particularly, we propose a Hierarchical Memory Network (HMN) with a bidirectional GRU (BiGRU) as the utterance reader and a BiGRU fusion layer for the interaction between historical utterances. For memory summarizing, we propose an Attention GRU (AGRU) where we utilize the attention weights to update the internal state of GRU. We further promote the AGRU to a bidirectional variant (BiAGRU) to balance the contextual information from recent memories and that from distant memories. We conduct experiments on two emotion conversation datasets with extensive analysis, demonstrating the efficacy of our AGHMN models.


title: On the M\"obius function in all short intervals

abstract: We show that, for the M\"obius function $\mu(n)$, we have $$ \sum_{x < n\leq x+x^{\theta}}\mu(n)=o(x^{\theta}) $$ for any $\theta>0.55$. This improves on a result of Ramachandra from 1976, which is valid for $\theta>7/12$. Ramachandra's result corresponded to Huxley's $7/12$ exponent for the prime number theorem in short intervals. The main new idea leading to the improvement is using Ramar\'e's identity to extract a small prime factor from the $n$-sum. The proof method also allows us to improve on an estimate of Zhan for the exponential sum of the M\"obius function as well as some results on multiplicative functions and almost primes in short intervals.


title: Grammar Compressed Sequences with Rank/Select Support

abstract: Sequence representations supporting not only direct access to their symbols, but also rank/select operations, are a fundamental building block in many compressed data structures. Several recent applications need to represent highly repetitive sequences, and classical statistical compression proves ineffective. We introduce, instead, grammar-based representations for repetitive sequences, which use up to 6% of the space needed by statistically compressed representations, and support direct access and rank/select operations within tens of microseconds. We demonstrate the impact of our structures in text indexing applications.


title: Comparing hyperbolic and extremal lengths for shortest curves

abstract: We give a linear upper bound for the extremal length in terms of the hyperbolic length for curves of a shortest pair of pants decomposition. As a consequence, we obtain results on the asymptotics of the maximum of the extremal length systole in moduli space for large genus. We also obtain upper bounds for the extremal length of shortest pair of pants decompositions and for the renormalized volume of 3-manifolds.


title: Global stability of fluid flows despite transient growth of energy

abstract: Verifying nonlinear stability of a laminar fluid flow against all perturbations is a central challenge in fluid dynamics. Past results rely on monotonic decrease of a perturbation energy or a similar quadratic generalized energy. None show stability for the many flows that seem to be stable despite these energies growing transiently. Here a broadly applicable method to verify global stability of such flows is presented. It uses polynomial optimization computations to construct non-quadratic Lyapunov functions that decrease monotonically. The method is used to verify global stability of 2D plane Couette flow at Reynolds numbers above the energy stability threshold found by Orr in 1907. This is the first global stability result for any flow that surpasses the energy method.


title: Eigenvectors from eigenvalues revisited

abstract: This is a remark on a recent post by P. Denton, S. Parke, T. Tao, X. Zhang, Eigenvectors from eigenvalues, arXiv:1908.03795v1


title: Note on eigenvectors from eigenvalues

abstract: Denton, Parke, Tao and Zhang gave a new method which determines eigenvectors from eigenvalues for Hermitian matrices with distinct eigenvalues. In this short note, we extend the above result to general Hermitian matrices.


title: Asymptotics of the nucleus ground-state and single-hole nature of the bound states of the nuclei

abstract: We consider nuclei composed of nucleons which interact via two-body potentials decreasing exponentially at infinity. Protons and neutrons are not distinguished in order to simplify notations. The basic result is the rigorous mathematical proof that the bound eigenstates of the nuclei belong to the subspace $P \mathscr{H}$ spanned by the states of a single hole in the ground state $\psi_{0}$ of the parent nucleus with an extra nucleon, as in the uncorrelated models. This follows from the exponential decay of $\psi_{0}$ when a nucleon is very far apart from the residual nucleus. We prove that the real difference from the uncorrelated models is that $P \mathscr{H}$ has an infinite dimension and contains generalized single hole states, distinguishable from the usual ones by the fact that one cannot assign a wave function to the hole. The bound eigenstates of the nuclei are just states of this kind. Some physical consequences are discussed, in particular the unexpected fact that the dynamical part of the single-hole Hamiltonian, although nonnull, does not affect the single-hole overlaps with the bound eigenstates. The decay of $\psi_{0}$ provides the asymptotic behaviours of many single-hole quantities, in particular the nuclear density matrix. Thus a by product of this paper is the rigorous proof of the method developed by Van Neck, Waroquier and Heyde to calculate overlaps and separation energies.


title: Machine-learning non-stationary noise out of gravitational wave detectors

abstract: Signal extraction out of background noise is a common challenge in high precision physics experiments, where the measurement output is often a continuous data stream. To improve the signal to noise ratio of the detection, witness sensors are often used to independently measure background noises and subtract them from the main signal. If the noise coupling is linear and stationary, optimal techniques already exist and are routinely implemented in many experiments. However, when the noise coupling is non-stationary, linear techniques often fail or are sub-optimal. Inspired by the properties of the background noise in gravitational wave detectors, this work develops a novel algorithm to efficiently characterize and remove non-stationary noise couplings, provided there exist witnesses of the noise source and of the modulation. In this work, the algorithm is described in its most general formulation, and its efficiency is demonstrated with examples from the data of the Advanced LIGO gravitational wave observatory, where we could obtain an improvement of the detector gravitational wave reach without introducing any bias on the source parameter estimation.

categories: gr-qc astro-ph.IM cs.LG physics.data-an physics.ins-det


title: Breakdown of Liesegang precipitation bands in a simplified fast reaction limit of the Keller-Rubinow model

abstract: We study solutions to the integral equation \[ \omega(x) = \Gamma - x^2 \int_{0}^1 K(\theta) \, H(\omega(x\theta)) \, \mathrm d \theta \] where $\Gamma>0$, $K$ is a weakly degenerate kernel satisfying, among other properties, $K(\theta) \sim k \, (1-\theta)^\sigma$ as $\theta \to 1$ for constants $k>0$ and $\sigma \in (0, \log_2 3 -1)$, $H$ denotes the Heaviside function, and $x \in [0,\infty)$. This equation arises from a reaction-diffusion equation describing Liesegang precipitation band patterns under certain simplifying assumptions. We argue that the integral equation is an analytically tractable paradigm for the clustering of precipitation rings observed in the full model. This problem is nontrivial as the right hand side fails a Lipschitz condition so that classical contraction mapping arguments do not apply. Our results are the following. Solutions to the integral equation, which initially feature a sequence of relatively open intervals on which $\omega$ is positive ("rings") or negative ("gaps") break down beyond a finite interval $[0,x^*]$ in one of two possible ways. Either the sequence of rings accumulates at $x^*$ ("non-degenerate breakdown") or the solution cannot be continued past one of its zeroes at all ("degenerate breakdown"). Moreover, we show that degenerate breakdown is possible within the class of kernels considered. Finally, we prove existence of generalized solutions which extend the integral equation past the point of breakdown.


title: Statistical Mechanics of the Directed 2-distance Minimal Dominating Set problem

abstract: The directed L-distance minimal dominating set (MDS) problem has wide practical applications in the fields of computer science and communication networks. Here, we study this problem from the perspective of purely theoretical interest. We only give results for an Erd$\acute{o}$s R$\acute{e}$nyi (ER) random graph and regular random graph, but this work can be extended to any type of networks. We develop spin glass theory to study the directed 2-distance MDS problem. First, we find that the belief propagation algorithm does not converge when the inverse temperature exceeds a threshold on either an ER random network or regular random network. Second, the entropy density of replica symmetric theory has a transition point at a finite inverse temperature on a regular random graph when the node degree exceeds 4 and on an ER random graph when the node degree exceeds 6.6; there is no entropy transition point (or $\beta=\infty$) in other circumstances. Third, the results of the replica symmetry (RS) theory are in perfect agreement with those of belief propagation (BP) algorithm while the results of the belief propagation decimation (BPD) algorithm are better than those of the greedy heuristic algorithm.


title: Shapelets for earthquake detection

abstract: This paper introduces EQShapelets (EarthQuake Shapelets) a time-series shape-based approach embedded in machine learning to autonomously detect earthquakes. It promises to overcome the challenges in the field of seismology related to automated detection and cataloging of earthquakes. EQShapelets are amplitude and phase-independent, i.e., their detection sensitivity is irrespective of the magnitude of the earthquake and the time of occurrence. They are also robust to noise and other spurious signals. The detection capability of EQShapelets is tested on one week of continuous seismic data provided by the Northern California Seismic Network (NCSN) obtained from a station in central California near the Calaveras Fault. EQShapelets combined with a Random Forest classifier, detected all of the cataloged earthquakes and 281 uncataloged events with lower false detection rate thus offering a better performance than autocorrelation and FAST algorithms. The primary advantage of EQShapelets over competing methods is the interpretability and insight it offers. Shape-based approaches are intuitive, visually meaningful and offers immediate insight into the problem domain that goes beyond their use in accurate detection. EQShapelets, if implemented at a large scale, can significantly reduce catalog completeness magnitudes and can serve as an effective tool for near real-time earthquake monitoring and cataloging.


title: Quillen metric for singular families of Riemann surfaces with cusps and compact perturbation theorem

abstract: We study the behavior of the Quillen metric for the family of Riemann surfaces with cusps when the additional cusps are created by degeneration. More precisely, in our previous paper, we've seen that the renormalization of the Quillen metric associated with a family of Riemann surfaces with cusps extends continuously over the locus of singular curves. The main result of this article shows that, modulo some explicit universal constant, this continuous extension coincides with the Quillen metric of the normalization of singular curves. This result shows that the Quillen metric is compatible with the adjunction of cusps. When this theorem is applied directly to the moduli space of curves, we obtain the compatibility of the Quillen metric with clutching morphisms in the moduli space of pointed stable curves. As one application, we obtain the compatibility between our definition of the analytic torsion and the definition of Takhtajan-Zograf using lengths of closed geodesics. As a consequence of the proof of the main theorem, we get an explicit relation in terms of Bott-Chern forms between the Quillen metric associated with a cusped metric and the Quillen metric associated with a metric on the compactified Riemann surface. This refines relative compact perturbation theorem we obtained before by pinning down the universal constant.


title: The homeomorphism group of the first uncountable ordinal

abstract: We show that the topology of pointwise convergence on scattered spaces is compatible with the group structure of their homeomorphism group. We then establish a few topological properties of the homeomorphism group of the first uncountable ordinal, such as amenability and Roelcke-precompactness.


title: From deformation theory of wheeled props to classification of Kontsevich formality maps

abstract: We study homotopy theory of the wheeled prop controlling Poisson structures on arbitrary formal graded finite-dimensional manifolds and prove, in particular, that Grothendieck-Teichmueller group acts on that wheeled prop faithfully and homotopy non-trivially. Next we apply this homotopy theory to the study of the deformation complex of an arbitrary Maxim Kontsevich formality map and compute the full cohomology group of that deformation complex in terms of the cohomology of a certain graph complex introduced earlier by Maxim Kontsevich in [K1] and studied by Thomas Willwacher in [W1].


title: Absolute pressure and gas species identification with an optically levitated rotor

abstract: The authors describe a novel variety of spinning-rotor vacuum gauge in which the rotor is a ${\sim}4.7{\text -}\mu$m-diameter silica microsphere, optically levitated. A rotating electrostatic field is used to apply torque to the permanent electric dipole moment of the silica microsphere and control its rotational degrees of freedom. When released from a driving field, the microsphere's angular velocity decays exponentially with a damping time inversely proportional to the residual gas pressure, and dependent on gas composition. The gauge is calibrated by measuring the rotor mass with electrostatic co-levitation, and assuming a spherical shape, confirmed separately, and uniform density. The gauge is cross-checked against a capacitance manometer by observing the torsional drag due to a number of different gas species. The techniques presented can be used to perform absolute vacuum measurements localized in space, owing to the small dimensions of the microsphere and the ability to translate the optical trap in three dimensions, as well as measurements in magnetic field environments. In addition, the dynamics of the microsphere, paired with a calibrated vacuum gauge, can be used to measure the effective molecular mass of a gas mixture without the need for ionization and at pressures up to approximately 1 mbar.


title: Tool-Supported Experiments for Continuously Collecting Data of Subjective Video Quality Assessments During Video Playback

abstract: The adequate use of documentation for communication is one challenge in requirements engineering (RE). In recent years, several researchers addressed this challenge by using videos as a communication mechanism. All of them concluded that this way of using videos has the potential to facilitate requirements communication. Nevertheless, software professionals are not directors and thus do not necessarily know what constitutes a good video. This lack of knowledge is one crucial reason why videos are still not an established communication mechanism in RE. When videos shall be established in the RE activities, practices, and techniques, requirements engineers have to acquire the necessary knowledge to produce and use good videos on their own at moderate costs, yet sufficient quality. In our research project ViViReq (see Acknowledgment), we aspire to bridge this knowledge gap about what constitutes a good video. Whether a video is good or not depends on its quality perceived by its viewers. However, video quality is a rather ill-defined concept due to numerous unspecified technical and subjective characteristics. As part of our research plan, we develop a quality model for videos inspired by the idea of Femmer and Vogelsang to define and evaluate the quality of videos as RE artifacts. In addition to evaluating videos, this quality model can be used to identify the relevant characteristics of videos for their specific purpose which can be further used to specify requirements, their criteria for satisfaction, and corresponding measures. Therefore, software professionals may use the quality model as guidance for producing and using videos.


title: Superpixel Soup: Monocular Dense 3D Reconstruction of a Complex Dynamic Scene

abstract: This work addresses the task of dense 3D reconstruction of a complex dynamic scene from images. The prevailing idea to solve this task is composed of a sequence of steps and is dependent on the success of several pipelines in its execution. To overcome such limitations with the existing algorithm, we propose a unified approach to solve this problem. We assume that a dynamic scene can be approximated by numerous piecewise planar surfaces, where each planar surface enjoys its own rigid motion, and the global change in the scene between two frames is as-rigid-as-possible (ARAP). Consequently, our model of a dynamic scene reduces to a soup of planar structures and rigid motion of these local planar structures. Using planar over-segmentation of the scene, we reduce this task to solving a "3D jigsaw puzzle" problem. Hence, the task boils down to correctly assemble each rigid piece to construct a 3D shape that complies with the geometry of the scene under the ARAP assumption. Further, we show that our approach provides an effective solution to the inherent scale-ambiguity in structure-from-motion under perspective projection. We provide extensive experimental results and evaluation on several benchmark datasets. Quantitative comparison with competing approaches shows state-of-the-art performance.


title: An inductive construction of minimal codes

abstract: We provide new families of minimal codes in any characteristic. Also, an inductive construction of minimal codes is presented.


title: Cosmological volume acceleration in dust epoch: using scaling solutions and variable cosmological term $\Lambda(t)$ within an anisotropic cosmological model

abstract: Under the premise that the current observations of the cosmic microwave background radiation set a very stringent limit to the anisotropy of the universe, we consider an anistropic model in the presence of a barotropic perfect fluid and a homogeneous scalar field, which transits to a flat FRW cosmology for late times in a dust epoch, presenting an accelerated volume expansion. Furtheremore, the scalar field is identified with a varying cosmological term via $V(\phi(t))=2\Lambda(t)$. Exact solutions to the EKG system are obtained by proposing an anisotropic extension of the scaling solutions scenario: $\rm\rho\sim \eta^{-n},\ \rho_\phi\sim \eta^{-m}$, with $\rm\eta^3=a_1a_2a_3$ the volume function of the anistropic model ($\rm a_1,\, a_2,\, a_3$ being the scale factors).


title: Quantum correlations of qubit-qutrit systems under general collective dephasing

abstract: Most studies of collective dephasing for bipartite as well as multipartite quantum systems focus on a very specific orientation of magnetic field, that is, z-orientation. However, in practical situations, there are always small fluctuations in stochastic field and it is necessary that more general orientations of fields should be considered. We extend this problem to qubit-qutrit systems and study correlation dynamics for entanglement and local quantum uncertainty for some specific quantum states. We find that certain quantum states exhibit freezing dynamics both for entanglement and local quantum uncertainty. We analyze the asymptotic states and find the conditions for having non-zero entanglement and local quantum uncertainty. Our results are relevant for ion-trap experiments and can be verified with current experimental setups.


title: Electromagnetic and gravitomagnetic structure and radii of nucleons

abstract: Taking into account the PDFs, obtained by different Collaborations, the momentum transfer dependence of GPDs of the nucleons is obtained. The calculated electromagnetic and gravitomagnetic form factors of nucleons are used for the description of different form factors and the nucleons elastic scattering in a wide energy and momentum transfer region with a minimum number of fitting parameters. The electromagnetic and gravitomagnetic radii of the nucleon are calculated using the obtained momentum transfer dependence of GPDs with different forms of PDFs obtained by different Collaborations. The comparison of the calculations, taking into account the PDFs obtained by different Collaborations, of mean square electromagnetic and gravitomagnetic radii of nucleons is made.


title: Remarks on mass dimension one fermions: The underlying aspects, bilinear forms, Spinor Classification and RIM decomposition

abstract: In the present essay we review the underlying physical information behind the first concrete example describing a mass dimension one fermion - namely Elko spinors. We start the program exploring the physical information by evaluating the Elko bilinear forms, both within the proper orthochronous Lorentz subgroup as well as within the VSR theory. As we shall see, such structures do not hold the right observance of the Fierz-Pauli-Kofink quadratic relations. Thus, by the aforementioned reasons, we develop a deformation of the Clifford algebra basis. Such protocol can be accomplished by taking precisely the right Elko dual structure during the construction of the bilinear forms related to these spinors. With the appropriated bilinear forms at hands, we search for a real physical interpretation in order to achieve a deeper understanding of such spinor fields. Aiming an interesting application, we present a relation concerning Elko spinors and the neutrino physics via the Heisenberg non-linear theory by means of a bijective linear map between Elko spinors and the so-called Restricted Inomata-McKinley (RIM) spinors. Thus, we describe some of its properties. Some interesting results concerning the construction of RIM-decomposable spinors emerge from such prescription.


title: AssemblyNet: A large ensemble of CNNs for 3D Whole Brain MRI Segmentation

abstract: Whole brain segmentation using deep learning (DL) is a very challenging task since the number of anatomical labels is very high compared to the number of available training images. To address this problem, previous DL methods proposed to use a single convolution neural network (CNN) or few independent CNNs. In this paper, we present a novel ensemble method based on a large number of CNNs processing different overlapping brain areas. Inspired by parliamentary decision-making systems, we propose a framework called AssemblyNet, made of two "assemblies" of U-Nets. Such a parliamentary system is capable of dealing with complex decisions, unseen problem and reaching a consensus quickly. AssemblyNet introduces sharing of knowledge among neighboring U-Nets, an "amendment" procedure made by the second assembly at higher-resolution to refine the decision taken by the first one, and a final decision obtained by majority voting. During our validation, AssemblyNet showed competitive performance compared to state-of-the-art methods such as U-Net, Joint label fusion and SLANT. Moreover, we investigated the scan-rescan consistency and the robustness to disease effects of our method. These experiences demonstrated the reliability of AssemblyNet. Finally, we showed the interest of using semi-supervised learning to improve the performance of our method.


title: SINet: Extreme Lightweight Portrait Segmentation Networks with Spatial Squeeze Modules and Information Blocking Decoder

abstract: Designing a lightweight and robust portrait segmentation algorithm is an important task for a wide range of face applications. However, the problem has been considered as a subset of the object segmentation problem and less handled in the semantic segmentation field. Obviously, portrait segmentation has its unique requirements. First, because the portrait segmentation is performed in the middle of a whole process of many real-world applications, it requires extremely lightweight models. Second, there has not been any public datasets in this domain that contain a sufficient number of images with unbiased statistics. To solve the first problem, we introduce the new extremely lightweight portrait segmentation model SINet, containing an information blocking decoder and spatial squeeze modules. The information blocking decoder uses confidence estimates to recover local spatial information without spoiling global consistency. The spatial squeeze module uses multiple receptive fields to cope with various sizes of consistency in the image. To tackle the second problem, we propose a simple method to create additional portrait segmentation data which can improve accuracy on the EG1800 dataset. In our qualitative and quantitative analysis on the EG1800 dataset, we show that our method outperforms various existing lightweight segmentation models. Our method reduces the number of parameters from 2.1M to 86.9K (around 95.9% reduction), while maintaining the accuracy under an 1% margin from the state-of-the-art portrait segmentation method. We also show our model is successfully executed on a real mobile device with 100.6 FPS. In addition, we demonstrate that our method can be used for general semantic segmentation on the Cityscapes dataset. The code and dataset are available in https://github.com/HYOJINPARK/ExtPortraitSeg .


title: Gradient Method for Continuous Influence Maximization with Budget-Saving Considerations

abstract: Continuous influence maximization (CIM) generalizes the original influence maximization by incorporating general marketing strategies: a marketing strategy mix is a vector $\boldsymbol x = (x_1,\dots,x_d)$ such that for each node $v$ in a social network, $v$ could be activated as a seed of diffusion with probability $h_v(\boldsymbol x)$, where $h_v$ is a strategy activation function satisfying DR-submodularity. CIM is the task of selecting a strategy mix $\boldsymbol x$ with constraint $\sum_i x_i \le k$ where $k$ is a budget constraint, such that the total number of activated nodes after the diffusion process, called influence spread and denoted as $g(\boldsymbol x)$, is maximized. In this paper, we extend CIM to consider budget saving, that is, each strategy mix $\boldsymbol x$ has a cost $c(\boldsymbol x)$ where $c$ is a convex cost function, we want to maximize the balanced sum $g(\boldsymbol x) + \lambda(k - c(\boldsymbol x))$ where $\lambda$ is a balance parameter, subject to the constraint of $c(\boldsymbol x) \le k$. We denote this problem as CIM-BS. The objective function of CIM-BS is neither monotone, nor DR-submodular or concave, and thus neither the greedy algorithm nor the standard result on gradient method could be directly applied. Our key innovation is the combination of the gradient method with reverse influence sampling to design algorithms that solve CIM-BS: For the general case, we give an algorithm that achieves $\left(\frac{1}{2}-\varepsilon\right)$-approximation, and for the case of independent strategy activations, we present an algorithm that achieves $\left(1-\frac{1}{e}-\varepsilon\right)$ approximation.


title: Safe Policies for Reinforcement Learning via Primal-Dual Methods

abstract: In this paper, we study the learning of safe policies in the setting of reinforcement learning problems. This is, we aim to control a Markov Decision Process (MDP) of which we do not know the transition probabilities, but we have access to sample trajectories through experience. We define safety as the agent remaining in a desired safe set with high probability during the operation time. We therefore consider a constrained MDP where the constraints are probabilistic. Since there is no straightforward way to optimize the policy with respect to the probabilistic constraint in a reinforcement learning framework, we propose an ergodic relaxation of the problem. The advantages of the proposed relaxation are threefold. (i) The safety guarantees are maintained in the case of episodic tasks and they are kept up to a given time horizon for continuing tasks. (ii) The constrained optimization problem despite its non-convexity has arbitrarily small duality gap if the parametrization of the policy is rich enough. (iii) The gradients of the Lagrangian associated with the safe-learning problem can be easily computed using standard policy gradient results and stochastic approximation tools. Leveraging these advantages, we establish that primal-dual algorithms are able to find policies that are safe and optimal. We test the proposed approach in a navigation task in a continuous domain. The numerical results show that our algorithm is capable of dynamically adapting the policy to the environment and the required safety levels.


title: An Accelerating Cosmological Model from a Parametrization of Hubble Parameter

abstract: In view of late-time cosmic acceleration, a dark energy cosmological model is revisited wherein Einstein's cosmological constant is considered as a candidate of dark energy. Exact solution of Einstein field equations (EFEs) is derived in a homogeneous isotropic background in classical general relativity. The solution procedure is adopted, in a model independent way (or the cosmological parametrization). A simple parametrization of the Hubble parameter (H) as a function of cosmic time `t' is considered which produces an exponential type of evolution of the scale factor (a) and also shows a negative value of deceleration parameter at the present time with a signature flip from early deceleration to late acceleration. Cosmological dynamics of the model obtained have been discussed illustratively for different phases of the evolution of the universe. The evolution of different cosmological parameters are shown graphically for at and closed cases of Friedmann-Lemaitre-Robertson-Walker (FLRW) space-time for the presented model (open case is incompatible to the present scenario). We have also constrained our model parameters with the updated (36 points) observational Hubble dataset.


title: Investigating Active Learning and Meta-Learning for Iterative Peptide Design

abstract: Often the development of novel functional peptides is not amenable to high throughput or purely computational screening methods. Peptides must be synthesized one at a time in a process that does not generate large amounts of data. One way this method can be improved is by ensuring that each experiment provides the best improvement in both peptide properties and predictive modeling accuracy. Here, we study the effectiveness of active learning, optimizing experiment order, and meta-learning, transferring knowledge between contexts, to reduce the number of experiments necessary to build a predictive model. We present a multi-task benchmark database of peptides designed to advance these methods for experimental design. Each task is binary classification of peptides represented as a sequence string. We find neither active learning method tested to be better than random choice. The meta-learning method Reptile was found to improve average accuracy across datasets. Combining meta-learning with active learning offers inconsistent benefits.


title: Reversible Computation in Wireless Communications

abstract: This chapter presents the pioneering work in applying reversible computation paradigms to wireless communications. These applications range from developing reversible hardware architectures for underwater acoustic communications to novel distributed optimisation procedures in large radio-frequency antenna arrays based on reversing Petri nets. Throughout the chapter, we discuss the rationale for introducing reversible computation in the domain of wireless communications, exploring the inherently reversible properties of communication channels and systems formed by devices in a wireless network.


title: On Universal Features for High-Dimensional Learning and Inference

abstract: We consider the problem of identifying universal low-dimensional features from high-dimensional data for inference tasks in settings involving learning. For such problems, we introduce natural notions of universality and we show a local equivalence among them. Our analysis is naturally expressed via information geometry, and represents a conceptually and computationally useful analysis. The development reveals the complementary roles of the singular value decomposition, Hirschfeld-Gebelein-R\'enyi maximal correlation, the canonical correlation and principle component analyses of Hotelling and Pearson, Tishby's information bottleneck, Wyner's common information, Ky Fan $k$-norms, and Brieman and Friedman's alternating conditional expectations algorithm. We further illustrate how this framework facilitates understanding and optimizing aspects of learning systems, including multinomial logistic (softmax) regression and the associated neural network architecture, matrix factorization methods for collaborative filtering and other applications, rank-constrained multivariate linear regression, and forms of semi-supervised learning.


title: All-Multiplicity Non-Planar MHV Amplitudes in sYM at Two Loops

abstract: We give a closed-form, prescriptive representation of all-multiplicity two-loop MHV amplitude integrands in fully-color-dressed (non-planar) maximally supersymmetric Yang-Mills theory.


title: OmniFold: A Method to Simultaneously Unfold All Observables

abstract: Collider data must be corrected for detector effects ("unfolded") to be compared with many theoretical calculations and measurements from other experiments. Unfolding is traditionally done for individual, binned observables without including all information relevant for characterizing the detector response. We introduce OmniFold, an unfolding method that iteratively reweights a simulated dataset, using machine learning to capitalize on all available information. Our approach is unbinned, works for arbitrarily high-dimensional data, and naturally incorporates information from the full phase space. We illustrate this technique on a realistic jet substructure example from the Large Hadron Collider and compare it to standard binned unfolding methods. This new paradigm enables the simultaneous measurement of all observables, including those not yet invented at the time of the analysis.


title: Open System Tensor Networks and Kramers' Crossover for Quantum Transport

abstract: Tensor networks are a powerful tool for many-body ground states with limited entanglement. These methods can nonetheless fail for certain time-dependent processes - such as quantum transport or quenches - where entanglement growth is linear in time. Matrix-product-state decompositions of the resulting out-of-equilibrium states require a bond dimension that grows exponentially, imposing a hard limit on simulation timescales. However, in the case of transport, if the reservoir modes of a closed system are arranged according to their scattering structure, the entanglement growth can be made logarithmic. Here, we apply this ansatz to open systems via extended reservoirs that have explicit relaxation. This enables transport calculations that can access steady states, time dynamics and noise, and periodic driving (e.g., Floquet states). We demonstrate the approach by calculating the transport characteristics of an open, interacting system. These results open a path to scalable and numerically systematic many-body transport calculations with tensor networks.


title: Improved global fit to Non-Standard neutrino Interactions using COHERENT energy and timing data

abstract: We perform a global fit to neutrino oscillation and coherent neutrino-nucleus scattering data, using both timing and energy information from the COHERENT experiment. The results are used to set model-independent bounds on four-fermion effective operators inducing non-standard neutral-current neutrino interactions. We quantify the allowed ranges for their Wilson coefficients, as well as the status of the LMA-D solution, for a wide class of new physics models with arbitrary ratios between the strength of the operators involving up and down quarks. Our results are presented for the COHERENT experiment alone, as well as in combination with the global data from oscillation experiments. We also quantify the dependence of our results for COHERENT with respect to the choice of quenching factor, nuclear form factor, and the treatment of the backgrounds.


title: Neutrino propagation hinders fast pairwise flavor conversions

abstract: Neutrino flavor conversions may dramatically affect the inner working of compact astrophysical objects as well as the synthesis of the heavier elements. We present the first sophisticated numerical solution of the neutrino flavor conversion within a (2+1+1) dimensional setup: we include the advective term in the neutrino equations of motion and track the flavor evolution in two spatial dimensions, one angular variable, and time. Notably, the advective term hinders the development of neutrino pairwise conversions, if the conditions favoring such conversions (i.e., crossings between the angular distributions of $\nu_e$ and $\bar\nu_e$ or a non-negligible flux of neutrinos traveling backward with respect to the main propagation direction) exist for time scales shorter than the typical time scale of the advective term. As a consequence, fast pairwise conversions can only occur when the conditions favoring flavor conversions are self-sustained and global, such as the ones induced by the lepton emission self-sustained asymmetry (LESA) in core-collapse supernovae. Our work highlights the major impact of the dynamical evolution of the neutrino field on the growth of flavor instabilities and the strong interplay between classical and quantum effects. Critical limitations of the linear stability analysis, used to predict neutrino flavor instabilities, are also pointed out.


title: Long-time asymptotics of solutions to the Keller-Rubinow model for Liesegang rings in the fast reaction limit

abstract: We consider the Keller--Rubinow model for Liesegang rings in one spatial dimension in the fast reaction limit as introduced by Hilhorst, van der Hout, Mimura, and Ohnishi in 2007. Numerical evidence suggests that solutions to this model converge, independent of the initial concentration, to a universal profile for large times in parabolic similarity coordinates. For the concentration function, the notion of convergence appears to be similar to attraction to a stable equilibrium point in phase space. The reaction term, however, is discontinuous so that it can only convergence in a much weaker, averaged sense. This also means that most of the traditional analytical tools for studying the long-time behavior fail on this problem. In this paper, we identify the candidate limit profile as the solution of a certain one-dimensional boundary value problem which can be solved explicitly. We distinguish two nontrivial regimes. In the first, the transitional regime, precipitation is restricted to a bounded region in space. We prove that the concentration converges to a single asymptotic profile. In the second, the supercritical regime, we show that the concentration converges to one of a one-parameter family of asymptotic profiles, selected by a solvability condition for the one-dimensional boundary value problem. Here, our convergence result is only conditional: we prove that if convergence happens, either pointwise for the concentration or in an averaged sense for the precipitation function, then the other field converges likewise; convergence in concentration is uniform, and the asymptotic profile is indeed the profile selected by the solvability condition. A careful numerical study suggests that the actual behavior of the equation is indeed the one suggested by the theorem.


title: Core-collapse, superluminous, and gamma-ray burst supernova host galaxy populations at low redshift: the importance of dwarf and starbursting galaxies

abstract: We present a comprehensive study of an unbiased sample of 150 nearby (<z> = 0.014) core-collapse supernova (CCSN) host galaxies drawn from the All-Sky Automated Survey for Supernovae (ASAS-SN) for direct comparison to the nearest LGRB and SLSN hosts. We use public imaging surveys to gather multi-wavelength photometry for all CCSN host galaxies and fit their spectral energy distributions (SEDs) to derive stellar masses and integrated star formation rates. CCSNe populate galaxies across a wide range of stellar masses, from blue and compact dwarf galaxies to large spiral galaxies. We find 33(+4,-4) per cent of CCSNe are in dwarf galaxies (M < 10^9 M_Sun) and 2(+2,-1) per cent are in dwarf starburst galaxies (sSFR > 10^-8 yr^-1). We reanalyse low-redshift SLSN and LGRB hosts from the literature (out to $z<0.3$) in a homogeneous way and compare against the CCSN host sample. The relative SLSN to CCSN supernova rate is increased in low-mass galaxies and at high specific star-formation rates. These parameters are strongly covariant and we cannot break the degeneracy between them with our current sample, although there is some evidence that both factors may play a role. Larger unbiased samples of CCSNe from projects such as ZTF and LSST will be needed to determine whether host-galaxy mass (a proxy for metallicity) or specific star-formation rate (a proxy for star-formation intensity and potential IMF variation) is more fundamental in driving the preference for SLSNe and LGRBs in unusual galaxy environments.


title: The Role of Quasar Radiative Feedback on Galaxy Formation during Cosmic Reionization

abstract: Recent observations have found that many $z\sim 6$ quasar fields lack galaxies. This unexpected lack of galaxies may potentially be explained by quasar radiation feedback. In this paper I present a suite of 3D radiative transfer cosmological simulations of quasar fields. I find that quasar radiation suppresses star formation in low mass galaxies, mainly by photo-dissociating their molecular hydrogen. Photo-heating also plays a role, but only after $\sim$100 Myr. However, galaxies which already have stellar mass above $10^5 M_\odot$ when the quasar turns on will not be suppressed significantly. Quasar radiative feedback suppresses the faint end of the galaxy luminosity function (LF) within $1$ pMpc, but to a far lesser degree than the field-to-field variation of the LF. My study also suggests that by using the number of bright galaxies ($M_{1500}<-16$) around quasars, we can potentially recover the underlying mass overdensity, which allows us to put reliable constraints on quasar environments.


title: Multivariate Time-series Analysis of Variable Objects in the Gaia Mission

abstract: In astronomy, we are witnessing an enormous increase in the number of source detections, precision, and diversity of measurements. Additionally, multi-epoch data is becoming the norm, making time-series analyses an important aspect of current astronomy. The Gaia mission is an outstanding example of a multi-epoch survey that provides measurements in a large diversity of domains, with its broad-band photometry; spectrophotometry in blue and red (used to derive astrophysical parameters); spectroscopy (employed to infer radial velocities, v sin(i), and other astrophysical parameters); and its extremely precise astrometry. Most of all that information is provided for sources covering the entire sky. Here, we present several properties related to the Gaia time series, such as the time sampling; the different types of measurements; the Gaia G, G BP and G RP-band photometry; and Gaia-inspired studies using the CORrelation-RAdial-VELocities data to assess the potential of the information on the radial velocity, the FWHM, and the contrast of the cross-correlation function. We also present techniques (which are used or are under development) that optimize the extraction of astrophysical information from the different instruments of Gaia, such as the principal component analysis and the multi-response regression. The detailed understanding of the behavior of the observed phenomena in the various measurement domains can lead to richer and more precise characterization of the Gaia data, including the definition of more informative attributes that serve as input to (our) machine-learning algorithms.


title: Gauge-invariant observables in gravity and electromagnetism: black hole backgrounds and null dressings

abstract: We address questions regarding construction and implications of gauge-invariant "dressed" observables in nontrivial background geometries such as that of a black hole. Formally, such observables can be constructed, e.g. by locating points with geodesics launched from infinity. However, practical complications arise in non-trivial geometries, and in particular for observables behind black hole horizons. Greater simplicity can be achieved by considering null constructions where the dressing lies along a null geodesic, or null surface such as a cone. We first investigate basic properties of these null dressings in the simpler context of electromagnetism. Since null constructions provide simple dressings for gauge-invariant observables inside black holes, they also allow us to investigate the question of compatibility of observables inside and outside black holes, and in particular the idea of black hole complementarity. While such observables in general have non-vanishing state-dependent commutators, the failure to commute does not appear particularly enhanced by the presence of the horizon.


title: Baryonic clues to the puzzling diversity of dwarf galaxy rotation curves

abstract: We use a compilation of disc galaxy rotation curves to assess the role of the luminous component ("baryons") in the rotation curve diversity problem. As in earlier work, we find that rotation curve shape correlates with baryonic surface density: high surface density galaxies have rapidly-rising rotation curves consistent with cuspy cold dark matter halos; slowly-rising rotation curves (characteristic of galaxies with inner mass deficits or "cores") occur only in low surface density galaxies. The correlation, however, seems too weak to be the main driver of the diversity. In addition, dwarf galaxies exhibit a clear trend, from "cuspy" systems where baryons are unimportant in the inner mass budget to "cored" galaxies where baryons actually dominate. This trend constrains the various scenarios proposed to explain the diversity, such as (i) baryonic inflows and outflows during galaxy formation; (ii) dark matter self-interactions; (iii) variations in the baryonic mass structure coupled to rotation velocities through the "mass discrepancy-acceleration relation" (MDAR); or (iv) non-circular motions in gaseous discs. Together with analytical modeling and cosmological hydrodynamical simulations, our analysis shows that each of these scenarios has promising features, but none seems to fully account for the observed diversity. The MDAR, in particular, is inconsistent with the observed trend between rotation curve shape and baryonic importance; either the trend is caused by systematic errors in the data or the MDAR does not apply. The origin of the dwarf galaxy rotation curve diversity and its relation to the structure of cold dark matter halos remains an open issue.


title: Automatic Differentiable Monte Carlo: Theory and Application

abstract: Differentiable programming has emerged as a key programming paradigm empowering rapid developments of deep learning while its applications to important computational methods such as Monte Carlo remain largely unexplored. Here we present the general theory enabling infinite-order automatic differentiation on expectations computed by Monte Carlo with unnormalized probability distributions, which we call "automatic differentiable Monte Carlo" (ADMC). By implementing ADMC algorithms on computational graphs, one can also leverage state-of-the-art machine learning frameworks and techniques to traditional Monte Carlo applications in statistics and physics. We illustrate the versatility of ADMC by showing some applications: fast search of phase transitions and accurately finding ground states of interacting many-body models in two dimensions. ADMC paves a promising way to innovate Monte Carlo in various aspects to achieve higher accuracy and efficiency, e.g. easing or solving the sign problem of quantum many-body models through ADMC.

categories: physics.comp-ph cond-mat.dis-nn cond-mat.stat-mech cond-mat.str-el cs.LG


title: Evidence of Rotational Fr\"ohlich Coupling in Polaronic Trions

abstract: Electrons commonly couple through Fr\"ohlich interactions with longitudinal optical phonons to form polarons. However, trions possess a finite angular momentum and should therefore couple instead to rotational optical phonons. This creates a polaronic trion whose binding energy is determined by the crystallographic orientation of the lattice. Here, we demonstrate theoretically within the Fr\"ohlich approach and experimentally by photoluminescence emission that the bare trion binding energy (20 meV) is significantly enhanced by the phonons at the interface between the two-dimensional semiconductor MoS$_2$ and the bulk transition metal oxide SrTiO$_3$. The low-temperature {binding energy} changes from 60 meV in [001]-oriented substrates to 90 meV for [111] orientation, as a result of the counter-intuitive interplay between the rotational axis of the MoS$_2$ trion and that of the SrTiO$_3$ phonon mode.


title: Aspects of the Normal State Resistivity of Cuprate Superconductors

abstract: Planar normal state resistivity data taken from three families of cuprate superconductors are compared with theoretical calculations from the recent extremely correlated Fermi liquid theory (ECFL). The two hole doped cuprate materials $LSCO$ and $BSLCO$ and the electron doped material $LCCO$ have yielded rich data sets at several densities $\delta$ and temperatures T, thereby enabling a systematic comparison with theory. The recent ECFL resistivity calculations for the highly correlated $t$-$t'$-$J$ model by us give the resistivity for a wide set of model parameters. After using X-ray diffraction and angle resolved photoemission data to fix parameters appearing in the theoretical resistivity, only one parameter, the magnitude of the hopping $t$, remains undetermined. For each data set, the slope of the experimental resistivity at a single temperature-density point is sufficient to determine $t$, and hence the resistivity on absolute scale at all remaining densities and temperatures. This procedure is shown to give a fair account of the entire data.


title: Searching for Sterile Neutrino with X-ray Intensity Mapping

abstract: The cosmological X-ray emission associated to the possible radiative decay of sterile neutrinos is composed by a collection of lines at different energies. For a given mass, each line corresponds to a given redshift. In this work, we cross correlate such line emission with catalogs of galaxies tracing the dark matter distribution at different redshifts. We derive observational prospects by correlating the X-ray sky that will be probed by the eROSITA and Athena missions with current and near future photometric and spectroscopic galaxy surveys. A relevant and unexplored fraction of the parameter space of sterile neutrinos can be probed by this technique.


title: Complementarity of Peculiar Velocity Surveys and Redshift Space Distortions for Testing Gravity

abstract: Peculiar-velocity surveys of the low-redshift universe have significant leverage to constrain the growth rate of cosmic structure and test gravity. Wide-field imaging surveys combined with multi-object spectrographs (e.g. ZTF2, LSST, DESI, 4MOST) can use Type Ia supernovae as informative tracers of the velocity field, reaching few percent constraints on the growth rate $f\sigma_8$ at $z\lesssim0.2$ where density tracers cannot do better than $\sim10\%$. Combining the high-redshift DESI survey mapping redshift space distortions with a low-redshift supernova peculiar velocity survey using LSST and DESI can determine the gravitational growth index to $\sigma(\gamma)\approx0.02$, testing general relativity. We study the characteristics needed for the peculiar velocity survey, and how its complementarity with clustering surveys improves when going from a $\Lambda$CDM model assumption to a $w_0$-$w_a$ cosmology.


title: Nonlinear Bell inequality for macroscopic measurements

abstract: The correspondence principle suggests that quantum systems grow classical when large. Classical systems cannot violate Bell inequalities. Yet agents given substantial control can violate Bell inequalities proven for large-scale systems. We consider agents who have little control, implementing only general operations suited to macroscopic experimentalists: preparing small-scale entanglement and measuring macroscopic properties while suffering from noise. That experimentalists so restricted can violate a Bell inequality appears unlikely, in light of earlier literature. Yet we prove a Bell inequality that such an agent can violate, even if experimental errors have variances that scale as the system size. A violation implies nonclassicality, given limitations on particles' interactions. A product of singlets violates the inequality; experimental tests are feasible for photons, solid-state systems, atoms, and trapped ions. Consistently with known results, violations of our Bell inequality cannot disprove local hidden-variables theories. By rejecting the disproof goal, we show, one can certify nonclassical correlations under reasonable experimental assumptions.


title: A QCD Debye mass in a broad temperature range

abstract: The Debye mass sets a scale for the screening of static charges and the scattering of fast charges within a gauge plasma. Inspired by its potential cosmological applications, we determine a QCD Debye mass at 2-loop order in a broad temperature range (1 GeV ... 10 TeV), demonstrating how quark mass thresholds get smoothly crossed. Along the way, integration-by-parts identities pertinent to massive loops at finite temperature are illuminated.


title: To beta or not to beta: can higher-order Jeans analysis break the mass-anisotropy degeneracy in simulated dwarfs?

abstract: We test a non-parametric higher-order Jeans analysis method, GravSphere, on 32 simulated dwarf galaxies comparable to classical Local Group dwarfs like Fornax. The galaxies are selected from the APOSTLE suite of cosmological hydrodynamics simulations with Cold Dark Matter (CDM) and Self-Interacting Dark Matter (SIDM) models, allowing us to investigate cusps and cores in density distributions. We find that, for CDM dwarfs, the recovered enclosed mass profiles have a bias of no more than 10 per cent, with a 50 per cent scatter in the inner regions and a 20 per cent scatter near the half-light radius, consistent with standard mass estimators. The density profiles are also recovered with a bias of no more than 10 per cent and a scatter of 30 per cent in the inner regions. For SIDM dwarfs, the mass and density profiles are recovered within our 95 per cent confidence intervals, but are biased towards cuspy dark matter distributions. This is mainly due to a lack of sufficient constraints from the data. We explore the sources of scatter in the accuracy of the recovered profiles and suggest a $\chi^2$ statistic to separate successful models from biased ones. Finally, we show that the uncertainties on the mass profiles obtained with GravSphere are smaller than those for comparable Jeans methods, and that they can be further improved if stronger priors, motivated by cosmological simulations, are placed on the velocity anisotropy. We conclude that GravSphere is a promising Jeans-based approach for modelling dark matter distributions in dwarf galaxies.


title: A Dynamical Survey of Stellar-Mass Black Holes in 50 Milky Way Globular Clusters

abstract: Recent numerical simulations of globular clusters (GCs) have shown that stellar-mass black holes (BHs) play a fundamental role in driving cluster evolution and shaping their present-day structure. Rapidly mass-segregating to the center of GCs, BHs act as a dynamical energy source via repeated super-elastic scattering, delaying onset of core collapse and limiting mass segregation for visible stars. While recent discoveries of BH candidates in Galactic and extragalactic GCs have further piqued interest in BH-mediated cluster dynamics, numerical models show that even if significant BH populations remain in today's GCs, they are typically in configurations that are not directly detectable. We demonstrated in Weatherford et al. (2018) that an anti-correlation between a suitable measure of mass segregation ($\Delta$) in observable stellar populations and the number of retained BHs in GC models can be applied to indirectly probe BH populations in real GCs. Here, we estimate the number and total mass of BHs presently retained in 50 Milky Way GCs from the ACS Globular Cluster Survey by measuring $\Delta$ between populations of main sequence stars, using correlations found between $\Delta$ and BH retention in the CMC Cluster Catalog models. We demonstrate that the range in $\Delta$'s distribution from our models matches that for observed GCs to a remarkable degree. Our results further provide the narrowest constraints to-date on the retained BH populations in the GCs analyzed. Of these 50 GCs, we identify NGCs 2808, 5927, 5986, 6101, and 6205 to presently contain especially large BH populations, each with total BH mass exceeding $10^3\,\rm{M_{\odot}}$.


title: Incompressibility of classical distributions

abstract: In blind compression of quantum states, a sender Alice is given a specimen of a quantum state $\rho$ drawn from a known ensemble (but without knowing what $\rho$ is), and she transmits sufficient quantum data to a receiver Bob so that he can decode a near perfect specimen of $\rho$. For many such states drawn iid from the ensemble, the asymptotically achievable rate is the number of qubits required to be transmitted per state. The Holevo information is a lower bound for the achievable rate, and is attained for pure state ensembles, or in the related scenario of entanglement-assisted visible compression of mixed states wherein Alice knows what state is drawn. In this paper, we prove a general and robust lower bound on the achievable rate for ensembles of classical states, which holds even in the least demanding setting when Alice and Bob share free entanglement and a constant per-copy error is allowed. We apply the bound to a specific ensemble of only two states and prove a near-maximal separation (saturating the dimension bound in leading order) between the best achievable rate and the Holevo information for constant error. This also implies that the ensemble is incompressible -- compression does not reduce the communication cost by much. Since the states are classical, the observed incompressibility is not fundamentally quantum mechanical. We lower bound the difference between the achievable rate and the Holevo information in terms of quantitative limitations to clone the specimen or to distinguish the two classical states.


title: A low-scale flavon model with a $Z_N$ symmetry

abstract: We propose a model that explains the fermion mass hierarchy by the Froggatt-Nielsen mechanism with a discrete $Z_N^F$ symmetry. As a concrete model, we study a supersymmetric model with a single flavon coupled to the minimal supersymmetric Standard Model. Flavon develops a TeV scale vacuum expectation value for realizing flavor hierarchy, an appropriate $\mu$-term and the electroweak scale, hence the model has a low cutoff scale. We demonstrate how the flavon is successfully stabilized together with the Higgs bosons in the model. The discrete flavor symmetry $Z_N^F$ controls not only the Standard Model fermion masses, but also the Higgs potential and a mass of the Higgsino which is a good candidate for dark matter. The hierarchy in the Higgs-flavon sector is determined in order to make the model anomaly-free and realize a stable electroweak vacuum. We show that this model can explain the fermion mass hierarchy, realistic Higgs-flavon potential and thermally produced dark matter at the same time. We discuss flavor violating processes induced by the light flavon which would be detected in future experiments.


title: A Scrambled Method of Moments

abstract: Quasi-Monte Carlo (qMC) methods are a powerful alternative to classical Monte-Carlo (MC) integration. Under certain conditions, they can approximate the desired integral at a faster rate than the usual Central Limit Theorem, resulting in more accurate estimates. This paper explores these methods in a simulation-based estimation setting with an emphasis on the scramble of Owen (1995). For cross-sections and short-panels, the resulting Scrambled Method of Moments simply replaces the random number generator with the scramble (available in most softwares) to reduce simulation noise. Scrambled Indirect Inference estimation is also considered. For time series, qMC may not apply directly because of a curse of dimensionality on the time dimension. A simple algorithm and a class of moments which circumvent this issue are described. Asymptotic results are given for each algorithm. Monte-Carlo examples illustrate these results in finite samples, including an income process with "lots of heterogeneity."


title: New constraints on Heavy Neutral Leptons from Super-Kamiokande data

abstract: Heavy neutral leptons are predicted in many extensions of the Standard Model with massive neutrinos. If kinematically accessible, they can be copiously produced from kaon and pion decays in atmospheric showers, and subsequently decay inside large neutrino detectors. We perform a search for these long-lived particles using Super-Kamiokande multi-GeV neutrino data and derive stringent limits on the mixing with electron, muon and tau neutrinos as a function of the long-lived particle mass. We also present the limits on the branching ratio versus lifetime plane, which are helpful in determining the constraints in non-minimal models where the heavy neutral leptons have new interactions with the Standard Model.


title: From Boundary Data to Bound States II: Scattering Angle to Dynamical Invariants (with Twist)

abstract: We recently introduced in [1910.03008] a "boundary-to-bound" dictionary between gravitational scattering data and observables for bound states of non-spinning bodies. In this paper, we elaborate further on this (holographic) map. We start by deriving the following -- remarkably simple -- formula relating the periastron advance to the scattering angle: $\Delta \Phi(J,{\cal E}) =\chi(J,{\cal E}) + \chi (-J,{\cal E})$, via analytic continuation in angular momentum and binding energy. Using explicit expressions from [1910.03008], we confirm its validity to all orders in the Post-Minkowskian (PM) expansion. Furthermore, we reconstruct the radial action for the bound state directly from the knowledge of the scattering angle. The radial action enables us to write compact expressions for dynamical invariants in terms of the deflection angle to all PM orders, which can also be written as a function of the PM-expanded amplitude. As an example, we reproduce our result in [1910.03008] for the periastron advance, and compute the radial and azimuthal frequencies and redshift variable to two-loops. Agreement is found in the overlap between PM and Post-Newtonian (PN) schemes. Last but not least, we initiate the study of our dictionary including spin. We demonstrate that the same relation between deflection angle and periastron advance applies for aligned-spin contributions, with $J$ the (canonical) total angular momentum. Explicit checks are performed to display perfect agreement using state-of-the-art PN results in the literature. Using the map between test- and two-body dynamics, we also compute the periastron advance up to quadratic order in the spin, to one-loop and to all orders in velocity. We conclude with a discussion on the generalized "impetus formula" for spinning bodies and black holes as "elementary particles".


title: The Orbit of WASP-12b is Decaying

abstract: WASP-12b is a transiting hot Jupiter on a 1.09-day orbit around a late-F star. Since the planet's discovery in 2008, the time interval between transits has been decreasing by $29\pm 2$ msec year$^{-1}$. This is a possible sign of orbital decay, although the previously available data left open the possibility that the planet's orbit is slightly eccentric and is undergoing apsidal precession. Here, we present new transit and occultation observations that provide more decisive evidence for orbital decay, which is favored over apsidal precession by a $\Delta\mathrm{BIC}$ of 22.3 or Bayes factor of 70,000. We also present new radial-velocity data that rule out the R{\o}mer effect as the cause of the period change. This makes WASP-12 the first planetary system for which we can be confident that the orbit is decaying. The decay timescale for the orbit is $P/\dot{P} = 3.25\pm 0.23$ Myr. Interpreting the decay as the result of tidal dissipation, the modified stellar tidal quality factor is $Q'_\star = 1.8 \times10^{5}$.