title: Using UDP for Internet Transport Evolution

abstract: The increasing use of middleboxes (e.g., NATs, firewalls) in the Internet has made it harder and harder to deploy new transport or higher layer protocols, or even extensions to existing ones. Current work to address this Internet transport ossification has led to renewed interest in UDP as an encapsulation for making novel transport protocols deployable in the Internet. Examples include Google's QUIC and the WebRTC data channel. The common assumption made by these approaches is that encapsulation over UDP works in the present Internet. This paper presents a measurement study to examine this assumption, and provides guidance for protocol design based on our measurements. The key question is "can we run new transport protocols for the Internet over UDP?" We find that the answer is largely "yes": UDP works on most networks, and impairments are generally confined to access networks. This allows relatively simple fallback strategies to work around it. Our answer is based on a twofold methodology. First, we use the RIPE Atlas platform to basically check UDP connectivity and first-packet latency. Second, we deploy copycat, a new tool for comparing TCP loss, latency, and throughput with UDP by generating TCP-shaped traffic with UDP headers.


title: Massive spinons in $S=1/2$ spin chains: spinon-pair operator representation

abstract: Spinons are among the generic excitations in one-dimensional spin systems, they can be massless or massive. The quantitative description of massive spinons poses a considerable challenge in spite of various variational approaches. We show that a representation in terms of hopping and Bogoliubov spinon processes, which we call "spinon-pair" operators, and their combination is possible. We refer to such a representation as second quantized form. Neglecting terms which change the number of spinons yields the variational results. Treating the bilinear and quartic terms by continuous unitary transformations leads to considerably improved results. Thus, we provide the proof-of-principle that systems displaying massive spinons as elementary excitations can be treated in second quantization based on spinon-pair representation.


title: Search for long-lived scalar particles in $B^+ \to K^+ \chi (\mu^+\mu^-)$ decays

abstract: A search for a long-lived scalar particle $\chi$ is performed, looking for the decay ${B^+ \to K^+ \chi}$ with ${\chi \to \mu^+\mu^-}$ in $pp$ collision data corresponding to an integrated luminosity of $3\, {\rm fb}^{-1}$, collected by the LHCb experiment at centre-of-mass energies of $\sqrt{s}=7$ and 8$\,$TeV. This new scalar particle, predicted by Hidden Sector models, is assumed to have a narrow width. The signal would manifest itself as an excess in the dimuon invariant mass distribution over the Standard Model background. No significant excess is observed in the accessible ranges of mass ${250 < m(\chi) < 4700\,\rm MeV/c^2}$ and lifetime ${0.1 < \tau(\chi) < 1000\,\rm ps}$. Upper limits on the branching fraction $\mathcal{B}(B^+ \to K^+ \chi (\mu^+\mu^-))$ at 95% confidence level are set as a function of $m(\chi)$ and $\tau(\chi)$, varying between $2\times10^{-10}$ and $10^{-7}$. These are the most stringent limits to date. The limits are interpreted in the context of a model with a light inflaton particle.


title: Planck-scale mass equidistribution of toral Laplace eigenfunctions

abstract: We study the small scale distribution of the $L^2$-mass of eigenfunctions of the Laplacian on the the two-dimensional flat torus. Given an orthonormal basis of eigenfunctions, Lester and Rudnick showed the existence of a density one subsequence whose $L^2$-mass equidistributes more-or-less down to the Planck scale. We give a more precise version of their result showing equidistribution holds down to a small power of log above Planck scale, and also showing that the $L^2$-mass fails to equidistribute at a slightly smaller power of log above the Planck scale. This article rests on a number of results about the proximity of lattice points on circles, much of it based on foundational work of Javier Cilleruelo.


title: Quantifying the degree of average contraction of Collatz orbits

abstract: We here elaborate on a quantitative argument to support the validity of the Collatz conjecture, also known as the (3x + 1) or Syracuse conjecture. The analysis is structured as follows. First, three distinct fixed points are found for the third iterate of the Collatz map, which hence organise in a period 3 orbit of the original map. These are 1, 2 and 4, the elements which define the unique attracting cycle, as hypothesised by Collatz. To carry out the calculation we write the positive integers in modulo 8 (mod8 ), obtain a closed analytical form for the associated map and determine the transitions that yield contracting or expanding iterates in the original, infinite-dimensional, space of positive integers. Then, we consider a Markov chain which runs on the reduced space of mod8 congruence classes of integers. The transition probabilities of the Markov chain are computed from the deterministic map, by employing a measure that is invariant for the map itself. Working in this setting, we demonstrate that the stationary distribution sampled by the stochastic system induces a contracting behaviour for the orbits of the deterministic map on the original space of the positive integers. Sampling the equilibrium distribution on the congruence classes mod8^m for any m, which amounts to arbitrarily reducing the degree of imposed coarse graining, returns an identical conclusion.


title: The Heavy Photon Search beamline and its performance

abstract: The Heavy Photon Search (HPS) is an experiment to search for a hidden sector photon, aka a heavy photon or dark photon, in fixed target electroproduction at the Thomas Jefferson National Accelerator Facility (JLab). The HPS experiment searches for the e$^+$e$^-$ decay of the heavy photon with bump hunt and detached vertex strategies using a compact, large acceptance forward spectrometer, consisting of a silicon microstrip detector (SVT) for tracking and vertexing, and a PbWO$_4$ electromagnetic calorimeter for energy measurement and fast triggering. To achieve large acceptance and good vertexing resolution, the first layer of silicon detectors is placed just 10 cm downstream of the target with the sensor edges only 500 $\mu$m above and below the beam. Placing the SVT in such close proximity to the beam puts stringent requirements on the beam profile and beam position stability. As part of an approved engineering run, HPS took data in 2015 and 2016 at 1.05 GeV and 2.3 GeV beam energies, respectively. This paper describes the beam line and its performance during that data taking.


title: Coverings of torus knots in $S^2\times S^1$ and universals

abstract: Let $t_{\alpha,\beta}\subset S^2\times S^1$ be an ordinary fiber of a Seifert fibering of $S^2\times S^1$ with two exceptional fibers of order $\alpha$. We show that any Seifert manifold with Euler number zero is a branched covering of $S^2\times S^1$ with branching $t_{\alpha,\beta}$ if $\alpha\geq3$. We compute the Seifert invariants of the Abelian covers of $S^2\times S^1$ branched along a $t_{\alpha,\beta}$. We also show that $t_{2,1}$, a non-trivial torus knot in $S^2\times S^1$, is not universal.


title: Logic-based Clustering and Learning for Time-Series Data

abstract: To effectively analyze and design cyberphysical systems (CPS), designers today have to combat the data deluge problem, i.e., the burden of processing intractably large amounts of data produced by complex models and experiments. In this work, we utilize monotonic Parametric Signal Temporal Logic (PSTL) to design features for unsupervised classification of time series data. This enables using off-the-shelf machine learning tools to automatically cluster similar traces with respect to a given PSTL formula. We demonstrate how this technique produces interpretable formulas that are amenable to analysis and understanding using a few representative examples. We illustrate this with case studies related to automotive engine testing, highway traffic analysis, and auto-grading massively open online courses.


title: Structure Preserving H-infinity Optimal PI Control

abstract: A multi-variable PI (proportional integrating) controller is proved to be optimal for an important class of control problems where performance is specified in terms of frequency weighted H-infinity norms. The problem class includes networked systems with a subsystem in each node and control action along each edge. For such systems, the optimal PI controller is decentralized in the sense that control action along a given network edge is entirely determined by states at nodes connected by that edge.


title: Relativistic Zitterbewegung in non-Hermitian photonic waveguide systems

abstract: Zitterbewegung (ZB) is a phenomenon in relativistic quantum systems where the electron wave packet exhibits a trembling or oscillating behavior during its motion, caused by its interaction or coupling with the negative energy state. To directly observe ZB in electronic systems is difficult, due to the challenges associated with the atomic scale wavelength of the electron. Photonic systems offer an alternative paradigm. We exploit the concept of pseudo parity-time (pseudo PT ) symmetry to study ZB in non-Hermitian quantum systems implemented as an experimentally feasible optical waveguide array. In particular, the non-Hermitian Hamiltonian is realized through evanescent coupling among the waveguides to form a one-dimensional lattice with periodic modulations in gain and loss along the guiding direction. As the modulation frequency is changed, we obtain a number of phenomena including periodically suppressed ZB trembling, spatial energy localization, and Hermitian-like ZB oscillations. We calculate phase diagrams indicating the emergence of different types of dynamical behaviors of the relativistic non-Hermitian quantum system in an experimentally justified parameter space. We provide numerical results and a physical analysis to explain the distinct dynamical behaviors revealed by the phase diagrams. Our findings provide a deeper understanding of both the relativistic ZB phenomenon and non-Hermitian pseudo-PT systems, with potential applications in controlling/harnessing light propagation in waveguide-based optical systems.


title: Quantum noise generated by local random Hamiltonians

abstract: We investigate the impact of a local random unitary noise on multipartite quantum states of an arbitrary dimension. We follow the dynamical approach, in which the single-particle unitaries are generated by local random Hamiltonians. Assuming short evolution time we derive a lower bound on the fidelity between an initial and the final state transformed by this type of noise. This result is based on averaging the Tamm-Mandelstam bound and holds for a wide class of distributions of random Hamiltonians fulfilling specific symmetry conditions. It is showed how the sensitivity of a given pure quantum state to the discussed type of noise depends on the properties of single-particle and bipartite reduced states.


title: Adaptive selection of sampling points for Uncertainty Quantification

abstract: We present a simple and robust strategy for the selection of sampling points in Uncertainty Quantification. The goal is to achieve the fastest possible convergence in the cumulative distribution function of a stochastic output of interest. We assume that the output of interest is the outcome of a computationally expensive nonlinear mapping of an input random variable, whose probability density function is known. We use a radial function basis to construct an accurate interpolant of the mapping. This strategy enables adding new sampling points one at a time, adaptively. This takes into full account the previous evaluations of the target nonlinear function. We present comparisons with a stochastic collocation method based on the Clenshaw-Curtis quadrature rule, and with an adaptive method based on hierarchical surplus, showing that the new method often results in a large computational saving.


title: Learning from Simulated and Unsupervised Images through Adversarial Training

abstract: With recent progress in graphics, it has become more tractable to train models on synthetic images, potentially avoiding the need for expensive annotations. However, learning from synthetic images may not achieve the desired performance due to a gap between synthetic and real image distributions. To reduce this gap, we propose Simulated+Unsupervised (S+U) learning, where the task is to learn a model to improve the realism of a simulator's output using unlabeled real data, while preserving the annotation information from the simulator. We develop a method for S+U learning that uses an adversarial network similar to Generative Adversarial Networks (GANs), but with synthetic images as inputs instead of random vectors. We make several key modifications to the standard GAN algorithm to preserve annotations, avoid artifacts, and stabilize training: (i) a 'self-regularization' term, (ii) a local adversarial loss, and (iii) updating the discriminator using a history of refined images. We show that this enables generation of highly realistic images, which we demonstrate both qualitatively and with a user study. We quantitatively evaluate the generated images by training models for gaze estimation and hand pose estimation. We show a significant improvement over using synthetic images, and achieve state-of-the-art results on the MPIIGaze dataset without any labeled real data.


title: Superscattering of pseudospin-1 wave in photonic lattice

abstract: We uncover a superscattering behavior of pseudospin-1 wave from weak scatterers in the subwavelength regime where the scatterer size is much smaller than wavelength. The phenomenon manifests itself as unusually strong scattering characterized by extraordinarily large values of the cross section even for arbitrarily weak scatterer strength. We establish analytically and numerically that the physical origin of superscattering is revival resonances, for which the conventional Born theory breaks down. The phenomenon can be experimentally tested using synthetic photonic systems.


title: The Rearrangement Number

abstract: How many permutations of the natural numbers are needed so that every conditionally convergent series of real numbers can be rearranged to no longer converge to the same sum? We define the \emph{rearrangement number}, a new cardinal characteristic of the continuum, as the answer to this question. We compare the rearrangement number with several natural variants, for example one obtained by requiring the rearranged series to still converge but to a new, finite limit. We also compare the rearrangement number with several well-studied cardinal characteristics of the continuum. We present some new forcing constructions designed to add permutations that rearrange series from the ground model in particular ways, thereby obtaining consistency results going beyond those that follow from comparisons with familiar cardinal characteristics. Finally, we deal briefly with some variants concerning rearrangements by a special sort of permutation and with rearranging some divergent series to become (conditionally) convergent.


title: AR Scorpii is a new white dwarf in the ejector state

abstract: Marsh et al. (2016) have recently reported the discovery of a radio-pulsing white dwarf in the cataclysmic variable AR Sco. The period of pulsations which are also seen in the optical and UV is about 117 seconds. High intensity of pulsing radiation and non-thermal character of its spectrum leave little room for doubt that the white dwarf in AR Sco operates as a spin-powered pulsar and, therefore, is in the ejector state. We show that this system is very much resembling a well-known object AE Aqr. In both systems the compact components are spin-powered and have relatively strong surface magnetic field of order of 100-500 MG. They originated due to accretion spin-up in the previous epoch during which the magnetic field of the white dwarf had substantially evolved being initially buried by the accreted matter and recovered to its initial value after the spin-up phase had ended.


title: Hidden symmetries of deformed oscillators

abstract: We associate with each simple Lie algebra a system of second-order differential equations invariant under a non-compact real form of the corresponding Lie group. In the limit of a contraction to a Schr\"odinger algebra, these equations reduce to a system of ordinary harmonic oscillators. We provide two clarifying examples of such deformed oscillators: one system invariant under SO(2,3) transformations, and another system featuring $G_{2(2)}$ symmetry. The construction of invariant actions requires adding semi-dynamical degrees of freedom; we illustrate the algorithm with the two examples mentioned.


title: Understanding Image and Text Simultaneously: a Dual Vision-Language Machine Comprehension Task

abstract: We introduce a new multi-modal task for computer systems, posed as a combined vision-language comprehension challenge: identifying the most suitable text describing a scene, given several similar options. Accomplishing the task entails demonstrating comprehension beyond just recognizing "keywords" (or key-phrases) and their corresponding visual concepts. Instead, it requires an alignment between the representations of the two modalities that achieves a visually-grounded "understanding" of various linguistic elements and their dependencies. This new task also admits an easy-to-compute and well-studied metric: the accuracy in detecting the true target among the decoys. The paper makes several contributions: an effective and extensible mechanism for generating decoys from (human-created) image captions; an instance of applying this mechanism, yielding a large-scale machine comprehension dataset (based on the COCO images and captions) that we make publicly available; human evaluation results on this dataset, informing a performance upper-bound; and several baseline and competitive learning approaches that illustrate the utility of the proposed task and dataset in advancing both image and language comprehension. We also show that, in a multi-task learning setting, the performance on the proposed task is positively correlated with the end-to-end task of image captioning.


title: The upper bound on the lowest mass halo

abstract: We explore the connection between galaxies and dark matter halos in the Milky Way (MW) and quantify the implications on properties of the dark matter particle and the phenomenology of low-mass galaxy formation. This is done through a probabilistic comparison of the luminosity function of MW dwarf satellite galaxies to models based on two suites of zoom-in simulations. One suite is dark-matter-only while the other includes a disk component, therefore we can quantify the effect of the MW's baryonic disk on our results. We apply numerous Stellar-Mass-Halo-Mass (SMHM) relations allowing for multiple complexities: scatter, a characteristic break scale, and subhalos which host no galaxy. In contrast to previous works we push the model/data comparison to the faintest dwarfs by modeling observational incompleteness, allowing us to draw three new conclusions. Firstly, we constrain the SMHM relation for $10^2<M_*/M_\odot<10^8$ galaxies, allowing us to bound the peak halo mass of the faintest MW satellite to $M_\mathrm{vir}>2.4\times10^8M_\odot$ ($1\sigma$). Secondly, by translating to a Warm Dark Matter (WDM) cosmology, we bound the thermal relic mass $m_\mathrm{WDM}>2.9$ keV at 95\% confidence, on a par with recent constraints from the Lyman-$\alpha$ forest. Lastly, we find that the observed number of ultra-faint MW dwarfs is in tension with the theoretical prediction that reionisation prevents galaxy formation in almost all $10^8M_\odot$ halos. This can be tested with the next generation of deep imaging surveys. To this end, we predict the likely number of detectable satellite galaxies in the Subaru/HSC survey and the LSST. Confronting these predictions with future observations will be amongst our strongest tests of WDM and the effect reionisation on low-mass systems.


title: Influence of Defects on Photoconductivity and Photocatalytic Activity of Nitrogen-Doped Titania

abstract: Samples of nitrogen-doped titanium dioxide (anatase, $0.2 \leq N \leq 1.0~wt\%$) prepared by the sol-gel method were investigated using X-band EPR spectroscopy, photoconductivity and photocatalysis measurements. $N\bullet$ and $NO\bullet$ paramagnetic defects in N-TiO2 have been observed; their concentrations and spin-Hamiltonian parameters were calculated. An increase both in the rate of the generation of free charge carriers and in the rate of photocatalysis was found in N-TiO2 in contrast with non-doped titania under visible light. The correlation of the density of the observed radicals with the photoconductivity and photocatalysis data is discussed.


title: Speciation: Goldschmidt's Heresy, Once Supported by Gould and Dawkins, is Again Reinstated

abstract: The view that the initiation of branching into two sympatric species may not require natural selection emerged in Victorian times (Fleeming Jenkin, George Romanes, William Bateson). In the 1980s paleontologist Steven Jay Gould gave a theoretical underpinning of this non-genic "chromosomal" view, thus reinstating Richard Goldschmidt's "heresy" of the 1930s. From modelling studies with computer-generated "biomorphs," zoologist Richard Dawkins also affirmed Goldschmidt, proclaiming the "evolution of evolvability." However, in the 1990s, while Gould and Dawkins were recanting, bioinformatic, biochemical and cytological studies were providing a deeper underpinning. In 2001 this came under attack from leaders in the field who favored Dawkins' genic emphasis. Now, with growing evidence from multiple sources, we can reinstate again Goldschmidt's view and clarify its nineteenth century roots.


title: SampleRNN: An Unconditional End-to-End Neural Audio Generation Model

abstract: In this paper we propose a novel model for unconditional audio generation based on generating one audio sample at a time. We show that our model, which profits from combining memory-less modules, namely autoregressive multilayer perceptrons, and stateful recurrent neural networks in a hierarchical structure is able to capture underlying sources of variations in the temporal sequences over very long time spans, on three datasets of different nature. Human evaluation on the generated samples indicate that our model is preferred over competing models. We also show how each component of the model contributes to the exhibited performance.


title: Convergence Rates for Greedy Kaczmarz Algorithms, and Faster Randomized Kaczmarz Rules Using the Orthogonality Graph

abstract: The Kaczmarz method is an iterative algorithm for solving systems of linear equalities and inequalities, that iteratively projects onto these constraints. Recently, Strohmer and Vershynin [J. Fourier Anal. Appl., 15(2):262-278, 2009] gave a non-asymptotic convergence rate analysis for this algorithm, spurring numerous extensions and generalizations of the Kaczmarz method. Rather than the randomized selection rule analyzed in that work, in this paper we instead discuss greedy and approximate greedy selection rules. We show that in some applications the computational costs of greedy and random selection are comparable, and that in many cases greedy selection rules give faster convergence rates than random selection rules. Further, we give the first multi-step analysis of Kaczmarz methods for a particular greedy rule, and propose a provably-faster randomized selection rule for matrices with many pairwise-orthogonal rows.


title: Poisson measure as a spectral measure of a family of commuting selfadjoint operators, connected with some moment problem

abstract: It is proved that the Poisson measure is a spectral measure of some family of commuting selfadjoint operators acting on a space constructed from some generalization of the moment problem.


title: Detecting and characterizing high frequency oscillations in epilepsy - A case study of big data analysis

abstract: We develop a framework to uncover and analyze dynamical anomalies from massive, nonlinear and non-stationary time series data. The framework consists of three steps: preprocessing of massive data sets to eliminate erroneous data segments, application of the empirical mode decomposition and Hilbert transform paradigm to obtain the fundamental components embedded in the time series at distinct time scales, and statistical/scaling analysis of the components. As a case study, we apply our framework to detecting and characterizing high frequency oscillations (HFOs) from a big database of rat EEG recordings. We find a striking phenomenon: HFOs exhibit on-off intermittency that can be quantified by algebraic scaling laws. Our framework can be generalized to big data-related problems in other fields such as large-scale sensor data and seismic data analysis.


title: Atom: Horizontally Scaling Strong Anonymity

abstract: Atom is an anonymous messaging system that protects against traffic-analysis attacks. Unlike many prior systems, each Atom server touches only a small fraction of the total messages routed through the network. As a result, the system's capacity scales near-linearly with the number of servers. At the same time, each Atom user benefits from "best possible" anonymity: a user is anonymous among all honest users of the system, against an active adversary who controls the entire network, a portion of the system's servers, and any number of malicious users. The architectural ideas behind Atom have been known in theory, but putting them into practice requires new techniques for (1) avoiding the reliance on heavy general-purpose multi-party computation protocols, (2) defeating active attacks by malicious servers at minimal performance cost, and (3) handling server failure and churn. Atom is most suitable for sending a large number of short messages, as in a microblogging application or a high-security communication bootstrapping ("dialing") for private messaging systems. We show that, on a heterogeneous network of 1,024 servers, Atom can transit a million Tweet-length messages in 28 minutes. This is over 23x faster than prior systems with similar privacy guarantees.


title: New and refined bounds for expected maxima of fractional Brownian motion

abstract: For the fractional Brownian motion $B^H$ with the Hurst parameter value $H$ in (0,1/2), we derive new upper and lower bounds for the difference between the expectations of the maximum of $B^H$ over [0,1] and the maximum of $B^H$ over the discrete set of values $ in^{-1},$ $i=1,\ldots, n.$ We use these results to improve our earlier upper bounds for the expectation of the maximum of $B^H$ over $[0,1]$ and derive new upper bounds for Pickands' constant.


title: "What is Relevant in a Text Document?": An Interpretable Machine Learning Approach

abstract: Text documents can be described by a number of abstract concepts such as semantic category, writing style, or sentiment. Machine learning (ML) models have been trained to automatically map documents to these abstract concepts, allowing to annotate very large text collections, more than could be processed by a human in a lifetime. Besides predicting the text's category very accurately, it is also highly desirable to understand how and why the categorization process takes place. In this paper, we demonstrate that such understanding can be achieved by tracing the classification decision back to individual words using layer-wise relevance propagation (LRP), a recently developed technique for explaining predictions of complex non-linear classifiers. We train two word-based ML models, a convolutional neural network (CNN) and a bag-of-words SVM classifier, on a topic categorization task and adapt the LRP method to decompose the predictions of these models onto words. Resulting scores indicate how much individual words contribute to the overall classification decision. This enables one to distill relevant information from text documents without an explicit semantic information extraction step. We further use the word-wise relevance scores for generating novel vector-based document representations which capture semantic information. Based on these document vectors, we introduce a measure of model explanatory power and show that, although the SVM and CNN models perform similarly in terms of classification accuracy, the latter exhibits a higher level of explainability which makes it more comprehensible for humans and potentially more useful for other applications.


title: Linear-time logics -- a coalgebraic perspective

abstract: We describe a general approach to deriving linear-time logics for a wide variety of state-based, quantitative systems, by modelling the latter as coalgebras whose type incorporates both branching and linear behaviour. Concretely, we define logics whose syntax is determined by the type of linear behaviour, and whose domain of truth values is determined by the type of branching behaviour, and we provide two semantics for them: a step-wise semantics akin to that of standard coalgebraic logics, and a path-based semantics akin to that of standard linear-time logics. The former semantics is useful for model checking, whereas the latter is the more natural semantics, as it measures the extent with which qualitative properties hold along computation paths from a given state. Our main result is the equivalence of the two semantics. We also provide a semantic characterisation of a notion of logical distance induced by these logics. Instances of our logics support reasoning about the possibility, likelihood or minimal cost of exhibiting a given linear-time property.


title: Chemistry of TMC-1 with multiply deuterated species and spin chemistry of H2, H2+, H3+ and their isotopologues

abstract: Deuterated species are unique and powerful tools in astronomy since they can probe the physical conditions, chemistry, and ionization level of various astrophysical media. Recent observations of several deuterated species along with some of their spin isomeric forms have rekindled the interest for more accurate studies on deuterium fractionation. This paper presents the first publicly available chemical network of multiply deuterated species along with spin chemistry implemented on the latest state-of-the-art gas-grain chemical code `NAUTILUS'. D/H ratios for all deuterated species observed at different positions of TMC-1 are compared with the results of our model, which considers multiply deuterated species along with the spin chemistry of light hydrogen bearing species H2, H2+, H3+ and their isotopologues. We also show the differences in the modeled abundances of non-deuterated species after the inclusion of deuteration and spin chemistry in the model. Finally, we present a list of potentially observable deuterated species in TMC-1 awaiting detection.


title: A State Space Approach for Piecewise-Linear Recurrent Neural Networks for Reconstructing Nonlinear Dynamics from Neural Measurements

abstract: The computational properties of neural systems are often thought to be implemented in terms of their network dynamics. Hence, recovering the system dynamics from experimentally observed neuronal time series, like multiple single-unit (MSU) recordings or neuroimaging data, is an important step toward understanding its computations. Ideally, one would not only seek a state space representation of the dynamics, but would wish to have access to its governing equations for in-depth analysis. Recurrent neural networks (RNNs) are a computationally powerful and dynamically universal formal framework which has been extensively studied from both the computational and the dynamical systems perspective. Here we develop a semi-analytical maximum-likelihood estimation scheme for piecewise-linear RNNs (PLRNNs) within the statistical framework of state space models, which accounts for noise in both the underlying latent dynamics and the observation process. The Expectation-Maximization algorithm is used to infer the latent state distribution, through a global Laplace approximation, and the PLRNN parameters iteratively. After validating the procedure on toy examples, the approach is applied to MSU recordings from the rodent anterior cingulate cortex obtained during performance of a classical working memory task, delayed alternation. A model with 5 states turned out to be sufficient to capture the essential computational dynamics underlying task performance, including stimulus-selective delay activity. The estimated models were rarely multi-stable, but rather were tuned to exhibit slow dynamics in the vicinity of a bifurcation point. In summary, the present work advances a semi-analytical (thus reasonably fast) maximum-likelihood estimation framework for PLRNNs that may enable to recover the relevant dynamics underlying observed neuronal time series, and directly link them to computational properties.


title: Curvature of Metrics on Semple Jet bundles

abstract: We discuss the Morse estimates for the curvature of several metrics on Semple weighted projective bundle over a projective variety. Following Demailly works on holomorphic Morse inequalities we show an analogue of his results along the Green-Griffiths conjecture for invariant jets.


title: BSEPACK User's Guide

abstract: This is the user manual for the software package BSEPACK (Bethe--Salpeter Eigenvalue Solver Package).


title: Milne's Equation revisited: Exact Asymptotic Solutions

abstract: We present novel approaches for solving Milne's equation, which was introduced in 1930 as an efficient numerical scheme for the Schr\"odinger equation. Milne's equation appears in a wide class of physical problems, ranging from astrophysics and cosmology, to quantum mechanics and quantum optics. We show how a third order linear differential equation is equivalent to Milne's non-linear equation, and can be used to accurately calculate Milne's amplitude and phase functions. We also introduce optimization schemes to achieve a convenient, fast, and accurate computation of wave functions using an economical parametrization. These new optimization procedures answer the long standing question of finding non-oscillatory solutions of Milne's equation. We apply them to long-range potentials and find numerically exact asymptotic solutions.


title: Automatic Interpretation of Unordered Point Cloud Data for UAV Navigation in Construction

abstract: The objective of this work is to develop a data processing system that can automatically generate waypoints for navigation of an unmanned aerial vehicle (UAV) to inspect surfaces of structures like buildings and bridges. The input includes data recorded by two 2D laser scanners, orthogonally mounted on the UAV, and an inertial measurement unit (IMU). To achieve the goal, algorithms are developed to process the data collected. They are separated into three major groups: (i) the data registration and filtering to generate a 3D model of the structure and control the density of point clouds for data completeness enhancement; (ii) the surface and obstacle detection to assist the UAV in monitoring tasks; and (iii) the waypoint generation to set the flight path. Experiments on different data sets show that the developed system is able to reconstruct a 3D point cloud of the structure, extract its surfaces and objects, and generate waypoints for the UAV to accomplish inspection tasks.


title: On the Generalized Wavelet-Galerkin Method

abstract: In the frame of the traditional wavelet-Galerkin method based on the compactly supported wavelets, it is important to calculate the so-called connection coefficients that are some integrals whose integrands involve products of wavelets, their derivatives as well as some known coefficients in considered differential equations. However, even for linear differential equations with non-constant coefficient, the computation of connect coefficients becomes rather time-consuming and often even impossible. In this paper, we propose a generalized wavelet-Galerkin method based on the compactly supported wavelets, which is computationally very efficient even for differential equations with non-constant coefficients, no matter linear or nonlinear problems. Some related mathematical theorems are proved, based on which the basic ideas of the generalized wavelet-Galerkin method are described in details. In addition, some examples are used to illustrate its validity and high efficiency. A nonlinear example shows that the generalized wavelet-Galerkin method is not only valid to solve nonlinear problems, but also possesses the ability to find new solutions of multi-solution problems. This method can be widely applied to various types of both linear and nonlinear differential equations in science and engineering.


title: Grain boundary-mediated nanopores in molybdenum disulfide grown by chemical vapor deposition

abstract: Molybdenum disulfide (MoS2) is a particularly interesting member of the family of two-dimensional (2D) materials due to its semiconducting and tunable electronic properties. Currently, the most reliable method for obtaining high-quality industrial scale amounts of 2D materials is chemical vapor deposition (CVD), which results in polycrystalline samples. As grain boundaries (GBs) are intrinsic defect lines within CVD-grown 2D materials, their atomic structure is of paramount importance. Here, through atomic-scale analysis of micrometer-long GBs, we show that covalently bound boundaries in 2D MoS2 tend to be decorated by nanopores. Such boundaries occur when differently oriented MoS2 grains merge during growth, whereas the overlap of grains leads to boundaries with bilayer areas. Our results suggest that the nanopore formation is related to stress release in areas with a high concentration of dislocation cores at the grain boundaries, and that the interlayer interaction leads to intrinsic rippling at the overlap regions. This provides insights for the controlled fabrication of large-scale MoS 2 samples with desired structural properties for applications.


title: Active rotational and translational microrheology beyond the linear spring regime

abstract: Active particle tracking microrheometers have the potential to perform accurate broad-band measurements of viscoelasticity within microscopic systems. Generally, their largest possible precision is limited by Brownian motion and low frequency changes to the system. The signal to noise ratio is usually improved by increasing the size of the driven motion compared to the Brownian as well as averaging over repeated measurements. New theory is presented here which gives the complex shear modulus when the motion of a spherical particle is driven by non-linear forces. In some scenarios error can be further reduced by applying a variable transformation which linearises the equation of motion. This allows normalisation which eliminates low frequency drift in the particle's equilibrium position. Using this method will easily increase the signal strength enough to significantly reduce the measurement time for the same error. Thus the method is more conducive to measuring viscoelasticity in slowly changing microscopic systems, such as a living cell.


title: Simultaneous Partial Inverses and Decoding Interleaved Reed-Solomon Codes

abstract: The paper introduces the simultaneous partial-inverse problem (SPI) for polynomials and develops its application to decoding interleaved Reed--Solomon codes beyond half the minimum distance. While closely related both to standard key equations and to well-known Pad\'e approximation problems, the SPI problem stands out in several respects. First, the SPI problem has a unique solution (up to a scale factor), which satisfies a natural degree bound. Second, the SPI problem can be transformed (monomialized) into an equivalent SPI problem where all moduli are monomials. Third, the SPI problem can be solved by an efficient algorithm of the Berlekamp--Massey type. Fourth, decoding interleaved Reed--Solomon codes (or subfield-evaluation codes) beyond half the minimum distance can be analyzed in terms of a partial-inverse condition for the error pattern: if that condition is satisfied, then the (true) error locator polynomial is the unique solution of a standard key equation and can be computed in many different ways, including the well-known multi-sequence Berlekamp--Massey algorithm and the SPI algorithm of this paper. Two of the best performance bounds from the literature (the Schmidt--Sidorenko--Bossert bound and the Roth--Vontobel bound) are generalized to hold for the partial-inverse condition and thus to apply to several different decoding algorithms.


title: Testing parity-violating physics from cosmic rotation power reconstruction

abstract: We study the reconstruction of the cosmic rotation power spectrum produced by parity-violating physics, with an eye to ongoing and near future cosmic microwave background (CMB) experiments such as BICEP Array, CMBS4, LiteBIRD and Simons Observatory. In addition to the inflationary gravitational waves and gravitational lensing, measurements of other various effects on CMB polarization open new window into the early Universe. One of these is anisotropies of the cosmic polarization rotation which probes the Chern-Simons term generally predicted by string theory. The anisotropies of the cosmic rotation are also generated by the primordial magnetism and in the Standard Model extention framework. The cosmic rotation anisotropies can be reconstructed as quadratic in CMB anisotropies. However, the power of the reconstructed cosmic rotation is a CMB four-point correlation and is not directly related to the cosmic-rotation power spectrum. Understanding all contributions in the four-point correlation is required to extract the cosmic rotation signal. Assuming a scale-invariant rotation spectrum motivated by the inflationary cosmic-rotation models, we employ simulation to quantify each contribution to the four-point correlation and find that (1) a secondary contraction of the trispectrum increases the total signal-to-noise, (2) a bias from the lensing-induced trispectrum is significant compared to the statistical errors in, e.g., LiteBIRD and CMBS4-like experiments, (3) the use of a realization-dependent estimator decreases the statistical errors by 10%-20%, depending on experimental specifications, and (4) other higher-order contributions are negligible at least for near future experiments.


title: Boltzmann-Langevin Approach to Pre-equilibrium Correlations in Nuclear Collisions

abstract: Correlations born before the onset of hydrodynamic flow can leave observable traces on the final state particles. Measurement of these correlations can yield important information on the isotropization and thermalization process. Starting from a Boltzmann-like kinetic theory in the presence of dynamic Langevin noise, we derive a new partial differential equation for the two-particle correlation function that respects the microscopic conservation laws. We illustrate how this equation can be used to study the effect of thermalization on long range correlations.


title: Human Action Attribute Learning From Video Data Using Low-Rank Representations

abstract: Representation of human actions as a sequence of human body movements or action attributes enables the development of models for human activity recognition and summarization. We present an extension of the low-rank representation (LRR) model, termed the clustering-aware structure-constrained low-rank representation (CS-LRR) model, for unsupervised learning of human action attributes from video data. Our model is based on the union-of-subspaces (UoS) framework, and integrates spectral clustering into the LRR optimization problem for better subspace clustering results. We lay out an efficient linear alternating direction method to solve the CS-LRR optimization problem. We also introduce a hierarchical subspace clustering approach, termed hierarchical CS-LRR, to learn the attributes without the need for a priori specification of their number. By visualizing and labeling these action attributes, the hierarchical model can be used to semantically summarize long video sequences of human actions at multiple resolutions. A human action or activity can also be uniquely represented as a sequence of transitions from one action attribute to another, which can then be used for human action recognition. We demonstrate the effectiveness of the proposed model for semantic summarization and action recognition through comprehensive experiments on five real-world human action datasets.


title: Interplay of excitation transport and atomic motion in flexible Rydberg aggregates

abstract: Strong resonant dipole-dipole interactions in flexible Rydberg aggregates enable the formation of exciton pulses, the interplay of atomic motion and electronic excitation transfer which feature high fidelity entanglement transport. Here, we demonstrate the coherent splitting of such pulses into two modes, which induce strongly different atomic motion, leading to clear signatures of nonadiabatic effects in atomic density profiles. The mechanism exploits local nonadiabatic effects at a conical intersection, turning them from a decoherence source into an asset. The conical intersection is a consequence of the exciton pulses moving along a linear Rydberg chain and approaching an additional linear, perpendicularly aligned Rydberg chain. The intersection provides a sensitive knob controlling the propagation direction and coherence properties of exciton pulses. We demonstrate that this scenario can be exploited as an exciton switch, controlling direction and coherence properties of the joint pulse on the second of the chains. Initially, we demonstrate the pulse splitting on planar aggregates with atomic motion one-dimensionally constrained and employing isotropic interactions. Subsequently, we confirm the splitting mechanism for a fully realistic scenario in which all spatial restrictions are removed and the full anisotropy of the dipole-dipole interactions is taken into account. Our results enable the experimental observation of non-adiabatic electronic dynamics and entanglement transport with Rydberg atoms. The conical intersection crossings are clearly evident, both in atomic mean position information and excited state spectra of the Rydberg system. This suggests flexible Rydberg aggregates as a test-bench for quantum chemical effects in experiments on much-inflated length scales. The fundamental ideas discussed here have general implications for excitons on a dynamic network.


title: Surveillance and Intervention of Infrastructure-Free Mobile Communications: A New Wireless Security Paradigm

abstract: Conventional wireless security assumes wireless communications are rightful and aims to protect them against malicious eavesdropping and jamming attacks. However, emerging infrastructure-free mobile communication networks are likely to be illegally used (e.g., by criminals or terrorists) but difficult to be monitored, thus imposing new challenges on the public security. To tackle this issue, this article presents a paradigm shift of wireless security to the surveillance and intervention of infrastructure-free suspicious and malicious wireless communications, by exploiting legitimate eavesdropping and jamming jointly. In particular, {\emph{proactive eavesdropping}} (via jamming) is proposed to intercept and decode information from suspicious communication links for the purpose of inferring their intentions and deciding further measures against them. {\emph{Cognitive jamming}} (via eavesdropping) is also proposed so as to disrupt, disable, and even spoof the targeted malicious wireless communications to achieve various intervention tasks.


title: Kinetic simulations of X-B and O-X-B mode conversion and its deterioration at high input power

abstract: Spherical tokamak plasmas are typically overdense and thus inaccessible to externally-injected microwaves in the electron cyclotron range. The electrostatic electron Bernstein wave (EBW), however, provides a method to access the plasma core for heating and diagnostic purposes. Understanding the details of the coupling process to electromagnetic waves is thus important both for the interpretation of microwave diagnostic data and for assessing the feasibility of EBW heating and current drive. While the coupling is reasonably well-understood in the linear regime, nonlinear physics arising from high input power has not been previously quantified. To tackle this problem, we have performed one- and two-dimensional fully kinetic particle-in-cell simulations of the two possible coupling mechanisms, namely X-B and O-X-B mode conversion. We find that the ion dynamics has a profound effect on the field structure in the nonlinear regime, as high amplitude short-scale oscillations of the longitudinal electric field are excited in the region below the high-density cut-off prior to the arrival of the EBW. We identify this effect as the instability of the X wave with respect to resonant scattering into an EBW and a lower-hybrid wave. We calculate the instability rate analytically and find this basic theory to be in reasonable agreement with our simulation results.


title: Prediction and Characterization of Multiple Extremal Paths in Continuously Monitored Qubits

abstract: We examine most-likely paths between initial and final states for diffusive quantum trajectories in continuously monitored pure-state qubits, obtained as extrema of a stochastic path integral. We demonstrate the possibility of "multipaths" in the dynamics of continuously-monitored qubit systems, wherein multiple most-likely paths travel between the same pre- and post-selected states over the same time interval. Most-likely paths are expressed as solutions to a Hamiltonian dynamical system. The onset of multipaths may be determined by analyzing the evolution of a Lagrange manifold in this phase space, and is mathematically analogous to the formation of caustics in ray optics or semiclassical physics. Additionally, we develop methods for finding optimal traversal times between states, or optimal final states given an initial state and evolution time; both give insight into the measurement dynamics of continuously-monitored quantum states. We apply our methods in two systems: a qubit with two non-commuting observables measured simultaneously, and a qubit measured in one observable while subject to Rabi drive. In the two-observable case we find multipaths due to caustics, bounded by a diverging Van-Vleck determinant, and their onset time. We also find multipaths generated by paths with different "winding numbers" around the Bloch sphere in both systems.


title: Optimal Resource Allocation for Cellular Networks with MATLAB Instructions

abstract: This report presents a more detailed description of the algorithm and simulations published in papers [1, 2]. It includes a step by step description of the algorithm and included the corresponding flow chart. In addition, detailed instructions of the MATLAB code used to simulate the proposed allocation algorithm in [1, 2] is presented. The report starts with a brief motivation of resource allocation problem in wireless networks. Then, some of the prior related work on the subject are mentioned. Finally, we provide the details instructions on MATLAB functions used in our algorithm. More rigorous analysis and proofs of the problem and algorithm are present in [1, 2] and further discussions presented in [3, 4].


title: Anatomy of Scholarly Information Behavior Patterns in the Wake of Academic Social Media Platforms

abstract: As more scholarly content is born digital or converted to a digital format, digital libraries are becoming increasingly vital to researchers seeking to leverage scholarly big data for scientific discovery. Although scholarly products are available in abundance-especially in environments created by the advent of social networking services-little is known about international scholarly information needs, information-seeking behavior, or information use. The purpose of this paper is to address these gaps via an in-depth analysis of the information needs and information-seeking behavior of researchers, both students and faculty, at two universities, one in the U.S. and the other in Qatar. Based on this analysis, the study identifies and describes new behavior patterns on the part of researchers as they engage in the information-seeking process. The analysis reveals that the use of academic social networks has notable effects on various scholarly activities. Further, this study identifies differences between students and faculty members in regard to their use of academic social networks, and it identifies differences between researchers according to discipline. Although the researchers who participated in the present study represent a range of disciplinary and cultural backgrounds, the study reports a number of similarities in terms of the researchers' scholarly activities.


title: Validating gravitational-wave detections: The Advanced LIGO hardware injection system

abstract: Hardware injections are simulated gravitational-wave signals added to the Laser Interferometer Gravitational-wave Observatory (LIGO). The detectors' test masses are physically displaced by an actuator in order to simulate the effects of a gravitational wave. The simulated signal initiates a control-system response which mimics that of a true gravitational wave. This provides an end-to-end test of LIGO's ability to observe gravitational waves. The gravitational-wave analyses used to detect and characterize signals are exercised with hardware injections. By looking for discrepancies between the injected and recovered signals, we are able to characterize the performance of analyses and the coupling of instrumental subsystems to the detectors' output channels. This paper describes the hardware injection system and the recovery of injected signals representing binary black hole mergers, a stochastic gravitational wave background, spinning neutron stars, and sine-Gaussians.


title: Infinite Volume of Noncommutative Black Hole Wrapped by Finite Surface

abstract: The volume of a black hole under noncommutative spacetime background is found to be infinite, in contradiction with the surface area of a black hole, or its Bekenstein-Hawking (BH) entropy, which is well-known to be finite. Our result rules out the possibility of interpreting the entropy of a black hole by counting the number of modes wrapped inside its surface if the final evaporation stage can be properly treated. It implies the statistical interpretation for the BH entropy can be independent of the volume, provided spacetime is noncommutative. The effect of radiation back reaction is found to be small and doesn't influence the above conclusion.


title: Spectral algorithms for tensor completion

abstract: In the tensor completion problem, one seeks to estimate a low-rank tensor based on a random sample of revealed entries. In terms of the required sample size, earlier work revealed a large gap between estimation with unbounded computational resources (using, for instance, tensor nuclear norm minimization) and polynomial-time algorithms. Among the latter, the best statistical guarantees have been proved, for third-order tensors, using the sixth level of the sum-of-squares (SOS) semidefinite programming hierarchy (Barak and Moitra, 2014). However, the SOS approach does not scale well to large problem instances. By contrast, spectral methods --- based on unfolding or matricizing the tensor --- are attractive for their low complexity, but have been believed to require a much larger sample size. This paper presents two main contributions. First, we propose a new unfolding-based method, which outperforms naive ones for symmetric $k$-th order tensors of rank $r$. For this result we make a study of singular space estimation for partially revealed matrices of large aspect ratio, which may be of independent interest. For third-order tensors, our algorithm matches the SOS method in terms of sample size (requiring about $rd^{3/2}$ revealed entries), subject to a worse rank condition ($r\ll d^{3/4}$ rather than $r\ll d^{3/2}$). We complement this result with a different spectral algorithm for third-order tensors in the overcomplete ($r\ge d$) regime. Under a random model, this second approach succeeds in estimating tensors of rank $d\le r \ll d^{3/2}$ from about $rd^{3/2}$ revealed entries.


title: Sequential nonparametric tests for a change in distribution: an application to detecting radiological anomalies

abstract: We propose a sequential nonparametric test for detecting a change in distribution, based on windowed Kolmogorov--Smirnov statistics. The approach is simple, robust, highly computationally efficient, easy to calibrate, and requires no parametric assumptions about the underlying null and alternative distributions. We show that both the false-alarm rate and the power of our procedure are amenable to rigorous analysis, and that the method outperforms existing sequential testing procedures in practice. We then apply the method to the problem of detecting radiological anomalies, using data collected from measurements of the background gamma-radiation spectrum on a large university campus. In this context, the proposed method leads to substantial improvements in time-to-detection for the kind of radiological anomalies of interest in law-enforcement and border-security applications.


title: Homotopical properties of the simplicial Maurer-Cartan functor

abstract: We consider the category whose objects are filtered, or complete, $L_\infty$-algebras and whose morphisms are $\infty$-morphisms which respect the filtrations. We then discuss the homotopical properties of the Getzler-Hinich simplicial Maurer-Cartan functor which associates to each filtered $L_\infty$-algebra a Kan simplicial set, or $\infty$-groupoid. In previous work with V. Dolgushev, we showed that this functor sends weak equivalences of filtered $L_\infty$-algebras to weak homotopy equivalences of simplicial sets. Here we sketch a proof of the fact that this functor also sends fibrations to Kan fibrations. To the best of our knowledge, only special cases of this result have previously appeared in the literature. As an application, we show how these facts concerning the simplicial Maurer--Cartan functor provide a simple $\infty$-categorical formulation of the Homotopy Transfer Theorem.


title: Large time asymptotics of solutions to the short-pulse equation

abstract: We consider the long-time behavior of solutions to the short-pulse equation. Using the method of testing by wave packets, we prove small data global existence and modified scattering.


title: Norm inflation for the generalized Boussinesq and Kawahara equations

abstract: We consider ill-posedness of the Cauchy problem for the generalized Boussinesq and Kawahara equations. We prove norm inflation with general initial data, an improvement over the ill-posedness results by Geba et al., Nonlinear Anal. 95 (2014), 404-413 for the generalized Boussinesq equations and by Kato, Adv. Differential Equations 16 (2011), no. 3-4, 257-287 for the Kawahara equation.


title: On the chiral covariant approach to $\rho\rho$ scattering

abstract: We examine in detail a recent work (D.~G\"ulmez, U.-G.~Mei\ss ner and J.~A.~Oller, Eur. Phys. J. C 77:460 (2017)), where improvements to make $\rho\rho$ scattering relativistically covariant are made. The paper has the remarkable conclusion that the $J=2$ state disappears with a potential which is much more attractive than for $J=0$, where a bound state is found. We trace this abnormal conclusion to the fact that an "on-shell" factorization of the potential is done in a region where this potential is singular and develops a large discontinuous and unphysical imaginary part. A method is developed, evaluating the loops with full $\rho$ propagators, and we show that they do not develop singularities and do not have an imaginary part below threshold. With this result for the loops we define an effective potential, which when used with the Bethe-Salpeter equation provides a state with $J=2$ around the energy of the $f_2(1270)$. In addition, the coupling of the state to $\rho\rho$ is evaluated and we find that this coupling and the $T$ matrix around the energy of the bound state are remarkably similar to those obtained with a drastic approximation used previously, in which the $q^2$ terms of the propagators of the exchanged $\rho$ mesons are dropped, once the cut-off in the $\rho\rho$ loop function is tuned to reproduce the bound state at the same energy.


title: Object Shape Approximation & Contour Adaptive Depth Image Coding for Virtual View Synthesis

abstract: A depth image provides partial geometric information of a 3D scene, namely the shapes of physical objects as observed from a particular viewpoint. This information is important when synthesizing images of different virtual camera viewpoints via depth-image-based rendering (DIBR). It has been shown that depth images can be efficiently coded using contour-adaptive codecs that preserve edge sharpness, resulting in visually pleasing DIBR-synthesized images. However, contours are typically losslessly coded as side information (SI), which is expensive if the object shapes are complex. In this paper, we pursue a new paradigm in depth image coding for color-plus-depth representation of a 3D scene: we pro-actively simplify object shapes in a depth and color image pair to reduce depth coding cost, at a penalty of a slight increase in synthesized view distortion. Specifically, we first mathematically derive a distortion upper-bound proxy for 3DSwIM---a quality metric tailored for DIBR-synthesized images. This proxy reduces interdependency among pixel rows in a block to ease optimization. We then approximate object contours via a dynamic programming (DP) algorithm to optimally trade off coding cost of contours using arithmetic edge coding (AEC) with our proposed view synthesis distortion proxy. We modify the depth and color images according to the approximated object contours in an inter-view consistent manner. These are then coded respectively using a contour-adaptive image codec based on graph Fourier transform (GFT) for edge preservation and HEVC intra. Experimental results show that by maintaining sharp but simplified object contours during contour-adaptive coding, for the same visual quality of DIBR-synthesized virtual views, our proposal can reduce depth image coding rate by up to 22% compared to alternative coding strategies such as HEVC intra.


title: Accelerating Lattice QCD Multigrid on GPUs Using Fine-Grained Parallelization

abstract: The past decade has witnessed a dramatic acceleration of lattice quantum chromodynamics calculations in nuclear and particle physics. This has been due to both significant progress in accelerating the iterative linear solvers using multi-grid algorithms, and due to the throughput improvements brought by GPUs. Deploying hierarchical algorithms optimally on GPUs is non-trivial owing to the lack of parallelism on the coarse grids, and as such, these advances have not proved multiplicative. Using the QUDA library, we demonstrate that by exposing all sources of parallelism that the underlying stencil problem possesses, and through appropriate mapping of this parallelism to the GPU architecture, we can achieve high efficiency even for the coarsest of grids. Results are presented for the Wilson-Clover discretization, where we demonstrate up to 10x speedup over present state-of-the-art GPU-accelerated methods on Titan. Finally, we look to the future, and consider the software implications of our findings.


title: Multiple simultaneous Gamma-ray emission sites in 3C 84

abstract: The mis-aligned active galactic nuclei (AGN), 3C\,84, has been observed approximately monthly at four frequencies simultaneously (22\,GHz -- 129\,GHz) on the Korean VLBI Network (KVN) since 2013 as part of the interferometric monitoring of $\gamma$-ray bright AGN (iMOGABA) program. 3C\,84 is known to have to sites of bright radio emission, with the first thought to be near the central super-massive black hole and the second in a slowly moving emission feature, several milliarcseconds south of the core. Analysis of this data suggests that the highly variable Gamma-ray emission originates near the SMBH and the slowly rising $\gamma$-ray emission originates in the slow moving feature.


title: Streaming GPU Singular Value and Dynamic Mode Decompositions

abstract: This work develops a parallelized algorithm to compute the dynamic mode decomposition (DMD) on a graphics processing unit using the streaming method of snapshots singular value decomposition. This allows the algorithm to operate efficiently on streaming data by avoiding redundant inner-products as new data becomes available. In addition, it is possible to leverage the native compressed format of many data streams, such as HD video and computational physics codes that are represented sparsely in the Fourier domain, to massively reduce data transfer from CPU to GPU and to enable sparse matrix multiplications. Taken together, these algorithms facilitate real-time streaming DMD on high-dimensional data streams. We demonstrate the proposed method on numerous high-dimensional data sets ranging from video background modeling to scientific computing applications, where DMD is becoming a mainstay algorithm. The computational framework is developed as an open-source library written in C++ with CUDA, and the algorithms may be generalized to include other DMD advances, such as compressed sensing DMD, multi resolution DMD, or DMD with control. Keywords: Singular value decomposition, dynamic mode decomposition, streaming computations, graphics processing unit, video background modeling, scientific computing.


title: Multiple Periodic Solutions for $\Gamma$-symmetric Newtonian Systems

abstract: The existence of periodic solutions in $\Gamma$-symmetric Newtonian systems $\ddot{x}=-\nabla f(x)$ can be effectively studied by means of the $(\Gamma\times O(2))$-equivariant gradient degree with values in the Euler ring $U(\Gamma\times O(2))$. In this paper, we show that in the case of $\Gamma$ being a finite group, the Euler ring $U(\Gamma\times O(2))$ and the related basic degrees are effectively computable using Euler ring homomorphisms, the Burnside ring $A(\Gamma\times O(2))$ and the reduced $(\Gamma\times O(2))$-degree with no free parameters. We present several examples of Newtonian systems with various symmetries, for which we show existence of multiple periodic solutions. We also provide exact value of the equivariant topological invariant for those problems.


title: On Physical Realizability of Nonlinear Quantum Stochastic Differential Equations

abstract: In this article we study physical realizability for a class of nonlinear quantum stochastic differential equations (QSDEs). Physical realizability is a property in which a QSDE corresponds to the dynamics of an open quantum system. We derive a sufficient and necessary condition for a nonlinear QSDE to be physically realizable.


title: Markov-Nash Equilibria in Mean-Field Games with Discounted Cost

abstract: In this paper, we consider discrete-time dynamic games of the mean-field type with a finite number $N$ of agents subject to an infinite-horizon discounted-cost optimality criterion. The state space of each agent is a locally compact Polish space. At each time, the agents are coupled through the empirical distribution of their states, which affects both the agents' individual costs and their state transition probabilities. We introduce a new solution concept of the Markov-Nash equilibrium, under which a policy is player-by-player optimal in the class of all Markov policies. Under mild assumptions, we demonstrate the existence of a mean-field equilibrium in the infinite-population limit $N \to \infty$, and then show that the policy obtained from the mean-field equilibrium is approximately Markov-Nash when the number of agents $N$ is sufficiently large.


title: MRFI Stream Arbitration: An Efficient Arbitration Scheme for NoC with Emerging Interconnect Technology

abstract: An improved version of stream arbitration based on multiband RF interconnect (MRFI) is proposed. Thanks to the simultaneous multiple channel transmitting/receiving feature of MRFI, dynamic bandwidth allocation is achieved in the proposed arbitration algorithm. With dynamic bandwidth allocation, MRFI based arbitration can guarantee 100% channel bandwidth utilization, which is a significant improvement compared with original RF-I based stream arbitration whose channel bandwidth utilization is only around 30%~50%.


title: Visualizing dispersive features in 2D image via minimum gradient method

abstract: We developed a minimum gradient based method to track ridge features in 2D image plot, which is a typical data representation in many momentum resolved spectroscopy experiments. Through both analytic formulation and numerical simulation, we compare this new method with existing DC (distribution curve) based and higher order derivative based analyses. We find that the new method has good noise resilience and enhanced contrast especially for weak intensity features, meanwhile preserves the quantitative local maxima information from the raw image. An algorithm is proposed to extract 1D ridge dispersion from the 2D image plot, whose quantitative application to angle-resolved photoemission spectroscopy measurements on high temperature superconductors is demonstrated.


title: On incomplete and synchronizing finite sets

abstract: This paper situates itself in the theory of variable length codes and of finite automata where the concepts of completeness and synchronization play a central role. In this theoretical setting, we investigate the problem of finding upper bounds to the minimal length of synchronizing words and incompletable words of a finite language X in terms of the length of the words of X. This problem is related to two well-known conjectures formulated by Cerny and Restivo, respectively. In particular, if Restivo's conjecture is true, our main result provides a quadratic bound for the minimal length of a synchronizing pair of any finite synchronizing complete code with respect to the maximal length of its words.


title: Semi-coherent Detection and Performance Analysis for Ambient Backscatter System

abstract: We study a novel communication mechanism, ambient backscatter, that utilizes radio frequency (RF) signals transmitted from an ambient source as both energy supply and information carrier to enable communications between low-power devices. Different from existing non-coherent schemes, we here design the semi-coherent detection, where channel parameters can be obtained from unknown data symbols and a few pilot symbols. We first derive the optimal detector for the complex Gaussian ambient RF signal from likelihood ratio test and compute the corresponding closed-form bit error rate (BER). To release the requirement for prior knowledge of the ambient RF signal, we next design a suboptimal energy detector with ambient RF signals being either the complex Gaussian or the phase shift keying (PSK). The corresponding detection thresholds, the analytical BER, and the outage probability are also obtained in closed-form. Interestingly, the complex Gaussian source would cause an error floor while the PSK source does not, which brings nontrivial indication of constellation design as opposed to the popular Gaussian-embedded literatures. Simulations are provided to corroborate the theoretical studies.


title: Cosmological Inference from Host-Selected Type Ia Supernova Samples

abstract: We compare two Type Ia supernova (SN Ia) samples that are drawn from a spectroscopically confirmed SN Ia sample: a host-selected sample in which SNe Ia are restricted to those that have a spectroscopic redshift from the host; and a broader, more traditional sample in which the redshift could come from either the SN or the host. The host-selected sample is representative of SN samples that will use the redshift of the host to infer the SN redshift, long after the SN has faded from view. We find that SNe Ia that are selected on the availability of a redshift from the host differ from SNe Ia that are from the broader sample. The former tend to be redder, have narrower light curves, live in more massive hosts, and tend to be at lower redshifts. We find that constraints on the equation of state of dark energy, $w$, and the matter density, $\Omega_M$, remain consistent between these two types of samples. Our results are important for ongoing and future supernova surveys, which unlike previous supernova surveys, will have limited real-time follow-up to spectroscopically classify the SNe they discover. Most of the redshifts in these surveys will come from the hosts.


title: Relative Rips Machine

abstract: We study isometric actions of finitely presented groups on $\mathbb{R}$-trees. In this paper, we develop a relative version of the Rips machine to study $\textit{pairs}$ of such actions. An important example of a $\textit{pair}$ is a group action on an $\mathbb{R}$-tree and a subgroup action on its minimal invariant subtree.


title: Triggering the formation of direct collapse black holes by their congeners

abstract: Direct collapse black holes (DCBHs) are excellent candidates as seeds of supermassive black holes (SMBHs) observed at $z \gsim 6$. The formation of a DCBH requires a strong external radiation field to suppress $\rm H_2$ formation and cooling in a collapsing gas cloud. Such strong field is not easily achieved by first stars or normal star-forming galaxies. Here we investigate a scenario in which the previously-formed DCBH can provide the necessary radiation field for the formation of additional ones. Using one-zone model and the simulated DCBH Spectral Energy Distributions (SEDs) filtered through absorbing gas initially having column density $N_{\rm H}$, we derive the critical field intensity, $J_{\rm LW}^{\rm crit}$, to suppress $\rm H_2$ formation and cooling. For the SED model with $N_{\rm H}=1.3\times10^{25}$ cm$^{-2}$, $8.0\times10^{24}$ cm$^{-2}$ and $5.0\times10^{24}$ cm$^{-2}$, we obtain $J_{\rm LW}^{\rm crit}\approx22$, 35 and 54, all much smaller than the critical field intensity for normal star-forming galaxies ($J_{\rm LW}^{\rm crit}\simgt 1000$). X-ray photons from previously-formed DCBHs build up a high-$z$ X-ray background (XRB) that may boost the $J_{\rm LW}^{\rm crit}$. However, we find that in the three SED models $J_{\rm LW}^{\rm crit}$ only increases to $\approx80$, 170 and 390 respectively even when $\dt{\rho}_\bullet$ reaches the maximum value allowed by the present-day XRB level ($0.22, 0.034, 0.006~M_\odot$yr$^{-1}$Mpc$^{-3}$), still much smaller than the galactic value. Although considering the XRB from first galaxies may further increase $J_{\rm LW}^{\rm crit}$, we conclude that our investigation supports a scenario in which DCBH may be more abundant than predicted by models only including galaxies as external radiation sources.


title: SU(N) Fermi liquid at finite temperature

abstract: We consider the thermodynamic potential $\Omega$ of an N component Fermi gas with a short range interaction obeying SU(N) symmetry. We analyze especially the part of $\Omega$ that depends on the temperature T non-analytically for small T . We examine the temperature range where one can observe this $T^4 ln T$ contribution and discuss how it can be extracted experimentally.


title: The dark side of solar photospheric G-band bright points

abstract: Bright small-scale magnetic elements found mainly in intergranular lanes at the solar surface are named bright points (BPs). They show high contrasts in Fraunhofer G-band observations and are described by nearly vertical slender flux tubes or sheets. A recent comparison between BP observations in the ultraviolet (UV) and visible spectral range recorded with the balloon-borne observatory SUNRISE and state-of-the-art magnetohydrodynamical (MHD) simulations revealed a kiloGauss magnetic field for 98% of the synthetic BPs. Here we address the opposite question, namely which fraction of pixels hosting kiloGauss fields coincides with an enhanced G-band brightness. We carried out 3D radiation MHD simulations for three magnetic activity levels (corresponding to the quiet Sun, weak and strong plage) and performed a full spectral line synthesis in the G-band. Only 7% of the kiloGauss pixels in our quiet-Sun simulation coincide with a brightness lower than the mean quiet-Sun intensity, while 23% of the pixels in the weak-plage simulation and even 49% in the strong-plage simulation are associated with a local darkening. Dark strong-field regions are preferentially found in the cores of larger flux patches that are rare in the quiet Sun, but more common in plage regions, often in the vertices of granulation cells. The significant brightness shortfall in the core of larger flux patches coincide with a slight magnetic field weakening. KiloGauss elements in the quiet Sun are on average brighter than similar features in plage regions. Almost all strong-field pixels display a more or less vertical magnetic field orientation. Hence in the quiet Sun, G-band BPs correspond almost one-to-one with kiloGauss elements. In weak plage the correspondence is still very good, but not perfect.


title: On the genus of the complete tripartite graph $K_{n,n,1}$

abstract: For even $n$ we prove that the genus of the complete tripartite graph $K_{n,n,1}$ is $\lceil (n-1) (n-2)/4 \rceil$. This is the least number of bridges needed to build a complete $n$-way road interchange where changing lanes is not allowed. Both the theoretical result, and the surprising link to modelling road intersections are new.


title: Information-Centric Wireless Networks with Virtualization and D2D Communications

abstract: Wireless network virtualization and information-centric networking (ICN) are two promising technologies for next generation wireless networks. Although some excellent works have focused on these two technologies, device-to-device (D2D) communications have not beeen investigated in information-centric virtualized cellular networks. Meanwhile, content caching in mobile devices has attracted great attentions due to the saved backhaul consumption or reduced transmission latency in D2D-assisted cellular networks. However, when it comes to the multi-operator scenario, the direct content sharing between different operators via D2D communications is typically infeasible. In this article, we propose a novel information-centric virtualized cellular network framework with D2D communications, enabling not only content caching in the air, but also inter-operator content sharing between mobile devices. Moreover, we describe the key components in the proposed framework, and present the interactions among them. In addition, we incorporate and formulate the content caching strategies in resource allocation optimization, to maximize the total utility of mobile virtual network operators (MVNOs) through caching popular contents in mobile devices. Simulations results demonstrate the effectiveness of the proposed framework and scheme with different system parameters.


title: Direct relation between the excess current noise spectral density and the complex conductivity of disordered conductor

abstract: We put forward a fundamental relationship between the normalized spectral density of excess current noise and the imaginary part of complex conductivity.


title: Radio Emission of Air Showers with Energy E$_0$ $\geq$ 10$^{19}$ eV by Yakutsk Array Data

abstract: In this paper, we present results obtained from the measurements of radio emission at frequency of 32 MHz with energy more than 10$^{19}$ eV. Generalized formula that describe lateral distribution and depends on main characteristic of the air showers: energy E$_0$ and depth of maximum X$_{max}$ was derived. The formula has a good agreement with data at average and large distances from shower axis. Employing the ratio of radio emission amplitude at distances 175 m and 725 m we determined the depth of maximum X$_{max}$ for air shower with energy 3.7$\cdot$10$^{19}$ eV, which in our case is equal to X$_{max}$ = 769$\pm$34g$\cdot$cm$^{-2}$.


title: Search for Neutrino Generated Air Shower Candidates with Energy More than 5$\cdot$10$^{18}$ eV and Zenith Angle $\theta$ $\geq$ 50$^\circ$

abstract: Neutrino air showers can be formed in any part of the atmosphere passing a long was in the matter due to its physical properties. In general, air showers produced by neutrinos are highly inclined and formed near the ground level, i.e. young showers. Therefore, one should expect a large number of peaks in the signal of such air showers [1, 2, 3]. The goal of our work is to search for air shower candidates produced by neutrino. For this purposes, we analyzed large amount of data from scintillation detectors with different area and energy thresholds [1, 2]. Preliminary analysis of Yakutsk array data indicated the absence of air shower produced by neutrino, but it does not mean that such air showers does not exist. It is going to need a further analysis of highly inclined showers. In order to do that improved methodology for recording and processing of air showers required.


title: Cross-Color Channel Perceptually Adaptive Quantization for HEVC

abstract: HEVC includes a Coding Unit (CU) level luminance-based perceptual quantization technique known as AdaptiveQP. AdaptiveQP perceptually adjusts the Quantization Parameter (QP) at the CU level based on the spatial activity of raw input video data in a luma Coding Block (CB). In this paper, we propose a novel cross-color channel adaptive quantization scheme which perceptually adjusts the CU level QP according to the spatial activity of raw input video data in the constituent luma and chroma CBs; i.e., the combined spatial activity across all three color channels (the Y, Cb and Cr channels). Our technique is evaluated in HM 16 with 4:4:4, 4:2:2 and 4:2:0 YCbCr JCT-VC test sequences. Both subjective and objective visual quality evaluations are undertaken during which we compare our method with AdaptiveQP. Our technique achieves considerable coding efficiency improvements, with maximum BD-Rate reductions of 15.9% (Y), 13.1% (Cr) and 16.1% (Cb) in addition to a maximum decoding time reduction of 11.0%.


title: A global Torelli theorem for singular symplectic varieties

abstract: We systematically study the moduli theory of symplectic varieties (in the sense of Beauville) which admit a resolution by an irreducible symplectic manifold. In particular, we prove an analog of Verbitsky's global Torelli theorem for the locally trivial deformations of such varieties. Verbitsky's work on ergodic complex structures replaces twistor lines as the essential global input. In so doing we extend many of the local deformation-theoretic results known in the smooth case to such (not-necessarily-projective) symplectic varieties. We deduce a number of applications to the birational geometry of symplectic manifolds, including some results on the classification of birational contractions of $K3^{[n]}$-type varieties.


title: Presence of Turbulent and Ordered Local Structure within ICME shock-sheath and its Contribution in Forbush decrease

abstract: The transient interplanetary disturbances evoke short time cosmic ray flux decrease which is known as Forbush decrease. The traditional model and understanding of Forbush decrease suggest that the substructure of interplanetary counterpart of coronal mass ejection (ICME) independently contributes in cosmic ray flux decrease. These substructures, shock-sheath and magnetic cloud (MC) manifest as classical two-step Forbush decrease. The recent work by Raghav et. al. (2016a) has shown multi-step decreases and recoveries within shock-sheath. However, this can not be explained by ideal shock-sheath barrier model. Further, they suggested that the local structures within the ICME's sub-structure (MC and shock sheath) could explain this deviation of FD profile from the classical FD. Therefore, present study attempts to investigate the cause of multi-step cosmic ray flux decrease and their respective recovery within shock-sheath in detail. 3D-hodogram method has been utilized to get more details of the local structures within the shock-sheath. A 3D-hodogram method unambiguously suggests the formation of small scale local structures within the ICME (shock-sheath and even in MC). Moreover, the method could differentiate the turbulent and ordered Interplanetary Magnetic Field (IMF) regions within the substructures of ICME. The study explicitly suggests that the turbulent and ordered IMF regions within shock-sheath do influence cosmic-ray variations uniquely.


title: A Base Camp for Scaling AI

abstract: Modern statistical machine learning (SML) methods share a major limitation with the early approaches to AI: there is no scalable way to adapt them to new domains. Human learning solves this in part by leveraging a rich, shared, updateable world model. Such scalability requires modularity: updating part of the world model should not impact unrelated parts. We have argued that such modularity will require both "correctability" (so that errors can be corrected without introducing new errors) and "interpretability" (so that we can understand what components need correcting). To achieve this, one could attempt to adapt state of the art SML systems to be interpretable and correctable; or one could see how far the simplest possible interpretable, correctable learning methods can take us, and try to control the limitations of SML methods by applying them only where needed. Here we focus on the latter approach and we investigate two main ideas: "Teacher Assisted Learning", which leverages crowd sourcing to learn language; and "Factored Dialog Learning", which factors the process of application development into roles where the language competencies needed are isolated, enabling non-experts to quickly create new applications. We test these ideas in an "Automated Personal Assistant" (APA) setting, with two scenarios: that of detecting user intent from a user-APA dialog; and that of creating a class of event reminder applications, where a non-expert "teacher" can then create specific apps. For the intent detection task, we use a dataset of a thousand labeled utterances from user dialogs with Cortana, and we show that our approach matches state of the art SML methods, but in addition provides full transparency: the whole (editable) model can be summarized on one human-readable page. For the reminder app task, we ran small user studies to verify the efficacy of the approach.


title: Criticality in Cell Differentiation

abstract: Cell differentiation is an important process in living organisms. Differentiation is mostly based on binary decisions with the progenitor cells choosing between two specific lineages. The differentiation dynamics have both deterministic and stochastic components. Several theoretical studies suggest that cell differentiation is a bifurcation phenomenon, well-known in dynamical systems theory. The bifurcation point has the character of a critical point with the system dynamics exhibiting specific features in its vicinity. These include the critical slowing down, rising variance and lag-1 autocorrelation function, strong correlations between the fluctuations of key variables and non-Gaussianity in the distribution of fluctuations. Recent experimental studies provide considerable support to the idea of criticality in cell differentiation and in other biological processes like the development of the fruit fly embryo. In this Review, an elementary introduction is given to the concept of criticality in cell differentiation. The correspondence between the signatures of criticality and experimental observations on blood cell differentiation in mice is further highlighted.


title: On component groups of Jacobians of quaternionic modular curves

abstract: We use a combinatorial result relating the discriminant of the cycle pairing on a weighted finite graph to the eigenvalues of its Laplacian to deduce a formula for the orders of component groups of Jacobians of modular curves arising from quaternion algebras over $\mathbb{F}_q(T)$ or $\mathbb{Q}$. Our formula over $\mathbb{Q}$ recovers a result of Jordan and Livn\'e.


title: DARN: a Deep Adversial Residual Network for Intrinsic Image Decomposition

abstract: We present a new deep supervised learning method for intrinsic decomposition of a single image into its albedo and shading components. Our contributions are based on a new fully convolutional neural network that estimates absolute albedo and shading jointly. Our solution relies on a single end-to-end deep sequence of residual blocks and a perceptually-motivated metric formed by two adversarially trained discriminators. As opposed to classical intrinsic image decomposition work, it is fully data-driven, hence does not require any physical priors like shading smoothness or albedo sparsity, nor does it rely on geometric information such as depth. Compared to recent deep learning techniques, we simplify the architecture, making it easier to build and train, and constrain it to generate a valid and reversible decomposition. We rediscuss and augment the set of quantitative metrics so as to account for the more challenging recovery of non scale-invariant quantities. We train and demonstrate our architecture on the publicly available MPI Sintel dataset and its intrinsic image decomposition, show attenuated overfitting issues and discuss generalizability to other data. Results show that our work outperforms the state of the art deep algorithms both on the qualitative and quantitative aspect.


title: Spaces of Sums of Powers and Real Rank Boundaries

abstract: We investigate properties of Waring decompositions of real homogeneous forms. We study the moduli of real decompositions, so-called Space of Sums of Powers, naturally included in the Variety of Sums of Powers. Explicit results are obtained for quaternary quadrics, relating the algebraic boundary of ${\rm SSP}$ to various loci in the Hilbert scheme of four points in $\mathbb{P}^3$. Further, we study the locus of general real forms whose real rank coincides with the complex rank. In case of quaternary quadrics the boundary of this locus is a degree forty hypersurface $J(\sigma_3(v_3(\mathbb{P}^3)),\tau(v_3(\mathbb{P}^3)))$.


title: Concentration inequalities for Poisson point processes with application to adaptive intensity estimation

abstract: We derive concentration inequalities for maxima of empirical processes associated with Poisson point processes. The proofs are based on a careful application of Ledoux's entropy method. We demonstrate the utility of the obtained concentration inequalities by application to adaptive intensity estimation.


title: LSE Precoders for Massive MIMO with Signal Constraints: Fundamental Limits

abstract: This paper proposes the nonlinear Least Square Error (LSE) precoders for multiuser MIMO broadcast channels. The output signals of LSE Precoders are limited to be chosen from a predefined set which let these precoders address several constraints such as peak power limitation, constant envelope transmission and discrete constellations. We study the large-system performance of these precoders via the replica method from statistical physics, and derive a closed-form expression for the asymptotic distortion. Our results demonstrate that an LSE precoder with the output peak-to-average power ratio of $3~{\rm dB}$ can track the performance of the Regularized Zero Forcing (RZF) precoder closely. As the peak-to-average power ratio reduces to one, the constant envelope precoder is recovered. The investigations depict that the performance of the RZF precoder is achieved by the constant envelope precoder with $20\%$ of more transmit antennas. For $M$-PSK constellations, our analysis gives a lower-bound on the asymptotic distortion which is tight for moderate antenna-to-user ratios and deviates as the ratio grows. We improve this bound by deriving the replica solution under one-step of replica symmetry breaking. Our numerical investigations for this case show that the bound is tight for antenna-to-user ratios less than $5$.


title: Long and Short Memory in Economics: Fractional-Order Difference and Differentiation

abstract: Long and short memory in economic processes is usually described by the so-called discrete fractional differencing and fractional integration. We prove that the discrete fractional differencing and integration are the Grunwald-Letnikov fractional differences of non-integer order d. Equations of ARIMA(p,d,q) and ARFIMA(p,d,q) models are the fractional-order difference equations with the Grunwald-Letnikov differences of order d. We prove that the long and short memory with power law should be described by the exact fractional-order differences, for which the Fourier transform demonstrates the power law exactly. The fractional differencing and the Grunwald-Letnikov fractional differences cannot give exact results for the long and short memory with power law, since the Fourier transform of these discrete operators satisfy the power law in the neighborhood of zero only. We prove that the economic processes with the continuous time long and short memory, which is characterized by the power law, should be described by the fractional differential equations.


title: On Garland's vanishing theorem for $\mathrm{SL}_n$

abstract: This is an expository paper on Garland's vanishing theorem specialized to the case when the linear algebraic group is $\mathrm{SL}_n$. Garland's theorem can be stated as a vanishing of the cohomology groups of certain finite simplicial complexes. The method of the proof is quite interesting on its own. It relates the vanishing of cohomology to the assertion that the minimal positive eigenvalue of a certain combinatorial laplacian is sufficiently large. Since the 1970's, this idea has found applications in a variety of problems in representation theory, group theory, and combinatorics, so the paper might be of interest to a wide audience. The paper is intended for non-specialists and graduate students.


title: Programmable Active Janus Droplets Driven by Water/Alcohol Phase Separation

abstract: We report the existence of self-propelled Janus droplets driven by phase separation, which are able to deliver cargo in a programmable manner. The self-propelling droplets are initially formed by a water/ethanol mixture in a squalane/monoolein solution, and evolve in up to three stages depending on ethanol concentration. In the first stage, the droplet propulsion is generated by Marangoni flow originating from the solubilization of ethanol in the oily phase. During this process the droplets absorb surfactant molecules; in combination with the continuous loss of ethanol this leads to a phase separation of the water/ethanol/monoolein mixture and the formation of Janus droplets, i.e. a water-rich droplet connected to an ethanol-rich droplet that is able to deliver cargo. We characterize the different evolution stages of self-propulsion by the flow field around the droplet that evolves from a weak pusher, over a neutral swimmer, to a dimer of neutral swimmers. Finally, we utilize this active system to deliver DNA as a cargo. Tuning the delay time before phase separation, by varying the chemical composition of the droplets, several different cargo delivery processes can be programmed.


title: Optical and Near-IR study of nova V2676 Oph 2012

abstract: We present optical spectrophotometric and near-infrared (NIR) photometric observations of the nova V2676 Oph covering the period from 2012 March 29 through 2015 May 8. The optical spectra and photometry of the nova have been taken from SMARTS and Asiago; the near-infrared photometry was obtained from SMARTS and Mt. Abu. The spectra are dominated by strong H I lines from the Balmer series, Fe II, N I and [O I] lines in the initial days, typical of an Fe II type nova. The measured FWHM for the H-beta and H-alpha lines was 800-1200 km/s. There was pronounced dust formation starting 90 days after the outburst. The J - K color was the largest among recent dust forming novae.


title: On weighted generalized functions associated with quadratic forms

abstract: In this article we consider certain types of weighted generalized functions associated with nondegenerate quadratic forms. Such functions and their derivatives are used for constructing fundamental solutions of iterated ultra-hyperbolic equations with Bessel operator and for constructing negative real powers of ultra-hyperbolic operators with Bessel operator.


title: Tradeoffs between microbial growth phases lead to frequency-dependent and non-transitive selection

abstract: Mutations in a microbial population can increase the frequency of a genotype not only by increasing its exponential growth rate, but also by decreasing its lag time or adjusting the yield (resource efficiency). The contribution of multiple life-history traits to selection is a critical question for evolutionary biology as we seek to predict the evolutionary fates of mutations. Here we use a model of microbial growth to show there are two distinct components of selection corresponding to the growth and lag phases, while the yield modulates their relative importance. The model predicts rich population dynamics when there are tradeoffs between phases: multiple strains can coexist or exhibit bistability due to frequency-dependent selection, and strains can engage in rock-paper-scissors interactions due to non-transitive selection. We characterize the environmental conditions and patterns of traits necessary to realize these phenomena, which we show to be readily accessible to experiments. Our results provide a theoretical framework for analyzing high-throughput measurements of microbial growth traits, especially interpreting the pleiotropy and correlations between traits across mutants. This work also highlights the need for more comprehensive measurements of selection in simple microbial systems, where the concept of an ordinary fitness landscape breaks down.