title: X-Ray Spectral Model from Clumpy Torus and Its Application to Circinus Galaxy

abstract: We construct an X-ray spectral model from the clumpy torus in an active galactic nucleus (AGN), designated as "XCLUMPY", utilizing the Monte Carlo simulation for Astrophysics and Cosmology framework (MONACO: Odaka et al. 2011, 2016). The adopted geometry of the torus is the same as that in Nenkova et al. (2008), who assume a power law distribution of clumps in the radial direction and a normal distribution in the elevation direction. We investigate the dependence of the X-ray continuum and Fe K$\alpha$ fluorescence line profile on the torus parameters. Our model is compared with other torus models: MYTorus model (Murphy & Yaqoob 2009), Ikeda model (Ikeda et al. 2009), and CTorus model (Liu & Li 2014). As an example, we also present the results applied to the broadband X-ray spectra of the Circinus galaxy observed with XMM-Newton, Suzaku, and NuSTAR. Our model can well reproduce the data, yielding a hydrogen column density along the equatorial plane $N_{\mathrm{H}}^{\mathrm{Equ}} = 9.08_{-0.08}^{+0.14} \times 10^{24}$ cm$^{-2}$, a torus angular width $\sigma = 14.7_{-0.39}^{+0.44}$ degree, and a 2--10 keV luminosity $\log L_{2-10}/\mathrm{erg \ s^{-1}} = 42.8$. These results are discussed in comparison with the observations in other wavelengths.


title: Rapid Luminosity Decline and Subsequent Reformation of the Innermost Dust Distribution in the Changing-look AGN Mrk 590

abstract: We examine the long-term optical/near-infrared (NIR) flux variability of a "changing-look" active galactic nucleus (AGN) Mrk 590 between 1998 and 2007. Multi-band multi-epoch optical/NIR photometry data from the SDSS Stripe 82 database and the Multicolor Active Galactic Nuclei Monitoring (MAGNUM) project reveal that Mrk 590 experienced a sudden luminosity decrease during the period from 2000 to 2001. Detection of dust reverberation lag signals between $V$- and $K$-band light curves obtained by the MAGNUM project during the faint state in $2003-2007$ suggests that the dust torus innermost radius $R_\text{dust}$ of Mrk 590 had become very small [$R_\text{dust} \simeq 32$~ light-days (lt-days)] by the year 2004 according to the aforementioned significant decrease in AGN luminosity. The $R_\text{dust}$ in the faint state is comparable to the H$\beta$ broad line region (BLR) radius of $R_{\text{H}\beta, \text{BLR}} \simeq 26$ lt-days measured by previous reverberation mapping observations during the bright state of Mrk 590 in $1990-1996$. These observations indicate that the innermost radius of the dust torus in Mrk 590 decreased rapidly after the AGN ultraviolet-optical luminosity drop, and that the replenishment time scale of the innermost dust distribution is less than 4 years, which is much shorter than the free fall time scale of BLR gas or dust clouds. We suggest that rapid replenishment of the innermost dust distribution can be accomplished either by new dust formation in radiatively-cooled BLR gas clouds or by new dust formation in the disk atmosphere and subsequent vertical wind from the dusty disk as a result of radiation pressure.


title: Atomic Modeling of Photoionization Fronts in Nitrogen Gas

abstract: Photoionization fronts play a dominant role in many astrophysical environments, but remain difficult to achieve in a laboratory experiment. Recent papers have suggested that experiments using a nitrogen medium held at ten atmospheres of pressure that is irradiated by a source with a radiation temperature of T$_{\rm R}\sim$ 100 eV can produce viable photoionization fronts. We present a suite of one-dimensional numerical simulations using the \helios\ multi-material radiation hydrodynamics code that models these conditions and the formation of a photoionization front. We study the effects of varying the atomic kinetics and radiative transfer model on the hydrodynamics and ionization state of the nitrogen gas, finding that more sophisticated physics, in particular a multi-angle long characteristic radiative transfer model and a collisional-radiative atomics model, dramatically changes the atomic kinetic evolution of the gas. A photoionization front is identified by computing the ratios between the photoionization rate, the electron impact ionization rate, and the total recombination rate. We find that due to the increased electron temperatures found using more advanced physics that photoionization fronts are likely to form in our nominal model. We report results of several parameter studies. In one of these, the nitrogen pressure is fixed at ten atmospheres and varies the source radiation temperature while another fixes the temperature at 100 eV and varied the nitrogen pressure. Lower nitrogen pressures increase the likelihood of generating a photoionization front while varying the peak source temperature has little effect.


title: How to measure the entropy of a mesoscopic system via thermoelectric transport

abstract: Entropy is a fundamental thermodynamic quantity indicative of the accessible degrees of freedom in a system. While it has been suggested that the entropy of a mesoscopic system can yield nontrivial information on emergence of exotic states, its measurement in such small electron-number system is a daunting task. Here we propose a method to extract the entropy of a Coulomb-blockaded mesoscopic system from transport measurements. We prove analytically and demonstrate numerically the applicability of the method to such a mesoscopic system of arbitrary spectrum and degeneracies. We then apply our procedure to measurements of thermoelectric response of a single quantum dot, and demonstrate how it can be used to deduce the entropy change across Coulomb-blockade valleys, resolving, along the way, a long-standing puzzle of the experimentally observed finite thermoelectric response at the apparent particle-hole symmetric point.


title: Testing $\Lambda$CDM With Dwarf Galaxy Morphology

abstract: The leading tensions to the collisionless cold dark matter (CDM) paradigm are the "small-scale controversies", discrepancies between observations at the dwarf-galactic scale and their simulational counterparts. In this work we consider methods to infer 3D morphological information on Local Group dwarf spheroidals, and test the fitness of CDM+hydrodynamics simulations to the observed galaxy shapes. We find that the subpopulation of dwarf galaxies with mass-to-light ratio $\gtrsim 100 M_\odot/L_\odot$ reflects an oblate morphology. This is discrepant with the dwarf galaxies with mass-to-light ratio $\lesssim 100 M_\odot/L_\odot$, which reflect prolate morphologies, and more importantly with simulations of CDM-sourced galaxies which are explicitly prolate. Although more simulations and data are called for, if evidence of oblate pressure-supported stellar distributions persists, we argue that an underlying oblate non-CDM dark matter halo may be required, and present this as motivation for future studies.


title: No Permanent Friends or Enemies: Tracking Relationships between Nations from News

abstract: Understanding the dynamics of international politics is important yet challenging for civilians. In this work, we explore unsupervised neural models to infer relations between nations from news articles. We extend existing models by incorporating shallow linguistics information and propose a new automatic evaluation metric that aligns relationship dynamics with manually annotated key events. As understanding international relations requires carefully analyzing complex relationships, we conduct in-person human evaluations with three groups of participants. Overall, humans prefer the outputs of our model and give insightful feedback that suggests future directions for human-centered models. Furthermore, our model reveals interesting regional differences in news coverage. For instance, with respect to US-China relations, Singaporean media focus more on "strengthening" and "purchasing", while US media focus more on "criticizing" and "denouncing".


title: The stochastic thin-film equation: existence of nonnegative martingale solutions

abstract: We consider the stochastic thin-film equation with colored Gaussian Stratonovich noise in one space dimension and establish the existence of nonnegative weak (martingale) solutions. The construction is based on a Trotter-Kato-type decomposition into a deterministic and a stochastic evolution, which yields an easy to implement numerical algorithm. Compared to previous work, no interface potential has to be included, the initial data and the solution can have de-wetted regions of positive measure, and the Trotter-Kato scheme allows for a simpler proof of existence than in case of It\^o noise.


title: The Vast Potential of Exoplanet Satellites for High-Energy Astrophysics

abstract: The photometric precision, monitoring baselines, and rapid, even sampling rates required by modern satellites designed for detecting the signal of transiting exoplanets are ideally suited to a large number of applications in high-energy astrophysics. I will exemplify this by discussing the results for active galactic nuclei from Kepler and summarizing other high-energy results from Kepler/K2. These precision instruments are currently underutilized for high-energy applications despite their great potential, due in part to complex systematics affecting the data. I will summarize these effects and mitigation approaches, and conclude by discussing how the recently launched Transiting Exoplanet Survey Satellite mission will differ from Kepler/K2 in ways significant to the high-energy community.


title: Interacting Edge States of Fermionic Symmetry-Protected Topological Phases in Two Dimensions

abstract: Recently, it has been found that there exist symmetry-protected topological phases of fermions, which have no realizations in non-interacting fermionic systems or bosonic models. We study the edge states of such an intrinsically interacting fermionic SPT phase in two spatial dimensions, protected by $\mathbb{Z}_4\times\mathbb{Z}_2^T$ symmetry. We model the edge Hilbert space by replacing the internal $\mathbb{Z}_4$ symmetry with a spatial translation symmetry, and design an exactly solvable Hamiltonian for the edge model. We show that at low-energy the edge can be described by a two-component Luttinger liquid, with nontrivial symmetry transformations that can only be realized in strongly interacting systems. We further demonstrate the symmetry-protected gaplessness under various perturbations, and the bulk-edge correspondence in the theory.


title: A Conditional Lower Bound on Graph Connectivity in MapReduce

abstract: MapReduce (and its open source implementation Hadoop) has become the de facto platform for processing large data sets. MapReduce offers a streamlined computational framework by interleaving sequential and parallel computation while hiding underlying system issues from the programmer. Due to the popularity of MapReduce, there have been attempts in the theoretical computer science community to understand the power and limitations of the MapReduce framework. In the most widely studied MapReduce models each machine has memory sub-linear in the input size to the problem, hence cannot see the entire input. This restriction places many limitations on algorithms that can be developed for the model; however, the current understanding of these restrictions is still limited. In this paper, our goal is to work towards understanding problems which do not admit efficient algorithms in the MapReduce model. We study the basic question of determining if a graph is connected or not. We concentrate on instances of this problem where an algorithm is to determine if a graph consists of a single cycle or two disconnected cycles. In this problem, locally every part of the graph is similar and the goal is to determine the global structure of the graph. We consider a natural class of algorithms that can store/process/transfer the information only in the form of paths and show that no randomized algorithm cannot answer the decision question in a sub-logarithmic number of rounds. Currently, there are no absolute super constant lower bounds on the number of rounds known for any problem in MapReduce. We introduce some of the first lower bounds for a natural graph problem, albeit for a restricted class of algorithms. We believe our result makes progress towards understanding the limitations of MapReduce.


title: Universal divergence of the Renyi entropy of a thinly sliced torus at the Ising fixed point

abstract: The entanglement entropy of a quantum critical system can provide new universal numbers that depend on the geometry of the entangling bipartition. We calculate a universal number called $\kappa$, which arises when a quantum critical system is embedded on a two-dimensional torus and bipartitioned into two cylinders. In the limit when one of the cylinders is a thin slice through the torus, $\kappa$ parameterizes a divergence that occurs in the entanglement entropy sub-leading to the area law. Using large-scale Monte Carlo simulations of an Ising model in 2+1 dimensions, we access the second Renyi entropy, and determine that, at the Wilson-Fisher (WF) fixed point, $\kappa_{2,\text{WF}} = 0.0174(5)$. This result is significantly different from its value for the Gaussian fixed point, known to be $\kappa_{2,\text{Gaussian}} \approx 0.0227998$.


title: Determining Phonon Mean Free Path Spectrum by Ballistic Phonon Resistance within a Nanoslot-Patterned Thin Film

abstract: At micro- to nano-scales, classical size effects in heat conduction play an important role in suppressing the thermal transport process. Such effects occur when the characteristic lengths become commensurate to the mean free paths (MFPs) of heat carriers that are mainly phonons for nonmetallic crystals. Beyond existing experimental efforts on thin films using laser-induced thermal gratings, this work provides the complete theoretical analysis for a new approach to extract the effective phonon MFP distribution for the in-plane heat conduction within a thin film or flake-like sample. In this approach, nanoslots are patterned on a suspended thin film. Phonons will transport ballistically through the neck region between adjacent nanoslots if the phonon MFPs are much longer than the neck width. The associated "ballistic thermal resistance" for varied neck dimensions can then be used to reconstruct the phonon MFP distribution within the film. The technique can be further extended to two-dimensional materials when the relaxation time approximation is reasonably accurate.


title: Mass Production of a Trigger Data Serializer ASIC for the Upgrade of the Muon Spectrometer at the ATLAS Experiment

abstract: The Trigger Data Serializer (TDS) is a custom ASIC designed for the upgrade of the innermost station of the endcap ATLAS Muon Spectrometer. It is a mixed-signal chip with two working modes that can handle up to 128 detector channels. A total of 6,000 TDS ASICs have been produced for detector operation. This paper discusses a custom automatic test platform we developed to provide quality control of the TDS ASICs. We introduce the design, test procedures, and results obtained from this TDS testing platform.


title: On singular moduli that are S-units

abstract: Recently Yu. Bilu, P. Habegger and L. K\"uhne proved that no singular modulus can be a unit in the ring of algebraic integers. In this paper we study for which sets S of prime numbers there is no singular modulus that is an S-units. Here we prove that when the set S contains only primes congruent to 1 modulo 3 then no singular modulus can be an S-unit. We then give some remarks on the general case and we study the norm factorizations of a special family of singular moduli.


title: RepGN:Object Detection with Relational Proposal Graph Network

abstract: Region based object detectors achieve the state-of-the-art performance, but few consider to model the relation of proposals. In this paper, we explore the idea of modeling the relationships among the proposals for object detection from the graph learning perspective. Specifically, we present relational proposal graph network (RepGN) which is defined on object proposals and the semantic and spatial relation modeled as the edge. By integrating our RepGN module into object detectors, the relation and context constraints will be introduced to the feature extraction of regions and bounding boxes regression and classification. Besides, we propose a novel graph-cut based pooling layer for hierarchical coarsening of the graph, which empowers the RepGN module to exploit the inter-regional correlation and scene description in a hierarchical manner. We perform extensive experiments on COCO object detection dataset and show promising results.


title: N$^{3}$LO predictions for the decay of the Higgs boson to bottom quarks

abstract: We present a fully-differential calculation of the $H\rightarrow b\overline{b}$ decay at next-to-next- to-next-to-leading order (N$^3$LO) accuracy. Our calculation considers diagrams in which the Higgs boson couples directly to the bottom quarks, i.e. the perturbative order we consider is $\mathcal{O}(\alpha_s^3y_b^2)$. In order to regulate the infrared divergences present at this order we use the Projection-to-Born technique coupled with N-jettiness slicing. After validating our methodology at next-to-next-to-leading order (NNLO) we present exclusive jet rates and differential distributions for jet observables at N3LO accuracy using the Durham jet algorithm in the Higgs rest frame.


title: $H \rightarrow b\overline{b}j$ at Next-to-Next-to-Leading Order Accuracy

abstract: We present the calculation of the decay $H \rightarrow b\overline{b}j$ at next-to-next-to-leading order (NNLO) accuracy. We consider contributions in which the Higgs boson couples directly to bottom quarks, i.e. our predictions are accurate to order $\mathcal{O}(\alpha_s^3 y_b^2)$ . We calculate the various components needed to construct the NNLO contribution, including an independent calculation of the two-loop amplitudes. We compare our results for the two-loop amplitudes to an existing calculation finding agreement. We present multiple checks on our two-loop expression using the known infrared factorization properties as the emitted gluon becomes soft or collinear. We use our results to construct a Monte Carlo implementation of $H \rightarrow b\overline{b}j$ and present jet rates and differential distributions in the Higgs rest frame using the Durham jet algorithm.


title: Constrained Restless Bandits for Dynamic Scheduling in Cyber-Physical Systems

abstract: This paper studies a class of constrained restless multi-armed bandits (CRMAB). The constraints are in the form of time varying set of actions (set of available arms). This variation can be either stochastic or semi-deterministic. Given a set of arms, a fixed number of them can be chosen to be played in each decision interval. The play of each arm yields a state dependent reward. The current states of arms are partially observable through binary feedback signals from arms that are played. The current availability of arms is fully observable. The objective is to maximize long term cumulative reward. The uncertainty about future availability of arms along with partial state information makes this objective challenging. Applications for CRMAB can be found in resource allocation in cyber-physical systems involving components with time varying availability. First, this optimization problem is analyzed using Whittle's index policy. To this end, a constrained restless single-armed bandit is studied. It is shown to admit a threshold-type optimal policy and is also indexable. An algorithm to compute Whittle's index is presented. An alternate solution method with lower complexity is also presented in the form of an online rollout policy. A detailed discussion on the complexity of both these schemes is also presented, which suggests that online rollout policy with short look ahead is simpler to implement than Whittle's index computation. Further, upper bounds on the value function are derived in order to estimate the degree of sub-optimality of various solutions. The simulation study compares the performance of Whittle's index, online rollout, myopic and modified Whittle's index policies.


title: VoteNet: A Deep Learning Label Fusion Method for Multi-Atlas Segmentation

abstract: Deep learning (DL) approaches are state-of-the-art for many medical image segmentation tasks. They offer a number of advantages: they can be trained for specific tasks, computations are fast at test time, and segmentation quality is typically high. In contrast, previously popular multi-atlas segmentation (MAS) methods are relatively slow (as they rely on costly registrations) and even though sophisticated label fusion strategies have been proposed, DL approaches generally outperform MAS. In this work, we propose a DL-based label fusion strategy (VoteNet) which locally selects a set of reliable atlases whose labels are then fused via plurality voting. Experiments on 3D brain MRI data show that by selecting a good initial atlas set MAS with VoteNet significantly outperforms a number of other label fusion strategies as well as a direct DL segmentation approach. We also provide an experimental analysis of the upper performance bound achievable by our method. While unlikely achievable in practice, this bound suggests room for further performance improvements. Lastly, to address the runtime disadvantage of standard MAS, all our results make use of a fast DL registration approach.


title: Harmonia: Near-Linear Scalability for Replicated Storage with In-Network Conflict Detection

abstract: Distributed storage employs replication to mask failures and improve availability. However, these systems typically exhibit a hard tradeoff between consistency and performance. Ensuring consistency introduces coordination overhead, and as a result the system throughput does not scale with the number of replicas. We present Harmonia, a replicated storage architecture that exploits the capability of new-generation programmable switches to obviate this tradeoff by providing near-linear scalability without sacrificing consistency. To achieve this goal, Harmonia detects read-write conflicts in the network, which enables any replica to serve reads for objects with no pending writes. Harmonia implements this functionality at line rate, thus imposing no performance overhead. We have implemented a prototype of Harmonia on a cluster of commodity servers connected by a Barefoot Tofino switch, and have integrated it with Redis. We demonstrate the generality of our approach by supporting a variety of replication protocols, including primary-backup, chain replication, Viewstamped Replication, and NOPaxos. Experimental results show that Harmonia improves the throughput of these protocols by up to 10X for a replication factor of 10, providing near-linear scalability up to the limit of our testbed.


title: Quark and Polyakov-loop correlations in effective models at zero and non-vanishing density

abstract: The aim of this work is to shed light on some lesser known aspects of Polyakov-loop--extended chiral models (namely the Polyakov-loop extended Nambu--Jona-Lasinio and Quark-Meson models), especially on the correlation of the quark sector with the Polyakov loop. We show that the ordering of chiral and Polyakov-loop transitions and their difference in temperature as seen in lattice QCD calculations could be realized with a critical scale of the Polyakov-loop potential that is larger than the one in pure gauge theory. The comparison of the results for the Polyakov-loop susceptibility obtained using the self-consistent medium-dependent quark mass with those obtained while keeping these masses at a fixed value allows to disentangle chiral-symmetry restoration and center-symmetry breaking effects. Furthermore, a confined chirally restored phase is identified by a plateau in the quark contribution to thermodynamics and by sigma and pion spectral functions that coincide but have a small width. We also discuss that, for some large chemical potential values, the explicit center-symmetry breaking is so strong that statistical deconfinement is realized at infinitely small temperatures. Both the missing sensitivity of the Polyakov loop to the quark mass, except at close to the chiral transition, and the Polyakov loop being zero at zero temperature at all chemical potentials, can be interpreted as indications of a missing mechanism which accounts for the quark back-reaction on the Polyakov loop.


title: Non-Stationary Polar Codes for Resistive Memories

abstract: Resistive memories are considered a promising memory technology enabling high storage densities with in-memory computing capabilities. However, the readout reliability of resistive memories is impaired due to the inevitable existence of wire resistance, resulting in the sneak path problem. Motivated by this problem, we study polar coding over channels with different reliability levels, termed non-stationary polar codes, and we propose a technique improving its bit error rate (BER) performance. We then apply the framework of non-stationary polar codes to the crossbar array and evaluate its BER performance under two modeling approaches, namely binary symmetric channels (BSCs) and binary asymmetric channels (BSCs). Finally, we propose a technique for biasing the proportion of high-resistance states in the crossbar array and show its advantage in reducing further the BER. Several simulations are carried out using a SPICE-like simulator, exhibiting significant reduction in BER.


title: Stochastically modeled weakly reversible reaction networks with a single linkage class

abstract: It has been known for nearly a decade that deterministically modeled reaction networks that are weakly reversible and consist of a single linkage class have trajectories that are bounded from both above and below by positive constants (so long as the initial condition has strictly positive components). It is conjectured that the stochastically modeled analogs of these systems are positive recurrent. We prove this conjecture in the affirmative under the following additional assumptions: (i) the system is binary and (ii) for each species, there is a complex (vertex in the associated reaction diagram) that is a multiple of that species. To show this result, a new proof technique is developed in which we study the recurrence properties of the $n-$step embedded discrete time Markov chain.


title: Photon-Dark Photon Conversions in Background Electromagnetic Fields

abstract: The mixing of photons with light pseudoscalars in the presence of external electromagnetic fields has been used extensively to search for axion-like-particles. A similar effect for dark photon propagating states is usually not considered due to the Landau-Yang theorem. We point out that mixing between photon and dark photon propagating states in background electromagnetic fields can indeed occur, in non-linear QED, through a four-photon vertex by integrating out the electron box diagram. Starting from the Schwinger Lagrangian, we derive the equations of motion for dark photons interacting with the Standard Model photon through gauge kinetic terms. We provide expressions for the perpendicular and parallel refractive indices in series expansions in the critical field strength, valid both in the strong and weak background field limits. We then consider mixing between the photon-dark photon propagating system in the presence of pure electric and magnetic background fields, and work out the probability of conversion when the background fields are homogeneous. We indicate outlines of the calculation in the inhomogeneous case, and finally express our results in the active-sterile basis, where we find that the mixing induced by background fields can lead to corrections to the tree-level mixing in the zero field limit that is usually considered to probe such systems. Our results may find applications for probing photon-dark photon conversions in the vicinity of neutron stars and in table-top petawatt laser experiments.


title: Scattered compact sets in continuous images of \v{C}ech-complete spaces

abstract: Assume hat a functionally Hausdorff space $X$ is a continuous image of a \v{C}ech complete space $P$ with Lindel\"of number $l(P)<\mathfrak c$. Then the following conditions are equivalent: (i) every compact subset of $X$ is scattered, (ii) for every continuous map $f:X\to Y$ to a functionally Hausdorff space $Y$ the image $f(X)$ has cardinality $|f(X)|\le \max\{l(P),\psi(Y)\}$, (iii) no continuous map $f:X\to[0,1]$ is surjective. Also we prove the equivalence of the conditions: (a) $\omega_1<\mathfrak b$, (b) a K-analytic space $X$ (with a unique non-isolated point) is countable if and only if every compact subset of $X$ is countable.


title: Rationally elliptic toric varieties

abstract: We give a characterization of all complete smooth toric varieties whose rational homotopy is of elliptic type. All such toric varieties of complex dimension not more than three are explicitly described.


title: On Acoustic Modeling for Broadband Beamforming

abstract: In this work, we describe limitations of the free-field propagation model for designing broadband beamformers for microphone arrays on a rigid surface. Towards this goal, we describe a general framework for quantifying the microphone array performance in a general wave-field by directly solving the acoustic wave equation. The model utilizes Finite-Element-Method (FEM) for evaluating the response of the microphone array surface to background 3D planar and spherical waves. The effectiveness of the framework is established by designing and evaluating a representative broadband beamformer under realistic acoustic conditions.


title: Intentional Computational Level Design

abstract: The procedural generation of levels and content in video games is a challenging AI problem. Often such generation relies on an intelligent way of evaluating the content being generated so that constraints are satisfied and/or objectives maximized. In this work, we address the problem of creating levels that are not only playable but also revolve around specific mechanics in the game. We use constrained evolutionary algorithms and quality-diversity algorithms to generate small sections of Super Mario Bros levels called scenes, using three different simulation approaches: Limited Agents, Punishing Model, and Mechanics Dimensions. All three approaches are able to create scenes that give opportunity for a player to encounter or use targeted mechanics with different properties. We conclude by discussing the advantages and disadvantages of each approach and compare them to each other.


title: The $x_i$-eigenvalue problem on some new fuzzy spheres

abstract: We study the eigenvalue equation for the 'Cartesian coordinates' observables $x_i$ on the fully $O(2)$-covariant fuzzy circle $\{S^1_\Lambda\}_{\Lambda\in\mathbb{N}}$ ($i=1,2$) and on the fully $O(3)$-covariant fuzzy 2-sphere $\{S^2_\Lambda\}_{\Lambda\in\mathbb{N}}$ ($i=1,2,3$) introduced in [G. Fiore, F. Pisacane, J. Geom. Phys. 132 (2018), 423-451]. We show that the spectrum and eigenvectors of $x_i$ fulfill a number of properties which are expected for $x_i$ to approximate well the corresponding coordinate operator of a quantum particle forced to stay on the unit sphere.


title: Measurement of a $^7$Li tune-out wavelength by phase-patterned atom interferometry

abstract: Atom interferometers typically use the total populations the interferometer's output ports as the signal, but finer spatial structure can contain useful information. We pattern a matter-wave phase profile onto an atomic sample. An interferometer translates the phase into a measurable pattern in the atomic density that we use perform the first direct precision measurement of the $^7$Li tune-out wavelength near 671 nm. Expressed as a detuning from the $|2S_{1/2},F{=}2\rangle{\rightarrow}|2P_{1/2},F'{=}2\rangle$ transition, we find 3329.3(1.4) MHz for the tensor-shifted tune out of the $|2S_{1/2},F{=}2,m_F{=}0\rangle$ state with $\sigma^\pm$ light polarization and 3310.1(4.9) MHz for the tune out of the the scalar polarizability. This technique may be generalized for directly sensing spatially varying phase profiles.


title: The physics of (ir)rational choice

abstract: Even though classic theories and models of discrete choice pose man as a rational being, it has been shown extensively that people persistently violate rationality in their actual choices. Recent models of decision-making take these violations often (partially) into account, however, a unified framework has not been established. Here we propose such a framework, inspired by the Ising model from statistical physics, and show that representing choice problems as a graph, together with a simple choice process, allows us to explain both rational decisions as well as violations of rationality.


title: Search for sub-solar mass ultracompact binaries in Advanced LIGO's second observing run

abstract: We present an Advanced LIGO and Advanced Virgo search for sub-solar mass ultracompact objects in data obtained during Advanced LIGO's second observing run. In contrast to a previous search of Advanced LIGO data from the first observing run, this search includes the effects of component spin on the gravitational waveform. We identify no viable gravitational wave candidates consistent with sub-solar mass ultracompact binaries with at least one component between 0.2 - 1.0 solar masses. We use the null result to constrain the binary merger rate of (0.2 solar mass, 0.2 solar mass) binaries to be less than 3.7 x 10^5 Gpc^-3 yr^-1 and the binary merger rate of (1.0 solar mass, 1.0 solar mass) binaries to be less than 5.2 x 10^3 Gpc^-3 yr^-1. Sub-solar mass ultracompact objects are not expected to form via known stellar evolution channels, though it has been suggested that primordial density fluctuations or particle dark matter with cooling mechanisms and/or nuclear interactions could form black holes with sub-solar masses. Assuming a particular primordial black hole formation model, we constrain a population of merging 0.2 solar mass black holes to account for less than 16% of the dark matter density and a population of merging 1.0 solar mass black holes to account for less than 2% of the dark matter density. We discuss how constraints on the merger rate and dark matter fraction may be extended to arbitrary black hole population models that predict sub-solar mass binaries.


title: Consistent SDNs through Network State Fuzzing

abstract: The conventional wisdom is that a software-defined network (SDN) operates under the premise that the logically centralized control plane has an accurate representation of the actual data plane state. Unfortunately, bugs, misconfigurations, faults or attacks can introduce inconsistencies that undermine correct operation. Previous work in this area, however, lacks a holistic methodology to tackle this problem and thus, addresses only certain parts of the problem. Yet, the consistency of the overall system is only as good as its least consistent part. Motivated by an analogy of network consistency checking with program testing, we propose to add active probe-based network state fuzzing to our consistency check repertoire. Hereby, our system, PAZZ, combines production traffic with active probes to periodically test if the actual forwarding path and decision elements (on the data plane) correspond to the expected ones (on the control plane). Our insight is that active traffic covers the inconsistency cases beyond the ones identified by passive traffic. PAZZ prototype was built and evaluated on topologies of varying scale and complexity. Our results show that PAZZ requires minimal network resources to detect persistent data plane faults through fuzzing and localize them quickly while outperforming baseline approaches.


title: Safety-margin-based design and redesign considering mixed epistemic model uncertainty and aleatory parameter uncertainty

abstract: At the initial design stage engineers often rely on low-fidelity models that have high epistemic uncertainty. Traditional safety-margin-based deterministic design resorts to testing (e.g. prototype experiment, evaluation of high-fidelity simulation, etc.) to reduce epistemic uncertainty and achieve targeted levels of safety. Testing is used to calibrate models and prescribe redesign when tests are not passed. After calibration, reduced epistemic model uncertainty can be leveraged through redesign to restore safety or improve design performance; however, redesign may be associated with substantial costs or delays. In this paper, a methodology is described for optimizing the safety-margin-based design, testing, and redesign process to allow the designer to tradeoff between the risk of future redesign and the possible performance and reliability benefits. The proposed methodology represents the epistemic model uncertainty with a Kriging surrogate and is applicable in a wide range of design problems. The method is illustrated on a cantilever beam bending example and then a sounding rocket example. It is shown that redesign acts as a type of quality control measure to prevent an undesirable initial design from being accepted as the final design. It is found that the optimal design/redesign strategy for maximizing expected design performance includes not only redesign to correct an initial design that is later revealed to be unsafe, but also redesign to improve performance when an initial design is later revealed to be too conservative (e.g. too heavy).


title: Sub-GeV Dark Matter Searches and Electric Field Studies for the LUX and LZ Experiments

abstract: The nature of dark matter (DM) remains a mystery since it has so far eluded detection in the laboratory. To that end, the Large Underground Xenon (LUX) experiment was built to directly observe the interaction of DM with xenon target nuclei. LUX acquired data from April 2013 to May 2016 at SURF in South Dakota, which led to publications of many world-leading exclusion limits that probe much of the unexplored DM parameter space. This manuscript describes two novel direct detection methods that used the first LUX dataset to place limits on sub-GeV DM. The Bremsstrahlung and Migdal effects consider electron recoils that accompany the standard DM-nucleus scattering, thereby extending the reach of the LUX detector to lower DM masses. The spin-independent DM-nucleon scattering was constrained for four different classes of mediators for DM particles with masses of 0.4-5 GeV/c$^{2}$. The detector conditions changed significantly before its final 332 live-days of data acquisition. The electric fields varied in a non-trivial non-symmetric manner, which triggered a need for a fully 3D model of the electric fields inside the LUX detector. The successful modeling of these electric fields, described herein, enabled a thorough understanding of the detector throughout its scientific program and strengthened its sensitivity to DM. The LUX-ZEPLIN (LZ) experiment is a next-generation xenon detector soon to start searching for DM. However, increasingly large noble liquid detectors like LZ are facing challenges with applications of high voltage (HV). The Xenon Breakdown Apparatus (XeBrA) at the Lawrence Berkeley National Laboratory was built to characterize the HV behavior of liquid xenon and liquid argon. Results from XeBrA will serve not only to improve our understanding of the physical processes involved in the breakdown but also to inform the future of noble liquid detector engineering.


title: Exploring the Limitations of Behavior Cloning for Autonomous Driving

abstract: Driving requires reacting to a wide variety of complex environment conditions and agent behaviors. Explicitly modeling each possible scenario is unrealistic. In contrast, imitation learning can, in theory, leverage data from large fleets of human-driven cars. Behavior cloning in particular has been successfully used to learn simple visuomotor policies end-to-end, but scaling to the full spectrum of driving behaviors remains an unsolved problem. In this paper, we propose a new benchmark to experimentally investigate the scalability and limitations of behavior cloning. We show that behavior cloning leads to state-of-the-art results, including in unseen environments, executing complex lateral and longitudinal maneuvers without these reactions being explicitly programmed. However, we confirm well-known limitations (due to dataset bias and overfitting), new generalization issues (due to dynamic objects and the lack of a causal model), and training instability requiring further research before behavior cloning can graduate to real-world driving. The code of the studied behavior cloning approaches can be found at https://github.com/felipecode/coiltraine .


title: Tuning Topological Superconductivity in Phase-Controlled Josephson Junctions with Rashba and Dresselhaus Spin-Orbit Coupling

abstract: Recently, topological superconductors based on Josephson junctions in two-dimensional electron gases with strong Rashba spin-orbit coupling have been proposed as attractive alternatives to wire-based setups. Here, we elucidate how phase-controlled Josephson junctions based on quantum wells with [001] growth direction and an arbitrary combination of Rashba and Dresselhaus spin-orbit coupling can also host Majorana bound states for a wide range of parameters as long as the magnetic field is oriented appropriately. Hence, Majorana bound states based on Josephson junctions can appear in a wide class of two-dimensional electron gases. We study the effect of spin-orbit coupling, the Zeeman energies, and the superconducting phase difference to create a full topological phase diagram and find the optimal stability region to observe Majorana bound states in narrow junctions. Surprisingly, for equal Rashba and Dresselhaus spin-orbit coupling, well localized Majorana bound states can appear only for phase differences $\phi\neq\pi$ as the topological gap protecting the Majorana bound states vanishes at $\phi=\pi$. Our results show that the ratio between Rashba and Dresselhaus spin-orbit coupling or the choice of the in-plane crystallographic axis along which the superconducting phase bias is applied offer additional tunable knobs to test Majorana bound states in these systems. Finally, we discuss signatures of Majorana bound states that could be probed experimentally by tunneling conductance measurements at the edge of the junction.


title: Observations of Solar Coronal Rain in Null Point Topologies

abstract: Coronal rain is the well-known phenomenon in which hot plasma high in the Sun's corona undergoes rapid cooling (from > 10^6 K to < 10^4 K), condenses, and falls to the surface. Coronal rain appears frequently in active region coronal loops and is very common in post-flare loops. This Letter presents discovery observations, which show that coronal rain is ubiquitous in the embedded bipole very near a coronal hole boundary. Our observed structures formed when the photospheric decay of active region leading sunspots resulted in a large parasitic polarity embedded in a background unipolar region. We observe coronal rain to appear within the legs of closed loops well under the fan surface, as well as preferentially near separatrices of the resulting coronal topology: the spine lines, null point, and fan surface. We analyze 3 events using SDO Atmospheric Imaging Assembly (AIA) observations in the 304, 171, and 211 {/AA} channels, as well as SDO Helioseismic and Magnetic Imager (HMI) magnetograms. The frequency of rain formation and the ease with which it is observed strongly suggests that this phenomenon is generally present in null-point topologies of this size scale. We argue that these rain events could be explained by the classic process of thermal nonequilibrium or via interchange reconnection at the null; it is also possible that both mechanisms are present. Further studies with higher spatial resolution data and MHD simulations will be required to determine the exact mechanism(s).


title: TTS Skins: Speaker Conversion via ASR

abstract: We present a fully convolutional wav-to-wav network for converting between speakers' voices, without relying on text. Our network is based on an encoder-decoder architecture, where the encoder is pre-trained for the task of Automatic Speech Recognition, and a multi-speaker waveform decoder is trained to reconstruct the original signal in an autoregressive manner. We train the network on narrated audiobooks, and demonstrate multi-voice TTS in those voices, by converting the voice of a TTS robot.


title: Some Aspect of Certain two Subclass of Analytic Functions with Negative Coefficients Defined by Rafid Operator

abstract: In this paper, we define the subclasses $R_{\mu,p}^{\delta}(\alpha;A,B)\ $ and $ P_{\mu,p}^{\delta}(\alpha;A,B)\ $ of analytic functions in the open unit disc of complex plain. Then the neighborhood properties, integral means inequalities and some results concerning the partial sums of the functions belonging to these two subclasses are discussed.


title: Inviscid limit of the active interface equations

abstract: We present a detailed solution of the active interface equations in the inviscid limit. The active interface equations were previously introduced as a toy model of membrane-protein systems: they describe a stochastic interface where growth is stimulated by inclusions which themselves move on the interface. In the inviscid limit, the equations reduce to a pair of coupled conservation laws. After discussing how the inviscid limit is obtained, we turn to the corresponding Riemann problem: the solution of the set of conservation laws with discontinuous initial condition. In particular, by considering two physically meaningful initial conditions, a giant trough and a giant peak in the interface, we elucidate the generation of shock waves and rarefaction fans in the system. Then, by combining several Riemann problems, we construct an oscillating solution of the active interface with periodic boundaries conditions. The existence of this oscillating state reflects the reciprocal coupling between the two conserved quantities in our system.


title: FPGA-accelerated machine learning inference as a service for particle physics computing

abstract: New heterogeneous computing paradigms on dedicated hardware with increased parallelization, such as Field Programmable Gate Arrays (FPGAs), offer exciting solutions with large potential gains. The growing applications of machine learning algorithms in particle physics for simulation, reconstruction, and analysis are naturally deployed on such platforms. We demonstrate that the acceleration of machine learning inference as a web service represents a heterogeneous computing solution for particle physics experiments that potentially requires minimal modification to the current computing model. As examples, we retrain the ResNet-50 convolutional neural network to demonstrate state-of-the-art performance for top quark jet tagging at the LHC and apply a ResNet-50 model with transfer learning for neutrino event classification. Using Project Brainwave by Microsoft to accelerate the ResNet-50 image classification model, we achieve average inference times of 60 (10) milliseconds with our experimental physics software framework using Brainwave as a cloud (edge or on-premises) service, representing an improvement by a factor of approximately 30 (175) in model inference latency over traditional CPU inference in current experimental hardware. A single FPGA service accessed by many CPUs achieves a throughput of 600--700 inferences per second using an image batch of one, comparable to large batch-size GPU throughput and significantly better than small batch-size GPU throughput. Deployed as an edge or cloud service for the particle physics computing model, coprocessor accelerators can have a higher duty cycle and are potentially much more cost-effective.

categories: physics.data-an hep-ex physics.comp-ph physics.ins-det


title: Fast state and trap rotation of a particle in an anisotropic potential

abstract: We study the dynamics of a quantum or classical particle in a two-dimensional rotating anisotropic harmonic potential. By a sequence of symplectic transformations for constant rotation velocity we find uncoupled normal generalized coordinates and conjugate momenta in which the Hamiltonian takes the form of two independent harmonic oscillators. The decomposition into normal-mode dynamics enables us to design fast trap-rotation processes to produce a rotated version of an arbitrary initial state, when the two normal frequencies are commensurate.


title: HEPCloud, an Elastic Hybrid HEP Facility using an Intelligent Decision Support System

abstract: HEPCloud is rapidly becoming the primary system for provisioning compute resources for all Fermilab-affiliated experiments. In order to reliably meet the peak demands of the next generation of High Energy Physics experiments, Fermilab must plan to elastically expand its computational capabilities to cover the forecasted need. Commercial cloud and allocation-based High Performance Computing (HPC) resources both have explicit and implicit costs that must be considered when deciding when to provision these resources, and at which scale. In order to support such provisioning in a manner consistent with organizational business rules and budget constraints, we have developed a modular intelligent decision support system (IDSS) to aid in the automatic provisioning of resources spanning multiple cloud providers, multiple HPC centers, and grid computing federations. In this paper, we discuss the goals and architecture of the HEPCloud Facility, the architecture of the IDSS, and our early experience in using the IDSS for automated facility expansion both at Fermi and Brookhaven National Laboratory.


title: Classifications of $\ell$-Zero-Sumfree Sets

abstract: The set of all $\ell$-zero-sumfree subsets of $\mathbb{Z}/n\mathbb{Z}$ is a simplicial complex denoted by $\Delta_{n,\ell}$ We create an algorithm via defining a set of integer partitions we call $(n,\ell)$-congruent partitions in order to compute this complex for moderately-sized parameters $n$ and $\ell$. We also theoretically determine $\Delta_{n,\ell}$ for several infinite families of parameters, and compute the intersection lattices and the characteristic polynomials of the corresponding coordinate subspace arrangements.


title: End-to-End Environmental Sound Classification using a 1D Convolutional Neural Network

abstract: In this paper, we present an end-to-end approach for environmental sound classification based on a 1D Convolution Neural Network (CNN) that learns a representation directly from the audio signal. Several convolutional layers are used to capture the signal's fine time structure and learn diverse filters that are relevant to the classification task. The proposed approach can deal with audio signals of any length as it splits the signal into overlapped frames using a sliding window. Different architectures considering several input sizes are evaluated, including the initialization of the first convolutional layer with a Gammatone filterbank that models the human auditory filter response in the cochlea. The performance of the proposed end-to-end approach in classifying environmental sounds was assessed on the UrbanSound8k dataset and the experimental results have shown that it achieves 89% of mean accuracy. Therefore, the propose approach outperforms most of the state-of-the-art approaches that use handcrafted features or 2D representations as input. Furthermore, the proposed approach has a small number of parameters compared to other architectures found in the literature, which reduces the amount of data required for training.


title: Physical Symmetries Embedded in Neural Networks

abstract: Neural networks are a central technique in machine learning. Recent years have seen a wave of interest in applying neural networks to physical systems for which the governing dynamics are known and expressed through differential equations. Two fundamental challenges facing the development of neural networks in physics applications is their lack of interpretability and their physics-agnostic design. The focus of the present work is to embed physical constraints into the structure of the neural network to address the second fundamental challenge. By constraining tunable parameters (such as weights and biases) and adding special layers to the network, the desired constraints are guaranteed to be satisfied without the need for explicit regularization terms. This is demonstrated on upervised and unsupervised networks for two basic symmetries: even/odd symmetry of a function and energy conservation. In the supervised case, the network with embedded constraints is shown to perform well on regression problems while simultaneously obeying the desired constraints whereas a traditional network fits the data but violates the underlying constraints. Finally, a new unsupervised neural network is proposed that guarantees energy conservation through an embedded symplectic structure. The symplectic neural network is used to solve a system of energy-conserving differential equations and out-performs an unsupervised, non-symplectic neural network.


title: Quantum-Assisted Clustering Algorithms for NISQ-Era Devices

abstract: In the NISQ-era of quantum computing, we should not expect to see quantum devices that provide an exponential improvement in runtime for practical problems, due to the lack of error correction and small number of qubits available. Nevertheless, these devices should be able to provide other performance improvements, particularly when combined with existing classical machines. In this article, we develop several hybrid quantum-classical clustering algorithms that can be employed as subroutines on small, NISQ-era devices. These new hybrid algorithms require a number of qubits that is at most logarithmic in the size of the data, provide performance improvement and/or runtime improvement over their classical counterparts, and do not require a black-box oracle. Consequently, we are able to provide a promising near-term application of NISQ-era devices.


title: Playgol: learning programs through play

abstract: Children learn though play. We introduce the analogous idea of learning programs through play. In this approach, a program induction system (the learner) is given a set of tasks and initial background knowledge. Before solving the tasks, the learner enters an unsupervised playing stage where it creates its own tasks to solve, tries to solve them, and saves any solutions (programs) to the background knowledge. After the playing stage is finished, the learner enters the supervised building stage where it tries to solve the user-supplied tasks and can reuse solutions learnt whilst playing. The idea is that playing allows the learner to discover reusable general programs on its own which can then help solve the user-supplied tasks. We claim that playing can improve learning performance. We show that playing can reduce the textual complexity of target concepts which in turn reduces the sample complexity of a learner. We implement our idea in Playgol, a new inductive logic programming system. We experimentally test our claim on two domains: robot planning and real-world string transformations. Our experimental results suggest that playing can substantially improve learning performance. We think that the idea of playing (or, more verbosely, unsupervised bootstrapping for supervised program induction) is an important contribution to the problem of developing program induction approaches that self-discover BK.


title: From GAN to WGAN

abstract: This paper explains the math behind a generative adversarial network (GAN) model and why it is hard to be trained. Wasserstein GAN is intended to improve GANs' training by adopting a smooth metric for measuring the distance between two probability distributions.


title: Latest results from RHIC + Progress on determining $\hat{q}L$ in RHI collisions using di-hadron correlations

abstract: Results from Relativistic Heavy Ion Collider Physics in 2018 and plans for the future at Brookhaven National Laboratory are presented.


title: Quantifying the search for solid Li-ion electrolyte materials by anion: a data-driven perspective

abstract: We compile data and machine learned models of solid Li-ion electrolyte performance to assess the state of materials discovery efforts and build new insights for future efforts. Candidate electrolyte materials must satisfy several requirements, chief among them fast ionic conductivity and robust electrochemical stability. Considering these two requirements, we find new evidence to suggest that optimization of the sulfides for fast ionic conductivity and wide electrochemical stability may be more likely than optimization of the oxides, and that the oft-overlooked chlorides and bromides may be particularly promising families for Li-ion electrolytes. We also find that the nitrides and phosphides appear to be the most promising material families for electrolytes stable against Li-metal anodes. Furthermore, the spread of the existing data in performance space suggests that fast conducting materials that are stable against both Li metal and a >4V cathode are exceedingly rare, and that a multiple-electrolyte architecture is a more likely path to successfully realizing a solid-state Li metal battery by approximately an order of magnitude or more. Our model is validated by its reproduction of well-known trends that have emerged from the limited existing data in recent years, namely that the electronegativity of the lattice anion correlates with ionic conductivity and electrochemical stability. In this work, we leverage the existing data to make solid electrolyte performance trends quantitative for the first time, building a roadmap to complement material discovery efforts around desired material performance.


title: Theory of capacities in fracional Sobolev spaces with variable exponents

abstract: In this paper we develop a capacities theory connected with the fractional Sobolev spaces with variable exponents. Two kinds of capacities are studied: Sobolev capacity and relative capacity. Basic properties of capacities, including monotonicity, outer capacity and several results, are studies. We prove that both capacities is a Choquet capacity and all borel sets are capacitable.


title: Distribution LMP-based Transactive Day-ahead Market with Variable Renewable Generation

abstract: The largescale penetration of variable renewable energy (VRE) and their generation uncertainties poses a major challenge for the distribution system operator (DSO) to efficiently determine the day-ahead real and reactive power distribution locational marginal prices (DLMPs) and their underlying components. In this paper, we propose a DLMP-based transactive day-ahead market (DAM) model, that in addition to energy and losses, determines prices for creating congestions and voltage violations under peak-load and large-scale stochastic VRE penetration conditions. To account for the VRE uncertainties and the effect of their large-scale penetration on the DLMP components and distributed energy resources' (DERs) schedules, we propose a novel data-driven probability efficient point (PEP) method that computes the optimal total VRE generation at different confidence (risk) levels to incorporate in the proposed transactive DAM model. We perform a wide range of simulation studies on a modified IEEE 69-node system to validate the proposed methods and demonstrate the effect of peak load conditions, large-scale VRE penetration, and inclusion of battery energy storage systems (BESS) on the resulting positive or negative real and reactive power DLMPs and their components.


title: Microwave-to-optical conversion via four-wave-mixing in a cold ytterbium ensemble

abstract: Interfacing superconducting qubits with optical photons requires noise-free microwave-to-optical transducers, a technology currently not realized at the single-photon level. We propose to use four-wave-mixing in an ensemble of cold ytterbium (Yb) atoms prepared in the metastable 'clock' state. The parametric process uses two high-lying Rydberg states for bidirectional conversion between a 10 GHz microwave photon and an optical photon in the telecommunication E-band. To avoid noise photons due to spontaneous emission, we consider continuous operation far detuned from the intermediate states. We use an input-output formalism to predict conversion efficiencies of $\approx50\%$ with bandwidths of $\approx100$ kHz.


title: Chiral perturbation theory in the environment with chiral imbalance

abstract: We discuss the description of chiral imbalance imprint in pion matter. It is based on the Chiral Perturbation Theory supplemented by medium induced chiral chemical potential in the covariant derivative. Such an implementation of chiral imbalance recently was explored in linear sigma model inspired by QCD. The relationship between two sigma models is examined. A possible experimental detection of chiral imbalance in charged pion decays inside the fireball is pointed out.


title: Equivariant mappings and invariant sets on Minkowski space

abstract: In this paper we introduce the systematic study of invariant functions and equivariant mappings defined on Minkowski space under the action of the Lorentz group. We adapt some known results from the orthogonal group acting on the Euclidean space to the Lorentz group acting on the Minkowski space. In addition, an algorithm is given to compute generators of the ring of functions that are invariant under an important class of Lorentz subgroups, namely when these are generated by involutions, which is also useful to compute equivariants. Furthermore, general results on invariant subspaces of the Minkowski space are presented, with a characterization of invariant lines and planes in the two lowest dimensions.


title: Landmark Proportional Subdistribution Hazards Models for Dynamic Prediction of Cumulative Incidence Functions

abstract: An individualized risk prediction model that dynamically updates the probability of a clinical event from a specific cause is valuable for physicians to be able to optimize personalized treatment strategies in real-time by incorporating all available information collected over the follow-up. However, this is more complex and challenging when competing risks are present, because it requires simultaneously updating the overall survival and the cumulative incidence functions (CIFs) while adjusting for the time-dependent covariates and time-varying covariate effects. In this study, we developed a landmark proportional subdistribution hazards (PSH) model and a more comprehensive supermodel by extending the landmark method to the Fine-Gray model. The performance of our models was assessed via simulations and through analysis of data from a multicenter clinical trial for breast cancer patients. Our proposed models have appealing advantages over other dynamic prediction models for data with competing risks. First, our models are robust against violations of the PSH assumption and can directly predict the conditional CIFs bypassing the estimation of overall survival and greatly simplify the prediction procedure. Second, our landmark PSH supermodel enables researchers to make predictions at a set of landmark points in one step. Third, the proposed models can easily incorporate various types of time-dependent information using existing standard software without computational burden.


title: Gravitational probes of ultra-light axions

abstract: The axion is a hypothetical, well-motivated dark-matter particle whose existence would explain the lack of charge-parity violation in the strong interaction. In addition to this original motivation, an `axiverse' of ultra-light axions (ULAs) with masses $10^{-33}\,{\rm eV}\lesssim m_{\rm a}\lesssim 10^{-10}\,{\rm eV}$ also emerges from string theory. Depending on the mass, such a ULA contributes to the dark-matter density, or alternatively, behaves like dark energy. At these masses, ULAs' classical wave-like properties are astronomically manifested, potentially mitigating observational tensions within the $\Lambda$CDM paradigm on local-group scales. ULAs also provide signatures on small scales such as suppression of structure, interference patterns and solitons to distinguish them from heavier dark matter candidates. Through their gravitational imprint, ULAs in the presently allowed parameter space furnish a host of observational tests to target in the next decade, altering standard predictions for microwave background anisotropies, galaxy clustering, Lyman-$\alpha$ absorption by neutral hydrogen along quasar sightlines, pulsar timing, and the black-hole mass spectrum.


title: Efficient $2^3S$ positronium production by stimulated decay from the $3^3P$ level

abstract: We investigate experimentally the possibility of enhancing the production of $2^3S$ positronium atoms by driving the $1^3S$-$3^3P$ and $3^3P$-$2^3S$ transitions, overcoming the natural branching ratio limitation of spontaneous decay from $3^3P$ to $2^3S$. The decay of $3^3P$ positronium atoms towards the $2^3S$ level has been effciently stimulated by a 1312.2nm broadband IR laser pulse. The dependence of the stimulating transition efficiency on the intensity of the IR pulse has been measured to find the optimal enhancement conditions. A maximum relative increase of $ \times (3.1 \pm 1.0) $ in the $2^3S$ production efficiency, with respect to the case where only spontaneous decay is present, was obtained.


title: Optimal approximation order of piecewise constants on convex partitions

abstract: We prove that the error of the best nonlinear $L_p$-approximation by piecewise constants on convex partitions is $\mathcal{O}\big(N^{-\frac{2}{d+1}}\big)$, where $N$ the number of cells, for all functions in the Sobolev space $W^2_q(\Omega)$ on a cube $\Omega\subset\mathbb{R}^d$, $d\geqslant 2$, as soon as $\frac{2}{d+1} + \frac{1}{p} - \frac{1}{q}\geqslant 0$. The approximation order $\mathcal{O}\big(N^{-\frac{2}{d+1}}\big)$ is achieved on a polyhedral partition obtained by anisotropic refinement of an adaptive dyadic partition. Further estimates of the approximation order from the above and below are given for various Sobolev and Sobolev-Slobodeckij spaces $W^r_q(\Omega)$ embedded in $L_p(\Omega)$, some of which also improve the standard estimate $\mathcal{O}\big(N^{-\frac 1d}\big)$ known to be optimal on isotropic partitions.


title: Performance of the Gemini Planet Imager Non-Redundant Mask and spectroscopy of two close-separation binaries HR 2690 and HD 142527

abstract: The Gemini Planet Imager (GPI) contains a 10-hole non-redundant mask (NRM), enabling interferometric resolution in complement to its coronagraphic capabilities. The NRM operates both in spectroscopic (integral field spectrograph, henceforth IFS) and polarimetric configurations. NRM observations were taken between 2013 and 2016 to characterize its performance. Most observations were taken in spectroscopic mode with the goal of obtaining precise astrometry and spectroscopy of faint companions to bright stars. We find a clear correlation between residual wavefront error measured by the AO system and the contrast sensitivity by comparing phase errors in observations of the same source, taken on different dates. We find a typical 5-$\sigma$ contrast sensitivity of $2-3~\times~10^{-3}$ at $\sim\lambda/D$. We explore the accuracy of spectral extraction of secondary components of binary systems by recovering the signal from a simulated source injected into several datasets. We outline data reduction procedures unique to GPI's IFS and describe a newly public data pipeline used for the presented analyses. We demonstrate recovery of astrometry and spectroscopy of two known companions to HR 2690 and HD 142527. NRM+polarimetry observations achieve differential visibility precision of $\sigma\sim0.4\%$ in the best case. We discuss its limitations on Gemini-S/GPI for resolving inner regions of protoplanetary disks and prospects for future upgrades. We summarize lessons learned in observing with NRM in spectroscopic and polarimetric modes.


title: DeepLocalization: Landmark-based Self-Localization with Deep Neural Networks

abstract: We address the problem of vehicle self-localization from multi-modal sensor information and a reference map. The map is generated off-line by extracting landmarks from the vehicle's field of view, while the measurements are collected similarly on the fly. Our goal is to determine the autonomous vehicle's pose from the landmark measurements and map landmarks. To learn this mapping, we propose DeepLocalization, a deep neural network that regresses the vehicle's translation and rotation parameters from unordered and dynamic input landmarks. The proposed network architecture is robust to changes of the dynamic environment and can cope with a small number of extracted landmarks. During the training process we rely on synthetically generated ground-truth. In our experiments, we evaluate two inference approaches in real-world scenarios. We show that DeepLocalization can be combined with regular GPS signals and filtering algorithms such as the extended Kalman filter. Our approach achieves state-of-the-art accuracy and is about ten times faster than the related work.


title: Interface engineering of graphene nanosheet reinforced ZrB$_2$ composites by tuning surface contacts

abstract: The mechanical properties of heterophase interfaces are critically important for the behaviour of graphene-reinforced composites. In this work, the structure, adhesion, cleavage and sliding of heterophase interfaces, formed between a ZrB$_2$ matrix and graphene nanosheets, are systematically investigated by density functional theory, and compared to available experimental data. We demonstrate that the surface chemistry of the ZrB$_2$ matrix material largely shapes the interface structures (of either Zr-C-Zr or B-C-B type) and the nature of the interfacial interaction. The Zr-C-Zr interfaces present strong chemical bonding and their response to mechanical stress is significantly influenced by graphene corrugation. In contrast B-C-B interfaces, interacting through the relatively weak $\pi$-$\pi$ stacking, show attributes similar to 2D materials heterostructures. Our theoretical results provide insights into the interface bonding mechanisms in graphene/ceramic composites, and emphasize the prospect for their design via interface engineering enabled by surface contacts.


title: Electrical control of magnetization in S/F/S junctions on a 3D topological insulator

abstract: Strong dependence of the Josephson energy on the magnetization orientation in Josephson junctions with ferromagnetic interlayers and spin-orbit coupling opens a way to control magnetization by Josephson current or Josephson phase. Here we investigate the perspectives of magnetization control in superconductor/ferromagnet/superconductor (S/F/S) Josephson junctions on the surface of a 3D topological insulator hosting Dirac quasiparticles. Due to the spin-momentum locking of these Dirac quasiparticles a strong dependence of the Josephson current-phase relation on the magnetization orientation is realized. It is demonstrated that this can lead to splitting of the ferromagnet's easy-axis in the voltage driven regime. We show that such a splitting can lead to stabilization of an unconventional four-fold degenerate ferromagnetic state.


title: Estimates and monotonicity for a heat flow of isometric G2-structures

abstract: Given a $7$-dimensional compact Riemannian manifold $\left( M,g\right) $ that admits $G_{2}$-structure, all the $G_{2}$-structures that are compatible with the metric $g$ are parametrized by unit sections of an octonion bundle over $M$. We define a natural energy functional on unit octonion sections and consider its associated heat flow. The critical points of this functional and flow precisely correspond to $G_{2}$-structures with divergence-free torsion. In this paper, we first derive estimates for derivatives of $V\left( t\right) $ along the flow and prove that the flow exists as long as the torsion remains bounded. We also prove a monotonicity formula and and an $\varepsilon $-regularity result for this flow. Finally, we show that within a metric class of $G_{2}$-structures that contains a torsion-free $G_{2}$-structure, under certain conditions, the flow will converge to a torsion-free $G_{2}$-structure.


title: Broadband lightweight flat lenses for longwave-infrared imaging

abstract: We experimentally demonstrate imaging in the longwave-infrared (LWIR) spectral band (8um to 12um) using a single polymer flat lens based upon multi-level diffractive optics. The device thickness is only 10{\mu}m, and chromatic aberrations are corrected over the entire LWIR band with one surface. Due to the drastic reduction in device thickness, we are able to utilize polymers with absorption in the LWIR, allowing for inexpensive manufacturing via imprint lithography. The weight of our lens is less than 100 times those of comparable refractive lenses. We fabricated and characterized two different flat lenses. Even with about 25% absorption losses, experiments show that our flat polymer lenses obtain good imaging with field of view of ~35degrees and angular resolution less than 0.013 degrees. The flat lenses were characterized with two different commercial LWIR image sensors. Finally, we show that by using lossless, higher-refractive-index materials like silicon, focusing efficiencies in excess of 70% can be achieved over the entire LWIR band. Our results firmly establish the potential for lightweight, ultra-thin, broadband lenses for high-quality imaging in the LWIR band.


title: Existence and stability analysis of solutions for a ultradian glucocorticoid rhythmicity and acute stress model

abstract: The hypothalamic pituitary adrenal (HPA) axis responds to physical and mental challenge to maintain homeostasis in part by controlling the body's cortisol level. Dysregulation of the HPA axis is implicated in numerous stress-related diseases. For a structured model of the HPA axis that includes the glucocorticoid receptor but does not take into account the system response delay, we analyze linear and non-linear stability of stationary solutions. For a second mathematical model that describes the mechanism of the HPA axis self-regulatory activities and takes into account a delay of system response, we prove existence of periodic solutions under certain assumptions on ranges of parameter values and analyze stability of these solutions with respect to the time delay value.


title: Self-Supervised Audio-Visual Co-Segmentation

abstract: Segmenting objects in images and separating sound sources in audio are challenging tasks, in part because traditional approaches require large amounts of labeled data. In this paper we develop a neural network model for visual object segmentation and sound source separation that learns from natural videos through self-supervision. The model is an extension of recently proposed work that maps image pixels to sounds. Here, we introduce a learning approach to disentangle concepts in the neural networks, and assign semantic categories to network feature channels to enable independent image segmentation and sound source separation after audio-visual training on videos. Our evaluations show that the disentangled model outperforms several baselines in semantic segmentation and sound source separation.


title: Semi-bandit Optimization in the Dispersed Setting

abstract: The goal of data-driven algorithm design is to obtain high-performing algorithms for specific application domains using machine learning and data. Across many fields in AI, science, and engineering, practitioners will often fix a family of parameterized algorithms and then optimize those parameters to obtain good performance on example instances from the application domain. In the online setting, we must choose algorithm parameters for each instance as they arrive, and our goal is to be competitive with the best fixed algorithm in hindsight. There are two major challenges in online data-driven algorithm design. First, it can be computationally expensive to evaluate the loss functions that map algorithm parameters to performance, which often require the learner to run a combinatorial algorithm to measure its performance. Second, the losses can be extremely volatile and have sharp discontinuities. However, we show that in many applications, evaluating the loss function for one algorithm choice can sometimes reveal the loss for a range of similar algorithms, essentially for free. We develop online optimization algorithms capable of using this kind of extra information by working in the semi-bandit feedback setting. Our algorithms achieve regret bounds that are essentially as good as algorithms under full-information feedback and are significantly more computationally efficient. We apply our semi-bandit results to obtain the first provable guarantees for data-driven algorithm design for linkage-based clustering and we improve the best regret bounds for designing greedy knapsack algorithms.


title: Decentralized and Parallel Primal and Dual Accelerated Methods for Stochastic Convex Programming Problems

abstract: We introduce primal and dual stochastic gradient oracle methods for decentralized convex optimization problems. Both for primal and dual oracles, the proposed methods are optimal in terms of the number of communication steps. However, for all classes of the objective, the optimality in terms of the number of oracle calls per node takes place only up to a logarithmic factor and the notion of smoothness. By using mini-batching technique, we show that the proposed methods with stochastic oracle can be additionally parallelized at each node. The considered algorithms can be applied to many data science problems and inverse problems.


title: An Inexact Interior-Point Lagrangian Decomposition Algorithm with Inexact Oracles

abstract: We develop a new inexact interior-point Lagrangian decomposition method to solve a wide range class of constrained composite convex optimization problems. Our method relies on four techniques: Lagrangian dual decomposition, self-concordant barrier smoothing, path-following, and proximal-Newton technique. It also allows one to approximately compute the solution of the primal subproblems (called the slave problems), which leads to inexact oracles (i.e., inexact gradients and Hessians) of the smoothed dual problem (called the master problem). The smoothed dual problem is nonsmooth, we propose to use an inexact proximal-Newton method to solve it. By appropriately controlling the inexact computation at both levels: the slave and master problems, we still estimate a polynomial-time iteration-complexity of our algorithm as in standard short-step interior-point methods. We also provide a strategy to recover primal solutions and establish complexity to achieve an approximate primal solution. We illustrate our method through two numerical examples on well-known models with both synthetic and real data and compare it with some existing state-of-the-art methods.


title: Intersection property and interaction decomposition

abstract: The decomposition into interaction subspaces is a hierarchical decomposition of the spaces of cylindrical functions of a finite product space, also called factor spaces. It is an important construction in graphical models and a standard way to prove the Hammersley-Clifford theorem that relates Markov fields to Gibbs fields and plays a central role in Kellerer's result for the linearized marginal problem. We define an intersection of sum property, or simply intersection property, and show that it characterizes collections of vector subspaces over a poset that can be hierarchically decomposed into direct sums, giving therefore a general setting for such construction to hold. We will call this generalization the interaction decomposition. The intersection property is the Bayesian intersection property when specified to factor spaces which, under this new perspective on the interaction decomposition, appears to be a structure property. An application is the extension of the decomposition into interaction subspaces for any product of any set.


title: A new approach in an analytical method for diffusion dynamics for the presence of delocalized sink in a potential well: Application to different potential curves

abstract: We provide a new approach to solve one dimension Fokker-Planck equation in the Laplace domain for the case where a particle is evolving in a potential energy curve in the presence of general delocalized sink. We also calculate rate constants in the presence of non-localized sink on different potential energy curves. In the previous method, we need to solve matrix equations to calculate rate constant but in our method, it is not required. We also calculate the rate constant by using the method to some known potential curve, viz, flat potential, linear potential, and parabolic potential.


title: Graph Element Networks: adaptive, structured computation and memory

abstract: We explore the use of graph neural networks (GNNs) to model spatial processes in which there is no a priori graphical structure. Similar to finite element analysis, we assign nodes of a GNN to spatial locations and use a computational process defined on the graph to model the relationship between an initial function defined over a space and a resulting function in the same space. We use GNNs as a computational substrate, and show that the locations of the nodes in space as well as their connectivity can be optimized to focus on the most complex parts of the space. Moreover, this representational strategy allows the learned input-output relationship to generalize over the size of the underlying space and run the same model at different levels of precision, trading computation for accuracy. We demonstrate this method on a traditional PDE problem, a physical prediction problem from robotics, and learning to predict scene images from novel viewpoints.


title: Genie: A Generator of Natural Language Semantic Parsers for Virtual Assistant Commands

abstract: To understand diverse natural language commands, virtual assistants today are trained with numerous labor-intensive, manually annotated sentences. This paper presents a methodology and the Genie toolkit that can handle new compound commands with significantly less manual effort. We advocate formalizing the capability of virtual assistants with a Virtual Assistant Programming Language (VAPL) and using a neural semantic parser to translate natural language into VAPL code. Genie needs only a small realistic set of input sentences for validating the neural model. Developers write templates to synthesize data; Genie uses crowdsourced paraphrases and data augmentation, along with the synthesized data, to train a semantic parser. We also propose design principles that make VAPL languages amenable to natural language translation. We apply these principles to revise ThingTalk, the language used by the Almond virtual assistant. We use Genie to build the first semantic parser that can support compound virtual assistants commands with unquoted free-form parameters. Genie achieves a 62% accuracy on realistic user inputs. We demonstrate Genie's generality by showing a 19% and 31% improvement over the previous state of the art on a music skill, aggregate functions, and access control.


title: A deep learning based solution for construction equipment detection: from development to deployment

abstract: This paper aims at providing researchers and engineering professionals with a practical and comprehensive deep learning based solution to detect construction equipment from the very first step of its development to the last one which is deployment. This paper focuses on the last step of deployment. The first phase of solution development, involved data preparation, model selection, model training, and model evaluation. The second phase of the study comprises of model optimization, application specific embedded system selection, and economic analysis. Several embedded systems were proposed and compared. The review of the results confirms superior real-time performance of the solutions with a consistent above 90% rate of accuracy. The current study validates the practicality of deep learning based object detection solutions for construction scenarios. Moreover, the detailed knowledge, presented in this study, can be employed for several purposes such as, safety monitoring, productivity assessments, and managerial decisions.


title: Triangular decomposition of character varieties

abstract: A marked surface is a compact oriented surface equipped with some pairwise disjoint arcs embedded in its boundary. In this paper, we extend the notion of character varieties to marked surfaces, in such a way that they have a nice behaviour for the operation of gluing two boundary arcs together. These stated character varieties are affine Poisson varieties which coincide with the Culler-Shalen character varieties when the surface is unmarked and are closely related to the Fock-Rosly and Alekseev-Kosmann-Malkin-Meinrenken constructions in the marked case. These Poisson varieties are the classical moduli spaces underlying stated skein algebras and share similar properties. In particular, stated character varieties admit triangular decompositions, associated to triangulations of the surface. We identify the Zariski tangent spaces of these varieties with some twisted groupoid cohomological groups and provide a generalization of Goldman's formula for the Poisson bracket of curve functions in terms of intersection form in homology.


title: Making Meaning: Semiotics Within Predictive Knowledge Architectures

abstract: Within Reinforcement Learning, there is a fledgling approach to conceptualizing the environment in terms of predictions. Central to this predictive approach is the assertion that it is possible to construct ontologies in terms of predictions about sensation, behaviour, and time---to categorize the world into entities which express all aspects of the world using only predictions. This construction of ontologies is integral to predictive approaches to machine knowledge where objects are described exclusively in terms of how they are perceived. In this paper, we ground the Pericean model of semiotics in terms of Reinforcement Learning Methods, describing Peirce's Three Categories in the notation of General Value Functions. Using the Peircean model of semiotics, we demonstrate that predictions alone are insufficient to construct an ontology; however, we identify predictions as being integral to the meaning-making process. Moreover, we discuss how predictive knowledge provides a particularly stable foundation for semiosis\textemdash the process of making meaning\textemdash and suggest a possible avenue of research to design algorithmic methods which construct semantics and meaning using predictions.


title: When is a Prediction Knowledge?

abstract: Within Reinforcement Learning, there is a growing collection of research which aims to express all of an agent's knowledge of the world through predictions about sensation, behaviour, and time. This work can be seen not only as a collection of architectural proposals, but also as the beginnings of a theory of machine knowledge in reinforcement learning. Recent work has expanded what can be expressed using predictions, and developed applications which use predictions to inform decision-making on a variety of synthetic and real-world problems. While promising, we here suggest that the notion of predictions as knowledge in reinforcement learning is as yet underdeveloped: some work explicitly refers to predictions as knowledge, what the requirements are for considering a prediction to be knowledge have yet to be well explored. This specification of the necessary and sufficient conditions of knowledge is important; even if claims about the nature of knowledge are left implicit in technical proposals, the underlying assumptions of such claims have consequences for the systems we design. These consequences manifest in both the way we choose to structure predictive knowledge architectures, and how we evaluate them. In this paper, we take a first step to formalizing predictive knowledge by discussing the relationship of predictive knowledge learning methods to existing theories of knowledge in epistemology. Specifically, we explore the relationships between Generalized Value Functions and epistemic notions of Justification and Truth.


title: Radical pair intersystem crossing: Quantum dynamics or incoherent kinetics?

abstract: Magnetic field effects on radical pair reactions arise due to the interplay of coherent electron spin dynamics and spin relaxation effects, a rigorous treatment of which requires the solution of the Liouville-von Neumann equation. However, it is often found that simple incoherent kinetic models of the radical pair singlet-triplet intersystem crossing provide an acceptable description of experimental measurements. In this paper we outline the theoretical basis for this incoherent kinetic description, elucidating its connection to exact quantum mechanics. We show in particular how the finite lifetime of the radical pair spin states, as well as any additional spin-state dephasing, leads to incoherent intersystem crossing. We arrive at simple expressions for the radical pair spin state interconversion rates to which the functional form proposed recently by Steiner et al. [J. Phys. Chem. C 122, 11701 (2018)] can be regarded as an approximation. We also test the kinetic master equation against exact quantum dynamical simulations for a model radical pair and for a series of $\text{PTZ}^{\bullet+}\text{-Ph}_\text{n}\text{-PDI}^{\bullet-}$ molecular wires.


title: Co-isometric weighted composition operators on Hilbert spaces of analytic functions

abstract: We obtain a necessary and sufficient condition for a weighted composition operator to be co-isometric on a general weighted Hardy space of analytic functions in the unit disk whose reproducing kernel has the usual natural form. This turns out to be equivalent to the property of being unitary. The result reveals a dichotomy identifying a specific family of weighted Hardy spaces as the only ones that support non-trivial operators of this kind.


title: Adaptive Huber Regression on Markov-dependent Data

abstract: High-dimensional linear regression has been intensively studied in the community of statistics in the last two decades. For the convenience of theoretical analyses, classical methods usually assume independent observations and sub-Gaussian-tailed errors. However, neither of them hold in many real high-dimensional time-series data. Recently [Sun, Zhou, Fan, 2019, J. Amer. Stat. Assoc., in press] proposed Adaptive Huber Regression (AHR) to address the issue of heavy-tailed errors. They discover that the robustification parameter of the Huber loss should adapt to the sample size, the dimensionality, and the moments of the heavy-tailed errors. We progress in a vertical direction and justify AHR on dependent observations. Specifically, we consider an important dependence structure -- Markov dependence. Our results show that the Markov dependence impacts on the adaption of the robustification parameter and the estimation of regression coefficients in the way that the sample size should be discounted by a factor depending on the spectral gap of the underlying Markov chain.


title: A Novel BiLevel Paradigm for Image-to-Image Translation

abstract: Image-to-image (I2I) translation is a pixel-level mapping that requires a large number of paired training data and often suffers from the problems of high diversity and strong category bias in image scenes. In order to tackle these problems, we propose a novel BiLevel (BiL) learning paradigm that alternates the learning of two models, respectively at an instance-specific (IS) and a general-purpose (GP) level. In each scene, the IS model learns to maintain the specific scene attributes. It is initialized by the GP model that learns from all the scenes to obtain the generalizable translation knowledge. This GP initialization gives the IS model an efficient starting point, thus enabling its fast adaptation to the new scene with scarce training data. We conduct extensive I2I translation experiments on human face and street view datasets. Quantitative results validate that our approach can significantly boost the performance of classical I2I translation models, such as PG2 and Pix2Pix. Our visualization results show both higher image quality and more appropriate instance-specific details, e.g., the translated image of a person looks more like that person in terms of identity.


title: Deep Learning for Power System Security Assessment

abstract: Security assessment is among the most fundamental functions of power system operator. The sheer complexity of power systems exceeding a few buses, however, makes it an extremely computationally demanding task. The emergence of deep learning methods that are able to handle immense amounts of data, and infer valuable information appears as a promising alternative. This paper has two main contributions. First, inspired by the remarkable performance of convolutional neural networks for image processing, we represent for the first time power system snapshots as 2-dimensional images, thus taking advantage of the wide range of deep learning methods available for image processing. Second, we train deep neural networks on a large database for the NESTA 162-bus system to assess both N-1 security and small-signal stability. We find that our approach is over 255 times faster than a standard small-signal stability assessment, and it can correctly determine unsafe points with over 99% accuracy.


title: Automatic Target Detection for Sparse Hyperspectral Images

abstract: In this work, a novel target detector for hyperspectral imagery is developed. The detector is independent on the unknown covariance matrix, behaves well in large dimensions, distributional free, invariant to atmospheric effects, and does not require a background dictionary to be constructed. Based on a modification of the robust principal component analysis (RPCA), a given hyperspectral image (HSI) is regarded as being made up of the sum of a low-rank background HSI and a sparse target HSI that contains the targets based on a pre-learned target dictionary specified by the user. The sparse component is directly used for the detection, that is, the targets are simply detected at the non-zero entries of the sparse target HSI. Hence, a novel target detector is developed, which is simply a sparse HSI generated automatically from the original HSI, but containing only the targets with the background is suppressed. The detector is evaluated on real experiments, and the results of which demonstrate its effectiveness for hyperspectral target detection especially when the targets are well matched to the surroundings.


title: Predict Future Sales using Ensembled Random Forests

abstract: This is a method report for the Kaggle data competition 'Predict future sales'. In this paper, we propose a rather simple approach to future sales predicting based on feature engineering, Random Forest Regressor and ensemble learning. Its performance turned out to exceed many of the conventional methods and get final score 0.88186, representing root mean squared error. As of this writing, our model ranked 5th on the leaderboard. (till 8.5.2018)


title: #Cyberbullying in the Digital Age: Exploring People's Opinions with Text Mining

abstract: This study used text mining to investigate people's insights about cyberbullying. English-language tweets were collected and analyzed by R software. Our analysis demonstrated three major themes: (a) the major actions that needed to be taken into consideration (e.g. guiding parents and teachers to cyberbullying prevention, funding schools to fight cyberbullying), (b) certain events that were important to people (e.g. the Michigan cyberbullying law), and (c) people's major concerns in this regard (e.g. mental health issues among students). Parents and teachers have an important role in educating, informing, warning, preventing, and protecting against cyberbullying behaviors. The frequency of negative sentiments was almost 2.45 times more than positive sentiments.


title: Towards Solving the Navier-Stokes Equation on Quantum Computers

abstract: In this paper, we explore the suitability of upcoming novel computing technologies, in particular adiabatic annealing based quantum computers, to solve fluid dynamics problems that form a critical component of several science and engineering applications. We start with simple flows with well-studied flow properties, and provide a framework to convert such systems to a form amenable for deployment on such quantum annealers. We analyze the solutions obtained both qualitatively and quantitatively as well as the sensitivities of the various solution selection schemes on the obtained solution.

categories: cs.NA math.NA physics.comp-ph physics.flu-dyn


title: A "rare'' plane set with Hausdorff dimension 2

abstract: We prove that for every at most countable family $\{f_k(x)\}$ of real functions on $[0,1)$ there is a single-valued real function $F(x)$, $x\in[0,1)$, such that the Hausdorff dimension of the graph $\Gamma$ of $F(x)$ equals 2, and for every $C\in\mathbb R$ and every $k$, the intersection of $\Gamma$ with the graph of the function $f_k(x)+C$ consists of at most one point. We also construct a family of functions of cardinality continuum and a function $F$ with similar properties.