title: Well-Founded Operators for Normal Hybrid MKNF Knowledge Bases

abstract: Hybrid MKNF knowledge bases have been considered one of the dominant approaches to combining open world ontology languages with closed world rule-based languages. Currently, the only known inference methods are based on the approach of guess-and-verify, while most modern SAT/ASP solvers are built under the DPLL architecture. The central impediment here is that it is not clear what constitutes a constraint propagator, a key component employed in any DPLL-based solver. In this paper, we address this problem by formulating the notion of unfounded sets for nondisjunctive hybrid MKNF knowledge bases, based on which we propose and study two new well-founded operators. We show that by employing a well-founded operator as a constraint propagator, a sound and complete DPLL search engine can be readily defined. We compare our approach with the operator based on the alternating fixpoint construction by Knorr et al [2011] and show that, when applied to arbitrary partial partitions, the new well-founded operators not only propagate more truth values but also circumvent the non-converging behavior of the latter. In addition, we study the possibility of simplifying a given hybrid MKNF knowledge base by employing a well-founded operator, and show that, out of the two operators proposed in this paper, the weaker one can be applied for this purpose and the stronger one cannot. These observations are useful in implementing a grounder for hybrid MKNF knowledge bases, which can be applied before the computation of MKNF models. The paper is under consideration for acceptance in TPLP.


title: Relaxation time and critical slowing down of a spin-torque oscillator

abstract: The relaxation phenomena of spin-torque oscillators consisting of nanostructured ferromagnets are interesting research targets in magnetism. A theoretical study on the relaxation time of a spin-torque oscillator from one self-oscillation state to another is investigated. By solving the Landau-Lifshitz-Gilbert equation both analytically and numerically, it is shown that the oscillator relaxes to the self-oscillation state exponentially within a few nanoseconds, except when magnetization is close to a critical point. The relaxation rate, which is an inverse of relaxation time, is proportional to the current. On the other hand, a critical slowing down appears near the critical point, where relaxation is inversely proportional to time, and the relaxation time becomes on the order of hundreds of nanoseconds. These conclusions are primarily obtained for a spin-torque oscillator consisting of a perpendicularly magnetized free layer and an in-plane magnetized pinned layer, and are further developed for application to arbitrary types of spin-torque oscillators.


title: Long-Term Memory Networks for Question Answering

abstract: Question answering is an important and difficult task in the natural language processing domain, because many basic natural language processing tasks can be cast into a question answering task. Several deep neural network architectures have been developed recently, which employ memory and inference components to memorize and reason over text information, and generate answers to questions. However, a major drawback of many such models is that they are capable of only generating single-word answers. In addition, they require large amount of training data to generate accurate answers. In this paper, we introduce the Long-Term Memory Network (LTMN), which incorporates both an external memory module and a Long Short-Term Memory (LSTM) module to comprehend the input data and generate multi-word answers. The LTMN model can be trained end-to-end using back-propagation and requires minimal supervision. We test our model on two synthetic data sets (based on Facebook's bAbI data set) and the real-world Stanford question answering data set, and show that it can achieve state-of-the-art performance.


title: Disorder and the neural representation of complex odors: smelling in the real world

abstract: Animals smelling in the real world use a small number of receptors to sense a vast number of natural molecular mixtures, and proceed to learn arbitrary associations between odors and valences. Here, we propose a new interpretation of how the architecture of olfactory circuits is adapted to meet these immense complementary challenges. First, the diffuse binding of receptors to many molecules compresses a vast odor space into a tiny receptor space, while preserving similarity. Next, lateral interactions "densify" and decorrelate the response, enhancing robustness to noise. Finally, disordered projections from the periphery to the central brain reconfigure the densely packed information into a format suitable for flexible learning of associations and valences. We test our theory empirically using data from Drosophila. Our theory suggests that the neural processing of olfactory information differs from the other senses in its fundamental use of disorder.


title: On Thompson's conjecture for finite simple exceptional groups of Lie type

abstract: Let $G$ be a finite group, $N(G)$ be the set of conjugacy classes of the group $G$. In the present paper it is proved $G\simeq L$ if $N(G)=N(L)$, where $G$ is a finite group with trivial center and $L$ is a finite simple group of exceptional Lie type or Tits group.


title: Controllability and Stabilizability Analysis of Signed Consensus Networks

abstract: Signed networks have been a topic of recent interest in the network control community as they allow studying antagonistic interactions in multi-agent systems. Although dynamical characteristics of signed networks have been well-studied, notions such as controllability and stabilizability for signed networks for protocols such as consensus are missing in the literature. Classically, graph automorphisms with respect to the input nodes have been used to characterize uncontrollability of consensus networks. In this paper, we show that in addition to the graph symmetry, the topological property of structural balance facilitates the derivation of analogous sufficient conditions for uncontrollability for signed networks. In particular, we provide an analysis which shows that a gauge transformation induced by structural balance allows symmetry arguments to hold for signed consensus networks. Lastly, we use fractional automorphisms to extend our observations to output controllability and stabilizability of signed networks.


title: Distributed Nash Equilibrium Seeking under Partial-Decision Information via the Alternating Direction Method of Multipliers

abstract: In this paper, we consider the problem of finding a Nash equilibrium in a multi-player game over generally connected networks. This model differs from a conventional setting in that players have partial information on the actions of their opponents and the communication graph is not necessarily the same as the players' cost dependency graph. We develop a relatively fast algorithm within the framework of inexact-ADMM, based on local information exchange between the players. We prove its convergence to Nash equilibrium for fixed step-sizes and analyze its convergence rate. Numerical simulations illustrate its benefits when compared to a consensus-based gradient type algorithm with diminishing step-sizes.


title: Unified framework to determine Gaussian states in continuous variable systems

abstract: Gaussian states are the backbone of quantum information protocols with continuous variable systems, whose power relies fundamentally on the entanglement between the different modes. In the case of global pure states, knowledge of the reduced states in a given bipartition of a multipartite quantum system bears information on the entanglement in such bipartition. For Gaussian states, the reduced states are also Gaussian, so there determination requires essentially the experimental determination of their covariance matrix. Here, we develop strategies to determine the covariance matrix of an arbitrary n-mode bosonic Gaussian state through measurement of the total phase acquired when appropriate metaplectic evolutions, associated with quadratic Hamiltonians, are applied. Simply one-mode metaplectic evolutions, such rotations, squeezing and shear transformations, in addition to a single two-mode rotation, allows to determine all the covariance matrix elements of a n-mode bosonic system. All the single-mode metaplectic evolutions are applied conditionally to a state in which an ancilla qubit is entangled with the n-mode system. The ancillary system provides, after measurement, the value of the total phase of each evolution. The proposed method is experimentally friendly to be implemented in the most currently used continuous variable systems.


title: Signed graphs and the freeness of the Weyl subarrangements of type $B_{\ell}$

abstract: A Weyl arrangement is the hyperplane arrangement defined by a root system. Arnold and Saito proved that every Weyl arrangement is free. The Weyl subarrangements of type $A_{\ell}$ are represented by simple graphs. Stanley gave a characterization of freeness for this type of arrangements in terms of thier graph. In addition, The Weyl subarrangements of type $B_{\ell}$ can be represented by signed graphs. A characterization of freeness for them is not known. However, characterizations of freeness for a few restricted classes are known. For instance, Edelman and Reiner characterized the freeness of the arrangements between type $ A_{\ell-1} $ and type $ B_{\ell} $. In this paper, we give a characterization of the freeness and supersolvability of the Weyl subarrangements of type $B_{\ell}$ under certain assumption.


title: An IMEX Finite Element Method for a linearized Cahn-Hilliard-Cook equation driven by the space derivative of a space-time white noise

abstract: We consider a model initial- and Dirichlet boundary- value problem for a linearized Cahn-Hilliard-Cook equation, in one space dimension, forced by the space derivative of a space-time white noise. First, we introduce a canvas problem the solution to which is a regular approximation of the mild solution to the problem and depends on a finite number of random variables. Then, fully-discrete approximations of the solution to the canvas problem are constructed using, for discretization in space, a Galerkin finite element method based on $H^2$ piecewise polynomials, and, for time-stepping, an implicit/explicit method. Finally, we derive a strong a priori estimate of the error approximating the mild solution to the problem by the canvas problem solution, and of the numerical approximation error of the solution to the canvas problem.


title: Load Balancing in the Non-Degenerate Slowdown Regime

abstract: We analyse Join-the-Shortest-Queue in a contemporary scaling regime known as the Non-Degenerate Slowdown regime. Join-the-Shortest-Queue (JSQ) is a classical load balancing policy for queueing systems with multiple parallel servers. Parallel server queueing systems are regularly analysed and dimensioned by diffusion approximations achieved in the Halfin-Whitt scaling regime. However, when jobs must be dispatched to a server upon arrival, we advocate the Non-Degenerate Slowdown regime (NDS) to compare different load-balancing rules. In this paper we identify novel diffusion approximation and timescale separation that provides insights into the performance of JSQ. We calculate the price of irrevocably dispatching jobs to servers and prove this to within 15% (in the NDS regime) of the rules that may manoeuvre jobs between servers. We also compare ours results for the JSQ policy with the NDS approximations of many modern load balancing policies such as Idle-Queue-First and Power-of-$d$-choices policies which act as low information proxies for the JSQ policy. Our analysis leads us to construct new rules that have identical performance to JSQ but require less communication overhead than power-of-2-choices.


title: Phonon Spectra in the Parent Superconducting Iron-tuned Telluride Fe$_{1+x}$Te from Inelastic Neutron Scattering and Ab Initio Calculations

abstract: We report inelastic neutron scattering measurements of phonon spectra in the parent superconductor iron-tuned chalcogenide Fe$_{1+x}$Te, for two different x contents (x $\leq$ 0.11), using neutron time-of-flight technique. Thermal neutron spectroscopy allowed to collect the low-temperature Stokes spectra over an extended Q-range, at 2, 40 and 120K - hence covering both the magnetic monoclinic and the paramagnetic tetragonal phases. Whereas cold-neutrons allowed to measure high-resolution anti-Stokes spectra at 140, 220 and 300K, thus covering the tetragonal phase. Our results evidence a spin-phonon coupling behaviour towards the observed noticeable temperature-dependent change of the Stokes spectra across the transition temperatures. On the other hand, the anti-Stokes spectra reveal a pronounced hardening of the low-energy, acoustic region, of the phonon spectrum, upon heating, indicating a strong anharmonicity and a subtle dependence of phonons on structural evolution within the tetragonal phase. Experimental results are accompanied by ab initio calculations of phonon spectra of the tetragonal stoichiometric phase for a comparison with the high-resolution anti-Stokes spectra. Calculations included different density functional methods. Spin polarization and van der Waals interaction, were either considered or neglected, individually or concomitantly, in order to study their respective effect on lattice dynamics description. Our results suggest that including van der Waals interaction has only a slight effect on phonon dynamics, however, phonon spectra are better described when spin polarization is included, in a cooperative way with van der Waals interactions.


title: Algorithms for zero-dimensional ideals using linear recurrent sequences

abstract: Inspired by Faug\`ere and Mou's sparse FGLM algorithm, we show how using linear recurrent multi-dimensional sequences can allow one to perform operations such as the primary decomposition of an ideal, by computing the annihilator of one or several such sequences.


title: Model Based Diagnosis of Multiple Observations with Implicit Hitting Sets

abstract: Model based diagnosis finds a growing range of practical applications, and significant performance-wise improvements have been achieved in recent years. Some of these improvements result from formulating the problem with maximum satisfiability (MaxSAT). Whereas recent work focuses on analyzing failing observations separately, it is also the case that in practical settings there may exist many failing observations. This paper first investigates the drawbacks of analyzing failing observations separately. It then shows that existing solutions do not scale for large systems. Finally, the paper proposes a novel approach for diagnosing systems with many failing observations. The proposed approach is based on implicit hitting sets and so is tightly related with the original seminal work on model based diagnosis. The experimental results demonstrate not only the importance of analyzing multiple observations simultaneously, but also the significance of the implicit hitting set approach.


title: Probing the cold and warm molecular gas in the Whirlpool Galaxy: Herschel SPIRE-FTS Observations of the central region of M51 (NGC 5194)

abstract: We present Herschel SPIRE-FTS intermediate-sampled mapping observations of the central ~8 kpc (~150") of M51, with a spatial resolution of 40". We detect 4 12CO transitions (J=4-3 to J=7-6) and the CI 3P2-3P1 and 3P1-3P0 transitions. We supplement these observations with ground based observations of 12CO J=1-0 to J=3-2 and perform a two-component non-LTE analysis. We find that the molecular gas in the nucleus and centre regions has a cool component (T_kin~10-20 K) with a moderate but poorly constrained density (n(H2)~10^3-10^6 cm^-3), as well as significant molecular gas in a warmer (T_kin~300-3000 K), lower density (n(H2)~10^1.6-10^2.5 cm^-3) component. We compare our CO line ratios and calculated densities along with ratios of CO to total infrared luminosity to a grid of photon dominated region (PDR) models and find that the cold molecular gas likely resides in PDRs with a field strength of G_o~10^2. The warm component likely requires an additional source of mechanical heating, from supernovae and stellar winds or possibly shocks produced in the strong spiral density wave. When compared to similar two-component models of other star-forming galaxies published as part of the Very Nearby Galaxies Survey (Arp 220, M82 and NGC 4038/39), M51 has the lowest density for the warm component, while having a warm gas mass fraction that is comparable to those of Arp 220 and M82, and significantly higher than that of NGC 4038/39.


title: A Tissue Engineered Model of Aging: Interdependence and Cooperative Effects in Failing Tissues

abstract: Aging remains a fundamental open problem in modern biology. Although there exist a number of theories on aging on the cellular scale, nearly nothing is known about how microscopic failures cascade to macroscopic failures of tissues, organs and ultimately the organism. The goal of this work is to bridge microscopic cell failure to macroscopic manifestations of aging. We use tissue engineered constructs to control the cellular-level damage and cell-cell distance in individual tissues to establish the role of complex interdependence and interactions between cells in aging tissues. We found that while microscopic mechanisms drive aging, the interdependency between cells plays a major role in tissue death, providing evidence on how cellular aging is connected to its higher systemic consequences.


title: Higher-order congruence relations on affine moment graphs: The subgeneric case

abstract: We study the structure algebra $\mathcal{Z}$ of the stable moment graph for the case of the affine root system $A_{1}$. The structure algebra $\mathcal{Z}$ is an algebra over a symmetric algebra and in particular, it is a module over a symmetric algebra. We study this module structure on $\mathcal{Z}$ and we construct a basis. By "setting $c$ equal to zero" in $\mathcal{Z}$, we obtain the module $\mathcal{Z}_{c=0}$. This module can be described in terms of the finite root system $A_{1}$ and we show that it is determined by a set of certain divisibility relations. These relations can be regarded as a generalization of ordinary moment graph relations that define sections of sheaves on moment graphs, and because of this we call them higher-order congruence relations.


title: A Generalised Seizure Prediction with Convolutional Neural Networks for Intracranial and Scalp Electroencephalogram Data Analysis

abstract: Seizure prediction has attracted a growing attention as one of the most challenging predictive data analysis efforts in order to improve the life of patients living with drug-resistant epilepsy and tonic seizures. Many outstanding works have been reporting great results in providing a sensible indirect (warning systems) or direct (interactive neural-stimulation) control over refractory seizures, some of which achieved high performance. However, many works put heavily handcraft feature extraction and/or carefully tailored feature engineering to each patient to achieve very high sensitivity and low false prediction rate for a particular dataset. This limits the benefit of their approaches if a different dataset is used. In this paper we apply Convolutional Neural Networks (CNNs) on different intracranial and scalp electroencephalogram (EEG) datasets and proposed a generalized retrospective and patient-specific seizure prediction method. We use Short-Time Fourier Transform (STFT) on 30-second EEG windows with 50% overlapping to extract information in both frequency and time domains. A standardization step is then applied on STFT components across the whole frequency range to prevent high frequencies features being influenced by those at lower frequencies. A convolutional neural network model is used for both feature extraction and classification to separate preictal segments from interictal ones. The proposed approach achieves sensitivity of 81.4%, 81.2%, 82.3% and false prediction rate (FPR) of 0.06/h, 0.16/h, 0.22/h on Freiburg Hospital intracranial EEG (iEEG) dataset, Children's Hospital of Boston-MIT scalp EEG (sEEG) dataset, and Kaggle American Epilepsy Society Seizure Prediction Challenge's dataset, respectively. Our prediction method is also statistically better than an unspecific random predictor for most of patients in all three datasets.


title: A Two-band Model for p-wave Superconductivity

abstract: In this paper we study the effects of hybridization in the superconducting properties of a two-band system. We consider the cases that these bands are formed by electronic orbitals with angular momentum, such that, the hybridization $V(\mathbf{k})$ among them can be symmetric or antisymmetric under inversion symmetry. We take into account only intra-band attractive interactions in the two bands and investigate the appearance of an induced inter-band pairing gap. We show that (inter-band) superconducting orderings are induced in the total absence of attractive interaction between the two bands, which turns out to be completely dependent on the hybridization between them. For the case of antisymmetric hybridization we show that the induced inter-band superconductivity has a p-wave symmetry.


title: Local Large deviation: A McMillian Theorem for Coloured Random Graph Processes

abstract: For a finite typed graph on $n$ nodes and with type law $\mu,$ we define the so-called spectral potential $\rho_{\lambda}(\,\cdot,\,\mu),$ of the graph.From the $\rho_{\lambda}(\,\cdot,\,\mu)$ we obtain Kullback action or the deviation function, $\mathcal{H}_{\lambda}(\pi\,\|\,\nu),$ with respect to an empirical pair measure, $\pi,$ as the Legendre dual. For the finite typed random graph conditioned to have an empirical link measure $\pi$ and empirical type measure $\mu$, we prove a Local large deviation principle (LLDP), with rate function $\mathcal{H}_{\lambda}(\pi\,\|\,\nu)$ and speed $n.$ We deduce from this LLDP, a full conditional large deviation principle and a weak variant of the classical McMillian Theorem for the typed random graphs. Given the typical empirical link measure, $\lambda\mu\otimes\mu,$ the number of typed random graphs is approximately equal $e^{n\|\lambda\mu\otimes\mu\|H\big(\lambda\mu\otimes\mu/\|\lambda\mu\otimes\mu\|\big)}.$ Note that we do not require any topological restrictions on the space of finite graphs for these LLDPs.


title: Impact of an introductory lab on students' understanding of measurement uncertainty

abstract: Physics lab courses are an essential part of the physics undergraduate curriculum. Learning goals for these classes often include the ability to interpret measurements and uncertainties. The Physics Measurement Questionnaire (PMQ) is an established open-response survey that probes students' understanding of measurement uncertainty along three dimensions: data collection, data analysis, and data comparison. It classifies students' reasoning into point-like and set-like paradigms, with the set-like paradigm more aligned with expert reasoning. In the context of a course transformation effort at the University of Colorado Boulder, we examine over 500 student responses to the PMQ both before and after instruction in the pre-transformed course. We describe changes in students' overall reasoning, measured by aggregating four probes of the PMQ. In particular, we observe large shifts towards set-like reasoning by the end of the course.


title: Student reasoning about measurement uncertainty in an introductory lab course

abstract: Proficiency with calculating, reporting, and understanding measurement uncertainty is a nationally recognized learning outcome for undergraduate physics lab courses. The Physics Measurement Questionnaire (PMQ) is a research-based assessment tool that measures such understanding. The PMQ was designed to characterize student reasoning into point or set paradigms, where the set paradigm is more aligned with expert reasoning. We analyzed over 500 student open-ended responses collected at the beginning and the end of a traditional introductory lab course at the University of Colorado Boulder. We discuss changes in students' understanding over a semester by analyzing pre-post shifts in student responses regarding data collection, data analysis, and data comparison.


title: An stratification of $B^4(2,K_C)$ over a general curve

abstract: For a general curve C of genus $g\geq 10$, we show that the Brill- Noether locus $B^4(2,K_C)$ contains irreducible sub-varieties $B_3\supset B_4 \supset \cdots \supset B_n$, where $B_n$ is of dimension $3g-10-n$ and $B_3$ is an irreducible component of the expected dimension $3g-13$.


title: No Tau? No Astronomy!

abstract: Since 2013 IceCube cascade showers sudden overabundance have shown a fast flavor change above 30-60 TeV up to PeV energy. This flavor change from dominant muon tracks at TeVs to shower events at higher energies, has been indebted to a new injection of a neutrino astronomy. However the recent published 54 neutrino HESE, high energy starting events, as well as the 38 external muon tracks made by trough going muon formed around the IceCube, none of them are pointing to any expected X-gamma or radio sources: no one in connection to GRB, no toward active BL Lac, neither to AGN source in Fermi catalog. No clear correlation with nearby mass distribution (Local Group), nor to galactic plane. Moreover there have not been any record (among a dozen of 200 TeV energetic events) of the expected double bang due to the tau neutrino birth and decay. An amazing and surprising unfair distribution in flavor versus expected democratic one. Finally there is not a complete consistence of the internal HESE event spectra and the external crossing muon track ones. Moreover the apparent sudden astrophysical neutrino flux rise at 60 TeV might be probably also suddenly cut at a few PeV in order to hide the (unobserved , yet) Glashow resonance peak at 6.3 PeV. A more mondane prompt charmed atmospheric neutrino component may explain most of the IceCube puzzles. If this near future, 2017-2018, it does not shine tau neutrino signals somewhere (by tau airshowers in AUGER, TA, ASHRA or double bang in IceCube) there are a list of consequences to face. These missing correlations and in particular the tau signature absence force us to claim : No Tau? No neutrino Astronomy.


title: Phase transitions in the $q$-coloring of random hypergraphs

abstract: We study in this paper the structure of solutions in the random hypergraph coloring problem and the phase transitions they undergo when the density of constraints is varied. Hypergraph coloring is a constraint satisfaction problem where each constraint includes $K$ variables that must be assigned one out of $q$ colors in such a way that there are no monochromatic constraints, i.e. there are at least two distinct colors in the set of variables belonging to every constraint. This problem generalizes naturally coloring of random graphs ($K=2$) and bicoloring of random hypergraphs ($q=2$), both of which were extensively studied in past works. The study of random hypergraph coloring gives us access to a case where both the size $q$ of the domain of the variables and the arity $K$ of the constraints can be varied at will. Our work provides explicit values and predictions for a number of phase transitions that were discovered in other constraint satisfaction problems but never evaluated before in hypergraph coloring. Among other cases we revisit the hypergraph bicoloring problem ($q=2$) where we find that for $K=3$ and $K=4$ the colorability threshold is not given by the one-step-replica-symmetry-breaking analysis as the latter is unstable towards more levels of replica symmetry breaking. We also unveil and discuss the coexistence of two different 1RSB solutions in the case of $q=2$, $K \ge 4$. Finally we present asymptotic expansions for the density of constraints at which various phase transitions occur, in the limit where $q$ and/or $K$ diverge.

categories: cond-mat.dis-nn cond-mat.stat-mech cs.DM math.PR


title: Dynamical pruning of rooted trees with applications to 1D ballistic annihilation

abstract: We introduce generalized dynamical pruning on rooted binary trees with edge lengths. The pruning removes parts of a tree $T$, starting from the leaves, according to a pruning function defined on subtrees within $T$. The generalized pruning encompasses a number of discrete and continuous pruning operations, including the tree erasure and Horton pruning. The main result is invariance of a finite critical binary Galton-Watson tree with exponential edge lengths with respect to the generalized dynamical pruning for an arbitrary admissible pruning function. The second part of the paper examines the continuum 1-D ballistic annihilation model $A+A \rightarrow \emptyset$ for a constant particle density and initial velocity that alternates between the values of $\pm$1. The model evolution is equivalent to a generalized dynamical pruning of the shock tree that represents dynamics of sinks (points of particle annihilation), with the pruning function equal to the total tree length. The shock tree is isometric to the level set tree of the model potential (integral of velocity). This equivalence allows us to construct a complete probabilistic description of the annihilation dynamics for the initial velocity that alternates between the values of $\pm$1 at the epochs of a stationary Poisson process. Finally, we discuss several real tree representations of the ballistic annihilation model, closely connected to the shock wave tree.


title: Coiling of an elastic beam inside a disk: A model for spider-capture silk

abstract: Motivated by recent experimental observations of capillary-induced spooling of fibers inside droplets both in spider capture silk and in synthetic systems, we investigate the behavior of a fiber packed in a drop. Using a simplified 2D model, we provide analytical predictions for the buckling threshold and the deep post-buckling asymptotic behavior. The threshold for spooling is found to be in particularly good agreement with experimental results. We further solve the Elastica equations for a fiber confined in a soft potential, and track the equilibrium paths using numerical continuation techniques. A wealth of different paths corresponding to different symmetries is uncovered, and their stability is finally discussed.


title: $LlogL$-integrability of the velocity gradient for Stokes system with drifts in $L_\infty (BMO^{-1})$

abstract: For any weak solution of the Stokes system with drifts in $L_\infty (BMO^{-1})$, we prove a reverse H\"older inequality and $LlogL$-higher integrability of the velocity gradients.


title: Observing future Solar Flare Neutrino in Hyper-KamioKande in Japan, Korea and in IceCube

abstract: The largest Solar flare have been recorded in gamma flash and hard spectra up to tens GeV energy. The present building and upgrade of Hyper-Kamiokande (HK) and IceCube (as well as Deep Core, PINGU) neutrino detectors do offer a novel way to largest trace of solar Flare: their sudden anti-neutrino (or neutrino) flare made by proton scattering and pion decays via Delta resonance production. These signals might be observable at largest flare by HK via soft spectra up to tens- hundred MeV energy and by IceCube-PINGU at higher, GeV up to hundred GeV, energies. We show the expected rate of signals for the most powerful solar flare occurred in recent decades.


title: Beam Energy Dependence of Jet-Quenching Effects in Au+Au Collisions at $\sqrt{s_{_{ \mathrm{NN}}}}$ = 7.7, 11.5, 14.5, 19.6, 27, 39, and 62.4 GeV

abstract: We report measurements of the nuclear modification factor, $R_{ \mathrm{CP}}$, for charged hadrons as well as identified $\pi^{+(-)}$, $K^{+(-)}$, and $p(\overline{p})$ for Au+Au collision energies of $\sqrt{s_{_{ \mathrm{NN}}}}$ = 7.7, 11.5, 14.5, 19.6, 27, 39, and 62.4 GeV. We observe a clear high-$p_{\mathrm{T}}$ net suppression in central collisions at 62.4 GeV for charged hadrons which evolves smoothly to a large net enhancement at lower energies. This trend is driven by the evolution of the pion spectra, but is also very similar for the kaon spectra. While the magnitude of the proton $R_{ \mathrm{CP}}$ at high $p_{\mathrm{T}}$ does depend on collision energy, neither the proton nor the anti-proton $R_{ \mathrm{CP}}$ at high $p_{\mathrm{T}}$ exhibit net suppression at any energy. A study of how the binary collision scaled high-$p_{\mathrm{T}}$ yield evolves with centrality reveals a non-monotonic shape that is consistent with the idea that jet-quenching is increasing faster than the combined phenomena that lead to enhancement.


title: Cooperative Kernels: GPU Multitasking for Blocking Algorithms (Extended Version)

abstract: There is growing interest in accelerating irregular data-parallel algorithms on GPUs. These algorithms are typically blocking, so they require fair scheduling. But GPU programming models (e.g.\ OpenCL) do not mandate fair scheduling, and GPU schedulers are unfair in practice. Current approaches avoid this issue by exploiting scheduling quirks of today's GPUs in a manner that does not allow the GPU to be shared with other workloads (such as graphics rendering tasks). We propose cooperative kernels, an extension to the traditional GPU programming model geared towards writing blocking algorithms. Workgroups of a cooperative kernel are fairly scheduled, and multitasking is supported via a small set of language extensions through which the kernel and scheduler cooperate. We describe a prototype implementation of a cooperative kernel framework implemented in OpenCL 2.0 and evaluate our approach by porting a set of blocking GPU applications to cooperative kernels and examining their performance under multitasking. Our prototype exploits no vendor-specific hardware, driver or compiler support, thus our results provide a lower-bound on the efficiency with which cooperative kernels can be implemented in practice.


title: Eigenvalues of the Thurston operator

abstract: Let $f:\hat{\mathbb C}\to \hat{\mathbb C}$ be a postcritically finite rational map. Let $\mathcal Q(\hat{\mathbb C})$ be the space of meromorphic quadratic differentials on $ \hat{\mathbb C}$ with simple poles. We study the set of eigenvalues of the pushforward operator $f_*:\mathcal Q(\hat{\mathbb C})\to \mathcal Q(\hat{\mathbb C})$. In particular, we show that when $f:\mathbb C \to \mathbb C$ is a unicritical polynomial of degree $D$ with periodic critical point, the eigenvalues of $f_*:\mathcal Q(\hat{\mathbb C})\to \mathcal Q(\hat{\mathbb C})$ are contained in the annulus $\bigl\{\frac{1}{4D}<|\lambda|<1\bigr\}$ and belong to $\frac{1}{D} \mathbb U$ where $\mathbb U$ is the group of algebraic units.


title: Injection-locked phonon lasing driven by an individual plasmonic metamolecule

abstract: Nanophotonic devices take advantage of geometry-dependent optical properties to confine and enhance the interaction of light with matter on small scales. By carefully patterning nanoscale geometries, coupling of responses across distinct physical domains can be achieved, such as photon and phonon modes in optomechanical crystals. Extreme optical confinement of plasmonics can localize strong cross-domain interactions into nanometer-scale volumes, smaller than many nanosystem's thermal and mechanical eigenmodes, opening new regimes for engineering reconfigurable, dynamic nanosystems. Here, we present a nanosystem based on an electrically actuated, mechanically coupled gap-plasmon metamolecule resonator. Integrated electrical actuation reveals distinct dispersive and reactive optomechanical coupling regimes of the metamolecules, and enables both broad spectral tuning of individual plasmonic resonances and selective optical transduction of nanomechanical modes. Strong thermomechanical backaction from a single point-like metamolecule yields optically driven self-oscillation, and we show that this metamolecule phonon laser can be injection-locked to enhance weak mechanical stimulii. This multi-domain nanosystem advances nanomechanical sensing and optical modulation, opens new avenues for studying nanoscale nonlinear optomechanics, and may serve as a building block for future reconfigurable "smart" metamaterials.


title: On the Compactness, Efficiency, and Representation of 3D Convolutional Networks: Brain Parcellation as a Pretext Task

abstract: Deep convolutional neural networks are powerful tools for learning visual representations from images. However, designing efficient deep architectures to analyse volumetric medical images remains challenging. This work investigates efficient and flexible elements of modern convolutional networks such as dilated convolution and residual connection. With these essential building blocks, we propose a high-resolution, compact convolutional network for volumetric image segmentation. To illustrate its efficiency of learning 3D representation from large-scale image data, the proposed network is validated with the challenging task of parcellating 155 neuroanatomical structures from brain MR images. Our experiments show that the proposed network architecture compares favourably with state-of-the-art volumetric segmentation networks while being an order of magnitude more compact. We consider the brain parcellation task as a pretext task for volumetric image segmentation; our trained network potentially provides a good starting point for transfer learning. Additionally, we show the feasibility of voxel-level uncertainty estimation using a sampling approximation through dropout.


title: The Way Forward - Closing Remarks at Quark Matter 2017

abstract: This contribution is a written version of my closing talk presented at the Quark Matter 2017 conference. Neither the talk nor this contribution to the conference proceedings is intended as a comprehensive summary. Rather, a brief discussion is presented of emerging themes and challenges in the field of relativistic heavy ion physics.


title: New non-linear modified massless Klein--Gordon equation

abstract: The massless Klein--Gordon equation on arbitrary curved backgrounds allows for solutions which develop "tails" inside the light cone and, therefore, do not strictly follow null geodesics as discovered by DeWitt and Brehme almost sixty years ago. A modification of the massless Klein--Gordon equation is presented which always exhibits null geodesic propagation of waves on arbitrary curved spacetimes. This new equation is derived from a Lagrangian which exhibits current--current interaction. Its non--linearity is due to a self--coupling term which is related to the quantum mechanical Bohm potential.


title: Ground state properties of 3d metals from self-consistent GW approach

abstract: Self consistent GW approach (scGW) has been applied to calculate the ground state properties (equilibrium Wigner-Seitz radius $S_{WZ}$ and bulk modulus $B$) of 3d transition metals Sc, Ti, V, Fe, Co, Ni, and Cu. The approach systematically underestimates $S_{WZ}$ with average relative deviation from the experimental data about 1% and it overestimates the calculated bulk modulus with relative error about 25%. It is shown that scGW is superior in accuracy as compared to the local density approximation (LDA) but it is less accurate than the generalized gradient approach (GGA) for the materials studied. If compared to the random phase approximation (RPA), scGW is slightly less accurate, but its error for the 3d metals looks more systematic. The systematic nature of the deviation from the experimental data suggests that the next order of the perturbation theory should allow one to reduce the error.


title: Capacity of Wireless Distributed Storage Systems with Broadcast Repair

abstract: In wireless distributed storage systems, storage nodes are connected by wireless channels, which are broadcast in nature. This paper exploits this unique feature to design an efficient repair mechanism, called broadcast repair, for wireless distributed storage systems in the presence of multiple-node failures. Due to the broadcast nature of wireless transmission, we advocate a new measure on repair performance called repair-transmission bandwidth. In contrast to repair bandwidth, which measures the average number of packets downloaded by a newcomer to replace a failed node, repair-transmission bandwidth measures the average number of packets transmitted by helper nodes per failed node. A fundamental study on the storage capacity of wireless distributed storage systems with broadcast repair is conducted by modeling the storage system as a multicast network and analyzing the minimum cut of the corresponding information flow graph. The fundamental tradeoff between storage efficiency and repair-transmission bandwidth is also obtained for functional repair. The performance of broadcast repair is compared both analytically and numerically with that of cooperative repair, the basic repair method for wired distributed storage systems with multiple-node failures. While cooperative repair is based on the idea of allowing newcomers to exchange packets, broadcast repair is based on the idea of allowing a helper to broadcast packets to all newcomers simultaneously. We show that broadcast repair outperforms cooperative repair, offering a better tradeoff between storage efficiency and repair-transmission bandwidth.


title: Comparing Jupiter interior structure models to Juno gravity measurements and the role of a dilute core

abstract: The Juno spacecraft has measured Jupiter's low-order, even gravitational moments, $J_2$--$J_8$, to an unprecedented precision, providing important constraints on the density profile and core mass of the planet. Here we report on a selection of interior models based on ab initio computer simulations of hydrogen-helium mixtures. We demonstrate that a dilute core, expanded to a significant fraction of the planet's radius, is helpful in reconciling the calculated $J_n$ with Juno's observations. Although model predictions are strongly affected by the chosen equation of state, the prediction of an enrichment of $Z$ in the deep, metallic envelope over that in the shallow, molecular envelope holds. We estimate Jupiter's core to contain an 7--25 Earth mass of heavy elements. We discuss the current difficulties in reconciling measured $J_n$ with the equations of state, and with theory for formation and evolution of the planet.


title: Utilizing Microcavities to Suppress Third-order Cascades in Fifth-order Raman Spectra

abstract: Nonlinear optical signals in the condensed phase are often accompanied by sequences of lower-order processes, known as cascades, which share the same phase matching and power dependence on the incoming fields and are thus hard to distinguish. The suppression of cascading in order to reveal the desired nonlinear signal has been a major challenge in multidimensional Raman spectroscopy, i.e., the $\chi^{(5)}$ signal being masked by cascading signals given by a product of two $\chi^{(3)}$ processes. Since cascading originates from the exchange of a virtual photon between molecules, it can be manipulated by performing the experiment in an optical microcavity. Using a quantum electrodynamical (QED) treatment we demonstrate that the $\chi^{(3)}$ cascading contributions can be greatly suppressed. By optimizing the cavity size and the incoming pulse directions, we show that up to $\sim$99.5\% suppression of the cascading signal is possible.


title: Constraining the HI-Halo Mass Relation From Galaxy Clustering

abstract: We study the dependence of galaxy clustering on atomic gas mass using a sample of $\sim$16,000 galaxies with redshift in the range of $0.0025<z<0.05$ and HI mass of $M_{\rm HI}>10^8M_{\odot}$, drawn from the 70% complete sample of the Arecibo Legacy Fast ALFA survey. We construct subsamples of galaxies with $M_{\rm HI}$ above different thresholds, and make volume-limited clustering measurements in terms of three statistics: the projected two-point correlation function, the projected cross-correlation function with respect to a reference sample selected from the Sloan Digital Sky Survey, and the redshift-space monopole moment. In contrast to previous studies, which found no/weak HI-mass dependence, we find both the clustering amplitude on scales above a few Mpc and the bias factors to increase significantly with increasing HI mass for subsamples with HI mass thresholds above $10^9M_{\odot}$. For HI mass thresholds below $10^9M_{\odot}$, while the measurements have large uncertainties caused by the limited survey volume and sample size, the inferred galaxy bias factors are systematically lower than the minimum halo bias factor from mass-selected halo samples. The simple halo model, in which galaxy content is only determined by halo mass, has difficulties in interpreting the clustering measurements of the HI-selected samples. We extend the simple model by including the halo formation time as an additional parameter. A model that puts HI-rich galaxies into halos that formed late can reproduce the clustering measurements reasonably well. We present the implications of our best-fitting model on the correlation of HI mass with halo mass and formation time, as well as the halo occupation distributions and HI mass functions for central and satellite galaxies. These results are compared with the predictions from semi-analytic galaxy formation models and hydrodynamic galaxy formation simulations.


title: Shared-memory Graph Truss Decomposition

abstract: We present PKT, a new shared-memory parallel algorithm and OpenMP implementation for the truss decomposition of large sparse graphs. A k-truss is a dense subgraph definition that can be considered a relaxation of a clique. Truss decomposition refers to a partitioning of all the edges in the graph based on their k-truss membership. The truss decomposition of a graph has many applications. We show that our new approach PKT consistently outperforms other truss decomposition approaches for a collection of large sparse graphs and on a 24-core shared-memory server. PKT is based on a recently proposed algorithm for k-core decomposition.


title: Nonequilibrium-induced enhancement of dynamical quantum coherence and entanglement of spin arrays

abstract: The random magnetic field produced by nuclear spins has long been viewed as the dominating source of decoherence in the quantum-dot based spins. Here we obtain in both exact and analytical manner the dynamics of spin qubits coupled to nuclear spin environments via the hyperfine interaction, going beyond the weak system-bath interaction and Markovian approximation. We predict that the detailed-balance breaking produced by chemical potential gradient in nuclear baths leads to the rapid oscillations of populations, quantum coherence and entanglement, which are absent in the conventional case (i.e., Overhauser noise). This is attributed to the nonequilibrium feature of the system as shown in the relation between the oscillation period and the chemical potential imbalance. Our results reveal the essentiality of nonequilibriumness with detailed-balance breaking for enhancing the dynamical coherence and entanglement of spin qubits. Moreover, our exact solution explicitly demonstrates that the non-Markovian bath comprised by nuclear spins can preserve the collective quantum state, due to the recovery of coherence. Finally, we propose an experiment using ultracold trapped ions to observe these nonequilibrium and memory effects.


title: Further results on the expected hitting time, the cover cost and the related invariants of graphs

abstract: A close relation between hitting times of the simple random walk on a graph, the Kirchhoff index, resistance-centrality, and related invariants of unicyclic graphs is displayed. Combining with the graph transformations and some other techniques, sharp upper and lower bounds on the cover cost (resp. reverse cover cost) of a vertex in an $n$-vertex unicyclic graph are determined. All the corresponding extremal graphs are identified, respectively.


title: Infinite families of 2-designs from GA_1(q) actions

abstract: Group action is a standard approach to obtain $t$-designs. In this approach, selecting a specific permutation group with a certain degree of transitivity or homogeneity and a proper set of base blocks is important for obtaining $t$-$(v, k, \lambda)$ designs with computable parameters $t, v, k$, and $\lambda$. The general affine group $\GA_1(q)$ is $2$-transitive on $\gf(q)$, and has relatively a small size. In this paper, we determine the parameters of a number of infinite families of $2$-designs obtained from the action of the group $\GA_1(q)$ on certain base blocks, and demonstrate that some of the $2$-designs give rise to linear codes with optimal or best parameters known. Open problems are also presented.


title: Stable cosmology in ghost-free quasidilaton theory

abstract: We present a novel cosmological solution in the framework of extended quasidilaton theory which underwent scrutiny recently. We only consider terms that do not generate the Boulware-Deser degree of freedom, hence the "ghost-free" quasidilaton theory, and show three new branches of cosmological evolution therein. One of the solutions passes the perturbative stability tests. This new solution exhibits a late time self-acceleration and all graviton polarizations acquire masses that converge to a constant in the asymptotic future. Moreover, all modes propagate at the speed of light. We propose that this solution can be used as a benchmark model for future phenomenological studies.


title: The hydrodynamic limit of a randomized load balancing network

abstract: Randomized load balancing networks arise in a variety of applications, and allow for efficient sharing of resources, while being relatively easy to implement. We consider a network of parallel queues in which incoming jobs with independent and identically distributed service times are assigned to the shortest queue among a randomly chosen subset of $d$ queues, and leave the network on completion of service. Prior work on dynamical properties of this model has focused on the case of exponential service distributions. In this work, we analyze the more realistic case of general service distributions. We first introduce a novel particle representation of the state of the network, and characterize the state dynamics via a sequence of interacting measure-valued stochastic processes. Under mild assumptions, we show that the sequence of scaled state processes converges, as the number of servers goes to infinity, to a hydrodynamic limit that is characterized as the unique solution to a countable system of coupled deterministic measure-valued equations. We also establish a propagation of chaos result that shows that finite collections of queues are asymptotically independent. The general framework developed here is potentially useful for analyzing a larger class of models arising in diverse fields including biology and materials science.


title: Highly excited exciton-polariton condensates

abstract: Exciton-polaritons are a coherent electron-hole-photon (e-h-p) system where condensation has been observed in semiconductor microcavities. In contrast to equilibrium Bose-Einstein condensation (BEC) for long lifetime systems, polariton condensates have a dynamical nonequilibrium feature owing to the similar physical structure that they have to semiconductor lasers. One of the distinguishing features of a condensate to a laser is the presence of strong coupling between the matter and photon fields. Irrespective of its equilibrium or nonequilibrium nature, exciton-polariton have been observed to maintain strong coupling. We show that by investigating high density regime of exciton-polariton condensates, the negative branch directly observed in photoluminescence. This is evidence that the present e-h-p system is still in the strong coupling regime, contrary to past results where the system reduced to standard lasing at high density.


title: Truncated $\mathcal{V}$-fractional Taylor's formula with applications

abstract: In this paper, we present and prove a new truncated $\mathcal{V}$-fractional Taylor's formula using the truncated $\mathcal{V}$-fractional variation of constants formula. In this sense, we present the truncated $\mathcal{V}$-fractional Taylor's remainder by means of $\mathcal{V}$-fractional integral, essential for analyzing and comparing the error, when approaching functions by polynomials. From these new results, some applications were made involving some inequalities, specifically, we generalize the Cauchy-Schwartz inequality.


title: Separating CMB Stokes Q and U polarization signals from Non-Gaussian Emissions

abstract: In this paper we estimate diffuse foreground minimized Cosmic Microwave Background (CMB) Stokes Q and U polarization maps based upon the fundamental concept of Gaussian nature of CMB and strong non-Gaussian nature of astrophysical polarized foregrounds using WMAP nine year published polarization maps. We excise regions of the sky that define position of the known point sources, regions that are strongly contaminated by either the detector noise or by the diffuse foregrounds or both, and then perform foreground minimizations over the surviving sky regions that constitute approximately $50\%$ of the full sky area. We critically evaluate performance of foreground minimizations in several ways and show that our foreground minimization method removes significant foregrounds from input maps. The cleaned Stokes \{Q, U\} polarization maps have less EE and BB power from relevant sky region compared to WMAP foreground-reduced Stokes \{Q, U\} polarization maps at different multipole ranges. We validate our methodology by performing detailed Monte Carlo simulations. The main driving machinery of our method is an internal-linear-combination (ILC) approach, however, unlike simple variance minimization performed in usual ILC method, the core of the method being dependent on the theoretically well motivated concept of Gaussianity of CMB polarization a direct connection is established between observations and models of inflation. Additionally, the method, like the usual ILC method, is independent on modeling uncertainties of polarized foregrounds. It will be useful to apply our method in future generation low-noise CMB polarization experiments.


title: Coherent Emission Mechanisms in Astrophysical Plasmas

abstract: Three known examples of coherent emission in radio astronomical sources are reviewed: plasma emission, electron cyclotron maser emission (ECME) and pulsar radio emission.


title: The totally nonnegative Grassmannian is a ball

abstract: We prove that three spaces of importance in topological combinatorics are homeomorphic to closed balls: the totally nonnegative Grassmannian, the compactification of the space of electrical networks, and the cyclically symmetric amplituhedron.


title: Lifting CDCL to Template-based Abstract Domains for Program Verification

abstract: The success of Conflict Driven Clause Learning (CDCL) for Boolean satisfiability has inspired adoption in other domains. We present a novel lifting of CDCL to program analysis called Abstract Conflict Driven Learning for Programs (ACDLP). ACDLP alternates between model search, which performs over-approximate deduction with constraint propagation, and conflict analysis, which performs under-approximate abduction with heuristic choice. We instantiate the model search and conflict analysis algorithms to an abstract domain of template polyhedra, strictly generalizing CDCL from the Boolean lattice to a richer lattice structure. Our template polyhedra can express intervals, octagons and restricted polyhedral constraints over program variables. We have imple- mented ACDLP for automatic bounded safety verification of C programs. We evaluate the performance of our analyser by comparing with CBMC, which uses CDCL, and Astree, a commercial abstract interpretation tool. We observe two orders of magnitude reduction in the number of decisions, propagations, and conflicts as well as a 1.5x speedup in runtime compared to CBMC. Compared to Astree, ACDLP solves twice as many benchmarks and has much higher precision. This is the first instantiation of CDCL with a template polyhedra abstract domain.


title: A multivariate integral representation on $\mathrm{GL}_2 \times \mathrm{GSp}_4$ inspired by the pullback formula

abstract: We give a two variable Rankin-Selberg integral inspired by consideration of Garrett's pullback formula. For a globally generic cusp form on $\mathrm{GL}_2\times \mathrm{GSp}_4$, the integral represents the product of the $\mathrm{Std}\times \mathrm{Spin}$ and $\mathbf{1} \times \mathrm{Std}$ $L$-functions. We prove a result concerning an Archimedean principal series representation in order to verify a case of Jiang's first-term identity relating certain non-Siegel Eisenstein series on symplectic groups. Using it, we obtain a new proof of a known result concerning possible poles of these $L$-functions.


title: Global well-posedness of the periodic cubic fourth order NLS in negative Sobolev spaces

abstract: We consider the Cauchy problem for the cubic fourth order nonlinear Schr\"odinger equation (4NLS) on the circle. In particular, we prove global well-posedness of the renormalized 4NLS in negative Sobolev spaces $H^s(\mathbb{T})$, $s > -\frac{1}{3}$, with enhanced uniqueness. The proof consists of two separate arguments. (i) We first prove global existence in $H^s(\mathbb{T})$, $s > -\frac{9}{20}$, via the short-time Fourier restriction norm method. By following the argument in Guo-Oh for the cubic NLS, this also leads to non-existence of solutions for the (non-renormalized) 4NLS in negative Sobolev spaces. (ii) We then prove enhanced uniqueness in $H^s(\mathbb{T})$, $s > -\frac{1}{3}$, by establishing an energy estimate for the difference of two solutions with the same initial condition. For this purpose, we perform an infinite iteration of normal form reductions on the $H^s$-energy functional, allowing us to introduce an infinite sequence of correction terms to the $H^s$-energy functional in the spirit of the $I$-method. In fact, the main novelty of this paper is this reduction of the $H^s$-energy functionals (for a single solution and for the difference of two solutions with the same initial condition) to sums of infinite series of multilinear terms of increasing degrees.


title: A Robust t-process Regression Model with Independent Errors

abstract: Gaussian process regression (GPR) model is well-known to be susceptible to outliers. Robust process regression models based on t-process or other heavy-tailed processes have been developed to address the problem. However, due to the nature of the current definition for heavy-tailed processes, the unknown process regression function and the random errors are always defined jointly and thus dependently. This definition, mainly owing to the dependence assumption involved, is not justified in many practical problems and thus limits the application of those robust approaches. It also results in a limitation of the theory of robust analysis. In this paper, we propose a new robust process regression model enabling independent random errors. An efficient estimation procedure is developed. Statistical properties, such as unbiasness and information consistency, are provided. Numerical studies show that the proposed method is robust against outliers and has a better performance in prediction compared with the existing models. We illustrate that the estimated random-effects are useful in detecting outlying curves.


title: Spin-polarized electric current in silicene nanoribbons induced by atomic adsorption

abstract: We investigate the nonequilibrium transport properties of a silicene armchair nanoribbon with a random distribution of adsorbed atoms in apex positions. A ferromagnetic insulator grown below the nanoribbon splits spin-up and spin-down electron bands and gives rise to a spin polarization of the conductance. The conductance vanishes when the Fermi energy matches the adatom levels due to the coupling of adatom localized states with the continuum spectra of the nanoribbon. This is the well-known Fano effect, resulting in a spin-dependent antiresonance in the conductance. The different antiresonance energies of spin-up and spin-down electrons give rise to a full spin polarization of the conductance in a broad energy window. This spin-dependent Fano effect opens the possibility to using it in spintronics as a tuneable source of polarized electrons.


title: Asymptotic stability of the stationary Navier-Stokes flows in Besov spaces

abstract: We discuss the asymptotic stability of stationary solutions to the incompressible Navier-Stokes equations on the whole space in Besov spaces with positive smoothness and low integrability. A critical estimate for the semigroup generated by the Laplacian with a perturbation is main ingredient of the argument.


title: Fano varieties with large Seshadri constants

abstract: We show that the set of Fano varieties (with arbitrary singularities) whose anticanonical divisors have large Seshadri constants satisfies certain weak and birational boundedness. We also classify singular Fano varieties of dimension $n$ whose anticanonical divisors have Seshadri constants at least $n$, generalizing an earlier result of Liu and the author.


title: Efficient Adjoint Computation for Wavelet and Convolution Operators

abstract: First-order optimization algorithms, often preferred for large problems, require the gradient of the differentiable terms in the objective function. These gradients often involve linear operators and their adjoints, which must be applied rapidly. We consider two example problems and derive methods for quickly evaluating the required adjoint operator. The first example is an image deblurring problem, where we must compute efficiently the adjoint of multi-stage wavelet reconstruction. Our formulation of the adjoint works for a variety of boundary conditions, which allows the formulation to generalize to a larger class of problems. The second example is a blind channel estimation problem taken from the optimization literature where we must compute the adjoint of the convolution of two signals. In each example, we show how the adjoint operator can be applied efficiently while leveraging existing software.


title: Option Pricing and Hedging for Discrete Time Autoregressive Hidden Markov Model

abstract: In this paper we solve the discrete time mean-variance hedging problem when asset returns follow a multivariate autoregressive hidden Markov model. Time dependent volatility and serial dependence are well established properties of financial time series and our model covers both. To illustrate the relevance of our proposed methodology, we first compare the proposed model with the well-known hidden Markov model via likelihood ratio tests and a novel goodness-of-fit test on the S\&P 500 daily returns. Secondly, we present out-of-sample hedging results on S\&P 500 vanilla options as well as a trading strategy based on theoretical prices, which we compare to simpler models including the classical Black-Scholes delta-hedging approach.


title: Some ergodic properties of metrics on hyperbolic groups

abstract: Let $\Gamma$ be a non-elementary Gromov-hyperbolic group, and $\partial \Gamma$ denote its Gromov boundary. We consider $\Gamma$-invariant proper $\delta$-hyperbolic, quasi-convex metric $d$ on $\Gamma$, and the associated Patterson-Sullivan measure class $[\nu]$ on $\partial^{(2)}\Gamma$, and its square $[\nu\times\nu]$ on $\partial^{(2)}\Gamma$ -- the space of distinct pairs of points on the boundary. We construct an analogue of a geodesic flow to study ergodicity properties of the $\Gamma$-actions on $(\partial\Gamma,\nu)$ and on $(\partial^{(2)}\Gamma,[\nu\times\nu])$.


title: SpArcFiRe: morphological selection effects due to reduced visibility of tightly winding arms in distant spiral galaxies

abstract: The Galaxy Zoo has provided morphological data on many galaxies. Several biases have been identified in the Galaxy Zoo data. Here we report on a newly discovered selection effect: astronomers interested in studying spiral galaxies may select a set of spiral galaxies based upon a threshold in spirality (the fraction of Galaxy Zoo humans who report seeing spiral structure). SpArcFiRe is an automated tool that decomposes a spiral galaxy into its constituent spiral arms, providing objective, quantitative data on their structure. SpArcFiRe measures the pitch angle of spiral arms. We have observed that when selecting a set of spiral galaxies based on a threshold on spirality, the pitch angle of spiral arms appear increase with redshift. We hypothesize that this is a selection effect: tightly-wound spiral arms become less visible as images degrade with increasing redshift, leading to fewer such galaxies being included in the sample at higher redshifts. We corroborate this hypothesis by artificially degrading images of nearby galaxies, then using a machine learning algorithm trained on Galaxy Zoo data to provide a spirality for each artificially degraded image. It correctly predicts that spirality decreases as image quality degrades. Thus, the mean pitch angle of those galaxies remaining above the spirality threshold is higher than those eliminated by the selection effect. This demonstrates that users who select samples of galaxies using a threshold of Galaxy Zoo votes must carefully consider the possibility of selection effects on morphological measures, even if the measure itself is believed to be objective and unbiased. Finally, we also perform an empirical sensitivity analysis to demonstrate that SpArcFiRe's output changes in a smooth and predictable fashion to changes in its internal algorithmic parameters.


title: Automatic Classification of Bright Retinal Lesions via Deep Network Features

abstract: The diabetic retinopathy is timely diagonalized through color eye fundus images by experienced ophthalmologists, in order to recognize potential retinal features and identify early-blindness cases. In this paper, it is proposed to extract deep features from the last fully-connected layer of, four different, pre-trained convolutional neural networks. These features are then feeded into a non-linear classifier to discriminate three-class diabetic cases, i.e., normal, exudates, and drusen. Averaged across 1113 color retinal images collected from six publicly available annotated datasets, the deep features approach perform better than the classical bag-of-words approach. The proposed approaches have an average accuracy between 91.23% and 92.00% with more than 13% improvement over the traditional state of art methods.


title: A flux calibration method for remote sensing satellites using stars

abstract: Star surveys and model analyses show that many stars have absolute stable fluxes as good as 3% in 0.3-35{\mu}m wavebands and about 1% in the visible wavebands. The relative flux calibrations between stars are better than 0.2%. Some stars have extremely stable fluxes and can be used as long term flux calibration sources. Stellar brightness is several orders of magnitude lower than most ground objects while the stars do not usually appear in remote sensing cameras, which makes the stars inappropriate for being calibration sources. The calibration method using stars discussed in this paper is through a mini-camera attached to remote sensing satellite. The mini-camera works at similar wavebands as the remote sensing cameras and it can observe the stars and the ground objects alternatively. High signal-to-noise ratio is achieved for the relatively faint stars through longer exposure time. Simultaneous precise cross-calibration is obtained as the mini-camera and remote sensing cameras look at the ground objects at the same time. The fluxes from the stars used as calibration standards are transferred to the remote sensing cameras through this procedure. Analysis shows that a 2% accurate calibration is possible.


title: Comparison of Flow Scheduling Policies for Mix of Regular and Deadline Traffic in Datacenter Environments

abstract: Datacenters are the main infrastructure on top of which cloud computing services are offered. Such infrastructure may be shared by a large number of tenants and applications generating a spectrum of datacenter traffic. Delay sensitive applications and applications with specific Service Level Agreements (SLAs), generate deadline constrained flows, while other applications initiate flows that are desired to be delivered as early as possible. As a result, datacenter traffic is a mix of two types of flows: deadline and regular. There are several scheduling policies for either traffic type with focus on minimizing completion times or deadline miss rate. In this report, we apply several scheduling policies to mix traffic scenario while varying the ratio of regular to deadline traffic. We consider FCFS (First Come First Serve), SRPT (Shortest Remaining Processing Time) and Fair Sharing as deadline agnostic approaches and a combination of Earliest Deadline First (EDF) with either FCFS or SRPT as deadline-aware schemes. In addition, for the latter, we consider both cases of prioritizing deadline traffic (Deadline First) and prioritizing regular traffic (Deadline Last). We study both light-tailed and heavy-tailed flow size distributions and measure mean, median and tail flow completion times (FCT) for regular flows along with Deadline Miss Rate (DMR) and average lateness for deadline flows. We also consider two operation regimes of lightly-loaded (low utilization) and heavily-loaded (high utilization). We find that performance of deadline-aware schemes is highly dependent on fraction of deadline traffic. With light-tailed flow sizes, we find that FCFS performs better in terms of tail times and average lateness while SRPT performs better in average times and deadline miss rate. For heavy-tailed flow sizes, except for tail times, SRPT performs better in all other metrics.


title: Improving the staggered grid Lagrangian hydrodynamics for modeling multi-material flows

abstract: In this work, we make two improvements on the staggered grid hydrodynamics (SGH) Lagrangian scheme for modeling 2-dimensional compressible multi-material flows on triangular mesh. The first improvement is the construction of a dynamic local remeshing scheme for preventing mesh distortion. The remeshing scheme is similar to many published algorithms except that it introduces some special operations for treating grids around multi-material interfaces. This makes the simulation of extremely deforming and topology-variable multi-material processes possible, such as the complete process of a heavy fluid dipping into a light fluid. The second improvement is the construction of an Euler-like flow on each edge of the mesh to count for the "edge-bending" effect, so as to mitigate the "checkerboard" oscillation that commonly exists in Lagrangian simulations, especially the triangular mesh based simulations. Several typical hydrodynamic problems are simulated by the improved staggered grid Lagrangian hydrodynamic method to test its performance.


title: A Nested Attention Neural Hybrid Model for Grammatical Error Correction

abstract: Grammatical error correction (GEC) systems strive to correct both global errors in word order and usage, and local errors in spelling and inflection. Further developing upon recent work on neural machine translation, we propose a new hybrid neural model with nested attention layers for GEC. Experiments show that the new model can effectively correct errors of both types by incorporating word and character-level information,and that the model significantly outperforms previous neural models for GEC as measured on the standard CoNLL-14 benchmark dataset. Further analysis also shows that the superiority of the proposed model can be largely attributed to the use of the nested attention mechanism, which has proven particularly effective in correcting local errors that involve small edits in orthography.


title: DDCCast: Meeting Point to Multipoint Transfer Deadlines Across Datacenters using ALAP Scheduling Policy

abstract: Large cloud companies manage dozens of datacenters across the globe connected using dedicated inter-datacenter networks. An important application of these networks is data replication which is done for purposes such as increased resiliency via making backup copies, getting data closer to users for reduced delay and WAN bandwidth usage, and global load balancing. These replications usually lead to network transfers with deadlines that determine the time prior to which all datacenters should have a copy of the data. Inter-datacenter networks have limited capacity and need be utilized efficiently to maximize performance. In this report, we focus on applications that transfer multiple copies of objects from one datacenter to several datacenters given deadline constraints. Existing solutions are either deadline agnostic, or only consider point-to-point transfers. We propose DDCCast, a simple yet effective deadline aware point to multipoint technique based on DCCast and using ALAP traffic allocation. DDCCast performs careful admission control using temporal planning, uses rate-allocation and rate-limiting to avoid congestion and sends traffic over forwarding trees that are carefully selected to reduce bandwidth usage and maximize deadline meet rate. We perform experiments confirming DDCCast's potential to reduce total bandwidth usage by up to $45\%$ while admitting up to $25\%$ more traffic into the network compared to existing solutions that guarantee deadlines.


title: More efficient second harmonic generation of whispering gallery modes by selective out-coupling

abstract: We demonstrate second harmonic generation (SHG) in an $x$-cut congruent lithium niobate (LN) whispering gallery mode resonator. We first show theoretically that independent control of the coupling of the pump and signal modes is optimal for high conversion rates. A scheme based on our earlier work in Ref. [1] is then implemented experimentally to verify this. Thereby we are able to improve on the efficiency of SHG by more than an order of magnitude by selectively out-coupling using a LN prism, utilizing the birefringence of it and the resonator in kind. We report 5.28%/mW efficiency for SHG from 1555.4 nm to 777.7 nm.


title: LoopInvGen: A Loop Invariant Generator based on Precondition Inference

abstract: We describe the LoopInvGen tool for generating loop invariants that can provably guarantee correctness of a program with respect to a given specification. LoopInvGen is an efficient implementation of the inference technique originally proposed in our earlier work on PIE (https://doi.org/10.1145/2908080.2908099). In contrast to existing techniques, LoopInvGen is not restricted to a fixed set of features -- atomic predicates that are composed together to build complex loop invariants. Instead, we start with no initial features, and use program synthesis techniques to grow the set on demand. This not only enables a less onerous and more expressive approach, but also appears to be significantly faster than the existing tools over the SyGuS-COMP 2018 benchmarks from the INV track.


title: Is QBism the Future of Quantum Physics?

abstract: The purpose of this book is to explain Quantum Bayesianism ('QBism') to "people without easy access to mathematical formulas and equations" (4-5). Qbism is an interpretation of quantum mechanics that "doesn't meddle with the technical aspects of the theory [but instead] reinterprets the fundamental terms of the theory and gives them new meaning" (3). The most important motivation for QBism, enthusiastically stated on the book's cover, is that QBism provides "a way past quantum theory's paradoxes and puzzles" such that much of the weirdness associated with quantum theory "dissolves under the lens of QBism". Could a non-technical book that almost fits in your pocket really succeed in resolving the notorious paradoxes of quantum theory? I believe the answer is: not in this case. There are three primary reasons. Firstly, the argument that QBism solves quantum paradoxes is not convincing. Secondly, it proves difficult to pin down exactly what QBism says. Thirdly, as a scientific theory QBism seems explanatorily inert, which is maybe why the topic of scientific explanation is not broached. Sections 2-4 below discuss each of these problems in turn. However, there is still a wealth of insights that make this book a worthwhile read, and I shall begin with those.


title: Integral Ricci curvature bounds for possibly collapsed spaces with Ricci curvature bounded from below

abstract: Assuming a lower bound on the Ricci curvature of a complete Riemannian manifold, for $p< 1/2$ we show the existence of bounds on the local $L^p$ norm of the Ricci curvature that depend only on the dimension and which improve with volume collapse.


title: Matrix-Based Characterization of the Motion and Wrench Uncertainties in Robotic Manipulators

abstract: Characterization of the uncertainty in robotic manipulators is the focus of this paper. Based on the random matrix theory (RMT), we propose uncertainty characterization schemes in which the uncertainty is modeled at the macro (system) level. This is different from the traditional approaches that model the uncertainty in the parametric space of micro (state) level. We show that perturbing the system matrices rather than the state of the system provides unique advantages especially for robotic manipulators. First, it requires only limited statistical information that becomes effective when dealing with complex systems where detailed information on their variability is not available. Second, the RMT-based models are aware of the system state and configuration that are significant factors affecting the level of uncertainty in system behavior. In this study, in addition to the motion uncertainty analysis that was first proposed in our earlier work, we also develop an RMT-based model for the quantification of the static wrench uncertainty in multi-agent cooperative systems. This model is aimed to be an alternative to the elaborate parametric formulation when only rough bounds are available on the system parameters. We discuss that how RMT-based model becomes advantageous when the complexity of the system increases. We perform experimental studies on a KUKA youBot arm to demonstrate the superiority of the RMT-based motion uncertainty models. We show that how these models outperform the traditional models built upon Gaussianity assumption in capturing real-system uncertainty and providing accurate bounds on the state estimation errors. In addition, to experimentally support our wrench uncertainty quantification model, we study the behavior of a cooperative system of mobile robots. It is shown that one can rely on less demanding RMT-based formulation and yet meets the acceptable accuracy.


title: Networked Fairness in Cake Cutting

abstract: We introduce a graphical framework for fair division in cake cutting, where comparisons between agents are limited by an underlying network structure. We generalize the classical fairness notions of envy-freeness and proportionality to this graphical setting. Given a simple undirected graph G, an allocation is envy-free on G if no agent envies any of her neighbor's share, and is proportional on G if every agent values her own share no less than the average among her neighbors, with respect to her own measure. These generalizations open new research directions in developing simple and efficient algorithms that can produce fair allocations under specific graph structures. On the algorithmic frontier, we first propose a moving-knife algorithm that outputs an envy-free allocation on trees. The algorithm is significantly simpler than the discrete and bounded envy-free algorithm recently designed by Aziz and Mackenzie for complete graphs. Next, we give a discrete and bounded algorithm for computing a proportional allocation on descendant graphs, a class of graphs by taking a rooted tree and connecting all its ancestor-descendant pairs.


title: Complete Call-by-Value Calculi of Control Operators II: Strong Termination

abstract: We provide characterization of the strong termination property of the CCV (complete call-by-value) lambda-mu calculus introduced in the first part of this series of the paper. The calculus is complete with respect to the standard continuation-passing style (CPS) semantics. The union-intersection type systems for the calculus is developed in the previous paper. We characterize the strong normalizability of terms of the calculus in terms of the CPS semantics and typeability.


title: The Formation of Massive Molecular Filaments and Massive Stars Triggered by a MHD Shock Wave

abstract: Recent observations suggest that intensive molecular cloud collision can trigger massive star/cluster formation. The most important physical process caused by the collision is a shock compression. In this paper, the influence of a shock wave on the evolution of a molecular cloud is studied numerically by using isothermal magnetohydrodynamics (MHD) simulations with the effect of self-gravity. Adaptive-mesh-refinement and sink particle techniques are used to follow long-time evolution of the shocked cloud. We find that the shock compression of turbulent inhomogeneous molecular cloud creates massive filaments, which lie perpendicularly to the background magnetic field as we have pointed out in a previous paper. The massive filament shows global collapse along the filament, which feeds a sink particle located at the collapse center. We observe high accretion rate dot{M}_acc > 10^{-4} M_sun/yr that is high enough to allow the formation of even O-type stars. The most massive sink particle achieves M>50 M_sun in a few times 10^5 yr after the onset of the filament collapse.


title: Rationally Equivalent Points On Hypersurfaces In ${\mathbb P}^n$

abstract: We prove a conjecture of Voisin that no two distinct points on a very general hypersurface of degree $2n$ in ${\mathbb P}^n$ are rationally equivalent.


title: Stanley sequences with odd character

abstract: Given a set of integers containing no 3-term arithmetic progressions, one constructs a Stanley sequence by choosing integers greedily without forming such a progression. Independent Stanley sequences are a "well-structured" class of Stanley sequences with two main parameters: the character $\lambda(A)$ and the repeat factor $\rho(A)$. Rolnick conjectured that for every $\lambda \in \mathbb{N}_0\backslash\{1, 3, 5, 9, 11, 15\}$, there exists an independent Stanley sequence $S(A)$ such that $\lambda(A) =\lambda$. This paper demonstrates that $\lambda(A) \not\in \{1, 3, 5, 9, 11, 15\}$ for any independent Stanley sequence $S(A)$.


title: A Tutorial on Thompson Sampling

abstract: Thompson sampling is an algorithm for online decision problems where actions are taken sequentially in a manner that must balance between exploiting what is known to maximize immediate performance and investing to accumulate new information that may improve future performance. The algorithm addresses a broad range of problems in a computationally efficient manner and is therefore enjoying wide use. This tutorial covers the algorithm and its application, illustrating concepts through a range of examples, including Bernoulli bandit problems, shortest path problems, product recommendation, assortment, active learning with neural networks, and reinforcement learning in Markov decision processes. Most of these problems involve complex information structures, where information revealed by taking an action informs beliefs about other actions. We will also discuss when and why Thompson sampling is or is not effective and relations to alternative algorithms.


title: A note on some variations of the $\gamma$-graph

abstract: For a graph $G$, the $\gamma$-graph of $G$, $G(\gamma)$, is the graph whose vertices correspond to the minimum dominating sets of $G$, and where two vertices of $G(\gamma)$ are adjacent if and only if their corresponding dominating sets in $G$ differ by exactly two adjacent vertices. In this paper, we present several variations of the $\gamma$-graph including those using identifying codes, locating-domination, total-domination, paired-domination, and the upper-domination number. For each, we show that for any graph $H$, there exist infinitely many graphs whose $\gamma$-graph variant is isomorphic to $H$.


title: PT-symmetry in nonlinear twisted multi-core fibers

abstract: We address propagation of light in nonlinear twisted multi-core fibers with alternating amplifying and absorbing cores that are arranged into the PT - symmetric structure. In this structure, the coupling strength between neighboring cores and global energy transport can be controlled not only by the nonlinearity strength, but also by gain and losses and by the fiber twisting rate. The threshold level of gain/losses, at which PT -symmetry breaking occurs, is a non-monotonic function of the fiber twisting rate and it can be reduced nearly to zero or, instead, notably increased just by changing this rate. Nonlinearity usually leads to the monotonic reduction of the symmetry breaking threshold in such fibers.


title: DroneCells: Improving 5G Spectral Efficiency using Drone-mounted Flying Base Stations

abstract: We study a cellular networking scenario, called DroneCells, where miniaturized base stations (BSs) are mounted on flying drones to serve mobile users. We propose that the drones never stop, and move continuously within the cell in a way that reduces the distance between the BS and the serving users, thus potentially improving the spectral efficiency of the network. By considering the practical mobility constraints of commercial drones, we design drone mobility algorithms to improve the spectral efficiency of DroneCells. As the optimal problem is NP-hard, we propose a range of practically realizable heuristics with varying complexity and performance. Simulations show that, using the existing consumer drones, the proposed algorithms can readily improve spectral efficiency by 34\% and the 5-percentile packet throughput by 50\% compared to the scenario, where drones hover over fixed locations. More significant gains can be expected with more agile drones in the future. A surprising outcome is that the drones need to fly only at minimal speeds to achieve these gains, avoiding any negative effect on drone battery lifetime. We further demonstrate that the optimal solution provides only modest improvements over the best heuristic algorithm, which employs Game Theory to make mobility decisions for drone BSs.


title: PHL 6625: A Minor Merger-Associated QSO Behind NGC 247

abstract: PHL 6625 is a luminous quasi-stellar object (QSO) at z = 0.3954 located behind the nearby galaxy NGC 247 (z = 0.0005). Hubble Space Telescope (HST) observations revealed an arc structure associated with it. We report on spectroscopic observations with the Very Large Telescope (VLT) and multiwavelength observations from the radio to the X-ray band for the system, suggesting that PHL 6625 and the arc are a close pair of merging galaxies, instead of a strong gravitational lens system. The QSO host galaxy is estimated to be (4-28) x 10^10 M_sun, and the mass of the companion galaxy of is estimated to be M_* = (6.8 +/- 2.4) x 10^9 M_sun, suggesting that this is a minor merger system. The QSO displays typical broad emission lines, from which a black hole mass of about (2-5) x 10^8 M_sun and an Eddington ratio of about 0.01-0.05 can be inferred. The system represents an interesting and rare case where a QSO is associated with an ongoing minor merger, analogous to Arp 142.


title: Weakly distance-regular digraphs whose attached association schemes are regular

abstract: In this paper, we study commutative weakly distance-regular digraphs whose attached association schemes are regular, and give a characterization of mixed arcs. As an application, we classify such digraphs of diameter 2.


title: Dissipative Stability Conditions for Linear Coupled Differential-Difference Systems via a Dynamical Constraints Approach

abstract: In this short note, we derive dissipative conditions with slack variables for a linear coupled differential-difference (CDDS) via constructing a Krasovskii functional. The approach can be interpreted as a generalization of the Finsler Lemma approach for standard LTI systems proposed previously in \cite{de2001stability}. We also show that the proposed slack variables scheme is equivalent to the approach based on directly substituting the system trajectory $\dot{\bm{x}}(t)$, similar to the case of LTI system.


title: The SuperWASP catalogue of 4963 RR Lyr stars: identification of 983 Blazhko effect candidates

abstract: We set out to compile a catalogue of RRab pulsating variables in the SuperWASP archive and identify candidate Blazhko effect objects within this catalogue. We analysed their light curves and power spectra for correlations in their common characteristics to further our understanding of the phenomenon. Pulsation periods were found for each SWASP RRab object using PDM techniques. Low frequency periodic signals detected in the CLEAN power spectra of RRab stars were matched with modulation sidebands and combined with pairs of sidebands to produce a list of candidate Blazhko periods. A novel technique was used in an attempt to identify Blazhko effect stars by comparing scatter at different parts of the folded light curve. Pulsation amplitudes were calculated based on phase folded light curves. The SuperWASP RRab catalogue consists of 4963 objects of which 3397 are previously unknown. We discovered 983 distinct candidates for Blazhko effect objects, 613 of these being previously unknown in the literature as RR Lyrae stars, and 894 are previously unknown to be Blazhko effect stars. Correlations were investigated between the scatter of points on the light curve, the periods and amplitudes of the objects' pulsations, and those of the Blazhko effect. A statistical analysis has been performed on a large population of Blazhko effect stars from the wide-field SuperWASP survey. No correlations were found between the Blazhko period and other parameters including the Blazhko amplitude, although we confirmed a lower rate of occurrence of the Blazhko effect in long pulsation period objects.


title: Microwave photonics with superconducting quantum circuits

abstract: In the past 20 years, impressive progress has been made both experimentally and theoretically in superconducting quantum circuits, which provide a platform for manipulating microwave photons. This emerging field of superconducting quantum microwave circuits has been driven by many new interesting phenomena in microwave photonics and quantum information processing. For instance, the interaction between superconducting quantum circuits and single microwave photons can reach the regimes of strong, ultra-strong, and even deep-strong coupling. Many higher-order effects, unusual and less familiar in traditional cavity quantum electrodynamics with natural atoms, have been experimentally observed, e.g., giant Kerr effects, multi-photon processes, and single-atom induced bistability of microwave photons. These developments may lead to improved understanding of the counterintuitive properties of quantum mechanics, and speed up applications ranging from microwave photonics to superconducting quantum information processing. In this article, we review experimental and theoretical progress in microwave photonics with superconducting quantum circuits. We hope that this global review can provide a useful roadmap for this rapidly developing field.


title: InferSpark: Statistical Inference at Scale

abstract: The Apache Spark stack has enabled fast large-scale data processing. Despite a rich library of statistical models and inference algorithms, it does not give domain users the ability to develop their own models. The emergence of probabilistic programming languages has showed the promise of developing sophisticated probabilistic models in a succinct and programmatic way. These frameworks have the potential of automatically generating inference algorithms for the user defined models and answering various statistical queries about the model. It is a perfect time to unite these two great directions to produce a programmable big data analysis framework. We thus propose, InferSpark, a probabilistic programming framework on top of Apache Spark. Efficient statistical inference can be easily implemented on this framework and inference process can leverage the distributed main memory processing power of Spark. This framework makes statistical inference on big data possible and speed up the penetration of probabilistic programming into the data engineering domain.


title: Testing Forecast Accuracy of Expectiles and Quantiles with the Extremal Consistent Loss Functions

abstract: Forecast evaluations aim to choose an accurate forecast for making decisions by using loss functions. However, different loss functions often generate different ranking results for forecasts, which complicates the task of comparisons. In this paper, we develop statistical tests for comparing performances of forecasting expectiles and quantiles of a random variable under consistent loss functions. The test statistics are constructed with the extremal consistent loss functions of Ehm et.al. (2016). The null hypothesis of the tests is that a benchmark forecast at least performs equally well as a competing one under all extremal consistent loss functions. It can be shown that if such a null holds, the benchmark will also perform at least equally well as the competitor under all consistent loss functions. Thus under the null, when different consistent loss functions are used, the result that the competitor does not outperform the benchmark will not be altered. We establish asymptotic properties of the proposed test statistics and propose to use the re-centered bootstrap to construct their empirical distributions. Through simulations, we show the proposed test statistics perform reasonably well. We then apply the proposed method on (1) re-examining abilities of some often-used predictors on forecasting risk premium of the S&P500 index; (2) comparing performances of experts' forecasts on annual growth of U.S. real gross domestic product; (3) evaluating performances of estimated daily value at risk of the S&P500 index.


title: Strain-tuning of edge magnetism in zigzag graphene nanoribbons

abstract: Using the determinant quantum Monte-Carlo method, we elucidate the strain tuning of edge magnetism in zigzag graphene nanoribbons. Our intensive numerical results show that a relatively weak Coulomb interaction may induce a ferromagnetic-like behaviour with a proper strain, and the edge magnetism can be enhanced greatly as the strain along the zigzag edge increases, which provides another way to control graphene magnetism even at room temperature.


title: Exhaustive search for sparse variable selection in linear regression

abstract: We propose a K-sparse exhaustive search (ES-K) method and a K-sparse approximate exhaustive search method (AES-K) for selecting variables in linear regression. With these methods, K-sparse combinations of variables are tested exhaustively assuming that the optimal combination of explanatory variables is K-sparse. By collecting the results of exhaustively computing ES-K, various approximate methods for selecting sparse variables can be summarized as density of states. With this density of states, we can compare different methods for selecting sparse variables such as relaxation and sampling. For large problems where the combinatorial explosion of explanatory variables is crucial, the AES-K method enables density of states to be effectively reconstructed by using the replica-exchange Monte Carlo method and the multiple histogram method. Applying the ES-K and AES-K methods to type Ia supernova data, we confirmed the conventional understanding in astronomy when an appropriate K is given beforehand. However, we found the difficulty to determine K from the data. Using virtual measurement and analysis, we argue that this is caused by data shortage.


title: Image Segmentation Algorithms Overview

abstract: The technology of image segmentation is widely used in medical image processing, face recognition pedestrian detection, etc. The current image segmentation techniques include region-based segmentation, edge detection segmentation, segmentation based on clustering, segmentation based on weakly-supervised learning in CNN, etc. This paper analyzes and summarizes these algorithms of image segmentation, and compares the advantages and disadvantages of different algorithms. Finally, we make a prediction of the development trend of image segmentation with the combination of these algorithms.


title: Structured H-infinity control of infinite dimensional systems

abstract: We develop a novel frequency-based H-infinity control method for a large class of infinite-dimensional Linear-Time-Invariant systems in transfer function form. Major benefits of our approach is that reduction or identification techniques are not needed thereby avoiding possible distortions. It can exploit either transfer function models or input/output frequency response data when available. It computes simple practical controllers of any size and structure. We use a non-smooth trust-region method to compute arbitrarily structured locally optimal H-infinity controllers for a frequency sampled approximation of the underlying infinite-dimensional H-infinity problem in such a way that exponential stability in closed-loop is guaranteed, and the optimal value of the approximation captures the value of the infinite-dimensional problem within a prior tolerance level. We demonstrate the versatility and practicality of our method on a variety of infinite-dimensional H-infinity synthesis problems, including distributed and boundary control of PDEs, control of dead time and delay systems, and using a rich testing set.