title: VIS3COS: I. survey overview and the role of environment and stellar mass on star formation

abstract: We present the VIMOS Spectroscopic Survey of a Supercluster in the COSMOS field (VIS3COS) at z~0.84. We use VIMOS high-resolution spectra (GG475 filter) to spectroscopically select 490 galaxies in and around the super-structure and an additional 481 galaxies in the line of sight. We present the redshift distribution, the catalog to be made public and the first results on the properties of individual galaxies and stacked spectra (3500{\AA}<\lambda<4200{\AA} rest-frame). We probe a wide range of densities and environments (from low-density field to clusters and rich groups). We find a decrease in the median star formation rate from low to high density environments in all bins of stellar mass and a sharp rise of the quenched fraction (from ~10% to ~40-60%) of intermediate stellar mass galaxies $(10<\log_{10}\left(M_\star/\mathrm{M_\odot}\right)<10.75)$ from filaments to clusters. The quenched fraction for massive galaxies shows little dependence on the environment being constant at ~30-40%. We find a break in the strength of the [OII] emission, with nearly constant line equivalent widths at lower densities (~-11{\AA}) and then a drop to ~-2.5{\AA} towards higher densities. The break in the [OII] line strength happens at similar densities ($\log_{10}(1+\delta)\sim0.0-0.5$) as the observed rise in the quenched fraction. Our results may provide further clues regarding the different environmental processes affecting galaxies with different stellar masses and highlight the advantages of a single dataset in the COSMOS field probing a wide range of stellar masses and environments. We hypothesize that quenching mechanisms are enhanced in high-density regions.


title: The Type I Superluminous Supernova PS16aqv: Lightcurve Complexity and Deep Limits on Radioactive Ejecta in a Fast Event

abstract: [Abridged] We present UV/optical observations of PS16aqv (SN 2016ard), a Type I superluminous supernova (SLSN-I) classified as part of our search for low-$z$ SLSNe. PS16aqv is a fast evolving SLSNe-I that reached a peak absolute magnitude of $M_{r} \approx -22.1$. The lightcurves exhibit a significant undulation at 30 rest-frame days after peak, with a behavior similar to undulations seen in the slowly fading SLSN-I SN 2015bn. This similarity strengthens the case that fast and slow SLSNe-I form a continuum with a common origin. At $\approx\!80$ days after peak, the lightcurves exhibit a transition to a slow decline, followed by significant subsequent steepening, indicative of a plateau phase or a second significant undulation. Deep limits at $\approx280$ days after peak imply a tight constraint on the nickel mass, $M_{\rm Ni} \lesssim 0.35$ M$_{\odot}$ (lower than for previous SLSNe-I), and indicate that some SLSNe-I do not produce significantly more nickel than normal Type Ic SNe. Using MOSFiT, we model the lightcurve with a magnetar central engine model and find $P_{\rm spin} \approx 0.9$ ms, $B \approx 1.5 \times 10^{14}$ G, and $M_{\rm ej} \approx 16$ M$_{\odot}$. The implied rapid spin-down time and large reservoir of available energy coupled with the high ejecta mass may account for the fast evolving lightcurve and slow spectroscopic evolution. We also study the location of PS16aqv in its host galaxy and find that it occurred at an offset of $2.46 \pm 0.21$ kpc from the central star-forming region. We find the host galaxy exhibits low metallicity and spatially varying extinction and star formation rate, with the explosion site exhibiting lower values than the central region. The complexity seen in the lightcurves of PS16aqv and other events highlights the importance of obtaining well-sampled lightcurves for exploring deviations from a uniform decline.


title: Approximate Data Depth Revisited

abstract: Halfspace depth and $\beta$-skeleton depth are two types of depth functions in nonparametric data analysis. The halfspace depth of a query point $q\in \mathbb{R}^d$ with respect to $S\subset\mathbb{R}^d$ is the minimum portion of the elements of $S$ which are contained in a halfspace which passes through $q$. For $\beta \geq 1$, the $\beta$-skeleton depth of $q$ with respect to $S$ is defined to be the total number of \emph{$\beta$-skeleton influence regions} that contain $q$, where each of these influence regions is the intersection of two hyperballs obtained from a pair of points in $S$. The $\beta$-skeleton depth introduces a family of depth functions that contain \emph{spherical depth} and \emph{lens depth} if $\beta=1$ and $\beta=2$, respectively. The main results of this paper include approximating the planar halfspace depth and $\beta$-skeleton depth using two different approximation methods. First, the halfspace depth is approximated by the $\beta$-skeleton depth values. For this method, two dissimilarity measures based on the concepts of \emph{fitting function} and \emph{Hamming distance} are defined to train the halfspace depth function by the $\beta$-skeleton depth values obtaining from a given data set. The goodness of this approximation is measured by a function of error values. Secondly, computing the planar $\beta$-skeleton depth is reduced to a combination of some range counting problems. Using existing results on range counting approximations, the planar $\beta$-skeleton depth of a query point is approximated in $O(n\;poly(1/\varepsilon,\log n))$, $\beta\geq 1$. Regarding the $\beta$-skeleton depth functions, it is also proved that this family of depth functions converge when $\beta \to \infty$. Finally, some experimental results are provided to support the proposed method of approximation and convergence of $\beta$-skeleton depth functions.


title: A combination theorem for Anosov subgroups

abstract: We prove an analogue of Klein combination theorem for Anosov subgroups by using a local-to-global principle for Morse quasigeodesics.


title: Testing Alignment of Node Attributes with Network Structure Through Label Propagation

abstract: Attributed network data is becoming increasingly common across fields, as we are often equipped with information about nodes in addition to their pairwise connectivity patterns. This extra information can manifest as a classification, or as a multidimensional vector of features. Recently developed methods that seek to extend community detection approaches to attributed networks have explored how to most effectively combine connectivity and attribute information to identify quality communities. These methods often rely on some assumption of the dependency relationships between attributes and connectivity. In this work, we seek to develop a statistical test to assess whether node attributes align with network connectivity. The objective is to quantitatively evaluate whether nodes with similar connectivity patterns also have similar attributes. To address this problem, we use a node sampling and label propagation approach. We apply our method to several synthetic examples that explore how network structure and attribute characteristics affect the empirical p-value computed by our method. Finally, we apply the test to a network generated from a single cell mass cytometry (CyTOF) dataset and show that our test can identify markers associated with distinct sub populations of single cells.


title: Algorithms for Estimating Trends in Global Temperature Volatility

abstract: Trends in terrestrial temperature variability are perhaps more relevant for species viability than trends in mean temperature. In this paper, we develop methodology for estimating such trends using multi-resolution climate data from polar orbiting weather satellites. We derive two novel algorithms for computation that are tailored for dense, gridded observations over both space and time. We evaluate our methods with a simulation that mimics these data's features and on a large, publicly available, global temperature dataset with the eventual goal of tracking trends in cloud reflectance temperature variability.


title: Analytical and numerical study of backreacting one-dimensional holographic superconductors in the presence of Born-Infeld electrodynamics

abstract: We analytically as well as numerically study the effects of Born-Infeld nonlinear electrodynamics on the properties of $(1+1)$-dimensional s-wave holographic superconductors. We relax the probe limit and further assume the scalar and gauge fields affect on the background spacetime. We thus explore the effects of backreaction on the condensation of the scalar hair. For the analytical method, we employ the Sturm-Liouville eigen value problem and for the numerical method, we employ the shooting method. We show that these methods are powerful enough to analyze the critical temperature and phase transition of the one dimensional holographic superconductor. We find out that increasing the backreaction as well as nonlinearity makes the condensation harder to form. In addition, this one-dimensional holographic superconductor faces with second order phase transition and their critical exponent has the mean field value $\beta={1}/{2}$.


title: Comparisons and challenges of modern neutrino scattering experiments (TENSIONS2016 report)

abstract: Over the last decade, there has been enormous effort to measure neutrino interaction cross sections important to oscillation experiments. However, a number of results from modern experiments appear to be in tension with each other, despite purporting to measure the same processes. The TENSIONS2016 workshop was held at University of Pittsburgh July 24-31, 2016 and was sponsored by the Pittsburgh High Energy Physics, Astronomy, and Cosmology Center (PITT-PACC). The focus was on bringing experimentalists from three experiments together to compare results in detail and try to find the source of tension by clarifying and comparing signal definitions and the analysis strategies used for each measurement. A set of comparisons between the measurements using a consistent set of models was also made. This paper summarizes the main conclusions of that work.


title: Quantum electrodynamics of a superconductor-insulator phase transition

abstract: A chain of Josephson junctions implements one of the simplest many-body models undergoing a superconductor-insulator (SI) quantum phase transition between states with zero and infinite resistance. Apart from zero resistance, the superconducting state is necessarily accompanied by a sound-like mode due to collective oscillations of the phase of the complex-valued order parameter. Exciting this phase mode results in transverse photons propagating along the chain. Surprisingly little is known about the fate of this mode upon entering the insulating state, where the order parameter's amplitude remains non-zero, but the phase ordering is "melted" by quantum fluctuations. Here we report momentum-resolved radio-frequency spectroscopy of collective modes in nanofabricated chains of Al/AlOx/Al tunnel junctions. We find that the phase mode survives remarkably far into the insulating regime, such that $\textrm{M}\Omega$-resistance chains carry $\textrm{GHz}$-frequency alternating currents as nearly ideal superconductors. The insulator reveals itself through broadening and random frequency shifts of collective mode resonances, originated from intrinsic interactions. By pushing the chain parameters deeper into the insulating state, we achieved propagation with the speed of light down to $8\times 10^5~\textrm{m/s}$ and the wave impedance up to $23~\textrm{k}\Omega$. The latter quantity exceeds the predicted critical impedance by an order of magnitude, which opens the problem of quantum electrodynamics of a Bose glass insulator for both theory and experiment. Notably, the effective fine structure constant of such a 1D vacuum exceeds a unity, promising transformative applications to quantum science and technology.


title: Theoretical considerations about heavy-ion fusion in potential scattering

abstract: We carefully compare the one-dimensional WKB barrier tunneling model, and the one-channel Sch\"odinger equation with a complex optical potential calculation of heavy-ion fusion, for a light and a heavy system. It is found that the major difference between the two approaches occurs around the critical energy, above which the effective potential for the grazing angular momentum ceases to exhibit a pocket. The value of this critical energy is shown to be strongly dependent on the nuclear potential at short distances, on the inside region of the Coulomb barrier, and this dependence is much more important for heavy systems. Therefore the nuclear fusion process is expected to provide information on the nuclear potential in this inner region. We compare calculations with available data to show that the results are consistent with this expectation.


title: White Dwarfs as Dark Matter Detectors

abstract: Dark matter that is capable of sufficiently heating a local region in a white dwarf will trigger runaway fusion and ignite a type Ia supernova. This was originally proposed in Graham et al. (2015) and used to constrain primordial black holes which transit and heat a white dwarf via dynamical friction. In this paper, we consider dark matter (DM) candidates that heat through the production of high-energy standard model (SM) particles, and show that such particles will efficiently thermalize the white dwarf medium and ignite supernovae. Based on the existence of long-lived white dwarfs and the observed supernovae rate, we derive new constraints on ultra-heavy DM which produce SM particles through DM-DM annihilations, DM decays, and DM-SM scattering interactions in the stellar medium. As a concrete example, we rule out supersymmetric Q-ball DM in parameter space complementary to terrestrial bounds. We put further constraints on DM that is captured by white dwarfs, considering the formation and self-gravitational collapse of a DM core which heats the star via decays and annihilations within the core. It is also intriguing that the DM-induced ignition discussed in this work provide an alternative mechanism of triggering supernovae from sub-Chandrasekhar, non-binary progenitors.


title: Non-stationary localized oscillations of an infinite string, with time-varying tension, lying on the Winkler foundation with a point elastic inhomogeneity

abstract: We consider non-stationary oscillations of an infinite string with time-varying tension. The string lies on the Winkler foundation with a point inhomogeneity (a concentrated spring of negative stiffness). In such a system with constant parameters (the string tension), under certain conditions a trapped mode of oscillation exists and is unique. Therefore, applying a non-stationary external excitation to this system can lead to the emergence of the string oscillations localized near the inhomogeneity. We provide an analytical description of non-stationary localized oscillations of the string with slowly time-varying tension using the asymptotic procedure based on successive application of two asymptotic methods, namely the method of stationary phase and the method of multiple scales. The obtained analytical results were verified by independent numerical calculations based on the finite difference method. The applicability of the analytical formulas was demonstrated for various types of external excitation and laws governing the varying tension. In particular, we have shown that in the case when the trapped mode frequency approaches zero, localized low-frequency oscillations with increasing amplitude precede the localized string buckling. The dependence of the amplitude of such oscillations on its frequency is more complicated in comparison with the case of a one degree of freedom system with time-varying stiffness.


title: Impact of geometry on many-body localization

abstract: The impact of geometry on many body localization is studied on simple, exemplary systems amenable to exact diagonalization treatment. The crossover between ergodic and MBL phase for uniform as well as quasi-random disorder is analyzed using statistics of energy levels. It is observed that the transition to many-body localized phase is correlated with the number of nearest coupled neighbors. The crossover from extended to localized systems is approximately described by the so called plasma model.


title: Fixed-PSNR Lossy Compression for Scientific Data

abstract: Error-controlled lossy compression has been studied for years because of extremely large volumes of data being produced by today's scientific simulations. None of existing lossy compressors, however, allow users to fix the peak signal-to-noise ratio (PSNR) during compression, although PSNR has been considered as one of the most significant indicators to assess compression quality. In this paper, we propose a novel technique providing a fixed-PSNR lossy compression for scientific data sets. We implement our proposed method based on the SZ lossy compression framework and release the code as an open-source toolkit. We evaluate our fixed-PSNR compressor on three real-world high-performance computing data sets. Experiments show that our solution has a high accuracy in controlling PSNR, with an average deviation of 0.1 ~ 5.0 dB on the tested data sets.


title: Poincar\'e gauge gravity: An overview

abstract: We review the basics and the current status of the Poincar\'e gauge theory of gravity. The general dynamical scheme of Poincar\'e gauge gravity (PG) is formulated, and its physical consequences are outlined. In particular, we discuss exact solutions with and without torsion, highlight the cosmological aspects, and consider the probing of the spacetime geometry.


title: $T\bar{T}$ Partition Function from Topological Gravity

abstract: The $T{\bar T}$ deformation of a relativistic two-dimensional theory results in a solvable gravitational theory. Deformed scattering amplitudes can be obtained from coupling the undeformed theory to the flat space Jackiw--Teitelboim (JT) gravity. We show that the JT description is applicable and useful also in finite volume. Namely, we calculate the torus partition function of an arbitrary matter theory coupled to the JT gravity, formulated in the first order (vierbein) formalism. The first order description provides a natural set of dynamical clocks and rods for this theory, analogous to the target space coordinates in string theory. These dynamical coordinates play the role of relational observables allowing to define a torus path integral for the JT gravity. The resulting partition function is one-loop exact and reproduces the $T\bar{T}$ deformed finite volume spectrum.


title: Auto-locking Waveguide Amplifier System for Lidar and Magnetometric Applications

abstract: We describe a compact waveguide amplifier system that is suitable for optically pumping rubidium magnetometers. The system consists of an auto-locking vacuum-sealed external cavity diode laser, a semiconductor tapered amplifier and a pulsing unit based on an acousto-optic modulator. The diode laser utilises optical feedback from an interference filter to narrow the linewidth of an inexpensive laser diode to ~500 kHz. This output is scannable over an 8 GHz range (at 780 nm) and can be locked without human intervention to any spectral marker in an expandable library of reference spectra, using the autolocking controller. The tapered amplifier amplifies the output from 50 mW up to 2 W with negligible distortions in the spectral quality. The system can operate at visible and near infrared wavelengths with MHz repetition rates. We demonstrate optical pumping of rubidium vapour with this system for magnetometric applications. The magnetometer detects the differential absorption of two orthogonally polarized components of a linearly polarized probe laser following optical pumping by a circularly polarized pump laser. The differential absorption signal is studied for a range of pulse lengths, pulse amplitudes and DC magnetic fields. Our results suggest that this laser system is suitable for optically pumping spin-exchange free magnetometers.


title: Political Discussion and Leanings on Twitter: the 2016 Italian Constitutional Referendum

abstract: The recent availability of large, high-resolution data sets of online human activity allowed for the study and characterization of the mechanisms shaping human interactions at an unprecedented level of accuracy. To this end, many efforts have been put forward to understand how people share and retrieve information when forging their opinion about a certain topic. Specifically, the detection of the political leaning of a person based on its online activity can support the forecasting of opinion trends in a given population. Here, we tackle this challenging task by combining complex networks theory and machine learning techniques. In particular, starting from a collection of more than 6 millions tweets, we characterize the structure and dynamics of the Italian online political debate about the constitutional referendum held in December 2016. We analyze the discussion pattern between different political communities and characterize the network of contacts therein. Moreover, we set up a procedure to infer the political leaning of Italian Twitter users, which allows us to accurately reconstruct the overall opinion trend given by official polls (Pearson's r=0.88) as well as to predict with good accuracy the final outcome of the referendum. Our study provides a large-scale examination of the Italian online political discussion through sentiment-analysis, thus setting a baseline for future studies on online political debate modeling.


title: Batch Normalization in the final layer of generative networks

abstract: Generative Networks have shown great promise in generating photo-realistic images. Despite this, the theory surrounding them is still an active research area. Much of the useful work with Generative networks rely on heuristics that tend to produce good results. One of these heuristics is the advice not to use Batch Normalization in the final layer of the generator network. Many of the state-of-the-art generative network architectures use this heuristic, but the reasons for doing so are inconsistent. This paper will show that this is not necessarily a good heuristic and that Batch Normalization can be beneficial in the final layer of the generator network either by placing it before the final non-linear activation, usually a $tanh$ or replacing the final $tanh$ activation altogether with Batch Normalization and clipping. We show that this can lead to the faster training of Generator networks by matching the generator to the mean and standard deviation of the target distribution's image colour values.


title: Consistent Iron Abundances Derived from Neutral and Singly-Ionized Iron Lines in Ultraviolet and Optical Spectra of Six Warm Metal-Poor Stars

abstract: Neutral Fe lines in metal-poor stars yield conflicting abundances depending on whether and how deviations from local thermodynamic equilibrium (LTE) are considered. We have collected new high resolution and high signal-to-noise ultraviolet (UV) spectra of three warm dwarf stars with [Fe/H] = -2.9 with the Space Telescope Imaging Spectrograph on the Hubble Space Telescope. We locate archival UV spectra for three other warm dwarfs with [Fe/H] = -3.3, -2.2, and -1.6, supplemented with optical spectra for all six stars. We calculate stellar parameters using methods that are largely independent of the spectra, adopting broadband photometry, color-temperature relations, Gaia parallaxes, and assumed masses. We use the LTE line analysis code MOOG to derive Fe abundances from hundreds of Fe I and Fe II lines with wavelengths from 2290 to 6430 Angstroms. The [Fe/H] ratios derived separately from Fe I and Fe II lines agree in all six stars, with [Fe II/H] - [Fe I/H] ranging from +0.00 +/- 0.07 to -0.12 +/- 0.09 dex, when strong lines and Fe I lines with lower excitation potential < 1.2 eV are excluded. This constrains the extent of any deviations from LTE that may occur within this parameter range. While our result confirms non-LTE calculations for some warm, metal-poor dwarfs, it may not be generalizable to more metal-poor dwarfs, where deviations from LTE are predicted to be larger. We also investigate trends of systematically lower abundances derived from Fe I lines in the Balmer continuum region (3100-3700 Angstroms), and we conclude that no proposed explanation for this effect can fully account for the observations presently available.


title: Investigations of optical pumping for magnetometry using an auto-locking laser system

abstract: We have developed a versatile pulsed laser system for high precision magnetometry. The operating wavelength of the system can be configured to optically pump alkali vapors such as rubidium and cesium. The laser system consists of an auto-locked, interference filter stabilized, external cavity diode laser (ECDL), a tapered waveguide amplifier, and a pulsing module. The auto-locking controller can be used by an untrained operator to stabilize the laser frequency with respect to a library of atomic, molecular, and solid-state spectral markers. The ECDL output can be amplified from 20 mW to 2 W in continuous wave (CW) mode. The pulsing module, which includes an acousto-optic modulator (AOM), can generate pulses with durations of 20 ns and repetition rates of several MHz. Accordingly, the laser system is well suited for applications such as gravimetry, magnetometry, and differential-absorption lidar. In this work, we focus on magnetometric applications and demonstrate that the laser source is suitable for optically pumping rubidium vapor. We also describe numerical simulations of optical pumping relevant to the rubidium D1 and D2 transitions at 795 nm and 780 nm respectively. These studies are relevant to the design and construction of a new generation of portable, rubidium, spin-exchange relaxation-free (SERF) magnetometers, capable of sensitivities of 1 fT Hz-1/2 1.


title: Dynamic Monopolies in Reversible Bootstrap Percolation

abstract: We study an extremal question for the (reversible) $r-$bootstrap percolation processes. Given a graph and an initial configuration where each vertex is active or inactive, in the $r-$bootstrap percolation process the following rule is applied in discrete-time rounds: each vertex gets active if it has at least $r$ active neighbors, and an active vertex stays active forever. In the reversible $r$-bootstrap percolation, each vertex gets active if it has at least $r$ active neighbors, and inactive otherwise. We consider the following question on the $d$-dimensional torus: how many vertices should be initially active so that the whole graph becomes active? Our results settle an open problem by Balister, Bollob\'as, Johnson, and Walters and generalize the results by Flocchini, Lodi, Luccio, Pagli, and Santoro.


title: A Conformal Collider for Holographic CFTs

abstract: We develop a formalism to study the implications of causality on OPE coefficients in conformal field theories with large central charge and a sparse spectrum of higher spin operators. The formalism has the interpretation of a new conformal collider-type experiment for these class of CFTs and hence it has the advantage of requiring knowledge only about CFT three-point functions. This is accomplished by considering the holographic null energy operator which was introduced in arXiv:1709.03597 as a generalization of the averaged null energy operator. Analyticity properties of correlators in the Regge limit imply that the holographic null energy operator is a positive operator in a subspace of the total CFT Hilbert space. Utilizing this positivity condition, we derive bounds on three-point functions $\langle TO_1O_2\rangle$ of the stress tensor with various operators for CFTs with large central charge and a sparse spectrum. After imposing these constraints, we also find that the operator product expansions of all primary operators in the Regge limit have certain universal properties. All of these results are consistent with the expectation that CFTs in this class, irrespective of their microscopic details, admit universal gravity-like holographic dual descriptions. Furthermore, this connection enables us to constrain various inflationary observables such as the amplitude of chiral gravity waves, non-gaussanity of gravity waves and tensor-to-scalar ratio.


title: A Comparative Study of Single-Constraint Routing in Wireless Mesh Networks Using Different Dynamic Programming Algorithms

abstract: Finding the shortest route in wireless mesh networks is an important aspect. Many techniques are used to solve this problem like dynamic programming, evolutionary algorithms, weighted-sum techniques, and others. In this paper, we use dynamic programming techniques to find the shortest path in wireless mesh networks due to their generality, reduction of complexity and facilitation of numerical computation, simplicity in incorporating constraints, and their conformity to the stochastic nature of some problems. The routing problem is a multi-objective optimization problem with some constraints such as path capacity and end-to-end delay. Single-constraint routing problems and solutions using Dijkstra, Bellman- Ford, and Floyd-Warshall algorithms are proposed in this work with a discussion on the difference between them. These algorithms find the shortest route through finding the optimal rate between two nodes in the wireless networks but with bounded end-to-end delay. The Dijkstra-based algorithm is especially favorable in terms of processing time. We also present a comparison between our proposed single-constraint Dijkstra-based routing algorithm and the mesh routing algorithm (MRA) existing in the literature to clarify the merits of the former.


title: More green space is related to less antidepressant prescription rates in the Netherlands: A Bayesian geoadditive quantile regression approach

abstract: Exposure to green space seems to be beneficial for self-reported mental health. In this study we used an objective health indicator, namely antidepressant prescription rates. Current studies rely exclusively upon mean regression models assuming linear associations. It is, however, plausible that the presence of green space is non-linearly related with different quantiles of the outcome antidepressant prescription rates. These restrictions may contribute to inconsistent findings. Our aim was to assess antidepressant prescription rates in relation to green space, and to analyze how the relationship varies non-linearly across different quantiles of antidepressant prescription rates. We used cross-sectional data for the year 2014 at a municipality level in the Netherlands. Ecological Bayesian geoadditive quantile regressions were fitted for the 15, 50, and 85 percent quantiles to estimate green space-prescription rate correlations, controlling for confounders. The results suggested that green space was overall inversely and non-linearly associated with antidepressant prescription rates. More important, the associations differed across the quantiles, although the variation was modest. Significant non-linearities were apparent: The associations were slightly positive in the lower quantile and strongly negative in the upper one. Our findings imply that an increased availability of green space within a municipality may contribute to a reduction in the number of antidepressant prescriptions dispensed. Green space is thus a central health and community asset, whilst a minimum level of 28 percent needs to be established for health gains. The highest effectiveness occurred at a municipality surface percentage higher than 79 percent. This inverse dose-dependent relation has important implications for setting future community-level health and planning policies.


title: The Role of Models and Megamodels at Runtime

abstract: In model-driven software development a multitude of interrelated models are used to systematically realize a software system. This results in a complex development process since the models and the relations between the models have to be managed. Similar problems appear when following a model-driven approach for managing software systems at runtime. A multitude of interrelated runtime models are employed simultaneously, and thus they have to be maintained at runtime. While for the development case megamodels have emerged to address the problem of managing models and relations, the problem is rather neglected for the case of runtime models by applying ad-hoc solutions. Therefore, we propose to utilize megamodel concepts for the case of multiple runtime models. Based on the current state of research, we present a categorization of runtime models and conceivable relations between them. The categorization describes the role of interrelated models at runtime and demonstrates that several approaches already employ multiple runtime models and relations. Then, we show how megamodel concepts help in organizing and utilizing runtime models and relations in a model-driven manner while supporting a high level of automation. Finally, the role of interrelated models and megamodels at runtime is discussed for self-adaptive software systems and exemplified by a case study.


title: Adaptation and Abstract Runtime Models

abstract: Runtime adaptability is often a crucial requirement for today's complex software systems. Several approaches use an architectural model as a runtime representation of a managed system for monitoring, reasoning and performing adaptation. To ease the causal connection between a system and a model, these models are often closely related to the implementation and at a rather low level of abstraction. This makes them as complex as the implementation and it impedes reusability and extensibility of autonomic managers. Moreover, the models often do not cover different concerns, like security or performance, and therefore they do not support several self-management capabilities at once. In this paper we propose a model-driven approach that provides multiple architectural runtime models at different levels of abstraction as a basis for adaptation. Each runtime model abstracts from the underlying system and platform leveraging reusability and extensibility of managers that work on these models. Moreover, each model focuses on a specific concern which simplifies the work of autonomic managers. The different models are maintained automatically at runtime using model-driven engineering techniques that also reduce development efforts. Our approach has been implemented for the broadly adopted Enterprise Java Beans component standard and its application is presented in a self-healing scenario requiring structural adaptation.


title: Robust Handling of Polysemy via Sparse Representations

abstract: Words are polysemous and multi-faceted, with many shades of meanings. We suggest that sparse distributed representations are more suitable than other, commonly used, (dense) representations to express these multiple facets, and present Category Builder, a working system that, as we show, makes use of sparse representations to support multi-faceted lexical representations. We argue that the set expansion task is well suited to study these meaning distinctions since a word may belong to multiple sets with a different reason for membership in each. We therefore exhibit the performance of Category Builder on this task, while showing that our representation captures at the same time analogy problems such as "the Ganga of Egypt" or "the Voldemort of Tolkien". Category Builder is shown to be a more expressive lexical representation and to outperform dense representations such as Word2Vec in some analogy classes despite being shown only two of the three input terms.


title: Measurement of the production cross section for single top quarks in association with W bosons in proton-proton collisions at $\sqrt{s}=$ 13 TeV

abstract: A measurement is presented of the associated production of a single top quark and a W boson in proton-proton collisions at $\sqrt{s}=$ 13 TeV by the CMS Collaboration at the CERN LHC. The data collected corresponds to an integrated luminosity of 35.9 fb$^{-1}$. The measurement is performed using events with one electron and one muon in the final state along with at least one jet originated from a bottom quark. A multivariate discriminant, exploiting the kinematic properties of the events, is used to separate the signal from the dominant $\mathrm{t\overline{t}}$ background. The measured cross section of 63.1 $\pm$ 1.8 (stat) $\pm$ 6.4 (syst) $\pm$ 2.1 (lumi) pb is in agreement with the standard model expectation.


title: Resolvent Estimates on Asymptotically Hyperbolic Spaces

abstract: We extend Vasy's results on semiclassical high energy estimates for the meromorphic continuation of the resolvent for asymptotically hyperbolic manifolds to metrics that are not necessarily even. Vasy's method gives the meromorphic continuation of the resolvent and high energy estimates in strips, assuming that the geodesic flow is non-trapping, without having to construct a parametrix. We prove that the size of the strip to which the resolvent can be extended meromorphically is the same obtained by Guillarmou.


title: Mass-improvement of the vector current in three-flavor QCD

abstract: We determine two improvement coefficients which are relevant to cancel mass-dependent cutoff effects in correlation functions with operator insertions of the non-singlet local QCD vector current. This determination is based on degenerate three-flavor QCD simulations of non-perturbatively O(a) improved Wilson fermions with tree-level improved gauge action. Employing a very robust strategy that has been pioneered in the quenched approximation leads to an accurate estimate of a counterterm cancelling dynamical quark cutoff effects linear in the trace of the quark mass matrix. To our knowledge this is the first time that such an effect has been determined systematically with large significance.


title: Ultrafast scanning electron microscope applied for studying the interaction between free electrons and optical near-fields of periodic nanostructures

abstract: In this paper we describe an ultrafast scanning electron microscope setup developed for the research of inelastic scattering of electrons at optical near-fields of periodic dielectric nanostructures. Electron emission from the Schottky cathode is controlled by ultraviolet femtosecond laser pulses. The electron pulse duration at the interaction site is characterized via cross-correlation of the electrons with an infrared laser pulse that excites a synchronous periodic near-field on the surface of a silicon nanostructure. The lower limit of 410 fs is found in the regime of a single electron per pulse. The role of pulse broadening due to Coulomb interaction in multielectron pulses is investigated. The setup is used to demonstrate an increase of the interaction distance between the electrons and the optical near-fields by introducing a pulse-front-tilt to the infrared laser beam. Further we show the dependence of the final electron spectra on the resonance condition between the phase velocity of the optical near-field and the electron propagation velocity. The resonance is controlled by adjusting the initial electron energy/velocity and by introducing a linear chirp to the structure period allowing to increase the final electron energy gain up to a demonstrated 3.8 keV.


title: Asset Price Bubbles: An Option-based Indicator

abstract: We construct a statistical indicator for the detection of short-term asset price bubbles based on the information content of bid and ask market quotes for plain vanilla put and call options. Our construction makes use of the martingale theory of asset price bubbles and the fact that such scenarios where the price for an asset exceeds its fundamental value can in principle be detected by analysis of the asymptotic behavior of the implied volatility surface. For extrapolating this implied volatility, we choose the SABR model, mainly because of its decent fit to real option market quotes for a broad range of maturities and its ease of calibration. As main theoretical result, we show that under lognormal SABR dynamics, we can compute a simple yet powerful closed-form martingale defect indicator by solving an ill-posed inverse calibration problem. In order to cope with the ill-posedness and to quantify the uncertainty which is inherent to such an indicator, we adopt a Bayesian statistical parameter estimation perspective. We probe the resulting posterior densities with a combination of optimization and adaptive Markov chain Monte Carlo methods, thus providing a full-blown uncertainty estimation of all the underlying parameters and the martingale defect indicator. Finally, we provide real-market tests of the proposed option-based indicator with focus on tech stocks due to increasing concerns about a tech bubble 2.0.


title: Schmidt gap in random spin chains

abstract: We numerically investigate the low-lying entanglement spectrum of the ground state of random one-dimensional spin chains obtained after partition of the chain into two equal halves. We consider two paradigmatic models: the spin-1/2 random transverse field Ising model, solved exactly, and the spin-1 random Heisenberg model, simulated using the density matrix renormalization group. In both cases we analyze the mean Schmidt gap, defined as the difference between the two largest eigenvalues of the reduced density matrix of one of the two partitions, averaged over many disorder realizations. We find that the Schmidt gap detects the critical point very well and scales with universal critical exponents.


title: Processing of missing data by neural networks

abstract: We propose a general, theoretically justified mechanism for processing missing data by neural networks. Our idea is to replace typical neuron's response in the first hidden layer by its expected value. This approach can be applied for various types of networks at minimal cost in their modification. Moreover, in contrast to recent approaches, it does not require complete data for training. Experimental results performed on different types of architectures show that our method gives better results than typical imputation strategies and other methods dedicated for incomplete data.


title: An efficient hybrid orbital representation for quantum Monte Carlo calculations

abstract: The scale and complexity of quantum system to which real-space quantum Monte Carlo (QMC) can be applied in part depends on the representation and memory usage of the trial wavefunction. B-splines, the computationally most efficient basis set, can have memory requirements exceeding the capacity of a single computational node. This situation has traditionally forced a difficult choice of either using slow internode communication or a potentially less accurate but smaller basis set such as Gaussians. Here, we introduce a hybrid representation of the single particle orbitals that combine a localized atomic basis set around atomic cores and B-splines in the interstitial regions to reduce the memory usage while retaining high speed of evaluation and either retaining or increasing overall accuracy. We present a benchmark calculation for NiO demonstrating a superior accuracy while using only one eighth the memory required for conventional B-splines. The hybrid orbital representation therefore expands the overall range of systems that can be practically studied with QMC.


title: Quantifying the abundance of faint, low-redshift satellite galaxies in the COSMOS survey

abstract: Faint dwarf satellite galaxies are important as tracers of small-scale structure, but remain poorly characterized outside the Local Group, due to the difficulty of identifying them consistently at larger distances. We review a recently proposed method for estimating the average satellite population around a given sample of nearby bright galaxies, using a combination of size and magnitude cuts (to select low-redshift dwarf galaxies preferentially) and clustering measurements (to estimate the fraction of true satellites in the cut sample). We test this method using the high-precision photometric redshift catalog of the COSMOS survey, exploring the effect of specific cuts on the clustering signal. The most effective of the size-magnitude cuts considered recover the clustering signal around low-redshift primaries (z < 0.15) with about two-thirds of the signal and 80\%\ of the signal-to-noise ratio obtainable using the full COSMOS photometric redshifts. These cuts are also fairly efficient, with more than one third of the selected objects being clustered satellites. We conclude that structural selection represents a useful tool in characterizing dwarf populations to fainter magnitudes and/or over larger areas than are feasible with spectroscopic surveys.} In reviewing the low-redshift content of the COSMOS field, we also note the existence of several dozen objects that appear resolved or partially resolved in the HST imaging, and are confirmed to be local (at distances of $\sim$250 Mpc or less) by their photometric or spectroscopic redshifts. This underlines the potential for future space-based surveys to reveal local populations of intrinsically faint galaxies through imaging alone.


title: Hyperuniform states generated by a critical friction field

abstract: Hyperuniform states are an efficient way to fill up space for disordered systems. In these states the particle distribution is disordered at the short scale but becomes increasingly uniform when looked at large scales. Hyperuniformity appears in several systems, in static or quasistatic regimes as well as close to transitions to absorbing states. Here, we show that a vibrated granular layer, at the critical point of the liquid-to-solid transition, displays dynamic hyperuniformity. Prior to the transition, patches of the solid phase form, with length scales and mean lifetimes that diverge critically at the transition point. When reducing the wavenumber, density fluctuations encounter increasingly more patches that block their propagation, resulting in a static structure factor that tends to zero for small wavenumbers at the critical point, which is a signature of hyperuniformity. A simple model demonstrates that this coupling of a density field to a highly fluctuating scalar friction field gives rise to dynamic hyperuniform states. Finally, we show that the structure factor detects better the emergence of hyperuniformity, compared to the particle number variance.


title: LECTOR Based Clock Gating for Low Power Multi-Stage Flip Flop Applications

abstract: Power dissipation in integrated circuits is one of the major concerns to the research community, at the verge when more number of transistors are integrated on a single chip. The substantial source of power dissipation in sequential elements of the integrated circuit is due to the fast switching of high frequency clock signals. These signals do not carry any information and are mainly intended to synchronize the operation of sequential components. This unnecessary switching of Clock, during the HOLD phase of either logic 1 or logic 0, may be eliminated using a technique, called Clock Gating. In this paper, we have incorporated a recent clock gating style called LECTOR based clock gating LB CG to drive multi stage architecture and simulated its performance using 90nm CMOS Predictive Technology Model PTM with a power supply of 1.1V at 18GHz clock frequency. A substantial savings in terms of average power in comparison to its non gated correspondent have been observed.


title: Learning to Collaborate for User-Controlled Privacy

abstract: It is becoming increasingly clear that users should own and control their data. Utility providers are also becoming more interested in guaranteeing data privacy. As such, users and utility providers should collaborate in data privacy, a paradigm that has not yet been developed in the privacy research community. We introduce this concept and present explicit architectures where the user controls what characteristics of the data she/he wants to share and what she/he wants to keep private. This is achieved by collaborative learning a sensitization function, either a deterministic or a stochastic one, that retains valuable information for the utility tasks but it also eliminates necessary information for the privacy ones. As illustration examples, we implement them using a plug-and-play approach, where no algorithm is changed at the system provider end, and an adversarial approach, where minor re-training of the privacy inferring engine is allowed. In both cases the learned sanitization function keeps the data in the original domain, thereby allowing the system to use the same algorithms it was using before for both original and privatized data. We show how we can maintain utility while fully protecting private information if the user chooses to do so, even when the first is harder than the second, as in the case here illustrated of identity detection while hiding gender.


title: Discovery of Nonlinear Multiscale Systems: Sampling Strategies and Embeddings

abstract: A major challenge in the study of dynamical systems is that of model discovery: turning data into models that are not just predictive, but provide insight into the nature of the underlying dynamical system that generated the data. This problem is made more difficult by the fact that many systems of interest exhibit diverse behaviors across multiple time scales. We introduce a number of data-driven strategies for discovering nonlinear multiscale dynamical systems and their embeddings from data. We consider two canonical cases: (i) systems for which we have full measurements of the governing variables, and (ii) systems for which we have incomplete measurements. For systems with full state measurements, we show that the recent sparse identification of nonlinear dynamical systems (SINDy) method can discover governing equations with relatively little data and introduce a sampling method that allows SINDy to scale efficiently to problems with multiple time scales. Specifically, we can discover distinct governing equations at slow and fast scales. For systems with incomplete observations, we show that the Hankel alternative view of Koopman (HAVOK) method, based on time-delay embedding coordinates, can be used to obtain a linear model and Koopman invariant measurement system that nearly perfectly captures the dynamics of nonlinear quasiperiodic systems. We introduce two strategies for using HAVOK on systems with multiple time scales. Together, our approaches provide a suite of mathematical strategies for reducing the data required to discover and model nonlinear multiscale systems.


title: Wasserstein Measure Coresets

abstract: The proliferation of large data sets and Bayesian inference techniques motivates demand for better data sparsification. Coresets provide a principled way of summarizing a large dataset via a smaller one that is guaranteed to match the performance of the full data set on specific problems. Classical coresets, however, neglect the underlying data distribution, which is often continuous. We address this oversight by introducing Wasserstein measure coresets, an extension of coresets which by definition takes into account generalization. Our formulation of the problem, which essentially consists in minimizing the Wasserstein distance, is solvable via stochastic gradient descent. This yields an algorithm which simply requires sample access to the data distribution and is able to handle large data streams in an online manner. We validate our construction for inference and clustering.


title: Assessing Health Care Interventions via an Interrupted Time Series Model: Study Power and Design Considerations

abstract: The delivery and assessment of quality health care is complex with many interacting and interdependent components. In terms of research design and statistical analysis, this complexity and interdependency makes it difficult to assess the true impact of interventions designed to improve patient health care outcomes. Interrupted time series (ITS) is a quasi-experimental design developed for inferring the effectiveness of a health policy intervention while accounting for temporal dependence within a single system or unit. Current standardized ITS methods do not simultaneously analyze data for several units, nor are there methods to test for the existence of a change point and to assess statistical power for study planning purposes in this context. To address this limitation we propose the `Robust Multiple ITS' (R-MITS) model, appropriate for multi-unit ITS data, that allows for inference regarding the estimation of a global change point across units in the presence of a potentially lagged (or anticipatory) treatment effect. Under the R-MITS model, one can formally test for the existence of a change point and estimate the time delay between the formal intervention implementation and the over-all-unit intervention effect. We conducted empirical simulation studies to assess the type one error rate of the testing procedure, power for detecting specified change-point alternatives, and accuracy of the proposed estimating methodology. R-MITS is illustrated by analyzing patient satisfaction data from a hospital that implemented and evaluated a new care delivery model in multiple units.


title: Quadrature Histograms in Maximum Likelihood Quantum State Tomography

abstract: Quantum state tomography aims to determine the quantum state of a system from measured data and is an essential tool for quantum information science. When dealing with continuous variable quantum states of light, tomography is often done by measuring the field amplitudes at different optical phases using homodyne detection. The quadrature-phase homodyne measurement outputs a continuous variable, so to reduce the computational cost of tomography, researchers often discretize the measurements. We show that this can be done without significantly degrading the fidelity between the estimated state and the true state. This paper studies different strategies for determining the histogram bin widths. We show that computation time can be significantly reduced with little loss in the fidelity of the estimated state when the measurement operators corresponding to each histogram bin are integrated over the bin width.


title: Effects of dissociation/recombination on the day-night temperature contrasts of ultra-hot Jupiters

abstract: Secondary eclipse observations of ultra-hot Jupiters have found evidence that hydrogen is dissociated on their daysides. Additionally, full-phase light curve observations of ultra-hot Jupiters show a smaller day-night emitted flux contrast than that expected from previous theory. Recently, it was proposed by Bell & Cowan (2018) that the heat intake to dissociate hydrogen and heat release due to recombination of dissociated hydrogen can affect the atmospheric circulation of ultra-hot Jupiters. In this work, we add cooling/heating due to dissociation/recombination into the analytic theory of Komacek & Showman (2016) and Zhang & Showman (2017) for the dayside-nightside temperature contrasts of hot Jupiters. We find that at high values of incident stellar flux, the day-night temperature contrast of ultra-hot Jupiters may decrease with increasing incident stellar flux due to dissociation/recombination, the opposite of that expected without including the effects of dissociation/recombination. We propose that a combination of a greater number of full-phase light curve observations of ultra-hot Jupiters and future General Circulation Models that include the effects of dissociation/recombination could determine in detail how the atmospheric circulation of ultra-hot Jupiters differs from that of cooler planets.


title: Computing Kantorovich-Wasserstein Distances on $d$-dimensional histograms using $(d+1)$-partite graphs

abstract: This paper presents a novel method to compute the exact Kantorovich-Wasserstein distance between a pair of $d$-dimensional histograms having $n$ bins each. We prove that this problem is equivalent to an uncapacitated minimum cost flow problem on a $(d+1)$-partite graph with $(d+1)n$ nodes and $dn^{\frac{d+1}{d}}$ arcs, whenever the cost is separable along the principal $d$-dimensional directions. We show numerically the benefits of our approach by computing the Kantorovich-Wasserstein distance of order 2 among two sets of instances: gray scale images and $d$-dimensional biomedical histograms. On these types of instances, our approach is competitive with state-of-the-art optimal transport algorithms.


title: On the Yang-Baxter Poisson algebra in non-ultralocal integrable systems

abstract: A common approach to the quantization of integrable models starts with the formal substitution of the Yang-Baxter Poisson algebra with its quantum version. However it is difficult to discern the presence of such an algebra for the so-called non-ultralocal models. The latter includes the class of non-linear sigma models which are most interesting from the point of view of applications. In this work, we investigate the emergence of the Yang-Baxter Poisson algebra in a non-ultralocal system which is related to integrable deformations of the Principal Chiral Field.


title: Sequential Learning of Principal Curves: Summarizing Data Streams on the Fly

abstract: When confronted with massive data streams, summarizing data with dimension reduction methods such as PCA raises theoretical and algorithmic pitfalls. Principal curves act as a nonlinear generalization of PCA and the present paper proposes a novel algorithm to automatically and sequentially learn principal curves from data streams. We show that our procedure is supported by regret bounds with optimal sublinear remainder terms. A greedy local search implementation (called \texttt{slpc}, for Sequential Learning Principal Curves) that incorporates both sleeping experts and multi-armed bandit ingredients is presented, along with its regret computation and performance on synthetic and real-life data.


title: Palatini formulation of $f(R,T)$ gravity theory, and its cosmological implications

abstract: We consider the Palatini formulation of $f(R,T)$ gravity theory, in which a nonminimal coupling between the Ricci scalar and the trace of the energy-momentum tensor is introduced, by considering the metric and the affine connection as independent field variables. The field equations and the equations of motion for massive test particles are derived, and we show that the independent connection can be expressed as the Levi-Civita connection of an auxiliary, energy-momentum trace dependent metric, related to the physical metric by a conformal transformation. Similarly to the metric case, the field equations impose the non-conservation of the energy-momentum tensor. We obtain the explicit form of the equations of motion for massive test particles in the case of a perfect fluid, and the expression of the extra-force, which is identical to the one obtained in the metric case. The thermodynamic interpretation of the theory is also briefly discussed. We investigate in detail the cosmological implications of the theory, and we obtain the generalized Friedmann equations of the $f(R,T)$ gravity in the Palatini formulation. Cosmological models with Lagrangians of the type $f=R-\alpha ^2/R+g(T)$ and $f=R+\alpha ^2R^2+g(T)$ are investigated. These models lead to evolution equations whose solutions describe accelerating Universes at late times.


title: Photon entanglement entropy as a probe of many-body correlations and fluctuations

abstract: Recent theoretical and experiments have explored the use of entangled photons as a spectroscopic probe of material systems. We develop here a theoretical description for entropy production in the scattering of an entangled biphoton state within an optical cavity. We develop this using perturbation theory by expanding the biphoton scattering matrix in terms of single-photon terms in which we introduce the photon-photon interaction via a complex coupling constant, $\xi$. We show that the von Neumann entropy provides a succinct measure of this interaction. We then develop a microscopic model and show that in the limit of fast fluctuations, the entanglement entropy vanishes whereas in the limit the coupling is homogeneous broadened, the entanglement entropy depends upon the magnitude of the fluctuations and reaches a maximum.


title: Visualization of Surface-acoustic-wave Potential by Transmission-mode Microwave Impedance Microscopy

abstract: Elastic waves propagating in piezoelectric materials are accompanied by a time-varying electric potential, which is of critical importance for acousto-electronic applications. The spatial mapping of such a potential at microwave frequencies is challenging since the characteristic length scale is determined by the acoustic wavelength of several micrometers. In this work, we report the visualization of surface acoustic waves (SAWs) on ferroelectric samples by transmission-mode microwave impedance microscopy (T-MIM). The SAW potential launched by the interdigital transducer is detected by the tip and demodulated by the microwave electronics as time-independent spatial patterns. Wave phenomena such as interference and diffraction are imaged and the results are in excellent agreement with the theoretical analysis. Our work opens up a new avenue to study various electromechanical systems in a spatially resolved manner.


title: Review of Cyber-Physical Attacks and Counter Defense Mechanisms for Advanced Metering Infrastructure in Smart Grid

abstract: The Advanced Metering Infrastructure (AMI) is a vital element in the current development of the smart grid. AMI technologies provide electric utilities with an effective way of continuous monitoring and remote control of smart grid components. However, owing to its increasing scale and cyber-physical nature, the AMI has been faced with security threats in both cyber and physical domains. This paper provides a comprehensive review of the crucial cyber-physical attacks and counter defense mechanisms in the AMI. First, two attack surfaces are surveyed in the AMI including the communication network and smart meters. The potential cyber-physical attacks are then reviewed for each attack surface. Next, the attack models and their cyber and physical impacts on the smart grid are studied for comparison. Counter defense mechanisms that help mitigate these security threats are discussed. Finally, several mathematical tools which may help in analysis and implementation of security solutions are summarized.


title: Efficient simulation of Gaussian Markov random fields by Chebyshev polynomial approximation

abstract: This paper presents an algorithm to simulate Gaussian random vectors whose precision matrix can be expressed as a polynomial of a sparse matrix. This situation arises in particular when simulating Gaussian Markov random fields obtained by the discretization by finite elements of the solutions of some stochastic partial derivative equations. The proposed algorithm uses a Chebyshev polynomial approximation to compute simulated vectors with a linear complexity. This method is asymptotically exact as the approximation order grows. Criteria based on tests of the statistical properties of the produced vectors are derived to determine minimal orders of approximation.


title: Surgery induces exact sequences in Lagrangian cobordisms

abstract: We prove that if L_0 and L_1 are exact branes intersecting in precisely one point, then there exists a fiber sequence in the infinity-category of Lagrangian cobordisms consisting of L_0, L_1, and a surgery of L_0 with L_1. By combining this with the exact functor from [Tan], we find analogues of results of Biran and Cornea in the wrapped and exact setting.


title: Combinatorial Properties of Metrically Homogeneous Graphs

abstract: Ramsey theory looks for regularities in large objects. Model theory studies algebraic structures as models of theories. The structural Ramsey theory combines these two fields and is concerned with Ramsey-type questions about certain model-theoretic structures. In 2005, Ne\v{s}et\v{r}il initiated a systematic study of the so-called Ramsey classes of finite structures. This thesis is a contribution to the programme; we find Ramsey expansions of the primitive 3-constrained classes from Cherlin's catalogue of metrically homogeneous graphs. A key ingradient is an explicit combinatorial algorithm to fill-in the missing distances in edge-labelled graphs to obtain structures from Cherlin's classes. This algorithm also implies the extension property for partial automorphisms (EPPA), another combinatorial property of classes of finite structures.


title: Incept-N: A Convolutional Neural Network based Classification Approach for Predicting Nationality from Facial Features

abstract: The nationality of a human being is a well-known identifying characteristic used for every major authentication purpose in every country. Albeit advances in the application of Artificial Intelligence and Computer Vision in different aspects, its contribution to this specific security procedure is yet to be cultivated. With a goal to successfully applying computer vision techniques to predict the nationality of a person based on his facial features, we have proposed this novel method and have achieved an average of 93.6% accuracy with very low misclassification rate.


title: Method G: Uncertainty Quantification for Distributed Data Problems using Generalized Fiducial Inference

abstract: It is not unusual for a data analyst to encounter data sets distributed across several computers. This can happen for reasons such as privacy concerns, efficiency of likelihood evaluations, or just the sheer size of the whole data set. This presents new challenges to statisticians as even computing simple summary statistics such as the median becomes computationally challenging. Furthermore, if other advanced statistical methods are desired, novel computational strategies are needed. In this paper we propose a new approach for distributed analysis of massive data that is suitable for generalized fiducial inference and is based on a careful implementation of a "divide and conquer" strategy combined with importance sampling. The proposed approach requires only small amount of communication between nodes, and is shown to be asymptotically equivalent to using the whole data set. Unlike most existing methods, the proposed approach produces uncertainty measures (such as confidence intervals) in addition to point estimates for parameters of interest. The proposed approach is also applied to the analysis of a large set of solar images.


title: Schr\"odinger formalism for a particle constrained to a surface in $\mathbb{R}_1^3$

abstract: In this work it is studied the Schr\"odinger equation for a non-relativistic particle restricted to move on a surface $S$ in a three-dimensional Minkowskian medium $\mathbb{R}_1^3$, i.e., the space $\mathbb{R}^3$ equipped with the metric $\text{diag}(-1,1,1)$. After establishing the consistency of the interpretative postulates for the new Schr\"odinger equation, namely the conservation of probability and the hermiticity of the new Hamiltonian built out of the Laplacian in $\mathbb{R}_1^3$, we investigate the confining potential formalism in the new effective geometry. Like in the well-known Euclidean case, it is found a geometry-induced potential acting on the dynamics $V_S = - \frac{\hbar^{2}}{2m} \left(\varepsilon H^2-K\right)$ which, besides the usual dependence on the mean ($H$) and Gaussian ($K$) curvatures of the surface, has the remarkable feature of a dependence on the signature of the induced metric of the surface: $\varepsilon= +1$ if the signature is $(-,+)$, and $\varepsilon=1$ if the signature is $(+,+)$. Applications to surfaces of revolution in $\mathbb{R}^3_1$ are examined, and we provide examples where the Schr\"odinger equation is exactly solvable. It is hoped that our formalism will prove useful in the modeling of novel materials such as hyperbolic metamaterials, which are characterized by a hyperbolic dispersion relation, in contrast to the usual spherical (elliptic) dispersion typically found in conventional materials.


title: Designing communication systems via iterative improvement: error correction coding with Bayes decoder and codebook optimized for source symbol error

abstract: In most error correction coding (ECC) frameworks, the typical error metric is the bit error rate (BER) which measures the number of bit errors. For this metric, the positions of the bits are not relevant to the decoding, and in many noise models, not relevant to the BER either. In many applications this is unsatisfactory as typically all bits are not equal and have different significance. We consider the problem of bit error correction and mitigation where bits in different positions have different importance. For error correction, we look at ECC from a Bayesian perspective and introduce Bayes estimators with general loss functions to take into account the bit significance. We propose ECC schemes that optimize this error metric. As the problem is highly nonlinear, traditional ECC construction techniques are not applicable. Using exhaustive search is cost prohibitive, and thus we use iterative improvement search techniques to find good codebooks. We optimize both general codebooks and linear codes. We provide numerical experiments to show that they can be superior to classical linear block codes such as Hamming codes and decoding methods such as minimum distance decoding. For error mitigation, we study the case where ECC is not possible or not desirable, but significance aware encoding of information is still beneficial in reducing the average error. We propose a novel number presentation format suitable for emerging storage media where the noise magnitude is unknown and possibly large and show that it has lower mean error than the traditional number format.


title: Efficient Online Portfolio with Logarithmic Regret

abstract: We study the decades-old problem of online portfolio management and propose the first algorithm with logarithmic regret that is not based on Cover's Universal Portfolio algorithm and admits much faster implementation. Specifically Universal Portfolio enjoys optimal regret $\mathcal{O}(N\ln T)$ for $N$ financial instruments over $T$ rounds, but requires log-concave sampling and has a large polynomial running time. Our algorithm, on the other hand, ensures a slightly larger but still logarithmic regret of $\mathcal{O}(N^2(\ln T)^4)$, and is based on the well-studied Online Mirror Descent framework with a novel regularizer that can be implemented via standard optimization methods in time $\mathcal{O}(TN^{2.5})$ per round. The regret of all other existing works is either polynomial in $T$ or has a potentially unbounded factor such as the inverse of the smallest price relative.


title: Can machine learning identify interesting mathematics? An exploration using empirically observed laws

abstract: We explore the possibility of using machine learning to identify interesting mathematical structures by using certain quantities that serve as fingerprints. In particular, we extract features from integer sequences using two empirical laws: Benford's law and Taylor's law and experiment with various classifiers to identify whether a sequence is, for example, nice, important, multiplicative, easy to compute or related to primes or palindromes.


title: Reduction of power grid fluctuations by communication between smart devices

abstract: The increase of electric demand and the progressive integration of renewable sources threatens the stability of the power grid. To solve this issue, several methods have been proposed to control the demand side instead of increasing the spinning reserve on the supply side. Here we focus on dynamic demand control (DDC), a method in which appliances can delay its scheduled operation if the electric frequency is outside a suitable range. We have recently shown that DDC effectively reduces small and medium-size frequency fluctuations but, due to the need of recovering pending tasks, the probability of large demand peaks, and hence large frequency fluctuations, may actually increase. Although these events are very rare they can potentially trigger a failure of the system and therefore strategies to avoid them have to be addressed. In this work, we introduce a new method including communication among DDC devices belonging to a given group, such that they can coordinate opposite actions to keep the group demand more stable. We show that for this method the amount of pending tasks decreases by a factor 10 while large frequency fluctuations are significantly reduced or even completely avoided.


title: DeepLogic: Towards End-to-End Differentiable Logical Reasoning

abstract: Combining machine learning with logic-based expert systems in order to get the best of both worlds are becoming increasingly popular. However, to what extent machine learning can already learn to reason over rule-based knowledge is still an open problem. In this paper, we explore how symbolic logic, defined as logic programs at a character level, is learned to be represented in a high-dimensional vector space using RNN-based iterative neural networks to perform reasoning. We create a new dataset that defines 12 classes of logic programs exemplifying increased level of complexity of logical reasoning and train the networks in an end-to-end fashion to learn whether a logic program entails a given query. We analyse how learning the inference algorithm gives rise to representations of atoms, literals and rules within logic programs and evaluate against increasing lengths of predicate and constant symbols as well as increasing steps of multi-hop reasoning.


title: Reachability Analysis for Spatial Concurrent Constraint Systems with Extrusion

abstract: Spatial concurrent constraint programming (SCCP) is an algebraic model of spatial modalities in constrained-based process calculi; it can be used to reason about spatial information distributed among the agents of a system. This work presents an executable rewriting logic semantics of SCCP with extrusion (i.e., process mobility) that uses rewriting modulo SMT, a novel technique that combines the power of term rewriting, matching algorithms, and SMT-solving. In this setting, constraints are encoded as formulas in a theory with a satisfaction relation decided by an SMT solver, while the topology of the spatial hierarchy is encoded as part of the term structure of symbolic states. By being executable, the rewriting logic specification offers support for the inherent symbolic and challenging task of reachability analysis in the constrained-based model. The approach is illustrated with examples about the automatic verification of fault-tolerance, consistency, and privacy in distributed spatial and hierarchical systems.


title: A case study of hurdle and generalized additive models in astronomy: the escape of ionizing radiation

abstract: The dark ages of the Universe end with the formation of the first generation of stars residing in primeval galaxies. These objects were the first to produce ultraviolet ionizing photons in a period when the cosmic gas changed from a neutral state to an ionized one, known as Epoch of Reionization (EoR). A pivotal aspect to comprehend the EoR is to probe the intertwined relationship between the fraction of ionizing photons capable to escape dark haloes, also known as the escape fraction ($f_{esc}$), and the physical properties of the galaxy. This work develops a sound statistical model suitable to account for such non-linear relationships and the non-Gaussian nature of $f_{esc}$. This model simultaneously estimates the probability that a given primordial galaxy starts the ionizing photon production and estimates the mean level of the $f_{esc}$ once it is triggered. The model was employed in the First Billion Years simulation suite, from which we show that the baryonic fraction and the rate of ionizing photons appear to have a larger impact on $f_{esc}$ than previously thought. A naive univariate analysis of the same problem would suggest smaller effects for these properties and a much larger impact for the specific star formation rate, which is lessened after accounting for other galaxy properties and non-linearities in the statistical model.


title: Structural discontinuity as generalized strain and Fourier transform for discrete-continuous systems

abstract: We consider a segmented structure, possibly connected with a continuous medium, as initially homogeneous, where discontinuities arise as localized strains induced by self-equilibrated localized actions. Under this formulation augmented by interface conditions, the linearized formulation remains valid. This approach eliminates the need for examining separate sections with subsequent conjugation. Only conditions related to the discontinuities should be satisfied, while the continuity in other respects preserves itself automatically. No obstacle remains for the continuous Fourier transform. For a uniform partitioning, the discrete transform is used together with the continuous one. We demonstrate the technique by obtaining the Floquet wave dispersive relations with their dependence upon interface stiffness. To this end, we briefly consider the flexural wave in the segmented beam on Winkler's foundation, the gravity wave in a plate (also segmented) on deep water and the Floquet-Rayleigh wave in such a plate on an elastic half-space. Besides, we present the wave equations developed for an elastic medium with discontinuities.


title: Resolving Force Indeterminacy in Contact Dynamics Using Compatibility Conditions

abstract: Contact dynamics (CD) is a powerful method to solve the dynamics of large systems of colliding rigid bodies. CD can be computationally more efficient than classical penalty-based discrete element methods (DEM) for simulating contact between stiff materials such as rock, glass, or engineering metals. However, by idealizing bodies as perfectly rigid, contact forces computed by CD can be non-unique due to indeterminacy in the contact network, which is a common occurence in dense granular flows. We propose a CD method that is designed to identify only the unique set of contact forces that would be predicted by a soft particle method, such as DEM, in the limit of large stiffness. The method involves applying an elastic compatibility condition to the contact forces, which maintains no-penetration constraints but filters out force distributions that could not have arisen from stiff elastic contacts. The method can be used as a post-processing step that could be integrated into existing CD codes with minimal effort. We demonstrate its efficacy in a variety of indeterminate problems, including some involving multiple materials, non-spherical shapes, and nonlinear contact constitutive laws.


title: Region-Based Classification of PolSAR Data Using Radial Basis Kernel Functions With Stochastic Distances

abstract: Region-based classification of PolSAR data can be effectively performed by seeking for the assignment that minimizes a distance between prototypes and segments. Silva et al (2013) used stochastic distances between complex multivariate Wishart models which, differently from other measures, are computationally tractable. In this work we assess the robustness of such approach with respect to errors in the training stage, and propose an extension that alleviates such problems. We introduce robustness in the process by incorporating a combination of radial basis kernel functions and stochastic distances with Support Vector Machines (SVM). We consider several stochastic distances between Wishart: Bhatacharyya, Kullback-Leibler, Chi-Square, R\'{e}nyi, and Hellinger. We perform two case studies with PolSAR images, both simulated and from actual sensors, and different classification scenarios to compare the performance of Minimum Distance and SVM classification frameworks. With this, we model the situation of imperfect training samples. We show that SVM with the proposed kernel functions achieves better performance with respect to Minimum Distance, at the expense of more computational resources and the need of parameter tuning. Code and data are provided for reproducibility.


title: Acceleration of Non-Linear Minimisation with PyTorch

abstract: I show that a software framework intended primarily for training of neural networks, PyTorch, is easily applied to a general function minimisation problem in science. The qualities of PyTorch of ease-of-use and very high efficiency are found to be applicable in this domain and lead to two orders of magnitude improvement in time-to-solution with very small software engineering effort.


title: Neural Architecture Search using Deep Neural Networks and Monte Carlo Tree Search

abstract: Neural Architecture Search (NAS) has shown great success in automating the design of neural networks, but the prohibitive amount of computations behind current NAS methods requires further investigations in improving the sample efficiency and the network evaluation cost to get better results in a shorter time. In this paper, we present a novel scalable Monte Carlo Tree Search (MCTS) based NAS agent, named AlphaX, to tackle these two aspects. AlphaX improves the search efficiency by adaptively balancing the exploration and exploitation at the state level, and by a Meta-Deep Neural Network (DNN) to predict network accuracies for biasing the search toward a promising region. To amortize the network evaluation cost, AlphaX accelerates MCTS rollouts with a distributed design and reduces the number of epochs in evaluating a network by transfer learning, which is guided with the tree structure in MCTS. In 12 GPU days and 1000 samples, AlphaX found an architecture that reaches 97.84\% top-1 accuracy on CIFAR-10, and 75.5\% top-1 accuracy on ImageNet, exceeding SOTA NAS methods in both the accuracy and sampling efficiency. Particularly, we also evaluate AlphaX on NASBench-101, a large scale NAS dataset; AlphaX is 3x and 2.8x more sample efficient than Random Search and Regularized Evolution in finding the global optimum. Finally, we show the searched architecture improves a variety of vision applications from Neural Style Transfer, to Image Captioning and Object Detection.


title: Overcoming catastrophic forgetting problem by weight consolidation and long-term memory

abstract: Sequential learning of multiple tasks in artificial neural networks using gradient descent leads to catastrophic forgetting, whereby previously learned knowledge is erased during learning of new, disjoint knowledge. Here, we propose a new approach to sequential learning which leverages the recent discovery of adversarial examples. We use adversarial subspaces from previous tasks to enable learning of new tasks with less interference. We apply our method to sequentially learning to classify digits 0, 1, 2 (task 1), 4, 5, 6, (task 2), and 7, 8, 9 (task 3) in MNIST (disjoint MNIST task). We compare and combine our Adversarial Direction (AD) method with the recently proposed Elastic Weight Consolidation (EWC) method for sequential learning. We train each task for 20 epochs, which yields good initial performance (99.24% correct task 1 performance). After training task 2, and then task 3, both plain gradient descent (PGD) and EWC largely forget task 1 (task 1 accuracy 32.95% for PGD and 41.02% for EWC), while our combined approach (AD+EWC) still achieves 94.53% correct on task 1. We obtain similar results with a much more difficult disjoint CIFAR10 task, which to our knowledge had not been attempted before (70.10% initial task 1 performance, 67.73% after learning tasks 2 and 3 for AD+EWC, while PGD and EWC both fall to chance level). Our results suggest that AD+EWC can provide better sequential learning performance than either PGD or EWC.


title: My camera can see through fences: A deep learning approach for image de-fencing

abstract: In recent times, the availability of inexpensive image capturing devices such as smartphones/tablets has led to an exponential increase in the number of images/videos captured. However, sometimes the amateur photographer is hindered by fences in the scene which have to be removed after the image has been captured. Conventional approaches to image de-fencing suffer from inaccurate and non-robust fence detection apart from being limited to processing images of only static occluded scenes. In this paper, we propose a semi-automated de-fencing algorithm using a video of the dynamic scene. We use convolutional neural networks for detecting fence pixels. We provide qualitative as well as quantitative comparison results with existing lattice detection algorithms on the existing PSU NRT data set and a proposed challenging fenced image dataset. The inverse problem of fence removal is solved using split Bregman technique assuming total variation of the de-fenced image as the regularization constraint.


title: Multi-view Sentence Representation Learning

abstract: Multi-view learning can provide self-supervision when different views are available of the same data. The distributional hypothesis provides another form of useful self-supervision from adjacent sentences which are plentiful in large unlabelled corpora. Motivated by the asymmetry in the two hemispheres of the human brain as well as the observation that different learning architectures tend to emphasise different aspects of sentence meaning, we create a unified multi-view sentence representation learning framework, in which, one view encodes the input sentence with a Recurrent Neural Network (RNN), and the other view encodes it with a simple linear model, and the training objective is to maximise the agreement specified by the adjacent context information between two views. We show that, after training, the vectors produced from our multi-view training provide improved representations over the single-view training, and the combination of different views gives further representational improvement and demonstrates solid transferability on standard downstream tasks.


title: Hausdorff dimension of divergent trajectories on homogeneous space

abstract: For one parameter subgroup action on a finite volume homogeneous space, we consider the set of points admitting divergent on average trajectories. We show that the Hausdorff dimension of this set is strictly less than the manifold dimension of the homogeneous space. As a corollary we know that the Hausdorff dimension of the set of points admitting divergent trajectories is not full, which proves a conjecture of Y. Cheung.


title: DVAE#: Discrete Variational Autoencoders with Relaxed Boltzmann Priors

abstract: Boltzmann machines are powerful distributions that have been shown to be an effective prior over binary latent variables in variational autoencoders (VAEs). However, previous methods for training discrete VAEs have used the evidence lower bound and not the tighter importance-weighted bound. We propose two approaches for relaxing Boltzmann machines to continuous distributions that permit training with importance-weighted bounds. These relaxations are based on generalized overlapping transformations and the Gaussian integral trick. Experiments on the MNIST and OMNIGLOT datasets show that these relaxations outperform previous discrete VAEs with Boltzmann priors. An implementation which reproduces these results is available at https://github.com/QuadrantAI/dvae .


title: Adhesive contact between a rigid body of arbitrary shape and a thin elastic coating

abstract: Application of the principle of energy balance to a rigid indenter in contact with elastic layer on a flat rigid substrate provides a very simple derivation of the detachment criterion which earlier has been obtained by much more complicated asymptotic analysis. This criterion allows calculating the adhesive strength of arbitrary contact of a flat-ended indenter, which occurs to be proportional to the area of the face of the indenter and does not depend on its shape. Similarly, the adhesive contact problem can be easily solved in the case of arbitrary three-dimensional shape.


title: Favorable Collisional Demixing of Ash and Fuel in Magnetized Inertial Fusion

abstract: Magnetized inertial fusion experiments are approaching regimes where the radial transport is dominated by collisions between magnetized ions, providing an opportunity to exploit effects usually associated with steady-state magnetic fusion. In particular, the low-density hotspot characteristic of magnetized liner inertial fusion results in diamagnetic and thermal frictions which can demix thermalized ash from fuel, accelerating the fusion reaction. For reactor regimes in which there is substantial burnup of the fuel, increases in the fusion energy yield on the order of 5% are possible.


title: Closed Walk Sampler: An Efficient Method for Estimating Eigenvalues of Large Graphs

abstract: Eigenvalues of a graph are of high interest in graph analytics for Big Data due to their relevance to many important properties of the graph including network resilience, community detection and the speed of viral propagation. Accurate computation of eigenvalues of extremely large graphs is usually not feasible due to the prohibitive computational and storage costs and also because full access to many social network graphs is often restricted to most researchers. In this paper, we present a series of new sampling algorithms which solve both of the above-mentioned problems and estimate the two largest eigenvalues of a large graph efficiently and with high accuracy. Unlike previous methods which try to extract a subgraph with the most influential nodes, our algorithms sample only a small portion of the large graph via a simple random walk, and arrive at estimates of the two largest eigenvalues by estimating the number of closed walks of a certain length. Our experimental results using real graphs show that our algorithms are substantially faster while also achieving significantly better accuracy on most graphs than the current state-of-the-art algorithms.


title: Odd characteristic classes in entire cyclic homology and equivariant loop space homology

abstract: Given a compact manifold $M$ and $g\in C^{\infty}(M,U(l;\mathbb{C}))$ we construct a Chern character $\mathrm{Ch}^-(g)$ which lives in the odd part of the equivariant (entire) cyclic Chen-normalized bar complex $\underline{\mathscr{C}}(\Omega_{\mathbb{T}}(M\times \mathbb{T}))$ of $M$, and which is mapped to the odd Bismut-Chern character under the equivariant Chen integral map. It is also shown that the assignment $g\mapsto \mathrm{Ch}^-(g)$ induces a well-defined group homomorphism from the $K^{-1}$ theory of $M$ to the odd homology group of $\underline{\mathscr{C}}(\Omega_{\mathbb{T}}(M\times \mathbb{T}))$


title: Efficient Atlasing and Search of Configuration Spaces of Point-Sets Constrained by Distance Intervals

abstract: For configurations of point-sets that are pairwise constrained by distance intervals, the EASAL software implements a suite of algorithms that characterize the structure and geometric properties of the configuration space. The algorithms generate, describe and explore these configuration spaces using generic rigidity properties, classical results for stratification of semi-algebraic sets, and new results for efficient sampling by convex parametrization. The paper reviews the key theoretical underpinnings, major algorithms and their implementation. The paper outlines the main applications such as the computation of free energy and kinetics of assembly of supramolecular structures or of clusters in colloidal and soft materials. In addition, the paper surveys select experimental results and comparisons.


title: Butterfly-Net: Optimal Function Representation Based on Convolutional Neural Networks

abstract: Deep networks, especially convolutional neural networks (CNNs), have been successfully applied in various areas of machine learning as well as to challenging problems in other scientific and engineering fields. This paper introduces Butterfly-Net, a low-complexity CNN with structured and sparse cross-channel connections, together with a Butterfly initialization strategy for a family of networks. Theoretical analysis of the approximation power of Butterfly-Net to the Fourier representation of input data shows that the error decays exponentially as the depth increases. Combining Butterfly-Net with a fully connected neural network, a large class of problems are proved to be well approximated with network complexity depending on the effective frequency bandwidth instead of the input dimension. Regular CNN is covered as a special case in our analysis. Numerical experiments validate the analytical results on the approximation of Fourier kernels and energy functionals of Poisson's equations. Moreover, all experiments support that training from Butterfly initialization outperforms training from random initialization. Also, adding the remaining cross-channel connections, although significantly increase the parameter number, does not much improve the post-training accuracy and is more sensitive to data distribution.


title: A Method for Predicting Nonequilibrium Thermal Expansion Using Steepest-Entropy-Ascent Quantum Thermodynamics

abstract: Steepest-entropy-ascent quantum thermodynamics (SEAQT) is an intriguing approach that describes equilibrium and dynamic processes in a self-consistent way. The applicability is limited to mainly gas phases because of a complex eigenstructure (eigenvalues and eigenfunctions) of solid or liquid phases. In this contribution, the SEAQT modeling is extended to a condensed phase by constructing a simplified eigenstructure (so-called pseudo-eigenstructure), and the applicability is demonstrated by calculating the thermal expansion of metallic silver in three cases: (a) at stable equilibrium, (b) along three irreversible paths from an initial nonequilibrium state to stable equilibrium, and (c) along an irreversible path between two stable equilibrium states. The SEAQT framework with an anharmonic pseudo-eigenstructure predicts reasonable values for equilibrium thermal expansion. For the irreversible cases considered, the SEAQT approach makes it possible to predict the time-dependence of lattice relaxations from the initial state to the final state.


title: The line shapes of the $Z_b(10610)$ and $Z_b(10650)$ in the elastic and inelastic channels revisited

abstract: The most recent experimental data for all measured production and decay channels of the bottomonium-like states $Z_b(10610)$ and $Z_b(10650)$ are analysed simultaneously using solutions of the Lippmann-Schwinger equations which respect constraints from unitarity and analyticity. The interaction potential in the open-bottom channels $B^{(*)}\bar{B}^{*}+\mbox{c.c.}$ contains short-range interactions as well as one-pion exchange. It is found that the long-range interaction does not affect the line shapes as long as only $S$ waves are considered. Meanwhile, the line shapes can be visibly modified once $D$ waves, mediated by the strong tensor forces from the pion exchange potentials, are included. However, in the fit they get balanced largely by a momentum dependent contact term that appears to be needed also to render the results for the line shapes independent of the cut-off. The resulting line shapes are found to be insensitive to various higher-order interactions included to verify stability of the results. Both $Z_b$ states are found to be described by the poles located on the unphysical Riemann sheets in the vicinity of the corresponding thresholds. In particular, the $Z_b(10610)$ state is associated with a virtual state residing just below the $B\bar{B}^{*}/\bar B{B}^{*}$ threshold while the $Z_b(10650)$ state most likely is a shallow state located just above the $B^*\bar{B}^{*}$ threshold.


title: Fisher Efficient Inference of Intractable Models

abstract: Maximum Likelihood Estimators (MLE) has many good properties. For example, the asymptotic variance of MLE solution attains equality of the asymptotic Cram{\'e}r-Rao lower bound (efficiency bound), which is the minimum possible variance for an unbiased estimator. However, obtaining such MLE solution requires calculating the likelihood function which may not be tractable due to the normalization term of the density model. In this paper, we derive a Discriminative Likelihood Estimator (DLE) from the Kullback-Leibler divergence minimization criterion implemented via density ratio estimation and a Stein operator. We study the problem of model inference using DLE. We prove its consistency and show that the asymptotic variance of its solution can attain the equality of the efficiency bound under mild regularity conditions. We also propose a dual formulation of DLE which can be easily optimized. Numerical studies validate our asymptotic theorems and we give an example where DLE successfully estimates an intractable model constructed using a pre-trained deep neural network.


title: Subspace Selection via DR-Submodular Maximization on Lattices

abstract: The subspace selection problem seeks a subspace that maximizes an objective function under some constraint. This problem includes several important machine learning problems such as the principal component analysis and sparse dictionary selection problem. Often, these problems can be solved by greedy algorithms. Here, we are interested in why these problems can be solved by greedy algorithms, and what classes of objective functions and constraints admit this property. To answer this question, we formulate the problems as optimization problems on lattices. Then, we introduce a new class of functions, directional DR-submodular functions, to characterize the approximability of problems. We see that the principal component analysis, sparse dictionary selection problem, and these generalizations have directional DR-submodularities. We show that, under several constraints, the directional DR-submodular function maximization problem can be solved efficiently with provable approximation factors.