title: Image-Guided Depth Sampling and Reconstruction

abstract: Depth acquisition, based on active illumination, is essential for autonomous and robotic navigation. LiDARs (Light Detection And Ranging) with mechanical, fixed, sampling templates are commonly used in today's autonomous vehicles. An emerging technology, based on solid-state depth sensors, with no mechanical parts, allows fast, adaptive, programmable scans. In this paper, we investigate the topic of adaptive, image-driven, sampling and reconstruction strategies. First, we formulate a piece-wise linear depth model with several tolerance parameters and estimate its validity for indoor and outdoor scenes. Our model and experiments predict that, in the optimal case, about 20-60 piece-wise linear structures can approximate well a depth map. This translates to a depth-to-image sampling ratio of about 1/1200. We propose a simple, generic, sampling and reconstruction algorithm, based on super-pixels. We reach a sampling rate which is still far from the optimal case. However, our sampling improves grid and random sampling, consistently, for a wide variety of reconstruction methods. Moreover, our proposed reconstruction achieves state-of-the-art results, compared to image-guided depth completion algorithms, reducing the required sampling rate by a factor of 3-4. A single-pixel depth camera built in our lab illustrates the concept.


title: A countable partition for singular flows, and its application on the entropy theory

abstract: In this paper, we construct a countable partition $\mathscr{A}$ for flows with hyperbolic singularities by introducing a new cross section at each singularity. Such partition forms a Kakutani tower in a neighborhood of the singularity, and turns out to have finite metric entropy for every invariant probability measure. Moreover, each element of $\mathscr{A}^\infty$ will stay in a scaled tubular neighborhood for arbitrarily long time. This new construction enables us to study entropy theory for singular flows away from homoclinic tangencies, and show that the entropy function is upper semi-continuous with respect to both invariant measures and the flows.


title: On Flying Backwards: Preventing Run-away of Small, Low-speed, Fixed-wing UAVs in Strong Winds

abstract: Small, low-speed fixed-wing Unmanned Aerial Vehicles (UAVs) operating autonomously, beyond-visual-line-of-sight (BVLOS) will inevitably encounter winds rising to levels near or exceeding the vehicles' nominal airspeed. In this paper, we develop a nonlinear lateral-directional path following guidance law with explicit consideration of online wind estimates. Energy efficient airspeed reference compensation logic is developed for excess wind scenarios (i.e. when the wind speed rises above the airspeed), enabling either mitigation, prevention, or over-powering of excess wind induced run-away from a given path. The developed guidance law is demonstrated on a representative small, low-speed test UAV in two flight experiments conducted in mountainous regions of Switzerland with strong, turbulent wind conditions, gusts reaching up to 13 meters per second. We demonstrate track-keeping errors of less than 1 meter consistently maintained during a representative duration of gusting, excess winds and a mean ground speed undershoot of 0.5 meters per second from the commanded minimum forward ground speed demonstrated in over 5 minutes of the showcased flight results.


title: Permutations avoiding a pattern of length three under Mallows distributions

abstract: We consider permutations avoiding a pattern of length three under the family of Mallows distributions. In particular, for any pattern $\tau\in S_3-\{321\}$, we obtain rather precise results on the asymptotic probability as $n\to\infty$ that a permutation $\sigma\in S_n$ under the Mallows distribution with parameter $q\in(0,1)$ avoids the pattern. By a duality between the parameters $q$ and $\frac1q$, we also obtain rather precise results on the above probability for $q>1$ and any pattern $\tau\in S_3-\{123\}$.


title: Slice Dirac operator over octonions

abstract: The slice Dirac operator over octonions is a slice counterpart of the Dirac operator over quaternions. It involves a new theory of stem functions, which is the extension from the commutative $ O(1) $ case to the non-commutative $ O(3) $ case. For functions in the kernel of the slice Dirac operator over octonions, we establish the representation formula, the Cauchy integral formula (and, more in general, the Cauchy-Pompeiu formula), and the Taylor as well as the Laurent series expansion formulas.


title: Simultaneous Clustering and Optimization for Evolving Datasets

abstract: Simultaneous clustering and optimization (SCO) has recently drawn much attention due to its wide range of practical applications. Many methods have been previously proposed to solve this problem and obtain the optimal model. However, when a dataset evolves over time, those existing methods have to update the model frequently to guarantee accuracy; such updating is computationally infeasible. In this paper, we propose a new formulation of SCO to handle evolving datasets. Specifically, we propose a new variant of the alternating direction method of multipliers (ADMM) to solve this problem efficiently. The guarantee of model accuracy is analyzed theoretically for two specific tasks: ridge regression and convex clustering. Extensive empirical studies confirm the effectiveness of our method.


title: A Convergence Result for Dirichlet Semigroups on Tubular Neighbourhoods and the Marginals of Conditional Brownian Motion

abstract: We investigate yet another approach to understand the limit behaviour of Brownian motion conditioned to stay within a tubular neighbourhood around a closed and connected submanifold of a Riemannian manifold. In this context, we identify a second order generator subject to Dirichlet conditions on the boundary of the tube and study its associated semigroups. After a suitable rescaling and renormalization procedure, we obtain convergence of these semigroups, both in $L^2$ and in Sobolev spaces of arbitrarily large index, to a limit semigroup, as the tube diameter tends to zero. As a byproduct, we conclude that the conditional Brownian motion converges in finite dimensional distributions to a limit process supported by the path space of the submanifold.


title: Reconstruction of the magnetic field for a Schr\"odinger operator in a cylindrical setting

abstract: In this thesis we consider a magnetic Schr\"odinger inverse problem over a compact domain contained in an infinite cylindrical manifold. We show that, under certain conditions on the electromagnetic potentials, we can recover the magnetic field from boundary measurements in a constructive way. A fundamental tool for this procedure is a global Carleman estimate for the magnetic Schr\"odinger operator. We prove this by conjugating the magnetic operator essentially into the Laplacian, and using the Carleman estimates for it proven by Kenig-Salo-Uhlmann in the anisotropic setting, see [KSU11a]. The conjugation is achieved through pseudodifferential operators over the cylinder, for which we develop the necessary results. The main motivations to attempt this question are the following results concerning the magnetic Schr\"odinger operator: first, the solution to the uniqueness problem in the cylindrical setting in [DSFKSU09], and, second, the reconstruction algorithm in the Euclidean setting from [Sal06]. We will also borrow ideas from the reconstruction of the electric potential in the cylindrical setting from [KSU11b]. These two new results answer partially the Carleman estimate problem (Question 4.3.) proposed in [Sal13] and the reconstruction for the magnetic Schr\"odinger operator mentioned in the introduction of [KSU11b]. To our knowledge, these are the first global Carleman estimates and reconstruction procedure for the magnetic Schr\"odinger operator available in the cylindrical setting.


title: A Sub-Gaussian estimate for Dirichlet Heat Kernels on Tubular Neighbourhoods and Tightness of Conditional Brownian Motion

abstract: We prove tightness of a family of path measures $\nu_{\varepsilon}$ on tubes $L(\varepsilon)$ of small diameters around a closed and connected submanifold $L$ of another Riemannian manifold $M$. Together with a convergence result for Dirichlet semigroups on tubular neighbourhoods, that implies weak convergence of the measures as the tube radius $\varepsilon$ tends to zero to a measure supported by the path space of the submanifold. As a consequence, we obtain weak convergence of the measures obtained by conditioning Brownian motion to stay within the tubes $L(\varepsilon)$ up to a finite time $T>0$, and we identify the limit measure.


title: Pairwise Multi-marginal Optimal Transport and Embedding for Earth Mover's Distance

abstract: We investigate the problem of pairwise multi-marginal optimal transport, that is, given a collection of probability distributions $\{P_\alpha\}$ on a Polish space $\mathcal{X}$, to find a coupling $\{X_\alpha\}$, $X_\alpha\sim P_\alpha$, such that $\mathbf{E}[c(X_\alpha,X_\beta)]\le r\inf_{X\sim P_\alpha,Y\sim P_\beta}\mathbf{E}[c(X,Y)]$ for all $\alpha,\beta$, where $c$ is a cost function and $r\ge1$. In other words, every pair $(X_\alpha,X_\beta)$ has an expected cost at most a factor of $r$ from its lowest possible value. This can be regarded as a locality sensitive hash function for probability distributions, and has applications such as robust and distributed computation of transport plans. It can also be considered as a bi-Lipschitz embedding of the collection of probability distributions into the space of random variables taking values on $\mathcal{X}$. For $c(x,y)=\Vert x-y\Vert_2^q$ on $\mathbb{R}^n$, where $q>0$, we show that a finite $r$ is attainable if and only if either $n=1$ or $0<q<1$. As $n\to\infty$, the growth rate of the smallest possible $r$ is exactly $\Theta(n^{q/2})$ if $0<q<1$. Hence, the metric space of probability distributions on $\mathbb{R}^n$ with finite $q$-th absolute moments, $0<q<1$, with the earth mover's distance (or 1-Wasserstein distance) with respect to the snowflake metric $c(x,y)=\Vert x-y\Vert_2^q$, is bi-Lipschitz embeddable into $L_1$ with distortion $O(n^{q/2})$. If we consider $c(x,y)=\Vert x-y\Vert_2$ (i.e., $q=1$) on the grid $[0..s]^n$ instead of $\mathbb{R}^n$, then $r=O(\sqrt{n}\log s)$ is attainable, which implies the embeddability of the space of probability distributions on $[0..s]^n$ into $L_1$ with distortion $O(\sqrt{n}\log s)$, and improves upon the $O(n\log s)$ result by Indyk and Thaper. The case of the discrete metric cost $c(x,y)=\mathbf{1}\{x\neq y\}$ and more general metric and ultrametric costs are also investigated.


title: Lorentz Breaking and SU(2)_L x U(1)_Y Gauge Invariance for Neutrino Decays

abstract: Conceivable Lorentz-violating effects in the neutrino sector remain a research area of great general interest, as they touch upon the very foundations on which the Standard Model and our general understanding of fundamental interactions is laid. Here, we investigate the relation of Lorentz violation in the neutrino sector in light of the fact that neutrinos and corresponding left-handed charged leptons form SU(2)_L doublets under the electroweak gauge group. Lorentz-violating effects thus cannot be fully separated from questions related to gauge invariance. The model dependence of the effective interaction Lagrangians used in various recent investigations is investigated, with a special emphasis on neutrino splitting, otherwise known as neutrino-pair Cerenkov radiation, NPCR, and vacuum pair emission (electron-positron-pair Cerenkov radiation, LPCR). We investigate two scenarios in which Lorentz violating effects do not necessarily also break electroweak gauge invarianceThe first of these involves a restricted set of gauge transformation, a subgroup of SU(2)_L x U(1)_Y, while in the second, differential Lorentz violation is exclusively introduced by the mixing of the neutrino flavor and mass eigenstates. Our study culminates in a model which fully preserves SU(2)_L x U(1)_Y gauge invariance, involves flavor-dependent Lorentz-breaking parameters, and still allows for NPCR and LPCR decays to proceed.


title: Existence and $L^{\infty}$-estimates for elliptic equations involving convolution

abstract: In this paper, with a fixed $p\in (1,+\infty)$ and a bounded domain $\Omega \subset \mathbb{R}^N$ whose boundary $\partial\Omega$ fulfills the $C^1$ regularity, we study a boundary value problem involving a nonlocal operator assigning to $u$ the convolution $\rho \ast E(u)$ of $\rho$ with $E(u)$, where $\rho$ is an integrable function on $\mathbb{R}^N$ and $E$ is an extension operator related to $\Omega$. Under verifiable conditions, we prove the existence of a (weak) solution to our problem by using the surjectivity theorem for pseudomonotone operators. Moreover, through a modified version of Moser iteration up to the boundary, we show that (any) weak solution to our problem is bounded.


title: Cosmological constraints in extended parameter space from the Planck 2018 Legacy release

abstract: We present new constraints on extended cosmological scenarios using the recent data from the Planck 2018 Legacy release. In addition to the six parameters of the standard LCDM model, we also simultaneously vary the dark energy equation of state, the neutrino mass, the neutrino effective number, the running of the spectral index and the lensing amplitude $A_L$. We confirm that a resolution of the Hubble tension is given by a dark energy equation of state with w<-1, ruling out quintessence models at high statistical significance. This solution is, however, not supported by BAO and Pantheon data. We find no evidence for evolving dark energy, i.e. $w_a \neq 0$. The neutrino effective number is always in agreement with the expectations of the standard model based on three active neutrinos. The running of the spectral index also is always consistent with zero. Despite the increase in the number of parameters, the $A_L$ lensing anomaly is still present at more than two standard deviations. The $A_L$ anomaly significantly affects the bounds on the neutrino mass that can be larger by a factor four with respect to those derived under standard $\Lambda$CDM. While the lensing data reduces the evidence for $A_L>1$, the inclusion of BAO and Pantheon increase its statistical significance.


title: Discrete Kontorovich-Lebedev transforms

abstract: Discrete analogs of the classical Kontorovich-Lebedev transforms are introduced and investigated. It involves series with the modified Bessel function or Macdonald function $K_{in}(x), x >0, n \in \mathbb{N}, i $ is the imaginary unit, and incomplete Bessel functions. Several expansions of suitable functions and sequences in terms of these series and integrals are established. As an application, a Dirichlet boundary value problem in the upper half-plane for inhomogeneous Helmholtz equation is solved.


title: Network Inference from Consensus Dynamics with Unknown Parameters

abstract: We explore the problem of inferring the graph Laplacian of a weighted, undirected network from snapshots of a single or multiple discrete-time consensus dynamics, subject to parameter uncertainty, taking place on the network. Specifically, we consider three problems in which we assume different levels of knowledge about the diffusion rates, observation times, and the input signal power of the dynamics. To solve these underdetermined problems, we propose a set of algorithms that leverage the spectral properties of the observed data and tools from convex optimization. Furthermore, we provide theoretical performance guarantees associated with these algorithms. We complement our theoretical work with numerical experiments, that demonstrate how our proposed methods outperform current state-of-the-art algorithms and showcase their effectiveness in recovering both synthetic and real-world networks.


title: Learning to Transport with Neural Networks

abstract: We compare several approaches to learn an Optimal Map, represented as a neural network, between probability distributions. The approaches fall into two categories: ``Heuristics'' and approaches with a more sound mathematical justification, motivated by the dual of the Kantorovitch problem. Among the algorithms we consider a novel approach involving dynamic flows and reductions of Optimal Transport to supervised learning.


title: Direct Observation of Infrared Plasmonic Fano Antiresonances by a Nanoscale Electron Probe

abstract: In this Letter, we exploit recent breakthroughs in monochromated aberration-corrected scanning transmission electron microscopy (STEM) to resolve infrared plasmonic Fano antiresonances in individual nanofabricated disk-rod dimers. Using a combination of electron energy-loss spectroscopy (EELS) and theoretical modeling, we investigate and characterize a subspace of the weak coupling regime between quasi-discrete and quasi-continuum localized surface plasmon resonances where infrared plasmonic Fano antiresonances appear. This work illustrates the capability of STEM instrumentation to experimentally observe nanoscale plasmonic responses that were previously the domain only of higher resolution infrared spectroscopies.


title: Experimental observation of coupled valley and spin Hall effect in p-doped WSe2 devices

abstract: Giant spin Hall effect (GSHE) has been observed in heavy metal materials such as Ta, Pt, and W, where spins are polarized in the surface plane and perpendicular to the charge current direction. Spins generated in these materials have successfully switched magnets with in-plane magnetic anisotropy (IMA) and perpendicular magnetic anisotropy (PMA) through spin orbit torque (SOT) mechanism. It is generally accepted that PMA magnets are preferred over IMA magnets in data storage applications owing to their large thermal stability even at ultra scaled dimensions. However, SOT switching of PMA magnets by conventional GSHE materials requires either a small external magnetic field, a local dipolar field, or introducing tilted anisotropy to break the symmetry with respect to the magnetization. To deterministically switch a PMA without any additional assistance, nonconventional GSHE materials that can generate spins with polarization perpendicular to the surfaces are needed. Several monolayer transition metal dichalcogenides (TMDs) have been predicted to generate such out of plane spins due to their 2D nature and unique band structures. Interestingly, opposite spins are locked to their respective sub-band in each K valley of the TMD valence band with substantially large energy splitting, which enables polarized spins to be accessible through electrical gating and spatially separated by electric field through the valley Hall effect (VHE). Therefore, spatial separation and accumulation of spins in these 2D TMDs are uniquely referred to as coupled valley and spin Hall effect. Here, we report an experiment of electrical generation of spin current with out of plane polarization in monolayer WSe2 and detection of spin signals through a nonlocal spin valve structure built on a lateral graphene spin diffusion channel that partially overlaps with WSe2.


title: Notes on the norm of pre-Schwarzian derivatives on bi-univalent functions of order $\alpha$

abstract: In the present paper we estimate the norm of the pre-Schwarzian derivative of bi-starlike functions of order $\alpha$ where $\alpha\in[0,1)$. Initially this problem was handled by Rahmatan et al. in [Bull Iran Math Soc {\bf43}: 1037-1043, 2017]. We pointed out that the proofs and bounds by Rahmatan et al. are incorrect and present correct proofs and bounds.


title: The distinct flavors of Zipf's law in the rank-size and in the size-distribution representations, and its maximum-likelihood fitting

abstract: In the last years, researchers have realized the difficulties of fitting power-law distributions properly. These difficulties are higher in Zipf's systems, due to the discreteness of the variables and to the existence of two representations for these systems, i.e., two versions about which is the random variable to fit. The discreteness implies that a power law in one of the representations is not a power law in the other, and vice versa. We generate synthetic power laws in both representations and apply a state-of-the-art fitting method (based on maximum-likelihood plus a goodness-of-fit test) for each of the two random variables. It is important to stress that the method does not fit the whole distribution, but the tail, understood as the part of a distribution above a cut-off that separates non-power-law behavior from power-law behavior. We find that, no matter which random variable is power-law distributed, the rank-size representation is not adequate for fitting, whereas the representation in terms of the distribution of sizes leads to the recovery of the simulated exponents, may be with some bias.


title: Sound Event Detection in Multichannel Audio using Convolutional Time-Frequency-Channel Squeeze and Excitation

abstract: In this study, we introduce a convolutional time-frequency-channel "Squeeze and Excitation" (tfc-SE) module to explicitly model inter-dependencies between the time-frequency domain and multiple channels. The tfc-SE module consists of two parts: tf-SE block and c-SE block which are designed to provide attention on time-frequency and channel domain, respectively, for adaptively recalibrating the input feature map. The proposed tfc-SE module, together with a popular Convolutional Recurrent Neural Network (CRNN) model, are evaluated on a multi-channel sound event detection task with overlapping audio sources: the training and test data are synthesized TUT Sound Events 2018 datasets, recorded with microphone arrays. We show that the tfc-SE module can be incorporated into the CRNN model at a small additional computational cost and bring significant improvements on sound event detection accuracy. We also perform detailed ablation studies by analyzing various factors that may influence the performance of the SE blocks. We show that with the best tfc-SE block, error rate (ER) decreases from 0.2538 to 0.2026, relative 20.17\% reduction of ER, and 5.72\% improvement of F1 score. The results indicate that the learned acoustic embeddings with the tfc-SE module efficiently strengthen time-frequency and channel-wise feature representations to improve the discriminative performance.


title: Grand Unified Theories in Renormalisable, Classically Scale Invariant Gravity

abstract: We analyze $SO(N)$ and $SU(N)$ gauge theories with scalars in adjoint and fundamental representations, coupled to renormalisable, classically scale invariant gravity. In the specific case of $SO(12),$ we show that the quantum field theory can be can be asymptotically free in all couplings (hence ultra-violet complete). For a region of parameter space, Dimensional Transmutation occurs, with the adjoint vacuum expectation value breaking $SO(12) \to SU(6)\otimes U(1)$ and producing a Low Energy Effective Theory having Einstein-Hilbert gravity. We verify that certain minima are locally stable and lie within the catchment basin of the ultraviolet fixed points.


title: Understanding the Nature of Mean-Field Semiclassical Light-Matter Dynamics: An Investigation of Energy Transfer, Electron-Electron Correlations, External Driving and Long-Time Detailed Balance

abstract: Semiclassical electrodynamics is an appealing approach for studying light-matter interactions, especially for realistic molecular systems. However, there is no unique semiclassical scheme. On the one hand, intermolecular interactions can be described instantaneously by static two-body interactions connecting different molecules plus a classical transverse E-field; we will call this Hamiltonian #I. On the other hand, intermolecular interactions can also be described as effects that are mediated exclusively through a classical one-body E-field without any quantum effects at all (assuming we ignore electronic exchange); we will call this Hamiltonian #II. Moreover, one can also mix these two Hamiltonians into a third, hybrid Hamiltonian, which preserves quantum electron-electron correlations for lower excitations but describes higher excitations in a mean-field way. To investigate which semiclassical scheme is most reliable for practical use, here we study the real-time dynamics of a pair of identical two-level systems (TLSs) undergoing either resonance energy transfer (RET) or collectively driven dynamics. While all approaches perform reasonably well when there is no strong external excitation, we find that no single approach is perfect for all conditions. Each method has its own distinct problems: Hamiltonian #I performs best for RET but behaves in a complicated manner for driven dynamics. Hamiltonian #II is always stable, but obviously fails for RET at short distances. One key finding is that, under externally driving, a full configuration interaction description of Hamiltonian #I strongly overestimates the long-time electronic energy, highlighting the not obvious fact that, if one plans to merge quantum molecules with classical light, a full, exact treatment of electron-electron correlations can actually lead to worse results than a simple mean-field treatment.


title: Multi-block Bregman proximal alternating linearized minimization and its application to orthogonal nonnegative matrix factorization

abstract: We introduce and analyze BPALM and A-BPALM, two multi-block proximal alternating linearized minimization algorithms using Bregman distances for solving structured nonconvex problems. The objective function is the sum of a multi-block relatively smooth function (i.e., relatively smooth by fixing all the blocks except one) and block separable (nonsmooth) nonconvex functions. It turns out that the sequences generated by our algorithms are subsequentially convergent to critical points of the objective function, while they are globally convergent under KL inequality assumption. Further, the rate of convergence is further analyzed for functions satisfying the {\L}ojasiewicz's gradient inequality. We apply this framework to orthogonal nonnegative matrix factorization (ONMF) that satisfies all of our assumptions and the related subproblems are solved in closed forms, where some preliminary numerical results is reported.


title: Deep Neural Network for Semantic-based Text Recognition in Images

abstract: State-of-the-art text spotting systems typically aim to detect isolated words or word-by-word text in images of natural scenes and ignore the semantic coherence within a region of text. However, when interpreted together, seemingly isolated words may be easier to recognize. On this basis, we propose a novel "semantic-based text recognition" (STR) deep learning model that reads text in images with the help of understanding context. STR consists of several modules. We introduce the Text Grouping and Arranging (TGA) algorithm to connect and order isolated text regions. A text-recognition network interprets isolated words. Benefiting from semantic information, a sequenceto-sequence network model efficiently corrects inaccurate and uncertain phrases produced earlier in the STR pipeline. We present experiments on two new distinct datasets that contain scanned catalog images of interior designs and photographs of protesters with hand-written signs, respectively. Our results show that our STR model outperforms a baseline method that uses state-of-the-art single-wordrecognition techniques on both datasets. STR yields a high accuracy rate of 90% on the catalog images and 71% on the more difficult protest images, suggesting its generality in recognizing text.


title: Optimistic planning for the near-optimal control of nonlinear switched discrete-time systems with stability guarantees

abstract: Originating in the artificial intelligence literature, optimistic planning (OP) is an algorithm that generates near-optimal control inputs for generic nonlinear discrete-time systems whose input set is finite. This technique is therefore relevant for the near-optimal control of nonlinear switched systems, for which the switching signal is the control. However, OP exhibits several limitations, which prevent its application in a standard control context. First, it requires the stage cost to take values in [0,1], an unnatural prerequisite as it excludes, for instance, quadratic stage costs. Second, it requires the cost function to be discounted. Third, it applies for reward maximization, and not cost minimization. In this paper, we modify OP to overcome these limitations, and we call the new algorithm OPmin. We then make stabilizability and detectability assumptions, under which we derive near-optimality guarantees for OPmin and we show that the obtained bound has major advantages compared to the bound originally given by OP. In addition, we prove that a system whose inputs are generated by OPmin in a receding-horizon fashion exhibits stability properties. As a result, OPmin provides a new tool for the near-optimal, stable control of nonlinear switched discrete-time systems for generic cost functions.


title: Programmable In-Network Security for Context-aware BYOD Policies

abstract: Bring Your Own Device (BYOD) has become the new norm in enterprise networks, but BYOD security remains a top concern. Context-aware security, which enforces access control based on dynamic runtime context, holds much promise. Recent work has developed SDN solutions to collect device context for network-wide access control in a central controller. However, the central controller poses a bottleneck that can become an attack target, and processing context changes at remote software has low agility. We present a new paradigm, programmable in-network security (Poise), which is enabled by the emergence of programmable switches. At the heart of Poise is a novel switch primitive, which can be programmed to support a wide range of context-aware policies in hardware. Users of Poise specify concise policies, and Poise compiles them into different instantiations of the security primitive in P4. Compared to centralized SDN defenses, Poise is resilient to control plane saturation attacks, and it dramatically increases defense agility.


title: Uncertainty in the Hot Hand Fallacy: Detecting Streaky Alternatives to Random Bernoulli Sequences

abstract: We study a class of permutation tests of the randomness of a collection of Bernoulli sequences and their application to analyses of the human tendency to perceive streaks of consecutive successes as overly representative of positive dependence - the hot hand fallacy. In particular, we study permutation tests of the null hypothesis of randomness (i.e., that trials are i.i.d.) based on test statistics that compare the proportion of successes that directly follow k consecutive successes with either the overall proportion of successes or the proportion of successes that directly follow k consecutive failures. We characterize the asymptotic distributions of these test statistics and their permutation distributions under randomness, under a set of general stationary processes, and under a class of Markov chain alternatives, which allow us to derive their local asymptotic power. The results are applied to evaluate the empirical support for the hot hand fallacy provided by four controlled basketball shooting experiments. We establish that substantially larger data sets are required to derive an informative measurement of the deviation from randomness in basketball shooting. In one experiment, for which we were able to obtain data, multiple testing procedures reveal that one shooter exhibits a shooting pattern significantly inconsistent with randomness - supplying strong evidence that basketball shooting is not random for all shooters all of the time. However, we find that the evidence against randomness in this experiment is limited to this shooter. Our results provide a mathematical and statistical foundation for the design and validation of experiments that directly compare deviations from randomness with human beliefs about deviations from randomness, and thereby constitute a direct test of the hot hand fallacy.


title: GraphBLAST: A High-Performance Linear Algebra-based Graph Framework on the GPU

abstract: High-performance implementations of graph algorithms are challenging to implement on new parallel hardware such as GPUs because of three challenges: (1) the difficulty of coming up with graph building blocks, (2) load imbalance on parallel hardware, and (3) graph problems having low arithmetic intensity. To address some of these challenges, GraphBLAS is an innovative, on-going effort by the graph analytics community to propose building blocks based on sparse linear algebra, which will allow graph algorithms to be expressed in a performant, succinct, composable and portable manner. In this paper, we examine the performance challenges of a linear-algebra-based approach to building graph frameworks and describe new design principles for overcoming these bottlenecks. Among the new design principles is exploiting input sparsity, which allows users to write graph algorithms without specifying push and pull direction. Exploiting output sparsity allows users to tell the backend which values of the output in a single vectorized computation they do not want computed. Load-balancing is an important feature for balancing work amongst parallel workers. We describe the important load-balancing features for handling graphs with different characteristics. The design principles described in this paper have been implemented in "GraphBLAST", the first high-performance linear algebra-based graph framework on NVIDIA GPUs that is open-source. The results show that on a single GPU, GraphBLAST has on average at least an order of magnitude speedup over previous GraphBLAS implementations SuiteSparse and GBTL, comparable performance to the fastest GPU hardwired primitives and shared-memory graph frameworks Ligra and Gunrock, and better performance than any other GPU graph framework, while offering a simpler and more concise programming model.


title: Defence Against the Modern Arts: the Curse of Statistics -- FRStat

abstract: For several decades, legal and scientific scholars have argued that conclusions from forensic examinations should be supported by statistical data and reported within a probabilistic framework. Multiple models have been proposed to quantify the probative value of forensic evidence. Unfortunately, several of these models rely on ad-hoc strategies that are not scientifically sound. The opacity of the technical jargon that is used to present these models and their results and the complexity of the techniques involved make it very difficult for the untrained user to separate the wheat from the chaff. This series of paper is intended to help forensic scientists and lawyers recognise issues in tools proposed to interpret the results of forensic examinations. This paper focuses on the tool proposed by the Latent Print Branch of the U.S. Defense Forensic Science Center (DFSC) and called FRStat. In this paper, I explore the compatibility of the results outputted by FRStat with the language used by the DFCS to report the conclusions of their fingerprint examinations, as well as the appropriateness of the statistical modelling underpinning the tool and the validation of its performance.


title: SDSS IV MaNGA -- Star-Formation Driven Biconical Outflows in the Local Univers

abstract: We present a sample of 48 nearby galaxies with central, biconical outflows identified by the Mapping Nearby Galaxies at APO (MaNGA) survey. All considered galaxies have star formation driven bi-conical central outflows (SFB), with no signs of AGN. We find that the SFB outflows require high central concentration of the star formation rate. This increases the gas velocity dispersion over the equilibrium limit and helps maintain the gas outflows. The central starbursts increase the metallicity, extinction, and the alpha/Fe ratio in the gas. Significant amount of young stellar population at the centers suggests that the SFBs are associated with the formation of young bulges in galaxies. More than 70% of SFB galaxies are group members or have companions with no prominent interaction, or show asymmetry of external isophotes. In 15% SFB cases stars and gas rotate in the opposite directions, which points at the gas infall from satellites as the primary reason for triggering the SFB phenomena.


title: On generic double shuffle relations, localized multiple polylogarithms and algebraic functions

abstract: We define subvarieties of $\mathcal{M}_{0,n}$ equipped with algebraic functions that are solutions to the generic double shuffle equations satisfied by multiple polylogarithms on $\mathcal{M}_{0,n}$.


title: Effect of Interim Adaptations in Group Sequential Designs

abstract: This manuscript investigates unconditional and conditional-on-stopping maximum likelihood estimators (MLEs), information measures and information loss associated with conditioning in group sequential designs (GSDs). The possibility of early stopping brings truncation to the distributional form of MLEs; sequentially, GSD decisions eliminate some events from the sample space. Multiple testing induces mixtures on the adapted sample space. Distributions of MLEs are mixtures of truncated distributions. Test statistics that are asymptotically normal without GSD, have asymptotic distributions, under GSD, that are non-normal mixtures of truncated normal distributions under local alternatives; under fixed alternatives, asymptotic distributions of test statistics are degenerate. Estimation of various statistical quantities such as information, information fractions, and confidence intervals should account for the effect of planned adaptations. Calculation of adapted information fractions requires substantial computational effort. Therefore, a new GSD is proposed in which stage-specific sample sizes are fully determined by desired operational characteristics, and calculation of information fractions is not needed.


title: Effects of quantum fluctuations on the dynamics of dipolar Bose polarons

abstract: We study the dynamics of dipolar Bose polarons in the presence of the normal and anomalous fluctuations using the time-dependent Hartree-Fock-Bogoliubov theory. The density profiles of the condensate, the anomalous component and the impurity are deeply analyzed. The time evolution of the width and the center-of-mass oscillation of such quantities is also highlighted. We calculate corrections due to quantum fluctuations and impurity to the chemical potential and the radius of the condensate and of the anomalous component in the weak coupling regime using the Thomas-Fermi approximation. Effects of the dipole-dipole interaction, impurity-host interaction and the anomalous fluctuation on the width and on the breathing frequencies of the impurity are discussed by variational and numerical means.


title: Almost everywhere uniqueness of blow-up limits for the lower dimensional obstacle problem

abstract: We answer a question left open in [Arch. Rat. Mech. Anal. 230 (1) (2018), 125-184] and [Arch. Rat. Mech. Anal. 230 (2) (2018), 783-784], by proving that the blow-up limit of minimizers $u$ of the lower dimensional obstacle problem is unique at generic point of the free-boundary.


title: Singularities and Genus of the k-Ellipse

abstract: A k-ellipse is a plane curve consisting of all points whose distances from k fixed foci sum to a constant. We determine the singularities and genus of its Zariski closure in the complex projective plane. The paper resolves an open problem stated by Nie, Parrilo and Sturmfels in 2008.


title: Toric ideals of Minkowski sums of unit simplices

abstract: In this paper, we discuss the toric ideals of Minkowski sums of unit simplices. More precisely, we prove that the toric ideal of Minkowski sum of unit simplices has a squarefree initial ideal and is generated by quadratic binomials. Moreover, we also prove that Minkowski sums of unit simplices have the integer decomposition property. Those results are a partial contribution to Oda conjecture and B{\o}gvad conjecture.


title: Covers of rational double points in mixed characteristic

abstract: We further the classification of rational surface singularities. Suppose $(S, \mathfrak{n}, \mathcal{k})$ is a strictly Henselian regular local ring of mixed characteristic $(0, p > 5)$. We classify functions $f$ for which $S/(f)$ has an isolated rational singularity at the maximal ideal $\mathfrak{n}$. The classification of such functions are used to show that if $(R, \mathfrak{m}, \mathcal{k})$ is an excellent, strictly Henselian, Gorenstein rational singularity of dimension $2$ and mixed characteristic $(0, p > 5)$, then there exists a split finite cover of $\mbox{Spec}(R)$ by a regular scheme. We give an application of our result to the study of $2$-dimensional BCM-regular singularities in mixed characteristic.


title: Automatic Playtesting for Game Parameter Tuning via Active Learning

abstract: Game designers use human playtesting to gather feedback about game design elements when iteratively improving a game. Playtesting, however, is expensive: human testers must be recruited, playtest results must be aggregated and interpreted, and changes to game designs must be extrapolated from these results. Can automated methods reduce this expense? We show how active learning techniques can formalize and automate a subset of playtesting goals. Specifically, we focus on the low-level parameter tuning required to balance a game once the mechanics have been chosen. Through a case study on a shoot-`em-up game we demonstrate the efficacy of active learning to reduce the amount of playtesting needed to choose the optimal set of game parameters for two classes of (formal) design objectives. This work opens the potential for additional methods to reduce the human burden of performing playtesting for a variety of relevant design concerns.


title: An algebra of distributions related to a star product with separation of variables

abstract: Given a star product with separation of variables $\star$ on a pseudo-K\"ahler manifold $M$ and a point $x_0 \in M$, we construct an associative algebra of formal distributions supported at $x_0$. We use this algebra to express the formal oscillatory exponents of a family of formal oscillatory integrals related to the star product $\star$.


title: Gradient sub-grid-scale model for relativistic MHD Large Eddy Simulations

abstract: MHD turbulence is likely to play an important role in several astrophysical scenarios where the magnetic Reynolds is very large. Numerically, these cases can be studied efficiently by means of Large Eddy Simulations, in which the computational resources are used to evolve the system only up to a finite grid size. The resolution is not fine enough to capture all the relevant small-scale physics at play, which is instead effectively modeled by a set of additional terms in the evolution equations, dubbed as sub-grid-scale model. Here we extend such approach, commonly used in non-relativistic/non-magnetic/incompressible fluid dynamics, applying the so-called gradient model to a general set of balance-law equations, that includes the relevant case in which a non-trivial inversion of conserved to primitive fields is needed. In particular, we focus on the relativistic compressible ideal MHD scenario, providing for the first time (and for any equation of state) all the additional sub-grid-scale terms. As an application, we consider box simulations of the relativistic Kelvin-Helmholtz instability, which is also the first mechanism responsible for the magnetic field amplification in binary neutron star mergers and cannot yet be fully captured by the finest-grid and longest simulations available. The performance of our model is numerically assessed by comparing it to the residuals arising from the filtering of high-resolution simulations. We find that the model can fit very well those residuals from resolutions a few times higher. Although the application shown here explicitly considers the Minkowski metric, it can be directly extended to general relativity, thus settling the basis to implement the gradient sub-grid model in a GRMHD binary merger. Our results suggest that this approach will be potentially able to unveil much better the small-scale dynamics achievable in full GRMHD simulations.


title: Automatic Game Design via Mechanic Generation

abstract: Game designs often center on the game mechanics---rules governing the logical evolution of the game. We seek to develop an intelligent system that generates computer games. As first steps towards this goal we present a composable and cross-domain representation for game mechanics that draws from AI planning action representations. We use a constraint solver to generate mechanics subject to design requirements on the form of those mechanics---what they do in the game. A planner takes a set of generated mechanics and tests whether those mechanics meet playability requirements---controlling how mechanics function in a game to affect player behavior. We demonstrate our system by modeling and generating mechanics in a role-playing game, platformer game, and combined role-playing-platformer game.


title: Explicit Characterization of Performance of a Class of Networked Linear Control Systems

abstract: We show that the steady-state variance as a performance measure for a class of networked linear control systems is expressible as the summation of a rational function over the Laplacian eigenvalues of the network graph. Moreover, we characterize the role of connectivity thresholds for the feedback (and observer) gain design of these networks. We use our framework to derive bounds and scaling laws for the performance of the dynamical network. Our approach generalizes and unifies the previous results on the performance measure of these networks for the case of arbitrary nodal dynamics. We bring extensions of our methodology for the case of decentralized observer-based output feedback as well as a class of composite networks. Numerous examples support our theoretical contributions.


title: Analytic and Numerical Study of Navier-Stokes Loop Equation in Turbulence

abstract: We developed analytic approach to the non-planar loop equation, which we derived in previous papers \cite{M19a},\cite{M19b},\cite{M19c}. We found quadratic integral equation for the vorticity distribution $\Omega(r)$ we introduced on a minimal surface. There are no corrections to the minimal surface though: it is still defined by mean external curvature equal to zero, for arbitrary non-planar loop. We also analyzed the loop equations with viscosity term in Navier-Stokes equations. This term creates boundary condition for $\Omega(r\in C)$. The leading viscosity correction term mixes the moments $\left< \Gamma^p \right>$ with $\left< \Gamma^{p-1} \right>$ resembling the bi-fractal behavior observed in \cite{S19} and explicitly breaking the time reversal symmetry. We also develop numerical approach to the loop equation with arbitrary curved loop and present \Mathematica notebook building triangulated minimal surface and then numerically solving these equations. As a result we obtain predictions for future numerical experiments which will compute vorticity distribution along the loop.


title: Monte-Carlo Tree Search for Simulation-based Strategy Analysis

abstract: Games are often designed to shape player behavior in a desired way; however, it can be unclear how design decisions affect the space of behaviors in a game. Designers usually explore this space through human playtesting, which can be time-consuming and of limited effectiveness in exhausting the space of possible behaviors. In this paper, we propose the use of automated planning agents to simulate humans of varying skill levels to generate game playthroughs. Metrics can then be gathered from these playthroughs to evaluate the current game design and identify its potential flaws. We demonstrate this technique in two games: the popular word game Scrabble and a collectible card game of our own design named Cardonomicon. Using these case studies, we show how using simulated agents to model humans of varying skill levels allows us to extract metrics to describe game balance (in the case of Scrabble) and highlight potential design flaws (in the case of Cardonomicon).


title: Testing Multispecies Coalescent Simulators using Summary Statistics

abstract: As genomic scale datasets motivate research on species tree inference, simulators of the multispecies coalescent (MSC) process are essential for the testing and evaluation of new inference methods. However, the simulators themselves must be tested to ensure they give valid samples from the coalescent process. In this work we develop several statistical tools using summary statistics to evaluate the fit of a simulated gene tree sample to the MSC model. Using these tests on samples from four published simulators, we uncover flaws in several. The tests are implemented as an R package, so that both developers and users will be able to easily check proper performance of future simulators.


title: ChemBO: Bayesian Optimization of Small Organic Molecules with Synthesizable Recommendations

abstract: In applications such as molecule design or drug discovery, it is desirable to have an algorithm which recommends new candidate molecules based on the results of past tests. These molecules first need to be synthesized and then tested for objective properties. We describe ChemBO, a Bayesian optimization framework for generating and optimizing organic molecules for desired molecular properties. While most existing data-driven methods for this problem do not account for sample efficiency or fail to enforce realistic constraints on synthesizability, our approach explores the synthesis graph in a sample-efficient way and produces synthesizable candidates. We implement ChemBO as a Gaussian process model and explore existing molecular kernels for it. Moreover, we propose a novel optimal-transport based distance and kernel that accounts for graphical information explicitly. In our experiments, we demonstrate the efficacy of the proposed approach on several molecular optimization problems.


title: Geometry and Generation of a New Graph Planarity Game

abstract: We introduce a new abstract graph game, Swap Planarity, where the goal is to reach a state without edge intersections and a move consists of swapping the locations of two vertices connected by an edge. We analyze this puzzle game using concepts from graph theory and graph drawing, computational geometry, and complexity. Furthermore, we specify quality criteria for puzzle instances, and describe a method to generate high-quality instances. We also report on experiments that show how well this generation process works.


title: Gevrey and formal Nilsson solutions of $A$-hypergeometric systems

abstract: We prove that the space of Gevrey solutions of an $A$--hypergeometric system along a coordinate subspace is contained in a space of formal Nilsson solutions. Moreover, under some additional conditions, both spaces are equal. In the process we prove some other results about formal Nilsson solutions.


title: Inference of visual field test performance from OCT volumes using deep learning

abstract: Visual field tests (VFT) are pivotal for glaucoma diagnosis and conducted regularly to monitor disease progression. Here we address the question to what degree aggregate VFT measurements such as Visual Field Index (VFI) and Mean Deviation (MD) can be inferred from Optical Coherence Tomography (OCT) scans of the Optic Nerve Head (ONH) or the macula. Accurate inference of VFT measurements from OCT could reduce examination time and cost. We propose a novel 3D Convolutional Neural Network (CNN) for this task and compare its accuracy with classical machine learning (ML) algorithms trained on common, segmentation-based OCT, features employed for glaucoma diagnostics. Peak accuracies were achieved on ONH scans when inferring VFI with a Pearson Correlation (PC) of 0.88$\pm$0.035 for the CNN and a significantly lower (p $<$ 0.01) PC of 0.74$\pm$0.090 for the best performing, classical ML algorithm - a Random Forest regressor. Estimation of MD was equally accurate with a PC of 0.88$\pm$0.023 on ONH scans for the CNN.


title: Restricted Linearized Augmented Lagrangian Method for Euler's Elastica Model

abstract: Euler's elastica model has been extensively studied and applied to image processing tasks. However, due to the high nonlinearity and nonconvexity of the involved curvature term, conventional algorithms suffer from slow convergence and high computational cost. Various fast algorithms have been proposed, among which, the augmented Lagrangian based ones are very popular in the community. However, parameter tuning might be very challenging for these methods. In this paper, a simple cutting-off strategy is introduced into the augmented Lagrangian based algorithms for minimizing the Euler's elastica energy, which leads to easy parameter tuning and fast convergence. The cutting-off strategy is based on an observation of inconsistency inside the augmented Lagrangian based algorithms. When the weighting parameter of the curvature term goes to zero, the energy functional boils down to the ROF model. So, a natural requirement is that its augmented Lagrangian based algorithms should also approach the augmented Lagrangian based algorithms formulated directly for solving the ROF model from the very beginning. Unfortunately, this is not the case for certain existing augmented Lagrangian based algorithms. The proposed cutting-off strategy helps to decouple the tricky dependence between the auxiliary splitting variables, so as to remove the observed inconsistency. Numerical experiments suggest that the proposed algorithm enjoys easier parameter-tuning, faster convergence and even higher quality of image restorations.


title: A note on generic Clifford algebras of binary cubic forms

abstract: We study the representation theoretic results of the binary cubic generic Clifford algebra $\mathcal C$, which is an Artin-Schelter regular algebra of global dimension five. In particular, we show that $\mathcal C$ is a PI algebra of PI degree three and compute its point variety and discriminant ideals. As a consequence, we give a necessary and sufficient condition on a binary cubic form $f$ for the associated Clifford algebra $\mathcal C_f$ to be an Azumaya algebra.


title: Ab Initio Bethe-Salpeter Equation Approach to Neutral Excitations in Molecules with Numeric Atom-Centered Orbitals

abstract: The Bethe-Salpeter equation (BSE) based on GW quasiparticle levels is a successful approach for calculating the optical gaps and spectra of solids and also for predicting the neutral excitations of small molecules. We here present an all-electron implementation of the GW+BSE formalism for molecules, using numeric atom-centered orbital (NAO) basis sets. We present benchmarks for low-lying excitation energies for a set of small organic molecules, denoted in the literature as "Thiel's set". Literature reference data based on Gaussian-type orbitals are reproduced to about one meV precision for the molecular benchmark set, when using the same GW quasiparticle energies and basis sets as the input to the BSE calculations. For valence correlation consistent NAO basis sets, as well as for standard NAO basis sets for ground state density-functional theory with extended augmentation functions, we demonstrate excellent convergence of the predicted low-lying excitations to the complete basis set limit. A simple and affordable augmented NAO basis set denoted "tier2+aug2" is recommended as a particularly efficient formulation for production calculations. We finally demonstrate that the same convergence properties also apply to linear-response time-dependent density functional theory within the NAO formalism.


title: Notes on k-rainbow independent domination in graphs

abstract: The $k$-rainbow independent domination number of a graph $G$, denoted $\gamma_{\rm rik}(G)$, is the cardinality of a smallest set consisting of two vertex-disjoint independent sets $V_1$ and $V_2$ for which every vertex in $V(G)\setminus (V_1\cup V_2)$ has neighbors in both $V_1$ and $V_2$. This domination invariant was proposed by {\v{S}}umenjak, Rall and Tepeh in (Applied Mathematics and Computation 333(15), 2018: 353-361), which allows to reduce the problem of computing the independent domination number of the generalized prism $G {\Box} K_k$ to an integer labeling problem on $G$. They proved a Nordhaus-Gaddum-type theorem: $5\leq \gamma_{\rm rik}(G)+\gamma_{\rm rik}(\overline{G})\leq n+3$ for every graph $G$ of order $n\geq 3$, where $\overline{G}$ is the complement of $G$. In this paper, we improve this result by showing that if $G$ is not isomorphic to the 5-cycle, then $5\leq \gamma_{\rm rik}(G)+\gamma_{\rm rik}(\overline{G})\leq n+2$. Moreover, we show that the problem of deciding whether a graph has a $k$-rainbow independent dominating function of a given weight is $\mathcal{NP}$-complete. Our results respond some open questions proposed by \v{S}umenjak, et al.


title: Hoffman's bound for hypergraphs

abstract: One of the best-known results in spectral graph theory is the inequality of Hoffman \[ \chi\left( G\right) \geq1-\frac{\lambda\left( G\right) }{\lambda_{\min }\left( G\right) }, \] where $\chi\left( G\right) $ is the chromatic number of a graph $G$ and $\lambda\left( G\right) ,$ $\lambda_{\min}\left( G\right) $ are the largest and the smallest eigenvalues of its adjacency matrix. In this note Hoffman's inequality is extended to weighted uniform $r$-graphs for every even $r$.


title: Testing Newton/GR, MoND and quantised inertia on wide binaries

abstract: Wide binary stars are within the low-acceleration regime in which galactic rotation curves deviate from Newtonian or general relativistic predictions. It has recently been observed that their rotation rates are similarly anomalous in a way that dark matter cannot explain, since it must be smooth on these small scales to fit galaxy rotation curves. Here, it is shown that Newtonian/GR models cannot predict these wide binaries since dark matter cannot be applied. It is also shown that MoND cannot predict these systems. However, a model which assumes that inertia is due to Unruh radiation made inhomogeneous in space by relativistic horizons (QI, quantised inertia) can predict these wide binaries, and it has the advantage of not needing an adjustable parameter.


title: Co-degrees resilience for perfect matchings in random hypergraphs

abstract: In this paper we prove an optimal co-degrees resilience property for the binomial $k$-uniform hypergraph model $H_{n,p}^k$ with respect to perfect matchings. That is, for a sufficiently large $n$ which is divisible by $k$, and $p\geq C_k\log_n/n$, we prove that with high probability every subgraph $H\subseteq H^k_{n,p}$ with minimum co-degree (meaning, the number of supersets every set of size $k-1$ is contained in) at least $(1/2+o(1))np$ contains a perfect matching.


title: Rydberg atom entanglements in the weak coupling regime

abstract: We present an entanglement scheme for Rydberg atoms using the van der Waals interaction phase induced by Ramsey-type pulsed interactions. This scheme realizes not only controlled phase operations between atoms at a distance larger than Rydberg blockade distance, but also various counter-intuitive entanglement examples, including two-atom entanglement in the presence of a closer third atom and $W$-state generation for partially-blockaded three atoms. Experimental realization is conducted with single rubidium atoms loaded in an array of optical tweezer dipole traps, to demonstrate the proposed entanglement generations and measurements.


title: The Omega spectrum for mod 2 KO-theory

abstract: The 8-periodic theory that comes from the KO-theory of the mod 2 Moore space is the same as the real first Morava K-theory obtained from the homotopy fixed points of the Z/(2) action on the first Morava K-theory. The first Morava K-theory, K(1), is just mod 2 KU-theory. We compute the homology Hopf algebras for the spaces in this Omega spectrum. There are a lot of maps into and out of these spaces and the spaces for KO- theory, KU-theory and the first Morava K-theory. For every one of these 98 maps (counting suspensions) there is a spectral sequence. We describe all 98 maps and spectral sequences. 48 of these maps involve our new spaces and 56 of the spectral sequences do. In addition, the maps on homotopy are all written down.


title: Galaxy evolution and radiative properties in the early Universe: multi-wavelength analysis in cosmological simulations

abstract: Recent observations have successfully detected UV or infrared flux from galaxies at the epoch of reionization. However, the origin of their radiative properties has not been fully understood yet. Combining cosmological hydrodynamic simulations and radiative transfer calculations, we present theoretical predictions of multi-wavelength radiative properties of the first galaxies at z=6-15. We find that most of the gas and dust are ejected from star-forming regions due to supernova (SN) feedback, which allows UV photons to escape. We show that the peak of SED rapidly shifts between UV and infrared wavelengths on a timescale of 100 Myr due to intermittent star formation and feedback. When dusty gas covers the star-forming regions, the galaxies become bright in the observed-frame sub-millimeter wavelengths. In addition, we find that the escape fraction of ionizing photons also changes between 1-40% at z>10. The mass fraction of HII region changes with the star formation history, resulting in the fluctuations of metal lines and Lyman-alpha line luminosities. In the starbursting phase of galaxies with the halo mass $\sim 10^{11}\,{\rm M_{\odot}}$ ($10^{12}\,{\rm M_{\odot}}$), the simulated galaxy has $L_{\rm [O\,III]} \sim 10^{42}\, (10^{43})\, {\rm erg\,s^{-1}}$, which is consistent with the observed star-forming galaxies at z>7. Our simulations suggest that deep [C II] observation with ALMA can trace the distribution of neutral gas extending over $\sim 20$ physical kpc. We also find that the luminosity ratio $L_{\rm [O\,III]}/L_{\rm [C\,II]}$ decreases with bolometric luminosity due to metal enrichment. Our simulations show that the combination of multi-wavelength observations by ALMA and JWST will be able to reveal the multi-phase ISM structure and the transition from starbursting to outflowing phases of high-z galaxies.


title: Shadow Detection for Ultrasound Images Using Unlabeled Data and Synthetic Shadows

abstract: Medical ultrasound is widely used technique for diagnosing internal organs. As common artifacts, shadows often appear in ultrasound images. Detecting such shadows is curious because they prevent accurate diagnosis. In this paper, we propose a novel shadow detection method based on auto-encoding structure. It once separates an input image into a shadow image and a content image using two decoders and combines them to reconstruct the input. To lead the network into separating the input, we inject synthetic shadows into the input and make the network to predict them as the shadow image. Since we know the rough shape of shadows as basic domain knowledge, we can generate plausible shadows. These processes are achieved by using only unlabeled data. Experiments on ultrasound images for fetal heart diagnosis shows the effectiveness of the method.


title: Coupling Effects in Multi-Stage Laser Wake-field Acceleration of Electrons

abstract: Staging laser wake-field acceleration is considered as a necessary technique for developing full-optical jitter-free electron accelerators. Splitting of the acceleration length into several technical parts with their lengths smaller than the dephasing length and with independent laser drivers allows generation of stable, reproducible acceleration fields. Temporal and spatial coupling of pre-accelerated electron bunches for their injection in the acceleration phase of a successive laser pulse wake field is the key part of the staging laser-driven acceleration. Here, characterization of the coupling is performed with dense, stable, a narrow energy band <3% and energy selectable electron beams with charges ~1.6 pC and energy ~10 MeV generated from a laser plasma cathode. Cumulative focusing of electron bunches in a low density pre-plasma, exhibiting the Budker- Bennett effect, is shown to result in the efficient injection of electrons even with a long distance between the injector and the booster in the laser pulse wake. Measured characteristics of electron beams modified by the booster wake field agree well with those obtained by multidimensional particle-in-cell simulations.


title: Graph Drawing with Morphing Partial Edges

abstract: A partial edge drawing (PED) of a graph is a variation of a node-link diagram. PED draws a link, which is a partial visual representation of an edge, and reduces visual clutter of the node-link diagram. However, more time is required to read a PED to infer undrawn parts. The authors propose a morphing edge drawing (MED), which is a PED that changes with time. In MED, links morph between partial and complete drawings; thus, a reduced load for estimation of undrawn parts in a PED is expected. Herein, a formalization of MED is shown based on a formalization of PED. Then, requirements for the scheduling of morphing are specified. The requirements inhibit morphing from crossing and shorten the overall time for morphing the edges. Moreover, an algorithm for a scheduling method implemented by the authors is illustrated and the effectiveness of PED from a reading time viewpoint is shown through an experimental evaluation.


title: Learning Compact Target-Oriented Feature Representations for Visual Tracking

abstract: Many state-of-the-art trackers usually resort to the pretrained convolutional neural network (CNN) model for correlation filtering, in which deep features could usually be redundant, noisy and less discriminative for some certain instances, and the tracking performance might thus be affected. To handle this problem, we propose a novel approach, which takes both advantages of good generalization of generative models and excellent discrimination of discriminative models, for visual tracking. In particular, we learn compact, discriminative and target-oriented feature representations using the Laplacian coding algorithm that exploits the dependence among the input local features in a discriminative correlation filter framework. The feature representations and the correlation filter are jointly learnt to enhance to each other via a fast solver which only has very slight computational burden on the tracking speed. Extensive experiments on three benchmark datasets demonstrate that this proposed framework clearly outperforms baseline trackers with a modest impact on the frame rate, and performs comparably against the state-of-the-art methods.


title: Impact of a Spinning Supermassive Black Hole on the Orbit and Gravitational Waves of a Nearby Compact Binary

abstract: Recent theoretical studies suggest that stellar-mass binary black holes (BBHs) would merge more efficiently due to the Kozai-Lidov mechanism if these binaries form in the vicinity of supermassive black holes (SMBHs). Since SMBHs are likely rotating rapidly, we continue our earlier study on the generalization of the Kozai-Lidov formalism to include the spin of the SMBH and study the evolution of a nearby BBH. We find that the eccentricity and orbital inclination of the BBH is significantly affected, because the spin (i) forces the orbital plane of the center-of-mass of the BBH around the SMBH to precess (the Lense-Thirring effect) and (ii) imposes an additional gravitomagnetic force on the BBH. As a result, the merger time of the BBH could be significantly different. We calculate the waveform from the BBH in one representative example and study its detectability by a milli-Hertz GW detector, such as the Laser Interferometer Space Antenna (LISA). We find that the signal is distinguishable from that in the case without spin. Our results imply that the BBHs in the LISA band could potentially be used to probe the spin of the SMBHs in galaxy centers.


title: Sensitivity Analysis of Treatment Effect to Unmeasured Confounding in Observational Studies with Survival and Competing Risks Outcomes

abstract: No unmeasured confounding is often assumed in estimating treatment effects in observational data when using approaches such as propensity scores and inverse probability weighting. However, in many such studies due to the limitation of the databases, collected confounders are not exhaustive, and it is crucial to examine the extent to which the resulting estimate is sensitive to the unmeasured confounders. We consider this problem for survival and competing risks data. Due to the complexity of models for such data, we adapt the simulated potential confounders approach of Carnegie et al. (2016), which provides a general tool for sensitivity analysis due to unmeasured confounding. More specifically, we specify one sensitivity parameter to quantify the association between an unmeasured confounder and the treatment assignment, and another set of parameters to quantify the association between the confounder and the time-to-event outcomes. By varying the magnitudes of the sensitivity parameters, we estimate the treatment effect of interest using the stochastic EM and the EM algorithms. We demonstrate the performance of our methods on simulated data, and apply them to a comparative effectiveness study in inflammatory bowel disease (IBD).


title: ReadsMap: a new tool for high precision mapping of DNAseq and RNAseq read sequences

abstract: There are currently plenty of programs available for mapping short sequences (reads) to a genome. Most of them, however, including such popular and actively developed programs as Bowtie, BWA, TopHat and many others, are based on Burrows-Wheeler Transform (BWT) algorithm. This approach is very effective for mapping high-homology reads, but runs into problems when mapping reads with high level of errors or SNP. Also it has problems with mapping RNASeq spliced reads (such as reads that aligning with gaps corresponding intron sequences), the kind that is essential for finding introns and alternative splicing gene isoforms. Meanwhile, finding intron positions is the most important task for determining the gene structure, and especially alternatively spliced variants of genes. In this paper, we propose a new algorithm that involves hashing reference genome. ReadsMap program, implementing such algorithm, demonstrate very high-accuracy mapping of large number of short reads to one or more genomic contigs. It is achieved mostly by better alignment of very short parts of reads separated by long introns with accounting information from mapping other reads containing the same intron inserted between bigger blocks. Availability and implementation: ReadsMap is implemented in C. It is incorporated in Fgenesh++ gene identification pipeline and is freely available to academic users at Softberry web server www.softberry.com.


title: Forecasting age distribution of death counts: An application to annuity pricing

abstract: We consider a compositional data analysis approach to forecasting the age distribution of death counts. Using the age-specific period life-table death counts in Australia obtained from the Human Mortality Database, the compositional data analysis approach produces more accurate one- to 20-step-ahead point and interval forecasts than Lee-Carter method, Hyndman-Ullah method, and two na\"{i}ve random walk methods. The improved forecast accuracy of period life-table death counts is of great interest to demographers for estimating survival probabilities and life expectancy, and to actuaries for determining temporary annuity prices for various ages and maturities. Although we focus on temporary annuity prices, we consider long-term contracts which make the annuity almost lifetime, in particular when the age at entry is sufficiently high.


title: Cross-lingual Text-independent Speaker Verification using Unsupervised Adversarial Discriminative Domain Adaptation

abstract: Speaker verification systems often degrade significantly when there is a language mismatch between training and testing data. Being able to improve cross-lingual speaker verification system using unlabeled data can greatly increase the robustness of the system and reduce human labeling costs. In this study, we introduce an unsupervised Adversarial Discriminative Domain Adaptation (ADDA) method to effectively learn an asymmetric mapping that adapts the target domain encoder to the source domain, where the target domain and source domain are speech data from different languages. ADDA, together with a popular Domain Adversarial Training (DAT) approach, are evaluated on a cross-lingual speaker verification task: the training data is in English from NIST SRE04-08, Mixer 6 and Switchboard, and the test data is in Chinese from AISHELL-I. We show that with the ADDA adaptation, Equal Error Rate (EER) of the x-vector system decreases from 9.331\% to 7.645\%, relatively 18.07\% reduction of EER, and 6.32\% reduction from DAT as well. Further data analysis of ADDA adapted speaker embedding shows that the learned speaker embeddings can perform well on speaker classification for the target domain data, and are less dependent with respect to the shift in language.


title: Characterizations of the Hardy space $\mathcal{H}_{FIO}^{1}(\mathbb{R}^{n})$ for Fourier integral operators

abstract: The Hardy spaces for Fourier integral operators $\mathcal{H}_{FIO}^{p}(\mathbb{R}^{n})$, for $1\leq p\leq \infty$, were introduced by Smith in [Smith,1998] and Hassell et al. in [Hassell-Portal-Rozendaal,2020]. In this article, we give several equivalent characterizations of $\mathcal{H}_{FIO}^{1}(\mathbb{R}^{n})$, for example in terms of Littlewood--Paley $g$ functions and maximal functions. This answers a question from [Rozendaal,2021]. We also give several applications of the characterizations.


title: Image to Video Domain Adaptation Using Web Supervision

abstract: Training deep neural networks typically requires large amounts of labeled data which may be scarce or expensive to obtain for a particular target domain. As an alternative, we can leverage webly-supervised data (i.e. results from a public search engine) which are relatively plentiful but may contain noisy results. In this work, we propose a novel two-stage approach to learn a video classifier using webly-supervised data. We argue that learning appearance features and then temporal features sequentially, rather than simultaneously, is an easier optimization for this task. We show this by first learning an image model from web images, which is used to initialize and train a video model. Our model applies domain adaptation to account for potential domain shift present between the source domain (webly-supervised data) and target domain and also accounts for noise by adding a novel attention component. We report results competitive with state-of-the-art for webly-supervised approaches on UCF-101 (while simplifying the training process) and also evaluate on Kinetics for comparison.


title: TopoTag: A Robust and Scalable Topological Fiducial Marker System

abstract: Fiducial markers have been playing an important role in augmented reality (AR), robot navigation, and general applications where the relative pose between a camera and an object is required. Here we introduce TopoTag, a robust and scalable topological fiducial marker system, which supports reliable and accurate pose estimation from a single image. TopoTag uses topological and geometrical information in marker detection to achieve higher robustness. Topological information is extensively used for 2D marker detection, and further corresponding geometrical information for ID decoding. Robust 3D pose estimation is achieved by taking advantage of all TopoTag vertices. Without sacrificing bits for higher recall and precision like previous systems, TopoTag can use full bits for ID encoding. TopoTag supports tens of thousands unique IDs and easily extends to millions of unique tags resulting in massive scalability. We collected a large test dataset including in total 169,713 images for evaluation, involving in-plane and out-of-plane rotation, image blur, different distances and various backgrounds, etc. Experiments on the dataset and real indoor and outdoor scene tests with a rolling shutter camera both show that TopoTag significantly outperforms previous fiducial marker systems in terms of various metrics, including detection accuracy, vertex jitter, pose jitter and accuracy, etc. In addition, TopoTag supports occlusion as long as the main tag topological structure is maintained and allows for flexible shape design where users can customize internal and external marker shapes. Code for our marker design/generation, marker detection, and dataset are available at http://herohuyongtao.github.io/research/publications/topo-tag/.


title: Minimal Surfaces in H^2xR

abstract: We give a fairly complete solution to the asymptotic Plateau Problem for minimal surfaces in H^2xR. In particular, we identify the collection of finite Jordan curves in the asymptotic cylinder which bounds a minimal surface in H^2xR.


title: Applying the Fokker--Planck equation to grating-based x-ray phase and dark-field imaging

abstract: X-ray imaging has conventionally relied upon attenuation to provide contrast. In recent years, two complementary modalities have been added; phase contrast and dark-field x-ray imaging, capturing weakly attenuating and sub-pixel sample structures respectively. These three modalities can be accessed using a crystal analyser, a grating interferometer or by looking at a directly-resolved grid, grating or speckle pattern. Grating and grid-based methods extract a differential phase signal by measuring how far a feature in the illumination has been shifted transversely due to the presence of a sample. The dark-field signal is extracted by measuring how the visibility of the structured illumination is decreased, typically due to the presence of sub-pixel structures in a sample. The strength of the dark-field signal may depend on the grating period, the pixel size and the set-up distances, and additional dark-field signal contributions may be seen as a result of strong phase effects or other factors. In this paper we show that the finite-difference form of the Fokker--Planck Equation can be applied to describe the drift (phase signal) and diffusion (dark-field signal) of the periodic or structured illumination used in phase contrast x-ray imaging with gratings, in order to better understand any cross-talk between attenuation, phase and dark-field x-ray signals. In future work, this mathematical description could be used as a basis for new approaches to the inverse problem of recovering both phase and dark-field information.


title: Fractional Newton-Raphson Method and Some Variants for the Solution of Non-linear Systems

abstract: The following document presents some novel numerical methods valid for one and several variables, which using the fractional derivative, allow to find solutions for some non-linear systems in the complex space using real initial conditions. The origin of these methods is the fractional Newton-Raphson method but unlike the latter, the orders of fractional derivatives proposed here are functions. In the first method, a function is used to guarantee an order of convergence (at least) quadratic, and in the others, a function is used to avoid the discontinuity that is generated when the fractional derivative of the constants is used, and with this, it is possible that the methods have at most an order of convergence (at least) linear.


title: V2S attack: building DNN-based voice conversion from automatic speaker verification

abstract: This paper presents a new voice impersonation attack using voice conversion (VC). Enrolling personal voices for automatic speaker verification (ASV) offers natural and flexible biometric authentication systems. Basically, the ASV systems do not include the users' voice data. However, if the ASV system is unexpectedly exposed and hacked by a malicious attacker, there is a risk that the attacker will use VC techniques to reproduce the enrolled user's voices. We name this the ``verification-to-synthesis (V2S) attack'' and propose VC training with the ASV and pre-trained automatic speech recognition (ASR) models and without the targeted speaker's voice data. The VC model reproduces the targeted speaker's individuality by deceiving the ASV model and restores phonetic property of an input voice by matching phonetic posteriorgrams predicted by the ASR model. The experimental evaluation compares converted voices between the proposed method that does not use the targeted speaker's voice data and the standard VC that uses the data. The experimental results demonstrate that the proposed method performs comparably to the existing VC methods that trained using a very small amount of parallel voice data.


title: The fault-tolerant cluster-sending problem

abstract: The development of fault-tolerant distributed systems that can tolerate Byzantine behavior has traditionally been focused on consensus protocols, which support fully-replicated designs. For the development of more sophisticated high-performance Byzantine distributed systems, more specialized fault-tolerant communication primitives are necessary, however. In this paper, we identify an essential communication primitive and study it in depth. In specifics, we formalize the cluster-sending problem, the problem of sending a message from one Byzantine cluster to another Byzantine cluster in a reliable manner. We not only formalize this fundamental problem, but also establish lower bounds on the complexity of this problem under crash failures and Byzantine failures. Furthermore, we develop practical cluster-sending protocols that meet these lower bounds and, hence, have optimal complexity. As such, our work provides a strong foundation for the further exploration of novel designs that address challenges encountered in fault-tolerant distributed systems.


title: A Deep Learning Approach for Tweet Classification and Rescue Scheduling for Effective Disaster Management

abstract: It is a challenging and complex task to acquire information from different regions of a disaster-affected area in a timely fashion. The extensive spread and reach of social media and networks allow people to share information in real-time. However, the processing of social media data and gathering of valuable information require a series of operations such as (1) processing each specific tweet for a text classification, (2) possible location determination of people needing help based on tweets, and (3) priority calculations of rescue tasks based on the classification of tweets. These are three primary challenges in developing an effective rescue scheduling operation using social media data. In this paper, first, we propose a deep learning model combining attention based Bi-directional Long Short-Term Memory (BLSTM) and Convolutional Neural Network (CNN) to classify the tweets under different categories. We use pre-trained crisis word vectors and global vectors for word representation (GLoVe) for capturing semantic meaning from tweets. Next, we perform feature engineering to create an auxiliary feature map which dramatically increases the model accuracy. In our experiments using real data sets from Hurricanes Harvey and Irma, it is observed that our proposed approach performs better compared to other classification methods based on Precision, Recall, F1-score, and Accuracy, and is highly effective to determine the correct priority of a tweet. Furthermore, to evaluate the effectiveness and robustness of the proposed classification model a merged dataset comprises of 4 different datasets from CrisisNLP and another 15 different disasters data from CrisisLex are used. Finally, we develop an adaptive multitask hybrid scheduling algorithm considering resource constraints to perform an effective rescue scheduling operation considering different rescue priorities.


title: Learning to Generalize to Unseen Tasks with Bilevel Optimization

abstract: Recent metric-based meta-learning approaches, which learn a metric space that generalizes well over combinatorial number of different classification tasks sampled from a task distribution, have been shown to be effective for few-shot classification tasks of unseen classes. They are often trained with episodic training where they iteratively train a common metric space that reduces distance between the class representatives and instances belonging to each class, over large number of episodes with random classes. However, this training is limited in that while the main target is the generalization to the classification of unseen classes during training, there is no explicit consideration of generalization during meta-training phase. To tackle this issue, we propose a simple yet effective meta-learning framework for metricbased approaches, which we refer to as learning to generalize (L2G), that explicitly constrains the learning on a sampled classification task to reduce the classification error on a randomly sampled unseen classification task with a bilevel optimization scheme. This explicit learning aimed toward generalization allows the model to obtain a metric that separates well between unseen classes. We validate our L2G framework on mini-ImageNet and tiered-ImageNet datasets with two base meta-learning few-shot classification models, Prototypical Networks and Relation Networks. The results show that L2G significantly improves the performance of the two methods over episodic training. Further visualization shows that L2G obtains a metric space that clusters and separates unseen classes well.


title: Revisiting consensus protocols through wait-free parallelization

abstract: The recent surge of blockchain systems has renewed the interest in traditional Byzantine fault-tolerant consensus protocols. Many such consensus protocols have a primary-backup design in which an assigned replica, the primary, is responsible for coordinating the consensus protocol. Although the primary-backup design leads to relatively simple and high performance consensus protocols, it places an unreasonable burden on a good primary and allows malicious primaries to substantially affect the system performance. In this paper, we propose a protocol-agnostic approach to improve the design of primary backup consensus protocols. At the core of our approach is a novel wait-free approach of running several instances of the underlying consensus protocol in parallel. To yield a high performance parallelized design, we present coordination-free techniques to order operations across parallel instances, deal with instance failures, and assign clients to specific instances. Consequently, the design we present is able to reduce the load on individual instances and primaries, while also reducing the adverse effects of any malicious replicas.


title: 3-uniform hypergraphs with few Berge paths of length three between any two vertices

abstract: Recently, Berge theta hypergraphs have received special attention due to the similarity with Berge even cycles. Let $r$-uniform Berge theta hypergraph $\Theta_{\ell,t}^{B}$ be the $r$-uniform hypergraph consisting of $t$ internally disjoint Berge paths of length $\ell$ with the same pair of endpoints. In this work, we determine the Tur\'{a}n number of $3$-uniform Berge theta hypergraph when $\ell=3$ and $t$ is relatively small. More precisely, we provide an explicit construction giving \begin{align*} \textup{ex}_{3}(n,\Theta_{3,217}^{B})=\Omega(n^{\frac{4}{3}}). \end{align*} This matches an earlier upper bound by He and Tait up to an absolute constant factor. The construction is algebraic, which is based on some equations over finite fields, and the parameter $t$ in our construction is much smaller than that in random algebraic construction. Our main technique is using the resultant of polynomials, which appears to be a powerful technique to eliminate variables.


title: Downlink Analysis of NOMA-enabled Cellular Networks with 3GPP-inspired User Ranking

abstract: This paper provides a comprehensive downlink analysis of non-orthogonal multiple access (NOMA) enabled cellular networks using tools from stochastic geometry. As a part of this analysis, we develop a novel 3GPP-inspired user ranking technique to construct a user cluster for the non-orthogonal transmission by grouping users from the cell center (CC) and cell edge (CE) regions. This technique allows to partition the users with distinct link qualities, which is imperative for harnessing NOMA performance gains. Our analysis is focused on the performance of a user cluster in the typical cell, which is significantly different from the standard stochastic geometry-based approach of analyzing the performance of the typical user. For this setting, we first derive the moments of the meta distributions for the CC and CE users under NOMA and orthogonal multiple access (OMA). Using this, we then derive the distributions of the transmission rates and mean packet delays under non-real time (NRT) and real-time (RT) service models, respectively, for both CC and CE users. Finally, we study two resource allocation (RA) techniques with the objective of maximizing the cell sum rate (CSR) under NRT service, and the sum effective capacity (SEC) under RT service. In addition to providing several useful design insights, our results demonstrate that NOMA provides improved rate region and higher CSR as compared to OMA. In addition, we also show that NOMA provides better SEC as compared to OMA for the higher user density.


title: Quasinormal modes of supersymmetric microstate geometries from the D1-D5 CFT

abstract: We revisit the study of the probe scalar quasinormal modes of a class of three-charge supersymmetric microstate geometries. We compute the real and imaginary parts of the quasinormal modes and show that in the parameter range when the geometries have large AdS region, the spectrum is precisely reproduced from a D1-D5 orbifold CFT analysis. The spectrum includes the slow decaying modes pointed out by Eperon, Reall, and Santos. We analyse in detail the nature of the quasinormal modes by studying the scalar wavefunction. We show that these modes correspond to slow leakage of excitation from AdS throat to infinity.


title: Quantum-enhanced least-square support vector machine: simplified quantum algorithm and sparse solutions

abstract: Quantum algorithms can enhance machine learning in different aspects. Here, we study quantum-enhanced least-square support vector machine (LS-SVM). Firstly, a novel quantum algorithm that uses continuous variable to assist matrix inversion is introduced to simplify the algorithm for quantum LS-SVM, while retaining exponential speed-up. Secondly, we propose a hybrid quantum-classical version for sparse solutions of LS-SVM. By encoding a large dataset into a quantum state, a much smaller transformed dataset can be extracted using quantum matrix toolbox, which is further processed in classical SVM. We also incorporate kernel methods into the above quantum algorithms, which uses both exponential growth Hilbert space of qubits and infinite dimensionality of continuous variable for quantum feature maps. The quantum LS-SVM exploits quantum properties to explore important themes for SVM such as sparsity and kernel methods, and stresses its quantum advantages ranging from speed-up to the potential capacity to solve classically difficult machine learning tasks.


title: Minimum Energy Analysis for Robust Gaussian Joint Source-Channel Coding with a Square-Law Profile

abstract: A distortion-noise profile is a function indicating the maximum allowed source distortion value for each noise level in the channel. In this paper, the minimum energy required to achieve a distortion noise profile is studied for Gaussian sources which are transmitted robustly over Gaussian channels. We provide improved lower and upper bounds for the minimum energy behavior of the square-law profile using a family of lower bounds and our proposed coding scheme.


title: New Venues in Formation and Detection for Primordial Black Hole Dark Matter

abstract: Primordial black holes (PBHs) are not as exotic as once thought and constitute a compelling non-particle dark matter (DM) candidate. We present a novel general PBH formation mechanism from scalar field fragmentation, which does not suffer from the inflaton potential fine-tuning that plagues many of the standard PBH formation models. We discuss how interactions of compact stars with very small sub-lunar/asteroid-size PBHs, which reside in the open window of parameter space where PBHs can constitute all of DM, allow for a slew of new astrophysical signatures that could shed light on PBH DM and are particularly interesting in the era of multi-messenger astronomy.


title: Multi-band optical variability of 3C 279 on diverse timescales

abstract: We have monitored the flat spectrum radio quasar, 3C 279, in the optical $B$, $V$, $R$ and $I$ passbands from 2018 February to 2018 July for 24 nights, with a total of 716 frames, to study flux, colour and spectral variability on diverse timescales. 3C\,279 was observed using seven different telescopes: two in India, two in Argentina, two in Bulgaria and one in Turkey to understand the nature of the source in optical regime. The source was found to be active during the whole monitoring period and displayed significant flux variations in $B$, $V$, $R$, and $I$ passbands. Variability amplitudes on intraday basis varied from 5.20\% to 17.9\%. A close inspection of variability patterns during our observation cycle reveals simultaneity among optical emissions from all passbands. During the complete monitoring period, progressive increase in the amplitude of variability with frequency was detected for our target. The amplitudes of variability in $B$, $V$, $R$ and $I$ passbands have been estimated to be 177\%, 172\%, 171\% and 158\%, respectively. Using the structure function technique, we found intraday timescales ranging from $\sim 23$ minutes to about 115 minutes. We also studied colour-magnitude relationship and found indications of mild bluer-when-brighter trend on shorter timescales. Spectral indices ranged from 2.3 to 3.0 with no clear trend on long term basis. We have also generated spectral energy distributions for 3C\,279 in optical $B$, $V$, $R$ and $I$ passbands for 17 nights. Finally, possible emission mechanisms causing variability in blazars are discussed briefly.


title: PERI: A Posit Enabled RISC-V Core

abstract: Owing to the failure of Dennard's scaling the last decade has seen a steep growth of prominent new paradigms leveraging opportunities in computer architecture. Two technologies of interest are Posit and RISC-V. Posit was introduced in mid-2017 as a viable alternative to IEEE 754-2008. Posit promises more accuracy, higher dynamic range, and fewer unused states along with simpler hardware designs as compared to IEEE 754-2008. RISC-V, on the other hand, provides a commercial-grade open-source ISA. It is not only elegant and simple but also highly extensible and customizable, thereby facilitating novel micro-architectural research and exploration. In this paper, we bring these two technologies together and propose the first Posit Enabled RISC-V core. The paper provides insights on how the current 'F' extension and the custom op-code space of RISC-V can be leveraged/modified to support Posit arithmetic. We also present implementation details of a parameterized and feature-complete Posit FPU which is integrated with the RISC-V compliant SHAKTI C-class core either as an execution unit or as an accelerator. To fully leverage the potential of Posit, we further enhance our Posit FPU, with minimal overheads, to support two different exponent sizes (with posit-size being 32-bits). This allows applications to switch from high-accuracy computation mode to a mode with higher dynamic-range at run-time. In the absence of viable software tool-chain to enable porting of applications in the Posit domain, we present a workaround on how certain applications can be modified minimally to exploit the existing RISC-V tool-chain. We also provide examples of applications which can perform better with Posit as compared to IEEE 754-2008. The proposed Posit FPU consumes 3507 slice LUTs and 1294 slice registers on an Artix-7-100T Xilinx FPGA while capable of operating at 100 MHz.


title: Dynamics of observables in a $q$-deformed harmonic oscillator

abstract: Chaos in classical systems has been studied in plenty over many years. Although the search for chaos in quantum systems has been an area of prominent research over the last few decades, the detailed analysis of many inherently chaotic quantum systems based on expectation values of dynamical variables has not been reported in the literature. In this paper, we extend the study of dynamical behaviour using expectation values of variables to a $q$-deformed harmonic oscillator. The system is found to be periodic, quasi-periodic or chaotic depending on the values of the deformation parameter $q$ and the deformed coherent amplitude $\alpha_{q}$, thus enabling us to explicitly classify the chaotic nature of the system on the basis of these parameters. The chaotic properties of the system are clearly illustrated through recurrence plots, power spectra, first-return-time distributions and Lyapunov exponents of the time series obtained for the expectation values of the dynamic variables.


title: Modular system for fluorescence-based single photon generation using retro-reflector

abstract: Apparatus for fluorescence-based single photon generation includes collection optics and various setups for characterization. Managing this system often reveals complexity in such a way that adjusting in a small region changes optimal alignments of others. We suggest here a modular system, where the optimal alignment is given to each compartment and tested independently. Based on this concept, we built a system for single photon generation with fluorescence center in hexagonal boron nitride nano-flake, advantageous for scaling up the number of single mode fiber output and a high degree of stability. The system allowed for a practical use of single photon stream extended over an hour with a uniform count rate of small fluctuation levels.