title: Far-field theory for trajectories of magnetic ellipsoids in rectangular and circular channels

abstract: We report a method to control the positions of ellipsoidal magnets in flowing channels of rectangular or circular cross section at low Reynolds number.A static uniform magnetic field is used to pin the particle orientation, and the particles move with translational drift velocities resulting from hydrodynamic interactions with the channel walls which can be described using Blake's image tensor.Building on his insights, we are able to present a far-field theory predicting the particle motion in rectangular channels, and validate the accuracy of the theory by comparing to numerical solutions using the boundary element method.We find that, by changing the direction of the applied magnetic field, the motion can be controlled so that particles move either to a curved focusing region or to the channel walls.We also use simulations to show that the particles are focused to a single line in a circular channel.Our results suggest ways to focus and segregate magnetic particles in lab-on-a-chip devices.


title: Inclusive production of vector quarkonia at the LHC

abstract: We discuss prompt production of $J/\psi$ mesons in proton-proton collisions at the LHC within NRQCD k$t$-factorization approach using Kimber-Martin-Ryskin (KMR) and Kutak-Stasto (KS) unintegrated gluon distribution (UGDF). We include both direct color-singlet production $gg \to J/\psi g$ as well as a feed-down from $\chi_{c} \to J/\psi \gamma$ and $\psi' \to J/\psi X$. The corresponding matrix elements for $gg \to J/\psi$, $gg \to \psi'$ and $gg \to \chi_{c}$ include parameters of the nonrelativistic space wave functions of the quarkonia at $r=0$, which are taken from potential models from the literature. We calculate the ratio of the corresponding cross sections for $\chi_{c2}/\chi_{c1}$. We compare our results with ATLAS experimental data. Differential distributions in rapidity and transverse momentum of $J/\psi$ and $\psi'$ are calculated and compared to experimental data of the ALICE and LHCb collaborations. We present results for three different values of energy 2.76 TeV, 7 TeV and 13 TeV.


title: Structural Transformations of Carbon and Boron Nitride Nanoscrolls at High Impact Collisions

abstract: The behavior of nanostructures under high strain-rate conditions has been object of theoretical and experimental investigations in recent years. For instance, it has been shown that carbon and boron nitride nanotubes can be unzipped into nanoribbons at high velocity impacts. However, the response of many nanostructures to high strain-rate conditions is still not completely understood. In this work we have investigated through fully atomistic reactive (ReaxFF) molecular dynamics (MD) simulations the mechanical behavior of carbon (CNS) and boron nitride nanoscrolls (BNS) colliding against solid targets at high velocities,. CNS (BNS) nanoscrolls are graphene (boron nitride) membranes rolled up into papyrus-like structures. Their open-ended topology leads to unique properties not found in close-ended analogues, such as nanotubes. Our results show that the collision products are mainly determined by impact velocities and by two impact angles, which define the position of the scroll (i) axis and (ii) open edge relative to the target. Our MD results showed that for appropriate velocities and orientations large-scale deformations and nanoscroll fracture can occur. We also observed unscrolling (scrolls going back to quasi-planar membranes), scroll unzipping into nanoribbons, and significant reconstruction due to breaking and/or formation of new chemical bonds. For particular edge orientations and velocities, conversion from open to close-ended topology is also possible, due to the fusion of nanoscroll walls.


title: Non-Linear Digital Self-Interference Cancellation for In-Band Full-Duplex Radios Using Neural Networks

abstract: Full-duplex systems require very strong self-interference cancellation in order to operate correctly and a significant part of the self-interference signal is due to non-linear effects created by various transceiver impairments. As such, linear cancellation alone is usually not sufficient and sophisticated non-linear cancellation algorithms have been proposed in the literature. In this work, we investigate the use of a neural network as an alternative to the traditional non-linear cancellation method that is based on polynomial basis functions. Measurement results from a full-duplex testbed demonstrate that a small and simple feed-forward neural network canceler works exceptionally well, as it can match the performance of the polynomial non-linear canceler with significantly lower computational complexity.


title: Realization of a Contactless Acoustic Levitation Motor via Doublet Mode Control and Autoresonance

abstract: This paper demonstrates analytically and experimentally an acoustic levitation motor which has the ability to levitate and rotate an object in the air without mechanical contact. To realize such a device two core methods are applied simultaneously; (i) resonance tracking with an Autoresonance feedback loop, (ii) generation of controlled structural traveling waves. The purpose of the first method is to achieve near-field acoustic levitation, which can levitate an object of a few kilograms. In this research, this is accomplished through high amplitude vibration of an aluminum annulus at ultrasonic frequencies (~30kHz). For high efficiency, the annulus is designed to have a very high Q value, and operating even slightly off resonance ceases levitation. Compounding this is the fact that the natural frequency constantly drifts as ambient conditions and loading change. To accommodate such a drift, and produce stable levitation automatically, a resonance tracking feedback loop is employed here. Simultaneously, the purpose of the second method is to achieve propulsion forces on the levitated object by propagating and controlling traveling waves in the aluminum annulus to create a thin layer of rotating air beneath the levitated body. Even though a single vibration mode can produce only standing waves, an axisymmetric structure possesses two modes per natural frequency, and excitation of a doublet mode pair can generate effective travelling waves. The present paper develops the theory behind the use of the Autoresonance feedback method for achieving constant levitation and propagating travelling waves in co- and counter rotating directions. It will be shown that all this can be accomplished with only single sensor. The result is a stable, repeatable and a highly controllable contactless acoustic levitation motor.


title: Magnetic Fields in Star-Forming Filaments in Different Environments

abstract: Cold, dense filaments, some appearing as infrared dark clouds, are the nurseries of stars. Tremendous progress in terms of temperature, density distribution and gas kinematics has been made in understanding the nature of these filaments. However, very little is known about the role played by magnetic fields in the evolution of these filaments. Here, I summarize the recent observational efforts and ongoing projects (POLSTAR survey) in this direction.


title: A Large Dimensional Study of Regularized Discriminant Analysis Classifiers

abstract: This article carries out a large dimensional analysis of standard regularized discriminant analysis classifiers designed on the assumption that data arise from a Gaussian mixture model with different means and covariances. The analysis relies on fundamental results from random matrix theory (RMT) when both the number of features and the cardinality of the training data within each class grow large at the same pace. Under mild assumptions, we show that the asymptotic classification error approaches a deterministic quantity that depends only on the means and covariances associated with each class as well as the problem dimensions. Such a result permits a better understanding of the performance of regularized discriminant analsysis, in practical large but finite dimensions, and can be used to determine and pre-estimate the optimal regularization parameter that minimizes the misclassification error probability. Despite being theoretically valid only for Gaussian data, our findings are shown to yield a high accuracy in predicting the performances achieved with real data sets drawn from the popular USPS data base, thereby making an interesting connection between theory and practice.


title: Heavy quarkonia in a contact interaction and an algebraic model: mass spectrum, decay constants, charge radii and elastic and transition form factors

abstract: For the flavor-singlet heavy quark system of bottomonia, we compute the masses of the ground state mesons in four different channels, namely, pseudo-scalar ($\eta_{b}(1S)$), vector ($\Upsilon(1S)$), scalar ($\chi_{b_0}(1P)$) and axial vector ($\chi_{b_{1}}(1P)$). We also calculate the weak decay constants of the $\eta_{b}(1S)$ and $\Upsilon(1S)$ as well as the charge radius of $\eta_{b}(1S)$. It complements our previous study of the corresponding charmonia systems: $\eta_c(1S)$, $J/\Psi(1S)$, $\chi_{c_0}(1P)$) and ($\chi_{c_{1}}(1P)$). The unified formalism for this analysis is provided by a symmetry-preserving Schwinger-Dyson equations treatment of a vector$\times$vector contact interaction. Whenever a comparison is possible, our results are in fairly good agreement with experimental data and model calculations based upon Schwinger-Dyson and Bethe-Salpeter equations involving sophisticated interaction kernels. Within the same framework, we also report the elastic and transition form factors to two photons for the pseudo-scalar channels $\eta_{c}(1S)$ and $\eta_{b}(1S)$ in addition to the elastic form factors for the vector mesons $J/\Psi$ and $\Upsilon$ for a wide range of photon momentum transfer squared ($Q^2$). For $\eta_{c}(1S)$ and $\eta_{b}(1S)$, we also provide predictions of an algebraic model which correlates remarkably well between the known infrared and ultraviolet limits of these form factors.


title: Protostellar Jets and Outflows in low-mass star formation

abstract: The driving mechanism of protostellar outflows and jets and their effects on the star formation process obtained from recent theoretical and numerical studies are described. Low-velocity outflows are driven by an outer region of the circumstellar disk, while high-velocity jets are driven near an inner edge of the disk. The disk angular momentum is effectively transferred by magnetic effects in the outflow and jet driving regions where the magnetic field is well coupled with neutral gas. On the other hand, in a high density gas region of the disk (or intermediate region), the magnetic field dissipates and is decoupled from neutral gas. Thus, in such a magnetically inactive region, no outward flow appears and the disk angular momentum is not effectively transferred by magnetic effects. Therefore, in the disk intermediate region, the disk surface density continues to increase and gravitational instability occurs and produce a non-axisymmetric (or spiral) structure. After spiral arms sufficiently develop, the disk angular momentum is transferred by gravitational torque and a large amount of the disk mass accretes onto the central protostar from the circumstellar disk. The episodic accretion induces time-variable high-velocity jets. The jets do not significantly contribute to a dynamical evolution of the protostar and circumstellar disk, while the low-velocity outflows can eject a large fraction of the infalling gas and determine the final stellar mass.


title: Pseudorandom States, Non-Cloning Theorems and Quantum Money

abstract: We propose the concept of pseudorandom states and study their constructions, properties, and applications. Under the assumption that quantum-secure one-way functions exist, we present concrete and efficient constructions of pseudorandom states. The non-cloning theorem plays a central role in our study---it motivates the proper definition and characterizes one of the important properties of pseudorandom quantum states. Namely, there is no efficient quantum algorithm that can create more copies of the state from a given number of pseudorandom states. As the main application, we prove that any family of pseudorandom states naturally gives rise to a private-key quantum money scheme.


title: Analyzing the Approximation Error of the Fast Graph Fourier Transform

abstract: The graph Fourier transform (GFT) is in general dense and requires O(n^2) time to compute and O(n^2) memory space to store. In this paper, we pursue our previous work on the approximate fast graph Fourier transform (FGFT). The FGFT is computed via a truncated Jacobi algorithm, and is defined as the product of J Givens rotations (very sparse orthogonal matrices). The truncation parameter, J, represents a trade-off between precision of the transform and time of computation (and storage space). We explore further this trade-off and study, on different types of graphs, how is the approximation error distributed along the spectrum.


title: Airy wave packets accelerating in space-time

abstract: Although diffractive spreading is an unavoidable feature of all wave phenomena, certain waveforms can attain propagation-invariance. A lesser-explored strategy for achieving optical selfsimilar propagation exploits the modification of the spatio-temporal field structure when observed in reference frames moving at relativistic speeds. For such an observer, it is predicted that the associated Lorentz boost can bring to a halt the axial dynamics of a wave packet of arbitrary profile. This phenomenon is particularly striking in the case of a self-accelerating beam -- such as an Airy beam -- whose peak normally undergoes a transverse displacement upon free-propagation. Here we synthesize an acceleration-free Airy wave packet that travels in a straight line by deforming its spatio-temporal spectrum to reproduce the impact of a Lorentz boost. The roles of the axial spatial coordinate and time are swapped, leading to `time-diffraction' manifested in self-acceleration observed in the propagating Airy wave-packet frame.


title: Active Tolerant Testing

abstract: In this work, we give the first algorithms for tolerant testing of nontrivial classes in the active model: estimating the distance of a target function to a hypothesis class C with respect to some arbitrary distribution D, using only a small number of label queries to a polynomial-sized pool of unlabeled examples drawn from D. Specifically, we show that for the class D of unions of d intervals on the line, we can estimate the error rate of the best hypothesis in the class to an additive error epsilon from only $O(\frac{1}{\epsilon^6}\log \frac{1}{\epsilon})$ label queries to an unlabeled pool of size $O(\frac{d}{\epsilon^2}\log \frac{1}{\epsilon})$. The key point here is the number of labels needed is independent of the VC-dimension of the class. This extends the work of Balcan et al. [2012] who solved the non-tolerant testing problem for this class (distinguishing the zero-error case from the case that the best hypothesis in the class has error greater than epsilon). We also consider the related problem of estimating the performance of a given learning algorithm A in this setting. That is, given a large pool of unlabeled examples drawn from distribution D, can we, from only a few label queries, estimate how well A would perform if the entire dataset were labeled? We focus on k-Nearest Neighbor style algorithms, and also show how our results can be applied to the problem of hyperparameter tuning (selecting the best value of k for the given learning problem).


title: Excursion Processes Associated with Elliptic Combinatorics

abstract: Researching elliptic analogues for equalities and formulas is a new trend in enumerative combinatorics which has followed the previous trend of studying $q$-analogues. Recently Schlosser proposed a lattice path model in the square lattice with a family of totally elliptic weight-functions including several complex parameters and discussed an elliptic extension of the binomial theorem. In the present paper, we introduce a family of discrete-time excursion processes on $\mathbb{Z}$ starting from the origin and returning to the origin in a given time duration $2T$ associated with Schlosser's elliptic combinatorics. The processes are inhomogeneous both in space and time and hence expected to provide new models in non-equilibrium statistical mechanics. By numerical calculation we show that the maximum likelihood trajectories on the spatio-temporal plane of the elliptic excursion processes and of their reduced trigonometric versions are not straight lines in general but are nontrivially curved depending on parameters. We analyze asymptotic probability laws in the long-term limit $T \to \infty$ for a simplified trigonometric version of excursion process. Emergence of nontrivial curves of trajectories in a large scale of space and time from the elementary elliptic weight-functions exhibits a new aspect of elliptic combinatorics.

categories: math-ph cond-mat.stat-mech math.MP math.PR nlin.SI


title: AGT relations for sheaves on surfaces

abstract: We consider a natural generalization of the Carlsson-Okounkov Ext operator on the K-theory groups of the moduli spaces of stable sheaves on a smooth projective surface. We compute the commutation relations between the Ext operator and the action of the deformed W-algebra on K-theory, which was developed in [18]. The conclusion is that the Ext operator is closely related with a vertex operator, thus giving a mathematical incarnation of the Alday-Gaiotto-Tachikawa correspondence for a general algebraic surface.


title: Echoes of Kerr-like wormholes

abstract: Structure at the horizon scale of black holes would give rise to echoes of the gravitational wave signal associated with the post-merger ringdown phase in binary coalescences. We study the waveform of echoes in static and stationary, traversable wormholes in which perturbations are governed by a symmetric effective potential. We argue that echoes are dominated by the wormhole quasinormal frequency nearest to the fundamental black hole frequency that controls the primary signal. We put forward an accurate method to construct the echoes waveform(s) from the primary signal and the quasinormal frequencies of the wormhole, which we characterize. We illustrate this in the static Damour-Solodukhin wormhole and in a new, rotating generalization that approximates a Kerr black hole outside the throat. Rotation gives rise to a potential with an intermediate plateau region that breaks the degeneracy of the quasinormal frequencies. Rotation also leads to late-time instabilities which, however, fade away for small angular momentum.


title: Symmetries in Two Dimensional Conformal Field Theories and Related Integrable Models

abstract: We present part of our investigations on two dimensional N=1 and N=2 superconformal field theories. As a direct generalization we consider the SU(2) coset models, in particular their renormalization group properties. A search and possible implementation of additional symmetries in related integrable models are also presented.


title: Interplanetary Scintillation studies with the Murchison Wide-field Array II: Properties of sub-arcsecond compact sources at low radio frequencies

abstract: We report the first astrophysical application of the technique of wide-field Interplanetary Scintillation (IPS) with the Murchison Widefield Array (MWA). This powerful technique allows us to identify and measure sub-arcsecond compact components in low-frequency radio sources across large areas of sky without the need for long-baseline interferometry or ionospheric calibration. We present the results of a five-minute observation of a 30x30 sq. deg MWA field at 162 MHz with 0.5 second time resolution. Of the 2550 continuum sources detected in this field, 302 (12 per cent) show rapid fluctuations caused by IPS. We find that at least 32% of bright low-frequency radio sources contain a sub-arcsec compact component that contributes over 40% of the total flux density. Perhaps surprisingly, peaked-spectrum radio sources are the dominant population among the strongly-scintillating, low-frequency sources in our sample. While gamma-ray AGN are generally compact, flat-spectrum radio sources at higher frequencies, the 162 MHz properties of many of the Fermi blazars in our field are consistent with a compact component embedded within more extended low-frequency emission. The detection of a known pulsar in our field shows that the wide-field IPS technique is at the threshold of sensitivity needed to detect new pulsars using image plane analysis, and scaling the current MWA sensitivity to that expected for SKA-low implies that large IPS-based pulsar searches will be feasible with SKA. Calibration strategies for the SKA require a better knowledge of the space density of compact sources at low radio frequencies, which IPS observations can now provide.


title: Universal gradient descent

abstract: In this book we collect many different and useful facts around gradient descent method. First of all we consider gradient descent with inexact oracle. We build a general model of optimized function that include composite optimization approach, level's methods, proximal methods etc. Then we investigate primal-dual properties of the gradient descent in general model set-up. At the end we generalize method to universal one.


title: Identifying Coherent Anomalies in Multi-Scale Spatio-Temporal Data using Markov Random Fields

abstract: Many physical processes involve spatio-temporal observations, which can be studied at different spatial and temporal scales. For example, rainfall data measured daily by rain gauges can be considered at daily, monthly or annual temporal scales, and local, grid-wise, region-wise or country-wise spatial scales. In this work, we focus on detection of anomalies in such multi-scale spatio-temporal data. We consider an anomaly as an event where the measured values over a spatio-temporally extended region are significantly different from their long-term means. However we aim to avoid setting any thresholds on the measured values and spatio-temporal sizes, because not only are thresholds subjective but also the long-term mean values often vary spatially and temporally. For this purpose we use spatio-Temporal Markov Random Field, where latent states indicate anomaly type (positive anomaly, negative anomaly, no anomaly/normal). Spatio-temporal coherence is maintained through suitable edge potentials. The model is extended to multiple spatio-temporal scales to achieve our second goal: anomalies at any scale should be defined both on the data at that scale, and also on anomalies at other scales. This allows us to trace an anomaly at a coarse scale to finer scales. For example, whether rainfall in a particular year is anomalous over a region should depend not only on the total volume of rainfall over the entire region, but also on whether there were such anomalies at the grid-scale, and the monthly scale. We use this approach to study rainfall anomalies over India -extremely diverse with respect to rainfall- for the period 1901-2011, and show its benefits over existing approaches.


title: Pearling Bifurcations in the strong Functionalized Cahn-Hilliard Free Energy

abstract: The Functionalized Cahn-Hilliard free energy supports phase separated morphologies of distinct codimension, including codimension-one bilayer and codimension-two filament morphologies. We characterize the linear stability of bilayer and filament morphologies associated to hypersurfaces within the strong functionalization scaling. In particular we show that the onset of the pearling instability, which triggers fast in-plane oscillations associated to bifurcation to higher codimensional morphology, is controlled by the functionalization parameters and the spatially constant value of the far-field chemical potential. Crucially, we show that onset of pearling is independent of the shape of the defining hypersurface.


title: Local noncentrosymmetric structure of Bi2Sr2CaCu2O8+y by X-ray magnetic circular dichroism at Cu K-edge XANES

abstract: The two-dimensional Bi2Sr2CaCu2O8+y (Bi2212), the most studied prototype cuprate superconductor, is a lamellar system made of a stack of two-dimensional corrugated CuO2 bilayers separated by Bi2O2+ySr2O2 layers. While the large majority of theories, proposed to interpret unconventional high Tc superconductivity in Bi2Sr2CaCu2O8+y, assume a centrosymmetric tetragonal CuO2 lattice for the [CuO2]Ca[CuO2] bilayer here we report new compelling results providing evidence for local noncentrosymmetric structure at the Cu site. We have measured polarized Cu K-edge XANES (x-ray absorption near edge structure) and the K-edge X-ray magnetic circular dichroism (XMCD) of a Bi2212 single crystal near optimum doping. The Cu K edge XMCD signal was measured at ID12 beamline of ESRF with the k-vector of x-ray beam parallel to c-axis i.e. with the electric field of x-ray beam E//ab, using a 17 T magnetic field parallel to the c-axis of a Bi2212 single crystal. Numerical simulations of the XMCD signal of Bi2212 by multiple scattering theory have shown agreement with the experimental XMCD signal only for the local structure with noncentrosymmetric Bb2b space group of Bi2Sr2CaCu2O8+y.


title: Time-Triggered Co-Scheduling of Computation and Communication with Jitter Requirements

abstract: The complexity of embedded application design is increasing with growing user demands. In particular, automotive embedded systems are highly complex in nature, and their functionality is realized by a set of periodic tasks. These tasks may have hard real-time requirements and communicate over an interconnect. The problem is to efficiently co-schedule task execution on cores and message transmission on the interconnect so that timing constraints are satisfied. Contemporary works typically deal with zero-jitter scheduling, which results in lower resource utilization, but has lower memory requirements. This article focuses on jitter-constrained scheduling that puts constraints on the tasks jitter, increasing schedulability over zero- jitter scheduling. The contributions of this article are: 1) Integer Linear Programming and Satisfiability Modulo Theory model exploiting problem-specific information to reduce the formulations complexity to schedule small applications. 2) A heuristic approach, employing three levels of scheduling scaling to real-world use-cases with 10000 tasks and messages. 3) An experimental evaluation of the proposed approaches on a case-study and on synthetic data sets showing the efficiency of both zero-jitter and jitter- constrained scheduling. It shows that up to 28% higher resource utilization can be achieved by having up to 10 times longer computation time with relaxed jitter requirements.


title: Counterfactual Explanations without Opening the Black Box: Automated Decisions and the GDPR

abstract: There has been much discussion of the right to explanation in the EU General Data Protection Regulation, and its existence, merits, and disadvantages. Implementing a right to explanation that opens the black box of algorithmic decision-making faces major legal and technical barriers. Explaining the functionality of complex algorithmic decision-making systems and their rationale in specific cases is a technically challenging problem. Some explanations may offer little meaningful information to data subjects, raising questions around their value. Explanations of automated decisions need not hinge on the general public understanding how algorithmic systems function. Even though such interpretability is of great importance and should be pursued, explanations can, in principle, be offered without opening the black box. Looking at explanations as a means to help a data subject act rather than merely understand, one could gauge the scope and content of explanations according to the specific goal or action they are intended to support. From the perspective of individuals affected by automated decision-making, we propose three aims for explanations: (1) to inform and help the individual understand why a particular decision was reached, (2) to provide grounds to contest the decision if the outcome is undesired, and (3) to understand what would need to change in order to receive a desired result in the future, based on the current decision-making model. We assess how each of these goals finds support in the GDPR. We suggest data controllers should offer a particular type of explanation, unconditional counterfactual explanations, to support these three aims. These counterfactual explanations describe the smallest change to the world that can be made to obtain a desirable outcome, or to arrive at the closest possible world, without needing to explain the internal logic of the system.


title: Minimal Exploration in Structured Stochastic Bandits

abstract: This paper introduces and addresses a wide class of stochastic bandit problems where the function mapping the arm to the corresponding reward exhibits some known structural properties. Most existing structures (e.g. linear, Lipschitz, unimodal, combinatorial, dueling, ...) are covered by our framework. We derive an asymptotic instance-specific regret lower bound for these problems, and develop OSSB, an algorithm whose regret matches this fundamental limit. OSSB is not based on the classical principle of "optimism in the face of uncertainty" or on Thompson sampling, and rather aims at matching the minimal exploration rates of sub-optimal arms as characterized in the derivation of the regret lower bound. We illustrate the efficiency of OSSB using numerical experiments in the case of the linear bandit problem and show that OSSB outperforms existing algorithms, including Thompson sampling.


title: Transitions from trees to cycles in adaptive flow networks

abstract: Transport networks are crucial to the functioning of natural and technological systems. Nature features transport networks that are adaptive over a vast range of parameters, thus providing an impressive level of robustness in supply. Theoretical and experimental studies have found that real-world transport networks exhibit both tree-like motifs and cycles. When the network is subject to load fluctuations, the presence of cyclic motifs may help to reduce flow fluctuations and, thus, render supply in the network more robust. While previous studies considered network topology via optimization principles, here, we take a dynamical systems approach and study a simple model of a flow network with dynamically adapting weights (conductances). We assume a spatially non-uniform distribution of rapidly fluctuating loads in the sinks and investigate what network configurations are dynamically stable. The network converges to a spatially non-uniform stable configuration composed of both cyclic and tree-like structures. Cyclic structures emerge locally in a transcritical bifurcation as the amplitude of the load fluctuations is increased. The resulting adaptive dynamics thus partitions the network into two distinct regions with cyclic and tree-like structures. The location of the boundary between these two regions is determined by the amplitude of the fluctuations. These findings may explain why natural transport networks display cyclic structures in the micro-vascular regions near terminal nodes, but tree-like features in the regions with larger veins.


title: Successive Cancellation Soft Output Detector For Uplink MU-MIMO Systems With One-bit ADCs

abstract: In this paper, we present a successive-cancellation-soft-output (SCSO) detector for an uplink multiuser multiple-input-multiple-output (MU-MIMO) system with one-bit analog-to-digital converters (ADCs). The proposed detector produces soft outputs (e.g., log-likelihood ratios (LLRs)) from one-bit quantized observations in a successive way: each user k's message is sequentially decoded from a channel decoder $k$ for $k=1,...,K$ in that order, and the previously decoded messages are exploited to improve the reliabilities of LLRs. Furthermore, we develop an efficient greedy algorithm to optimize a decoding order. Via simulation results, we demonstrate that the proposed ordered SCSO detector outperforms the other detectors for the coded MU-MIMO systems with one-bit ADCs.


title: Dark Energy Survey Year 1 Results: A Precise H0 Measurement from DES Y1, BAO, and D/H Data

abstract: We combine Dark Energy Survey Year 1 clustering and weak lensing data with Baryon Acoustic Oscillations (BAO) and Big Bang Nucleosynthesis (BBN) experiments to constrain the Hubble constant. Assuming a flat $\Lambda$CDM model with minimal neutrino mass ($\sum m_\nu = 0.06$ eV) we find $H_0=67.2^{+1.2}_{-1.0}$ km/s/Mpc (68% CL). This result is completely independent of Hubble constant measurements based on the distance ladder, Cosmic Microwave Background (CMB) anisotropies (both temperature and polarization), and strong lensing constraints. There are now five data sets that: a) have no shared observational systematics; and b) each constrain the Hubble constant with a few percent level precision. We compare these five independent measurements, and find that, as a set, the differences between them are significant at the $2.1\sigma$ level ($\chi^2/dof=20.1/11$, probability to exceed=4%). This difference is low enough that we consider the data sets statistically consistent with each other. The best fit Hubble constant obtained by combining all five data sets is $H_0 = 69.1^{+0.4}_{-0.6}$ km/s/Mpc.


title: Building Data-driven Models with Microstructural Images: Generalization and Interpretability

abstract: As data-driven methods rise in popularity in materials science applications, a key question is how these machine learning models can be used to understand microstructure. Given the importance of process-structure-property relations throughout materials science, it seems logical that models that can leverage microstructural data would be more capable of predicting property information. While there have been some recent attempts to use convolutional neural networks to understand microstructural images, these early studies have focused only on which featurizations yield the highest machine learning model accuracy for a single data set. This paper explores the use of convolutional neural networks for classifying microstructure with a more holistic set of objectives in mind: generalization between data sets, number of features required, and interpretability.


title: The Price of Information in Combinatorial Optimization

abstract: Consider a network design application where we wish to lay down a minimum-cost spanning tree in a given graph; however, we only have stochastic information about the edge costs. To learn the precise cost of any edge, we have to conduct a study that incurs a price. Our goal is to find a spanning tree while minimizing the disutility, which is the sum of the tree cost and the total price that we spend on the studies. In a different application, each edge gives a stochastic reward value. Our goal is to find a spanning tree while maximizing the utility, which is the tree reward minus the prices that we pay. Situations such as the above two often arise in practice where we wish to find a good solution to an optimization problem, but we start with only some partial knowledge about the parameters of the problem. The missing information can be found only after paying a probing price, which we call the price of information. What strategy should we adopt to optimize our expected utility/disutility? A classical example of the above setting is Weitzman's "Pandora's box" problem where we are given probability distributions on values of $n$ independent random variables. The goal is to choose a single variable with a large value, but we can find the actual outcomes only after paying a price. Our work is a generalization of this model to other combinatorial optimization problems such as matching, set cover, facility location, and prize-collecting Steiner tree. We give a technique that reduces such problems to their non-price counterparts, and use it to design exact/approximation algorithms to optimize our utility/disutility. Our techniques extend to situations where there are additional constraints on what parameters can be probed or when we can simultaneously probe a subset of the parameters.


title: Tunable magnetization relaxation of Fe_{2}Cr_{1-x}Co_{x}Si half-metallic Heusler alloys by band structure engineering

abstract: We report a systematic investigation on the magnetization relaxation properties of iron-based half-metallic Heusler alloy Fe$_{2}$Cr$_{1-x}$Co_${x}$Si (FCCS) thin films using broadband angular-resolved ferromagnetic resonance. Band structure engineering through Co doping (x) demonstrated by first-principles calculations is shown to tune the intrinsic magnetic damping over an order of magnitude, namely 0.01-0.0008. Notably, the intrinsic damping constants for samples with high Co concentration are among the lowest reported for Heusler alloys and even comparable to magnetic insulator yttrium iron garnet. Furthermore, a significant reduction of both isotropic and anisotropic contributions of extrinsic damping of the FCCS alloys was found in the FCCS films with x=0.5-0.75, which is of particular importance for applications. These results demonstrate a practical recipe to tailor functional magnetization for Heusler alloy-based spintronics at room temperature


title: ALMA + VLT observations of a Damped Lyman-{\alpha} absorbing galaxy: Massive, wide CO emission, gas-rich but with very low SFR

abstract: We are undertaking an ALMA survey of molecular gas in galaxies selected for their strong HI absorption, so-called DLA/sub-DLA galaxies. Here we report CO(2-1) detection from a DLA galaxy at z = 0.716. We also present optical and near-infrared spectra of the galaxy revealing [OII], H{\alpha} and [NII] emission lines shifted by ~170 km/s relative to the DLA, and providing an oxygen abundance 3.2 times solar, similar to the absorption metallicity. We report low unobscured SFR ~1 Msun/yr given the large reservoir of molecular gas, and also modest obscured SFR=4.5(+4.4,-2.6) Msun/yr based on far-IR and sub-mm data. We determine mass components of the galaxy: log[M*/Msun] = 10.80(+0.07,-0.14), log[M mol-gas/Msun] = 10.37 +/-0.04, and log[M dust/Msun] = 8.45(+0.10,-0.30). Surprisingly, this HI absorption-selected galaxy has no equivalent objects in CO surveys of flux-selected samples. The galaxy falls off current scaling relations for the SFR to molecular gas mass and CO Tully-Fisher relation. Detailed comparison of kinematical components of the absorbing, ionized and molecular gas, combined with their spatial distribution, suggests that part of the CO gas is both kinematically and spatially de-coupled from the main galaxy. It is thus possible that a major star burst in the past could explain the wide CO profile as well as the low SFR. Support for this also comes from the SED favouring an instantaneous burst of age ~0.5 Gyr. Our survey will establish whether flux-selected surveys of molecular gas are missing a key stage in the evolution of galaxies and their conversion of gas to stars.


title: Proposal for the Detection of Magnetic Monopoles in Spin Ice via Nanoscale Magnetometry

abstract: We present a proposal for applying nanoscale magnetometry to the search for magnetic monopoles in the spin ice materials holmium and dysprosium titanate. Employing Monte Carlo simulations of the dipolar spin ice model, we find that when cooled to below $1.5\,$K these materials exhibit a sufficiently low monopole density to enable the direct observation of magnetic fields from individual monopoles. At these temperatures we demonstrate that noise spectroscopy can capture the intrinsic fluctuations associated with monopole dynamics, allowing one to isolate the qualitative effects associated with both the Coulomb interaction between monopoles and the topological constraints implied by Dirac strings. We describe in detail three different nanoscale magnetometry platforms (muon spin rotation, nitrogen vacancy defects, and nanoSQUID arrays) that can be used to detect monopoles in these experiments, and analyze the advantages of each.


title: Revisiting the Gluon Spectrum in the Boost-Invariant Glasma from a Semi-Analytic Approach

abstract: In high energy heavy-ion collisions, the degrees of freedom at the very early stage can be effectively represented by strong classical gluonic fields within the Color Glass Condensate framework. As the system expands, the strong gluonic fields eventually become weak such that an equivalent description using the gluonic particle degrees of freedom starts to become valid. We revisit the spectrum of these gluonic particles by solving the classical Yang-Mills equations semi-analytically with the solutions having the form of power series expansions in the proper time. We propose a different formula for the gluon spectrum which is consistent with energy density during the whole time evolution. We find that the chromo-electric fields have larger contributions to the gluon spectrum than the chromo-magnetic fields do. Furthermore, the large momentum modes take less time to reach the weak-field regime while smaller momentum modes take more time. The resulting functional form of the gluon spectrum is exponential in nature and the spectrum is close to a thermal distrubtion with effective temperatures around $0.6$ to $0.9\, Q_s$ late in the Glasma evolution. The sensitiveness of the gluon spectrum to the infrared and the ultraviolet cut-offs are discussed.


title: Photobiological effects at Earth's surface following a 50 pc Supernova

abstract: We investigated the potential biological impacts at Earth's surface of stratospheric O3 depletion caused by nearby supernovae known to have occurred about 2.5 and 8 million years ago at about 50 pc distance. New and previously published atmospheric chemistry modeling results were combined with radiative transfer modeling to determine changes in surface-level Solar irradiance and biological responses. We find that UVB irradiance is increased by a factor of 1.1 to 2.8, with large variation in latitude, and seasonally at high latitude regions. Changes in UVA and PAR (visible light) are much smaller. DNA damage (in vitro) is increased by factors similar to UVB, while other biological impacts (erythema, skin cancer, cataracts, marine phytoplankton photosynthesis inhibition, and plant damage) are increased by smaller amounts. We conclude that biological impacts due to increased UV irradiance in this SN case are not mass-extinction level, but might be expected to contribute to changes in species abundances; this result fits well with species turnover observed around the Pliocene-Pleistocene boundary.


title: Values of quadratic forms on quasicrystals and a related problem

abstract: In this paper we study the set of values of quadratic form at points of a cut and project set. We will establish conditions which ensure that the set of values is dense. Our methods involve homogeneous dynamics and we will prove a orbit closure classification type result in this setting. This result has additional applications. In particular, we use it to study the set of integral values of a system consisting of a quadratic form and several linear forms.


title: Torsion of rational elliptic curves over the maximal abelian extension of Q

abstract: Let $E$ be an elliptic curve defined over $\mathbb{Q}$, and let $\mathbb{Q}^{ab}$ be the maximal abelian extension of $\mathbb{Q}$. In this article we classify the groups that can arise as $E(\mathbb{Q}^{ab})_{\text{tors}}$ up to isomorphism. The method illustrates techniques for finding explicit models of modular curves of mixed level structure. Moreover we provide an explicit algorithm to compute $E(\mathbb{Q}^{ab})_{\text{tors}}$ for any elliptic curve $E/\mathbb{Q}$.


title: Combinatorial cost: a coarse setting

abstract: The main inspiration for this paper is a paper by Elek where he introduces combinatorial cost for graph sequences. We show that having cost equal to 1 and hyperfiniteness are coarse invariants. We also show `cost-1' for box spaces behaves multiplicatively when taking subgroups. We show that graph sequences coming from Farber sequences of a group have property A if and only if the group is amenable. The same is true for hyperfiniteness. This generalises a theorem by Elek. Furthermore we optimise this result when Farber sequences are replaced by sofic approximations. In doing so we introduce a new concept: property almost-A.


title: Intelligent Parameter Tuning in Optimization-based Iterative CT Reconstruction via Deep Reinforcement Learning

abstract: A number of image-processing problems can be formulated as optimization problems. The objective function typically contains several terms specifically designed for different purposes. Parameters in front of these terms are used to control the relative weights among them. It is of critical importance to tune these parameters, as quality of the solution depends on their values. Tuning parameter is a relatively straightforward task for a human, as one can intelligently determine the direction of parameter adjustment based on the solution quality. Yet manual parameter tuning is not only tedious in many cases, but becomes impractical when a number of parameters exist in a problem. Aiming at solving this problem, this paper proposes an approach that employs deep reinforcement learning to train a system that can automatically adjust parameters in a human-like manner. We demonstrate our idea in an example problem of optimization-based iterative CT reconstruction with a pixel-wise total-variation regularization term. We set up a parameter tuning policy network (PTPN), which maps an CT image patch to an output that specifies the direction and amplitude by which the parameter at the patch center is adjusted. We train the PTPN via an end-to-end reinforcement learning procedure. We demonstrate that under the guidance of the trained PTPN for parameter tuning at each pixel, reconstructed CT images attain quality similar or better than in those reconstructed with manually tuned parameters.


title: Performance Analysis for Massive MIMO Downlink with Low Complexity Approximate Zero-Forcing Precoding

abstract: Zero-forcing (ZF) precoding plays an important role for massive MIMO downlink due to its near optimal performance. However, the high computation cost of the involved matrix inversion hinders its application. In this paper, we adopt the first order Neumann series (NS) for a low-complexity approximation. By introducing a relaxation parameter jointly with one selected user's interference to others into the precondition matrix, we propose the identity-plus-column NS (ICNS) method. By further exploiting the multi-user diversity gain via choosing the user with the largest interference to others, the ordered ICNS method is also proposed. Moreover, the sum-rate approximations of the proposed ICNS method and the competitive existing identity matrix based NS (INS) method are derived in closed-form, based on which the performance loss of ICNS due to inversion approximation compared with ideal ZF and its performance gain over INS are explicitly analyzed for three typical massive MIMO scenarios. Finally, simulations verify our analytical results and also show that the proposed two designs achieve better performance-complexity tradeoff than ideal ZF and existing low-complexity ZF precodings for practical large antenna number, correlated channels and not-so-small loading factor.


title: Ringel duality for certain strongly quasi-hereditary algebras

abstract: We introduce quasi-hereditary endomorphism algebras defined over a new class of finite dimensional monomial algebras with a special ideal structure. The main result is a uniform formula describing the Ringel duals of these quasi-hereditary algebras. As special cases, we obtain a Ringel-duality formula for a family of strongly quasi-hereditary algebras arising from a type A configuration of projective lines in a rational, projective surface as recently introduced by Hille and Ploog, for certain Auslander-Dlab-Ringel algebras, and for Eiriksson and Sauter's nilpotent quiver algebras when the quiver has no sinks and no sources. We also recover Tan's result that the Auslander algebras of self-injective Nakayama algebras are Ringel self-dual.


title: Inferences on the Higgs Boson and Axion Masses through a Maximum Entropy Principle

abstract: The Maximum Entropy Principle (MEP) is a method that can be used to infer the value of an unknown quantity in a set of probability functions. In this work we review two applications of MEP: one giving a precise inference of the Higgs boson mass value; and the other one allowing to infer the mass of the axion. In particular, for the axion we assume that it has a decay channel into pairs of neutrinos, in addition to the decay into two photons. The Shannon entropy associated to an initial ensemble of axions decaying into photons and neutrinos is then built for maximization.


title: Snow Queen is Evil and Beautiful: Experimental Evidence for Probabilistic Contextuality in Human Choices

abstract: We present unambiguous experimental evidence for (quantum-like) probabilistic contextuality in psychology. All previous attempts to find contextuality in a psychological experiment were unsuccessful because of the gross violations of marginal selectivity in behavioral data, making the traditional mathematical tests developed in quantum mechanics inapplicable. In our crowdsourcing experiment respondents were making two simple choices: of one of two characters in a story (The Snow Queen by Hans Christian Andersen), and of one of two characteristics, such as Kind and Evil, so that the character and the characteristic chosen matched the story line. The formal structure of the experiment imitated that of the Einstein-Podolsky-Rosen paradigm in the Bohm-Bell version. Marginal selectivity was violated, indicating that the two choices were directly influencing each other, but the application of a mathematical test developed in the Contextuality-by-Default theory, extending the traditional quantum-mechanical test, indicated a strong presence of contextuality proper, not reducible to direct influences.


title: Competition and Complexity in Amphiphilic Polymer Morphology

abstract: We analyze the competitive evolution of codimension one and two morphologies within the $H^{-1}$ gradient flow of the strong Functionalized Cahn-Hilliard equation. On a slow time scale a sharp hypersurface reduction yields a degenerate Mullins-Sekerka evolution for both codimension one and two hypersurfaces, leading to a geometric flow that depends locally on curvatures couples to the dynamic value of the spatially constant far-field chemical potential. Both codimension one and two morphologies admit two classes of bifurcations, one leads to pearling, a short-wavelength in-plane modulation of interfacial width, the other flips motion by curvature to the locally-ill posed motion against curvature, which leads to fingering instabilities. We present a bifurcation diagram for the morphological competition, and compare our results quantitatively to simulations of the full system and qualitatively to simulations of self-consistent mean field models and laboratory experiments; illuminating the role of the pearling bifurcation in the development of complex network morphologies.


title: Super RSK correspondence with symmetry

abstract: Super RSK correspondence is a bijective correspondence between superbiwords and pairs of semistandard supertableaux. Such a bijection was given by Bonetti, Senato and Venezia, via an insertion algorithm closely related to Schensted insertion. Notably, the symmetry property satisfied by the classical RSK bijection holds only in special cases under this bijection. We present a new super RSK bijection, based on the mixed insertion process defined by Haiman, where the symmetry property holds in complete generality.


title: On the variance of radio interferometric calibration solutions: Quality-based Weighting Schemes

abstract: This paper investigates the possibility of improving radio interferometric images using an algorithm inspired by an optical method known as "lucky imaging", which would give more weight to the best-calibrated visibilities used to make a given image. A fundamental relationship between the statistics of interferometric calibration solution residuals and those of the image-plane pixels is derived in this paper. This relationship allows us to understand and describe the statistical properties of the residual image. In this framework, the noise-map can be described as the Fourier transform of the covariance between residual visibilities in a new differential Fourier plane. Image-plane artefacts can be seen as one realisation of the pixel covariance distribution, which can be estimated from the antenna gain statistics. Based on this relationship, we propose a means of improving images made with calibrated visibilities using weighting schemes. This improvement would occur after calibration, but before imaging - it is thus ideally used between major iterations of self-calibration loops. Applying the weighting scheme to simulated data improves the noise level in the final image at negligible computational cost.


title: The direct relation between the duration of magnetic reconnection and the evolution of GOES light curves in solar flares

abstract: GOES soft X-ray light curves are used to measure the timing and duration of solar flare emission. The timing and duration of the magnetic reconnection and subsequent energy release which drives solar flares are unknown, though the light curves are presumably related. It is therefore critical to understand the physics which connects the two: how does the time scale of reconnection produce an observed GOES light curve? In this work, we model the formation and expansion of an arcade of loops with a hydrodynamic model, which we then use to synthesize GOES light curves. We calculate the FWHM and the e-folding decay time of the light curves and compare them to the separation of the centroids of the two ribbons which the arcade spans, which is representative of the size scale of the loops. We reproduce a linear relation between the two, as found observationally in previous work. We show that this demonstrates a direct connection between the duration of energy release and the evolution of these light curves. We also show that the cooling processes of individual loops comprising the flare arcade directly affect the measured time scales. From the clear consistency between the observed and modeled linearity, we conclude that the primary factors that control the flare time scales are the duration of reconnection and the loop lengths.


title: The MUSE Hubble Ultra Deep Field Survey. IX. Evolution of galaxy merger fraction since z~6

abstract: We provide, for the first time, robust observational constraints on the galaxy major merger fraction up to $z\approx 6$ using spectroscopic close pair counts. Deep Multi Unit Spectroscopic Explorer (MUSE) observations in the Hubble Ultra Deep Field (HUDF) and Hubble Deep Field South (HDF-S) are used to identify 113 secure close pairs of galaxies among a parent sample of 1801 galaxies spread over a large redshift range ($0.2<z<6$) and stellar masses ($10^7-10^{11} M_\odot$), thus probing about 12 Gyr of galaxy evolution. Stellar masses are estimated from spectral energy distribution (SED) fitting over the extensive UV-to-NIR HST photometry available in these deep Hubble fields, adding Spitzer IRAC bands to better constrain masses for high-redshift ($z\geqslant 3$) galaxies. These stellar masses are used to isolate a sample of 54 major close pairs with a galaxy mass ratio limit of 1:6. Among this sample, 23 pairs are identified at high redshift ($z\geqslant 3$) through their Ly$\alpha$ emission. The sample of major close pairs is divided into five redshift intervals in order to probe the evolution of the merger fraction with cosmic time. Our estimates are in very good agreement with previous close pair counts with a constant increase of the merger fraction up to $z\approx 3$ where it reaches a maximum of 20%. At higher redshift, we show that the fraction slowly decreases down to about 10% at $z\approx6$. The sample is further divided into two ranges of stellar masses using either a constant separation limit of $10^{9.5} M_\odot$ or the median value of stellar mass computed in each redshift bin. Overall, the major close pair fraction for low-mass and massive galaxies follows the same trend. These new, homogeneous, and robust estimates of the major merger fraction since $z\approx6$ are in good agreement with recent predictions of cosmological numerical simulations.


title: Nonlinear quantum Langevin equations for bosonic modes in solid-state systems

abstract: Based on the experimental evidence that impurities contribute to the dissipation properties of solid-state open quantum systems, we provide here a description in terms of nonlinear quantum Langevin equations of the role played by two-level systems in the dynamics of a bosonic degree of freedom. Our starting point is represented by the description of the system/environment coupling in terms of coupling to two separate reservoirs, modelling the interaction with external bosonic modes and two level systems, respectively. Furthermore, we show how this model represents a specific example of a class of open quantum systems that can be described by nonlinear quantum Langevin equations. Our analysis offers a potential explanation of the parametric effects recently observed in circuit-QED cavity optomechanics experiments.


title: Non Uniform On Chip Power Delivery Network Synthesis Methodology

abstract: In this paper, we proposed a non-uniform power delivery network (PDN) synthesis methodology. It first constructs initial PDN using uniform approach. Then preliminary power integrity analysis is performed to derive IR-safe candidate window. Congestion map is obtained based global route congestion estimation. A self-adaptive non-uniform PDN synthesis is then performed to globally and locally optimize PDN over selected regions. The PDN synthesis is congestion-driven and IR- guarded. Experimental results show significant timing important in trade-off small PDN length reduction with no EM/IR impact. We further explored potential power savings using our non-uniform PDN synthesis methodology.


title: Unveiling the Link between Complexity and Symmetry: Statistical Asymmetry

abstract: The concept of complexity appears in virtually all areas of knowledge. Its intuitive meaning shares similarities across fields, but disagreements between its details hinders a general definition, leading to a plethora of proposed measurements. While each might be appropriated to the problems it addresses, the lack of an underlying fundamental principle prevents the development of a unified theory. Here it is shown that the statistics of the amount of symmetry broken by systems can be used as such unifying principle. A general methodology is outlined and explicit expressions are given for cases in which it can capture the behavior of the two main groups of complexities currently in use. The presented results demonstrate that statistical asymmetry is an appropriate foundation for characterizing the general concept of complexity.


title: Fractional Brownian motion with zero Hurst parameter: a rough volatility viewpoint

abstract: Rough volatility models are becoming increasingly popular in quantitative finance. In this framework, one considers that the behavior of the log-volatility process of a financial asset is close to that of a fractional Brownian motion with Hurst parameter around 0.1. Motivated by this, we wish to define a natural and relevant limit for the fractional Brownian motion when $H$ goes to zero. We show that once properly normalized, the fractional Brownian motion converges to a Gaussian random distribution which is very close to a log-correlated random field.


title: On Ordinal Invariants in Well Quasi Orders and Finite Antichain Orders

abstract: We investigate the ordinal invariants height, length, and width of well quasi orders (WQO), with particular emphasis on width, an invariant of interest for the larger class of orders with finite antichain condition (FAC). We show that the width in the class of FAC orders is completely determined by the width in the class of WQOs, in the sense that if we know how to calculate the width of any WQO then we have a procedure to calculate the width of any given FAC order. We show how the width of WQO orders obtained via some classical constructions can sometimes be computed in a compositional way. In particular, this allows proving that every ordinal can be obtained as the width of some WQO poset. One of the difficult questions is to give a complete formula for the width of Cartesian products of WQOs. Even the width of the product of two ordinals is only known through a complex recursive formula. Although we have not given a complete answer to this question we have advanced the state of knowledge by considering some more complex special cases and in particular by calculating the width of certain products containing three factors. In the course of writing the paper we have discovered that some of the relevant literature was written on cross-purposes and some of the notions re-discovered several times. Therefore we also use the occasion to give a unified presentation of the known results.


title: A counterexample to Stein's Equi-n-square Conjecture

abstract: In 1975 Stein conjectured that in every $n\times n$ array filled with the numbers $1, \dots, n$ with every number occuring exactly $n$ times, there is a partial transversal of size $n-1$. In this note we show that this conjecture is false by constructing such arrays without partial transverals of size $n-\frac{1}{42}\ln n$.


title: A completely integrable flow of star-shaped curves on the light cone in Lorentzian $\mathbb{R}^4$

abstract: In this paper we prove that the space of differential invariants for curves with arc-length parameter in the light cone of Lorentzian $\mathbb{R}^4$, invariants under the centro-affine action of the Lorentzian group, is Poisson equivalent to the space of conformal differential invariants for curves in the M\"obius sphere. We use this relation to find realizations of solutions of a complexly coupled system of KdV equations as flows of curves in the cone.


title: Search for new physics in events with a leptonically decaying Z boson and a large transverse momentum imbalance in proton-proton collisions at $\sqrt{s} $ = 13 TeV

abstract: A search for new physics in events with a Z boson produced in association with large missing transverse momentum at the LHC is presented. The search is based on the 2016 data sample of proton-proton collisions recorded with the CMS experiment at $\sqrt{s} = $ 13 TeV, corresponding to an integrated luminosity of 35.9 fb$^{-1}$. The results of this search are interpreted in terms of a simplified model of dark matter production via spin-0 or spin-1 mediators, a scenario with a standard-model-like Higgs boson produced in association with the Z boson and decaying invisibly, a model of unparticle production, and a model with large extra spatial dimensions. No significant deviations from the background expectations are found, and limits are set on relevant model parameters, significantly extending the results previously achieved in this channel.


title: Statistical Challenges in Modeling Big Brain Signals

abstract: Brain signal data are inherently big: massive in amount, complex in structure, and high in dimensions. These characteristics impose great challenges for statistical inference and learning. Here we review several key challenges, discuss possible solutions, and highlight future research directions.


title: Block-modified Wishart matrices: the easy case

abstract: Associated to any complex Wishart matrix $W$ of parameters $(dn,dm)$ and any linear map $\varphi:M_n(\mathbb C)\to M_n(\mathbb C)$ is the "block-modified" matrix $\tilde{W}=(id\otimes\varphi)W$. Following some previous work with Nechita, we study here the asymptotic $*$-distribution of $\tilde{W}$, in the $d\to\infty$ limit, in the case where the modification map $\varphi$ is "easy", or more generally super-easy, in the quantum algebra/representation theory sense. Under suitable assumptions on $\varphi$ we obtain in this way a compound free Poisson law.


title: Quantum algebra from generalized q-Hermite polynomials

abstract: In this paper, we discuss new results related to the generalized discrete $q$-Hermite II polynomials $ \tilde h_{n,\alpha}(x;q)$, introduced by Mezlini et al. in 2014. Our aim is to give a continuous orthogonality relation, a $q$-integral representation and an evaluation at unity of the Poisson kernel, for these polynomials. Furthermore, we introduce $q$-Schr\"{o}dinger operators and give the raising and lowering operator algebra corresponding to these polynomials. Our results generate a new explicit realization of the quantum algebra $\mathsf{su}_{q}(1, 1)$, using the generators associated with a $q$-deformed generalized para-Bose oscillator.


title: Spherical Collapse Models with Clustered Dark Energy

abstract: We investigate the clustering effect of dark energy (DE) in the formation of galaxy clusters using the spherical collapse model. Assuming a fully clustered DE component, the spherical overdense region is treated as an isolated system which conserves the energy separately for both matter and DE inside the spherical region. Then, by introducing a parameter $r$ to characterize the degree of DE clustering, which is defined by the nonlinear density contrast ratio of matter to DE at turnaround in the recollapsing process, i.e. $r\equiv \nld_{\de,\ta}/\nld_{\m,\ta}$, we are able to uniquely determine the spherical collapsing process and hence obtain the virialized overdensity $\Dvir$ through a proper virialization scheme. Estimation of the virialized overdensities from current observation on galaxy clusters suggests that $0.5 < r < 0.8$ at $1\sigma$ level for the clustered DE with $w < -0.9$. Also, we compare our method to the linear perturbation theory that deals with the growth of DE perturbation at early times. While both results are consistent with each other, our method is practically simple and it shows that the collapse process is rather independent of initial DE perturbation and its evolution at early times.


title: Hierarchical Representations for Efficient Architecture Search

abstract: We explore efficient neural architecture search methods and show that a simple yet powerful evolutionary algorithm can discover new architectures with excellent performance. Our approach combines a novel hierarchical genetic representation scheme that imitates the modularized design pattern commonly adopted by human experts, and an expressive search space that supports complex topologies. Our algorithm efficiently discovers architectures that outperform a large number of manually designed models for image classification, obtaining top-1 error of 3.6% on CIFAR-10 and 20.3% when transferred to ImageNet, which is competitive with the best existing neural architecture search approaches. We also present results using random search, achieving 0.3% less top-1 accuracy on CIFAR-10 and 0.1% less on ImageNet whilst reducing the search time from 36 hours down to 1 hour.


title: Geostatistical inference in the presence of geomasking: a composite-likelihood approach

abstract: In almost any geostatistical analysis, one of the underlying, often implicit, modelling assump- tions is that the spatial locations, where measurements are taken, are recorded without error. In this study we develop geostatistical inference when this assumption is not valid. This is often the case when, for example, individual address information is randomly altered to provide pri- vacy protection or imprecisions are induced by geocoding processes and measurement devices. Our objective is to develop a method of inference based on the composite likelihood that over- comes the inherent computational limits of the full likelihood method as set out in Fanshawe and Diggle (2011). Through a simulation study, we then compare the performance of our proposed approach with an N-weighted least squares estimation procedure, based on a corrected version of the empirical variogram. Our results indicate that the composite-likelihood approach outper- forms the latter, leading to smaller root-mean-square-errors in the parameter estimates. Finally, we illustrate an application of our method to analyse data on malnutrition from a Demographic and Health Survey conducted in Senegal in 2011, where locations were randomly perturbed to protect the privacy of respondents.


title: Signatures of $\alpha$ clustering in ultra-relativistic collisions with light nuclei

abstract: We explore possible observable signatures of $\alpha$ clustering of light nuclei in ultra-relativistic nuclear collisions involving ${}^{7,9}$Be, ${}^{12}$C, and ${}^{16}$O. The clustering leads to specific spatial correlations of the nucleon distributions in the ground state, which are manifest in the earliest stage of the ultra-high energy reaction. The formed initial state of the fireball is sensitive to these correlations, and the effect influences, after the collective evolution of the system, the hadron production in the final stage. Specifically, we study effects on the harmonic flow in collisions of light clustered nuclei with a heavy target (${}^{208}$Pb), showing that measures of the elliptic flow are sensitive to clusterization in ${}^{7,9}$Be, whereas triangular flow is sensitive to clusterization in ${}^{12}$C and ${}^{16}$O. Specific predictions are made for model collisions at the CERN SPS energies. In another exploratory development we also examine the proton-beryllium collisions, where the $3/2^-$ ground state of ${}^{7,9}$Be nucleus is polarized by an external magnetic field. Clusterization leads to multiplicity distributions of participant nucleons which depend on the orientation of the polarization with respect to the collision axis, as well as on the magnetic number of the state. The obtained effects on multiplicities reach a factor of a few for collisions with a large number of participant nucleons.


title: Sampling and multilevel coarsening algorithms for fast matrix approximations

abstract: This paper addresses matrix approximation problems for matrices that are large, sparse and/or that are representations of large graphs. To tackle these problems, we consider algorithms that are based primarily on coarsening techniques, possibly combined with random sampling. A multilevel coarsening technique is proposed which utilizes a hypergraph associated with the data matrix and a graph coarsening strategy based on column matching. Theoretical results are established that characterize the quality of the dimension reduction achieved by a coarsening step, when a proper column matching strategy is employed. We consider a number of standard applications of this technique as well as a few new ones. Among the standard applications we first consider the problem of computing the partial SVD for which a combination of sampling and coarsening yields significantly improved SVD results relative to sampling alone. We also consider the Column subset selection problem, a popular low rank approximation method used in data related applications, and show how multilevel coarsening can be adapted for this problem. Similarly, we consider the problem of graph sparsification and show how coarsening techniques can be employed to solve it. Numerical experiments illustrate the performances of the methods in various applications.


title: Testing the photon-number statistics of a quantum key distribution light source

abstract: A commonly held tenet is that lasers well above threshold emit photons in a coherent state, which follow a Poissonian statistics when measured in photon number. This feature is often exploited to build quantum-based random number generators or to derive the secure key rate of quantum key distribution systems. Hence the photon number distribution of the light source can directly impact the randomness and the security distilled from such devices. Here, we propose a method based on measuring correlation functions to experimentally characterise a light source's photon statistics and use it in the estimation of a quantum key distribution system's key rate. This promises to be a useful tool for the certification of quantum-related technologies.


title: Data, Depth, and Design: Learning Reliable Models for Skin Lesion Analysis

abstract: Deep learning fostered a leap ahead in automated skin lesion analysis in the last two years. Those models are expensive to train and difficult to parameterize. Objective: We investigate methodological issues for designing and evaluating deep learning models for skin lesion analysis. We explore 10 choices faced by researchers: use of transfer learning, model architecture, train dataset, image resolution, type of data augmentation, input normalization, use of segmentation, duration of training, additional use of SVMs, and test data augmentation. Methods: We perform two full factorial experiments, for five different test datasets, resulting in 2560 exhaustive trials in our main experiment, and 1280 trials in our assessment of transfer learning. We analyze both with multi-way ANOVA. We use the exhaustive trials to simulate sequential decisions and ensembles, with and without the use of privileged information from the test set. Results -- main experiment: Amount of train data has disproportionate influence, explaining almost half the variation in performance. Of the other factors, test data augmentation and input resolution are the most influential. Deeper models, when combined, with extra data, also help. -- transfer experiment: Transfer learning is critical, its absence brings huge performance penalties. -- simulations: Ensembles of models are the best option to provide reliable results with limited resources, without using privileged information and sacrificing methodological rigor. Conclusions and Significance: Advancing research on automated skin lesion analysis requires curating larger public datasets. Indirect use of privileged information from the test set to design the models is a subtle, but frequent methodological mistake that leads to overoptimistic results. Ensembles of models are a cost-effective alternative to the expensive full-factorial and to the unstable sequential designs.


title: Sustainable Green Networking: Exploiting Degrees of Freedom towards Energy-Efficient 5G Systems

abstract: The carbon footprint concern in the development and deployment of 5G new radio systems has drawn the attention to several stakeholders. In this article, we analyze the critical power consuming component of all candidate 5G system architectures-the power amplifier (PA)-and propose PA-centric resource management solutions for green 5G communications. We discuss the impact of ongoing trends in cellular communications on sustainable green networking and analyze two communications architectures that allow exploiting the extra degrees-of-freedom (DoF) from multi-antenna and massive antenna deployments: small cells/distributed antenna network and massive MIMO. For small cell systems with a moderate number of antennas, we propose a peak to average power ratio-aware resource allocation scheme for joint orthogonal frequency and space division multiple access. For massive MIMO systems, we develop a highly parallel recurrent neural network for energy-efficient precoding. Simulation results for representative 5G deployment scenarios demonstrate an energy efficiency improvement of one order of magnitude or higher with respect to current state-of-the-art solutions.


title: Optimizing S-shaped utility and implications for risk management

abstract: We consider market players with tail-risk-seeking behaviour as exemplified by the S-shaped utility introduced by Kahneman and Tversky. We argue that risk measures such as value at risk (VaR) and expected shortfall (ES) are ineffective in constraining such players. We show that, in many standard market models, product design aimed at utility maximization is not constrained at all by VaR or ES bounds: the maximized utility corresponding to the optimal payoff is the same with or without ES constraints. By contrast we show that, in reasonable markets, risk management constraints based on a second more conventional concave utility function can reduce the maximum S-shaped utility that can be achieved by the investor, even if the constraining utility function is only rather modestly concave. It follows that product designs leading to unbounded S-shaped utilities will lead to unbounded negative expected constraining utilities when measured with such conventional utility functions. To prove these latter results we solve a general problem of optimizing an investor expected utility under risk management constraints where both investor and risk manager have conventional concave utility functions, but the investor has limited liability. We illustrate our results throughout with the example of the Black--Scholes option market. These results are particularly important given the historical role of VaR and that ES was endorsed by the Basel committee in 2012--2013.


title: Engineering reaction-diffusion networks with properties of neural tissue

abstract: We present an experimental system of networks of coupled non-linear chemical reactors, which we theoretically model within a reaction-diffusion framework. The networks consist of patterned arrays of diffusively coupled nanoliter-scale reactors containing the Belousov- Zhabotinsky (BZ) reaction. Microfluidic fabrication techniques are developed that provide the ability to vary the network topology, the reactor coupling strength and offer the freedom to choose whether an arbitrary reactor is inhibitory or excitatory coupled to its neighbor. This versatile experimental and theoretical framework can be used to create a wide variety of chemical networks. Here we design, construct and characterize chemical networks that achieve the complexity of central pattern generators (CPGs), which are found in the autonomic nervous system of a variety of organisms.


title: Pole, Pericenter, and Nodes of the Interstellar Minor Body A/2017 U1

abstract: Here, we explore some peculiar orbital features of the recently discovered asteroid A/2017 U1, which is a clear outlier when considering the average value of the eccentricity of known hyperbolic comets. As for the orientation of its orbit in space, the orbital plane of A/2017 U1 seems to be away from any obvious clusters present for this population. The orbital nodes of A/2017 U1 are well away from the paths of the planets of the Solar System and the Sun. All these orbital properties appear to confirm A/2017 U1 as the first known interstellar asteroid.


title: Bases for cluster algebras from orbifolds with one marked point

abstract: We generalize the construction of the bangle, band and the bracelet bases for cluster algebras from orbifolds to the case where there is only one marked point on the boundary.


title: SMA and ALMA Studies of Disk- and Planet Formation around Low-mass Protostars

abstract: We report our current SMA and ALMA studies of disk and planet formation around protostars. We have revealed that $r \gtrsim$100 AU scale disks in Keplerian rotation are ubiquitous around Class I sources. These Class I Keplerian disks are often embedded in rotating and infalling protostellar envelopes. The infalling speeds of the protostellar envelopes are typically $\sim$ 3-times smaller than the free-fall velocities, and the rotational profiles follow the $r^{-1}$ profile, that is, rotation with the conserved specific angular momentum. Our latest high-resolution ($\sim$0$\farcs$5) ALMA studies, as well as the other studies in the literature, have unveiled that $r \sim$100-AU scale Keplerian disks are also present in several Class 0 protostars, while in the other Class 0 sources the inferred upper limits of the Keplerian disks are very small ($r \lessim$20 AU). Our recent data analyses of the ALMA long baseline data of the Class I-II source HL Tau have revealed gaps in molecular gas as well as in dust in the surrounding disk, suggesting the presence of sub-Jovian planets in the disk. These results imply that disk and planet formation should be completed in the protostellar stage.


title: Sufficient Conditions for the Controllability of Wave Equations with a Transmission Condition at the Interface

abstract: We consider waves travelling in two different mediums each endowed with a different constant speed of propagation. At the interface between the two mediums, the refraction of the rays of the optic geometry is described by the Snell-Descartes law. We provide sufficient conditions on the geometry of the mediums and on the speed of propagation for the boundary controllability.


title: Attacking Binarized Neural Networks

abstract: Neural networks with low-precision weights and activations offer compelling efficiency advantages over their full-precision equivalents. The two most frequently discussed benefits of quantization are reduced memory consumption, and a faster forward pass when implemented with efficient bitwise operations. We propose a third benefit of very low-precision neural networks: improved robustness against some adversarial attacks, and in the worst case, performance that is on par with full-precision models. We focus on the very low-precision case where weights and activations are both quantized to $\pm$1, and note that stochastically quantizing weights in just one layer can sharply reduce the impact of iterative attacks. We observe that non-scaled binary neural networks exhibit a similar effect to the original defensive distillation procedure that led to gradient masking, and a false notion of security. We address this by conducting both black-box and white-box experiments with binary models that do not artificially mask gradients.


title: Non-Euclidean geometry, nontrivial topology and quantum vacuum effects

abstract: Space out of a topological defect of the Abrikosov-Nielsen-Olesen vortex type is locally flat but non-Euclidean. If a spinor field is quantized in such a space, then a variety of quantum effects is induced in the vacuum. Basing on the continuum model for long-wavelength electronic excitations, originating in the tight-binding approximation for the nearest neighbor interaction of atoms in the crystal lattice, we consider quantum ground state effects in monolayer structures warped into nanocones by a disclination; the nonzero size of the disclination is taken into account, and a boundary condition at the edge of the disclination is chosen to ensure self-adjointness of the Dirac-Weyl Hamiltonian operator. In the case of carbon nanocones, we find circumstances when the quantum ground state effects are independent of the boundary parameter and the disclination size.

categories: hep-th cond-mat.mtrl-sci cond-mat.str-el math-ph math.MP


title: Atomic and electronic structure of doped Si(111)(2root3x2root3)R30(degrees)-Sn interfaces

abstract: The hole doped Si(111)(2root3x2root3)R30(degrees)-Sn interface exhibits a symmetry-breaking insulator-insulator transition below 100 K that appears to be triggered by electron tunneling into the empty surface-state bands. No such transition is seen in electron-doped systems. To elucidate the nature and driving force of this phenomenon, the structure of the interface must be resolved. Here we report on an extensive experimental and theoretical study, including scanning tunneling microscopy and spectroscopy (STM/STS), dynamical low-energy electron diffraction (LEED) analysis, and density functional theory (DFT) calculations, to elucidate the structure of this interface. We consider six different structure models, three of which have been proposed before, and conclude that only two of them can account for the majority of experimental data. One of them is the model according to Tornevik et al. [C. Tornevik et al., Phys. Rev. B 44, 13144 (1991)] with a total Sn coverage of 14/12 monolayers (ML). The other is the 'revised trimer model' with a total Sn coverage of 13/12 ML, introduced in this work. These two models are very difficult to discriminate on the basis of DFT or LEED alone, but STS data clearly point toward the T\"ornevik model as the most viable candidate among the models considered here. The STS data also provide additional insights regarding the electron-injection-driven phase transformation as well as the critical role of valence band holes in this process. Similar processes may occur at other metal/semiconductor interfaces, provided they are non-metallic and can be doped. This could open up a new pathway toward the creation of novel surface phases with potentially very interesting and desirable electronic properties.


title: An Autonomous Programmable Actuator and Shape Reconfigurable Structures using Bistability and Shape Memory Polymers

abstract: Autonomous deployment and shape reconfiguration of structures is a crucial field of research in space exploration with emerging applications in the automotive, building and biomedical industries. Challenges in achieving autonomy include: bulky energy sources, imprecise deployment, jamming of components and lack of structural integrity. Leveraging advances in the fields of shape memory polymers, bistability and 3D multi-material printing, we present a 3D printed programmable actuator that enables the autonomous deployment and shape reconfiguration of structures activated though surrounding temperature change. Using a shape memory polymer as the temperature controllable energy source and a bistable mechanism as the linear actuator and force amplifier, the structures achieve precise geometric activation and quantifiable load bearing capacity. The proposed unit actuator integrates these two components and is designed to be assembled into larger deployable and shape reconfigurable structures. First, we demonstrate that the activation of the unit actuator can be sequenced by tailoring each shape memory polymer to a different activation time. Next, by changing the configuration of the actuator, we demonstrate an initially flat surface that transforms into a pyramid or a hyperbolic paraboloid, thus demonstrating a multi-state structure. Load bearing capability is demonstrated for both during activation and in the operating state.


title: Transient Co-orbitals of Venus: An Update

abstract: Venus has no known satellites, but has four known co-orbitals: (322756) 2001 CK32, 2002 VE68, 2012 XE133, and 2013 ND15. Here, we present numerical evidence suggesting that 2015 WZ12 is a possible Venus co-orbital; it might have been until recently a transient Trojan. Follow-up observations of this target in the near future will be difficult, though.


title: Phase Diagram of Stability for Massive Scalars in Anti-de Sitter Spacetime

abstract: We diagram the behavior of 5-dimensional anti-de Sitter spacetime against horizon formation in the gravitational collapse of a scalar field, treating the scalar field mass and width of initial data as free parameters, which we call the stability phase diagram. We find that the class of stable initial data becomes larger and shifts to smaller widths as the field mass increases. In addition to classifying initial data as stable or unstable, we identify two other classes based on nonperturbative behavior. The class of metastable initial data forms a horizon over longer time scales than suggested by the lowest order perturbation theory at computationally accessible amplitudes, and irregular initial data can exhibit non-monotonic and possibly chaotic behavior in the horizon formation times. Our results include evidence for chaotic behavior even in the collapse of a massless scalar field.


title: A Unified View of Piecewise Linear Neural Network Verification

abstract: The success of Deep Learning and its potential use in many safety-critical applications has motivated research on formal verification of Neural Network (NN) models. Despite the reputation of learned NN models to behave as black boxes and the theoretical hardness of proving their properties, researchers have been successful in verifying some classes of models by exploiting their piecewise linear structure and taking insights from formal methods such as Satisifiability Modulo Theory. These methods are however still far from scaling to realistic neural networks. To facilitate progress on this crucial area, we make two key contributions. First, we present a unified framework that encompasses previous methods. This analysis results in the identification of new methods that combine the strengths of multiple existing approaches, accomplishing a speedup of two orders of magnitude compared to the previous state of the art. Second, we propose a new data set of benchmarks which includes a collection of previously released testcases. We use the benchmark to provide the first experimental comparison of existing algorithms and identify the factors impacting the hardness of verification problems.


title: Series for $1/\pi$ of signature 20

abstract: Properties of theta functions and Eisenstein series dating to Jacobi and Ramanujan are used to deduce differential equations associated with McKay Thompson series of level 20. These equations induce expansions for modular forms of level 20 in terms of modular functions.The theory of singular values is applied to derive expansions for $1/\pi$ of signature $20$ analogous to those formulated by Ramanujan.


title: Almost instant brain atlas segmentation for large-scale studies

abstract: Large scale studies of group differences in healthy controls and patients and screenings for early stage disease prevention programs require processing and analysis of extensive multisubject datasets. Complexity of the task increases even further when segmenting structural MRI of the brain into an atlas with more than 50 regions. Current automatic approaches are time-consuming and hardly scalable; they often involve many error prone intermediate steps and don't utilize other available modalities. To alleviate these problems, we propose a feedforward fully convolutional neural network trained on the output produced by the state of the art models. Incredible speed due to available powerful GPUs neural network makes this analysis much easier and faster (from $>10$ hours to a minute). The proposed model is more than two orders of magnitudes faster than the state of the art and yet as accurate. We have evaluated the network's performance by comparing it with the state of the art in the task of differentiating region volumes of healthy controls and patients with schizophrenia on a dataset with 311 subjects. This comparison provides a strong evidence that speed did not harm the accuracy. The overall quality may also be increased by utilizing multi-modal datasets (not an easy task for other models) by simple adding more modalities as an input. Our model will be useful in large-scale studies as well as in clinical care solutions, where it can significantly reduce delay between the patient screening and the result.


title: Likelihood Analysis of the Sub-GUT MSSM in Light of LHC 13-TeV Data

abstract: We describe a likelihood analysis using MasterCode of variants of the MSSM in which the soft supersymmetry-breaking parameters are assumed to have universal values at some scale $M_{in}$ below the supersymmetric grand unification scale $M_{GUT}$, as can occur in mirage mediation and other models. In addition to $M_{in}$, such `sub-GUT' models have the 4 parameters of the CMSSM, namely a common gaugino mass $m_{1/2}$, a common soft supersymmetry-breaking scalar mass $m_0$, a common trilinear mixing parameter $A$ and the ratio of MSSM Higgs vevs $\tan\beta$, assuming that the Higgs mixing parameter $\mu > 0$. We take into account constraints on strongly- and electroweakly-interacting sparticles from $\sim 36$/fb of LHC data at 13 TeV and the LUX and 2017 PICO, XENON1T and PandaX-II searches for dark matter scattering, in addition to the previous LHC and dark matter constraints as well as full sets of flavour and electroweak constraints. We find a preference for $M_{in} \sim 10^5$ to $10^9$ GeV, with $M_{in} \sim M_{GUT}$ disfavoured by $\Delta \chi^2 \sim 3$ due to the ${\rm BR}(B_{s, d} \to \mu^+\mu^-)$ constraint. The lower limits on strongly-interacting sparticles are largely determined by LHC searches, and similar to those in the CMSSM. We find a preference for the LSP to be a Bino or Higgsino with $\tilde{\chi^0_1} \sim 1$ TeV, with annihilation via heavy Higgs bosons $H/A$ and stop coannihilation, or chargino coannihilation, bringing the cold dark matter density into the cosmological range. We find that spin-independent dark matter scattering is likely to be within reach of the planned LUX-Zeplin and XENONnT experiments. We probe the impact of the $(g-2)_\mu$ constraint, finding similar results whether or not it is included.


title: Level $17$ Ramanujan-Sato series

abstract: Two level 17 modular functions $$ r=q^{2}\prod_{n=1}^{\infty}(1-q^{n})^{\left(\frac{n}{17}\right)},\quad s=q^{2}\prod_{n=1}^{\infty}\frac{(1-q^{17n})^{3}}{(1-q^{n})^{3}}, $$ are used to construct a new class of Ramanujan-Sato series for $1/\pi$. The expansions are induced by modular identities similar to those level of 5 and 13 appearing in Ramanujan's Notebooks. A complete list of rational and quadratic series corresponding to singular values of the parameters is derived.


title: Locally stationary spatio-temporal interpolation of Argo profiling float data

abstract: Argo floats measure seawater temperature and salinity in the upper 2,000 m of the global ocean. Statistical analysis of the resulting spatio-temporal dataset is challenging due to its nonstationary structure and large size. We propose mapping these data using locally stationary Gaussian process regression where covariance parameter estimation and spatio-temporal prediction are carried out in a moving-window fashion. This yields computationally tractable nonstationary anomaly fields without the need to explicitly model the nonstationary covariance structure. We also investigate Student-$t$ distributed fine-scale variation as a means to account for non-Gaussian heavy tails in ocean temperature data. Cross-validation studies comparing the proposed approach with the existing state-of-the-art demonstrate clear improvements in point predictions and show that accounting for the nonstationarity and non-Gaussianity is crucial for obtaining well-calibrated uncertainties. This approach also provides data-driven local estimates of the spatial and temporal dependence scales for the global ocean which are of scientific interest in their own right.


title: Topology of polyhedral products over simplicial multiwedges

abstract: We prove that certain conditions on multigraded Betti numbers of a simplicial complex $K$ imply existence of a higher Massey product in cohomology of a moment-angle-complex $\mathcal Z_K$, which contains a unique element (a strictly defined product). Using the simplicial multiwedge construction, we find a family $\mathcal{F}$ of polyhedral products being smooth closed manifolds such that for any $l,r\geq 2$ there exists an $l$-connected manifold $M\in\mathcal F$ with a nontrivial strictly defined $r$-fold Massey product in $H^{*}(M)$. As an application to homological algebra, we determine a wide class of triangulated spheres $K$ such that a nontrivial higher Massey product of any order may exist in Koszul homology of their Stanley--Reisner rings. As an application to rational homotopy theory, we establish a combinatorial criterion for a simple graph $\Gamma$ to provide a (rationally) formal generalized moment-angle manifold $\mathcal Z_{P}^{J}=(D^{2j_{i}},S^{2j_{i}-1})^{\partial P^*}$, $J=(j_{1},\ldots,j_m)$ over a graph-associahedron $P=P_{\Gamma}$ and compute all the diffeomorphism types of formal moment-angle manifolds over graph-associahedra.


title: Early prediction of the duration of protests using probabilistic Latent Dirichlet Allocation and Decision Trees

abstract: Protests and agitations are an integral part of every democratic civil society. In recent years, South Africa has seen a large increase in its protests. The objective of this paper is to provide an early prediction of the duration of protests from its free flowing English text description. Free flowing descriptions of the protests help us in capturing its various nuances such as multiple causes, courses of actions etc. Next we use a combination of unsupervised learning (topic modeling) and supervised learning (decision trees) to predict the duration of the protests. Our results show a high degree (close to 90%) of accuracy in early prediction of the duration of protests.We expect the work to help police and other security services in planning and managing their resources in better handling protests in future.


title: An overview of DarkBit, the GAMBIT dark matter module

abstract: In this conference paper, I give an overview of the capabilities of DarkBit, a module of the GAMBIT global fitting code that calculates a range of dark matter observables and corresponding experimental likelihood functions. Included in the code are limits from the dark matter relic density, multiple direct detection experiments, and indirect searches in gamma-rays and neutrinos. I discuss the capabilities of the code, and then present recent results of GAMBIT scans of the parameter space of the minimal supersymmetric standard model, with a focus on sensitivities of future dark matter searches to the current best fit regions.


title: Fixing a Broken ELBO

abstract: Recent work in unsupervised representation learning has focused on learning deep directed latent-variable models. Fitting these models by maximizing the marginal likelihood or evidence is typically intractable, thus a common approximation is to maximize the evidence lower bound (ELBO) instead. However, maximum likelihood training (whether exact or approximate) does not necessarily result in a good latent representation, as we demonstrate both theoretically and empirically. In particular, we derive variational lower and upper bounds on the mutual information between the input and the latent variable, and use these bounds to derive a rate-distortion curve that characterizes the tradeoff between compression and reconstruction accuracy. Using this framework, we demonstrate that there is a family of models with identical ELBO, but different quantitative and qualitative characteristics. Our framework also suggests a simple new method to ensure that latent variable models with powerful stochastic decoders do not ignore their latent code.