title: An Empirical Comparison of Developer Retention in the RubyGems and npm Software Ecosystems

abstract: Software ecosystems can be viewed as socio-technical networks consisting of technical components (software packages) and social components (communities of developers) that maintain the technical components. Ecosystems evolve over time through socio-technical changes that may greatly impact the ecosystem's sustainability. Social changes like developer turnover may lead to technical degradation. This motivates the need to identify those factors leading to developer abandonment, in order to automate the process of identifying developers with high abandonment risk. This paper compares such factors for two software package ecosystems, RubyGems and npm. We analyse the evolution of their packages hosted on GitHub, considering development activity in terms of commits, and social interaction with other developers in terms of comments associated to commits, issues or pull requests. We analyse this socio-technical activity for more than 30k and 60k developers for RubyGems and npm respectively. We use survival analysis to identify which factors coincide with a lower survival probability. Our results reveal that developers with a higher probability to abandon an ecosystem: do not engage in discussions with other developers; do not have strong social and technical activity intensity; communicate or commit less frequently; and do not participate to both technical and social activities for long periods of time. Such observations could be used to automate the identification of developers with a high probability of abandoning the ecosystem and, as such, reduce the risks associated to knowledge loss.


title: Boolean perspectives of idioms and the Boyle derivative

abstract: We are concerned with the boolean or more general with the complemented properties of idioms (complete upper-continuous modular lattices). In [Simmons&Cantor] the author introduces a device which captures in some informal speaking how far the idiom is from be complemented, this device is the Cantor-Bendixson derivative. There exists another device that captures some boolean properties, the so-called Boyle-derivative, this derivative is an operator on the assembly (the frame of nuclei) of the idiom. The Boyle-derivative has its origins in module theory. In this investigation we produce an idiomatic analysis of the boolean properties of any idiom using the Boyle-derivative, we give conditions on a nucleus $j$ such that $[j, tp]$ is a complete boolean algebra. We also explore some properties of nuclei $j$ such that $A_{j}$ is a complemented idiom.


title: Multilayer Spectral Graph Clustering via Convex Layer Aggregation: Theory and Algorithms

abstract: Multilayer graphs are commonly used for representing different relations between entities and handling heterogeneous data processing tasks. Non-standard multilayer graph clustering methods are needed for assigning clusters to a common multilayer node set and for combining information from each layer. This paper presents a multilayer spectral graph clustering (SGC) framework that performs convex layer aggregation. Under a multilayer signal plus noise model, we provide a phase transition analysis of clustering reliability. Moreover, we use the phase transition criterion to propose a multilayer iterative model order selection algorithm (MIMOSA) for multilayer SGC, which features automated cluster assignment and layer weight adaptation, and provides statistical clustering reliability guarantees. Numerical simulations on synthetic multilayer graphs verify the phase transition analysis, and experiments on real-world multilayer graphs show that MIMOSA is competitive or better than other clustering methods.


title: Real Time Analytics: Algorithms and Systems

abstract: Velocity is one of the 4 Vs commonly used to characterize Big Data. In this regard, Forrester remarked the following in Q3 2014: "The high velocity, white-water flow of data from innumerable real-time data sources such as market data, Internet of Things, mobile, sensors, click-stream, and even transactions remain largely unnavigated by most firms. The opportunity to leverage streaming analytics has never been greater." Example use cases of streaming analytics include, but not limited to: (a) visualization of business metrics in real-time (b) facilitating highly personalized experiences (c) expediting response during emergencies. Streaming analytics is extensively used in a wide variety of domains such as healthcare, e-commerce, financial services, telecommunications, energy and utilities, manufacturing, government and transportation. In this tutorial, we shall present an in-depth overview of streaming analytics - applications, algorithms and platforms - landscape. We shall walk through how the field has evolved over the last decade and then discuss the current challenges - the impact of the other three Vs, viz., Volume, Variety and Veracity, on Big Data streaming analytics. The tutorial is intended for both researchers and practitioners in the industry. We shall also present state-of-the-affairs of streaming analytics at Twitter.


title: Kinematic interpretation of the Study quadric's ambient space

abstract: It is well known that real points of the Study quadric (sliced along a 3-dimensional generator space) correspond to displacements of the Euclidean 3-space. But we still lack of a kinematic meaning for the points of the ambient 7-dimensional projective space $P^7$. This paper gives one possible interpretation in terms of displacements of the Euclidean 4-space. From this point of view we also discuss the extended inverse kinematic map, motions corresponding to straight lines in $P^7$ and linear complexes of SE(3)-displacements. Moreover we present an application of this interpretation in the context of interactive motion design.


title: Quantum simulation of rainbow gravity by nonlocal nonlinearity

abstract: Rainbow gravity modifies general relativity by introducing an energy dependent metric, which is expected to have a role in the quantum theory of black holes and in quantum gravity at Planck energy scale. We show that rainbow gravity can be simulated in the laboratory by nonlinear waves in nonlocal media, as those occurring in Bose-condensed gases and nonlinear optics. We reveal that at a classical level, a nonlocal nonlinear Schr\"odinger equation may emulate the curved space time in proximity of a rotating black hole as dictated by the rainbow gravity scenario. We also demonstrate that a fully quantized analysis is possible. By the positive $\mathcal{P}$-representation, we study superradiance and show that the instability of a black-hole and the existence of an event horizon are inhibited by an energy dependent metric. Our results open the way to a number of fascinating experimental tests of quantum gravity theories and quantum field theory in curved manifolds, and also demonstrate that these theories may be novel tools for open problems in nonlinear quantum physics.


title: On the equivalence between different canonical forms of F(R) theory of gravity

abstract: Classical equivalence between Jordan's and Einstein's frame counterparts of F(R) theory of gravity has recently been questioned, since the two produce different Noether symmetries, which couldn't be translated back and forth using transformation relations. Here we add the Hamiltonian constraint equation, which is essentially the time-time component of Einstein's equation, through a Lagrange multiplier to the existence condition for Noether symmetry to show that all the three different canonical structures of F(R) theory of gravity, including the one which follows from Lagrange multiplier technique, admit each and every available symmetry independently. This establishes classical equivalence amongst all the three.


title: Risk Constrained Trading Strategies for Stochastic Generation with a Single-Price Balancing Market

abstract: Due to the limited predictability of wind power and other stochastic generation, trading this energy in competitive electricity markets is challenging. This paper derives revenue-maximising and risk-constrained strategies for stochastic generators participating in electricity markets with a single-price balancing mechanism. Starting from the optimal---and impractical---strategy of offering zero or nominal power, which exposes the participant to potentially large imbalance costs, we develop a number of strategies that control risk by hedging against penalising balancing prices in favour of rewarding ones. Trading strategies are formulated in a probabilistic framework in order to address asymmetry in balancing prices. The large-scale communication of system information characteristic of modern power systems is utilised to inputs for electricity price forecasts and probabilistic system length forecasts. A case study using data from the GB market in the UK is presented and the ability of the proposed strategies to increase revenue and reduce risk is demonstrated and analysed.


title: A combinatorial method for connecting BHV spaces representing different numbers of taxa

abstract: The phylogenetic tree space introduced by Billera, Holmes, and Vogtmann (BHV tree space) is a CAT(0) continuous space that represents trees with edge weights with an intrinsic geodesic distance measure. The geodesic distance measure unique to BHV tree space is well known to be computable in polynomial time, which makes it a potentially powerful tool for optimization problems in phylogenetics and phylogenomics. Specifically, there is significant interest in comparing and combining phylogenetic trees. For example, BHV tree space has been shown to be potentially useful in tree summary and consensus methods, which require combining trees with different number of leaves. Yet an open problem is to transition between BHV tree spaces of different maximal dimension, where each maximal dimension corresponds to the complete set of edge-weighted trees with a fixed number of leaves. We show a combinatorial method to transition between copies of BHV tree spaces in which trees with different numbers of taxa can be studied, derived from its topological structure and geometric properties. This method removes obstacles for embedding problems such as supertree and consensus methods in the BHV treespace framework.


title: Casimir Effect in the Rainbow Einstein's Universe

abstract: In the present paper we investigate the effects caused by the modification of the dispersion relation obtained by solving the Klein-Gordon equation in the closed Einstein's universe in the context of rainbow's gravity models. Thus, we analyse how the quantum vacuum fluctuations of the scalar field are modified when compared with the results obtained in the usual General Relativity scenario. The regularization, and consequently the renormalization, of the vacuum energy is performed adopting the Epstein-Hurwitz and Riemann's zeta functions.


title: What is the best planar cavity for maximizing coherent exciton-photon coupling

abstract: We compare alternative planar cavity structures for strong exciton$-$photon coupling, where the conventional distributed Bragg reflector (DBR) and three unconventional types of cavity mirrors$-$ air/GaAs DBR, Tamm $-$ plasmon mirror and sub$-$wavelength grating mirror. We design and optimize the planar cavities built with each type of mirror at one side or both sides for maximum vacuum field strength. We discuss the trade$-$off between performance and fabrication difficulty for each cavity structure. We show that cavities with sub$-$wavelength grating mirrors allow simultaneously strongest field and high cavity quality. The optimization principles and techniques developed in this work will guide the cavity design for research and applications of matter$-$light coupled semiconductors, especially new material systems that require greater flexibility in the choice of cavity materials and cavity fabrication procedures.


title: Protecting Genomic Privacy by a Sequence-Similarity Based Obfuscation Method

abstract: In the post-genomic era, large-scale personal DNA sequences are produced and collected for genetic medical diagnoses and new drug discovery, which, however, simultaneously poses serious challenges to the protection of personal genomic privacy. Existing genomic privacy-protection methods are either time-consuming or with low accuracy. To tackle these problems, this paper proposes a sequence similarity-based obfuscation method, namely IterMegaBLAST, for fast and reliable protection of personal genomic privacy. Specifically, given a randomly selected sequence from a dataset of DNA sequences, we first use MegaBLAST to find its most similar sequence from the dataset. These two aligned sequences form a cluster, for which an obfuscated sequence was generated via a DNA generalization lattice scheme. These procedures are iteratively performed until all of the sequences in the dataset are clustered and their obfuscated sequences are generated. Experimental results on two benchmark datasets demonstrate that under the same degree of anonymity, IterMegaBLAST significantly outperforms existing state-of-the-art approaches in terms of both utility accuracy and time complexity.


title: A GRB and Broad-lined Type Ic Supernova from a Single Central Engine

abstract: Unusually high velocities (< ~0.1c) and correspondingly high kinetic energies have been observed in a subset of Type Ic supernovae (so-called "broad-lined Ic" supernovae; SNe Ic-BL), prompting a search for a central engine model capable of generating such energetic explosions. A clue to the explosion mechanism may lie in the fact that all supernovae that accompany long-duration gamma-ray bursts belong to the SN Ic-BL class. Using a combination of two-dimensional relativistic hydrodynamics and radiation transport calculations, we demonstrate that the central engine responsible for long gamma-ray bursts can also trigger a SN Ic-BL. We find that a reasonable gamma-ray burst engine injected into a stripped Wolf-Rayet progenitor produces a relativistic jet with energy ~10^51 ergs, as well as a SN whose synthetic light curves and spectra are fully consistent with observed SNe Ic-BL during the photospheric phase. As a result of the jet's asymmetric energy injection, the SN spectra and light curves depend on viewing angle. The impact of viewing angle on the spectrum is particularly pronounced at early times, while the viewing angle dependence for the light curves (~10% variation in bolometric luminosity) persists throughout the photospheric phase.


title: Early formation of planetary building blocks inferred from Pb isotopic ages of chondrules

abstract: The most abundant components of primitive meteorites (chondrites) are millimeter-sized glassy spherical chondrules formed by transient melting events in the solar protoplanetary disk. Using Pb-Pb dates of 22 individual chondrules, we show that primary production of chondrules in the early solar system was restricted to the first million years after formation of the Sun and that these existing chondrules were recycled for the remaining lifetime of the protoplanetary disk. This is consistent with a primary chondrule formation episode during the early high-mass accretion phase of the protoplanetary disk that transitions into a longer period of chondrule reworking. An abundance of chondrules at early times provides the precursor material required to drive the efficient and rapid formation of planetary objects via chondrule accretion.


title: Using JAGS for Bayesian Cognitive Diagnosis Modeling: A Tutorial

abstract: In this article, the JAGS software program is systematically introduced to fit common Bayesian cognitive diagnosis models (CDMs), including the deterministic inputs, noisy "and" gate (DINA) model, the deterministic inputs, noisy "or" gate (DINO) model, the linear logistic model, the reduced reparameterized unified model (rRUM), and the log-linear CDM (LCDM). The unstructured latent structural model and the higher-order latent structural model are both introduced. We also show how to extend those models to consider the polytomous attributes, the testlet effect, and the longitudinal diagnosis. Finally, an empirical example is presented as a tutorial to illustrate how to use the JAGS codes in R.


title: Radiation mediated shocks in gamma-ray bursts: Pair creation

abstract: Sub-photospheric shock dissipation is one of the main proposed mechanisms for producing the prompt gamma-ray burst (GRB) emission. Such shocks are mediated by scattering of radiation. We introduce a time dependent, special relativistic code which dynamically couples Monte Carlo radiative transfer to the flow hydrodynamics. The code also self-consistently implements electron-positron pair production and annihilation. We simulate shocks with properties relevant for GRBs and study the steady-state solutions, which are accurate deep below the jet photosphere. The shock generates a power-law photon spectrum through the first-order Fermi mechanism, extending upwards from the typical upstream photon energy. Strong shocks (for which the downstream pressure is much larger than the upstream pressure) have rising $\nu F_\nu$ shock spectra. The spectrum extends up to $\epsilon_{max} \equiv E_{max}/m_e c^2 \sim v^2$ for non-relativistic shocks, where $m_e$ is the electron rest mass and $v$ is the relative speed between the upstream and downstream in units of the speed of light $c$. For mildly relativistic shocks the power law softens at $\epsilon \gtrsim 10^{-1}$ due to Klein-Nishina effects, and shocks with $v\gamma \gtrsim 1$, where $\gamma \equiv (1-v^2)^{-1/2}$, produce electron-positron pairs. As an example, a strong shock with $v\gamma = 3$ and a photon-to-proton ratio of $n_\gamma/n_p = 2 \times 10^5$ has a peak pair-to-proton ratio of $Z_\pm \approx 225$. The main effect of pairs in a steady-state shock is to decrease its spatial width by a factor of $\sim Z_\pm$. The post-shock spectrum thermalizes in the downstream. In absence of emission and absorption processes, kinetic equilibrium at temperature $\theta_d \equiv kT_d/m_e c^2 \approx \epsilon_d/3$ is reached at an optical depth of $\tau \gg \theta_d^{-1}$ behind the shock, where $\epsilon_d$ is the average downstream photon energy.


title: Generation of high-fidelity quantum control methods for multi-level systems

abstract: In recent decades there has been a rapid development of methods to experimentally control individual quantum systems. A broad range of quantum control methods has been developed for two-level systems, however the complexity of multi-level quantum systems make the development of analogous control methods extremely challenging. Here, we exploit the equivalence between multi-level systems with SU(2) symmetry and spin-1/2 systems to develop a technique for generating new robust, high-fidelity, multi-level control methods. As a demonstration of this technique, we develop new adiabatic and composite multi-level quantum control methods and experimentally realise these methods using an $^{171}$Yb$^+$ ion system. We measure the average infidelity of the process in both cases to be around $10^{-4}$, demonstrating that this technique can be used to develop high-fidelity multi-level quantum control methods and can, for example, be applied to a wide range of quantum computing protocols including implementations below the fault-tolerant threshold in trapped ions.


title: Anomaly Detection in Multivariate Non-stationary Time Series for Automatic DBMS Diagnosis

abstract: Anomaly detection in database management systems (DBMSs) is difficult because of increasing number of statistics (stat) and event metrics in big data system. In this paper, I propose an automatic DBMS diagnosis system that detects anomaly periods with abnormal DB stat metrics and finds causal events in the periods. Reconstruction error from deep autoencoder and statistical process control approach are applied to detect time period with anomalies. Related events are found using time series similarity measures between events and abnormal stat metrics. After training deep autoencoder with DBMS metric data, efficacy of anomaly detection is investigated from other DBMSs containing anomalies. Experiment results show effectiveness of proposed model, especially, batch temporal normalization layer. Proposed model is used for publishing automatic DBMS diagnosis reports in order to determine DBMS configuration and SQL tuning.


title: Perron-Frobenius theory for kernels and Crump-Mode-Jagers processes with macro-individuals

abstract: Perron-Frobenius theory developed for irreducible non-negative kernels deals with so-called $R$-positive recurrent kernels. If kernel $M$ is $R$-positive recurrent, then the main result determines the limit of the scaled kernel iterations $R^nM^n$ as $n\to\infty$. In the Nummelin's monograph this important result is proven using a regeneration method whose major focus is on $M$ having an atom. In the special case when $M=P$ is a stochastic kernel with an atom, the regeneration method has an elegant explaination in terms of an associated split chain. In this paper we give a new probabilistic interpretation of the general regeneration method in terms of multi-type Galton-Watson processes producing clusters of particles. Treating clusters as macro-individuals, we arrive at a single-type Crump-Mode-Jagers process with a naturally embedded renewal structure.


title: TensorFlow Estimators: Managing Simplicity vs. Flexibility in High-Level Machine Learning Frameworks

abstract: We present a framework for specifying, training, evaluating, and deploying machine learning models. Our focus is on simplifying cutting edge machine learning for practitioners in order to bring such technologies into production. Recognizing the fast evolution of the field of deep learning, we make no attempt to capture the design space of all possible model architectures in a domain- specific language (DSL) or similar configuration language. We allow users to write code to define their models, but provide abstractions that guide develop- ers to write models in ways conducive to productionization. We also provide a unifying Estimator interface, making it possible to write downstream infrastructure (e.g. distributed training, hyperparameter tuning) independent of the model implementation. We balance the competing demands for flexibility and simplicity by offering APIs at different levels of abstraction, making common model architectures available out of the box, while providing a library of utilities designed to speed up experimentation with model architectures. To make out of the box models flexible and usable across a wide range of problems, these canned Estimators are parameterized not only over traditional hyperparameters, but also using feature columns, a declarative specification describing how to interpret input data. We discuss our experience in using this framework in re- search and production environments, and show the impact on code health, maintainability, and development speed.


title: Distributed rank-1 dictionary learning: Towards fast and scalable solutions for fMRI big data analytics

abstract: The use of functional brain imaging for research and diagnosis has benefitted greatly from the recent advancements in neuroimaging technologies, as well as the explosive growth in size and availability of fMRI data. While it has been shown in literature that using multiple and large scale fMRI datasets can improve reproducibility and lead to new discoveries, the computational and informatics systems supporting the analysis and visualization of such fMRI big data are extremely limited and largely under-discussed. We propose to address these shortcomings in this work, based on previous success in using dictionary learning method for functional network decomposition studies on fMRI data. We presented a distributed dictionary learning framework based on rank-1 matrix decomposition with sparseness constraint (D-r1DL framework). The framework was implemented using the Spark distributed computing engine and deployed on three different processing units: an in-house server, in-house high performance clusters, and the Amazon Elastic Compute Cloud (EC2) service. The whole analysis pipeline was integrated with our neuroinformatics system for data management, user input/output, and real-time visualization. Performance and accuracy of D-r1DL on both individual and group-wise fMRI Human Connectome Project (HCP) dataset shows that the proposed framework is highly scalable. The resulting group-wise functional network decompositions are highly accurate, and the fast processing time confirm this claim. In addition, D-r1DL can provide real-time user feedback and results visualization which are vital for large-scale data analysis.


title: Extractor-Based Time-Space Lower Bounds for Learning

abstract: A matrix $M: A \times X \rightarrow \{-1,1\}$ corresponds to the following learning problem: An unknown element $x \in X$ is chosen uniformly at random. A learner tries to learn $x$ from a stream of samples, $(a_1, b_1), (a_2, b_2) \ldots$, where for every $i$, $a_i \in A$ is chosen uniformly at random and $b_i = M(a_i,x)$. Assume that $k,\ell, r$ are such that any submatrix of $M$ of at least $2^{-k} \cdot |A|$ rows and at least $2^{-\ell} \cdot |X|$ columns, has a bias of at most $2^{-r}$. We show that any learning algorithm for the learning problem corresponding to $M$ requires either a memory of size at least $\Omega\left(k \cdot \ell \right)$, or at least $2^{\Omega(r)}$ samples. The result holds even if the learner has an exponentially small success probability (of $2^{-\Omega(r)}$). In particular, this shows that for a large class of learning problems, any learning algorithm requires either a memory of size at least $\Omega\left((\log |X|) \cdot (\log |A|)\right)$ or an exponential number of samples, achieving a tight $\Omega\left((\log |X|) \cdot (\log |A|)\right)$ lower bound on the size of the memory, rather than a bound of $\Omega\left(\min\left\{(\log |X|)^2,(\log |A|)^2\right\}\right)$ obtained in previous works [R17,MM17b]. Moreover, our result implies all previous memory-samples lower bounds, as well as a number of new applications. Our proof builds on [R17] that gave a general technique for proving memory-samples lower bounds.


title: Time-Space Tradeoffs for Learning from Small Test Spaces: Learning Low Degree Polynomial Functions

abstract: We develop an extension of recently developed methods for obtaining time-space tradeoff lower bounds for problems of learning from random test samples to handle the situation where the space of tests is signficantly smaller than the space of inputs, a class of learning problems that is not handled by prior work. This extension is based on a measure of how matrices amplify the 2-norms of probability distributions that is more refined than the 2-norms of these matrices. As applications that follow from our new technique, we show that any algorithm that learns $m$-variate homogeneous polynomial functions of degree at most $d$ over $\mathbb{F}_2$ from evaluations on randomly chosen inputs either requires space $\Omega(mn)$ or $2^{\Omega(m)}$ time where $n=m^{\Theta(d)}$ is the dimension of the space of such functions. These bounds are asymptotically optimal since they match the tradeoffs achieved by natural learning algorithms for the problems.


title: Comodule Algebras and 2-Cocycles over the (Braided) Drinfeld Double

abstract: We show that for dually paired bialgebras, every comodule algebra over one of the paired bialgebras gives a comodule algebra over their Drinfeld double via a crossed product construction. These constructions generalize to working with bialgebra objects in a braided monoidal category of modules over a quasitriangular Hopf algebra. Hence two ways to provide comodule algebras over the braided Drinfeld double (the double bosonization) are provided. Furthermore, a map of second Hopf algebra cohomology spaces is constructed. It takes a pair of 2-cocycles over dually paired Hopf algebras and produces a 2-cocycle over their Drinfeld double. This construction also has an analogue for braided Drinfeld doubles.


title: Polynomial functors in manifold calculus

abstract: Let M be a smooth manifold, and let O(M) be the poset of open subsets of M. Manifold calculus, due to Goodwillie and Weiss, is a calculus of functors suitable for studying contravariant functors (cofunctors) F: O(M)--> Top from O(M) to the category of spaces. Weiss showed that polynomial cofunctors of degree <= k are determined by their values on O_k(M), where O_k(M) is the full subposet of O(M) whose objects are open subsets diffeomorphic to the disjoint union of at most k balls. Afterwards Pryor showed that one can replace O_k(M) by more general subposets and still recover the same notion of polynomial cofunctor. In this paper, we generalize these results to cofunctors from O(M) to any simplicial model category C. If conf(k, M) stands for the unordered configuration space of k points in M, we also show that the category of homogeneous cofunctors O(M) --> C of degree k is weakly equivalent to the category of linear cofunctors O(conf(k, M)) --> C provided that C has a zero object. Using a completely different approach, we also show that if C is a general model category and F: O_k(M) --> C is an isotopy cofunctor, then the homotopy right Kan extension of F along the inclusion O_k(M) --> O(M) is also an isotopy cofunctor.


title: Enabling Narrow(est) IWA Coronagraphy with STIS BAR5 and BAR10 Occulters

abstract: The Space Telescope Imaging Spectrograph's (STIS) BAR5 coronagraphic occulter was designed to provide high-contrast, visible-light, imaging in close (>= 0.15") angular proximity to bright point-sources. This is the smallest inner working angle (IWA) possible with HST's suite of coronagraphically augmented instruments through its mission lifetime. The STIS BAR5 image plane occulter, however, was damaged (bent and deformed) pre-launch and had not been enabled for GO science use following the installation of the instrument in 1997, during HST servicing mission SM2. With the success of the HST GO 12923 program, discussed herein, we explored and verified the functionality and utility of the BAR5 occulter. Thus, despite its physical damage, with updates to the knowledge of the aperture mask metrology and target pointing requirements, a robust determination of achievable raw and PSF-subtracted stellocentric image contrasts and fidelity was conducted. We also investigated, and herein report on, the use of the BAR10 rounded corners as narrow-angle occulters and compare IWA vs. contrast performance for the BAR5, BAR10, and Wedge occulters. With that, we provide recommendations for the most efficacious BAR5 and BAR10 use on-orbit in support of GO science. With color-matched PSF-template subtracted coronagraphy, inclusive of a small (+/- 1/4 pixel) 3-point cross-bar dithering strategy we recommend, we find BAR5 can deliver effective ~ 0.2" IWA image contrast of ~ 4 x 10^-5 pixel^-1 to ~ 1 x 10^-8 pixel^-1 at 2". With the pointing updates (to the PDB SIAF.dat file and/or implemented through APT) that we identified, and with observing strategies we explored, we recommend the use of STIS BAR5 coronagraphy as a fully "supported" capability for unique GO science.


title: The Impact of Baryons on the Large-Scale Structure of the Universe

abstract: Numerical simulations play an important role in current astronomy researches. Previous dark-matter-only simulations have represented the large-scale structure of the Universe. However, nowadays, hydro-dynamical simulations with baryonic models, which can directly present realistic galaxies, may twist these results from dark-matter-only simulations. In this chapter, we mainly focus on these three statistical methods: power spectrum, two-point correlation function and halo mass function, which are normally used to characterize the large-scale structure of the Universe. We review how these baryon processes influence the cosmology structures from very large scale to quasi-linear and non-linear scales by comparing dark-matter-only simulations with their hydro-dynamical counterparts. At last, we make a brief discussion on the impacts coming from different baryon models and simulation codes.


title: Embracing a new era of highly efficient and productive quantum Monte Carlo simulations

abstract: QMCPACK has enabled cutting-edge materials research on supercomputers for over a decade. It scales nearly ideally but has low single-node efficiency due to the physics-based abstractions using array-of-structures objects, causing inefficient vectorization. We present a systematic approach to transform QMCPACK to better exploit the new hardware features of modern CPUs in portable and maintainable ways. We develop miniapps for fast prototyping and optimizations. We implement new containers in structure-of-arrays data layout to facilitate vectorizations by the compilers. Further speedup and smaller memory-footprints are obtained by computing data on the fly with the vectorized routines and expanding single-precision use. All these are seamlessly incorporated in production QMCPACK. We demonstrate upto 4.5x speedups on recent Intel processors and IBM Blue Gene/Q for representative workloads. Energy consumption is reduced significantly commensurate to the speedup factor. Memory-footprints are reduced by up-to 3.8x, opening the possibility to solve much larger problems of future.


title: Numerical Simulations of Supernova Remnant Evolution in a Cloudy Interstellar Medium

abstract: The mixed morphology class of supernova remnants has centrally peaked X-ray emission along with a shell-like morphology in radio emission. White & Long proposed that these remnants are evolving in a cloudy medium wherein the clouds are evaporated via thermal conduction once being overrun by the expanding shock. Their analytical model made detailed predictions regarding temperature, density and emission profiles as well as shock evolution. We present numerical hydrodynamical models in 2D and 3D including thermal conduction, testing the White & Long model and presenting results for the evolution and emission from remnants evolving in a cloudy medium. We find that, while certain general results of the White & Long model hold, such as the way the remnants expand and the flattening of the X-ray surface brightness distribution, in detail there are substantial differences. In particular we find that the X-ray luminosity is dominated by emission from shocked cloud gas early on, leading to a bright peak which then declines and flattens as evaporation becomes more important. In addition, the effects of thermal conduction on the intercloud gas, which is not included in the White & Long model, are important and lead to further flattening of the X-ray brightness profile as well as lower X-ray emission temperatures.


title: A Review of Self-Exciting Spatio-Temporal Point Processes and Their Applications

abstract: Self-exciting spatio-temporal point process models predict the rate of events as a function of space, time, and the previous history of events. These models naturally capture triggering and clustering behavior, and have been widely used in fields where spatio-temporal clustering of events is observed, such as earthquake modeling, infectious disease, and crime. In the past several decades, advances have been made in estimation, inference, simulation, and diagnostic tools for self-exciting point process models. In this review, I describe the basic theory, survey related estimation and inference techniques from each field, highlight several key applications, and suggest directions for future research.


title: DM-PhyClus: A Bayesian phylogenetic algorithm for infectious disease transmission cluster inference

abstract: Background. Conventional phylogenetic clustering approaches rely on arbitrary cutpoints applied a posteriori to phylogenetic estimates. Although in practice, Bayesian and bootstrap-based clustering tend to lead to similar estimates, they often produce conflicting measures of confidence in clusters. The current study proposes a new Bayesian phylogenetic clustering algorithm, which we refer to as DM-PhyClus, that identifies sets of sequences resulting from quick transmission chains, thus yielding easily-interpretable clusters, without using any ad hoc distance or confidence requirement. Results. Simulations reveal that DM-PhyClus can outperform conventional clustering methods, as well as the Gap procedure, a pure distance-based algorithm, in terms of mean cluster recovery. We apply DM-PhyClus to a sample of real HIV-1 sequences, producing a set of clusters whose inference is in line with the conclusions of a previous thorough analysis. Conclusions. DM-PhyClus, by eliminating the need for cutpoints and producing sensible inference for cluster configurations, can facilitate transmission cluster detection. Future efforts to reduce incidence of infectious diseases, like HIV-1, will need reliable estimates of transmission clusters. It follows that algorithms like DM-PhyClus could serve to better inform public health strategies.


title: $B^0_s \to \ell^+ \ell^- \gamma$ as a Test of Lepton Flavor Universality

abstract: We discuss a number of strategies to reduce the $\mathcal B(B^0_s \to \ell^{+} \ell^{-} \gamma)$ theoretical error, and make such a measurement a new probe of the interactions that are interesting in the light of present-day flavor discrepancies. In particular, for low di-lepton invariant mass we propose to exploit the close parenthood between $\mathcal B(B^0_s \to \ell^{+} \ell^{-} \gamma)$ and the measured $\mathcal B(B^0_s \to \phi (\to K^+ K^-) \gamma)$. For high $q^2$, conversely, we exploit the fact that the decay is dominated by two form-factor combinations, plus contributions from broad charmonium that we model accordingly. We construct the ratio $R_\gamma$, akin to $R_K$ and likewise sensitive to lepton-universality violation. Provided the two rates in this ratio are integrated in a suitable region that minimises bremsstrahlung contributions while maximising statistics, the ratio is very close to unity and the form-factor dependence cancels to an extent that makes it a new valuable probe of lepton-universality violating contributions in the effective Hamiltonian. We finally speculate on additional ideas to extract short-distance information from resonance regions, which are theoretically interesting but statistically limited at present.


title: The Kontsevich-Rosenberg principle for bi-symplectic forms

abstract: In this expository note, we explain the so-called Van den Bergh functor, which enables the formalization of the Kontsevich-Rosenberg principle, whereby a structure on an associative algebra has geometric meaning if it induces standard geometric structures on its representation spaces. Crawley-Boevey, Etingof and Ginzburg proved that bi-symplectic forms satisfy this principle; this implies that bi-symplectic algebras can be regarded as noncommutative symplectic manifolds. In this note, we use the Van den Bergh functor to give an alternative proof.


title: When rule-based models need to count

abstract: Rule-based modelers dislike direct enumeration of cases when more efficient means of enumeration are available. We present an extension of the Kappa language which attaches to agents a notion of level. We detail two encodings that are more concise than the former practice.


title: Density Functional Theory of doped superfluid liquid helium and nanodroplets

abstract: During the last decade, density function theory (DFT) in its static and dynamic time dependent forms, has emerged as a powerful tool to describe the structure and dynamics of doped liquid helium and droplets. In this review, we summarize the activity carried out in this field within the DFT framework since the publication of the previous review article on this subject [M. Barranco et al., J. Low Temp. Phys. 142, 1 (2006)]. Furthermore, a comprehensive presentation of the actual implementations of helium DFT is given, which have not been discussed in the individual articles or are scattered in the existing literature. This is an Accepted Manuscript of an article published on August 2, 2017 by Taylor & Francis Group in Int. Rev. Phys. Chem. 36, 621 (2017), available online: http://dx.doi.org/10.1080/0144235X.2017.1351672


title: On Hilbert's 8th Problem

abstract: A Hadamard factorization of the Riemann Xi-function is constructed to characterize the zeros of the zeta function.


title: Cheryl's Birthday

abstract: We present four logic puzzles and after that their solutions. Joseph Yeo designed 'Cheryl's Birthday'. Mike Hartley came up with a novel solution for 'One Hundred Prisoners and a Light Bulb'. Jonathan Welton designed 'A Blind Guess' and 'Abby's Birthday'. Hans van Ditmarsch and Barteld Kooi authored the puzzlebook 'One Hundred Prisoners and a Light Bulb' that contains other knowledge puzzles, and that can also be found on the webpage http://personal.us.es/hvd/lightbulb.html dedicated to the book.


title: Electro-optical properties of phosphorene quantum dots

abstract: We study electronic and optical properties of single layer phosphorene quantum dots with various shapes, sizes, and edge types (including disordered edges) subjected to an external electric field normal to the structure plane. Compared to graphene quantum dots, in phosphorene clusters of similar shape and size there is a set of edge states with energies dispersed at around the Fermi level. These states make the majority of phosphorene quantum dots metallic and enrich the phosphorene absorption gap with low-energy absorption peaks tunable by the electric field. The presence of the edge states dispersed at around the Fermi level is a characteristic feature that is independent of the edge morphology and roughness.


title: Constructing elliptic curves from Galois representations

abstract: Given a non-isotrivial elliptic curve over an arithmetic surface, one obtains a lisse $\ell$-adic sheaf of rank two over the surface. This lisse sheaf has a number of straightforward properties: cyclotomic determinant, finite ramification, rational traces of Frobenius, and somewhere not potentially good reduction. We prove that any lisse sheaf of rank two possessing these properties comes from an elliptic curve.


title: Which Encoding is the Best for Text Classification in Chinese, English, Japanese and Korean?

abstract: This article offers an empirical study on the different ways of encoding Chinese, Japanese, Korean (CJK) and English languages for text classification. Different encoding levels are studied, including UTF-8 bytes, characters, words, romanized characters and romanized words. For all encoding levels, whenever applicable, we provide comparisons with linear models, fastText and convolutional networks. For convolutional networks, we compare between encoding mechanisms using character glyph images, one-hot (or one-of-n) encoding, and embedding. In total there are 473 models, using 14 large-scale text classification datasets in 4 languages including Chinese, English, Japanese and Korean. Some conclusions from these results include that byte-level one-hot encoding based on UTF-8 consistently produces competitive results for convolutional networks, that word-level n-grams linear models are competitive even without perfect word segmentation, and that fastText provides the best result using character-level n-gram encoding but can overfit when the features are overly rich.


title: The universal connection for principal bundles over homogeneous spaces and twistor space of coadjoint orbits

abstract: Given a holomorphic principal bundle $Q\, \longrightarrow\, X$, the universal space of holomorphic connections is a torsor $C_1(Q)$ for $\text{ad} Q \otimes T^*X$ such that the pullback of $Q$ to $C_1(Q)$ has a tautological holomorphic connection. When $X\,=\, G/P$, where $P$ is a parabolic subgroup of a complex simple group $G$, and $Q$ is the frame bundle of an ample line bundle, we show that $C_1(Q)$ may be identified with $G/L$, where $L\, \subset\, P$ is a Levi factor. We use this identification to construct the twistor space associated to a natural hyper-K\"ahler metric on $T^*(G/P)$, recovering Biquard's description of this twistor space, but employing only finite-dimensional, Lie-theoretic means.


title: Gramian Tensor Decomposition via Semidefinite Programming

abstract: In this paper we examine a symmetric tensor decomposition problem, the Gramian decomposition, posed as a rank minimization problem. We study the relaxation of the problem and consider cases when the relaxed solution is a solution to the original problem. In some instances of tensor rank and order, we prove generically that the solution to the relaxation will be optimal in the original. In other cases, we present interesting examples and approaches that demonstrate the intricacy of this problem.


title: Learning Visual Importance for Graphic Designs and Data Visualizations

abstract: Knowing where people look and click on visual designs can provide clues about how the designs are perceived, and where the most important or relevant content lies. The most important content of a visual design can be used for effective summarization or to facilitate retrieval from a database. We present automated models that predict the relative importance of different elements in data visualizations and graphic designs. Our models are neural networks trained on human clicks and importance annotations on hundreds of designs. We collected a new dataset of crowdsourced importance, and analyzed the predictions of our models with respect to ground truth importance and human eye movements. We demonstrate how such predictions of importance can be used for automatic design retargeting and thumbnailing. User studies with hundreds of MTurk participants validate that, with limited post-processing, our importance-driven applications are on par with, or outperform, current state-of-the-art methods, including natural image saliency. We also provide a demonstration of how our importance predictions can be built into interactive design tools to offer immediate feedback during the design process.


title: Theory and simulation of anode spots in low pressure plasmas

abstract: When electrodes are biased above the plasma potential, electrons accelerated through the associated electron sheath can dramatically increase the ionization rate of neutrals near the electrode surface. It has previously been observed that if the ionization rate is great enough, a double layer separates a luminous high-potential plasma attached to the electrode surface (called an anode spot or fireball) from the bulk plasma. Here, results of the first 2D particle-in-cell simulations of anode spot formation are presented along with a theoretical model describing the formation process. It is found that ionization leads to the buildup of an ion-rich layer adjacent to the electrode, forming a narrow potential well near the electrode surface that traps electrons born from ionization. Anode spot onset occurs when a quasineutral region is established in the potential well and the density in this region becomes large enough to violate the steady-state Langmuir condition, which is a balance between electron and ion fluxes across the double layer. A model for steady-state properties of the anode spot is also presented, which predicts values for the anode spot size, double layer potential drop, and form of the sheath at the electrode by considering particle, power, and current balance. These predictions are found to be consistent with the presented simulation and previous experiments.


title: Online unit clustering in higher dimensions

abstract: We revisit the online Unit Clustering and Unit Covering problems in higher dimensions: Given a set of $n$ points in a metric space, that arrive one by one, Unit Clustering asks to partition the points into the minimum number of clusters (subsets) of diameter at most one; while Unit Covering asks to cover all points by the minimum number of balls of unit radius. In this paper, we work in $\mathbb{R}^d$ using the $L_\infty$ norm. We show that the competitive ratio of any online algorithm (deterministic or randomized) for Unit Clustering must depend on the dimension $d$. We also give a randomized online algorithm with competitive ratio $O(d^2)$ for Unit Clustering of integer points (i.e., points in $\mathbb{Z}^d$, $d\in \mathbb{N}$, under $L_{\infty}$ norm). We show that the competitive ratio of any deterministic online algorithm for Unit Covering is at least $2^d$. This ratio is the best possible, as it can be attained by a simple deterministic algorithm that assigns points to a predefined set of unit cubes. We complement these results with some additional lower bounds for related problems in higher dimensions.


title: Gradient-enhanced kriging for high-dimensional problems

abstract: Surrogate models provide a low computational cost alternative to evaluating expensive functions. The construction of accurate surrogate models with large numbers of independent variables is currently prohibitive because it requires a large number of function evaluations. Gradient-enhanced kriging has the potential to reduce the number of function evaluations for the desired accuracy when efficient gradient computation, such as an adjoint method, is available. However, current gradient-enhanced kriging methods do not scale well with the number of sampling points due to the rapid growth in the size of the correlation matrix where new information is added for each sampling point in each direction of the design space. They do not scale well with the number of independent variables either due to the increase in the number of hyperparameters that needs to be estimated. To address this issue, we develop a new gradient-enhanced surrogate model approach that drastically reduced the number of hyperparameters through the use of the partial-least squares method that maintains accuracy. In addition, this method is able to control the size of the correlation matrix by adding only relevant points defined through the information provided by the partial-least squares method. To validate our method, we compare the global accuracy of the proposed method with conventional kriging surrogate models on two analytic functions with up to 100 dimensions, as well as engineering problems of varied complexity with up to 15 dimensions. We show that the proposed method requires fewer sampling points than conventional methods to obtain the desired accuracy, or provides more accuracy for a fixed budget of sampling points. In some cases, we get over 3 times more accurate models than a bench of surrogate models from the literature, and also over 3200 times faster than standard gradient-enhanced kriging models.


title: Internet of Tangible Things (IoTT): Challenges and Opportunities for Tangible Interaction with IoT

abstract: In the Internet of Things era, an increasing number of household devices and everyday objects are able to send to and retrieve information from the Internet, offering innovative services to the user. However, most of these devices provide only smartphone or web interfaces to control the IoT object properties and functions. As a result, generally, the interaction is disconnected from the physical world, decreasing the user experience and increasing the risk of isolating the user in digital bubbles. We argue that tangible interaction can counteract this trend and this paper discusses the potential benefits and the still open challenges of tangible interaction applied to the Internet of Things. To underline this need, we introduce the term Internet of Tangible Things. In the article, after an analysis of current open challenges for Human-Computer Interaction in IoT, we summarize current trends in tangible interaction and extrapolate eight tangible interaction properties that could be exploited for designing novel interactions with IoT objects. Through a systematic literature review of tangible interaction applied to IoT, we show what has been already explored in the systems that pioneered the field and the future explorations that still have to be conducted.


title: Stable rank for crossed products by actions of finite groups on C*-algebras

abstract: Let $G$ be a finite group, $A$ a unital separable finite simple nuclear C*-algebra, and $\alpha$ an action of $G$ on $A$. Assume that $A$ absorbs the Jiang-Su algebra $\mathcal{Z}$, the extremal boundary of the trace space of $A$ is compact and finite dimensional and that $\alpha$ fixes any tracial state of $A$. Then tsr$(A \rtimes_\alpha G) = 1$. In particular, when $A$ has a unique tracial state, we conclude it without above conditons on a tracial state space of $A$.


abstract: This is the Proceedings of the 2017 ICML Workshop on Human Interpretability in Machine Learning (WHI 2017), which was held in Sydney, Australia, August 10, 2017. Invited speakers were Tony Jebara, Pang Wei Koh, and David Sontag.


title: Evolutionary advantage of directional symmetry breaking in self-replicating polymers

abstract: Due to the asymmetric nature of the nucleotides, the extant informational biomolecule, DNA, is constrained to replicate unidirectionally on a template. As a product of molecular evolution that sought to maximize replicative potential, DNA's unidirectional replication poses a mystery since symmetric bidirectional self-replicators obviously would replicate faster than unidirectional self-replicators and hence would have been evolutionarily more successful. Here we carefully examine the physico-chemical requirements for evolutionarily successful primordial self-replicators and theoretically show that at low monomer concentrations that possibly prevailed in the primordial oceans, asymmetric unidirectional self-replicators would have an evolutionary advantage over bidirectional self-replicators. The competing requirements of low and high kinetic barriers for formation and long lifetime of inter-strand bonds respectively are simultaneously satisfied through asymmetric kinetic influence of inter-strand bonds, resulting in evolutionarily successful unidirectional self-replicators.


title: A discriminative view of MRF pre-processing algorithms

abstract: While Markov Random Fields (MRFs) are widely used in computer vision, they present a quite challenging inference problem. MRF inference can be accelerated by pre-processing techniques like Dead End Elimination (DEE) or QPBO-based approaches which compute the optimal labeling of a subset of variables. These techniques are guaranteed to never wrongly label a variable but they often leave a large number of variables unlabeled. We address this shortcoming by interpreting pre-processing as a classification problem, which allows us to trade off false positives (i.e., giving a variable an incorrect label) versus false negatives (i.e., failing to label a variable). We describe an efficient discriminative rule that finds optimal solutions for a subset of variables. Our technique provides both per-instance and worst-case guarantees concerning the quality of the solution. Empirical studies were conducted over several benchmark datasets. We obtain a speedup factor of 2 to 12 over expansion moves without preprocessing, and on difficult non-submodular energy functions produce slightly lower energy.


title: A differential Least Squares Deconvolution method for high precision spectroscopy of stars and exoplanets I. Application to obliquity measurements of HARPS observations of HD189733b

abstract: High precision measurements of stellar spectroscopic line profiles and their changes over time contain very valuable information about the physics of the stellar photosphere (stellar activity) and can be used to characterize extrasolar planets via the Rossiter-McLaughlin effect or from reflected light from the planet. In this paper we present a new method for measuring small changes in the mean line profile of a spectrum by performing what we call differential Least Squares Deconvolution (dLSD). The method consists in finding the convolution function (or kernel) required to transform a high signal-to-noise ratio template of the star into each observed spectrum. Compared to similar techniques, the method presented here does not require any assumptions on the template spectrum (eg. no line-list or cross-correlation mask required). We show that our implementation of dLSD is able to perform -at least- as good as other techniques by applying it to star-planet obliquity measurements of exoplanet HD183799 during its transit. Among other things, the method should enable model independent detection of light reflected by an exoplanet.


title: H\"older continuity of the integrated density of states for Extended Harper's Model with Liouville frequency

abstract: In this paper, we study the non-self-dual extended Harper's model with a Liouville frequency. Based on the work of \cite{SY}, we show that the integrated density of states (IDS for short) of the model is $\frac{1}{2}$-H$\ddot{\text{o}}$lder continuous. As an application, we also obtain the Carleson homogeneity of the spectrum.


title: Variants of the Riemann zeta function

abstract: We construct variants of the Riemann zeta function with convenient properties and make conjectures about their dynamics; some of the conjectures are based on an analogy with the dynamical system of zeta. More specifically, we study the family of functions $V_z: s \mapsto \zeta(s) \exp (zs)$. We observe convergence of $V_z$ fixed points along nearly logarithmic spirals with initial points at zeta fixed points and centered upon Riemann zeros. We can approximate these spirals numerically, so they might afford a means to study the geometry of the relationship of zeta fixed points to Riemann zeros.


title: Properties of size-dependent models having quasiperiodic boundary conditions

abstract: Boundary conditions effects are explored for size-dependent models in thermal equilibrium. Scalar and fermionic models are used for $D=1+3$ (films), $D=1+2$ (hollow cylinder) and $D=1+1$ (ring). For all models a minimal length is found, below which no thermally-induced phase transition occurs. Using quasiperiodic boundary condition controlled by a contour parameter $\theta$ ($\theta=0$ is a periodic boundary condition and $\theta=1$ is an antiperiodic condition) it results that the minimal length depends directly on the value of $\theta$. It is also argued that this parameter can be associated to an Aharonov-Bohm phase.


title: Regular versus singular order of contact on pseudoconvex hypersurfaces

abstract: The singular and regular type of a point on a real hypersurface $\mathcal H$ in $\mathbb C^n$ are shown to agree when the regular type is strictly less than 4. If $\mathcal H$ is pseudoconvex, we show they agree when the regular type is 4. A non-pseudoconvex example is given where the regular type is 4 and the singular type is infinite.


title: Tensor-Structured Coupled Cluster Theory

abstract: We derive and implement a new way of solving coupled cluster equations with lower computational scaling. Our method is based on decomposition of both amplitudes and two electron integrals, using a combination of tensor hypercontraction and canonical polyadic decomposition. While the original theory scales as $O(N^6)$ with respect to the number of basis functions, we demonstrate numerically that we achieve sub-millihartree difference from the original theory with $O(N^4)$ scaling. This is accomplished by solving directly for the factors that decompose the cluster operator. The proposed scheme is quite general and can be easily extended to other many-body methods.


title: Inheritance of Convexity for the $\mathcal{P}_{\min}$-Restricted Game

abstract: We consider restricted games on weighted graphs associated with minimum partitions. We replace in the classical definition of Myerson restricted game the connected components of any subgraph by the subcomponents corresponding to a minimum partition. This minimum partition $\mathcal{P}_{\min}$ is induced by the deletion of the minimum weight edges. We provide a characterization of the graphs satisfying inheritance of convexity from the underlying game to the restricted game associated with Pmin. Moreover, we prove that these graphs can be recognized in polynomial time.


title: Power packet transferability via symbol propagation matrix

abstract: Power packet is a unit of electric power transferred by a power pulse with an information tag. In Shannon's information theory, messages are represented by symbol sequences in a digitized manner. Referring to this formulation, we define symbols in power packetization as a minimum unit of power transferred by a tagged pulse. Here, power is digitized and quantized. In this paper, we consider packetized power in networks for a finite duration, giving symbols and their energies to the networks. A network structure is defined using a graph whose nodes represent routers, sources, and destinations. First, we introduce symbol propagation matrix (SPM) in which symbols are transferred at links during unit times. Packetized power is described as a network flow in a spatio-temporal structure. Then, we study the problem of selecting an SPM in terms of transferability, that is, the possibility to represent given energies at sources and destinations during the finite duration. To select an SPM, we consider a network flow problem of packetized power. The problem is formulated as an M-convex submodular flow problem which is known as generalization of the minimum cost flow problem and solvable. Finally, through examples, we verify that this formulation provides reasonable packetized power.


title: Randomly coloring graphs of bounded treewidth

abstract: We consider the problem of sampling a proper $k$-coloring of a graph of maximal degree $\Delta$ uniformly at random. We describe a new Markov chain for sampling colorings, and show that it mixes rapidly on graphs of bounded treewidth if $k\geq(1+\epsilon)\Delta$, for any $\epsilon>0$.


title: DirectDM: a tool for dark matter direct detection

abstract: We provide a Mathematica package, DirectDM, that takes as input the Wilson coefficients of the relativistic effective theory describing the interactions of dark matter with quarks, gluons and photons, and matches it onto an effective theory describing the interactions of dark matter with neutrons and protons. The nonperturbative matching is performed at leading order in a chiral expansion. The one-loop QCD and QED renormalization-group evolution from the electroweak scale down to the hadronic scale, as well as finite corrections at the heavy quark thresholds are taken into account. We also provide an interface with the package DMFormFactor so that, starting from the relativistic effective theory, one can directly obtain the event rates for direct detection experiments.


title: Heavy-Quark Hybrid Mass Splittings: Hyperfine and "Ultrafine"

abstract: It is argued that the heavy-quark limit of QCD requires a certain combination of hyperfine mass splittings in heavy-quark hybrid-meson multiplets to be unusually small. This observation will assist in the exploration of the heavy-quark hybrid spectrum at facilities such as PANDA. Alternatively, a large measured value for this mass splitting indicates that at least one member of the multiplet must contain significant light-quark degrees of freedom.


title: Extreme infrared variables from UKIDSS - II. an end-of-survey catalogue of eruptive YSOs and unusual stars

abstract: We present a catalogue of 618 high amplitude infrared variable stars (1 < DeltaK < 5 mag) detected by the two widely separated epochs of 2.2 um data in the UKIDSS Galactic plane survey, from searches covering 1470 deg^2. Most were discovered by a search of all fields at 30 < l < 230 deg. Sources include new dusty Mira variables, three new CV candidates, a blazar and a peculiar source that may be an interacting binary system. However, c. 60 per cent are YSOs, based on spatial association with star forming regions at distances ranging from 300 pc to over 10 kpc. This confirms our initial result in Contreras Pena et al.(Paper I) that YSOs dominate the high amplitude infrared variable sky in the Galactic disc. It is also supported by recently published VVV results at 295 < l < 350 deg. The spectral energy distributions of the YSOs indicate class I or flat spectrum systems in most cases, as in the VVV sample. A large number of variable YSOs are associated with the Cygnus X complex and other groups are associated with the North America/Pelican nebula, the Gemini OB1 molecular cloud, the Rosette complex, the Cone nebula, the W51 star forming region and the S86 and S236 HII regions. Most of the YSO variability is likely due to variable/episodic accretion on timescales of years, albeit usually less extreme than classical FUors and EXors. Luminosities at the 2010 WISE epoch range from c. 0.1 Lsun to 10^3 Lsun but only rarely exceed 10^2.5 Lsun.


title: Generative Adversarial Network-based Synthesis of Visible Faces from Polarimetric Thermal Faces

abstract: The large domain discrepancy between faces captured in polarimetric (or conventional) thermal and visible domain makes cross-domain face recognition quite a challenging problem for both human-examiners and computer vision algorithms. Previous approaches utilize a two-step procedure (visible feature estimation and visible image reconstruction) to synthesize the visible image given the corresponding polarimetric thermal image. However, these are regarded as two disjoint steps and hence may hinder the performance of visible face reconstruction. We argue that joint optimization would be a better way to reconstruct more photo-realistic images for both computer vision algorithms and human-examiners to examine. To this end, this paper proposes a Generative Adversarial Network-based Visible Face Synthesis (GAN-VFS) method to synthesize more photo-realistic visible face images from their corresponding polarimetric images. To ensure that the encoded visible-features contain more semantically meaningful information in reconstructing the visible face image, a guidance sub-network is involved into the training procedure. To achieve photo realistic property while preserving discriminative characteristics for the reconstructed outputs, an identity loss combined with the perceptual loss are optimized in the framework. Multiple experiments evaluated on different experimental protocols demonstrate that the proposed method achieves state-of-the-art performance.


title: Investigation of Turbulent Mixing and Local Reaction Rates on Deflagration to Detonation Transition in Methane-Oxygen

abstract: In the current study, the influence of turbulent mixing and local reaction rates on deflagration to detonation transition (DDT) was investigated using a state-of-the-art large eddy simulation (LES) strategy. Specifically, detonation attenuation by a porous medium, and the subsequent re-initiation for methane-oxygen, a moderately unstable mixture, was considered. The purpose of the investigation was to validate the numerical strategy with previous experimental observations, and to determine what specific roles turbulent mixing and shock compression have on flame acceleration during the final stages of DDT. The modelling procedure adopted was a grid-within-a-grid approach: The compressible linear eddy model for large eddy simulation (CLEM-LES). It was found that average turbulent fluctuations greater than the laminar flame speed were required in order to maintain wave velocities above the Chapman-Jouguet (CJ)-deflagration velocity threshold, a precursor requirement for DDT to occur. It was also found that sufficient turbulent burning on the flame surface was required in order to drive pressure waves to sufficiently strengthen the leading shock wave, locally, in order to trigger auto-ignition hot spots in the wave front. These local explosion events, which were found to burn out through turbulent surface reactions, drive transverse pressure waves outward. Upon subsequent shock reflections or interactions of the transverse waves, new local explosion events occurred, which further strengthened the adjacent leading shock wave above the CJ-detonation speed. Eventually, through this process, the wave sustained the CJ-detonation speed, on average, through the cyclic mechanism of local explosion events followed by turbulent surface reactions. Finally, combustion of the flame acceleration process was found to lie predominantly within the thin-reaction zones regime.


title: Probing the Physics of Narrow Line Regions in Active Galaxies IV: Full data release of The Siding Spring Southern Seyfert Spectroscopic Snapshot Survey (S7)

abstract: We present the second and final data release of the Siding Spring Southern Seyfert Spectroscopic Snapshot Survey (S7). Data are presented for 63 new galaxies not included in the first data release, and we provide 2D emission-line fitting products for the full S7 sample of 131 galaxies. The S7 uses the WiFeS instrument on the ANU 2.3m telescope to obtain spectra with a spectral resolution of R = 7000 in the red (540 - 700 nm) and R = 3000 in the blue (350 - 570 nm), over an integral field of 25x38 arcsec^2 with 1x1 arcsec^2 spatial pixels. The S7 contains both the largest sample of active galaxies and the highest spectral resolution of any comparable integral field survey to date. The emission-line fitting products include line fluxes, velocities and velocity dispersions across the WiFeS field of view, and an artificial neural network has been used to determine the optimal number of Gaussian kinematic components for emission-lines in each spaxel. Broad Balmer lines are subtracted from the spectra of nuclear spatial pixels in Seyfert 1 galaxies before fitting the narrow lines. We bin nuclear spectra and measure reddening-corrected nuclear fluxes of strong narrow lines for each galaxy. The nuclear spectra are classified on optical diagnostic diagrams, where the strength of the coronal line [FeVII]6087 is shown to be correlated with [OIII]/Hbeta. Maps revealing gas excitation and kinematics are included for the entire sample, and we provide notes on the newly observed objects.


title: Overdense microwave plasma heating in the CNT stellarator

abstract: Overdense plasmas have been attained with 2.45 GHz microwave heating in the low-field, low-aspect-ratio CNT stellarator. Densities higher than four times the ordinary (O) mode cutoff density were measured with 8 kW of power injected in the O-mode and, alternatively, with 6.5 kW in the extraordinary (X) mode. The temperature profiles peak at the plasma edge. This was ascribed to collisional damping of the X-mode at the upper hybrid resonant layer. The X-mode reaches that location by tunneling, mode-conversions or after polarization-scrambling reflections off the wall and in-vessel coils, regardless of the initial launch being in O- or X-mode. This interpretation was confirmed by full-wave numerical simulations. Also, as the CNT plasma is not completely ionized at these low microwave power levels, electron density was shown to increase with power. A dependence on magnetic field strength was also observed (for O-mode launch) and discussed.


title: Data driven modal decompositions: analysis and enhancements

abstract: The Dynamic Mode Decomposition (DMD) is a tool of trade in computational data driven analysis of fluid flows. More generally, it is a computational device for Koopman spectral analysis of nonlinear dynamical systems, with a plethora of applications in applied sciences and engineering. Its exceptional performance triggered developments of several modifications that make the DMD an attractive method in data driven framework. This work offers further improvements of the DMD to make it more reliable, and to enhance its functionality. In particular, data driven formula for the residuals allows selection of the Ritz pairs, thus providing more precise spectral information of the underlying Koopman operator, and the well-known technique of refining the Ritz vectors is adapted to data driven scenarios. Further, the DMD is formulated in a more general setting of weighted inner product spaces, and the consequences for numerical computation are discussed in detail. Numerical experiments are used to illustrate the advantages of the proposed method, designated as DDMD_RRR (Refined Rayleigh Ritz Data Driven Modal Decomposition).


title: Isolation effects in a system of two mutually communicating identical patches

abstract: Starting from the Fisher-Kolmogorov-Petrovskii-Piskunov equation (FKPP) we model the dynamic of a diffusive system with two mutually communicating identical patches and isolated of the remaining matrix. For this system we find the minimal size of each fragment in the explicit form and compare with the explicit results for similar problems found in the literature. From this comparison emerges an unexpected result that for a same set of the parameters, the isolated system studied in this work with size L, can be better or worst than the non isolated systems with the same size L, uniquely depending on the parameter $a_{0}$ (internal conditions of the patches). Due to the fact that this result is unexpected we propose a experimental verification.


title: The dependence of the IR-radio correlation on the metallicity

abstract: We have compiled a sample of 26 metal-poor galaxies with 12 + log(O/H) < 8.1 with both infrared continuum and 1.4 GHz radio continuum data. By comparing to galaxies at higher metallicity, we have investigated the dependence on the metallicity of the IR-radio relationship at 24 um, 70 um, 100 um and 160 um bands as well as the integrated FIR luminosity. It is found that metal-poor galaxies have on average lower qIR than metal-rich ones with larger offsets at longer IR wavelengths, from -0.06 dex in q24um to -0.6 dex in q160um. The qIR of all galaxies as a whole at 160 um show positive trends with the metallicity and IR-to-FUV ratio, and negative trends with the IR color, while those at lower IR wavelengths show weaker correlations. We proposed a mechanism that invokes combined effects of low obscured-SFR/total-SFR fraction and warm dust temperature at low metallicity to interpret the above behavior of qIR, with the former reducing the IR radiation and the latter further reducing the IR emission at longer IR wavelength. Other mechanisms that are related to the radio emission including the enhanced magnetic field strength and increased thermal radio contribution are unable to reconcile the IR-wavelength-dependent differences of qIR between metal-poor and metal- rich galaxies. In contrast to qIR, the mean total-SFR/radio ratio of metal-poor galaxies is the same as the metal-rich one, indicating the 1.4 GHz radio emission is still an effective tracer of SFRs at low metallicity.


title: Statistics of Deep Generated Images

abstract: Here, we explore the low-level statistics of images generated by state-of-the-art deep generative models. First, Variational auto-encoder (VAE~\cite{kingma2013auto}), Wasserstein generative adversarial network (WGAN~\cite{arjovsky2017wasserstein}) and deep convolutional generative adversarial network (DCGAN~\cite{radford2015unsupervised}) are trained on the ImageNet dataset and a large set of cartoon frames from animations. Then, for images generated by these models as well as natural scenes and cartoons, statistics including mean power spectrum, the number of connected components in a given image area, distribution of random filter responses, and contrast distribution are computed. Our analyses on training images support current findings on scale invariance, non-Gaussianity, and Weibull contrast distribution of natural scenes. We find that although similar results hold over cartoon images, there is still a significant difference between statistics of natural scenes and images generated by VAE, DCGAN and WGAN models. In particular, generated images do not have scale invariant mean power spectrum magnitude, which indicates existence of extra structures in these images. Inspecting how well the statistics of deep generated images match the known statistical properties of natural images, such as scale invariance, non-Gaussianity, and Weibull contrast distribution, can a) reveal the degree to which deep learning models capture the essence of the natural scenes, b) provide a new dimension to evaluate models, and c) allow possible improvement of image generative models (e.g., via defining new loss functions).


title: Criteria of existence for bounce solutions in false vacuum decay with gravity

abstract: The bounce solutions of self-interacting scalar fields coupled to gravity are studied using a semi-classical approach. We found that bounce solutions have a maximum required barrier curvature, in addition to the known minimum required barrier curvature. In particular, as the maximum barrier curvature is approached, the scale factor of the well-known Colemen-De Luccia (CDL) bounce solutions become divergent. Unlike the CDL or its more general oscillating bounce counterparts, this cannot be considered as a subset of the Hawking-Turok solution.


title: On the number of proper paths between vertices in edge-colored hypercubes

abstract: Given an integer $1\leq j <n$, define the $(j)$-coloring of a $n$-dimensional hypercube $H_{n}$ to be the $2$-coloring of the edges of $H_{n}$ in which all edges in dimension $i$, $1\leq i \leq j$, have color $1$ and all other edges have color $2$. Cheng et al. [Proper distance in edge-colored hypercubes, Applied Mathematics and Computation 313 (2017) 384-391.] determined the number of distinct shortest properly colored paths between a pair of vertices for the $(1)$-colored hypercubes. It is natural to consider the number for $(j)$-coloring, $j\geq 2$. In this note, we determine the number of different shortest proper paths in $(j)$-colored hypercubes for arbitrary $j$.


title: Universal Function Approximation by Deep Neural Nets with Bounded Width and ReLU Activations

abstract: This article concerns the expressive power of depth in neural nets with ReLU activations and bounded width. We are particularly interested in the following questions: what is the minimal width $w_{\text{min}}(d)$ so that ReLU nets of width $w_{\text{min}}(d)$ (and arbitrary depth) can approximate any continuous function on the unit cube $[0,1]^d$ aribitrarily well? For ReLU nets near this minimal width, what can one say about the depth necessary to approximate a given function? Our approach to this paper is based on the observation that, due to the convexity of the ReLU activation, ReLU nets are particularly well-suited for representing convex functions. In particular, we prove that ReLU nets with width $d+1$ can approximate any continuous convex function of $d$ variables arbitrarily well. These results then give quantitative depth estimates for the rate of approximation of any continuous scalar function on the $d$-dimensional cube $[0,1]^d$ by ReLU nets with width $d+3.$

categories: stat.ML cs.CG cs.LG math.FA math.ST stat.TH


title: A kind of conditional connectivity of transposition networks generated by $k$-trees

abstract: For a graph $G = (V, E)$, a subset $F\subset V(G)$ is called an $R_k$-vertex-cut of $G$ if $G -F$ is disconnected and each vertex $u \in V(G)- F$ has at least $k$ neighbors in $G -F$. The $R_k$-vertex-connectivity of $G$, denoted by $\kappa^k(G)$, is the cardinality of the minimum $R_k$-vertex-cut of $G$, which is a refined measure for the fault tolerance of network $G$. In this paper, we study $\kappa^2$ for Cayley graphs generated by $k$-trees. Let $Sym(n)$ be the symmetric group on $\{1, 2, \cdots ,n\}$ and $\mathcal{T}$ be a set of transpositions of $Sym(n)$. Let $G(\mathcal{T})$ be the graph on $n$ vertices $\{1, 2, . . . ,n\}$ such that there is an edge $ij$ in $G(\mathcal{T})$ if and only if the transposition $ij\in \mathcal{T}$. The graph $G(\mathcal{T})$ is called the transposition generating graph of $\mathcal{T}$. We denote by $Cay(Sym(n),\mathcal{T})$ the Cayley graph generated by $G(\mathcal{T})$. The Cayley graph $Cay(Sym(n),\mathcal{T})$ is denoted by $T_kG_n$ if $G(\mathcal{T})$ is a $k$-tree. We determine $\kappa^2(T_kG_n)$ in this work. The trees are $1$-trees, and the complete graph on $n$ vertices is a $n-1$-tree. Thus, in this sense, this work is a generalization of the such results on Cayley graphs generated by transposition generating trees and the complete-transposition graphs.


title: Generalized Entropy Agglomeration

abstract: Entropy Agglomeration (EA) is a hierarchical clustering algorithm introduced in 2013. Here, we generalize it to define Generalized Entropy Agglomeration (GEA) that can work with multiset blocks and blocks with rational occurrence numbers. We also introduce a numerical categorization procedure to apply GEA to numerical datasets. The software REBUS 2.0 is published with these capabilities: http://fidaner.wordpress.com/science/rebus2


title: Human Skin Detection Using RGB, HSV and YCbCr Color Models

abstract: Human Skin detection deals with the recognition of skin-colored pixels and regions in a given image. Skin color is often used in human skin detection because it is invariant to orientation and size and is fast to process. A new human skin detection algorithm is proposed in this paper. The three main parameters for recognizing a skin pixel are RGB (Red, Green, Blue), HSV (Hue, Saturation, Value) and YCbCr (Luminance, Chrominance) color models. The objective of proposed algorithm is to improve the recognition of skin pixels in given images. The algorithm not only considers individual ranges of the three color parameters but also takes into ac- count combinational ranges which provide greater accuracy in recognizing the skin area in a given image.


title: Global well-posedness for the 2D Boussinesq equations with a velocity damping term

abstract: In this paper, we prove global well-posedness of smooth solutions to the two-dimensional incompressible Boussinesq equations with only a velocity damping term when the initial data is close to an nontrivial equilibrium state $(0,x_2)$. As a by-product, under this equilibrium state, our result gives a positive answer to the question proposed by [ACWX] (see P.3597).


title: What Actions are Needed for Understanding Human Actions in Videos?

abstract: What is the right way to reason about human activities? What directions forward are most promising? In this work, we analyze the current state of human activity understanding in videos. The goal of this paper is to examine datasets, evaluation metrics, algorithms, and potential future directions. We look at the qualitative attributes that define activities such as pose variability, brevity, and density. The experiments consider multiple state-of-the-art algorithms and multiple datasets. The results demonstrate that while there is inherent ambiguity in the temporal extent of activities, current datasets still permit effective benchmarking. We discover that fine-grained understanding of objects and pose when combined with temporal reasoning is likely to yield substantial improvements in algorithmic accuracy. We present the many kinds of information that will be needed to achieve substantial gains in activity understanding: objects, verbs, intent, and sequential reasoning. The software and additional information will be made available to provide other researchers detailed diagnostics to understand their own algorithms.


title: High-resolution and Fully-programmable Microwave-shaper

abstract: We here propose and demonstrate a point-by-point programmable broadband microwave spectrum processor with high-resolution up to tens of MHz. We achieve this by bandwidth-minified mapping a programmable optical spectrum processor, which has much larger bandwidth and lower frequency resolution, to a microwave one, and make sure the mapping is a similarity transformation. As a comparison, the traditional optical-to-microwave mapping based on super-heterodyne is bandwidth-preserved and identical. Here the optical spectrum processor is firstly sliced by a high-quality-factor optical resonator, which has periodic transmission peak with MHz-level bandwidth, and then is mapped to microwave domain by optical-frequency-comb-assisted multi-heterodyne. We demonstrate the high frequency resolution and full function reconfigurability experimentally. In a numerical example, we show that the group delay variation of optical spectrum processor could be greatly enlarged after mapping. The resolution improvement and group delay magnification distinguish our proposal significantly from previous spectrum mapping from optics to microwave.


title: On degenerations of projective varieties to complexity-one T-varieties

abstract: Let $R$ be a positively graded finitely generated $\textbf{k}$-domain with Krull dimension $d+1$. We show that there is a homogeneous valuation $\mathfrak{v}: R \setminus \{0\} \to \mathbb{Z}^d$ of rank $d$ such that the associated graded $\text{gr}_\mathfrak{v}(R)$ is finitely generated. This then implies that any polarized $d$-dimensional projective variety $X$ has a flat deformation over $\mathbb{A}^1$, with reduced and irreducible fibers, to a polarized projective complexity-one $T$-variety (i.e. a variety with a faithful action of a $(d-1)$-dimensional torus $T$). As an application we conclude that any $d$-dimensional complex smooth projective variety $X$ equipped with an integral K\"ahler form has a proper $(d-1)$-dimensional Hamiltonian torus action on an open dense subset that extends continuously to all of $X$.


title: Pressure induced switching in ferroelectrics: on the junction between physics and electrochemistry

abstract: Pressure-induced polarization switching in ferroelectric thin films has emerged as a powerful method for domain patterning, allowing to create predefined domain patterns on free surfaces and under thin conductive top electrodes. However, the mechanisms for pressure induced polarization switching in ferroelectrics remain highly controversial, with flexoelectricity, polarization rotation and suppression, and bulk and surface electrochemical processes all being potentially relevant. Here we classify possible pressure induced switching mechanisms, perform elementary estimates, and study in depth using phase-field modelling. We show that magnitudes of these effects are remarkably close, and give rise to complex switching diagrams as a function of pressure and film thickness with non-trivial topology or switchable and non-switchable regions.


title: Current partition at zero-line intersection of quantum anomalous Hall topologies

abstract: At the interface between two-dimensional materials with different topologies, topologically protected one-dimensional states (also named as zero-line modes) arise. Here, we focus on the quantum anomalous Hall effect based zero-line modes formed at the interface between regimes with different Chern numbers. We find that, these zero-line modes are chiral and unilaterally conductive due to the breaking of time-reversal invariance. For a beam splitter consisting of two intersecting zero lines, the chirality ensures that current can only be injected from two of the four terminals. Our numerical results further show that, in the absence of contact resistance, the (anti-)clockwise partitions of currents from these two terminals are the same owing to the current conservation, which effectively simplifies the partition laws. We find that the partition is robust against relative shift of Fermi energy, but can be effectively adjusted by tuning the relative magnetization strengths at different regimes or relative angles between zero lines.


title: Sparse operator compression of higher-order elliptic operators with rough coefficients

abstract: We introduce the sparse operator compression to compress a self-adjoint higher-order elliptic operator with rough coefficients and various boundary conditions. The operator compression is achieved by using localized basis functions, which are energy-minimizing functions on local patches. On a regular mesh with mesh size $h$, the localized basis functions have supports of diameter $O(h\log(1/h))$ and give optimal compression rate of the solution operator. We show that by using localized basis functions with supports of diameter $O(h\log(1/h))$, our method achieves the optimal compression rate of the solution operator. From the perspective of the generalized finite element method to solve elliptic equations, the localized basis functions have the optimal convergence rate $O(h^k)$ for a $(2k)$th-order elliptic problem in the energy norm. From the perspective of the sparse PCA, our results show that a large set of Mat\'{e}rn covariance functions can be approximated by a rank-$n$ operator with a localized basis and with the optimal accuracy.


title: Neural Vector Spaces for Unsupervised Information Retrieval

abstract: We propose the Neural Vector Space Model (NVSM), a method that learns representations of documents in an unsupervised manner for news article retrieval. In the NVSM paradigm, we learn low-dimensional representations of words and documents from scratch using gradient descent and rank documents according to their similarity with query representations that are composed from word representations. We show that NVSM performs better at document ranking than existing latent semantic vector space methods. The addition of NVSM to a mixture of lexical language models and a state-of-the-art baseline vector space model yields a statistically significant increase in retrieval effectiveness. Consequently, NVSM adds a complementary relevance signal. Next to semantic matching, we find that NVSM performs well in cases where lexical matching is needed. NVSM learns a notion of term specificity directly from the document collection without feature engineering. We also show that NVSM learns regularities related to Luhn significance. Finally, we give advice on how to deploy NVSM in situations where model selection (e.g., cross-validation) is infeasible. We find that an unsupervised ensemble of multiple models trained with different hyperparameter values performs better than a single cross-validated model. Therefore, NVSM can safely be used for ranking documents without supervised relevance judgments.


title: Ellipsoidal Prediction Regions for Multivariate Uncertainty Characterization

abstract: While substantial advances are observed in probabilistic forecasting for power system operation and electricity market applications, most approaches are still developed in a univariate framework. This prevents from informing about the interdependence structure among locations, lead times and variables of interest. Such dependencies are key in a large share of operational problems involving renewable power generation, load and electricity prices for instance. The few methods that account for dependencies translate to sampling scenarios based on given marginals and dependence structures. However, for classes of decision-making problems based on robust, interval chance-constrained optimization, necessary inputs take the form of polyhedra or ellipsoids. Consequently, we propose a systematic framework to readily generate and evaluate ellipsoidal prediction regions, with predefined probability and minimum volume. A skill score is proposed for quantitative assessment of the quality of prediction ellipsoids. A set of experiments is used to illustrate the discrimination ability of the proposed scoring rule for misspecification of ellipsoidal prediction regions. Application results based on three datasets with wind, PV power and electricity prices, allow us to assess the skill of the resulting ellipsoidal prediction regions, in terms of calibration, sharpness and overall skill.


title: Dynamics of interacting fermions under spin-orbit coupling in an optical lattice clock

abstract: Quantum statistics and symmetrization dictate that identical fermions do not interact via s-wave collisions. However, in the presence of spin-orbit coupling (SOC), fermions prepared in identical internal states with distinct momenta become distinguishable. The resulting strongly interacting system can exhibit exotic topological and pairing behaviors, many of which are yet to be observed in condensed matter systems. Ultracold atomic gases offer a promising pathway for simulating these rich phenomena. Two recent experiments reported the observation of single atom SOC in optical lattice clocks (OLCs) based on alkaline-earth atoms. In these works encoding the effective spin degree of freedom in the long-lived electronic clock states significantly reduced the detrimental effects of spontaneous emission and heating that have thus far hindered the study of interacting SOC with alkali atoms. Beyond first studies of interacting SOC with alkali atoms in a bulk gas and with two particles in a lattice, here we enter a new regime of many-body interacting SOC in an OLC. Using clock spectroscopy, we observe the precession of the collective magnetization and the emergence of spin locking effects arising from an interplay between p-wave and SOC-induced exchange interactions. The many-body dynamics are well captured by a collective XXZ spin model, which describes a broad class of condensed matter systems ranging from superconductors to quantum magnets. Furthermore, our work will aid in the design of next-generation OLCs by offering a route for avoiding the observed large density shifts caused by SOC-induced exchange interactions.


title: Jackknife multiplier bootstrap: finite sample approximations to the $U$-process supremum with applications

abstract: This paper is concerned with finite sample approximations to the supremum of a non-degenerate $U$-process of a general order indexed by a function class. We are primarily interested in situations where the function class as well as the underlying distribution change with the sample size, and the $U$-process itself is not weakly convergent as a process. Such situations arise in a variety of modern statistical problems. We first consider Gaussian approximations, namely, approximate the $U$-process supremum by the supremum of a Gaussian process, and derive coupling and Kolmogorov distance bounds. Such Gaussian approximations are, however, not often directly applicable in statistical problems since the covariance function of the approximating Gaussian process is unknown. This motivates us to study bootstrap-type approximations to the $U$-process supremum. We propose a novel jackknife multiplier bootstrap (JMB) tailored to the $U$-process, and derive coupling and Kolmogorov distance bounds for the proposed JMB method. All these results are non-asymptotic, and established under fairly general conditions on function classes and underlying distributions. Key technical tools in the proofs are new local maximal inequalities for $U$-processes, which may be useful in other problems. We also discuss applications of the general approximation results to testing for qualitative features of nonparametric functions based on generalized local $U$-processes.


title: Bubble Magnetometry of Nanoparticle Heterogeneity and Interaction

abstract: Bubbles have a rich history as transducers in particle-physics experiments. In a solid-state analogue, we use bubble domains in nanomagnetic films to measure magnetic nanoparticles. This technique can determine the magnetic orientation of a single nanoparticle in a fraction of a second, and generate a full hysteresis loop in a few seconds, which is much faster than any other reported technique. We achieve this unprecedented speed by tuning the nanomagnetic properties of the films, including the Dzyaloshinskii-Moriya interaction, in the first application of topological protection from the skyrmion state to a nanoparticle sensor. We demonstrate the technique on iron/nickel nanorods and iron oxide nanoparticles, which delineate a wide range of properties and applications. Bubble magnetometry enables the first measurement with high throughput for statistical analysis of the magnetic hysteresis of dispersed nanoparticles, and the first direct measurement of a transition from superparamagnetic behavior as single nanoparticles to collective behavior in nanoscale agglomerates. These results demonstrate a breakthrough capability for measuring the heterogeneity and interaction of magnetic nanoparticles.


title: A method to adjust the impedance of the transmission line in a Multi-Strip Multi-Gap Resistive Plate Counter

abstract: While in a triggered experiment the matching of the RPC transmission line impedance with the one of the front-end electronics is less critical, for a trigger-less data recording this becomes mandatory. As expected, impedance matching is not straightforward when other requirements in terms of time and position resolutions, efficiency and granularity, have to be fulfilled in the same time. A method and the very first results obtained with a RPC prototype built based on it, presented in this paper, show that the impedance matching, independent of its granularity, can be achieved using an innovative architecture of the RPC.


title: Generalized Fr\'echet Bounds for Cell Entries in Multidimensional Contingency Tables

abstract: We consider the lattice, $\mathcal{L}$, of all subsets of a multidimensional contingency table and establish the properties of monotonicity and supermodularity for the marginalization function, $n(\cdot)$, on $\mathcal{L}$. We derive from the supermodularity of $n(\cdot)$ some generalized Fr\'echet inequalities complementing and extending inequalities of Dobra and Fienberg. Further, we construct new monotonic and supermodular functions from $n(\cdot)$, and we remark on the connection between supermodularity and some correlation inequalities for probability distributions on lattices. We also apply an inequality of Ky Fan to derive a new approach to Fr\'echet inequalities for multidimensional contingency tables.


title: Recent Trends in Deep Learning Based Natural Language Processing

abstract: Deep learning methods employ multiple processing layers to learn hierarchical representations of data and have produced state-of-the-art results in many domains. Recently, a variety of model designs and methods have blossomed in the context of natural language processing (NLP). In this paper, we review significant deep learning related models and methods that have been employed for numerous NLP tasks and provide a walk-through of their evolution. We also summarize, compare and contrast the various models and put forward a detailed understanding of the past, present and future of deep learning in NLP.