title: Natural Gas Flow Solutions with Guarantees: A Monotone Operator Theory Approach

abstract: We consider balanced flows in a natural gas transmission network and discuss computationally hard problems such as establishing if solution of the underlying nonlinear gas flow equations exists, if it is unique, and finding the solution. Particular topologies, e.g. trees, are known to be easy to solve based on a variational description of the gas flow equations, but these approaches do not generalize. In this paper, we show that the gas flow problem can be solved efficiently using the tools of monotone operator theory, provided that we look for solution within certain monotonicity domains. We characterize a family of monotonicity domains, described in terms of Linear Matrix Inequalities (LMI) in the state variables, each containing at most one solution. We also develop an efficient algorithm to choose a particular monotonicity domain, for which the LMI based condition simplifies to a bound on the flows. Performance of the technique is illustrated on exemplary gas networks.


title: A global existence result for a Keller-Segel type system with supercritical initial data

abstract: We consider a parabolic-elliptic Keller-Segel type system, which is related to a simplified model of chemotaxis. Concerning the maximal range of existence of solutions, there are essentially two kinds of results: either global existence in time for general subcritical ($\|\rho_0\|_1<8\pi$) initial data, or blow--up in finite time for suitably chosen supercritical ($\|\rho_0\|_1>8\pi$) initial data with concentration around finitely many points. As a matter of fact there are no results claiming the existence of global solutions in the supercritical case. We solve this problem here and prove that, for a particular set of initial data which share large supercritical masses, the corresponding solution is global and uniformly bounded.


title: Toolbox for continuous variable entanglement production and measurement using spontaneous parametric down conversion

abstract: We provide a toolbox for continuous variables quantum state engineering and characterization of biphoton states produced by spontaneous parametric down conversion in a transverse pump configuration. We show that the control of the pump beam's incidence spot and angle corresponds to phase space displacements of conjugate collective continuous variables of the biphoton. In particular, we illustrate with numerical simulations on a semiconductor device how this technique can be used to engineer and characterize arbitrary states of the frequency and time degrees of freedom.


title: Prominence Visibility in Hinode/XRT images

abstract: In this paper we study the soft X-ray (SXR) signatures of one particular prominence. The X-ray observations used here were made by the Hinode/XRT instrument using two different filters. Both of them have a pronounced peak of the response function around 10 A. One of them has a secondary smaller peak around 170 A, which leads to a contamination of SXR images. The observed darkening in both of these filters has a very large vertical extension. The position and shape of the darkening corresponds nicely with the prominence structure seen in SDO/AIA images. First we have investigated the possibility that the darkening is caused by X-ray absorption. But detailed calculations of the optical thickness in this spectral range show clearly that this effect is completely negligible. Therefore the alternative is the presence of an extended region with a large emissivity deficit which can be caused by the presence of cool prominence plasmas within otherwise hot corona. To reproduce the observed darkening one needs a very large extension along the line-of-sight of the region amounting to around 10$^5$ km. We interpret this region as the prominence spine, which is also consistent with SDO/AIA observations in EUV.


title: On skew polynomial codes and lattices from quotients of cyclic division algebras

abstract: We propose a variation of Construction A of lattices from linear codes defined using the quotient $\Lambda/\mathfrak p\Lambda$ of some order $\Lambda$ inside a cyclic division $F$-algebra, for $\mathfrak p$ a prime ideal of a number field $F$. To obtain codes over this quotient, we first give an isomorphism between $\Lambda/\mathfrak p\Lambda$ and a ring of skew polynomials. We then discuss definitions and basic properties of skew polynomial codes, which are needed for Construction A, but also explore further properties of the dual of such codes. We conclude by providing an application to space-time coding, which is the original motivation to consider cyclic division $F$-algebras as a starting point for this variation of Construction A.


title: Open-point topological games and productivity of dense-separable property

abstract: In this note we study the open-point topological games in order to analyze the least upper bound for density of dense subsets of a topological space. This way we may also analyze the behavior of such cardinal invariants in taking products of spaces. Various related cardinal equalities and inequalities are given. As an application we take a look at Banach spaces with the property (CSP) which can be formulated by stating that each weak-star dense linear subspace of the dual is weak-star separable.


title: A Convergent Gradient Descent Algorithm for Rank Minimization and Semidefinite Programming from Random Linear Measurements

abstract: We propose a simple, scalable, and fast gradient descent algorithm to optimize a nonconvex objective for the rank minimization problem and a closely related family of semidefinite programs. With $O(r^3 \kappa^2 n \log n)$ random measurements of a positive semidefinite $n \times n$ matrix of rank $r$ and condition number $\kappa$, our method is guaranteed to converge linearly to the global optimum.


title: Proper $G_a$-actions on $\C^4$ preserving a coordinate

abstract: We prove that the actions mentioned in the title are translations. We show also that for certain $G_a$-actions on affine fourfolds the quotient of the action is automatically affine and describe the geometric structure of such quotients.


title: Colorings, determinants and Alexander polynomials for spatial graphs

abstract: A {\em balanced} spatial graph has an integer weight on each edge, so that the directed sum of the weights at each vertex is zero. We describe the Alexander module and polynomial for balanced spatial graphs (originally due to Kinoshita \cite{ki}), and examine their behavior under some common operations on the graph. We use the Alexander module to define the determinant and $p$-colorings of a balanced spatial graph, and provide examples. We show that the determinant of a spatial graph determines for which $p$ the graph is $p$-colorable, and that a $p$-coloring of a graph corresponds to a representation of the fundamental group of its complement into a metacyclic group $\Gamma(p,m,k)$. We finish by proving some properties of the Alexander polynomial.


title: The Einstein-Hilbert functional and the Sasaki-Futaki invariant

abstract: We show that the Einstein-Hilbert functional, as a functional on the space of Reeb vector fields, detects the vanishing Sasaki-Futaki invariant. In particular, this provides an obstruction to the existence of a constant scalar curvature Sasakian metric. As an application we prove that K-semistable polarized Sasaki manifold has vanishing Sasaki-Futaki invariant. We then apply this result to show that under the right conditions on the Sasaki join manifolds of [7] a polarized Sasaki manifold is K-semistable if only if it has constant scalar curvature.


title: Some new lacunary $f$-statistical $A$-convergent sequence spaces of order $\alpha$

abstract: We study the concept of density for sets of natural numbers in some lacunary $A$-convergent sequence spaces. Also we are trying to investigate some relation between the ordinary convergence and module statistical convergence for evey unbounded modulus function. Morever we also study some results on the newly defined lacunary $f$-statistically $A$-convergent sequence spaces with respect to some Musielak-Orlicz function.


title: JExtract: An Eclipse Plug-in for Recommending Automated Extract Method Refactorings

abstract: Although Extract Method is a key refactoring for improving program comprehension, refactoring tools for such purpose are often underused. To address this shortcoming, we present JExtract, a recommendation system based on structural similarity that identifies Extract Method refactoring opportunities that are directly automated by IDE-based refactoring tools. Our evaluation suggests that JExtract is far more effective (w.r.t. recall and precision) to identify misplaced code in methods than JDeodorant, a state-of-the-art tool


title: On cycle-supermagic labelings of the disconnected graphs

abstract: In this paper we formulate cycle-supermagic labelings for the disjoint union of isomorphic copies of different families of graphs. We also prove that disjoint union of non isomorphic copies of fans and ladders are cycle-supermagic.


title: Modelling air pollution abatement in deep street canyons by means of air scrubbers

abstract: Deep street canyons are characterized by weak ventilation and recirculation of air. In such environment, the exposure to particulate matter and other air pollutants is enhanced, with a consequent worsening of both safety and health. The main solution adopted by the international community is aimed at the reduction of the emissions. In this theoretical study, we test a new solution: the removal of air pollutants close to their sources by a network of Air Pollution Abatement (APA) devices. The APA technology depletes gaseous and particulate air pollutants by a portable and low-consuming scrubbing system, that mimics the processes of wet and dry deposition. We estimate the potential pollutant abatement efficacy of a single absorber by Computational Fluid Dynamics (CFD) method. The presence of the scrubber effectively creates an additional sink at the bottom of the canyon, accelerating its cleaning process by up to 70%, when an almost perfect scrubber (90% efficiency) is simulated. The efficacy of absorber is not proportional to its internal abatement efficiency, but it increases rapidly at low efficiencies, then tends to saturate. In the particular configuration of the canyon we choose (aspect ratio of 3) the upwind corner is the most favourable for the absorber application. In the downwind corner the maximum pollutant abatement is -24%, while in the upwind corner the maximum abatement is twice as much at -51%. The efficacy of the absorber increases much faster at low efficiencies: at 25% efficiency, the abatement is already about half that obtained with the 90% efficiency. These first results, suggest strategies for the real-world application of a network of absorbers, and motivates further theoretical study to better characterize the details of the air flow inside and around the absorber.


title: Enhanced superconductivity and superconductor to insulator transition in nano-crystalline molybdenum thin films

abstract: Disorder driven superconductor to insulator transition via intermediate metallic regime is reported in nano-crystalline thin films of molybdenum. The nano-structured thin films have been deposited at room temperature using DC magnetron sputtering at different argon pressures. The grain size has been tuned using deposition pressure as the sole control parameter. A variation of particle sizes, room temperature resistivity and superconducting transition has been studied as a function of deposition pressure. The nano-crystalline molybdenum thin films are found to have large carrier concentration but very low mobility and electronic mean free path. Hall and conductivity measurements have been used to understand the effect of disorder on the carrier density and mobilities. Ioffe-Regel parameter is shown to correlate with the continuous metal-insulator transition in our samples.


title: The shear-free condition and constant-mean-curvature hyperboloidal initial data

abstract: We consider the Einstein-Maxwell-fluid constraint equations, and make use of the conformal method to construct and parametrize constant-mean-curvature hyperboloidal initial data sets that satisfy the shear-free condition. This condition is known to be necessary in order that a spacetime development admit a regular conformal boundary at future null infinity. We work with initial data sets in a variety of regularity classes, primarily considering those data sets whose geometries are weakly asymptotically hyperbolic, as defined in [arXiv:1506.03399]. These metrics are $C^{1,1}$ conformally compact, but not necessarily $C^2$ conformally compact. In order to ensure that the data sets we construct are indeed shear-free, we make use of the conformally covariant traceless Hessian introduced in [arXiv:1506.03399]. We furthermore construct a class of initial data sets with weakly asymptotically hyerbolic metrics that may be only $C^{0,1}$ conformally compact; these data sets are insufficiently regular to make sense of the shear-free condition.


title: Noise-induced synchronization in bidirectionally coupled type-I neurons

abstract: We present here some studies on noise-induced order and synchronous firing in a system of bidirectionally coupled generic type-I neurons. We find that transitions from unsynchronized to completely synchronized states occur beyond a critical value of noise strength that has a clear functional dependence on neuronal coupling strength and input values. For an inhibitory-excitatory (IE) synaptic coupling, the approach to a partially synchronized state is shown to vary qualitatively depending on whether the input is less or more than a critical value. We find that introduction of noise can cause a delay in the bifurcation of the firing pattern of the excitatory neuron for IE coupling.


title: Comments on the paper Growth and Characterization of Pure and Thiourea-Doped L-Alanine Single Crystals for NLO Devices

abstract: we show that the pure and thiourea-doped glycyl-L-alanine hydrochloride (GLAH and TU-GLAH) crystals claimed to have been grown by the authors of the title paper (Journal of Russian Laser Research, 34, 346-350 (2013)) are dubious crystals.


title: Fundamental representations of quantum affine superalgebras and R-matrices

abstract: We study a certain family of finite-dimensional simple representations over quantum affine superalgebras associated to general linear Lie superalgebras, the so-called fundamental representations: the denominators of rational $R$-matrices between two fundamental representations are computed; a cyclicity (and so simplicity) condition on tensor products of fundamental representations is proved.


title: Beltrami States in 2D Electron Magnetohydrodynamics

abstract: In this paper, the Hamiltonian formulations along with the Poisson brackets for two-dimensional (2D) electron magnetohydrodynamics (EMHD) flows are developed. These formulations are used to deduce the Beltrami states for 2D EMHD flows. In the massless electron limit, the EMHD Beltrami states reduce to the force-free states, though there is no force-free Beltrami state in the general EMHD case.


title: Gaia FGK Benchmark Stars: Effective temperatures and surface gravities

abstract: Large Galactic stellar surveys and new generations of stellar atmosphere models and spectral line formation computations need to be subjected to careful calibration and validation and to benchmark tests. We focus on cool stars and aim at establishing a sample of 34 Gaia FGK Benchmark Stars with a range of different metallicities. The goal was to determine the effective temperature and the surface gravity independently from spectroscopy and atmospheric models as far as possible. Fundamental determinations of Teff and logg were obtained in a systematic way from a compilation of angular diameter measurements and bolometric fluxes, and from a homogeneous mass determination based on stellar evolution models. The derived parameters were compared to recent spectroscopic and photometric determinations and to gravity estimates based on seismic data. Most of the adopted diameter measurements have formal uncertainties around 1%, which translate into uncertainties in effective temperature of 0.5%. The measurements of bolometric flux seem to be accurate to 5% or better, which contributes about 1% or less to the uncertainties in effective temperature. The comparisons of parameter determinations with the literature show in general good agreements with a few exceptions, most notably for the coolest stars and for metal-poor stars. The sample consists of 29 FGK-type stars and 5 M giants. Among the FGK stars, 21 have reliable parameters suitable for testing, validation, or calibration purposes. For four stars, future adjustments of the fundamental Teff are required, and for five stars the logg determination needs to be improved. Future extensions of the sample of Gaia FGK Benchmark Stars are required to fill gaps in parameter space, and we include a list of suggested candidates.


title: Graph-based compression of dynamic 3D point cloud sequences

abstract: This paper addresses the problem of compression of 3D point cloud sequences that are characterized by moving 3D positions and color attributes. As temporally successive point cloud frames are similar, motion estimation is key to effective compression of these sequences. It however remains a challenging problem as the point cloud frames have varying numbers of points without explicit correspondence information. We represent the time-varying geometry of these sequences with a set of graphs, and consider 3D positions and color attributes of the points clouds as signals on the vertices of the graphs. We then cast motion estimation as a feature matching problem between successive graphs. The motion is estimated on a sparse set of representative vertices using new spectral graph wavelet descriptors. A dense motion field is eventually interpolated by solving a graph-based regularization problem. The estimated motion is finally used for removing the temporal redundancy in the predictive coding of the 3D positions and the color characteristics of the point cloud sequences. Experimental results demonstrate that our method is able to accurately estimate the motion between consecutive frames. Moreover, motion estimation is shown to bring significant improvement in terms of the overall compression performance of the sequence. To the best of our knowledge, this is the first paper that exploits both the spatial correlation inside each frame (through the graph) and the temporal correlation between the frames (through the motion estimation) to compress the color and the geometry of 3D point cloud sequences in an efficient way.


title: Harbourne constants and conic configurations on the projective plane

abstract: In this note we exhibit the so-called Harbourne constants which capture and measure the Bounded Negativity on various birational models of an algebraic surface. We show an estimation for Harbourne constants for conic configurations on the complex projective plane with transversal intersection points.


title: The 12 prophets dataset

abstract: The "Ajeijadinho 3D" project is an initiative supported by the University of S\~ao Paulo (Museum of Science and Dean of Culture and Extension), which involves the 3D digitization of art works of Brazilian sculptor Antonio Francisco Lisboa, better known as Aleijadinho. The project made use of advanced acquisition and processing of 3D meshes for preservation and dissemination of the cultural heritage. The dissemination occurs through a Web portal, so that the population has the opportunity to meet the art works in detail using 3D visualization and interaction. The portal address is http://www.aleijadinho3d.icmc.usp.br. The 3D acquisitions were conducted over a week at the end of July 2013 in the cities of Ouro Preto, MG, Brazil and Congonhas do Campo, MG, Brazil. The scanning was done with a special equipment supplied by company Leica Geosystems, which allowed the work to take place at distances between 10 and 30 meters, defining a non-invasive procedure, simplified logistics, and without the need for preparation or isolation of the sites. In Ouro Preto, we digitized the churches of Francisco of Assis, Our Lady of Carmo, and Our Lady of Mercy; in Congonhas do Campo we scanned the entire Sanctuary of Bom Jesus de Matosinhos and his 12 prophets. Once scanned, the art works went through a long process of preparation, which required careful handling of meshes done by experts from the University of S\~ao Paulo in partnership with company Imprimate.


title: On enforcing maximum principles and achieving element-wise species balance for advection-diffusion-reaction equations under the finite element method

abstract: We present a robust computational framework for advective-diffusive-reactive systems that satisfies maximum principles, the non-negative constraint, and element-wise species balance property. The proposed methodology is valid on general computational grids, can handle heterogeneous anisotropic media, and provides accurate numerical solutions even for very high P\'eclet numbers. The significant contribution of this paper is to incorporate advection (which makes the spatial part of the differential operator non-self-adjoint) into the non-negative computational framework, and overcome numerical challenges associated with advection. We employ low-order mixed finite element formulations based on least-squares formalism, and enforce explicit constraints on the discrete problem to meet the desired properties. The resulting constrained discrete problem belongs to convex quadratic programming for which a unique solution exists. Maximum principles and the non-negative constraint give rise to bound constraints while element-wise species balance gives rise to equality constraints. The resulting convex quadratic programming problems are solved using an interior-point algorithm. Several numerical results pertaining to advection-dominated problems are presented to illustrate the robustness, convergence, and the overall performance of the proposed computational framework.


title: Approximate Inference with the Variational Holder Bound

abstract: We introduce the Variational Holder (VH) bound as an alternative to Variational Bayes (VB) for approximate Bayesian inference. Unlike VB which typically involves maximization of a non-convex lower bound with respect to the variational parameters, the VH bound involves minimization of a convex upper bound to the intractable integral with respect to the variational parameters. Minimization of the VH bound is a convex optimization problem; hence the VH method can be applied using off-the-shelf convex optimization algorithms and the approximation error of the VH bound can also be analyzed using tools from convex optimization literature. We present experiments on the task of integrating a truncated multivariate Gaussian distribution and compare our method to VB, EP and a state-of-the-art numerical integration method for this problem.


title: Robust Bayesian inference via coarsening

abstract: The standard approach to Bayesian inference is based on the assumption that the distribution of the data belongs to the chosen model class. However, even a small violation of this assumption can have a large impact on the outcome of a Bayesian procedure. We introduce a simple, coherent approach to Bayesian inference that improves robustness to perturbations from the model: rather than condition on the data exactly, one conditions on a neighborhood of the empirical distribution. When using neighborhoods based on relative entropy estimates, the resulting "coarsened" posterior can be approximated by simply tempering the likelihood---that is, by raising it to a fractional power---thus, inference is often easily implemented with standard methods, and one can even obtain analytical solutions when using conjugate priors. Some theoretical properties are derived, and we illustrate the approach with real and simulated data, using mixture models, autoregressive models of unknown order, and variable selection in linear regression.


title: GRINS: A Multiphysics Framework Based on the libMesh Finite Element Library

abstract: The progression of scientific computing resources has enabled the numerical approximation of mathematical models describing complex physical phenomena. A significant portion of researcher time is typically dedicated to the development of software to compute the numerical solutions. This work describes a flexible C++ software framework, built on the libMesh finite element library, designed to alleviate developer burden and provide easy access to modern computational algorithms, including quantity-of-interest-driven parallel adaptive mesh refinement on unstructured grids and adjoint-based sensitivities. Other software environments are highlighted and the current work motivated; in particular, the present work is an attempt to balance software infrastructure and user flexibility. The applicable class of problems and design of the software components is discussed in detail. Several examples demonstrate the effectiveness of the design, including applications that incorporate uncertainty. Current and planned developments are discussed.


title: Stationary properties of maximum entropy random walks

abstract: Maximum entropy (maxEnt) inference of state probabilities using state-dependent constraints is popular in the study of complex systems. In stochastic dynamical systems, the effect of state space topology and path-dependent constraints on the inferred state probabilities is unknown. To that end, we derive the transition probabilities and the stationary distribution of a maximum {\it path} entropy Markov process subject to state- and path-dependent constraints. The stationary distribution reflects a competition between path multiplicity and imposed constraints and is significantly different from the Boltzmann distribution. We illustrate our results with a particle diffusing on an energy landscape. Connections with the path integral approach to diffusion are discussed.


title: On the braking index of the unusual high-B rotation-powered pulsar PSR J1846-0258

abstract: PSR J1846-0258 is an object which straddles the boundary between magnetars and rotation powered pulsars. Though behaving for many years as a rotation-powered pulsar, in 2006, it exhibited distinctly magnetar-like behavior - emitting several short hard X-ray bursts, and a flux increase. Here we report on 7 years of post-outburst timing observations of PSR J1846-0258 using the Rossi X-ray Timing Explorer and the Swift X-ray Telescope. We measure the braking index over the post-magnetar outburst period to be $n=2.19\pm0.03$. This represents a change of $\Delta n=-0.46\pm0.03$ or a 14.5$\;\sigma$ difference from the pre-outburst braking index of $n=2.65\pm0.01$, which itself was measured over a span of 6.5 yr. So large and long-lived a change to a pulsar braking index is unprecedented and poses a significant challenge to models of pulsar spin-down.


title: Study of the 5p3/2 -> 6p3/2 electric dipole forbidden transition in atomic rubidium using optical-optical double resonance spectroscopy

abstract: Direct evidence of excitation of the 5p3/2 -> 6p3/2 electric dipole forbidden transition in atomic rubidium is presented. The experiments were performed in a room temperature rubidium cell with continuous wave extended cavity diode lasers. Optical-optical double resonance spectroscopy with counterpropagating beams allows the detection of the non-dipole transition free of Doppler broadening. The 5p3/2 state is prepared by excitation with a laser locked to the maximum F cyclic transition of the D2 line, and the forbidden transition is produced by excitation with a 911 nm laser. Production of the forbidden transition is monitored by detection of the 420 nm fluorescence that results from decay of the 6p3/2 state. Spectra with three narrow lines (~ 13 MHz FWHM) with the characteristic F - 1, F and F + 1 splitting of the 6p3/2 hyperfine structure in both rubidium isotopes were obtained. The results are in very good agreement with a direct calculation that takes into account the 5s -> 5p3/2 preparation dynamics, the 5p3/2 -> 6p3/2 non-dipole excitation geometry and the 6p3/2 -> 5s1/2 decay. The comparison also shows that the electric dipole forbidden transition is a very sensitive probe of the preparation dynamics.


title: A New Phenomenon: Sub-Tg, Solid-State, Plasticity-Induced Bonding in Polymers

abstract: Polymer self-adhesion due to the interdiffusion of macromolecules has been an active area of research for several decades [70, 43, 62, 42, 72, 73, 41]. Here, we report a new phenomenon of sub-Tg, solid-state, plasticity-induced bonding; where amorphous polymeric films were bonded together in a period of time on the order of a second in the solid-state at ambient temperatures nearly 60 K below their glass transition temperature (Tg) by subjecting them to active plastic deformation. Despite the glassy regime, the bulk plastic deformation triggered the requisite molecular mobility of the polymer chains, causing interpenetration across the interfaces held in contact. Quantitative levels of adhesion and the morphologies of the fractured interfaces validated the sub-Tg, plasticity-induced, molecular mobilization causing bonding. No-bonding outcomes (i) during the compression of films in a near hydrostatic setting (which inhibited plastic flow) and (ii) between an 'elastic' and a 'plastic' film further established the explicit role of plastic deformation in this newly reported sub-Tg solid-state bonding .


title: The computational complexity of calculating partition functions of optimal medians with Hamming distance

abstract: In this paper, we show that calculating the partition function of optimal medians of binary strings with Hamming distance is \#P-complete for several weight functions. The case when the weight function is the factorial function has application in bioinformatics. In that case, the partition function counts the most parsimonious evolutionary scenarios on a star tree under several models in bioinformatics. The results are extended to binary trees and we show that it is also \#P-complete to calculate the most parsimonious evolutionary scenarios on an arbitrary binary tree under the substitution model of biological sequences and under the Single Cut-or-Join model for genome rearrangements.


title: Entangling superconducting qubits in a multi-cavity system

abstract: Important tasks in cavity quantum electrodynamics include the generation and control of quantum states of spatially-separated particles distributed in different cavities. An interesting question in this context is how to prepare entanglement among particles located in different cavities, which are important for large-scale quantum information processing. We here consider a multi-cavity system where cavities are coupled to a superconducting (SC) qubit and each cavity hosts many SC qubits. We show that all intra-cavity SC qubits plus the coupler SC qubit can be prepared in an entangled Greenberger-Horne-Zeilinger (GHZ) state, by using a single operation and without the need of measurements. The GHZ state is created without exciting the cavity modes; thus greatly suppressing the decoherence caused by the cavity-photon decay and the effect of unwanted inter-cavity crosstalk on the operation. We also introduce two simple methods for entangling the intra-cavity SC qubits in a GHZ state. As an example, our numerical simulations show that it is feasible, with current circuit-QED technology, to prepare high-fidelity GHZ states, for up to nine SC qubits by using SC qubits distributed in two cavities. This proposal can in principle be used to implement a GHZ state for {\it an arbitrary number} of SC qubits distributed in multiple cavities. The proposal is quite general and can be applied to a wide range of physical systems, with the intra-cavity qubits being either atoms, NV centers, quantum dots, or various SC qubits.


title: Effective Hamiltonian for a Half-filled Asymmetric Ionic Hubbard Chain with Alternating On-site Interaction

abstract: We derive an effective spin Hamiltonian for the one-dimensional half-filled asymmetric ionic Hubbard model with alternating on-site interaction in the limit of strong repulsion. It is shown that the effective Hamiltonian is that of a spin $S=1/2$ anisotropic $XXZ$ Heisenberg chain with alternating next-nearest-neighbor and three-spin couplings in the presence of a uniform and a staggered magnetic field.


title: Simplicity and Finiteness of Discrete Spectrum of the Benjamin-Ono Scattering Operator

abstract: A spectral analysis is done on the $L$ operator of the Lax pair for the Benjamin-Ono equation. Simplicity and finiteness of the discrete spectrum are established as are needed for the Fokas and Ablowitz inverse scattering transform scheme. A crucial step in the simplicity proof is the discovery of a new identity connecting the $L^2$ norm of the eigenvector to its inner product with the scattering potential. The proof for finiteness is an extension of the ideas involved in the Birman-Schwinger bound for Schr\"odinger operators.


title: Edge states in honeycomb structures

abstract: An edge state is a time-harmonic solution of a conservative wave system, e.g. Schroedinger, Maxwell, which is propagating (plane-wave-like) parallel to, and localized transverse to, a line-defect or "edge". Topologically protected edge states are edge states which are stable against spatially localized (even strong) deformations of the edge. First studied in the context of the quantum Hall effect, protected edge states have attracted huge interest due to their role in the field of topological insulators. Theoretical understanding of topological protection has mainly come from discrete (tight-binding) models and direct numerical simulation. In this paper we consider a rich family of continuum PDE models for which we rigorously study regimes where topologically protected edge states exist. Our model is a class of Schroedinger operators on $\mathbb{R}^2$ with a background 2D honeycomb potential perturbed by an "edge-potential". The edge potential is a domain-wall interpolation, transverse to a prescribed "rational" edge, between two distinct periodic structures. General conditions are given for the bifurcation of a branch of topologically protected edge states from Dirac points of the background honeycomb structure. The bifurcation is seeded by the zero mode of a 1D effective Dirac operator. A key condition is a spectral no-fold condition for the prescribed edge. We then use this result to prove the existence of topologically protected edge states along zigzag edges of certain honeycomb structures. Our results are consistent with the physics literature and appear to be the first rigorous results on the existence of topologically protected edge states for continuum 2D PDE systems describing waves in a non-trivial periodic medium. We also show that the family of Hamiltonians we study contains cases where zigzag edge states exist, but which are not topologically protected.

categories: math-ph cond-mat.mes-hall math.AP math.MP quant-ph


title: The Extreme Value Machine

abstract: It is often desirable to be able to recognize when inputs to a recognition function learned in a supervised manner correspond to classes unseen at training time. With this ability, new class labels could be assigned to these inputs by a human operator, allowing them to be incorporated into the recognition function --- ideally under an efficient incremental update mechanism. While good algorithms that assume inputs from a fixed set of classes exist, e.g., artificial neural networks and kernel machines, it is not immediately obvious how to extend them to perform incremental learning in the presence of unknown query classes. Existing algorithms take little to no distributional information into account when learning recognition functions and lack a strong theoretical foundation. We address this gap by formulating a novel, theoretically sound classifier --- the Extreme Value Machine (EVM). The EVM has a well-grounded interpretation derived from statistical Extreme Value Theory (EVT), and is the first classifier to be able to perform nonlinear kernel-free variable bandwidth incremental learning. Compared to other classifiers in the same deep network derived feature space, the EVM is accurate and efficient on an established benchmark partition of the ImageNet dataset.


title: Syntactic categories for Nori motives

abstract: We give a new construction, based on categorical logic, of Nori's $\mathbb Q$-linear abelian category of mixed motives associated to a cohomology or homology functor with values in finite-dimensional vector spaces over $\mathbb Q$. This new construction makes sense for infinite-dimensional vector spaces as well, so that it associates a $\mathbb Q$-linear abelian category of mixed motives to any (co)homology functor, not only Betti homology (as Nori had done) but also, for instance, $\ell$-adic, $p$-adic or motivic cohomology. We prove that the $\mathbb Q$-linear abelian categories of mixed motives associated to different (co)homology functors are equivalent if and only a family (of logical nature) of explicit properties is shared by these different functors. The problem of the existence of a universal cohomology theory and of the equivalence of the information encoded by the different classical cohomology functors thus reduces to that of checking these explicit conditions.


title: Secure Degrees of Freedom of One-hop Wireless Networks with No Eavesdropper CSIT

abstract: We consider three channel models: the wiretap channel with $M$ helpers, the $K$-user multiple access wiretap channel, and the $K$-user interference channel with an external eavesdropper, when no eavesdropper's channel state information (CSI) is available at the transmitters. In each case, we establish the optimal sum secure degrees of freedom (s.d.o.f.) by providing achievable schemes and matching converses. We show that the unavailability of the eavesdropper's CSIT does not reduce the s.d.o.f. of the wiretap channel with helpers. However, there is loss in s.d.o.f. for both the multiple access wiretap channel and the interference channel with an external eavesdropper. In particular, we show that in the absence of eavesdropper's CSIT, the $K$-user multiple access wiretap channel reduces to a wiretap channel with $(K-1)$ helpers from a sum s.d.o.f. perspective, and the optimal sum s.d.o.f. reduces from $\frac{K(K-1)}{K(K-1)+1}$ to $\frac{K-1}{K}$. For the interference channel with an external eavesdropper, the optimal sum s.d.o.f. decreases from $\frac{K(K-1)}{2K-1}$ to $\frac{K-1}{2}$ in the absence of the eavesdropper's CSIT. Our results show that the lack of eavesdropper's CSIT does not have a significant impact on the optimal s.d.o.f. for any of the three channel models, especially when the number of users is large. This implies that physical layer security can be made robust to the unavailability of eavesdropper CSIT at high signal to noise ratio (SNR) regimes by careful modification of the achievable schemes as demonstrated in this paper.


title: On Liftings of Projective Indecomposable $G_{(1)}$-Modules

abstract: Let $G$ be a simple simply connected algebraic group over an algebraically closed field $k$ of characteristic $p$, with Frobenius kernel $G_{(1)}$. It is known that when $p\ge 2h-2$, where $h$ is the Coxeter number of $G$, the projective indecomposable $G_{(1)}$-modules (PIMs) lift to $G$, and this has been conjectured to hold in all characteristics. In this paper, we explore the lifting problem via extensions of algebraic groups, following the work of Parshall and Scott who in turn build upon ideas due to Donkin. We prove various results which augment this approach, and as an application demonstrate that the PIMs lift to $G_{(1)}H$, for particular closed subgroups $H \le G$ which contain a maximal torus of $G$.


title: Optomechanical Bell test

abstract: Photons of a laser beam driving the upper motional sideband of an optomechanical cavity can decay into photon-phonon pairs by means of an optomechanical parametric process. The phononic state can subsequently be mapped to a photonic state by exciting the lower sideband, hence creating photon-photon pairs out of an optomechanical system. We here demonstrate that these pairs can violate a Bell inequality when they are measured with photon counting techniques preceded by small displacement operations in the phase space. Our results show how to detect non-local light-mechanical correlations and open the way for device-independent quantum optomechanics where non-classical light-mechanical correlations can be certified without assumptions on the state structure or on the measurement devices.


title: Convergence of Sequential Quasi-Monte Carlo Smoothing Algorithms

abstract: Gerber and Chopin (2015) recently introduced Sequential quasi-Monte Carlo (SQMC) algorithms as an efficient way to perform filtering in state-space models. The basic idea is to replace random variables with low-discrepancy point sets, so as to obtain faster convergence than with standard particle filtering. Gerber and Chopin (2015) describe briefly several ways to extend SQMC to smoothing, but do not provide supporting theory for this extension. We discuss more thoroughly how smoothing may be performed within SQMC, and derive convergence results for the so-obtained smoothing algorithms. We consider in particular SQMC equivalents of forward smoothing and forward filtering backward sampling, which are the most well-known smoothing techniques. As a preliminary step, we provide a generalization of the classical result of Hlawka and M\"uck (1972) on the transformation of QMC point sets into low discrepancy point sets with respect to non uniform distributions. As a corollary of the latter, we note that we can slightly weaken the assumptions to prove the consistency of SQMC.


title: Improved Sampling Algorithms in Lattice QCD

abstract: Reverse Monte Carlo (RMC) is an algorithm that incorporates stochastic modification of the action as part of the process that updates the fields in a Monte Carlo simulation. Such update moves have the potential of lowering or eliminating potential barriers that may cause inefficiencies in exploring the field configuration space. The highly successful Cluster algorithms for spin systems can be derived from the RMC framework. In this work we apply RMC ideas to pure gauge theory, aiming to alleviate the critical slowing down observed in the topological charge evolution as well as other long distance observables. We present various formulations of the basic idea and report on our numerical experiments with these algorithms.


title: Adsorption of soft particles at fluid interfaces

abstract: Soft particles can be better emulsifiers than hard particles because they stretch at fluid interfaces. This deformation can increase adsorption energies by orders of magnitude relative to rigid particles. The deformation of a particle at an interface is governed by a competition of bulk elasticity and surface tension. When particles are partially wet by the two liquids, deformation is localized within a material-dependent distance $L$ from the contact line. At the contact line, the particle morphology is given by a balance of surface tensions. When the particle radius $R \ll L$, the particle adopts a lenticular shape identical to that of an adsorbed fluid droplet. Particle deformations can be elastic or plastic, depending on the relative values of the Young modulus, $E$, and yield stress, $\sigma_p$. When surface tensions favour complete spreading of the particles at the interface, plastic deformation can lead to unusual fried-egg morphologies. When deformable particles have surface properties that are very similar to one liquid phase, adsorption can be extremely sensitive to small changes of their affinity for the other liquid phase. These findings have implications for the adsorption of microgel particles at fluid interfaces and the performance of stimuli-responsive Pickering emulsions.


title: Hadamard well-posedness of the gravity water waves system

abstract: We consider in this article the system of (pure) gravity water waves in any dimension and in fluid domains with general bottoms. The unique solvability of the problem was established by Alazard-Burq-Zuily [Invent. Math, 198 (2014), no. 1, 71--163] at a low regularity level where the initial surface is $C^{\frac{3}{2}+}$ in terms of Sobolev embeddings, which allows the existence of free surfaces with unbounded curvature. Our result states that the solutions obtained above depend continuously on initial data in the strong topology where they are constructed. This completes a well-posedness result in the sense of Hadamard.


title: Invariant Theory of Artin-Schelter Regular Algebras: A survey

abstract: This is a survey of results that extend notions of the classical invariant theory of linear actions by finite groups on $k[x_1, \dots, x_n]$ to the setting of finite group or Hopf algebra $H$ actions on an Artin-Schelter regular algebra $A$. We investigate when $A^H$ is AS regular, or AS Gorenstein, or a "complete intersection" in a sense that is defined. Directions of related research are explored briefly.


title: A compactness result for energy-minimizing harmonic maps with rough domain metric

abstract: In 1996, Shi generalized the epsilon-regularity theorem of Schoen and Uhlenbeck to energy-minimizing harmonic maps from a domain equipped with a bounded measurable Riemannian metric. In the present work we prove a compactness result for such energy-minimizing maps. As an application, we combine our result with Shi's theorem to give an improved bound on the Hausdorff dimension of the singular set, assuming that the map has bounded energy at all scales. This last assumption can be removed when the target manifold is simply-connected.


title: A Tracing of the Fractional Temperature Field

abstract: This note is devoted to a study of $L^q$-tracing of the fractional temperature field $u(t,x)$ -- the weak solution of the fractional heat equation $(\partial_t+(-\Delta_x)^\alpha)u(t,x)=g(t,x)$ in $L^p(\mathbb R^{1+n}_+)$ subject to the initial temperature $u(0,x)=f(x)$ in $L^p(\mathbb R^n)$.


title: Power Partitions

abstract: In 1918, Hardy and Ramanujan published a seminal paper which included an asymptotic formula for the partition function. In their paper, they also claim without proof an asymptotic equivalence for $p^k(n)$, the number of partitions of a number $n$ into $k$-th powers. In this paper, we provide an asymptotic formula for $p^k(n)$, using the Hardy-Littlewood Circle Method. We also provide a formula for the difference function $p^k(n+1)-p^k(n)$. As a necessary step in the proof, we obtain a non-trivial bound on exponential sums of the form $\sum_{m=1}^q e(\frac{am^k}{q})$.


title: Percolation on Hypergraphs with Four-Edges

abstract: We study percolation on self-dual hypergraphs that contain hyperedges with four bounding vertices, or "four-edges", using three different generators, each containing bonds or sites with three distinct probabilities $p$, $r$, and $t$ connecting the four vertices. We find explicit values of these probabilities that satisfy the self-duality conditions discussed by Bollob\'as and Riordan. This demonstrates that explicit solutions of the self-duality conditions can be found using generators containing bonds and sites with independent probabilities. These solutions also provide new examples of lattices where exact percolation critical points are known. One of the generators exhibits three distinct criticality solutions ($p$, $r$, $t$). We carry out Monte-Carlo simulations of two of the generators on two different hypergraphs to confirm the critical values. For the case of the hypergraph and uniform generator studied by Wierman et al., we also determine the threshold $p = 0.441374\pm 0.000001$, which falls within the tight bounds that they derived. Furthermore, we consider a generator in which all or none of the vertices can connect, and find a soluble inhomogeneous percolation system that interpolates between site percolation on the union-jack lattice and bond percolation on the square lattice.


title: Cusp and $b_1$ growth for ball quotients and maps onto $\mathbb{Z}$ with finitely generated kernel

abstract: Let $M = \mathbb{B}^2 / \Gamma$ be a smooth ball quotient of finite volume with first betti number $b_1(M)$ and let $\mathcal{E}(M) \ge 0$ be the number of cusps (i.e., topological ends) of $M$. We study the growth rates that are possible in towers of finite-sheeted coverings of $M$. In particular, $b_1$ and $\mathcal{E}$ have little to do with one another, in contrast with the well-understood cases of hyperbolic $2$- and $3$-manifolds. We also discuss growth of $b_1$ for congruence arithmetic lattices acting on $\mathbb{B}^2$ and $\mathbb{B}^3$. Along the way, we provide an explicit example of a lattice in $\mathrm{PU}(2, 1)$ admitting a homomorphism onto $\mathbb{Z}$ with finitely generated kernel. Moreover, we show that any cocompact arithmetic lattice $\Gamma \subset \mathrm{PU}(n, 1)$ of simplest type contains a finite index subgroup with this property.


title: Geodesics in the Engel group with a sub-Lorentzian metric

abstract: Let E be the Engel group and D be a rank 2 bracket generating left invariant distribution with a Lorentzian metric, which is a nondegenerate metric of index 1. In this paper, we first study some properties of horizontal curves on E. Second, we prove that time-like normal geodesics are locally maximizers in the Engel group, and calculate the explicit expression of non-space-like geodesics.


title: A thorny path of field theory: from triviality to interaction and confinement

abstract: Summation of the perturbation series for the Gell-Mann--Low function \beta(g) of \phi^4 theory leads to the asymptotics \beta(g)=\beta_\infty g^\alpha at g\to\infty, where \alpha\approx 1 for space dimensions d=2,3,4. The natural hypothesis arises, that asymptotic behavior is \beta(g) \sim g for all d. Consideration of the "toy" zero-dimensional model confirms the hypothesis and reveals the origin of this result: it is related with a zero of a certain functional integral. This mechanism remains valid for arbitrary space dimensionality d. The same result for the asymptotics is obtained for explicitly accepted lattice regularization, while the use of high-temperature expansions allows to calculate the whole \beta-function. As a result, the \beta-function of four-dimensional \phi^4 theory is appeared to be non-alternating and has a linear asymptotics at infinity. The analogous situation is valid for QED. According to the Bogoliubov and Shirkov classification, it means possibility to construct the continuous theory with finite interaction at large distances. This conclusion is in visible contradiction with the lattice results indicating triviality of \phi^4 theory. This contradiction is resolved by a special character of renormalizability in \phi^4 theory: to obtain the continuous renormalized theory, there is no need to eliminate a lattice from the bare theory. In fact, such kind of renormalizability is not accidental and can be understood in the framework of Wilson's many-parameter renormalization group. Application of these ideas to QCD shows that Wilson's theory of confinement is not purely illustrative, but has a direct relation to a real situation. As a result, the problem of analytical proof of confinement and a mass gap can be considered as solved, at least on the physical level of rigor.

categories: hep-ph cond-mat.stat-mech hep-lat hep-th


title: A simple application of FIC to model selection

abstract: We have recently proposed a new information-based approach to model selection, the Frequentist Information Criterion (FIC), that reconciles information-based and frequentist inference. The purpose of this current paper is to provide a simple example of the application of this criterion and a demonstration of the natural emergence of model complexities with both AIC-like ($N^0$) and BIC-like ($\log N$) scaling with observation number $N$. The application developed is deliberately simplified to make the analysis analytically tractable.


title: Higgs inflation in Gauss-Bonnet braneworld

abstract: The measured masses of the Higgs boson and top quark indicate that the effective potential of the standard model either develops an unstable electroweak vacuum or stands stable all the way up to the Planck scale. In the latter case in which the top quark mass is about $2\sigma$ below its present central value, the Higgs boson can be the inflaton with the help of a large nonminimal coupling to curvature in four dimensions. We propose a scenario in which the Higgs boson can be the inflaton in a five-dimensional Gauss-Bonnet braneworld model to solve both the unitarity and stability problems which usually plague Higgs inflation. We find that in order for Higgs inflation to happen successfully in the Gauss-Bonnet regime, the extra dimension scale must appear roughly in the range between the TeV scale and the instability scale of standard model. At the tree level, our model can give rise to a naturally small nonminimal coupling $\xi\sim\mathcal{O}(1)$ for the Higgs quartic coupling $\lambda\sim\mathcal{O}(0.1)$ if the extra dimension scale lies at the TeV scale. At the loop level, the inflationary predictions at the tree level are preserved. Our model can be confronted with future experiments and observations from both particle physics and cosmology.


title: Baryon Number as the Fourth Color

abstract: We propose an extension of the Standard Model in which baryon number is promoted to be part of a non-Abelian gauge symmetry at high energies. Specifically, we consider the gauge group SU(4) x SU(2)_L x U(1)_X, where the SU(4) unifies baryon number and color. This symmetry is spontaneously broken down to the Standard Model gauge group at a scale which can be as low as a few TeV. The SU(4) structure implies that each SM quark comes along with an uncolored quark partner, the lightest of which is stabilized by the generalized baryon number symmetry and can play the role of dark matter. We explore circumstances under which one can realize a model of asymmetric dark matter whose relic abundance is connected to the observed baryon asymmetry, and discuss unique signatures that can be searched for at the LHC.


title: Brief Follow-up on Recent Studies of Theia's Accretion

abstract: Kaib & Cowan (2015) recently used terrestrial planet formation simulations to conclude that the moon-forming impactor (Theia) had only a ~5% or less chance of having the same oxygen isotope composition as Earth, while Mastrobuono-Battisti et al. (2015) used seemingly similar simulations and methods to arrive at a higher value of ~20% or more. Here we derive the results of both papers from a single set of simulations. Compared to Kaib & Cowan (2015), the analysis of Mastrobuono-Battisti et al. (2015) systematically yields more massive Theia analogs and imposes flatter isotopic gradients across the original protoplanetary disk. Both of these effects diminish isotopic differences between Earth and Theia analogs. While it is notoriously difficult to produce systems resembling our actual terrestrial planets, the analysis of Kaib & Cowan (2015) more often selects and analyzes Earth and Mars analogs at orbital locations near the real planets. Given this, we conclude that the greater isotopic differences between Earth and Theia found in Kaib & Cowan (2015) better reflect the predictions of terrestrial planet formation models. Finally, although simulation uncertainties and a terrestrial contribution to Moon formation enhance the fraction of Theia analogs consistent with the canonical giant impact hypothesis, this fraction still remains in the 5-8% range.


title: Probing flavor models with Ge-76-based experiments on neutrinoless double-beta decay

abstract: The physics impact of a staged approach for double-beta decay experiments based on Ge-76 is studied. The scenario considered relies on realistic time schedules envisioned by the GERDA and the MAJORANA collaborations, which are jointly working towards the realization of a future larger scale Ge-76 experiment. Intermediate stages of the experiments are conceived to perform quasi background-free measurements, and different data sets can be reliably combined to maximize the physics outcome. The sensitivity for such a global analysis is presented, with focus on how neutrino flavor models can be probed already with preliminary phases of the experiments. The synergy between theory and experiment yields strong benefits for both sides: the model predictions can be used to sensibly plan the experimental stages, and results from intermediate stages can be used to constrain whole groups of theoretical scenarios. This strategy clearly generates added value to the experimental efforts, while at the same time it allows to achieve valuable physics results as early as possible.


title: Amplitudes for massive vector and scalar bosons in spontaneously-broken gauge theory from the CHY representation

abstract: In the formulation of Cachazo, He, and Yuan, tree-level amplitudes for massless particles in gauge theory and gravity can be expressed as rational functions of the Lorentz invariants $k_a \cdot k_b$, $\epsilon_a \cdot k_b$, and $\epsilon_a \cdot \epsilon_b$, valid in any number of spacetime dimensions. We use dimensional reduction of higher-dimensional amplitudes of particles with internal momentum $\kappa$ to obtain amplitudes for massive particles in lower dimensions. In the case of gauge theory, we argue that these massive amplitudes belong to a theory in which the gauge symmetry is spontaneously broken by an adjoint Higgs field. Consequently, we show that tree-level $n$-point amplitudes containing massive vector and scalar bosons in this theory can be obtained by simply replacing $k_a \cdot k_b$ with $k_a \cdot k_b - \kappa_a \kappa_b $ in the corresponding massless amplitudes, where the masses of the particles are given by $|\kappa_a|$.


title: Evidence of Halo Assembly Bias in Massive Clusters

abstract: We present significant evidence of halo assembly bias for SDSS redMaPPer galaxy clusters in the redshift range $[0.1, 0.33]$. By dividing the 8,648 clusters into two subsamples based on the average member galaxy separation from the cluster center, we first show that the two subsamples have very similar halo mass of $M_{\rm 200m}\simeq 1.9\times 10^{14}~h^{-1}M_\odot$ based on the weak lensing signals at small radii $R<\sim 10~h^{-1}{\rm Mpc}$. However, their halo bias inferred from both the large-scale weak lensing and the projected auto-correlation functions differs by a factor of $\sim$1.5, which is a signature of assembly bias. The same bias hypothesis for the two subsamples is excluded at 2.5$\sigma$ in the weak lensing and 4.4$\sigma$ in the auto-correlation data, respectively. This result could bring a significant impact on both galaxy evolution and precision cosmology.


title: Real space Eliashberg approach to charge order of electrons coupled to dynamic antiferromagnetic fluctuations

abstract: We study charge ordered solutions for fermions on a square lattice interacting with dynamic antiferromagnetic fluctuations. Our approach is based on real space Eliashberg equations which are solved self-consistently. We first show that the antiferromagnetic fluctuations can induce arc features in the spectral functions, as spectral weight is suppressed at the hot spots; however, no real pseudogap is generated. At low temperature spontaneous charge order with a $d$-form factor can be stabilized for certain parameters. As long as the interacting Fermi surfaces possesses hot spots, the ordering wave vector corresponds to the diagonal connection of the hot spots, similar to the non-self-consistent case. Tendencies towards observed axial order only appear in situations without hot spots.


title: Integration Rules for Scattering Equations

abstract: As described by Cachazo, He and Yuan, scattering amplitudes in many quantum field theories can be represented as integrals that are fully localized on solutions to the so-called scattering equations. Because the number of solutions to the scattering equations grows quite rapidly, the contour of integration involves contributions from many isolated components. In this paper, we provide a simple, combinatorial rule that immediately provides the result of integration against the scattering equation constraints for any M\"obius-invariant integrand involving only simple poles. These rules have a simple diagrammatic interpretation that makes the evaluation of any such integrand immediate. Finally, we explain how these rules are related to the computation of amplitudes in the field theory limit of string theory.


title: The evolution of lossy compression

abstract: In complex environments, there are costs to both ignorance and perception. An organism needs to track fitness-relevant information about its world, but the more information it tracks, the more resources it must devote to memory and processing. Rate-distortion theory shows that, when errors are allowed, remarkably efficient internal representations can be found by biologically-plausible hill-climbing mechanisms. We identify two regimes: a high-fidelity regime where perceptual costs scale logarithmically with environmental complexity, and a low-fidelity regime where perceptual costs are, remarkably, independent of the environment. When environmental complexity is rising, Darwinian evolution should drive organisms to the threshold between the high- and low-fidelity regimes. Organisms that code efficiently will find themselves able to make, just barely, the most subtle distinctions in their environment.

categories: q-bio.NC cs.IT math.IT nlin.AO physics.soc-ph q-bio.PE


title: Renormalisation Group Corrections to Neutrino Mass Sum Rules

abstract: Neutrino mass sum rules are an important class of predictions in flavour models relating the Majorana phases to the neutrino masses. This leads, for instance, to enormous restrictions on the effective mass as probed in experiments on neutrinoless double beta decay. While up to now these sum rules have in practically all cases been taken to hold exactly, we will go here beyond that. After a discussion of the types of corrections that could possibly appear and elucidating on the theory behind neutrino mass sum rules, we estimate and explicitly compute the impact of radiative corrections, as these appear in general and thus hold for whole groups of models. We discuss all neutrino mass sum rules currently present in the literature, which together have realisations in more than 50 explicit neutrino flavour models. We find that, while the effect of the renormalisation group running can be visible, the qualitative features do not change. This finding strongly backs up the solidity of the predictions derived in the literature, and it thus marks a very important step in deriving testable and reliable predictions from neutrino flavour models.


title: Astrophysical hydrodynamics with a high-order discontinuous Galerkin scheme and adaptive mesh refinement

abstract: Solving the Euler equations of ideal hydrodynamics as accurately and efficiently as possible is a key requirement in many astrophysical simulations. It is therefore important to continuously advance the numerical methods implemented in current astrophysical codes, especially also in light of evolving computer technology, which favours certain computational approaches over others. Here we introduce the new adaptive mesh refinement (AMR) code TENET, which employs a high order discontinuous Galerkin (DG) scheme for hydrodynamics. The Euler equations in this method are solved in a weak formulation with a polynomial basis by means of explicit Runge-Kutta time integration and Gauss-Legendre quadrature. This approach offers significant advantages over commonly employed second order finite volume (FV) solvers. In particular, the higher order capability renders it computationally more efficient, in the sense that the same precision can be obtained at significantly less computational cost. Also, the DG scheme inherently conserves angular momentum in regions where no limiting takes place, and it typically produces much smaller numerical diffusion and advection errors than a FV approach. A further advantage lies in a more natural handling of AMR refinement boundaries, where a fall-back to first order can be avoided. Finally, DG requires no wide stencils at high order, and offers an improved data locality and a focus on local computations, which is favourable for current and upcoming highly parallel supercomputers. We describe the formulation and implementation details of our new code, and demonstrate its performance and accuracy with a set of two- and three-dimensional test problems. The results confirm that DG schemes have a high potential for astrophysical applications.


title: Discovering Uncolored Naturalness in Exotic Higgs Decays

abstract: Solutions to the hierarchy problem usually require top partners. In standard SUSY or composite Higgs theories, the partners carry SM color and are becoming increasingly constrained by LHC searches. However, theories like Folded SUSY (FS), Twin Higgs (TH) and Quirky Little Higgs (QLH) introduce uncolored top partners, which can be SM singlets or carry electroweak charge. Their small production cross section left doubt as to whether the LHC can effectively probe such scenarios. Typically, these partners are charged under their own mirror color gauge group. In FS and QLH, the absence of light mirror matter allows glueballs to form at the bottom of the mirror spectrum. This is also the case in some TH realizations. The Higgs can decay to these mirror glueballs, with the glueballs decaying into SM particles with potentially observable lifetimes. We undertake the first detailed study of this glueball signature and quantitatively demonstrate the discovery potential of uncolored naturalness via exotic Higgs decays at the LHC and a potential future 100 TeV collider. Our findings indicate that mirror glueballs are the smoking gun signature of natural FS and QLH type theories, in analogy to tree-level Higgs coupling shifts for the TH. We show that glueball masses in the $\sim$ 10-60 GeV mass range are theoretically preferred. Careful treatment of lifetime, mirror-hadronization and nonperturbative uncertainties is required to perform meaningful collider studies. We outline several new search strategies for exotic Higgs decays of the form $h\to XX \to 4f$ at the LHC, with $X$ having lifetimes in the $10 \mu m$ to $km$ range. We find that FS stops can be probed with masses up to 600 (1100) GeV at the LHC with 300 (3000) /fb of data, and TH top partners could be accessible with masses up to 900 (1500) GeV. This makes exotic Higgs decays the prime discovery channel for uncolored naturalness at the LHC.


title: Lord of the Rings: A Kinematic Distance to Circinus X-1 from a Giant X-Ray Light Echo

abstract: Circinus X-1 exhibited a bright X-ray flare in late 2013. Follow-up observations with Chandra and XMM-Newton from 40 to 80 days after the flare reveal a bright X-ray light echo in the form of four well-defined rings with radii from 5 to 13 arcminutes, growing in radius with time. The large fluence of the flare and the large column density of interstellar dust towards Circinus X-1 make this the largest and brightest set of rings from an X-ray light echo observed to date. By deconvolving the radial intensity profile of the echo with the MAXI X-ray lightcurve of the flare we reconstruct the dust distribution towards Circinus X-1 into four distinct dust concentrations. By comparing the peak in scattering intensity with the peak intensity in CO maps of molecular clouds from the Mopra Southern Galactic Plane CO Survey we identify the two innermost rings with clouds at radial velocity ~ -74 km/s and ~ -81 km/s, respectively. We identify a prominent band of foreground photoelectric absorption with a lane of CO gas at ~ -32 km/s. From the association of the rings with individual CO clouds we determine the kinematic distance to Circinus X-1 to be $D_{Cir X-1} = 9.4^{+0.8}_{-1.0}$ kpc. This distance rules out earlier claims of a distance around 4 kpc, implies that Circinus X-1 is a frequent super-Eddington source, and places a lower limit of $\Gamma \gtrsim 22$ on the Lorentz factor and an upper limit of $\theta_{jet} \lesssim 3^{\circ}$ on the jet viewing angle.


title: KAT-7 Science Verification: Cold Gas, Star Formation, and Substructure in the Nearby Antlia Cluster

abstract: The Antlia Cluster is a nearby, dynamically young structure, and its proximity provides a valuable opportunity for detailed study of galaxy and group accretion onto clusters. We present a deep HI mosaic completed as part of spectral line commissioning of the Karoo Array Telescope (KAT-7), and identify infrared counterparts from the WISE extended source catalog to study neutral atomic gas content and star formation within the cluster. We detect 37 cluster members out to a radius of ~0.9 Mpc with M_HI > 5x10^7 M_Sun. Of these, 35 are new HI detections, 27 do not have previous spectroscopic redshift measurements, and one is the Compton thick Seyfert II, NGC 3281, which we detect in HI absorption. The HI galaxies lie beyond the X-ray emitting region 200 kpc from the cluster center and have experienced ram pressure stripping out to at least 600 kpc. At larger radii, they are distributed asymmetrically suggesting accretion from surrounding filaments. Combining HI with optical redshifts, we perform a detailed dynamical analysis of the internal substructure, identify large infalling groups, and present the first compilation of the large scale distribution of HI, and star forming galaxies within the cluster. We find that elliptical galaxy NGC 3268 is at the center of the oldest substructure and argue that NGC 3258 and its companion population are more recent arrivals. Through the presence of HI and on-going star formation, we rank substructures with respect to their relative time since accretion onto Antlia.


title: The radial profile and flattening of the Milky Way's stellar halo to $\rm 80~$kpc from the SEGUE K-giant Survey

abstract: We characterise the radial density, metallicity and flattening profile of the Milky Way's stellar halo, based on the large sample of 1757 spectroscopically confirmed giant stars from SDSS/SEGUE-2 after excising stars that were algorithmically attributed to apparent halo substructure (including the Sagittarius stream). Compared to BHB stars or RR Lyrae, giants are more readily understood tracers of the overall halo star population, with less bias in age or metallicity. The well-characterized selection function of the sample enables forward modelling of those data, based on ellipsoidal stellar density models, $\nu_* (R,z)$, with Einasto profiles and (broken) power laws for their radial dependence, combined with a model for the metallicity gradient and the flattening profile. Among models with constant flattening, these data are reasonably well fit by an Einasto profile of $n=3.1\pm 0.5$ with an effective radius $\rm r_{eff} = 15\pm2~$kpc and a flattening of $q=0.7\pm 0.02$; or comparably well by an equally flattened broken power-law, with radial slopes of $\alpha_{in}=2.1\pm 0.3$ and $\alpha_{out}=3.8\pm 0.1$, with a break-radius of $r_{break}=18\pm1$~kpc; this is largely consistent with earlier work. We find a modest, but significant metallicity gradient within the "outer" stellar halo, $\rm [Fe/H]$ decreasing outward. If we allow for a variable flattening $q = f(r_{GC} )$, we find the distribution of halo giants to be considerably more flattened at small radii, $q({\rm 10~kpc})\sim 0.57$, compared to $q(>30{\rm kpc})\sim 0.8$. Remarkably, the data are then very well fit by a single power-law of index $\rm \sim 4.2\pm0.1$ of the variable $r_q\equiv\sqrt{R^2+(z/q(r))^2}$. In this simple and better fitting model, there is a break in flattening at $\sim 20$~kpc, instead of a break in the radial density function.


title: Reliability of the one-crossing approximation in describing the Mott transition

abstract: We assess the reliability of the one-crossing approximation (OCA) approach in quantitative description of the Mott transition in the framework of the dynamical mean field theory (DMFT). The OCA approach has been applied in the conjunction with DMFT to a number of heavy-fermion, actinide, transition metal compounds, and nanoscale systems. However, several recent studies in the framework of impurity models pointed out to serious deficiencies of OCA and raised questions regarding its reliability. Here we consider a single band Hubbard model on the Bethe lattice at finite temperatures and compare the results of OCA to those of a numerically exact quantum Monte Carlo (QMC) method. The temperature-local repulsion U phase diagram for the particle-hole symmetric case obtained by OCA is in good agreement with that of QMC, with the metal-insulator transition captured very well. We find, however, that the insulator to metal transition is shifted to higher values of U and, simultaneously, correlations in the metallic phase are significantly overestimated. This counter-intuitive behavior is due to simultaneous underestimations of the Kondo scale in the metallic phase and the size of the insulating gap. We trace the underestimation of the insulating gap to that of the second moment of the high-frequency expansion of the impurity spectral density. Calculations for the system away from the particle-hole symmetric case are also presented and discussed.


title: "Scheme-theoretic images" of morphisms of stacks

abstract: We give criteria for certain morphisms from an algebraic stack to a (not necessarily algebraic) stack to admit an (appropriately defined) scheme-theoretic image. We apply our criteria to show that certain natural moduli stacks of local Galois representations are algebraic (or Ind-algebraic) stacks.


title: Depolarizing channel parameter estimation using noisy initial states

abstract: We consider estimating the parameter associated with the qubit depolarizing channel when the available initial states that might be employed are mixed. We use quantum Fisher information as a measure of the accuracy of estimation to compare protocols which use collections of qubits in product states to one in which the qubits are in a correlated state. We show that, for certain parameter values and initial states, the correlated state protocol can yield a greater accuracy per channel invocation than the product state protocols. We show that, for some parameters and initial states, using more than two qubits and channel invocations is advantageous. These results stand in contrast to the known optimal case that uses pure initial states and a single channel invocation on a pair of entangled qubits.


title: The square sieve and the large sieve with square moduli

abstract: We give a short alternative proof using Heath-Brown's square sieve of a bound of the author for the large sieve with square moduli.


title: Typologies of the Popular Science Web Video

abstract: The creation of popular science web videos on the Internet has increased in recent years. The diversity of formats, genres, and producers makes it difficult to formulate a universal definition of science web videos since not every producer considers him- or herself to be a science communicator in an institutional sense, and professionalism and success on video platforms no longer depend exclusively on technical excellence or production costs. Entertainment, content quality, and authenticity have become the keys to community building and success. The democratization of science video production allows a new variety of genres, styles, and forms. This article provides a first overview of the typologies and characteristics of popular science web videos. To avoid a misleading identification of science web videos with institutionally produced videos, we steer clear of the term science communication video, since many of the actual producers are not even familiar with the academic discussion on science communication, and since the subject matter does not depend on political or educational strategies. A content analysis of 200 videos from 100 online video channels was conducted. Several factors such as narrative strategies, video editing techniques, and design tendencies with regard to cinematography, the number of shots, the kind of montage used, and even the spread use of sound design and special FX point to an increasing professionalism among science communicators independent of institutional or personal commitments: in general, it can be said that supposed amateurs are creating the visual language of science video communication. This study represents an important step in understanding the essence of current popular science web videos and provides an evidence-based definition as a helpful cornerstone for further studies on science communication within this kind of new media.


title: Confirmation of Wide-Field Signatures in Redshifted 21 cm Power Spectra

abstract: We confirm our recent prediction of the "pitchfork" foreground signature in power spectra of high-redshift 21 cm measurements where the interferometer is sensitive to large-scale structure on all baselines. This is due to the inherent response of a wide-field instrument and is characterized by enhanced power from foreground emission in Fourier modes adjacent to those considered to be the most sensitive to the cosmological H I signal. In our recent paper, many signatures from the simulation that predicted this feature were validated against Murchison Widefield Array (MWA) data, but this key pitchfork signature was close to the noise level. In this paper, we improve the data sensitivity through the coherent averaging of 12 independent snapshots with identical instrument settings and provide the first confirmation of the prediction with a signal-to-noise ratio > 10. This wide-field effect can be mitigated by careful antenna designs that suppress sensitivity near the horizon. Simple models for antenna apertures that have been proposed for future instruments such as the Hydrogen Epoch of Reionization Array and the Square Kilometre Array indicate they should suppress foreground leakage from the pitchfork by ~40 dB relative to the MWA and significantly increase the likelihood of cosmological signal detection in these critical Fourier modes in the three-dimensional power spectrum.


title: The final-state problem for the cubic-quintic NLS with non-vanishing boundary conditions

abstract: We construct solutions with prescribed scattering state to the cubic-quintic NLS $$ (i\partial_t+\Delta)\psi=\alpha_1 \psi-\alpha_{3}\vert \psi\vert^2 \psi+\alpha_5\vert \psi\vert^4 \psi $$ in three spatial dimensions in the class of solutions with $|\psi(x)|\to c >0$ as $|x|\to\infty$. This models disturbances in an infinite expanse of (quantum) fluid in its quiescent state --- the limiting modulus $c$ corresponds to a local minimum in the energy density. Our arguments build on work of Gustafson, Nakanishi, and Tsai on the (defocusing) Gross--Pitaevskii equation. The presence of an energy-critical nonlinearity and changes in the geometry of the energy functional add several new complexities. One new ingredient in our argument is a demonstration that solutions of such (perturbed) energy-critical equations exhibit continuous dependence on the initial data with respect to the \emph{weak} topology on $H^1_x$.


title: Asymptotics with a positive cosmological constant: II. Linear fields on de Sitter space-time

abstract: Linearized gravitational waves in de Sitter space-time are analyzed in detail to obtain guidance for constructing the theory of gravitational radiation in presence of a positive cosmological constant in full, nonlinear general relativity. Specifically: i) In the exact theory, the intrinsic geometry of $\scri$ is often assumed to be conformally flat in order to reduce the asymptotic symmetry group from $\Diff$ to the de Sitter group. Our {results show explicitly} that this condition is physically unreasonable; ii) We obtain expressions of energy-momentum and angular momentum fluxes carried by gravitational waves in terms of fields defined at $\scrip$; iii) We argue that, although energy of linearized gravitational waves can be arbitrarily negative in general, gravitational waves emitted by physically reasonable sources carry positive energy; and, finally iv) We demonstrate that the flux formulas reduce to the familiar ones in Minkowski space-time in spite of the fact that the limit $\Lambda \to 0$ is discontinuous (since, in particular, $\scri$ changes its space-like character to null in the limit).


title: High Energy Collisions of Black Holes Numerically Revisited

abstract: We use fully nonlinear numerical relativity techniques to study high energy head-on collision of nonspinning, equal-mass black holes to estimate the maximum gravitational radiation emitted by these systems. Our simulations include improvements in the construction of initial data, subsequent full numerical evolutions, and the computation of waveforms at infinity. The new initial data significantly reduces the spurious radiation content, allowing for initial speeds much closer to the speed of light, i.e. $v\sim0.99c$. Using these new techniques, We estimate the maximum radiated energy from head-on collisions to be $E_{\text{max}}/M_{\text{ADM}}=0.13\pm0.01$. This value differs from the second-order perturbative $(0.164)$ and zero-frequency-limit $(0.17)$ analytic computations, but is close to those obtained by thermodynamic arguments $(0.134)$ and by previous numerical estimates $(0.14\pm0.03)$.


title: Network Coding Over SATCOM: Lessons Learned

abstract: Satellite networks provide unique challenges that can restrict users' quality of service. For example, high packet erasure rates and large latencies can cause significant disruptions to applications such as video streaming or voice-over-IP. Network coding is one promising technique that has been shown to help improve performance, especially in these environments. However, implementing any form of network code can be challenging. This paper will use an example of a generation-based network code and a sliding-window network code to help highlight the benefits and drawbacks of using one over the other. In-order packet delivery delay, as well as network efficiency, will be used as metrics to help differentiate between the two approaches. Furthermore, lessoned learned during the course of our research will be provided in an attempt to help the reader understand when and where network coding provides its benefits.


title: CO2 Forest: Improved Random Forest by Continuous Optimization of Oblique Splits

abstract: We propose a novel algorithm for optimizing multivariate linear threshold functions as split functions of decision trees to create improved Random Forest classifiers. Standard tree induction methods resort to sampling and exhaustive search to find good univariate split functions. In contrast, our method computes a linear combination of the features at each node, and optimizes the parameters of the linear combination (oblique) split functions by adopting a variant of latent variable SVM formulation. We develop a convex-concave upper bound on the classification loss for a one-level decision tree, and optimize the bound by stochastic gradient descent at each internal node of the tree. Forests of up to 1000 Continuously Optimized Oblique (CO2) decision trees are created, which significantly outperform Random Forest with univariate splits and previous techniques for constructing oblique trees. Experimental results are reported on multi-class classification benchmarks and on Labeled Faces in the Wild (LFW) dataset.


title: A novel ZnO piezoelectric microcantilever energy scavenger: Fabrication and characterization

abstract: This novel piezoelectric zinc oxide (ZnO) thin film microcantilever was fabricated by micromachining technique. To release the cantilever, wet anisotropic etching of Silicon (Si) was performed by tetramethyl ammonium hydroxide (TMAH). The transverse piezoelectric coefficient d31 of the ZnO film, obtained from the deflection of the cantilever with influence of applied voltage, was calculated as 3.32 pC/N. The observed dynamic characterization of the novel piezoelectric microcantilever had linear response with the applied driving voltage. The obtained values of Young Modulus and Hardness were 208 +4 GPa and 4.84+ 0.1 GPa respectively. This inexpensive novel method provides additional design flexibility to fabricate vibrational energy harvesters. The easy steps of fabrication and cost effectiveness of this method may be preferred it over DRIE. The voltage induced due to deformation of ZnO cantilever were measured ~230mV. This microcantilever energy scavenger may be used to power the nano devices and sensors for medical and agricultural applications as a replacement of traditional bulky batteries.


title: Complexity of a Disjoint Matching Problem on Bipartite Graphs

abstract: We consider the following question: given an $(X,Y)$-bigraph $G$ and a set $S \subset X$, does $G$ contain two disjoint matchings $M_1$ and $M_2$ such that $M_1$ saturates $X$ and $M_2$ saturates $S$? When $|S|\geq |X|-1$, this question is solvable by finding an appropriate factor of the graph. In contrast, we show that when $S$ is allowed to be an arbitrary subset of $X$, the problem is NP-hard.


title: Structured Training for Neural Network Transition-Based Parsing

abstract: We present structured perceptron training for neural network transition-based dependency parsing. We learn the neural network representation using a gold corpus augmented by a large number of automatically parsed sentences. Given this fixed network representation, we learn a final layer using the structured perceptron with beam-search decoding. On the Penn Treebank, our parser reaches 94.26% unlabeled and 92.41% labeled attachment accuracy, which to our knowledge is the best accuracy on Stanford Dependencies to date. We also provide in-depth ablative analysis to determine which aspects of our model provide the largest gains in accuracy.


title: Evaluation of the Number of Different Genomes on Medium and Identification of Known Genomes Using Composition Spectra Approach

abstract: The article presents the theoretical foundations of the algorithm for calculating the number of different genomes in the medium under study and of two algorithms for determining the presence of a particular (known) genome in this medium. The approach is based on the analysis of the compositional spectra of subsequently sequenced samples of the medium. The theoretical estimations required for the implementation of the algorithms are obtained.


title: The interplay between resonant enhancement and quantum path dynamics in harmonic generation in helium

abstract: We present a theoretical study of the influence of resonant enhancement on quantum path dynamics in the generation of harmonics above and below the ionization threshold in helium. By varying the wavelength and intensity of the driving field from 425 nm to 500 nm and from 30 TW/cm${^2}$ to 140 TW/cm${^2}$, respectively, we identify enhancements of harmonics 7, 9, and 11 that correspond to multiphoton resonances between the ground state and the Stark shifted $1s2p$, $1s3p$, and $1s4p$ excited states. A time-frequency analysis of the emission shows that both the short and long quantum path contributions to the harmonic yield are enhanced through these bound state resonances. We analyze the sub-cycle time structure of the 9th harmonic yield in the vicinity of the resonances and find that on resonance the long trajectory contribution is phase shifted by approximately $\pi/4$. Finally, we compare the single atom and the macroscopic response of a helium gas and find that while the sub-cycle time profiles are slightly distorted by propagation effects, the phase shift of the long trajectory contribution is still recognizable.


title: The Lerch zeta function III. Polylogarithms and special values

abstract: This paper studies algebraic and analytic structures associated with the Lerch zeta function, extending the complex variables viewpoint taken in part II. The Lerch transcendent $\Phi(s, z, c)$ is obtained from the Lerch zeta function $\zeta(s, a, c)$ by the change of variable $z=e^{2 \pi i a}$. We show that it analytically continues to a maximal domain of holomorphy in three complex variables $(s, z, c)$ as a multivalued function defined over the base manifold ${\bf C} \times P^1({\bf C} \smallsetminus \{0, 1, \infty\}) \times ({\bf C}\smallsetminus {\bf Z})$. and compute the monodromy functions defining the multivaluedness. For positive integer values s=m and c=1 this function is closely related to the classical m-th order polylogarithm $Li_m(z)$ We study its behavior as a function of two variables $(z, c)$ for special values where s=m is an integer. For $m \ge 1$ it gives a one-parameter deformation of the polylogarithm, and satisfies a linear ODE with coefficients depending on c, of order m+1 of Fuchsian type. We determine its (m+1)-dimensional monodromy representation, which is a non-algebraic deformation of the monodromy of $Li_m(z).$


title: Private Graphon Estimation for Sparse Graphs

abstract: We design algorithms for fitting a high-dimensional statistical model to a large, sparse network without revealing sensitive information of individual members. Given a sparse input graph $G$, our algorithms output a node-differentially-private nonparametric block model approximation. By node-differentially-private, we mean that our output hides the insertion or removal of a vertex and all its adjacent edges. If $G$ is an instance of the network obtained from a generative nonparametric model defined in terms of a graphon $W$, our model guarantees consistency, in the sense that as the number of vertices tends to infinity, the output of our algorithm converges to $W$ in an appropriate version of the $L_2$ norm. In particular, this means we can estimate the sizes of all multi-way cuts in $G$. Our results hold as long as $W$ is bounded, the average degree of $G$ grows at least like the log of the number of vertices, and the number of blocks goes to infinity at an appropriate rate. We give explicit error bounds in terms of the parameters of the model; in several settings, our bounds improve on or match known nonprivate results.


title: Big Step Greedy Algorithm for Maximum Coverage Problem

abstract: This paper proposes a greedy heuristic named as Big step greedy heuristic and investigates the application of Big step greedy heuristic for maximum k-coverage problem. Greedy algorithms construct the solution in multiple steps, the classical greedy algorithm for maximum k-coverage problem, in each step selects one set that contains the greatest number of uncovered elements. The Big step greedy heuristic, in each step selects p (1 <= p <= k) sets such that the union of selected p sets contains the greatest number of uncovered elements by evaluating all possible p-combinations of given sets. When p=k Big step greedy algorithm behaves like exact algorithm that computes optimal solution by evaluating all possible k-combinations of given sets. When p=1 it behaves like the classical greedy algorithm.


title: A Uniform Contribution of Core-Collapse and Type Ia Supernovae to the Chemical Enrichment Pattern in the Outskirts of the Virgo Cluster

abstract: We present the first measurements of the abundances of $\alpha$-elements (Mg, Si, and S) extending out to beyond the virial radius of a cluster of galaxies. Our results, based on Suzaku Key Project observations of the Virgo Cluster, show that the chemical composition of the intra-cluster medium is consistent with being constant on large scales, with a flat distribution of the Si/Fe, S/Fe, and Mg/Fe ratios as a function of radius and azimuth out to 1.4 Mpc (1.3 $r_{200}$). Chemical enrichment of the intergalactic medium due solely to core collapse supernovae (SNcc) is excluded with very high significance; instead, the measured metal abundance ratios are generally consistent with the Solar value. The uniform metal abundance ratios observed today are likely the result of an early phase of enrichment and mixing, with both SNcc and type Ia supernovae (SNIa) contributing to the metal budget during the period of peak star formation activity at redshifts of 2-3. We estimate the ratio between the number of SNIa and the total number of supernovae enriching the intergalactic medium to be between 12-37%, broadly consistent with the metal abundance patterns in our own Galaxy or with the SNIa contribution estimated for the cluster cores.


title: Abstract Model Repair

abstract: Given a Kripke structure M and CTL formula $\varphi$, where M does not satisfy $\varphi$, the problem of Model Repair is to obtain a new model M' such that M' satisfies $\varphi$. Moreover, the changes made to M to derive M' should be minimum with respect to all such M'. As in model checking, state explosion can make it virtually impossible to carry out model repair on models with infinite or even large state spaces. In this paper, we present a framework for model repair that uses abstraction refinement to tackle state explosion. Our framework aims to repair Kripke Structure models based on a Kripke Modal Transition System abstraction and a 3-valued semantics for CTL. We introduce an abstract-model-repair algorithm for which we prove soundness and semi-completeness, and we study its complexity class. Moreover, a prototype implementation is presented to illustrate the practical utility of abstract-model-repair on an Automatic Door Opener system model and a model of the Andrew File System 1 protocol.


title: A Type-Theoretic Approach to Structural Resolution

abstract: Structural resolution (or S-resolution) is a newly proposed alternative to SLD-resolution that allows a systematic separation of derivations into term-matching and unification steps. Productive logic programs are those for which term-matching reduction on any query must terminate. For productive programs with coinductive meaning, finite term-rewriting reductions can be seen as measures of observation in an infinite derivation. Ability of handling corecursion in a productive way is an attractive computational feature of S-resolution. In this paper, we make first steps towards a better conceptual understanding of operational properties of S-resolution as compared to SLD-resolution. To this aim, we propose a type system for the analysis of both SLD-resolution and S-resolution. We formulate S-resolution and SLD-resolution as reduction systems, and show their soundness relative to the type system. One of the central methods of this paper is realizability transformation, which makes logic programs productive and non-overlapping. We show that S-resolution and SLD-resolution are only equivalent for programs with these two properties.


title: Fractional dynamics on networks: Emergence of anomalous diffusion and L\'evy flights

abstract: We introduce a formalism of fractional diffusion on networks based on a fractional Laplacian matrix that can be constructed directly from the eigenvalues and eigenvectors of the Laplacian matrix. This fractional approach allows random walks with long-range dynamics providing a general framework for anomalous diffusion and navigation, and inducing dynamically the small-world property on any network. We obtained exact results for the stationary probability distribution, the average fractional return probability and a global time, showing that the efficiency to navigate the network is greater if we use a fractional random walk in comparison to a normal random walk. For the case of a ring, we obtain exact analytical results showing that the fractional transition and return probabilities follow a long-range power-law decay, leading to the emergence of L\'evy flights on networks. Our general fractional diffusion formalism applies to regular, random and complex networks and can be implemented from the spectral properties of the Laplacian matrix, providing an important tool to analyze anomalous diffusion on networks.


title: Direct Measure of Giant Magnetocaloric Entropy Contributions in Ni-Mn-In

abstract: Off-stoichiometric alloys based on Ni 2 MnIn have drawn attention due to the coupled first order magnetic and structural transformations, and the large magnetocaloric entropy associated with the transformations. Here we describe calorimetric and magnetic studies of four compositions. The results provide a direct measure of entropy changes contributions including at the first-order phase transitions, and thereby a determination of the maximum field-induced entropy change corresponding to the giant magnetocaloric effect. We find a large excess entropy change, attributed to magneto-elastic coupling, but only in compositions with no ferromagnetic order in the high-temperature austenite phase. Furthermore, a molecular field model corresponding to antiferromagnetism of the low-temperature phases is in good agreement, and nearly independent of composition, despite significant differences in overall magnetic response of these materials.


title: A Model-Based Approach for Analog Spatio-Temporal Dynamic Forecasting

abstract: Analog forecasting has been applied in a variety of fields for predicting future states of complex nonlinear systems that require flexible forecasting methods. Past analog methods have almost exclu- sively been used in an empirical framework without the structure of a model-based approach. We propose a Bayesian model framework for analog forecasting, building upon previous analog methods but accounting for parameter uncertainty. Thus, unlike traditional analog forecasting methods, the use of Bayesian modeling allows one to rigorously quantify uncertainty to obtain realistic posterior predictive distributions. The model is applied to the long-lead time forecasting of mid-May averaged soil moisture anomalies in Iowa over a high-resolution grid of spatial locations. Sea Surface Tem- perature (SST) is used to find past time periods with similar trajectories to the current pre-forecast period. The analog model is developed on projection coefficients from a basis expansion of the soil moisture and SST fields. Separate models are constructed for locations falling in each Iowa Crop Reporting District (CRD) and the forecasting ability of the proposed model is compared against a variety of alternative methods and metrics.