title: Optimal Combination of Image Denoisers

abstract: Given a set of image denoisers, each having a different denoising capability, is there a provably optimal way of combining these denoisers to produce an overall better result? An answer to this question is fundamental to designing an ensemble of weak estimators for complex scenes. In this paper, we present an optimal combination scheme by leveraging deep neural networks and convex optimization. The proposed framework, called the Consensus Neural Network (CsNet), introduces three new concepts in image denoising: (1) A provably optimal procedure to combine the denoised outputs via convex optimization; (2) A deep neural network to estimate the mean squared error (MSE) of denoised images without needing the ground truths; (3) An image boosting procedure using a deep neural network to improve contrast and to recover lost details of the combined images. Experimental results show that CsNet can consistently improve denoising performance for both deterministic and neural network denoisers.


title: The high brightness temperature of B0529+483 revealed by RadioAstron and implications for interstellar scattering

abstract: The high brightness temperatures, $T_\mathrm{b}\gtrsim 10^{13}$ K, detected in several active galactic nuclei by RadioAstron space VLBI observations challenge theoretical limits. Refractive scattering by the interstellar medium may affect such measurements. We quantify the scattering properties and the sub-mas scale source parameters for the quasar B0529+483. Using RadioAstron correlated flux density measurements at 1.7, 4.8, and 22 GHz on projected baselines up to 240,000 km we find two characteristic angular scales in the quasar core, about 100 $\mu$as and 10 $\mu$as. Some indications of scattering substructure are found. Very high brightness temperatures, $T_\mathrm{b}\geq 10^{13}$ K, are estimated at 4.8 GHz and 22 GHz even taking into account the refractive scattering. Our findings suggest a clear dominance of the particle energy density over the magnetic field energy density in the core of this quasar.


title: Electron heating in the laser and static electric and magnetic fields

abstract: A 2D slab approximation of the interactions of electrons with intense linearly polarized laser radiation and static electric and magnetic fields is widely used for both numerical simulations and simplified semi-analytical models. It is shown that in this case electron dynamics can be conveniently described in the framework of 3/2 dimensional Hamiltonian approach. The electron acceleration beyond a standard ponderomotive scaling, caused by the synergistic effects of the laser and static electro-magnetic fields, is due to an onset of stochastic electron motion.


title: Current-carrying states in FFLO superconductors

abstract: We show that nonuniform superconductors of the Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) type conduct electric current in the way which is very different from the usual case. We discuss both equilibrium and nonequilibrium properties using a modified Ginzburg-Landau formalism. Among the novel features are the existence of two different critical currents and two distinct stable states able to carry a given current, the possibility of superconducting domain walls, and also a spontaneous supercurrent in a ring geometry.


title: An Upper Bound for the Depth of Some Classes of Polyhedra

abstract: K. Borsuk in the seventies introduced the notions of capacity and depth of compacta together with some relevant problems. In this paper, first, we introduce the concepts of the (strong) capacity and the (strong) depth of an object in an arbitrary category. Then in the category of groups, we compute the (strong) capacity and the (strong) depth of some well-known groups. Finally, we find an upper bound for the depth of some classes of finite polyhedra which generalizes a result of D. Kolodziejczyk in this subject.


title: The two-mass contribution to the three-loop pure singlet operator matrix element

abstract: We present the two-mass QCD contributions to the pure singlet operator matrix element at three loop order in x-space. These terms are relevant for calculating the structure function $F_2(x,Q^2)$ at $O(\alpha_s^3)$ as well as for the matching relations in the variable flavor number scheme and the heavy quark distribution functions at the same order. The result for the operator matrix element is given in terms of generalized iterated integrals that include square root letters in the alphabet, depending also on the mass ratio through the main argument. Numerical results are presented.


title: The Higgs boson in the minimal 3-3-1 model

abstract: In this paper we present the mass matrices and mass eigenstates for the CP-even neutral scalars in the minimal 331 model (m331) and its self-interactions, showing that the m331 automatically reproduces the Higgs potential of the Standard Model. We also present a method to generate numerical solutions for the quarks and leptons masses and their mixings, which we apply to study FCNC processes, being to calculate the contributions of all exotic neutral particles of the m331 to the mass differences in meson oscillations.


title: Techniques for proving Asynchronous Convergence results for Markov Chain Monte Carlo methods

abstract: Markov Chain Monte Carlo (MCMC) methods such as Gibbs sampling are finding widespread use in applied statistics and machine learning. These often lead to difficult computational problems, which are increasingly being solved on parallel and distributed systems such as compute clusters. Recent work has proposed running iterative algorithms such as gradient descent and MCMC in parallel asynchronously for increased performance, with good empirical results in certain problems. Unfortunately, for MCMC this parallelization technique requires new convergence theory, as it has been explicitly demonstrated to lead to divergence on some examples. Recent theory on Asynchronous Gibbs sampling describes why these algorithms can fail, and provides a way to alter them to make them converge. In this article, we describe how to apply this theory in a generic setting, to understand the asynchronous behavior of any MCMC algorithm, including those implemented using parameter servers, and those not based on Gibbs sampling.


title: Quarnet inference rules for level-1 networks

abstract: An important problem in phylogenetics is the construction of phylogenetic trees. One way to approach this problem, known as the supertree method, involves inferring a phylogenetic tree with leaves consisting of a set $X$ of species from a collection of trees, each having leaf-set some subset of $X$. In the 1980's characterizations, certain inference rules were given for when a collection of 4-leaved trees, one for each 4-element subset of $X$, can all be simultaneously displayed by a single supertree with leaf-set $X$. Recently, it has become of interest to extend such results to phylogenetic networks. These are a generalization of phylogenetic trees which can be used to represent reticulate evolution (where species can come together to form a new species). It has been shown that a certain type of phylogenetic network, called a level-1 network, can essentially be constructed from 4-leaved trees. However, the problem of providing appropriate inference rules for such networks remains unresolved. Here we show that by considering 4-leaved networks, called quarnets, as opposed to 4-leaved trees, it is possible to provide such rules. In particular, we show that these rules can be used to characterize when a collection of quarnets, one for each 4-element subset of $X$, can all be simultaneously displayed by a level-1 network with leaf-set $X$. The rules are an intriguing mixture of tree inference rules, and an inference rule for building up a cyclic ordering of $X$ from orderings on subsets of $X$ of size 4. This opens up several new directions of research for inferring phylogenetic networks from smaller ones, which could yield new algorithms for solving the supernetwork problem in phylogenetics.


title: Learning SO(3) Equivariant Representations with Spherical CNNs

abstract: We address the problem of 3D rotation equivariance in convolutional neural networks. 3D rotations have been a challenging nuisance in 3D classification tasks requiring higher capacity and extended data augmentation in order to tackle it. We model 3D data with multi-valued spherical functions and we propose a novel spherical convolutional network that implements exact convolutions on the sphere by realizing them in the spherical harmonic domain. Resulting filters have local symmetry and are localized by enforcing smooth spectra. We apply a novel pooling on the spectral domain and our operations are independent of the underlying spherical resolution throughout the network. We show that networks with much lower capacity and without requiring data augmentation can exhibit performance comparable to the state of the art in standard retrieval and classification benchmarks.


title: N=2 Moduli of AdS4 vacua: A fine-print study

abstract: We analyze the moduli spaces near N=4 supersymmetric AdS4 vacua of string theory by combining and comparing various approaches: (a) the known exact solutions of Type IIB string theory with localized 5-brane sourcees, (b) the holographically dual 3d quiver gauge theories; (c) gauged supergravity, and (d) the representations of the superconformal algebra Osp(4|4). Short multiplets containing the marginal N=2 deformations transform in the (2;0), (0;2) or (1;1) representations of the R-symmetry group SU(2)_H x SU(2)_C. The first two are classified by the chiral rings of the Higgs and Coulomb branches, while the latter contain mixed-branch operators. We identify the origin of these moduli in string theory, matching in particular the operators of the chiral rings with open strings on the magnetized 5-brane sources. Our results provide new evidence for the underlying holographic duality. The existence of a large number of bound-state moduli highlights the limitations of effective supergravity.


title: Towards a New Perturbation Theory of Molecular Liquids

abstract: First steps towards developing a new perturbation theory for molecular liquids are taken. By choosing a new form of splitting the site-site potential functions between molecules, we will get a set of atomic fluids as the reference system with known structure and thermodynamics. The perturbative part of the potential function is then expanded up to two terms. The excess Helmholtz free energy of the system is then obtained for three computable contributions. The derivation shows that the excess Helmholtz free energy has nothing to do with intra-atomic potentials, all contributions come merely from inter-atomic potentials. Then it is applied to compute the thermodynamics of two systems; hard sphere chain and carbon dioxide molecular fluids. The results compared with the computer simulation data show that the theory works well at low densities.


title: An Equation for Predicting Binding Strengths of Metal Cations to Protein of Human Serum Transferrin

abstract: Because human serum transferrin (hTF) exists freely in serum, it is a potential target for cancer treatment drugs and in curing iron-overloaded conditions in patients via long-term transfusion therapy. The understanding of the interactions between hTF and metal ions is very important for biological, pharmalogical, toxicological, and other protein engineering purposes. In this paper, a simple linear free energy correlation is proposed to predict the binding strength between hTF protein and metal cations. The stability constants for a family of metal-hTF complexes can be correlated to the non-solvation energies and the radii of cations. The binding strength is determined by both the physical properties (charge and size or ionic radius) and chemical properties (non-solvation energy) of a given cation. The binding strengths of either divalent and trivalent metals can then be predicted systematically.


title: Frustrated magnetism in tetragonal CoSe, analogue to superconducting FeSe

abstract: Recently synthesized metastable tetragonal CoSe, isostructural to the FeSe superconductor, offers a new avenue for investigating systems in close proximity to the iron-based superconductors. We present magnetic and transport property measurements on powders and single crystals of CoSe. High field magnetic susceptibility measurements indicate a suppression of the previously reported 10 K ferromagnetic transition with the magnetic susceptibility exhibiting time-dependence below the proposed transition. Dynamic scaling analysis of the time-dependence yields a critical relaxation time of $\tau^{*} = 0.064 \pm 0.008 $ s which in turn yields an activation energy of $E_{a}^{*}$ = 14.84 $\pm$ 0.59 K and an ideal glass temperature $T_{0}^{*}$ = 8.91 $\pm$ 0.09 K from Vogel-Fulcher analysis. No transition is observed in resistivity and specific heat measurements, but both measurements indicate that CoSe is metallic. These results are interpreted on the basis of CoSe exhibiting frustrated magnetic ordering arising from competing magnetic interactions. Arrott analysis of single crystal magnetic susceptibility has indicated that the magnetic moments lie in the $ab$-plane so frustration may arise from intralayer magnetic fluctuations as well as interlayer coupling. The results have implications for understanding the superconductivity in the iron chalcogenide systems as well as utilizing CoSe as a host for chemical and physical manipulation to tune and explore emergent phenomena within an expanding new class of transition metal chalcogenides.


title: A simple hysteretic constitutive model for unsaturated flow

abstract: In this paper we present a constitutive model to describe unsaturated flow that considers the hysteresis phenomena. This constitutive model provides simple mathematical expressions for both saturation and hydraulic conductivity curves, and a relationship between permeability and porosity. The model is based on the assumption that the porous media can be represented by a bundle of capillary tubes with throats or "ink-bottles" and a fractal pore size distribution. Under these hypotheses, hysteretic curves are obtained for saturation and relative hydraulic conductivity in terms of pressure head. However, a non-hysteretic relationship is obtained when relative hydraulic conductivity is expressed as a function of saturation. The proposed relationship between permeability and porosity is similar to the well-known Kozeny-Carman equation but depends on the fractal dimension. The performance of the constitutive model is tested against different sets of experimental data and previous models. In all of the cases the proposed expressions fit fairly well the experimental data and predicts values of permeability and hydraulic conductivity better than others models.


title: Hamiltonian Renormalisation II. Renormalisation Flow of 1+1 dimensional free scalar fields: Derivation

abstract: In the companion paper we motivated a renormalisation flow on Osterwalder-Schrader data (OS-data) consisting of 1. a Hilbert space, 2. a cyclic vacuum and 3. a Hamiltonian annihilating that vacuum. As the name suggests, the motivation was via the OS reconstruction theorem which allows to reconstruct the OS data from an OS measure satisfying (a subset of) the OS axioms, in particular reflection positivity. The guiding principle was to map the usual Wilsonian path integral renormalisation flow onto a flow of the corresponding OS data. We showed that this induced flow on the OS data has an unwanted feature which disqualifies the associated coarse grained Hamiltonians from being the projections of a continuum Hamiltonian onto vectors in the coarse grained Hilbert space. This motivated the definition of a direct Hamiltonian renormalisation flow which follows the guiding principle but does not suffer from the afore mentioned caveat. In order to test our proposal, we apply it to the only known completely solvable model, namely the case of free scalar quantum fields. In this paper we focus on the Klein Gordon field in two spacetime dimensions and illustrate the difference between the path integral induced and direct Hamiltonian flow. Generalisations to more general models in higher dimensions will be discussed in our companion papers.

categories: gr-qc hep-lat hep-th math-ph math.MP quant-ph


title: The stellar orbit distribution in present-day galaxies inferred from the CALIFA survey

abstract: Galaxy formation entails the hierarchical assembly of mass, along with the condensation of baryons and the ensuing, self-regulating star formation. The stars form a collisionless system whose orbit distribution retains dynamical memory that can constrain a galaxy's formation history. The ordered-rotation dominated orbits with near maximum circularity $\lambda_z \simeq1$ and the random-motion dominated orbits with low circularity $\lambda_z \simeq0$ are called kinematically cold and kinematically hot, respectively. The fraction of stars on `cold' orbits, compared to the fraction of stars on `hot' orbits, speaks directly to the quiescence or violence of the galaxies' formation histories. Here we present such orbit distributions, derived from stellar kinematic maps via orbit-based modelling for a well defined, large sample of 300 nearby galaxies. The sample, drawn from the CALIFA survey, includes the main morphological galaxy types and spans the total stellar mass range from $10^{8.7}$ to $10^{11.9}$ solar masses. Our analysis derives the orbit-circularity distribution as a function of galaxy mass, $p(\lambda_z~|~M_\star)$, and its volume-averaged total distribution, $p(\lambda_z)$. We find that across most of the considered mass range and across morphological types, there are more stars on `warm' orbits defined as $0.25\le \lambda_z \le 0.8$ than on either `cold' or `hot' orbits. This orbit-based "Hubble diagram" provides a benchmark for galaxy formation simulations in a cosmological context.


title: Phonological (un)certainty weights lexical activation

abstract: Spoken word recognition involves at least two basic computations. First is matching acoustic input to phonological categories (e.g. /b/, /p/, /d/). Second is activating words consistent with those phonological categories. Here we test the hypothesis that the listener's probability distribution over lexical items is weighted by the outcome of both computations: uncertainty about phonological discretisation and the frequency of the selected word(s). To test this, we record neural responses in auditory cortex using magnetoencephalography, and model this activity as a function of the size and relative activation of lexical candidates. Our findings indicate that towards the beginning of a word, the processing system indeed weights lexical candidates by both phonological certainty and lexical frequency; however, later into the word, activation is weighted by frequency alone.


title: Quantitative unique continuation for a parabolic equation

abstract: We address the quantitative uniqueness properties of the solutions of the parabolic equation $ \partial_t u - \Delta u = w_j (x,t) \partial_j u + v(x,t) u $ where $v$ and $w$ are bounded. We prove that for solutions $u$, the order of vanishing is bounded by $C(\Vert v\Vert_{L^\infty}^{2/3}+\Vert w\Vert_{L^\infty}^2)$ matching the upper bound previously established in the elliptic case. in the elliptic case.


title: Hyperuniformity with no fine tuning in sheared sedimenting suspensions

abstract: Particle suspensions, present in many natural and industrial settings, typically contain aggregates or other microstructures that can complicate macroscopic flow behaviors and damage processing equipment. Recent work found that applying uniform periodic shear near a critical transition can reduce fluctuations in the particle concentration across all length scales, leading to a hyperuniform state. However, this strategy for homogenization requires fine tuning of the strain amplitude. Here we show that in a model of sedimenting particles under periodic shear, there is a well-defined regime at low sedimentation speed where hyperuniform scaling automatically occurs. Our simulations and theoretical arguments show that the homogenization extends up to a finite lengthscale that diverges as the sedimentation speed approaches zero.


title: Modified hoop conjecture in expanding spacetimes and primordial black hole production in FRW universe

abstract: According to a variant of the hoop conjecture, if we localize two particles within the Schwarzschild radius corresponding to their center of mass energy, then a black hole will form. Despite a large body of work on the formation of primordial black holes, so far this conjecture has not been generalized to expanding spacetimes. We derive a formula which gives the distance within which two particles must be localized to give a black hole, and which crucially depends on the expansion rate of the background space. In the limit of a very slow expansion, we recover the flat spacetime case. In the opposite limit of the large expansion rate when the inverse Hubble radius is smaller than the Schwarzschild radius of a "would be" black hole, the new critical distance between two particles that can make a black hole becomes equal to the particle horizon, which is just a requirement that the particles are in a causal contact. This behavior also nicely illustrates why the Big Bang singularity is not a black hole. We then use our formula to calculate the number density, energy density and production rate of black holes produced in collisions of particles. We find that though black holes might be numerous at high temperatures, they never dominate over the background radiation below the Planck temperature.


title: Hopping in the crowd to unveil network topology

abstract: We introduce a nonlinear operator to model diffusion on a complex undirected network under crowded conditions. We show that the asymptotic distribution of diffusing agents is a nonlinear function of the nodes' degree and saturates to a constant value for sufficiently large connectivities, at variance with standard diffusion in the absence of excluded-volume effects. Building on this observation, we define and solve an inverse problem, aimed at reconstructing the a priori unknown connectivity distribution. The method gathers all the necessary information by repeating a limited number of independent measurements of the asymptotic density at a single node that can be chosen randomly. The technique is successfully tested against both synthetic and real data, and shown to estimate with great accuracy also the total number of nodes.


title: MWC349A and B Are Not Gravitationally Bound: New Evidence

abstract: The age and evolutionary status of MWC349A, the unique emission line star with maser and laser radiation in hydrogen recombination lines, remain unknown because the spectrum of the star is veiled by bright emission from the ionized disk and wind. The major argument for this massive ($>$10 M_sun) star being evolved is its association with a close-by (2.4 arcsec) companion, MWC349B, whose B0III spectrum implies an age of a few Myrs. However, newly obtained high-resolution spectra of MWC349B reveal a difference $\approx$ 35 km/s in the radial velocities of the two stars, which makes their being gravitationally bound highly improbable. An estimate of the relative proper motion of the two stars seems to confirm this conclusion. This reopens the previously suggested possibility that MWC349A is a young massive star in a region of active star formation close to the Cyg OB2 association. MWC349B, which moves with a speed $\geq$ 35 km/s relative to Cyg OB2, may be a runaway star from this association.


title: Endomorphisms of regular rooted trees induced by the action of polynomials on the ring $\mathbb Z_d$ of $d$-adic integers

abstract: We show that every polynomial in $\mathbb Z[x]$ defines an endomorphism of the $d$-ary rooted tree induced by its action on the ring $\mathbb Z_d$ of $d$-adic integers. The sections of this endomorphism also turn out to be induced by polynomials in $\mathbb Z[x]$ of the same degree. In the case of permutational polynomials acting on $\mathbb Z_d$ by bijections the induced endomorphisms are automorphisms of the tree. In the case of $\mathbb Z_2$ such polynomials were completely characterized by Rivest. As our main application we utilize the result of Rivest to derive the condition on the coefficients of a permutational polynomial $f(x)\in\mathbb Z[x]$ that is necessary and sufficient for $f$ to induce a level transitive automorphism of the binary tree, which is equivalent to the ergodicity of the action of $f(x)$ on $\mathbb Z_2$ with respect to the normalized Haar measure.


title: A Hint of a Blue Axion Isocurvature Spectrum?

abstract: It is known that if the Peccei-Quinn symmetry breaking field is displaced from its minimum during inflation, the axion isocurvature spectrum is generically strongly blue tilted with a break transition to a flat spectrum. We fit this spectrum (incorporated into the "vanilla" $\Lambda$-CDM cosmological model) to the Planck and BOSS DR11 data and find a mild hint for the presence of axionic blue-tilted isocurvature perturbations. We find the best fit parameter region is consistent with all of the dark matter being composed of QCD axions in the context of inflationary cosmology with an expansion rate of order $10^{8}$ GeV, the axion decay constant of order $10^{13}$ GeV, and the initial misalignment angle of order unity. Intriguingly, isocurvature with a spectral break may at least partially explain the low-$\ell$ vs. high-$\ell$ anomalies seen in the CMB data.


title: Numerical Simulations of Turbulent Molecular Clouds Regulated by Radiation Feedback Forces II: Radiation-Gas Interactions and Outflows

abstract: Momentum deposition by radiation pressure from young, massive stars may help to destroy molecular clouds and unbind stellar clusters by driving large-scale outflows. We extend our previous numerical radiation hydrodynamic study of turbulent, star-forming clouds to analyze the detailed interaction between non-ionizing UV radiation and the cloud material. Our simulations trace the evolution of gas and star particles through self-gravitating collapse, star formation, and cloud destruction via radiation-driven outflows. These models are idealized in that we include only radiation feedback and adopt an isothermal equation of state. Turbulence creates a structure of dense filaments and large holes through which radiation escapes, such that only ~50% of the radiation is (cumulatively) absorbed by the end of star formation. The surface density distribution of gas by mass as seen by the central cluster is roughly lognormal with sigma_ln(Sigma) = 1.3-1.7, similar to the externally-projected surface density distribution. This allows low surface density regions to be driven outwards to nearly 10 times their initial escape speed v_esc. Although the velocity distribution of outflows is broadened by the lognormal surface density distribution, the overall efficiency of momentum injection to the gas cloud is reduced because much of the radiation escapes. The mean outflow velocity is approximately twice the escape speed from the initial cloud radius. Our results are also informative for understanding galactic-scale wind driving by radiation, in particular the relationship between velocity and surface density for individual outflow structures, and the resulting velocity and mass distributions arising from turbulent sources.


title: Dark matter component decaying after recombination: Sensitivity to BAO and RSD probes

abstract: It has been recently suggested~\cite{Berezhiani:2015yta} that a subdominant fraction of dark matter decaying after recombination may alleviate tension between high-redshift (CMB anisotropy) and low-redshift (Hubble constant, cluster counts) measurements. In this report, we continue our previous study~\cite{Chudaykin:2016yfk} of the decaying dark matter (DDM) model adding all available recent baryon acoustic oscillation (BAO) and redshift space distortions (RSD) measurements. We find, that the BAO/RSD measurements generically prefer the standard $\Lambda$CDM and combined with other cosmological measurements impose an upper limit on the DDM fraction at the level of $\sim$\,5\,\%, strengthening by a factor of 1.5 limits obtained in \cite{Chudaykin:2016yfk} mostly from CMB data. However, the numbers vary from one analysis to other based on the same Baryon Oscillation Spectroscopic Survey (BOSS) Data Release 12 (DR12) galaxy sample. Overall, the model with a few percent DDM fraction provides a better fit to the combined cosmological data as compared to the $\Lambda$CDM: the cluster counting and direct measurements of the Hubble parameter are responsible for that. The improvement can be as large as 1.5\,$\sigma$ and grows to 3.3\,$\sigma$ when the CMB lensing power amplitude ${\rm A_L}$ is introduced as a free fitting parameter.


title: Schur's theory for partial projective representations

abstract: This article focuses on those aspects about partial actions of groups which are related to Schur's theory on projective representations. It provides an exhaustive description of the partial Schur multiplier, and this result is achieved by introducing the concept of a second partial cohomology group relative to an ideal, together with an appropriate analogue of a central extension. In addition, the new framework is proved to be consistent with the earlier notion of cohomology over partial modules.


title: Information Gathering with Peers: Submodular Optimization with Peer-Prediction Constraints

abstract: We study a problem of optimal information gathering from multiple data providers that need to be incentivized to provide accurate information. This problem arises in many real world applications that rely on crowdsourced data sets, but where the process of obtaining data is costly. A notable example of such a scenario is crowd sensing. To this end, we formulate the problem of optimal information gathering as maximization of a submodular function under a budget constraint, where the budget represents the total expected payment to data providers. Contrary to the existing approaches, we base our payments on incentives for accuracy and truthfulness, in particular, {\em peer-prediction} methods that score each of the selected data providers against its best peer, while ensuring that the minimum expected payment is above a given threshold. We first show that the problem at hand is hard to approximate within a constant factor that is not dependent on the properties of the payment function. However, for given topological and analytical properties of the instance, we construct two greedy algorithms, respectively called PPCGreedy and PPCGreedyIter, and establish theoretical bounds on their performance w.r.t. the optimal solution. Finally, we evaluate our methods using a realistic crowd sensing testbed.


title: Mechanisms of Coulomb dissociation processes

abstract: The Coulomb dissociation is studied of the $^8L$i nuclei on Pb target at energies 40.3 and 69.5 MeV/nucleon, in the experiment NSCL #03038. The $^{6,7}$Li, $^{4,6}$He, and $^2$H fragments were identified. The resonance decay and the direct breakup reactions were observed. The data give experimental evidence that the Coulomb dissociation is a two-step process. The projectile in the approaching phase is braked down and the valence neutron gets forced oscillation. The increasing Coulonb force holds the projectile and brings through the excited states up to the closest approach point. There released, the projectile is trapped into a primary, highly excited state in the continuum. The lifetime of the primary state depends on the multipolarity of the deformed projectile. At intermediate energy the collision is a sudden reaction, the valence neutron may stay --- during the impact --- in the forward or backward hemisphere of the nucleus, and the orbit of the valence neutron gets prolate or oblate, dipole or multipole deformation and the nucleus gets single-particle or collective excitation. The primary excited states, developing will decay prompt or delayed into the same reaction channel, resulting in the two decay mechanisms.


title: Measurements and Monte-Carlo simulations of the particle self-shielding effect of B4C grains in neutron shielding concrete

abstract: A combined measurement and Monte-Carlo simulation study was carried out in order to characterize the particle self-shielding effect of B4C grains in neutron shielding concrete. Several batches of a specialized neutron shielding concrete, with varying B4C grain sizes, were exposed to a 2 {\AA} neutron beam at the R2D2 test beamline at the Institute for Energy Technology located in Kjeller, Norway. The direct and scattered neutrons were detected with a neutron detector placed behind the concrete blocks and the results were compared to Geant4 simulations. The particle self-shielding effect was included in the Geant4 simulations by calculating effective neutron cross-sections during the Monte-Carlo simulation process. It is shown that this method well reproduces the measured results. Our results show that shielding calculations for low-energy neutrons using such materials would lead to an underestimate of the shielding required for a certain design scenario if the particle self-shielding effect is not included in the calculations.


title: The clumpy absorber in the high-mass X-ray binary Vela X-1

abstract: Bright and eclipsing, the high-mass X-ray binary Vela X-1 offers a unique opportunity to study accretion onto a neutron star from clumpy winds of O/B stars and to disentangle the complex accretion geometry of these systems. In Chandra-HETGS spectroscopy at orbital phase ~0.25, when our line of sight towards the source does not pass through the large-scale accretion structure such as the accretion wake, we observe changes in overall spectral shape on timescales of a few kiloseconds. This spectral variability is, at least in part, caused by changes in overall absorption and we show that such strongly variable absorption cannot be caused by unperturbed clumpy winds of O/B stars. We detect line features from high and low ionization species of silicon, magnesium and neon whose strengths and presence depend on the overall level of absorption. They imply a co-existence of cool and hot gas phases in the system that we interpret as a highly variable, structured accretion flow close to the compact object such as has been recently seen in simulations of wind accretion in high-mass X-ray binaries.


title: Learning to Organize Knowledge and Answer Questions with N-Gram Machines

abstract: Though deep neural networks have great success in natural language processing, they are limited at more knowledge intensive AI tasks, such as open-domain Question Answering (QA). Existing end-to-end deep QA models need to process the entire text after observing the question, and therefore their complexity in responding a question is linear in the text size. This is prohibitive for practical tasks such as QA from Wikipedia, a novel, or the Web. We propose to solve this scalability issue by using symbolic meaning representations, which can be indexed and retrieved efficiently with complexity that is independent of the text size. We apply our approach, called the N-Gram Machine (NGM), to three representative tasks. First as proof-of-concept, we demonstrate that NGM successfully solves the bAbI tasks of synthetic text. Second, we show that NGM scales to large corpus by experimenting on "life-long bAbI", a special version of bAbI that contains millions of sentences. Lastly on the WikiMovies dataset, we use NGM to induce latent structure (i.e. schema) and answer questions from natural language Wikipedia text, with only QA pairs as weak supervision.


title: Tidal shear and the consistency of microscopic Lagrangian halo approaches

abstract: We delineate the conditions under which the consistency relation for the non-Gaussian bias and the universality of the halo mass function hold in the context of microscopic Lagrangian descriptions of halos. The former is valid provided that the collapse barrier depends only on the physical fields (instead of fields normalized by their variance for example) and explicitly includes the effect of {\it all} physical fields such as the tidal shear. The latter holds provided that the response of the halo number density to a long-wavelength density fluctuation is equivalent to the response induced by shifting the spherical collapse threshold. Our results apply to any Lagrangian halo bias prescription. Effective "moving" barriers, which are ubiquitous in the literature, do not generally satisfy the consistency relation. Microscopic barriers including the tidal shear lead to two additional, second-order Lagrangian bias parameters which ensure that the consistency relation is satisfied. We provide analytic expressions for them.


title: Principal manifold estimation via model complexity selection

abstract: We propose a framework of principal manifolds to model high-dimensional data. This framework is based on Sobolev spaces and designed to model data of any intrinsic dimension. It includes principal component analysis and principal curve algorithm as special cases. We propose a novel method for model complexity selection to avoid overfitting, eliminate the effects of outliers, and improve the computation speed. Additionally, we propose a method for identifying the interiors of circle-like curves and cylinder/ball-like surfaces. The proposed approach is compared to existing methods by simulations and applied to estimate tumor surfaces and interiors in a lung cancer study.


title: Energy dependence of non-local potentials

abstract: Recently a variety of studies have shown the importance of including non-locality in the description of reactions. The goal of this work is to revisit the phenomenological approach to determining non-local optical potentials from elastic scattering. We perform a $\chi^2$ analysis of neutron elastic scattering data off $^{40}$Ca, $^{90}$Zr and $^{208}$Pb at energies $E \approx 5-40$ MeV, assuming a Perey and Buck or Tian, Pang, and Ma non-local form for the optical potential. We introduce energy and asymmetry dependencies in the imaginary part of the potential and refit the data to obtain a global parameterization. Independently of the starting point in the minimization procedure, an energy dependence in the imaginary depth is required for a good description of the data across the included energy range. We present two parameterizations, both of which represent an improvement over the original potentials for the fitted nuclei as well as for other nuclei not included in our fit. Our results show that, even when including the standard Gaussian non-locality in optical potentials, a significant energy dependence is required to describe elastic-scattering data.


title: Methodology study of machine learning for the neutron star equation of state

abstract: We discuss a methodology of machine learning to deduce the neutron star equation of state from a set of mass-radius observational data. We propose an efficient procedure to deal with a mapping from finite data points with observational errors onto an equation of state. We generate training data and optimize the neural network. Using independent validation data (mock observational data) we confirm that the equation of state is correctly reconstructed with precision surpassing observational errors. We finally discuss the relation between our method and Bayesian analysis with an emphasis put on generality of our method for underdetermined problems.


title: On continuous self-maps and homeomorphisms of the Golomb space

abstract: The Golomb space $\mathbb N_\tau$ is the set $\mathbb N$ of positive integers endowed with the topology $\tau$ generated by the base consisting of arithmetic progressions $\{a+bn\}_{n=0}^\infty$ with coprime $a,b$. We prove that the Golomb space $\mathbb N_\tau$ has continuum many continuous self-maps, contains a countable disjoint family of infinite closed connected subsets, the set $\Pi$ of prime numbers is a dense metrizable subspace of $\mathbb N_\tau$, and each homeomorphism $h$ of $\mathbb N_\tau$ has the following properties: $h(1)=1$, $h(\Pi)=\Pi$ and $\Pi_{h(x)}=h(\Pi_x)$ for all $x\in\mathbb N$. Here by $\Pi_x$ we denote the set of prime divisors of $x$.


title: Hyperreflexivity constants of the bounded $n$-cocycle spaces of group algebras and C$^*$-algebras

abstract: We introduced the concept of strong property $(\mathbb{B})$ with a constant for Banach algebras and, by applying certain analysis on the Fourier algebra of a unit circle, we show that all C$^*$-algebras and group algebras have the strong property $(\mathbb{B})$ with a constant given by $288\pi(1+\sqrt{2})$. We then use this result to find a concrete upper bound for the hyperreflexivity constant of $\mathcal{C}^n(A,X)$, the space of bounded $n$-cocycles from $A$ into $X$, where $A$ is a C$^*$-algebra or the group algebra of a group with an open subgroup of polynomial growth and $X$ is a Banach $A$-bimodule for which $\mathcal{H}{n+1}(A,X)$ is a Banach space. As another application, we show that for a locally compact amenable group $G$ and $1<p<\infty$, the space $CV_P(G)$ of convolution operators on $L^p(G)$ are hyperreflexive with a constant given by $288\pi(1+\sqrt{2})$. This is the generalization of a well-known result of E. Christiensen for $p=2$.


title: First principles second harmonic generation of transition metal dichalcogenides and boron nitride alloys: from monolayers and nanotubes to Haeckelites and Schwarzites

abstract: In order to shed light on the second harmonic generation (SHG) of new 2-D systems, first principles methods are used to calculate the second order susceptibility \chi(2) for different types of layered alloys such as monolayers of transition metal dichalcogenide (TMD) alloys, TMD Haeckelite alloys, nanotubes of TMD alloys, hexagonal boron nitride (h-BN) systems which include BxNyCz alloys, BN and BNC2 nanotubes, BxNxCy Haeckelites and BN Schwarzites (porous BN). It is found that the tungsten based alloys possess higher \chi(2) than Mo based at high photon energies, but at low energies, one type of MoSSe dominates. The hypothetical TMD Haeckelites NbSSe and Nb0.5Ta0.5S2 reveal the highest \chi(2) of all the calculated structures. Zigzag TMD alloy nanotubes show higher \chi(2) as the diameter is reduced and approximate to the monolayer for big diameters. BNC alloys exhibit a higher \chi(2) than the h-BN monolayer and are comparable to TMD alloys, except for one case which doubles its intensity. The BN tubes show an increase of \chi(2) as the diameter decreases, similarly to the TMD nanotubes. BxNxCy Haeckelites possess a very high \chi(2) and may shed light on the role of extended defects in nonlinear optical properties. One of the BN Schwarzites exhibits a higher \chi(2) than already known 3-D materials.


title: Political Polarization in Social Media: Analysis of the "Twitter Political Field" in Japan

abstract: There is an ongoing debate about whether the Internet is like a public sphere or an echo chamber. Among many forms of social media, Twitter is one of the most crucial online places for political debate. Most of the previous studies focus on the formal structure of the Twitter political field, such as its homophilic tendency, or otherwise limit the analysis to a few topics. In order to explore whether Twitter functions as an echo chamber in general, however, we have to investigate not only the structure but also the contents of Twitter's political field. Accordingly, we conducted both large-scale social network analysis and natural language processing. We firstly applied a community detection method to the reciprocal following network in Twitter and found five politically distinct communities in the field. We further examined dominant topics discussed therein by employing a topic model in analyzing the content of the tweets, and we found that a topic related to xenophobia is circulated solely in right-wing communities. To our knowledge, this is the first study to address echo chambers in Japanese Twitter political field and to examine the formal structure and the contents of tweets with the combination of large-scale social network analysis and natural language processing.


title: Wing Loss for Robust Facial Landmark Localisation with Convolutional Neural Networks

abstract: We present a new loss function, namely Wing loss, for robust facial landmark localisation with Convolutional Neural Networks (CNNs). We first compare and analyse different loss functions including L2, L1 and smooth L1. The analysis of these loss functions suggests that, for the training of a CNN-based localisation model, more attention should be paid to small and medium range errors. To this end, we design a piece-wise loss function. The new loss amplifies the impact of errors from the interval (-w, w) by switching from L1 loss to a modified logarithm function. To address the problem of under-representation of samples with large out-of-plane head rotations in the training set, we propose a simple but effective boosting strategy, referred to as pose-based data balancing. In particular, we deal with the data imbalance problem by duplicating the minority training samples and perturbing them by injecting random image rotation, bounding box translation and other data augmentation approaches. Last, the proposed approach is extended to create a two-stage framework for robust facial landmark localisation. The experimental results obtained on AFLW and 300W demonstrate the merits of the Wing loss function, and prove the superiority of the proposed method over the state-of-the-art approaches.


title: Algorithms and Data Structures to Accelerate Network Analysis

abstract: As the sheer amount of computer generated data continues to grow exponentially, new bottlenecks are unveiled that require rethinking our traditional software and hardware architectures. In this paper we present five algorithms and data structures (long queue emulation, lockless bimodal queues, tail early dropping, LFN tables, and multiresolution priority queues) designed to optimize the process of analyzing network traffic. We integrated these optimizations on R-Scope, a high performance network appliance that runs the Bro network analyzer, and present benchmarks showcasing performance speed ups of 5X at traffic rates of 10 Gbps.


title: Twisted Orlicz algebras, I

abstract: Let G be a locally compact group, let $\Omega:G\times G\to \mathbb{C}^*$ be a 2-cocycle, and let $\Phi$ be a Young function. In this paper, we consider the Orlicz space $L^\Phi(G)$ and investigate its algebraic property under the twisted convolution $\circledast$ coming from $\Omega$. We find sufficient conditions under which $(L^\Phi(G),\circledast)$ becomes a Banach algebra or a Banach $*$-algebra; we call it a {\it twisted Orlicz algebra}. Furthermore, we study its harmonic analysis properties, such as symmetry, existence of functional calculus, regularity, and having Wiener property, mostly for the case when $G$ is a compactly generated group of polynomial growth. We apply our methods to several important classes of polynomial as well as subexponential weights and demonstrate that our results could be applied to variety of cases.


title: Deep supervised learning using local errors

abstract: Error backpropagation is a highly effective mechanism for learning high-quality hierarchical features in deep networks. Updating the features or weights in one layer, however, requires waiting for the propagation of error signals from higher layers. Learning using delayed and non-local errors makes it hard to reconcile backpropagation with the learning mechanisms observed in biological neural networks as it requires the neurons to maintain a memory of the input long enough until the higher-layer errors arrive. In this paper, we propose an alternative learning mechanism where errors are generated locally in each layer using fixed, random auxiliary classifiers. Lower layers could thus be trained independently of higher layers and training could either proceed layer by layer, or simultaneously in all layers using local error information. We address biological plausibility concerns such as weight symmetry requirements and show that the proposed learning mechanism based on fixed, broad, and random tuning of each neuron to the classification categories outperforms the biologically-motivated feedback alignment learning technique on the MNIST, CIFAR10, and SVHN datasets, approaching the performance of standard backpropagation. Our approach highlights a potential biological mechanism for the supervised, or task-dependent, learning of feature hierarchies. In addition, we show that it is well suited for learning deep networks in custom hardware where it can drastically reduce memory traffic and data communication overheads.


title: Weak amenability of weighted Orlicz algebras

abstract: Let G be a locally compact abelian group, $\omega:G\to (0,\infty)$ be a weight, and ($\Phi$,$\Psi$) be a complementary pair of strictly increasing continuous Young functions. We show that for the weighted Orlicz algebra $L^\Phi_\omega(G)$, the weak amenability is obtained under conditions similar to the one considered by Y. Zhang for weighted group algebras. Our methods can be applied to various families of weighted Orlicz algebras, including weighted $L^p$-spaces.


title: Improving particle filter performance by smoothing observations

abstract: This article shows that increasing the observation variance at small scales can reduce the ensemble size required to avoid collapse in particle filtering of spatially-extended dynamics and improve the resulting uncertainty quantification at large scales. Particle filter weights depend on how well ensemble members agree with observations, and collapse occurs when a few ensemble members receive most of the weight. Collapse causes catastrophic variance underestimation. Increasing small-scale variance in the observation error model reduces the incidence of collapse by de-emphasizing small-scale differences between the ensemble members and the observations. Doing so smooths the posterior mean, though it does not smooth the individual ensemble members. Two options for implementing the proposed observation error model are described. Taking discretized elliptic differential operators as an observation error covariance matrix provides the desired property of a spectrum that grows in the approach to small scales. This choice also introduces structure exploitable by scalable computation techniques, including multigrid solvers and multiresolution approximations to the corresponding integral operator. Alternatively the observations can be smoothed and then assimilated under the assumption of independent errors, which is equivalent to assuming large errors at small scales. The method is demonstrated on a linear stochastic partial differential equation, where it significantly reduces the occurrence of particle filter collapse while maintaining accuracy. It also improves continuous ranked probability scores by as much as 25%, indicating that the weighted ensemble more accurately represents the true distribution. The method is compatible with other techniques for improving the performance of particle filters.


title: Shot noise of charge and spin transport in a junction with a precessing molecular spin

abstract: Magnetic molecules and nanomagnets can be used to influence the electronic transport in mesoscopic junction. In a magnetic field the precessional motion leads to resonances in the dc- and ac-transport properties of a nanocontact, in which the electrons are coupled to the precession. Quantities like the dc-conductance or the ac-response provide valuable information like the level structure and the coupling parameters. Here, we address the current noise properties of such contacts. This encompasses the charge current and spin-torque shot noise, which both show a step-like behavior as functions of bias voltage and magnetic field. The charge current noise shows pronounced dips around the steps, which we trace back to interference effects of electron in quasienergy levels coupled by the molecular spin precession. We show that some components of the noise of the spin-torque currents are directly related to the Gilbert damping and, hence, are experimentally accessible. Our results show that the noise characteristics allow to investigate in more detail the coherence of spin transport in contacts containing magnetic molecules.


title: Backward induction in presence of cycles

abstract: For the classical backward induction algorithm, the input is an arbitrary $n$-person positional game with perfect information modeled by a finite acyclic directed graph (digraph) and the output is a profile $(x_1, \ldots, x_n)$ of pure positional strategies that form some special subgame perfect Nash equilibrium. We extend this algorithm to work with digraphs that may have directed cycles. Each digraph admits a unique partition into strongly connected components, which will be treated as the outcomes of the game. Such a game will be called a {\em deterministic graphical multistage}(DGMS) game. If we identify the outcomes corresponding to all strongly connected components, except terminal positions, we obtain the so-called {\em deterministic graphical}(DG) games, which are frequent in the literature. The outcomes of a DG game are all terminal positions and one special outcome $c$ that is assigned to all infinite plays. We modify the backward induction procedure to adapt it for the DGMS games. However, by doing so, we lose two important properties: the modified algorithm always outputs a {\em Nash equilibrium} (NE) only when $n = 2$ and, even in this case, this NE may be not {\em subgame perfect}. (Yet, in the zero-sum case it is.) The lack of these two properties is not a fault of the algorithm, just (subgame perfect) Nash equilibria in pure positional strategies may fail to exist in the considered game. {\bf Keywords:} deterministic graphical (multistage) game, game in normal and in positional form, saddle point, Nash equilibrium, Nash-solvability, game form, positional structure, directed graph, digraph, directed cycle, acyclic digraph.


title: Scalable Recollections for Continual Lifelong Learning

abstract: Given the recent success of Deep Learning applied to a variety of single tasks, it is natural to consider more human-realistic settings. Perhaps the most difficult of these settings is that of continual lifelong learning, where the model must learn online over a continuous stream of non-stationary data. A successful continual lifelong learning system must have three key capabilities: it must learn and adapt over time, it must not forget what it has learned, and it must be efficient in both training time and memory. Recent techniques have focused their efforts primarily on the first two capabilities while questions of efficiency remain largely unexplored. In this paper, we consider the problem of efficient and effective storage of experiences over very large time-frames. In particular we consider the case where typical experiences are O(n) bits and memories are limited to O(k) bits for k << n. We present a novel scalable architecture and training algorithm in this challenging domain and provide an extensive evaluation of its performance. Our results show that we can achieve considerable gains on top of state-of-the-art methods such as GEM.


title: Multiplicity of self-adjoint realisations of the (2+1)-fermionic model of Ter-Martirosyan-- Skornyakov type

abstract: We reconstruct the whole family of self-adjoint Hamiltonians of Ter-Martirosyan-- Skornyakov type for a system of two identical fermions coupled with a third particle of different nature through an interaction of zero range. We proceed through an operator-theoretic approach based on the self-adjoint extension theory of Krein, Visik, and Birman. We identify the explicit Krein-Visik-Birman extension parameter as an operator on the "space of charges" for this model (the "Krein space") and we come to formulate a sharp conjecture on the dimensionality of its kernel. Based on our conjecture, for which we also discuss an amount of evidence, we explain the emergence of a multiplicity of extensions in a suitable regime of masses and we reproduce for the first time the previous partial constructions obtained by means of an alternative quadratic form approach.


title: Addressing Expensive Multi-objective Games with Postponed Preference Articulation via Memetic Co-evolution

abstract: This paper presents algorithmic and empirical contributions demonstrating that the convergence characteristics of a co-evolutionary approach to tackle Multi-Objective Games (MOGs) with postponed preference articulation can often be hampered due to the possible emergence of the so-called Red Queen effect. Accordingly, it is hypothesized that the convergence characteristics can be significantly improved through the incorporation of memetics (local solution refinements as a form of lifelong learning), as a promising means of mitigating (or at least suppressing) the Red Queen phenomenon by providing a guiding hand to the purely genetic mechanisms of co-evolution. Our practical motivation is to address MOGs of a time-sensitive nature that are characterized by computationally expensive evaluations, wherein there is a natural need to reduce the total number of true function evaluations consumed in achieving good quality solutions. To this end, we propose novel enhancements to co-evolutionary approaches for tackling MOGs, such that memetic local refinements can be efficiently applied on evolved candidate strategies by searching on computationally cheap surrogate payoff landscapes (that preserve postponed preference conditions). The efficacy of the proposal is demonstrated on a suite of test MOGs that have been designed.


title: Image Registration of Very Large Images via Genetic Programming

abstract: Image registration (IR) is a fundamental task in image processing for matching two or more images of the same scene taken at different times, from different viewpoints and/or by different sensors. Due to the enormous diversity of IR applications, automatic IR remains a challenging problem to this day. A wide range of techniques has been developed for various data types and problems. However, they might not handle effectively very large images, which give rise usually to more complex transformations, e.g., deformations and various other distortions. In this paper we present a genetic programming (GP)-based approach for IR, which could offer a significant advantage in dealing with very large images, as it does not make any prior assumptions about the transformation model. Thus, by incorporating certain generic building blocks into the proposed GP framework, we hope to realize a large set of specialized transformations that should yield accurate registration of very large images.


title: A Two-Phase Genetic Algorithm for Image Registration

abstract: Image Registration (IR) is the process of aligning two (or more) images of the same scene taken at different times, different viewpoints and/or by different sensors. It is an important, crucial step in various image analysis tasks where multiple data sources are integrated/fused, in order to extract high-level information. Registration methods usually assume a relevant transformation model for a given problem domain. The goal is to search for the "optimal" instance of the transformation model assumed with respect to a similarity measure in question. In this paper we present a novel genetic algorithm (GA)-based approach for IR. Since GA performs effective search in various optimization problems, it could prove useful also for IR. Indeed, various GAs have been proposed for IR. However, most of them assume certain constraints, which simplify the transformation model, restrict the search space or make additional preprocessing requirements. In contrast, we present a generalized GA-based solution for an almost fully affine transformation model, which achieves competitive results without such limitations using a two-phase method and a multi-objective optimization (MOO) approach. We present good results for multiple dataset and demonstrate the robustness of our method in the presence of noisy data.


title: Genetic Algorithm-Based Solver for Very Large Multiple Jigsaw Puzzles of Unknown Dimensions and Piece Orientation

abstract: In this paper we propose the first genetic algorithm (GA)-based solver for jigsaw puzzles of unknown puzzle dimensions and unknown piece location and orientation. Our solver uses a novel crossover technique, and sets a new state-of-the-art in terms of the puzzle sizes solved and the accuracy obtained. The results are significantly improved, even when compared to previous solvers assuming known puzzle dimensions. Moreover, the solver successfully contends with a mixed bag of multiple puzzle pieces, assembling simultaneously all puzzles.


title: An Automatic Solver for Very Large Jigsaw Puzzles Using Genetic Algorithms

abstract: In this paper we propose the first effective genetic algorithm (GA)-based jigsaw puzzle solver. We introduce a novel crossover procedure that merges two "parent" solutions to an improved "child" configuration by detecting, extracting, and combining correctly assembled puzzle segments. The solver proposed exhibits state-of-the-art performance, as far as handling previously attempted puzzles more accurately and efficiently, as well puzzle sizes that have not been attempted before. The extended experimental results provided in this paper include, among others, a thorough inspection of up to 30,745-piece puzzles (compared to previous attempts on 22,755-piece puzzles), using a considerably faster concurrent implementation of the algorithm. Furthermore, we explore the impact of different phases of the novel crossover operator by experimenting with several variants of the GA. Finally, we compare different fitness functions and their effect on the overall results of the GA-based solver.


title: A Generalized Genetic Algorithm-Based Solver for Very Large Jigsaw Puzzles of Complex Types

abstract: In this paper we introduce new types of square-piece jigsaw puzzles, where in addition to the unknown location and orientation of each piece, a piece might also need to be flipped. These puzzles, which are associated with a number of real world problems, are considerably harder, from a computational standpoint. Specifically, we present a novel generalized genetic algorithm (GA)-based solver that can handle puzzle pieces of unknown location and orientation (Type 2 puzzles) and (two-sided) puzzle pieces of unknown location, orientation, and face (Type 4 puzzles). To the best of our knowledge, our solver provides a new state-of-the-art, solving previously attempted puzzles faster and far more accurately, handling puzzle sizes that have never been attempted before, and assembling the newly introduced two-sided puzzles automatically and effectively. This paper also presents, among other results, the most extensive set of experimental results, compiled as of yet, on Type 2 puzzles.


title: A Genetic Algorithm-Based Solver for Very Large Jigsaw Puzzles

abstract: In this paper we propose the first effective automated, genetic algorithm (GA)-based jigsaw puzzle solver. We introduce a novel procedure of merging two "parent" solutions to an improved "child" solution by detecting, extracting, and combining correctly assembled puzzle segments. The solver proposed exhibits state-of-the-art performance solving previously attempted puzzles faster and far more accurately, and also puzzles of size never before attempted. Other contributions include the creation of a benchmark of large images, previously unavailable. We share the data sets and all of our results for future testing and comparative evaluation of jigsaw puzzle solvers.


title: The ALMA Phasing System: A Beamforming Capability for Ultra-High-Resolution Science at (Sub)Millimeter Wavelengths

abstract: The Atacama Millimeter/submillimeter Array (ALMA) Phasing Project (APP) has developed and deployed the hardware and software necessary to coherently sum the signals of individual ALMA antennas and record the aggregate sum in Very Long Baseline Interferometry (VLBI) Data Exchange Format. These beamforming capabilities allow the ALMA array to collectively function as the equivalent of a single large aperture and participate in global VLBI arrays. The inclusion of phased ALMA in current VLBI networks operating at (sub)millimeter wavelengths provides an order of magnitude improvement in sensitivity, as well as enhancements in u-v coverage and north-south angular resolution. The availability of a phased ALMA enables a wide range of new ultra-high angular resolution science applications, including the resolution of supermassive black holes on event horizon scales and studies of the launch and collimation of astrophysical jets. It also provides a high-sensitivity aperture that may be used for investigations such as pulsar searches at high frequencies. This paper provides an overview of the ALMA Phasing System design, implementation, and performance characteristics.


title: Approximate Gradient Coding via Sparse Random Graphs

abstract: Distributed algorithms are often beset by the straggler effect, where the slowest compute nodes in the system dictate the overall running time. Coding-theoretic techniques have been recently proposed to mitigate stragglers via algorithmic redundancy. Prior work in coded computation and gradient coding has mainly focused on exact recovery of the desired output. However, slightly inexact solutions can be acceptable in applications that are robust to noise, such as model training via gradient-based algorithms. In this work, we present computationally simple gradient codes based on sparse graphs that guarantee fast and approximately accurate distributed computation. We demonstrate that sacrificing a small amount of accuracy can significantly increase algorithmic robustness to stragglers.

categories: stat.ML cs.DC cs.IT cs.LG math.IT stat.CO


title: Separable discrete functions: recognition and sufficient conditions

abstract: A discrete function of $n$ variables is a mapping $g : X_1 \times \ldots \times X_n \rightarrow A$, where $X_1, \ldots, X_n$, and $A$ are arbitrary finite sets. Function $g$ is called {\em separable} if there exist $n$ functions $g_i : X_i \rightarrow A$ for $i = 1, \ldots, n$, such that for every input $x_1, \ldots ,x_n$ the function $g(x_1, \ldots, x_n)$ takes one of the values $g_1(x_1), \ldots ,g_n(x_n)$. Given a discrete function $g$, it is an interesting problem to ask whether $g$ is separable or not. Although this seems to be a very basic problem concerning discrete functions, the complexity of recognition of separable discrete functions of $n$ variables is known only for $n=2$. In this paper we will show that a slightly more general recognition problem, when $g$ is not fully but only partially defined, is NP-complete for $n \geq 3$. We will then use this result to show that the recognition of fully defined separable discrete functions is NP-complete for $n \geq 4$. The case $n = 2$ is well-studied in the context of game theory, where (separable) discrete functions of $n$ variables are referred to as (assignable) $n$-person game forms. There is a known sufficient condition for assignability (separability) of two-person game forms (discrete functions of two variables) called (weak) total tightness of a game form. This property can be tested in polynomial time, and can be easily generalized both to higher dimension and to partially defined functions. We will prove in this paper that weak total tightness implies separability for (partially defined) discrete functions of $n$ variables for any $n$, thus generalizing the above result known for $n=2$. Keywords: separable discrete functions, totally tight and assignable game forms


title: Internal wave energy flux from density perturbations in nonlinear stratifications

abstract: Internal gravity wave energy contributes significantly to the energy budget of the oceans, affecting mixing and the thermohaline circulation. Hence it is important to determine the internal wave energy flux $\mathbf{J} = p \, \mathbf{v}$, where $p$ is the pressure perturbation field and $\mathbf{v}$ is the velocity perturbation field. However, the pressure perturbation field is not directly accessible in laboratory or field observations. Previously, a Green's function based method was developed to calculate the instantaneous energy flux field from a measured density perturbation field $\rho(x,z,t)$, given a $constant$ buoyancy frequency $N$. Here we present methods for computing the instantaneous energy flux $\mathbf{J}(x,z,t)$ for a spatially varying $N(z)$, as in the oceans where $N(z)$ typically decreases by two orders of magnitude from shallow water to the deep ocean. Analytic methods are presented for computing $\mathbf{J}(x,z,t)$ from a density perturbation field for $N(z)$ varying linearly with $z$ and for $N^2(z)$ varying as $\tanh (z)$. To generalize this approach to $arbitrary$ $N(z)$, we present a computational method for obtaining $\mathbf{J}(x,z,t)$. The results for $\mathbf{J}(x,z,t)$ for the different cases agree well with results from direct numerical simulations of the Navier-Stokes equations. Our computational method can be applied to any density perturbation data using the MATLAB graphical user interface EnergyFlux.


title: Designing Coalition-Proof Reverse Auctions over Continuous Goods

abstract: This paper investigates reverse auctions that involve continuous values of different types of goods, general nonconvex constraints, and second stage costs. We seek to design the payment rules and conditions under which coalitions of participants cannot influence the auction outcome in order to obtain higher collective utility. Under the incentive-compatible Vickrey-Clarke-Groves mechanism, we show that coalition-proof outcomes are achieved if the submitted bids are convex and the constraint sets are of a polymatroid-type. These conditions, however, do not capture the complexity of the general class of reverse auctions under consideration. By relaxing the property of incentive-compatibility, we investigate further payment rules that are coalition-proof without any extra conditions on the submitted bids and the constraint sets. Since calculating the payments directly for these mechanisms is computationally difficult for auctions involving many participants, we present two computationally efficient methods. Our results are verified with several case studies based on electricity market data.


title: Grain boundary phases in bcc metals

abstract: We report a computational discovery of novel grain boundary structures and multiple grain boundary phases in elemental bcc tungsten. While grain boundary structures created by the \gamma-surface method as a union of two perfect half crystals have been studied extensively, it is known that the method has limitations and does not always predict the correct ground states. Here, we use a newly developed computational tool, based on evolutionary algorithms, to perform a grand-canonical search of a high-angle symmetric tilt boundary in tungsten, and we find new ground states and multiple phases that cannot be described using the conventional structural unit model. We use MD simulations to demonstrate that the new structures can coexist at finite temperature in a closed system, confirming these are examples of different GB phases. The new ground state is confirmed by first-principles calculations.


title: Generalized surface tension bounds in vacuum decay

abstract: Coleman and De Luccia (CDL) showed that gravitational effects can prevent the decay by bubble nucleation of a Minkowski or AdS false vacuum. In their thin-wall approximation this happens whenever the surface tension in the bubble wall exceeds an upper bound proportional to the difference of the square roots of the true and false vacuum energy densities. Recently it was shown that there is another type of thin-wall regime that differs from that of CDL in that the radius of curvature grows substantially as one moves through the wall. Not only does the CDL derivation of the bound fail in this case, but also its very formulation becomes ambiguous because the surface tension is not well-defined. We propose a definition of the surface tension and show that it obeys a bound similar in form to that of the CDL case. We then show that both thin-wall bounds are special cases of a more general bound that is satisfied for all bounce solutions with Minkowski or AdS false vacua. We discuss the limit where the parameters of the theory attain critical values and the bound is saturated. The bounce solution then disappears and a static planar domain wall solution appears in its stead. The scalar field potential then is of the form expected in supergravity, but this is only guaranteed along the trajectory in field space traced out by the bounce.


title: Stability of Fixed Points and Chaos in Fractional Systems

abstract: In this paper we propose a method to define the range of stability of fixed points for a variety of discrete fractional systems of the order $0 < \alpha <2$. The method is tested on various forms of fractional generalizations of the standard and logistic maps. Based on our analysis we make a conjecture that chaos is impossible in the corresponding continuous fractional systems.


title: Excitation Backprop for RNNs

abstract: Deep models are state-of-the-art for many vision tasks including video action recognition and video captioning. Models are trained to caption or classify activity in videos, but little is known about the evidence used to make such decisions. Grounding decisions made by deep networks has been studied in spatial visual content, giving more insight into model predictions for images. However, such studies are relatively lacking for models of spatiotemporal visual content - videos. In this work, we devise a formulation that simultaneously grounds evidence in space and time, in a single pass, using top-down saliency. We visualize the spatiotemporal cues that contribute to a deep model's classification/captioning output using the model's internal representation. Based on these spatiotemporal cues, we are able to localize segments within a video that correspond with a specific action, or phrase from a caption, without explicitly optimizing/training for these tasks.


title: Machine Learning Approaches for Traffic Volume Forecasting: A Case Study of the Moroccan Highway Network

abstract: In this paper, we aim to illustrate different approaches we followed while developing a forecasting tool for highway traffic in Morocco. Two main approaches were adopted: Statistical Analysis as a step of data exploration and data wrangling. Therefore, a beta model is carried out for a better understanding of traffic behavior. Next, we moved to Machine Learning where we worked with a bunch of algorithms such as Random Forest, Artificial Neural Networks, Extra Trees, etc. yet, we were convinced that this field of study is still considered under state of the art models, so, we were also covering an application of Long Short-Term Memory Neural Networks.


title: Local melting attracts grain boundaries in colloidal polycrystals

abstract: We find that laser-induced local melting attracts and deforms grain boundaries in 2D colloidal crystals. When a melted region in contact with the edge of a crystal grain recrystallizes, it deforms the grain boundary --- this attraction is driven by the multiplicity of deformed grain boundary configurations. Furthermore, the attraction provides a method to fabricate artificial colloidal crystal grains of arbitrary shape, enabling new experimental studies of grain boundary dynamics and ultimately hinting at a novel approach for fabricating materials with designer microstructures.


title: Operation of a LAr-TPC equipped with a multilayer LEM charge readout

abstract: A novel detector for the ionization signal in a single phase LAr-TPC, based on the adoption of a multilayer Large Electron Multiplier (LEM) replacing the traditional anodic wire arrays, has been experimented in the ICARINO test facility at the INFN Laboratories in Legnaro. Cosmic muon tracks were detected allowing the measurement of energy deposition and a first determination of the signal to noise ratio. The analysis of the recorded events demonstrated the 3D reconstruction capability of ionizing events in this device in liquid Argon, collecting a fraction of about 90% of the ionization signal with signal to noise ratio similar to that measured with more traditional wire chambers


title: Leave no Trace: Learning to Reset for Safe and Autonomous Reinforcement Learning

abstract: Deep reinforcement learning algorithms can learn complex behavioral skills, but real-world application of these methods requires a large amount of experience to be collected by the agent. In practical settings, such as robotics, this involves repeatedly attempting a task, resetting the environment between each attempt. However, not all tasks are easily or automatically reversible. In practice, this learning process requires extensive human intervention. In this work, we propose an autonomous method for safe and efficient reinforcement learning that simultaneously learns a forward and reset policy, with the reset policy resetting the environment for a subsequent attempt. By learning a value function for the reset policy, we can automatically determine when the forward policy is about to enter a non-reversible state, providing for uncertainty-aware safety aborts. Our experiments illustrate that proper use of the reset policy can greatly reduce the number of manual resets required to learn a task, can reduce the number of unsafe actions that lead to non-reversible states, and can automatically induce a curriculum.


title: Exact alignment recovery for correlated Erd\H{o}s-R\'enyi graphs

abstract: We consider the problem of perfectly recovering the vertex correspondence between two correlated Erd\H{o}s-R\'enyi (ER) graphs on the same vertex set. The correspondence between the vertices can be obscured by randomly permuting the vertex labels of one of the graphs. We determine the information-theoretic threshold for exact recovery, i.e. the conditions under which the entire vertex correspondence can be correctly recovered given unbounded computational resources.


title: Inducing skyrmions in ultrathin Fe films by hydrogen exposure

abstract: Magnetic skyrmions are localized nanometer-sized spin configurations with particle_like properties, which are envisioned to be used as bits in next_generation information technology. An essential step towards future skyrmion-based applications is to engineer key magnetic parameters for developing and stabilizing individual magnetic skyrmions. Here we demonstrate the tuning of the non_collinear magnetic state of an Fe double layer on an Ir111 substrate by loading the sample with atomic hydrogen. By using spin_polarized scanning tunneling microscopy, we discover that the hydrogenated system supports the formation of skyrmions in external magnetic fields, while the pristine Fe double layer does not. Based on ab initio calculations, we attribute this effect to the tuning of the Heisenberg exchange and the Dzyaloshinsky_Moriya interactions due to hydrogenation. In addition to interface engineering, hydrogenation of thin magnetic films offers a unique pathway to design and optimize the skyrmionic states in low_dimensional magnetic materials.


title: New convergence analysis of a primal-dual algorithm with large stepsizes

abstract: We consider a primal-dual algorithm for minimizing $f(x)+h\square l(Ax)$ with Fr\'echet differentiable $f$ and $l^*$. This primal-dual algorithm has two names in literature: Primal-Dual Fixed-Point algorithm based on the Proximity Operator (PDFP$^2$O) and Proximal Alternating Predictor-Corrector (PAPC). In this paper, we prove its convergence under a weaker condition on the stepsizes than existing ones. With additional assumptions, we show its linear convergence. In addition, we show that this condition (the upper bound of the stepsize) is tight and can not be weakened. This result also recovers a recently proposed positive-indefinite linearized augmented Lagrangian method. In addition, we apply this result to a decentralized consensus algorithm PG-EXTRA and derive the weakest convergence condition.


title: Measuring Territorial Control in Civil Wars Using Hidden Markov Models: A Data Informatics-Based Approach

abstract: Territorial control is a key aspect shaping the dynamics of civil war. Despite its importance, we lack data on territorial control that are fine-grained enough to account for subnational spatio-temporal variation and that cover a large set of conflicts. To resolve this issue, we propose a theoretical model of the relationship between territorial control and tactical choice in civil war and outline how Hidden Markov Models (HMMs) are suitable to capture theoretical intuitions and estimate levels of territorial control. We discuss challenges of using HMMs in this application and mitigation strategies for future work.


title: Learning Aggregated Transmission Propagation Networks for Haze Removal and Beyond

abstract: Single image dehazing is an important low-level vision task with many applications. Early researches have investigated different kinds of visual priors to address this problem. However, they may fail when their assumptions are not valid on specific images. Recent deep networks also achieve relatively good performance in this task. But unfortunately, due to the disappreciation of rich physical rules in hazes, large amounts of data are required for their training. More importantly, they may still fail when there exist completely different haze distributions in testing images. By considering the collaborations of these two perspectives, this paper designs a novel residual architecture to aggregate both prior (i.e., domain knowledge) and data (i.e., haze distribution) information to propagate transmissions for scene radiance estimation. We further present a variational energy based perspective to investigate the intrinsic propagation behavior of our aggregated deep model. In this way, we actually bridge the gap between prior driven models and data driven networks and leverage advantages but avoid limitations of previous dehazing approaches. A lightweight learning framework is proposed to train our propagation network. Finally, by introducing a taskaware image separation formulation with a flexible optimization scheme, we extend the proposed model for more challenging vision tasks, such as underwater image enhancement and single image rain removal. Experiments on both synthetic and realworld images demonstrate the effectiveness and efficiency of the proposed framework.


title: MinimalRNN: Toward More Interpretable and Trainable Recurrent Neural Networks

abstract: We introduce MinimalRNN, a new recurrent neural network architecture that achieves comparable performance as the popular gated RNNs with a simplified structure. It employs minimal updates within RNN, which not only leads to efficient learning and testing but more importantly better interpretability and trainability. We demonstrate that by endorsing the more restrictive update rule, MinimalRNN learns disentangled RNN states. We further examine the learning dynamics of different RNN structures using input-output Jacobians, and show that MinimalRNN is able to capture longer range dependencies than existing RNN architectures.


title: Removing instability of Polonyi-Starobinsky supergravity by adding FI term

abstract: Polonyi-Starobinsky (PS) supergravity is the N=1 supergravity model of Starobinsky inflation with spontaneous supersymmetry breaking (after inflation) due to Polonyi superfield, and inflaton belonging to a massive vector supermultiplet. The PS model is used for an explicit realization of the (super-heavy) gravitino dark matter scenario in cosmology. We find a potential instability in this model, and offer a mechanism for its removal by adding a Fayet-Iliopoulos (FI) term.


title: Mitigating Read-disturbance Errors in STT-RAM Caches by Using Data Compression

abstract: Due to its high density and close-to-SRAM read latency, spin transfer torque RAM (STT-RAM) is considered one of the most-promising emerging memory technologies for designing large last level caches (LLCs). However, in deep sub-micron region, STT-RAM shows read-disturbance error (RDE) whereby a read operation may modify the stored data value and this presents a severe threat to performance and reliability of STT-RAM caches. In this paper, we present a technique, named SHIELD, to mitigate RDE in STT-RAM LLCs. SHIELD uses data compression to reduce number of read operations from STT-RAM blocks to avoid RDE and also to reduce the number of bits written to cache during both write and restore operations. Experimental results have shown that SHIELD provides significant improvement in performance and energy efficiency. SHIELD consumes smaller energy than two previous RDE-mitigation techniques, namely high-current restore required read (HCRR, also called restore-after-read) and low-current long latency read (LCLL) and even an ideal RDE-free STT-RAM cache.


title: Faraday instability and subthreshold Faraday waves: surface waves emitted by walkers

abstract: A walker is a fluid entity comprising a bouncing droplet coupled to the waves that it generates at the surface of a vibrated bath. Thanks to this coupling, walkers exhibit a series of wave-particle features formerly thought to be exclusive to the quantum realm. In this paper, we derive a model of the Faraday surface waves generated by an impact upon a vertically vibrated liquid surface. We then particularise this theoretical framework to the case of forcing slightly below the Faraday instability threshold. Among others, this theory yields a rationale for the dependence of the wave amplitude to the phase of impact, as well as the characteristic timescale and length scale of viscous damping. The theory is validated with experiments of bead impact on a vibrated bath. We finally discuss implications of these results for the analogy between walkers and quantum particles.


title: Enumeration of Some Closed Knight Paths

abstract: The aim of the paper is to enumerate all closed knight paths of length n over a square board of size n+1. The closed knight paths of length 4, 6 and 8 are classified up to equivalence. We determine that there are exactly 3 equivalence classes of closed knight paths of length 4, exactly 25 equivalence classes of closed knight paths of length 6 and exactly 478 equivalence classes of closed knight paths of length 8.


title: Tree-Structured Boosting: Connections Between Gradient Boosted Stumps and Full Decision Trees

abstract: Additive models, such as produced by gradient boosting, and full interaction models, such as classification and regression trees (CART), are widely used algorithms that have been investigated largely in isolation. We show that these models exist along a spectrum, revealing never-before-known connections between these two approaches. This paper introduces a novel technique called tree-structured boosting for creating a single decision tree, and shows that this method can produce models equivalent to CART or gradient boosted stumps at the extremes by varying a single parameter. Although tree-structured boosting is designed primarily to provide both the model interpretability and predictive performance needed for high-stake applications like medicine, it also can produce decision trees represented by hybrid models between CART and boosted stumps that can outperform either of these approaches.


title: Co-attending Free-form Regions and Detections with Multi-modal Multiplicative Feature Embedding for Visual Question Answering

abstract: Recently, the Visual Question Answering (VQA) task has gained increasing attention in artificial intelligence. Existing VQA methods mainly adopt the visual attention mechanism to associate the input question with corresponding image regions for effective question answering. The free-form region based and the detection-based visual attention mechanisms are mostly investigated, with the former ones attending free-form image regions and the latter ones attending pre-specified detection-box regions. We argue that the two attention mechanisms are able to provide complementary information and should be effectively integrated to better solve the VQA problem. In this paper, we propose a novel deep neural network for VQA that integrates both attention mechanisms. Our proposed framework effectively fuses features from free-form image regions, detection boxes, and question representations via a multi-modal multiplicative feature embedding scheme to jointly attend question-related free-form image regions and detection boxes for more accurate question answering. The proposed method is extensively evaluated on two publicly available datasets, COCO-QA and VQA, and outperforms state-of-the-art approaches. Source code is available at https://github.com/lupantech/dual-mfa-vqa.


title: Prediction Scores as a Window into Classifier Behavior

abstract: Most multi-class classifiers make their prediction for a test sample by scoring the classes and selecting the one with the highest score. Analyzing these prediction scores is useful to understand the classifier behavior and to assess its reliability. We present an interactive visualization that facilitates per-class analysis of these scores. Our system, called Classilist, enables relating these scores to the classification correctness and to the underlying samples and their features. We illustrate how such analysis reveals varying behavior of different classifiers. Classilist is available for use online, along with source code, video tutorials, and plugins for R, RapidMiner, and KNIME at https://katehara.github.io/classilist-site/.


title: Thermodynamic formalism and the entropy at infinity of the geodesic flow

abstract: In this paper we study the ergodic theory and thermodynamic formalism of the geodesic flow on non-compact pinched negatively curved manifolds. We consider two notions of entropy at infinity, the topological and the measure theoretic entropy at infinity. We prove that both notions coincide and we relate this quantity with the upper semicontinuity of the entropy map. This relationship was already studied in [RV] for geometrically finite manifolds; in this paper we extend such results to arbitrary pinched negatively curved manifolds. We obtain several applications to the existence and stability of equilibrium states. Our approach has the advantage to include certain uniformly continous potentials, a situation that can not be studied through the methods from [PPS]. Of particular importance are potentials that vanish at infinity and those that satisfy a critical gap condition, that is, strongly positive recurrent potentials. We also prove some equidistribution and counting results for closed geodesics (beyond the geometrically finite case). In particular we obtain a version of the prime geodesic theorem for strongly positive recurrent manifolds.


title: Short proofs for generalizations of the Lov\'asz Local Lemma: Shearer's condition and cluster expansion

abstract: The Lov\'asz Local Lemma is a seminal result in probabilistic combinatorics. It gives a sufficient condition on a probability space and a collection of events for the existence of an outcome that simultaneously avoids all of those events. Over the years, more general conditions have been discovered under which the conclusion of the lemma continues to hold. In this note we provide short proofs of two of those more general results: Shearer's lemma and the cluster expansion lemma, in their "lopsided" form. We conclude by using the cluster expansion lemma to prove that the symmetric form of the local lemma holds with probabilities bounded by $1/ed$, rather than the bound $1/e(d+1)$ required by the traditional proofs.


title: MorphNet: Fast & Simple Resource-Constrained Structure Learning of Deep Networks

abstract: We present MorphNet, an approach to automate the design of neural network structures. MorphNet iteratively shrinks and expands a network, shrinking via a resource-weighted sparsifying regularizer on activations and expanding via a uniform multiplicative factor on all layers. In contrast to previous approaches, our method is scalable to large networks, adaptable to specific resource constraints (e.g. the number of floating-point operations per inference), and capable of increasing the network's performance. When applied to standard network architectures on a wide variety of datasets, our approach discovers novel structures in each domain, obtaining higher performance while respecting the resource constraint.


title: Relieving the Tension between Weak Lensing and Cosmic Microwave Background with Interacting Dark Matter and Dark Energy Models

abstract: We constrain interacting dark matter and dark energy (IDMDE) models using a 450-degree-square cosmic shear data from the Kilo Degree Survey (KiDS) and the angular power spectra from Planck's latest cosmic microwave background measurements. We revisit the discordance problem in the standard Lambda cold dark matter ($\Lambda$CDM) model between weak lensing and Planck datasets and extend the discussion by introducing interacting dark sectors. The IDMDE models are found to be able to alleviate the discordance between KiDS and Planck as previously inferred from the $\Lambda$CDM model, and moderately favored by a combination of the two datasets.


title: Scalable Relaxations of Sparse Packing Constraints: Optimal Biocontrol in Predator-Prey Network

abstract: Cascades represent rapid changes in networks. A cascading phenomenon of ecological and economic impact is the spread of invasive species in geographic landscapes. The most promising management strategy is often biocontrol, which entails introducing a natural predator able to control the invading population, a setting that can be treated as two interacting cascades of predator and prey populations. We formulate and study a nonlinear problem of optimal biocontrol: optimally seeding the predator cascade over time to minimize the harmful prey population. Recurring budgets, which typically face conservation organizations, naturally leads to sparse constraints which make the problem amenable to approximation algorithms. Available methods based on continuous relaxations scale poorly, to remedy this we develop a novel and scalable randomized algorithm based on a width relaxation, applicable to a broad class of combinatorial optimization problems. We evaluate our contributions in the context of biocontrol for the insect pest Hemlock Wolly Adelgid (HWA) in eastern North America. Our algorithm outperforms competing methods in terms of scalability and solution quality, and finds near optimal strategies for the control of the HWA for fine-grained networks -- an important problem in computational sustainability.


title: Whole planet coupling between climate, mantle, and core: Implications for the evolution of rocky planets

abstract: Earth's climate, mantle, and core interact over geologic timescales. Climate influences whether plate tectonics can take place on a planet, with cool climates being favorable for plate tectonics because they enhance stresses in the lithosphere, suppress plate boundary annealing, and promote hydration and weakening of the lithosphere. Plate tectonics plays a vital role in the long-term carbon cycle, which helps to maintain a temperate climate. Plate tectonics provides long-term cooling of the core, which is vital for generating a magnetic field, and the magnetic field is capable of shielding atmospheric volatiles from the solar wind. Coupling between climate, mantle, and core can potentially explain the divergent evolution of Earth and Venus. As Venus lies too close to the sun for liquid water to exist, there is no long-term carbon cycle and thus an extremely hot climate. Therefore plate tectonics cannot operate and a long-lived core dynamo cannot be sustained due to insufficient core cooling. On planets within the habitable zone where liquid water is possible, a wide range of evolutionary scenarios can take place depending on initial atmospheric composition, bulk volatile content, or the timing of when plate tectonics initiates, among other factors. Many of these evolutionary trajectories would render the planet uninhabitable. However, there is still significant uncertainty over the nature of the coupling between climate, mantle, and core. Future work is needed to constrain potential evolutionary scenarios and the likelihood of an Earth-like evolution.


title: Effects of magnetic dopants in (Li$_{0.8}$M$_{0.2}$OH)FeSe (M = Fe, Mn, Co): a density-functional theory study using band unfolding technique

abstract: The effects of Fe dopants in (Li$_{0.8}$Fe$_{0.2}$OH)FeSe on the electronic band structure are investigated by band unfolding ($k$-projection) technique based on first-principles supercell calculations. Doping 20\% Fe into the LiOH layers has significant effects on the band structure, that is, the Fe impurities doping electrons to the FeSe layers not only shift the Fermi level, but also induce substantial changes in the profile of bands around the Fermi level. However, the magnetic ordering in the dopants has minor effects on the band structure due to the fact that there is only a weak bonding between the LiOH and FeSe layers. Electronic bands for the surface FeSe layer show noticeable differences from those for inner layers in both the location of the Fermi level and details of the bands near the high symmetry points. The band structure for the surface FeSe layer where the Fe atoms are in checkerboard antiferromagnetic order is consistent with angle-resolved photoemission spectroscopy results. Mn and Co have similar doping effects on the band structure of (LiOH)FeSe.


title: Reduction of total-cost and average-cost MDPs with weakly continuous transition probabilities to discounted MDPs

abstract: This note describes sufficient conditions under which total-cost and average-cost Markov decision processes (MDPs) with general state and action spaces, and with weakly continuous transition probabilities, can be reduced to discounted MDPs. For undiscounted problems, these reductions imply the validity of optimality equations and the existence of stationary optimal policies. The reductions also provide methods for computing optimal policies. The results are applied to a capacitated inventory control problem with fixed costs and lost sales.