title: Stop memorizing: A data-dependent regularization framework for intrinsic pattern learning

abstract: Deep neural networks (DNNs) typically have enough capacity to fit random data by brute force even when conventional data-dependent regularizations focusing on the geometry of the features are imposed. We find out that the reason for this is the inconsistency between the enforced geometry and the standard softmax cross entropy loss. To resolve this, we propose a new framework for data-dependent DNN regularization, the Geometrically-Regularized-Self-Validating neural Networks (GRSVNet). During training, the geometry enforced on one batch of features is simultaneously validated on a separate batch using a validation loss consistent with the geometry. We study a particular case of GRSVNet, the Orthogonal-Low-rank Embedding (OLE)-GRSVNet, which is capable of producing highly discriminative features residing in orthogonal low-rank subspaces. Numerical experiments show that OLE-GRSVNet outperforms DNNs with conventional regularization when trained on real data. More importantly, unlike conventional DNNs, OLE-GRSVNet refuses to memorize random data or random labels, suggesting it only learns intrinsic patterns by reducing the memorizing capacity of the baseline DNN.


title: A $q$-extension of a partial differential equation and the Hahn polynomials

abstract: Using the theory of analytic functions of several complex variables, we prove that if an analytic function in several variables satisfies a system of $q$-partial differential equations, then, it can be expanded in terms of the product of the homogeneous Hahn polynomials. Some non-trivial applications of this expansion theorem to $q$-series are discussed.


title: Dynamics and level statistics of interacting fermions in the Lowest Landau Level

abstract: We consider the unitary dynamics of interacting fermions in the lowest Landau level, on spherical and toroidal geometries. The dynamics are driven by the interaction Hamiltonian which, viewed in the basis of single-particle Landau orbitals, contains correlated pair hopping terms in addition to static repulsion. This setting and this type of Hamiltonian has a significant history in numerical studies of fractional quantum Hall (FQH) physics, but the many-body quantum dynamics generated by such correlated hopping has not been explored in detail. We focus on initial states containing all the fermions in one block of orbitals. We characterize in detail how the fermionic liquid spreads out starting from such a state. We identify and explain differences with regular (single-particle) hopping Hamiltonians. Such differences are seen, e.g. in the entanglement dynamics, in that some initial block states are frozen or near-frozen, and in density gradients persisting in long-time equilibrated states. Examining the level spacing statistics, we show that the most common Hamiltonians used in FQH physics are not integrable, and explain that GOE statistics (level statistics corresponding to the Gaussian Orthogonal Ensemble) can appear in many cases despite the lack of time-reversal symmetry.


title: Distributed Computation in Node-Capacitated Networks

abstract: In this paper, we study distributed graph algorithms in networks in which the nodes have a limited communication capacity. Many distributed systems are built on top of an underlying networking infrastructure, for example by using a virtual communication topology known as an overlay network. Although this underlying network might allow each node to directly communicate with a large number of other nodes, the amount of communication that a node can perform in a fixed amount of time is typically much more limited. We introduce the Node-Capacitated Clique model as an abstract communication model, which allows us to study the effect of nodes having limited communication capacity on the complexity of distributed graph computations. In this model, the $n$ nodes of a network are connected as a clique and communicate in synchronous rounds. In each round, every node can exchange messages of $O(\log n)$ bits with at most $O(\log n)$ other nodes. When solving a graph problem, the input graph $G$ is defined on the same set of $n$ nodes, where each node knows which other nodes are its neighbors in $G$. To initiate research on the Node-Capacitated Clique model, we present distributed algorithms for the Minimum Spanning Tree (MST), BFS Tree, Maximal Independent Set, Maximal Matching, and Vertex Coloring problems. We show that even with only $O(\log n)$ concurrent interactions per node, the MST problem can still be solved in polylogarithmic time. In all other cases, the runtime of our algorithms depends linearly on the arboricity of $G$, which is a constant for many important graph families such as planar graphs.


title: Cross-domain attribute representation based on convolutional neural network

abstract: In the problem of domain transfer learning, we learn a model for the predic-tion in a target domain from the data of both some source domains and the target domain, where the target domain is in lack of labels while the source domain has sufficient labels. Besides the instances of the data, recently the attributes of data shared across domains are also explored and proven to be very helpful to leverage the information of different domains. In this paper, we propose a novel learning framework for domain-transfer learning based on both instances and attributes. We proposed to embed the attributes of dif-ferent domains by a shared convolutional neural network (CNN), learn a domain-independent CNN model to represent the information shared by dif-ferent domains by matching across domains, and a domain-specific CNN model to represent the information of each domain. The concatenation of the three CNN model outputs is used to predict the class label. An iterative algo-rithm based on gradient descent method is developed to learn the parameters of the model. The experiments over benchmark datasets show the advantage of the proposed model.


title: Quadrature Strategies for Constructing Polynomial Approximations

abstract: Finding suitable points for multivariate polynomial interpolation and approximation is a challenging task. Yet, despite this challenge, there has been tremendous research dedicated to this singular cause. In this paper, we begin by reviewing classical methods for finding suitable quadrature points for polynomial approximation in both the univariate and multivariate setting. Then, we categorize recent advances into those that propose a new sampling approach and those centered on an optimization strategy. The sampling approaches yield a favorable discretization of the domain, while the optimization methods pick a subset of the discretized samples that minimize certain objectives. While not all strategies follow this two-stage approach, most do. Sampling techniques covered include subsampling quadratures, Christoffel, induced and Monte Carlo methods. Optimization methods discussed range from linear programming ideas and Newton's method to greedy procedures from numerical linear algebra. Our exposition is aided by examples that implement some of the aforementioned strategies.


title: General solutions for nonlinear differential equations: a rule-based self-learning approach using deep reinforcement learning

abstract: A universal rule-based self-learning approach using deep reinforcement learning (DRL) is proposed for the first time to solve nonlinear ordinary differential equations and partial differential equations. The solver consists of a deep neural network-structured actor that outputs candidate solutions, and a critic derived only from physical rules (governing equations and boundary and initial conditions). Solutions in discretized time are treated as multiple tasks sharing the same governing equation, and the current step parameters provide an ideal initialization for the next owing to the temporal continuity of the solutions, which shows a transfer learning characteristic and indicates that the DRL solver has captured the intrinsic nature of the equation. The approach is verified through solving the Schr\"odinger, Navier-Stokes, Burgers', Van der Pol, and Lorenz equations and an equation of motion. The results indicate that the approach gives solutions with high accuracy, and the solution process promises to get faster.


title: Stochastic effects in autoimmune dynamics

abstract: Among various possible causes of autoimmune disease, an important role is played by infections that can result in a breakdown of immune tolerance, primarily through the mechanism of "molecular mimicry". In this paper we propose and analyse a stochastic model of immune response to a viral infection and subsequent autoimmunity, with account for the populations of T cells with different activation thresholds, regulatory T cells, and cytokines. We show analytically and numerically how stochasticity can result in sustained oscillations around deterministically stable steady states, and we also investigate stochastic dynamics in the regime of bi-stability. These results provide a possible explanation for experimentally observed variations in the progression of autoimmune disease. Computations of the variance of stochastic fluctuations provide practically important insights into how the size of these fluctuations depends on various biological parameters, and this also gives a headway for comparison with experimental data on variation in the observed numbers of T cells and organ cells affected by infection.


title: On the stochastic Lie algebra

abstract: We study the structure of the Lie algebra $\mathfrak{s}(n,\mathbb R)$ corresponding to the so-called stochastic Lie group $\mathcal{S} (n,\mathbb R)$. We obtain the Levi decomposition of the Lie algebra, classify Levi factor and classify the representation of the factor in $\mathbb{R}^n$. We discuss isomorphism of $\mathcal{S}(n,\mathbb R)$ with the group of invertible affine maps ${\it Aff}(n-1,\mathbb R)$. We prove that $\mathfrak s(n, \mathbb R)$ is generated by two generic elements.


title: Multitaper Spectral Estimation HDP-HMMs for EEG Sleep Inference

abstract: Electroencephalographic (EEG) monitoring of neural activity is widely used for sleep disorder diagnostics and research. The standard of care is to manually classify 30-second epochs of EEG time-domain traces into 5 discrete sleep stages. Unfortunately, this scoring process is subjective and time-consuming, and the defined stages do not capture the heterogeneous landscape of healthy and clinical neural dynamics. This motivates the search for a data-driven and principled way to identify the number and composition of salient, reoccurring brain states present during sleep. To this end, we propose a Hierarchical Dirichlet Process Hidden Markov Model (HDP-HMM), combined with wide-sense stationary (WSS) time series spectral estimation to construct a generative model for personalized subject sleep states. In addition, we employ multitaper spectral estimation to further reduce the large variance of the spectral estimates inherent to finite-length EEG measurements. By applying our method to both simulated and human sleep data, we arrive at three main results: 1) a Bayesian nonparametric automated algorithm that recovers general temporal dynamics of sleep, 2) identification of subject-specific "microstates" within canonical sleep stages, and 3) discovery of stage-dependent sub-oscillations with shared spectral signatures across subjects.


title: Enhanced Pricing and Management of Bundled Insurance Risks with Dependence-aware Prediction using Pair Copula Construction

abstract: We propose a dependence-aware predictive modeling framework for multivariate risks stemmed from an insurance contract with bundling features - an important type of policy increasingly offered by major insurance companies. The bundling feature naturally leads to longitudinal measurements of multiple insurance risks, and correct pricing and management of such risks is of fundamental interest to financial stability of the macroeconomy. We build a novel predictive model that fully captures the dependence among the multivariate repeated risk measurements. Specifically, the longitudinal measurement of each individual risk is first modeled using pair copula construction with a D-vine structure, and the multiple D-vines are then integrated by a flexible copula. The proposed model provides a unified modeling framework for multivariate longitudinal data that can accommodate different scales of measurements, including continuous, discrete, and mixed observations, and thus can be potentially useful for various economic studies. A computationally efficient sequential method is proposed for model estimation and inference, and its performance is investigated both theoretically and via simulation studies. In the application, we examine multivariate bundled risks in multi-peril property insurance using proprietary data from a commercial property insurance provider. The proposed model is found to provide improved decision making for several key insurance operations. For underwriting, we show that the experience rate priced by the proposed model leads to a 9% lift in the insurer's net revenue. For reinsurance, we show that the insurer underestimates the risk of the retained insurance portfolio by 10% when ignoring the dependence among bundled insurance risks.


title: Listening to the Universe through Indirect Detection

abstract: Indirect detection is the search for the particle nature of dark matter with astrophysical probes. Manifestly, it exists right at the intersection of particle physics and astrophysics, and the discovery potential for dark matter can be greatly extended using insights from both disciplines. This thesis provides an exploration of this philosophy. On the one hand, I will show how astrophysical observations of dark matter, through its gravitational interaction, can be exploited to determine the most promising locations on the sky to observe a particle dark matter signal. On the other, I demonstrate that refined theoretical calculations of the expected dark matter interactions can be used disentangle signals from astrophysical backgrounds. Both of these approaches will be discussed in the context of general searches, but also applied to the case of an excess of photons observed at the center of the Milky Way. This galactic center excess represents both the challenges and joys of indirect detection. Initially thought to be a signal of annihilating dark matter at the center of our own galaxy, it now appears more likely to be associated with a population of millisecond pulsars. Yet these pulsars were completely unanticipated, and highlight that indirect detection can lead to many new insights about the universe, hopefully one day including the particle nature of dark matter.


title: The nature of correlations in the insulating states of twisted bilayer graphene

abstract: The recently observed superconductivity in twisted bilayer graphene emerges from insulating states believed to arise from electronic correlations. While there have been many proposals to explain the insulating behaviour, the commensurability at which these states appear suggests that they are Mott insulators. Here we focus on the insulating states with $\pm 2$ electrons or holes with respect to the charge neutrality point. We show that the theoretical expectations for the Mott insulating states are not compatible with the experimentally observed dependence on temperature and magnetic field if, as frequently assumed, only the correlations between electrons on the same site are included. We argue that the inclusion of non-local (inter-site) correlations in the treatment of the Hubbard model can bring the predictions for the magnetic and temperature dependencies of the Mott transition to an agreement with experiments and have consequences for the critical interactions, the size of the gap, and possible pseudogap physics. The importance of the inter-site correlations to explain the experimental observations indicates that the observed insulating gap is not the one between the Hubbard bands and that antiferromagnetic-like correlations play a key role in the Mott transition.


title: Self-assembly of a drop pattern from a two-dimensional grid of nanometric metallic filaments

abstract: We report experiments, modeling and numerical simulations of the self--assembly of particle patterns obtained from a nanometric metallic square grid. Initially, nickel filaments of rectangular cross section are patterned on a SiO$_2$ flat surface, and then they are melted by laser irradiation with $\sim 20$ ns pulses. During this time, the liquefied metal dewets the substrate, leading to a linear array of drops along each side of the squares. The experimental data provides a series of SEM images of the resultant morphology as a function of the number of laser pulses or cumulative liquid lifetime. These data are analyzed in terms of fluid mechanical models that account for mass conservation and consider flow evolution with the aim to predict the final number of drops resulting from each side of the square. The aspect ratio, $\delta$, between the square sides' lengths and their widths is an essential parameter of the problem. Our models allow us to predict the $\delta$-intervals within which a certain final number of drops are expected. The comparison with experimental data shows a good agreement with the model that explicitly considers the Stokes flow developed in the filaments neck region that lead to breakup points. Also, numerical simulations, that solve the Navier-Stokes equations along with slip boundary condition at the contact lines, are implemented to describe the dynamics of the problem.


title: Cloud-Cloud Collision Induced Star Formation in IRAS 18223-1243

abstract: In the direction of l = 17.6 - 19 deg, the star-forming sites Sh 2-53 and IRAS 18223-1243 are prominently observed, and seem to be physically detached from each other. Sh 2-53 has been investigated at the junction of the molecular filaments, while a larger-scale environment of IRAS 18223-1243 remains unexplored. The goal of this paper is to investigate the star formation processes in the IRAS site (area ~0.4 deg x 0.4 deg). Based on the GRS 13CO line data, two molecular clouds, peaking at velocities of 45 and 51 km/s, are found. In the position-velocity plots, a relatively weak 13CO emission is detected at intermediate velocities (i.e. 47.5 - 49.5 km/s) between these two clouds, illustrating a link between two parallel elongated velocity structures. These clouds are physically connected in both space and velocity. The MAGPIS data at 20 cm trace free-free continuum emission toward the IRAS 18223-1243 source. Using the Spitzer and UKIDSS photometric data, we have identified infrared-excess young stellar objects (YSOs), and have observed their groups toward the intersection zones of the clouds. IRAS 18223-1243 is also spatially seen at an interface of the clouds. Considering these observational findings, we propose the onset of the collision of two clouds in the IRAS site about 1 Myr ago, which triggered the birth of massive star(s) and the YSO groups. A non-uniform distribution of the GPIPS H-band starlight mean polarization angles is also observed toward the colliding interfaces, indicating the impact of the collision on the magnetic field morphology.


title: Unitarity constraints on general scalar couplings with SARAH

abstract: We present an update of the Mathematica package SARAH to calculate unitarity constraints in BSM models. The new functions can perform an analytical and numerical calculation of the two-particle scattering matrix of (uncoloured) scalars. We do not make use of the simplifying assumption of a very large scattering energy, but include all contributions which could become important at small energies above the weak scale. This allows us to constrain trilinear scalar couplings. However, it can also modify (weakening or strengthening) the constraints on quartic couplings, which we show via the example of a singlet extended Standard Model.


title: Flux-driven Josephson parametric amplifier for sub-GHz frequencies fabricated with side-wall passivated spacer junction technology

abstract: We present experimental results on a Josephson parametric amplifier tailored for readout of ultra-sensitive thermal microwave detectors. In particular, we discuss the impact of fabrication details on the performance. We show that the small volume of deposited dielectric materials enabled by the side-wall passivated spacer niobium junction technology leads to robust operation across a wide range of operating temperatures up to 1.5 K. The flux-pumped amplifier has gain in excess of 20 dB in three-wave mixing and its center frequency is tunable between 540 MHz and 640 MHz. At 600 MHz, the amplifier adds 105 mK $\pm$ 9 mK of noise, as determined with the hot/cold source method. Phase-sensitive amplification is demonstrated with the device.


title: The structure of the space of ergodic measures of transitive partially hyperbolic sets

abstract: We provide examples of transitive partially hyperbolic dynamics (specific but paradigmatic examples of homoclinic classes) which blend different types of hyperbolicity in the one-dimensional center direction. These homoclinic classes have two disjoint parts: an "exposed" piece which is poorly homoclinically related with the rest and a "core" with rich homoclinic relations. There is an associated natural division of the space of ergodic measures which are either supported on the exposed piece or on the core. We describe the topology of these two parts and show that they glue along nonhyperbolic measures. Measures of maximal entropy are discussed in more detail. We present examples where the measure of maximal entropy is nonhyperbolic. We also present examples where the measure of maximal entropy is unique and nonhyperbolic, however in this case the dynamics is nontransitive.


title: Unitarity constraints in triplet extensions beyond the large s limit

abstract: Triplet extensions are attractive alternatives to the standard model (SM) of particle physics. While models with only one triplet are highly constrained by electroweak precision observables, this is not necessarily the case once several triplets are present as in the Georgi-Machacek model. As in all other BSM models, the parameter space of triplet extensions is constrained by the condition that perturbative unitarity is not violated. For this purpose, limits on the eigenvalues of the scalar $2 \to 2$ scattering matrix are set. It is very common in the BSM literature that the scattering matrix is calculated under one crucial assumption: the scattering energy s is so large that only point interactions involving quartic couplings provide non-negligible contributions. However, it is not given that this approximation is always valid - in fact, diagrams involving propagators can play an important role. We discuss at the examples of (i) the SM model extended by a real triplet, (ii) the Y = 1 triplet extension of the SM, and (iii) the Georgi-Machacek model, how the tree-level unitarity constraints are affected once the large s approximation is given up. For all models we find that the impact of (effective) cubic couplings can be crucial.


title: Improved unitarity constraints in Two-Higgs-Doublet-Models

abstract: Two-Higgs-Doublet-Models (THDMs) are among the simplest extensions of the standard model and are intensively studied in the literature. Using on-shell parameters such as the masses of the additional scalars as input, corresponds often to large quartic couplings in the underlying Lagrangian. Therefore, it is important to check if these couplings are for instance in agreement with perturbative unitarity. The common approach for doing this check is to consider the two-particle scattering matrix of scalars in the large centre-of-mass energy limit where only point interactions contribute. We show that this is not always a valid approximation: the full calculation including all tree-level contributions at finite energy can lead to much more stringent constraints. We show how the allowed regions in the parameter space are affected. In particular, the light Higgs window with a second Higgs below 125 GeV completely closes for large values of the $Z_2$ breaking parameter $|M_{12}|$. We also compare against the loop corrected constraints, which use also the large $\sqrt{s}$ approximation, and find that (effective) cubic couplings are often more important than radiative corrections.


title: Blended Conditional Gradients: the unconditioning of conditional gradients

abstract: We present a blended conditional gradient approach for minimizing a smooth convex function over a polytope P, combining the Frank--Wolfe algorithm (also called conditional gradient) with gradient-based steps, different from away steps and pairwise steps, but still achieving linear convergence for strongly convex functions, along with good practical performance. Our approach retains all favorable properties of conditional gradient algorithms, notably avoidance of projections onto P and maintenance of iterates as sparse convex combinations of a limited number of extreme points of P. The algorithm is lazy, making use of inexpensive inexact solutions of the linear programming subproblem that characterizes the conditional gradient approach. It decreases measures of optimality (primal and dual gaps) rapidly, both in the number of iterations and in wall-clock time, outperforming even the lazy conditional gradient algorithms of [arXiv:1410.8816]. We also present a streamlined version of the algorithm for the probability simplex.


title: Path length distribution in two-dimensional causal sets

abstract: We study the distribution of maximal-chain lengths between two elements of a causal set and its relationship with the embeddability of the causal set in a region of flat spacetime. We start with causal sets obtained from uniformly distributed points in Minkowski space. After some general considerations we focus on the 2-dimensional case and derive a recursion relation for the expected number of maximal chains $n_k$ as a function of their length $k$ and the total number of points $N$ between the maximal and minimal elements. By studying these theoretical distributions as well as ones generated from simulated sprinklings in Minkowski space we identify two features, the most probable path length or peak of the distribution $k_0$ and its width $\Delta$, which can be used both to provide a measure of the embeddability of the causal set as a uniform distribution of points in Minkowski space and to determine its dimensionality, if the causal set is manifoldlike in that sense. We end with a few simple examples of $n_k$ distributions for non-manifoldlike causal sets.


title: Transport coefficients of hot magnetized QCD matter beyond the lowest Landau level approximation

abstract: In this article, shear viscosity, bulk viscosity, and thermal conductivity of a hot QCD medium have been studied in the presence of a strong magnetic field. To model the hot magnetized QCD matter, an extended quasi-particle description of the hot QCD equation of state in the presence of the magnetic field has been adopted. The effects of higher Landau levels on the temperature dependence of viscous coefficients (bulk and shear viscosities) and thermal conductivity have been obtained by considering the $1\rightarrow 2$ processes in the presence of the strong magnetic field. An effective covariant kinetic theory has been set up in (1+1)-dimensional that includes mean field contributions in terms of quasi-particle dispersions and magnetic field to describe the Landau level dynamics of quarks. The sensitivity of these parameters to the magnitude of the magnetic field has also been explored. Both the magnetic field and mean field contributions have seen to play a significant role in obtaining the temperature behaviour of the transport coefficients of hot QCD medium.


title: Efficient Topological Materials Discovery Using Symmetry Indicators

abstract: Although the richness of spatial symmetries has led to a rapidly expanding inventory of possible topological crystalline (TC) phases of electrons, physical realizations have been slow to materialize due to the practical difficulty to ascertaining band topology in realistic calculations. Here, we integrate the recently established theory of symmetry indicators of band topology into first-principle band-structure calculations, and test it on a databases of previously synthesized crystals. The combined algorithm is found to efficiently unearth topological materials and predict topological properties like protected surface states. On applying our algorithm to just 8 out of the 230 space groups, we already discover numerous materials candidates displaying a diversity of topological phenomena, which are simultaneously captured in a single sweep. The list includes recently proposed classes of TC insulators that had no previous materials realization as well as other topological phases, including: (i) a screw-protected 3D TC insulator, \b{eta}-MoTe2, with gapped surfaces except for 1D helical "hinge" states; (ii) a rotation-protected TC insulator BiBr with coexisting surface Dirac cones and hinge states; (iii) non-centrosymmetric Z2 topological insulators undetectable using the well-established parity criterion, AgXO (X=Na,K,Rb); (iv) a Dirac semimetal MgBi2O6; (v) a Dirac nodal-line semimetal AgF2; and (vi) a metal with three-fold degenerate band crossing near the Fermi energy, AuLiMgSn. Our work showcases how the recent theoretical insights on the fundamentals of band structures can aid in the practical goal of discovering new topological materials.


title: Theory of Interacting Cavity Rydberg Polaritons

abstract: Photonic materials are an emerging platform to explore quantum matter and quantum dynamics. The development of Rydberg electromagnetically induced transparency provided a clear route to strong interactions between individual optical photons. In conjunction with carefully designed optical resonators, it is now possible to achieve extraordinary control of the properties of individual photons, introducing tunable gauge fields whilst imbuing the photons with mass and embedding them on curved spatial manifolds. Building on work formalizing Rydberg-mediated interactions between propagating photons, we develop a theory of interacting Rydberg polaritons in multimode optical resonators, where the strong interactions are married with tunable single-particle properties to build and probe exotic matter. In the presence of strong coupling between the resonator field and a Rydberg-dressed atomic ensemble, a quasiparticle called the "cavity Rydberg polariton" emerges. We investigate its properties, finding that it inherits both the fast dynamics of its photonic constituents and the strong interactions of its atomic constituents. We develop tools to properly renormalize the interactions when polaritons approach each other, and investigate the impact of atomic motion on the coherence of multi-mode polaritons, showing that most channels for atom-polariton cross-thermalization are strongly suppressed. Finally, we propose to harness the repeated diffraction and refocusing of the optical resonator to realize interactions which are local in momentum space. This work points the way to efficient modeling of polaritonic quantum materials in properly renormalized strongly interacting effective theories, thereby enabling experimental studies of photonic fractional quantum Hall fluids and crystals, plus photonic quantum information processors and repeaters.


title: Effects of viral and cytokine delays on dynamics of autoimmunity

abstract: A major contribution to the onset and development of autoimmune disease is known to come from infections. An important practical problem is identifying the precise mechanism by which the breakdown of immune tolerance as a result of immune response to infection leads to autoimmunity. In this paper, we develop a mathematical model of immune response to a viral infection, which includes T cells with different activation thresholds, regulatory T cells (Tregs), and~a cytokine mediating immune dynamics. Particular emphasis is made on the role of time delays associated with the processes of infection and mounting the immune response. Stability analysis of various steady states of the model allows us to identify parameter regions associated with different types of immune behaviour, such as, normal clearance of infection, chronic infection, and autoimmune dynamics. Numerical simulations are used to illustrate different dynamical regimes, and to identify basins of attraction of different dynamical states. An important result of the analysis is that not only the parameters of the system, but also the initial level of infection and the initial state of the immune system determine the progress and outcome of the dynamics.


title: Efficient Exploration of Gradient Space for Online Learning to Rank

abstract: Online learning to rank (OL2R) optimizes the utility of returned search results based on implicit feedback gathered directly from users. To improve the estimates, OL2R algorithms examine one or more exploratory gradient directions and update the current ranker if a proposed one is preferred by users via an interleaved test. In this paper, we accelerate the online learning process by efficient exploration in the gradient space. Our algorithm, named as Null Space Gradient Descent, reduces the exploration space to only the \emph{null space} of recent poorly performing gradients. This prevents the algorithm from repeatedly exploring directions that have been discouraged by the most recent interactions with users. To improve sensitivity of the resulting interleaved test, we selectively construct candidate rankers to maximize the chance that they can be differentiated by candidate ranking documents in the current query; and we use historically difficult queries to identify the best ranker when tie occurs in comparing the rankers. Extensive experimental comparisons with the state-of-the-art OL2R algorithms on several public benchmarks confirmed the effectiveness of our proposal algorithm, especially in its fast learning convergence and promising ranking quality at an early stage.


title: The fraction of ionizing radiation from massive stars that escapes to the intergalactic medium

abstract: The part played by stars in the ionization of the intergalactic medium remains an open question. A key issue is the proportion of the stellar ionizing radiation that escapes the galaxies in which it is produced. Spectroscopy of gamma-ray burst afterglows can be used to determine the neutral hydrogen column-density in their host galaxies and hence the opacity to extreme ultra-violet radiation along the lines-of-sight to the bursts. Thus, making the reasonable assumption that long-duration GRB locations are representative of the sites of massive stars that dominate EUV production, one can calculate an average escape fraction of ionizing radiation in a way that is independent of galaxy size, luminosity or underlying spectrum. Here we present a sample of NH measures for 138 GRBs in the range 1.6<z<6.7 and use it to establish an average escape fraction at the Lyman limit of <fesc>~0.005, with a 98% confidence upper limit of ~0.015. This analysis suggests that stars provide a small contribution to the ionizing radiation budget of the IGM at z<5, where the bulk of the bursts lie. At higher redshifts, z>5, firm conclusions are limited by the small size of the GRB sample, but any decline in average HI column-density seems to be modest. We also find no indication of a significant correlation of NH with galaxy UV luminosity or host stellar mass, for the subset of events for which these are available. We discuss in some detail a number of selection effects and potential biases. Drawing on a range of evidence we argue that such effects, while not negligible, are unlikely to produce systematic errors of more than a factor ~2, and so would not affect the primary conclusions. Given that many GRB hosts are low metallicity, high specific star-formation rate, dwarf galaxies, these results present a particular problem for the hypothesis that such galaxies dominated the reionization of the universe.


title: Mixup-Based Acoustic Scene Classification Using Multi-Channel Convolutional Neural Network

abstract: Audio scene classification, the problem of predicting class labels of audio scenes, has drawn lots of attention during the last several years. However, it remains challenging and falls short of accuracy and efficiency. Recently, Convolutional Neural Network (CNN)-based methods have achieved better performance with comparison to the traditional methods. Nevertheless, conventional single channel CNN may fail to consider the fact that additional cues may be embedded in the multi-channel recordings. In this paper, we explore the use of Multi-channel CNN for the classification task, which aims to extract features from different channels in an end-to-end manner. We conduct the evaluation compared with the conventional CNN and traditional Gaussian Mixture Model-based methods. Moreover, to improve the classification accuracy further, this paper explores the using of mixup method. In brief, mixup trains the neural network on linear combinations of pairs of the representation of audio scene examples and their labels. By employing the mixup approach for data argumentation, the novel model can provide higher prediction accuracy and robustness in contrast with previous models, while the generalization error can also be reduced on the evaluation data.


title: Reentrant cluster glass and stability of ferromagnetism in Ga2MnCo Heusler alloy

abstract: We present here a detailed investigation into the magnetic ordering of full Heusler alloy Ga$_2$MnCo using dc, ac magnetization measurements, neutron diffraction and neutron depolarization experiments. Crystal structure at room temperature was first confirmed to be L2$_1$ using the highly intense synchrotron X-ray diffraction (XRD) technique. Temperature dependent magnetization reveals that Ga$_2$MnCo enters a ferromagnetic (FM) state at $T_C = $154 K, characterized by a sharp increase in magnetization and a plateau-like region hereafter. As the temperature is decreased further, a sharp drop in magnetization is observed at $T_f$ = 50 K, hinting towards an antiferromagnetic (AFM) phase change. Neutron diffraction (ND) recorded over the range of temperature from 6 to 300 K, provides combined information regarding crystal as well as magnetic structure. Accordingly, an increase in the intensity of the ND pattern is seen at 150 K, signaling onset of long range FM order. However, there is no sign of appearance of superlattice reflections corresponding to the AFM phase, in the patterns recorded below 50 K. An unusual discontinuity in the unit cell volume is seen around $T_f$ indicating a coupling of this second transition with the contraction of the lattice. Attempts to unravel this interesting magnetic behaviour using ac susceptibility measurements lead to the existence of glassy magnetism below $T_f$. Systematic analysis of the susceptibility results along with neutron depolarization measurement, identifies the low temperature phase as a reentrant cluster glass.


title: Global asymptotic stability of bifurcating, positive equilibria of p-Laplacian boundary value problems with p-concave nonlinearities

abstract: We consider the parabolic, initial value problem $$ v_t =\Delta_p(v)+\lambda g(x,v)\phi_p(v), \quad \text{in $\Omega \times (0,\infty),$} $$ \[ v =0, \text{in $\partial\Omega \times (0,\infty),$}\tag{IVP} v =v_0\ge0, \text{in $\Omega \times \{0\},$} \] where $\Omega$ is a bounded domain in ${\mathbb R}^N$, for some integer $N\ge1$, with smooth boundary $\partial\Omega$, $\phi_p(s):=|s|^{p-1} {\rm sgn}s$, $s\in{\mathbb R}$, $\Delta_p$ denotes the $p$-Laplacian, with $p>\max\{2,N\}$, $v_0\in C^0(\overline{\Omega})$, and $\lambda>0$. The function $g:\overline{\Omega } \times [0,\infty)\to(0,\infty)$ is $C^0$ and, for each $x\in\overline{\Omega }$, the function $g(x,\cdot):[0,\infty)\to(0,\infty)$ is Lipschitz continuous and strictly decreasing. Clearly, (IVP) has the trivial solution $v\equiv0$, for all $\lambda>0$. In addition, there exists $0<\lambda_{\rm min}(g)<\lambda_{\rm max}(g)$ ($\lambda_{\rm max}(g)$ may be $\infty$) such that: $(a)$ if $\lambda\not\in(\lambda_{\rm min}(g),\lambda_{\rm max}(g))$ then (IVP) has no non-trivial, positive equilibrium; $(b)$ if $\lambda\in(\lambda_{\rm min}(g),\lambda_{\rm max}(g))$ then (IVP) has a unique, non-trivial, positive equilibrium $e_\lambda\in W_0^{1,p}(\Omega)$. We prove the following results on the positive solutions of (IVP): $(a)$ if $0<\lambda<\lambda_{\rm min}(g)$ then the trivial solution is globally asymptotically stable; $(b)$ if $\lambda_{\rm min}(g)<\lambda<\lambda_{\rm max}(g)$ then $e_\lambda$ is globally asymptotically stable; $(c)$ if $\lambda_{\rm max}(g)<\lambda$ then any non-trivial solution blows up in finite time.


title: Large Deviations Properties of Maximum Entropy Markov Chains from Spike Trains

abstract: We consider the maximum entropy Markov chain inference approach to characterize the collective statistics of neuronal spike trains, focusing on the statistical properties of the inferred model. We review large deviations techniques useful in this context to describe properties of accuracy and convergence in terms of sampling size. We use these results to study the statistical fluctuation of correlations, distinguishability and irreversibility of maximum entropy Markov chains. We illustrate these applications using simple examples where the large deviation rate function is explicitly obtained for maximum entropy models of relevance in this field.


title: Four-Dimensional Topological Insulators with Nodal-Line Boundary States

abstract: Conventional topological insulators and superconductors have topologically protected nodal points on their boundaries, and the recent interests in nodal-line semimetals only concerned bulk band structures. Here, we present a novel four-dimensional topological insulator protected by an anti-unitary reflection symmetry, whose boundary band has a single $PT$-symmetric nodal line with double topological charges. Inspired by the recent experimental realization of the four-dimensional quantum Hall effect, we also propose a cold-atom system which realizes the novel topological insulator with tunable parameters as extra dimensions.


title: GANE: A Generative Adversarial Network Embedding

abstract: Network embedding has become a hot research topic recently which can provide low-dimensional feature representations for many machine learning applications. Current work focuses on either (1) whether the embedding is designed as an unsupervised learning task by explicitly preserving the structural connectivity in the network, or (2) whether the embedding is a by-product during the supervised learning of a specific discriminative task in a deep neural network. In this paper, we focus on bridging the gap of the two lines of the research. We propose to adapt the Generative Adversarial model to perform network embedding, in which the generator is trying to generate vertex pairs, while the discriminator tries to distinguish the generated vertex pairs from real connections (edges) in the network. Wasserstein-1 distance is adopted to train the generator to gain better stability. We develop three variations of models, including GANE which applies cosine similarity, GANE-O1 which preserves the first-order proximity, and GANE-O2 which tries to preserves the second-order proximity of the network in the low-dimensional embedded vector space. We later prove that GANE-O2 has the same objective function as GANE-O1 when negative sampling is applied to simplify the training process in GANE-O2. Experiments with real-world network datasets demonstrate that our models constantly outperform state-of-the-art solutions with significant improvements on precision in link prediction, as well as on visualizations and accuracy in clustering tasks.


title: Machine learning with force-field inspired descriptors for materials: fast screening and mapping energy landscape

abstract: We present a complete set of chemo-structural descriptors to significantly extend the applicability of machine-learning (ML) in material screening and mapping energy landscape for multicomponent systems. These new descriptors allow differentiating between structural prototypes, which is not possible using the commonly used chemical-only descriptors. Specifically, we demonstrate that the combination of pairwise radial, nearest neighbor, bond-angle, dihedral-angle and core-charge distributions plays an important role in predicting formation energies, bandgaps, static refractive indices, magnetic properties, and modulus of elasticity for three-dimensional (3D) materials as well as exfoliation energies of two-dimensional (2D) layered materials. The training data consists of 24549 bulk and 616 monolayer materials taken from JARVIS-DFT database. We obtained very accurate ML models using gradient boosting algorithm. Then we use the trained models to discover exfoliable 2D-layered materials satisfying specific property requirements. Additionally, we integrate our formation energy ML model with a genetic algorithm for structure search to verify if the ML model reproduces the DFT convex hull. This verification establishes a more stringent evaluation metric for the ML model than what commonly used in data sciences. Our learnt model is publicly available on the JARVIS-ML website (https://www.ctcms.nist.gov/jarvisml ) property predictions of generalized materials.


title: Transversal special parabolic points in the graph of a polynomial obtained under Viro's patchworking

abstract: In this article we focus on the study of special parabolic points in surfaces arising as graphs of polynomials, we give a theorem of Viro's patchworking type to build families of real polynomials in two variables with a prescribed number of special parabolic points in their graphs. We use this result to build a family of degree d real polynomials in two variables with $(d-4)(2d-9)$ special parabolic points in its graph. This brings the number of special parabolic points closer to the upper bound of $(d-2)(5d-12)$ when $d \geq 13$, which is the best known up until now.


title: Reconfiguring structured light beams using nonlinear metasurfaces

abstract: Ultra-compact, low-loss, fast, and reconfigurable optical components, enabling manipulation of light by light, could open numerous opportunities for controlling light on the nanoscale. Nanostructured all-dielectric metasurfaces have been shown to enable extensive control of amplitude and phase of light in the linear optical regime. Among other functionalities, they offer unique opportunities for shaping the wave front of light to introduce the orbital angular momentum (OAM) to a beam. Such structured light beams bring a new degree of freedom for applications ranging from spectroscopy and micromanipulation to classical and quantum optical communications. To date, reconfigurability or tuning of the optical properties of all-dielectric metasurfaces have been achieved mechanically, thermally, electrically or optically, using phase-change or nonlinear optical materials. However, a majority of demonstrated tuning approaches are either slow or require high optical powers. Arsenic trisulfide (As$_2$S$_3$) chalcogenide glass offering ultra-fast and large $\chi^{(3)}$ nonlinearity as well as a low two-photon absorption coefficient in the near and mid-wave infrared spectral range, could provide a new platform for the realization of fast and relatively low-intensity reconfigurable metasurfaces. Here, we design and experimentally demonstrate an As$_2$S$_3$ chalcogenide glass based metasurface that enables reshaping of a conventional Hermite-Gaussian beam with no OAM into an OAM beam at low-intensity levels, while preserves the original beam's amplitude and phase characteristics at high-intensity levels. The proposed metasurface could find applications for a new generation of optical communication systems and optical signal processing.


title: Ground-based optical transmission spectroscopy of the small, rocky exoplanet GJ 1132b

abstract: Terrestrial Solar System planets either have high mean molecular weight atmospheres, as with Venus, Mars, and Earth, or no atmosphere at all, as with Mercury. We do not have sufficient observational information to know if this is typical of terrestrial planets or a phenomenon unique to the Solar System. The bulk of atmospheric exoplanet studies have focused on hot Jupiters and Neptunes, but recent discoveries of small, rocky exoplanets transiting small, nearby stars provide targets that are amenable to atmospheric study. GJ 1132b has a radius of 1.2 Earth radii and a mass of 1.6 Earth masses, and orbits an M-dwarf 12 parsecs away from the Solar System. We present results from five transits of GJ 1132b taken with the Magellan Clay Telescope and the LDSS3C multi-object spectrograph. We jointly fit our five data sets when determining the best-fit transit parameters both for the white light curve and wavelength-binned light curves. We bin the light curves into 20 nm wavelength bands to construct the transmission spectrum. Our results disfavor a clear, 10x solar metallicity atmosphere at 3.7 sigma confidence and a 10% H2O, 90% H2 atmosphere at 3.5 sigma confidence. Our data are consistent with a featureless spectrum, implying that GJ 1132b has a high mean molecular weight atmosphere or no atmosphere at all, though we do not account for the possible presence of aerosols. This result is in agreement with theoretical work which suggests that a planet of GJ 1132b's mass and insolation should not be able to retain a H2 envelope.


title: New explicit solution to the N-Queens Problem and its relation to the Millennium Problem

abstract: Using modular arithmetic of the ring $\mathbb{Z}_{n+1}$ we obtain a new short solution to the problem of existence of at least one solution to the $N$-Queens problem on an $N \times N$ chessboard. It was proved, that these solutions can be represented as the Queen function with the width fewer or equal to 3. It is shown, that this estimate could not be reduced. A necessary and sufficient condition of being a composition of solutions a solution is found. Based on the obtained results we formulate a conjecture about the width of the representation of arbitrary solution. If this conjecture is valid, it entails solvability of the $N$-Queens completion in polynomial time. The connection between the $N$-Queens completion and the Millennium $P$ vs $NP$ Problem is found by the group of mathematicians from Scotland in August 2017.


title: Higher Codimensional Alpha Invariants and Characterization of Projective Spaces

abstract: We generalize the definition of alpha invariant to arbitrary codimension. We also give a lower bound of these alpha invariants for K-semistable Q-Fano varieties and show that we can characterize projective spaces among all K-semistable Fano manifolds in terms of higher codimensional alpha invariants. Our results demonstrate the relation between alpha invariants of any codimension and volumes of Fano manifolds in the characterization of projective spaces.


title: Positive and Unlabeled Learning through Negative Selection and Imbalance-aware Classification

abstract: Motivated by applications in protein function prediction, we consider a challenging supervised classification setting in which positive labels are scarce and there are no explicit negative labels. The learning algorithm must thus select which unlabeled examples to use as negative training points, possibly ending up with an unbalanced learning problem. We address these issues by proposing an algorithm that combines active learning (for selecting negative examples) with imbalance-aware learning (for mitigating the label imbalance). In our experiments we observe that these two techniques operate synergistically, outperforming state-of-the-art methods on standard protein function prediction benchmarks.


title: Origin of magnetoresistance suppression in thin $\gamma$-MoTe$_2$

abstract: We use both classical magnetotransport and quantum oscillation measurements to study the thickness evolution of the extremely large magnetoresistance (XMR) material and type-II Weyl semimetal candidate, $\gamma$-MoTe$_2$, protected from oxidation. We find that the magnetoresistance is systematically suppressed with reduced thickness. This occurs concomitantly with both a decrease in carrier mobility and increase in electron-hole imbalance. We model the two effects separately and conclude that the XMR effect is more sensitive to the former.


title: Optimal Power Control for Fading Channels with Arbitrary Input Distributions and Delay-Sensitive Traffic

abstract: This paper presents the optimal power control policies maximizing the effective capacity achieved with arbitrary input distributions subject to an average power constraint and quality of service (QoS) requirements. The analysis leads to simplified expressions for the optimal power control strategies in the low power regime and two limiting cases, i.e., extremely stringent QoS constraints and vanishing QoS constraints. In the low power regime, the energy efficiency (EE) performance with the constant-power scheme is also determined by characterizing both the minimum energy per bit and wideband slope for arbitrary input signaling and general fading distributions. Subsequently, the results are specialized to Nakagami-m and Rician fading channels. Also, tradeoff between the effective capacity and EE is studied by determining the optimal power control scheme that maximizes the effective capacity subject to constraints on the minimum required EE and average transmission power. Circuit power consumption is explicitly considered in the EE formulation. Through numerical results, the performance comparison between constant-power and optimal power control schemes for different signal constellations and Gaussian signals is carried out. The impact of QoS constraints, input distributions, fading severity, and average transmit power level on the proposed power control schemes, maximum achievable effective capacity and EE is evaluated.


title: Parameter-Independent Quark Mass Relation in the U(3)$\times$U(3)$'$ Model

abstract: Recently, we have proposed a quark mass matrix model based on U(3)$\times$U(3)$'$ family symmetry, in which up- and down-quark mass matrices $M_u$ and $M_d$ are described only by complex parameters $a_u $ and $a_d $, respectively. When we use charged lepton masses as additional input values, we can successfully obtain predictions for quark masses and Cabibbo-Kobayashi-Maskawa mixing. Since we have only one complex parameter $a_q$ for each mass matrix $M_q$, we can obtain a parameter-independent mass relation by using three equations for ${\rm Tr}[H_q]$, ${\rm Tr}[H_q H_q]$ and ${\rm det}H_q$, where $H_q \equiv M_q M_q^\dagger$ ($q=u, d$). In this paper, we investigate its parameter-independent feature of the quark mass relation in the model.


title: A new construction of the degree of maximal monotone maps

abstract: The inclusion equations of the type $f \in T ( x)$ where $T: X \to 2^{X^{\ast}}$ is a maximal monotone map, are extensively studied in nonlinear analysis. In this paper, we present a new construction of the degree of maximal monotone maps of the form $T: Y \to 2^{X^*}$, where $Y$ is a locally uniformly convex and separable Banach space continuously embedded in the uniformly convex Banach space $X$. The advantage of the new construction lies in the remarkable simplicity it offers for calculation of degree in comparison with the classical one suggested by F. Browder. We prove a few classical theorems in convex analysis through the suggested degree.


title: A Partially Inexact Alternating Direction Method of Multipliers and its Iteration-Complexity Analysis

abstract: This paper proposes a partially inexact alternating direction method of multipliers for computing approximate solution of a linearly constrained convex optimization problem. This method allows its first subproblem to be solved inexactly using a relative approximate criterion, whereas a proximal term is added to its second subproblem in order to simplify it. A stepsize parameter is included in the updating rule of the Lagrangian multiplier to improve its computational performance. Pointwise and ergodic interation-complexity bounds for the proposed method are established. To the best of our knowledge, this is the first time that complexity results for an inexact ADMM with relative error criteria has been analyzed. Some preliminary numerical experiments are reported to illustrate the advantages of the new method.


title: Reconstruction of training samples from loss functions

abstract: This paper presents a new mathematical framework to analyze the loss functions of deep neural networks with ReLU functions. Furthermore, as as application of this theory, we prove that the loss functions can reconstruct the inputs of the training samples up to scalar multiplication (as vectors) and can provide the number of layers and nodes of the deep neural network. Namely, if we have all input and output of a loss function (or equivalently all possible learning process), for all input of each training sample $x_i \in \mathbb{R}^n$, we can obtain vectors $x'_i\in \mathbb{R}^n$ satisfying $x_i=c_ix'_i$ for some $c_i \neq 0$. To prove theorem, we introduce the notion of virtual polynomials, which are polynomials written as the output of a node in a deep neural network. Using virtual polynomials, we find an algebraic structure for the loss surfaces, called semi-algebraic sets. We analyze these loss surfaces from the algebro-geometric point of view. Factorization of polynomials is one of the most standard ideas in algebra. Hence, we express the factorization of the virtual polynomials in terms of their active paths. This framework can be applied to the leakage problem in the training of deep neural networks. The main theorem in this paper indicates that there are many risks associated with the training of deep neural networks. For example, if we have N (the dimension of weight space) + 1 nonsmooth points on the loss surface, which are sufficiently close to each other, we can obtain the input of training sample up to scalar multiplication. We also point out that the structures of the loss surfaces depend on the shape of the deep neural network and not on the training samples.


title: Degree conditions for embedding trees

abstract: We conjecture that every $n$-vertex graph of minimum degree at least $\frac k2$ and maximum degree at least $2k$ contains all trees with $k$ edges as subgraphs. We prove an approximate version of this conjecture for trees of bounded degree and dense host graphs. Our work also has implications on the Erd\H os--S\'os conjecture and the $\frac 23$-conjecture. We prove an approximate version of both conjectures for bounded degree trees and dense host graphs.


title: Scanner: Efficient Video Analysis at Scale

abstract: A growing number of visual computing applications depend on the analysis of large video collections. The challenge is that scaling applications to operate on these datasets requires efficient systems for pixel data access and parallel processing across large numbers of machines. Few programmers have the capability to operate efficiently at these scales, limiting the field's ability to explore new applications that leverage big video data. In response, we have created Scanner, a system for productive and efficient video analysis at scale. Scanner organizes video collections as tables in a data store optimized for sampling frames from compressed video, and executes pixel processing computations, expressed as dataflow graphs, on these frames. Scanner schedules video analysis applications expressed using these abstractions onto heterogeneous throughput computing hardware, such as multi-core CPUs, GPUs, and media processing ASICs, for high-throughput pixel processing. We demonstrate the productivity of Scanner by authoring a variety of video processing applications including the synthesis of stereo VR video streams from multi-camera rigs, markerless 3D human pose reconstruction from video, and data-mining big video datasets such as hundreds of feature-length films or over 70,000 hours of TV news. These applications achieve near-expert performance on a single machine and scale efficiently to hundreds of machines, enabling formerly long-running big video data analysis tasks to be carried out in minutes to hours.


title: Improved Sentence Modeling using Suffix Bidirectional LSTM

abstract: Recurrent neural networks have become ubiquitous in computing representations of sequential data, especially textual data in natural language processing. In particular, Bidirectional LSTMs are at the heart of several neural models achieving state-of-the-art performance in a wide variety of tasks in NLP. However, BiLSTMs are known to suffer from sequential bias - the contextual representation of a token is heavily influenced by tokens close to it in a sentence. We propose a general and effective improvement to the BiLSTM model which encodes each suffix and prefix of a sequence of tokens in both forward and reverse directions. We call our model Suffix Bidirectional LSTM or SuBiLSTM. This introduces an alternate bias that favors long range dependencies. We apply SuBiLSTMs to several tasks that require sentence modeling. We demonstrate that using SuBiLSTM instead of a BiLSTM in existing models leads to improvements in performance in learning general sentence representations, text classification, textual entailment and paraphrase detection. Using SuBiLSTM we achieve new state-of-the-art results for fine-grained sentiment classification and question classification.


title: A Realistic U(2) Model of Flavor

abstract: We propose a simple $U(2)$ model of flavor compatible with an $SU(5)$ GUT structure. All hierarchies in fermion masses and mixings arise from powers of two small parameters that control the $U(2)$ breaking. In contrast to previous $U(2)$ models this setup can be realized without supersymmetry and provides an excellent fit to all SM flavor observables including neutrinos. We also consider a variant of this model based on a $D_6 \times U(1)_F$ flavor symmetry, which closely resembles the $U(2)$ structure, but allows for Majorana neutrino masses from the Weinberg operator. Remarkably, in this case one naturally obtains large mixing in the lepton sector from small mixing in the quark sector. The model also offers a natural option for addressing the Strong CP Problem and Dark Matter by identifying the Goldstone boson of the $U(1)_F$ factor as the QCD axion.


title: Gravitino Decay in High Scale Supersymmetry with R-parity Violation

abstract: We consider the effects of R-parity violation due to the inclusion of a bilinear $\mu^\prime L H_u$ superpotential term in high scale supersymmetric models with an EeV scale gravitino as dark matter. Although the typical phenomenological limits on this coupling (e.g. due to lepton number violation and the preservation of the baryon asymmetry) are relaxed when the supersymmetric mass spectrum is assumed to be heavy (in excess of the inflationary scale of $3 \times 10^{13}$ GeV), the requirement that the gravitino be sufficiently long-lived so as to account for the observed dark matter density, leads to a relatively strong bound on $\mu^\prime \lesssim 20$ GeV. The dominant decay channels for the longitudinal component of the gravitino are $Z \nu, W^\pm l^\mp$, and $h\nu$. To avoid an excess neutrino signal in IceCube, our limit on $\mu'$ is then strengthened to $\mu' \lesssim 50$ keV. When the bound is saturated, we find that there is a potentially detectable flux of mono-chromatic neutrinos with EeV energies.


title: The impact of binaural white noise with oscillations of 100 to 750hz in the short-term visual working memory and the reactivity of alpha and beta cerebral waves

abstract: According to some researchers, noise is typically conceived as a detrimental factor in cognitive performance affecting perception, decision making, and motor function. However, in recent studies it is associated with white noise with concentration and calm, therefore, this research seeks to establish the impact of binaural white noise on the performance of short-term visual and working memory, the alpha - beta brain activity, attention - meditation, through the use of two auditory stimuli with frequency ranges of (100 to 450hz) and (100 to 750hz). This study was conducted in the city of Montes Claros, the Republic of Brazil, where seven participants were evaluated (n = 7) with an average age of 36.71, and two age groups (GP1) 21 to 30 and (GP2) 41 50, with university studies. Within the experimental process, the short-term visual memory tests were performed using the cognitive assessment battery CAB of CogniFit, and the recording of brain activities through the use of monopolar electroencephalogram and the eSense algorithms. With the results obtained and through the use of statistical tests, we can infer that the binaural white noise with oscillations of 100 to 750 Hz contributed to the performance of visual work memory in the short term


title: Thermal corrections to the gluon magnetic Debye mass

abstract: We compute the gluon polarization tensor in a thermo-magnetic environment in the strong magnetic field limit at zero and high temperature. The magnetic field effects are introduced using Schwinger's proper time method. Thermal effects are computed in the HTL approximation. At zero temperature, we reproduce the well-known result whereby for a non-vanishing quark mass, the polarization tensor reduces to the parallel structure and its coefficient develops an imaginary part corresponding to the threshold for quark-antiquark pair production. This coefficient is infrared finite and simplifies considerably when the quark mass vanishes. Keeping always the field strength as the largest energy scale, in the high temperature regime we analyze two complementary hierarchies of scales: $q^2\ll m_f^2\ll T^2$ and $m_f^2\ll q^2\ll T^2$. In the latter, we show that the polarization tensor is infrared finite as $m_f$ goes to zero. In the former, we discuss the thermal corrections to the magnetic Debye mass.


title: Curves containing all points of a finite projective Galois plane

abstract: In the projective plane PG(2,q) over a finite field of order q, a Tallini curve is a plane irreducible (algebraic) curve of (minimum) degree q+2 containing all points of PG(2,q). Such curves were investigated by G. Tallini in 1961, and by Homma and Kim in 2013. Our results concern the automorphism groups, the Weierstrass semigroups, the Hasse-Witt invariants, and quotient curves of the Tallini curves.


title: Accurate Kernel Learning for Linear Gaussian Markov Processes using a Scalable Likelihood Computation

abstract: We report an exact likelihood computation for Linear Gaussian Markov processes that is more scalable than existing algorithms for complex models and sparsely sampled signals. Better scaling is achieved through elimination of repeated computations in the Kalman likelihood, and by using the diagonalized form of the state transition equation. Using this efficient computation, we study the accuracy of kernel learning using maximum likelihood and the posterior mean in a simulation experiment. The posterior mean with a reference prior is more accurate for complex models and sparse sampling. Because of its lower computation load, the maximum likelihood estimator is an attractive option for more densely sampled signals and lower order models. We confirm estimator behavior in experimental data through their application to speleothem data.


title: Organic chemistry in a CO2 rich early Earth atmosphere

abstract: The emergence of life on the Earth has required a prior organic chemistry leading to the formation of prebiotic molecules. The origin and the evolution of the organic matter on the early Earth is not yet firmly understood. Several hypothesis, possibly complementary, are considered. They can be divided in two categories: endogenous and exogenous sources. In this work we investigate the contribution of a specific endogenous source: the organic chemistry occurring in the ionosphere of the early Earth where the significant VUV contribution of the young Sun involved an efficient formation of reactive species. We address the issue whether this chemistry can lead to the formation of complex organic compounds with CO2 as only source of carbon in an early atmosphere made of N2, CO2 and H2, by mimicking experimentally this type of chemistry using a low pressure plasma reactor. By analyzing the gaseous phase composition, we strictly identified the formation of H2O, NH3, N2O and C2N2. The formation of a solid organic phase is also observed, confirming the possibility to trigger organic chemistry in the upper atmosphere of the early Earth. The identification of Nitrogen-bearing chemical functions in the solid highlights the possibility for an efficient ionospheric chemistry to provide prebiotic material on the early Earth.


title: Highly Ordered and Pinched Magnetic Fields in the Class 0 Proto-Binary System L1448 IRS 2

abstract: We have carried out polarimetric observations with the Atacama Large Millimeter/submillimeter Array (ALMA) toward the Class 0 protostellar system L1448 IRS 2, which is a proto-binary embedded within a flattened, rotating structure, and for which a hint of a central disk has been suggested, but whose magnetic fields are aligned with the bipolar outflow on the cloud core scale. Our high sensitivity and high resolution ($\sim 100$ au) observations show a clear hourglass magnetic field morphology centered on the protostellar system, but the central pattern is consistent with a toroidal field indicative of a circumstellar disk, although other interpretations are also possible, including field lines dragged by an equatorial accretion flow into a configuration parallel to the midplane. If a relatively large disk does exist, it would suggest that the magnetic braking catastrophe is averted in this system, not through a large misalignment between the magnetic and rotation axes, but rather through some other mechanisms, such as non-ideal magneto-hydrodynamic effects and/or turbulence. We have also found a relationship of decreasing polarization fractions with intensities and the various slopes of this relationship can be understood as multiple polarization mechanisms and/or depolarization from a changing field morphology. In addition, we found a prominent clumpy depolarization strip crossing the center perpendicular to the bipolar outflow. Moreover, a rough estimate of the magnetic field strength indicates that the field is strong enough to hinder formation of a rotationally supported disk, which is inconsistent with the feature of a central toroidal field.


title: GumBolt: Extending Gumbel trick to Boltzmann priors

abstract: Boltzmann machines (BMs) are appealing candidates for powerful priors in variational autoencoders (VAEs), as they are capable of capturing nontrivial and multi-modal distributions over discrete variables. However, non-differentiability of the discrete units prohibits using the reparameterization trick, essential for low-noise back propagation. The Gumbel trick resolves this problem in a consistent way by relaxing the variables and distributions, but it is incompatible with BM priors. Here, we propose the GumBolt, a model that extends the Gumbel trick to BM priors in VAEs. GumBolt is significantly simpler than the recently proposed methods with BM prior and outperforms them by a considerable margin. It achieves state-of-the-art performance on permutation invariant MNIST and OMNIGLOT datasets in the scope of models with only discrete latent variables. Moreover, the performance can be further improved by allowing multi-sampled (importance-weighted) estimation of log-likelihood in training, which was not possible with previous models.


title: A Gaia DR 2 and VLT/FLAMES search for new satellites of the LMC

abstract: A wealth of tiny galaxies populates the surroundings of the Milky Way. Some of these objects might have their origin as former satellites of the Large Magellanic Cloud. Examples of the importance of understanding how many systems are genuine satellites of the MW or the LMC are the implications that the number and mass function of satellites have for dark matter theories and the treatment of baryonic physics in simulations of structure formation. Here we aim at deriving the bulk motions and estimates of the internal velocity dispersion and metallicity properties in four recently discovered distant southern dwarf galaxy candidates, Columba I, Reticulum III, Phoenix II and Horologium II. We combine Gaia DR2 astrometric measurements, photometry and new FLAMES spectroscopic data in the region of the CaII triplet lines; such combination is essential for finding potential member stars in these low luminosity systems. We find very likely member stars in all four satellites and are able to determine (or place limits on) the systems bulk motions and average internal properties. The systems are found to be very metal-poor, in agreement with dwarf galaxies and dwarf galaxy candidates of similar luminosity. The only system that we can place firmly in the category of dwarf galaxies is Phx II, given its resolved large velocity dispersion ($9.5_{-4.4}^{+6.8}$ km/s) and intrinsic metallicity spread (0.33 dex). Also for Col I we measure a clear metallicity spread. The orbital pole of Phx II is well constrained and close to that of the LMC, suggesting a prior association. The uncertainty on the orbital poles of the other systems are presently very large, so that an association cannot be excluded, apart from Col I. Using the numbers of potential former satellites of the LMC identified here and in the literature, we obtain for the LMC a dark matter mass of M$_{200}=1.9_{-0.9}^{+1.3}\times10^{11}$ M$_\odot$.


title: Resilient entanglement gates for trapped ions

abstract: Constructing a large scale ion trap quantum processor will require entangling gate operations that are robust in the presence of noise and experimental imperfection. We experimentally demonstrate how a new type of M{\o}lmer-S{\o}rensen gate protects against infidelity caused by heating of the motional mode used during the gate. Furthermore, we show how the same technique simultaneously provides significant protection against slow fluctuations and mis-sets in the secular frequency. Since this parameter sensitivity is worsened in cases where the ions are not ground state cooled, our method provides a path towards relaxing ion cooling requirements in practical realisations of quantum computing and simulation.


title: Quantitative Spectroscopy of Supergiants in the Local Group Dwarf Galaxy IC 1613: Metallicity and Distance

abstract: We present the spectral analysis of 21 blue supergiant stars of spectral type late B to early A within the Local Group dwarf galaxy IC 1613 based on VLT-FORS2 low resolution spectra. Combining our results with studies of early B type blue supergiants we report a wide bi-modal distribution of metallicities with two peaks around [Z] $\sim -0.50$ dex and [Z] $\sim -0.85$ dex. The bi-modal distribution correlates with spatial location, when compared with column densities of neutral hydrogen in IC 1613. While the low [Z] objects appear in regions of relatively high ISM HI column densities or close to them, the high [Z] supergiants are found in the central HI hole almost devoid of hydrogen. This suggests varied chemical evolution histories of the young stellar populations in IC 1613. Utilizing the Flux-Weighted Gravity-Luminosity Relation (FGLR), we determine IC 1613's distance modulus as $m - M$ = 24.39 $\pm$ 0.11 mag. This value is in agreement within previous distance measurements using the near-infrared period-luminosity relationship of Cepheids and the tip of the red giant branch.


title: Model-Driven Engineering of Self-Adaptive Software with EUREMA

abstract: The development of self-adaptive software requires the engineering of an adaptation engine that controls the underlying adaptable software by feedback loops. The engine often describes the adaptation by runtime models representing the adaptable software and by activities such as analysis and planning that use these models. To systematically address the interplay between runtime models and adaptation activities, runtime megamodels have been proposed. A runtime megamodel is a specific model capturing runtime models and adaptation activities. In this article, we go one step further and present an executable modeling language for ExecUtable RuntimE MegAmodels (EUREMA) that eases the development of adaptation engines by following a model-driven engineering approach. We provide a domain-specific modeling language and a runtime interpreter for adaptation engines, in particular feedback loops. Megamodels are kept alive at runtime and by interpreting them, they are directly executed to run feedback loops. Additionally, they can be dynamically adjusted to adapt feedback loops. Thus, EUREMA supports development by making feedback loops explicit at a higher level of abstraction and it enables solutions where multiple feedback loops interact or operate on top of each other and self-adaptation co-exists with offline adaptation for evolution.


title: A Testing Scheme for Self-Adaptive Software Systems with Architectural Runtime Models

abstract: Self-adaptive software systems (SASS) are equipped with feedback loops to adapt autonomously to changes of the software or environment. In established fields, such as embedded software, sophisticated approaches have been developed to systematically study feedback loops early during the development. In order to cover the particularities of feedback, techniques like one-way and in-the-loop simulation and testing have been included. However, a related approach to systematically test SASS is currently lacking. In this paper we therefore propose a systematic testing scheme for SASS that allows engineers to test the feedback loops early in the development by exploiting architectural runtime models. These models that are available early in the development are commonly used by the activities of a feedback loop at runtime and they provide a suitable high-level abstraction to describe test inputs as well as expected test results. We further outline our ideas with some initial evaluation results by means of a small case study.


title: The Sub-Geometric Phases in Density Matrix

abstract: In this letter, the generalization of geometric phase in density matrix is presented, we show that the extended sub-geometric phase have unified expression whatever in adiabatic or nonadiabatic procedure, the relations between them and the usual Berry phase or Aharonov-Anandan phase are established. We also demonstrated the influence of sub-geometric phases on the physical observables. Finally, our treatment is naturally used to investigate the geometric phase in mixed state.


title: ALVEC: Auto-scaling by Lotka Volterra Elastic Cloud: A QoS aware Non Linear Dynamical Allocation Model

abstract: Elasticity in resource allocation is still a relevant problem in cloud computing. There are many academic and white papers which have investigated the problem and offered solutions.\textbf{Unfortunately, there are scant evidence of determining scaling number dynamically.} Elasticity is defined as the ability to adapt with the changing workloads by provisioning and de-provisioning Cloud resources. We propose ALVEC, a novel model of resource allocation in Cloud data centers, inspired by population dynamics and Mathematical Biology, which addresses dynamic allocation by auto-tuning model parameters. The proposed model, governed by a coupled differential equation known as Lotka Volterra (LV), fares better in Service level agreement (SLA) management and Quality of Services (QoS). We show evidence of true elasticity, in theory and empirical comparisons. Additionally, ALVEC is able to predict the future load and allocate VM's accordingly. The proposed model, ALVEC is the first example of unsupervised resource allocation scheme.


title: A Study on Phase-Field Models for Brittle Fracture

abstract: In the phase-field modeling of brittle fracture, anisotropic constitutive assumptions for the degradation of stored elastic energy due to fracture are crucial to preventing cracking in compression and obtaining physically sound numerical solutions. Three energy decomposition models, the spectral decomposition, the volumetric-deviatoric split, and an improved volumetric-deviatoric split, and their effects on the performance of the phase-field modeling are studied. Meanwhile, anisotropic degradation of stiffness may lead to a small energy remaining on crack surfaces, which violates crack boundary conditions and can cause unphysical crack openings and propagation. A simple yet effective treatment for this is proposed: define a critically damaged zone with a threshold parameter and then degrade both the active and passive energies in the zone. A dynamic mesh adaptation finite element method is employed for the numerical solution of the corresponding elasticity system. Four examples, including two benchmark ones, one with complex crack systems, and one based on an experimental setting, are considered. Numerical results show that the spectral decomposition and improved volumetric-deviatoric split models, together with the improvement treatment of crack boundary conditions, can lead to crack propagation results that are comparable with the existing computational and experimental results. It is also shown that the numerical results are not very sensitive to the parameter defining the critically damaged zone.


title: Generators of invariant linear system on tropical curves for finite isometry group

abstract: For a tropical curve and a finite subgroup of the isometry group of the tropical curve, we prove, extending the work by Haase, Musiker and Yu, that the invariant part of the complete linear system associated to a invariant effective divisor on the tropical curve is finitely generated.


title: Gaia GraL II - Gaia DR2 Gravitational Lens Systems: The Known Multiply Imaged Quasars

abstract: Thanks to its spatial resolution the ESA/Gaia space mission offers a unique opportunity to discover new multiply-imaged quasars and to study the already known lensed systems at sub-milliarcsecond astrometric precisions. In this paper, we address the detection of the known multiply-imaged quasars from the Gaia Data Release 2 and determine the astrometric and photometric properties of the individually detected images found in the Gaia DR2 catalogue. We have compiled an exhaustive list of quasar gravitational lenses from the literature to search for counterparts in the Gaia Data Release 2. We then analyze the astrometric and photometric properties of these Gaia's detections. To highlight the tremendous potential of Gaia at the sub-milliarcsecond level we finally perform a simple Bayesian modeling of the well-known gravitational lens system HE0435-1223, using Gaia Data Release 2 and HST astrometry. From 478 known multiply imaged quasars, 200 have at least one image found in the Gaia Data Release 2. Among the 41 known quadruply-imaged quasars of the list, 26 have at least one image in the Gaia Data Release 2, 12 of which are fully detected, 6 have three counterparts, 7 have two and 1 has only one. As expected, the modeling of HE0435-1223 shows that the model parameters are significantly better constrained when using Gaia astrometry compared to HST astrometry, in particular the relative positions of the background quasar source and the centroid of the deflector. The Gaia sub-milliarcsecond astrometry also significantly reduces the parameter correlations. Besides providing an up-to-date list of multiply imaged quasars and their detection in the Gaia DR2, this paper shows that more complex modeling scenarios will certainly benefit from Gaia sub-milliarcsecond astrometry.


title: Prediction in Projection: A new paradigm in delay-coordinate reconstruction

abstract: Delay-coordinate embedding is a powerful, time-tested mathematical framework for reconstructing the dynamics of a system from a series of scalar observations. Most of the associated theory and heuristics are overly stringent for real-world data, however, and real-time use is out of the question due to the expert human intuition needed to use these heuristics correctly. The approach outlined in this thesis represents a paradigm shift away from that traditional approach. I argue that perfect reconstructions are not only unnecessary for the purposes of delay-coordinate based forecasting, but that they can often be less effective than reduced-order versions of those same models. I demonstrate this using a range of low- and high-dimensional dynamical systems, showing that forecast models that employ imperfect reconstructions of the dynamics---i.e., models that are not necessarily true embeddings---can produce surprisingly accurate predictions of the future state of these systems. I develop a theoretical framework for understanding why this is so. This framework, which combines information theory and computational topology, also allows one to quantify the amount of predictive structure in a given time series, and even to choose which forecast method will be the most effective for those data.


title: Dark matter and long-lived particles at the LHC

abstract: While the paradigm of a weakly interacting massive particle (WIMP) has guided our search strategies for dark matter in the past decades, their null-results have stimulated growing interest in alternative explanations pointing towards non-standard signatures. In this article we discuss the phenomenology of dark matter models that predict long-lived particle at the LHC. We focus on models with a $Z_2$-odd dark sector where - in decreasing order of the dark matter coupling - a coannihilation, conversion-driven freeze-out or superWIMP/freeze-in scenario could be realized.


title: The Stochastic Axion Scenario

abstract: For the minimal QCD axion model it is generally believed that overproduction of dark matter constrains the axion mass to be above a certain threshold, or at least that the initial misalignment angle must be tuned if the mass is below that threshold. We demonstrate that this is incorrect. During inflation, if the Hubble scale is low, the axion tends toward an equilibrium. This means the minimal QCD axion can naturally give the observed dark matter abundance in the entire lower part of the mass range, down to masses $\sim 10^{-12}$ eV (or $f_a$ up to almost the Planck scale). The axion abundance is generated by quantum fluctuations of the field during inflation. This mechanism generates cold dark matter with negligible isocurvature perturbations. In addition to the QCD axion, this mechanism can also generate a cosmological abundance of axion-like particles and other light fields.


title: RR Lyrae Variable Stars in the Crater II Dwarf Galaxy

abstract: We report the detection of RR Lyrae variable stars in Crater II, a recently discovered large and diffuse satellite dwarf galaxy of the Milky Way (MW). Based on B, V time-series photometry obtained with the Korea Microlensing Telescope Network (KMTNet) 1.6-m telescope at CTIO, we identified 83 ab-type and 13 c-type pulsators by fitting template light curves. The detected RR Lyrae stars are centrally concentrated, which ensures that most of them are members of Crater II. In terms of the distribution of RRab stars in the period-amplitude diagram, Crater II is clearly different from ultra-faint dwarf (UFD) galaxies, but very similar to the two classical MW dwarf spheroidal (dSph) galaxies Draco and Carina with Oosterhoff-intermediate (Oo-int) properties. Combined with the mean period of ab-type variables (<Pab> = 0.631+-0.004 d) and the c-type fraction (~0.14) in Crater II, this suggests an Oo-int classification for Crater II and implies that its nature is more like a dSph rather than a UFD. We also estimated the mean metallicity, reddening, and distance of Crater II, from the photometric and pulsation properties of the RR Lyrae stars. The stellar population model we have constructed indicates that Crater II is dominated by an old population, but is relatively younger than the oldest globular clusters in the MW. With a lack of high-amplitude short-period RRab stars, Crater II, like most of the other less massive dSphs, is probably not a surviving counterpart of the major building blocks of the MW halo.


title: Characterizing Circumgalactic Gas around Massive Ellipticals at z ~ 0.4 I. Initial Results

abstract: We present a new Hubble Space Telescope (HST) Cosmic Origins Spectrograph (COS) absorption-line survey to study halo gas around 16 luminous red galaxies (LRGs) at z=0.21-0.55. The LRGs are selected uniformly with stellar mass Mstar>1e11 Msun and no prior knowledge of the presence/absence of any absorption features. Based on observations of the full Lyman series, we obtain accurate measurements of neutral hydrogen column density N(HI) and find that high-N(HI) gas is common in these massive quiescent halos with a median of <log N(HI)> = 16.6 at projected distances d<~160 kpc. We measure a mean covering fraction of optically-thick gas with log N(HI)>~17.2 of <kappa>LLS=0.44^{+0.12}_{-0.11} at d<~160 kpc and <kappa>LLS=0.71^{+0.11}_{-0.20} at d<~100 kpc. The line-of-sight velocity separations between the HI absorbing gas and LRGs are characterized by a mean and dispersion of <v_{gas-gal}>=29 km/s and \sigma_v_{gas-gal}=171 km/s. Combining COS FUV and ground-based echelle spectra provides an expanded spectral coverage for multiple ionic transitions, from low-ionization MgII and SiII, to intermediate ionization SiIII and CIII, and to high-ionization OVI absorption lines. We find that intermediate ions probed by CIII and SiIII are the most prominent UV metal lines in LRG halos with a mean covering fraction of <kappa(CIII)>_{0.1}=0.75^{+0.08}_{-0.13} for W(977)>=0.1 Ang at d<160 kpc, comparable to what is seen for CIII in L* and sub-L* star-forming and red galaxies but exceeding MgII or OVI in quiescent halos. The COS-LRG survey shows that massive quiescent halos contain widespread chemically-enriched cool gas and that little distinction between LRG and star-forming halos is found in their HI and CIII content.


title: A descriptive set theorist's proof of the pointwise ergodic theorem

abstract: We give a short combinatorial proof of the classical pointwise ergodic theorem for probability measure preserving $\mathbb{Z}$-actions. Our approach reduces the theorem to a tiling problem: tightly tile each orbit by intervals with desired averages. This tiling problem is easy to solve for $\mathbb{Z}$ with intervals as tiles. However, it would be interesting to find other classes of groups and sequences of tiles for which this can be done, since then our approach would yield a pointwise ergodic theorem for such classes.


title: Variable broad lines and outflow in the weak blazar PBC J2333.9-2343

abstract: PBC J2333.9-2343 is a peculiar active nucleus with two giant radio lobes and a weak blazar-like nucleus at their center. In the present work we show new optical, UV, and X-ray data taken from the San Pedro M\'artir telescope, the New Technology Telescope, NTT/EFOSC2, and the Swift/XRT satellite. The source is highly variable at all frequencies, in particular the strongest variations are found in the broad H$\alpha$ component with a flux increase of 61$\pm$4 per cent between 2009 and 2016, following the X-ray flux increase of 62$\pm$6 per cent between 2010 and 2016. We also detected a broad H$\beta$ component in 2016, making the optical classification change from type 1.9 to type 1.8 in one year. We have also detected a broad component of the [OIII]$\lambda$5007 line, which is blue-shifted and of high velocity, suggesting an origin from a highly disturbed medium, possibly an outflow. The line flux variability and broad widths are indicative of a jet that is, at least in part, responsible for the ionization of the BLR and NLR.


title: Energetic $\gamma$-rays from TeV scale dark matter annihilation resummed

abstract: The annihilation cross section of TeV scale dark matter particles $\chi^0$ with electroweak charges into photons is affected by large quantum corrections due to Sudakov logarithms and the Sommerfeld effect. We calculate the semi-inclusive photon energy spectrum in $\chi^0\chi^0\to \gamma+X$ in the vicinity of the maximal photon energy $E_\gamma = m_\chi$ with NLL' accuracy in an all-order summation of the electroweak perturbative expansion adopting the pure wino model. This results in the most precise theoretical prediction of the annihilation rate for $\gamma$-ray telescopes with photon energy resolution of parametric order $m_W^2/m_\chi$ for photons with TeV energies.


title: Do Diffusion Protocols Govern Cascade Growth?

abstract: Large cascades can develop in online social networks as people share information with one another. Though simple reshare cascades have been studied extensively, the full range of cascading behaviors on social media is much more diverse. Here we study how diffusion protocols, or the social exchanges that enable information transmission, affect cascade growth, analogous to the way communication protocols define how information is transmitted from one point to another. Studying 98 of the largest information cascades on Facebook, we find a wide range of diffusion protocols - from cascading reshares of images, which use a simple protocol of tapping a single button for propagation, to the ALS Ice Bucket Challenge, whose diffusion protocol involved individuals creating and posting a video, and then nominating specific others to do the same. We find recurring classes of diffusion protocols, and identify two key counterbalancing factors in the construction of these protocols, with implications for a cascade's growth: the effort required to participate in the cascade, and the social cost of staying on the sidelines. Protocols requiring greater individual effort slow down a cascade's propagation, while those imposing a greater social cost of not participating increase the cascade's adoption likelihood. The predictability of transmission also varies with protocol. But regardless of mechanism, the cascades in our analysis all have a similar reproduction number ($\approx$ 1.8), meaning that lower rates of exposure can be offset with higher per-exposure rates of adoption. Last, we show how a cascade's structure can not only differentiate these protocols, but also be modeled through branching processes. Together, these findings provide a framework for understanding how a wide variety of information cascades can achieve substantial adoption across a network.


title: Vortices in a generalized Maxwell-Higgs model with visible and hidden sectors

abstract: We investigate the presence of vortices in generalized Maxwell-Higgs models with a hidden sector. The model engenders $U(1)\times U(1)$ symmetry, in a manner that the sectors are coupled via the visible magnetic permeability depending only on the hidden scalar field. We develop a first order framework in which the hidden sector decouples from the visible one. We illustrate the results with two specific examples, that give rise to the presence of vortices with internal structure.


title: Deep and rapid observations of strong-lensing galaxy clusters within the sky localisation of GW170814

abstract: We present observations of two strong-lensing galaxy clusters located within the $90$ per cent credible sky localization maps released following LIGO-Virgo's discovery of the binary black hole (BH-BH) gravitational wave (GW) source GW170814. Our objectives were (1) to search for candidate electromagnetic (EM) counterparts to GW170814 under the hypothesis that it was strongly-lensed, and thus more distant and less massive than inferred by LIGO-Virgo, and (2) to demonstrate the feasibility of rapid target of opportunity observations to search for faint lensed transient point sources in crowded cluster cores located within GW sky localizations. Commencing $20$ hours after discovery, and continuing over $12$ nights, we observed Abell 3084 ($z=0.22$) and SMACSJ0304.3$-$4401 ($z=0.46$) with GMOS on the Gemini-South telescope, and Abell 3084 with MUSE on ESO's Very Large Telescope. We detect no candidate EM counterparts in these data. Calibration of our photometric analysis methods using simulations yield $5\sigma$ detection limits for transients in difference images of the cores of these clusters of $i=25$. This is the most sensitive photometric search to date for counterparts to GW sources, and rules out the possibility that GW170814 was lensed by these clusters with a kilonova-like EM counterpart. Based on the detector frame masses of the compact objects, and assuming that at least one Neutron Star (NS) is required in the merging system to produce a kilonova-like counterpart, implies that GW170814 was neither a NS-NS nor NS-BH merger at $z>8$ lensed by either of these clusters. Also, in the first ever emission line search for counterparts to GW sources, we detected no lines down to a $5\sigma$ detection limit of $5\times10^{-17}{\rm erg\,s^{-1}\,cm^{-2}}$.