title: Reduced Order Optimal Control of the Convective FitzHugh-Nagumo Equation

abstract: In this paper, we compare three model order reduction methods: the proper orthogonal decomposition (POD), discrete empirical interpolation method (DEIM) and dynamic mode decomposition (DMD) for the optimal control of the convective FitzHugh-Nagumo (FHN) equations. The convective FHN equations consists of the semi-linear activator and the linear inhibitor equations, modeling blood coagulation in moving excitable media. The semilinear activator equation leads to a non-convex optimal control problem (OCP). The most commonly used method in reduced optimal control is POD. We use DEIM and DMD to approximate efficiently the nonlinear terms in reduced order models. We compare the accuracy and computational times of three reduced-order optimal control solutions with the full order discontinuous Galerkin finite element solution of the convection dominated FHN equations with terminal controls. Numerical results show that POD is the most accurate whereas POD-DMD is the fastest.


title: Nonlinear Model and its Inverse of an Audio System

abstract: This computer science master thesis aims at modelling the nonlinearities of a loudspeaker. A piecewise linear approximation is initially explored and then we present a nonlinear Volterra model to simulate the behavior of the system. The general theory of continuous and discrete Volterra series is summarised. A Normalized Least Mean Square algorithm is used to determine the Volterra series to third order. We also present as inverted system which is trained with the same algorithm. Training data for the models were collected measuring a physical speaker using a laser interferometer. Results indicate a decrease in Mean Squared Error compared to the linear model with a dependency on the particular test signal, the order and the parameters of the model.


title: Offsets between member galaxies and dark matter in clusters: a test with the Illustris simulation

abstract: Dark matter with a non-zero self-interacting cross section ($\sigma_{\rm SIDM}$) has been posited as a solution to a number of outstanding astrophysical mysteries. Many studies of merging galaxy clusters have given constraints on $\sigma_{\rm SIDM}$ based on the spatial offset between the member galaxy population and the dark matter distribution. Assuming $\sigma_{\rm SIDM} = 0$, how likely is it for us to see the galaxy-DM offset values observed in merging clusters of galaxies? To answer this question, we formulate a hypothesis test using data from Illustris, a $\Lambda$CDM cosmological simulation. We select 43 Illustris clusters and their galaxy members at z~0 and examine the accuracy of commonly used galaxy summary statistics, including kernel-density-estimation (KDE) luminosity peak, KDE number density peak, shrinking aperture, centroid and the location of the brightest cluster galaxy (BCG). We use the dark-matter particles to reproduce commonly adopted methods to identify dark-matter peaks based on gravitational lensing cluster maps. By analysing each cluster in 768 projections, we determine the optimistic noise floor in the measurements of the galaxy-DM offsets. We find that the choice of the galaxy summary statistics affects the inferred offset values substantially, with the BCG and the luminosity peak giving the tightest 68-th percentile offset levels, $\lesssim$ 4 kpc and $\lesssim$ 32 kpc, respectively. Shrinking aperture, number density and centroid give a large offset scatter of about 50-100 kpc at the 68-th percentile level, even for clusters with only one dominant mass component. Out of the 15 reported offsets from observed merging clusters that we examined, 13 of them are consistent with Illustris unrelaxed cluster offsets at the 2-sigma (95-th percentile) level, i.e. consistent with the hypothesis that $\Lambda$CDM is the true underlying physical model.


title: 1 to 2.4 micron Near-IR spectrum of the Giant Planet $\beta$ Pictoris b obtained with the Gemini Planet Imager

abstract: Using the Gemini Planet Imager (GPI) located at Gemini South, we measured the near-infrared (1.0-2.4 micron) spectrum of the planetary companion to the nearby, young star $\beta$ Pictoris. We compare the spectrum obtained with currently published model grids and with known substellar objects and present the best matching models as well as the best matching observed objects. Comparing the empirical measurement of the bolometric luminosity to evolutionary models, we find a mass of $12.9\pm0.2$ $\mathcal{M}_\mathrm{Jup}$, an effective temperature of $1724\pm15$ K, a radius of $1.46\pm0.01$ $\mathcal{R}_\mathrm{Jup}$, and a surface gravity of $\log g = 4.18\pm0.01$ [dex] (cgs). The stated uncertainties are statistical errors only, and do not incorporate any uncertainty on the evolutionary models. Using atmospheric models, we find an effective temperature of $1700-1800$ K and a surface gravity of $\log g = 3.5$-$4.0$ [dex] depending upon model. These values agree well with other publications and with "hot-start" predictions from planetary evolution models. Further, we find that the spectrum of $\beta$ Pic b best matches a low-surface gravity L2$\pm$1 brown dwarf. Finally comparing the spectrum to field brown dwarfs we find the the spectrum best matches 2MASS J04062677-381210 and 2MASS J03552337+1133437.


title: Accounting for Chromatic Atmospheric Effects on Barycentric Corrections

abstract: Atmospheric effects on stellar radial velocity measurements for exoplanet discovery and characterization have not yet been fully investigated for extreme precision levels. We carry out calculations to determine the wavelength dependence of barycentric corrections across optical wavelengths, due to the ubiquitous variations in air mass during observations. We demonstrate that radial velocity errors of at least several cm/s can be incurred if the wavelength dependence is not included in the photon-weighted barycentric corrections. A minimum of four wavelength channels across optical spectra (380-680 nm) are required to account for this effect at the 10 cm/s level, with polynomial fits of the barycentric corrections applied to cover all wavelengths. Additional channels may be required in poor observing conditions or to avoid strong telluric absorption features. Furthermore, consistent flux sampling on the order of seconds throughout the observation is necessary to ensure that accurate photon weights are obtained. Finally, we describe how a multiple-channel exposure meter will be implemented in the EXtreme PREcision Spectrograph (EXPRES).


title: How Dark Matter Came to Matter

abstract: The history of the dark matter problem can be traced back to at least the 1930s, but it was not until the early 1970s that the issue of 'missing matter' was widely recognized as problematic. In the latter period, previously separate issues involving missing mass were brought together in a single anomaly. We argue that reference to a straightforward 'accumulation of evidence' alone is inadequate to comprehend this episode. Rather, the rise of cosmological research, the accompanying renewed interest in the theory of relativity and changes in the manpower division of astronomy in the 1960s are key to understanding how dark matter came to matter. At the same time, this story may also enlighten us on the methodological dimensions of past practices of physics and cosmology.

categories: astro-ph.CO gr-qc hep-ph physics.hist-ph


title: Insulators and metals with topological order and discrete symmetry breaking

abstract: Numerous experiments have reported discrete symmetry breaking in the high temperature pseudogap phase of the hole-doped cuprates, including breaking of one or more of lattice rotation, inversion, or time-reversal symmetries. In the absence of translational symmetry breaking or topological order, these conventional order parameters cannot explain the gap in the charged fermion excitation spectrum in the anti-nodal region. Zhao et al. (1601.01688) and Jeong et al. (arXiv:1701.06485) have also reported inversion and time-reversal symmetry breaking in insulating Sr2IrO4 similar to that in the metallic cuprates, but co-existing with Neel order. We extend an earlier theory of topological order in insulators and metals, in which the topological order combines naturally with the breaking of these conventional discrete symmetries. We find translationally-invariant states with topological order co-existing with both Ising-nematic order and spontaneous charge currents. The link between the discrete broken symmetries and the topological-order-induced pseudogap explains why the broken symmetries do not survive in the confining phases without a pseudogap at large doping. Our theory also connects to the O(3) non-linear sigma model and CP1 descriptions of quantum fluctuations of the Neel order. In this framework, the optimal doping criticality of the cuprates is primarily associated with the loss of topological order.


title: A stellar census of the nearby, young 32 Orionis group

abstract: The 32 Orionis group was discovered almost a decade ago and despite the fact that it represents the first northern, young (age ~ 25 Myr) stellar aggregate within 100 pc of the Sun ($d \simeq 93$ pc), a comprehensive survey for members and detailed characterisation of the group has yet to be performed. We present the first large-scale spectroscopic survey for new (predominantly M-type) members of the group after combining kinematic and photometric data to select candidates with Galactic space motion and positions in colour-magnitude space consistent with membership. We identify 30 new members, increasing the number of known 32 Ori group members by a factor of three and bringing the total number of identified members to 46, spanning spectral types B5 to L1. We also identify the lithium depletion boundary (LDB) of the group, i.e. the luminosity at which lithium remains unburnt in a coeval population. We estimate the age of the 32 Ori group independently using both isochronal fitting and LDB analyses and find it is essentially coeval with the {\beta} Pictoris moving group, with an age of $24\pm4$ Myr. Finally, we have also searched for circumstellar disc hosts utilising the AllWISE catalogue. Although we find no evidence for warm, dusty discs, we identify several stars with excess emission in the WISE W4-band at 22 {\mu}m. Based on the limited number of W4 detections we estimate a debris disc fraction of $32^{+12}_{-8}$ per cent for the 32 Ori group.


title: ZFIRE: The Evolution of the Stellar Mass Tully-Fisher Relation to Redshift 2.0 < Z < 2.5 with MOSFIRE

abstract: Using observations made with MOSFIRE on Keck I as part of the ZFIRE survey, we present the stellar mass Tully-Fisher relation at 2.0 < z < 2.5. The sample was drawn from a stellar mass limited, Ks-band selected catalog from ZFOURGE over the CANDELS area in the COSMOS field. We model the shear of the Halpha emission line to derive rotational velocities at 2.2X the scale radius of an exponential disk (V2.2). We correct for the blurring effect of a two-dimensional PSF and the fact that the MOSFIRE PSF is better approximated by a Moffat than a Gaussian, which is more typically assumed for natural seeing. We find for the Tully-Fisher relation at 2.0 < z < 2.5 that logV2.2 =(2.18 +/- 0.051)+(0.193 +/- 0.108)(logM/Msun - 10) and infer an evolution of the zeropoint of Delta M/Msun = -0.25 +/- 0.16 dex or Delta M/Msun = -0.39 +/- 0.21 dex compared to z = 0 when adopting a fixed slope of 0.29 or 1/4.5, respectively. We also derive the alternative kinematic estimator S0.5, with a best-fit relation logS0.5 =(2.06 +/- 0.032)+(0.211 +/- 0.086)(logM/Msun - 10), and infer an evolution of Delta M/Msun= -0.45 +/- 0.13 dex compared to z < 1.2 if we adopt a fixed slope. We investigate and review various systematics, ranging from PSF effects, projection effects, systematics related to stellar mass derivation, selection biases and slope. We find that discrepancies between the various literature values are reduced when taking these into account. Our observations correspond well with the gradual evolution predicted by semi-analytic models.


title: Broken boost invariance in the Glasma via finite nuclei thickness

abstract: We simulate the creation and evolution of non-boost-invariant Glasma in the early stages of heavy ion collisions within the color glass condensate framework. This is accomplished by extending the McLerran-Venugopalan model to include a parameter for the Lorentz-contracted but finite width of the nucleus in the beam direction. We determine the rapidity profile of the Glasma energy density, which shows deviations from the boost-invariant result. Varying the parameters both broad and narrow profiles can be produced. We compare our results to experimental data from RHIC and find surprising agreement.


title: Bulk Connectedness and Boundary Entanglement

abstract: We prove, for any state in a conformal field theory defined on a set of boundary manifolds with corresponding classical holographic bulk geometry, that for any bipartition of the boundary into two non-clopen sets, the density matrix cannot be a tensor product of the reduced density matrices on each region of the bipartition. In particular, there must be entanglement across the bipartition surface. We extend this no-go theorem to general, arbitrary partitions of the boundary manifolds into non-clopen parts, proving that the density matrix cannot be a tensor product. This result gives a necessary condition for states to potentially correspond to holographic duals.


title: A Differential Abundance Analysis of Very Metal-Poor Stars

abstract: We have performed a differential, line-by-line, chemical abundance analysis, ultimately relative to the Sun, of nine very metal-poor main sequence halo stars, near [Fe/H]=$-$2 dex. Our abundances range from $-2.66\leq\mathrm{[Fe/H]}\leq-1.40$ dex with conservative uncertainties of 0.07 dex. We find an average [$\alpha$/Fe]$=0.34\pm0.09$ dex, typical of the Milky Way. While our spectroscopic atmosphere parameters provide good agreement with HST parallaxes, there is significant disagreement with temperature and gravity parameters indicated by observed colors and theoretical isochrones. Although a systematic underestimate of the stellar temperature by a few hundred degrees could explain this difference, it is not supported by current effective temperature studies and would create large uncertainties in the abundance determinations. Both 1D and $\langle$3D$\rangle$ hydrodynamical models combined with separate 1D non-LTE effects do not yet account for the atmospheres of real metal-poor MS stars, but a fully 3D non-LTE treatment may be able to explain the ionization imbalance found in this work.


title: $\texttt{galkin}$: a new compilation of the Milky Way rotation curve data

abstract: We present $\texttt{galkin}$, a novel compilation of kinematic measurements tracing the rotation curve of our Galaxy, together with a tool to treat the data. The compilation is optimised to Galactocentric radii between 3 and 20 kpc and includes the kinematics of gas, stars and masers in a total of 2780 measurements carefully collected from almost four decades of literature. A simple, user-friendly tool is provided to select, treat and retrieve the data of all source references considered. This tool is especially designed to facilitate the use of kinematic data in dynamical studies of the Milky Way with various applications ranging from dark matter constraints to tests of modified gravity.


title: The NuSTAR Hard X-ray Survey of the Norma Arm Region

abstract: We present a catalog of hard X-ray sources in a square-degree region surveyed by NuSTAR in the direction of the Norma spiral arm. This survey has a total exposure time of 1.7 Ms, and typical and maximum exposure depths of 50 ks and 1 Ms, respectively. In the area of deepest coverage, sensitivity limits of $5\times10^{-14}$ and $4\times10^{-14}$ erg s$^{-1}$ cm$^{-2}$ in the 3-10 and 10-20 keV bands, respectively, are reached. Twenty-eight sources are firmly detected and ten are detected with low significance; eight of the 38 sources are expected to be active galactic nuclei. The three brightest sources were previously identified as a low-mass X-ray binary, high-mass X-ray binary, and pulsar wind nebula. Based on their X-ray properties and multi-wavelength counterparts, we identify the likely nature of the other sources as two colliding wind binaries, three pulsar wind nebulae, a black hole binary, and a plurality of cataclysmic variables (CVs). The CV candidates in the Norma region have plasma temperatures of $\approx$10-20 keV, consistent with the Galactic Ridge X-ray emission spectrum but lower than temperatures of CVs near the Galactic Center. This temperature difference may indicate that the Norma region has a lower fraction of intermediate polars relative to other types of CVs compared to the Galactic Center. The NuSTAR log$N$-log$S$ distribution in the 10-20 keV band is consistent with the distribution measured by Chandra at 2-10 keV if the average source spectrum is assumed to be a thermal model with $kT\approx15$~keV, as observed for the CV candidates.


title: Self-similar accretion in thin disks around near-extremal black holes

abstract: Near-maximally spinning black holes display conformal symmetry in their near-horizon region, which is therefore the locus of critical phenomena. In this paper, we revisit the Novikov-Thorne accretion thin disk model and find a new self-similar radiation-dominated solution in the extremely high spin regime. Motivated by the self-consistency of the model, we require that matter flows at the sound speed at the innermost stable circular orbit (ISCO). We observe that, when the disk pressure is dominated by radiation at the ISCO, which occurs for the best-fitting Novikov-Thorne model of GRS 1915+105, the Shakura-Sunyaev viscosity parameter can be expressed in terms of the spin, mass accretion rate and radiative efficiency. We quantitatively describe how the exact thin disk solution approaches the self-similar solution in the vicinity of the ISCO and for increasing spins.


title: Stellar disruption events support the existence of the black hole event horizon

abstract: Many black hole (BH) candidates have been discovered in X-ray binaries and in the nuclei of galaxies. The prediction of Einstein's general relativity is that BHs have an event horizon --- a one-way membrane through which particles fall into the BH but cannot exit. However, except for the very few nearby supermassive BH candidates, our telescopes are unable to resolve and provide a direct proof of the event horizon. Here, we propose a novel observation that supports the existence of event horizons around supermassive BH candidates heavier than 10^{7.5} solar masses. Instead of an event horizon, if the BH candidate has a hard surface, when a star falls onto the surface, the shocked baryonic gas will form a radiation pressure supported envelope that shines at the Eddington luminosity for an extended period of time from months to years. We show that such emission has already been ruled out by the Pan-STARRS1 3pi survey if supermassive BH candidates have a hard surface at radius larger than 1 + 10^{-4.4} times the Schwarzschild radius. Future observations by LSST should be able to improve the limit to 1 + 10^{-6}.


title: Subleading Soft Graviton Theorem for Loop Amplitudes

abstract: Superstring field theory gives expressions for heterotic and type II string loop amplitudes that are free from ultraviolet and infrared divergences when the number of non-compact space-time dimensions is five or more. We prove the subleading soft graviton theorem in these theories to all orders in perturbation theory for S-matrix elements of arbitrary number of finite energy external states but only one external soft graviton. We also prove the leading soft graviton theorem for arbitrary number of finite energy external states and arbitrary number of soft gravitons. Since our analysis is based on general properties of one particle irreducible effective action, the results are valid in any theory of quantum gravity that gives finite result for the S-matrix order by order in perturbation theory without violating general coordinate invariance.


title: Unitary null energy condition violation in $P(X)$ cosmologies

abstract: A non-singular cosmological bounce in the Einstein frame can only take place if the Null Energy Condition (NEC) is violated. We explore situations where a single scalar field drives the NEC violation and derive the constraints imposed by demanding tree level unitarity on a cosmological background. We then focus on the explicit constraints that arise in P(X) theories and show that constraints from perturbative unitarity make it impossible for the NEC violation to occur within the region of validity of the effective field theory without also involving irrelevant operators that arise at a higher scale that would enter from integrating out more massive degrees of freedom. Within the context of P(X) theories we show that including such operators allows for a bounce that does not manifestly violate tree level unitarity, but at the price of either imposing a shift symmetry or involving technically unnatural small operator coefficients within the low-energy effective field theory.


title: Nanopteron solutions of diatomic Fermi-Pasta-Ulam-Tsingou lattices with small mass-ratio

abstract: Consider an infinite chain of masses, each connected to its nearest neighbors by a (nonlinear) spring. This is a Fermi-Pasta-Ulam-Tsingou lattice. We prove the existence of traveling waves in the setting where the masses alternate in size. In particular we address the limit where the mass ratio tends to zero. The problem is inherently singular and we find that the traveling waves are not true solitary waves but rather "nanopterons", which is to say, waves which asymptotic at spatial infinity to very small amplitude periodic waves. Moreover, we can only find solutions when the mass ratio lies in a certain open set. The difficulties in the problem all revolve around understanding Jost solutions of a nonlocal Schr\"odinger operator in its semi-classical limit.


title: Decidability and Independence of Conjugacy Problems in Finitely Presented Monoids

abstract: There have been several attempts to extend the notion of conjugacy from groups to monoids. The aim of this paper is study the decidability and independence of conjugacy problems for three of these notions (which we will denote by $\sim_p$, $\sim_o$, and $\sim_c$) in certain classes of finitely presented monoids. We will show that in the class of polycyclic monoids, $p$-conjugacy is "almost" transitive, $\sim_c$ is strictly included in $\sim_p$, and the $p$- and $c$-conjugacy problems are decidable with linear compexity. For other classes of monoids, the situation is more complicated. We show that there exists a monoid $M$ defined by a finite complete presentation such that the $c$-conjugacy problem for $M$ is undecidable, and that for finitely presented monoids, the $c$-conjugacy problem and the word problem are independent, as are the $c$-conjugacy and $p$-conjugacy problems.


title: Excess entropy and energy feedback from within cluster cores up to r$_{200}$

abstract: We estimate the "non-gravitational" entropy-injection profiles, $\Delta K$, and the resultant energy feedback profiles, $\Delta E$, of the intracluster medium for 17 clusters using their Planck SZ and ROSAT X-Ray observations, spanning a large radial range from $0.2r_{500}$ up to $r_{200}$. The feedback profiles are estimated by comparing the observed entropy, at fixed gas mass shells, with theoretical entropy profiles predicted from non-radiative hydrodynamic simulations. We include non-thermal pressure and gas clumping in our analysis. The inclusion of non-thermal pressure and clumping results in changing the estimates for $r_{500}$ and $r_{200}$ by 10\%-20\%. When clumpiness is not considered it leads to an under-estimation of $\Delta K\approx300$ keV cm$^2$ at $r_{500}$ and $\Delta K\approx1100$ keV cm$^2$ at $r_{200}$. On the other hand, neglecting non-thermal pressure results in an over-estimation of $\Delta K\approx 100$ keV cm$^2$ at $r_{500}$ and under-estimation of $\Delta K\approx450$ keV cm$^2$ at $r_{200}$. For the estimated feedback energy, we find that ignoring clumping leads to an under-estimation of energy per particle $\Delta E\approx1$ keV at $r_{500}$ and $\Delta E\approx1.5$ keV at $r_{200}$. Similarly, neglect of the non-thermal pressure results in an over-estimation of $\Delta E\approx0.5$ keV at $r_{500}$ and under-estimation of $\Delta E\approx0.25$ keV at $r_{200}$. We find entropy floor of $\Delta K\approx300$ keV cm$^2$ is ruled out at $\approx3\sigma$ throughout the entire radial range and $\Delta E\approx1$ keV at more than 3$\sigma$ beyond $r_{500}$, strongly constraining ICM pre-heating scenarios. We also demonstrate robustness of results w.r.t sample selection, X-Ray analysis procedures, entropy modeling etc.


title: Evolutions of Gowdy, Brill and Teukolsky initial data on a smooth lattice

abstract: Numerical results, based on a lattice method for computational general relativity, will be presented for Cauchy evolution of initial data for the Brill, Teukolsky and polarised Gowdy space-times. The simple objective of this paper is to demonstrate that the lattice method can, at least for these space-times, match results obtained from contemporary methods. Some of the issues addressed in this paper include the handling of axisymmetric instabilities (in the Brill space-time) and an implementation of a Sommerfeld radiation condition for the Brill and Teukolsky space-times. It will be shown that the lattice method performs particularly well in regard to the passage of the waves through the outer boundary. Questions concerning multiple black-holes, mesh refinement and long term stability will not be discussed here but may form the basis of future work.


title: Image Subtraction Reduction of Open Clusters M35 & NGC 2158 In The K2 Campaign-0 Super-Stamp

abstract: Observations were made of the open clusters M35 and NGC 2158 during the initial K2 campaign (C0). Reducing these data to high-precision photometric time-series is challenging due to the wide point spread function (PSF) and the blending of stellar light in such dense regions. We developed an image-subtraction-based K2 reduction pipeline that is applicable to both crowded and sparse stellar fields. We applied our pipeline to the data-rich C0 K2 super-stamp, containing the two open clusters, as well as to the neighboring postage stamps. In this paper, we present our image subtraction reduction pipeline and demonstrate that this technique achieves ultra-high photometric precision for sources in the C0 super-stamp. We extract the raw light curves of 3960 stars taken from the UCAC4 and EPIC catalogs and de-trend them for systematic effects. We compare our photometric results with the prior reductions published in the literature. For detrended, TFA-corrected sources in the 12--12.25 $\rm K_{p}$ magnitude range, we achieve a best 6.5 hour window running rms of 35 ppm falling to 100 ppm for fainter stars in the 14--14.25 $ \rm K_{p}$ magnitude range. For stars with $\rm K_{p}> 14$, our detrended and 6.5 hour binned light curves achieve the highest photometric precision. Moreover, all our TFA-corrected sources have higher precision on all time scales investigated. This work represents the first published image subtraction analysis of a K2 super-stamp. This method will be particularly useful for analyzing the Galactic bulge observations carried out during K2 campaign 9. The raw light curves and the final results of our detrending processes are publicly available at \url{http://k2.hatsurveys.org/archive/}.


title: MPC Validation and Aggregation of Unit Vectors

abstract: When dealing with privatized data, it is important to be able to protect against malformed user inputs. This becomes difficult in MPC systems as each server should not contain enough information to know what values any user has submitted. In this paper, we implement an MPC technique to verify blinded user inputs are unit vectors. In addition, we introduce a BGW circuit which can securely aggregate the blinded inputs while only releasing the result when it is above a public threshold. These distributed techniques take as input a unit vector. While this initially seems limiting compared to real number input, it is quite powerful for cases such as selecting from a list of options, indicating a location from a set of possibilities, or any system which uses one-hot encoding.


title: Noise-resilient preparation of quantum many-body ground states

abstract: Certain quantum many-body ground states can be prepared by a large-depth quantum circuit consisting of geometrically local gates. In the presence of noise, local expectation values deviate from the correct value at most by an amount comparable to the noise rate. This happens if the action of the noiseless circuit, restricted to certain subsystems, rapidly mixes local observables up to a small correction. The encoding circuit for the surface code is given as an example.


title: Force-gradient sensitive Kelvin probe force microscopy by dissipative electrostatic force modulation

abstract: We report a Kelvin probe force microscopy (KPFM) implementation using the dissipation signal of a frequency modulation atomic force microscopy that is capable of detecting the gradient of electrostatic force rather than electrostatic force. It features a simple implementation and faster scanning as it requires no low frequency modulation. We show that applying a coherent ac voltage with two times the cantilever oscillation frequency induces the dissipation signal proportional to the electrostatic force gradient which depends on the effective dc bias voltage including the contact potential difference. We demonstrate the KPFM images of a MoS$_2$ flake taken with the present method is in quantitative agreement with that taken with the frequency modulated Kelvin probe force microscopy technique.


title: Weighted Random Walk Sampling for Multi-Relational Recommendation

abstract: In the information overloaded web, personalized recommender systems are essential tools to help users find most relevant information. The most heavily-used recommendation frameworks assume user interactions that are characterized by a single relation. However, for many tasks, such as recommendation in social networks, user-item interactions must be modeled as a complex network of multiple relations, not only a single relation. Recently research on multi-relational factorization and hybrid recommender models has shown that using extended meta-paths to capture additional information about both users and items in the network can enhance the accuracy of recommendations in such networks. Most of this work is focused on unweighted heterogeneous networks, and to apply these techniques, weighted relations must be simplified into binary ones. However, information associated with weighted edges, such as user ratings, which may be crucial for recommendation, are lost in such binarization. In this paper, we explore a random walk sampling method in which the frequency of edge sampling is a function of edge weight, and apply this generate extended meta-paths in weighted heterogeneous networks. With this sampling technique, we demonstrate improved performance on multiple data sets both in terms of recommendation accuracy and model generation efficiency.


title: Context-Sensitive Super-Resolution for Fast Fetal Magnetic Resonance Imaging

abstract: 3D Magnetic Resonance Imaging (MRI) is often a trade-off between fast but low-resolution image acquisition and highly detailed but slow image acquisition. Fast imaging is required for targets that move to avoid motion artefacts. This is in particular difficult for fetal MRI. Spatially independent upsampling techniques, which are the state-of-the-art to address this problem, are error prone and disregard contextual information. In this paper we propose a context-sensitive upsampling method based on a residual convolutional neural network model that learns organ specific appearance and adopts semantically to input data allowing for the generation of high resolution images with sharp edges and fine scale detail. By making contextual decisions about appearance and shape, present in different parts of an image, we gain a maximum of structural detail at a similar contrast as provided by high-resolution data. We experiment on $145$ fetal scans and show that our approach yields an increased PSNR of $1.25$ $dB$ when applied to under-sampled fetal data \emph{cf.} baseline upsampling. Furthermore, our method yields an increased PSNR of $1.73$ $dB$ when utilizing under-sampled fetal data to perform brain volume reconstruction on motion corrupted captured data.


title: Huygens' principle and Dirac-Weyl equation

abstract: We investigate the validity of Huygens' principle for forward propagation in the massless Dirac-Weyl equation. The principle holds for odd space dimension n, while it is invalid for even n. We explicitly solve the cases n=1,2 and 3 and discuss generic $n$. We compare with the massless Klein-Gordon equation and comment on possible generalizations and applications.


title: Managing the Public to Manage Data: Citizen Science and Astronomy

abstract: Citizen science projects recruit members of the public as volunteers to process and produce datasets. These datasets must win the trust of the scientific community. The task of securing credibility involves, in part, applying standard scientific procedures to clean these datasets. However, effective management of volunteer behavior also makes a significant contribution to enhancing data quality. Through a case study of Galaxy Zoo, a citizen science project set up to generate datasets based on volunteer classifications of galaxy morphologies, this paper explores how those involved in running the project manage volunteers. The paper focuses on how methods for crediting volunteer contributions motivate volunteers to provide higher quality contributions and to behave in a way that better corresponds to statistical assumptions made when combining volunteer contributions into datasets. These methods have made a significant contribution to the success of the project in securing trust in these datasets, which have been well used by other scientists. Implications for practice are then presented for citizen science projects, providing a list of considerations to guide choices regarding how to credit volunteer contributions to improve the quality and trustworthiness of citizen science-produced datasets.


title: Growth of values of binary quadratic forms and Conway rivers

abstract: We study the growth of the values of binary quadratic forms $Q$ on a binary planar tree as it was described by Conway. We show that the corresponding Lyapunov exponents $\Lambda_Q(x)$ as a function of the path determined by $x\in \mathbb RP^1$ are twice the values of the corresponding exponents for the growth of Markov numbers \cite{SV}, except for the paths corresponding to the Conway rivers, when $\Lambda_Q(x)=0.$ The relation with Galois results about continued fraction expansions for quadratic irrationals is explained and interpreted geometrically.


title: A description length approach to determining the number of k-means clusters

abstract: We present an asymptotic criterion to determine the optimal number of clusters in k-means. We consider k-means as data compression, and propose to adopt the number of clusters that minimizes the estimated description length after compression. Here we report two types of compression ratio based on two ways to quantify the description length of data after compression. This approach further offers a way to evaluate whether clusters obtained with k-means have a hierarchical structure by examining whether multi-stage compression can further reduce the description length. We applied our criteria to determine the number of clusters to synthetic data and empirical neuroimaging data to observe the behavior of the criteria across different types of data set and suitability of the two types of criteria for different datasets. We found that our method can offer reasonable clustering results that are useful for dimension reduction. While our numerical results revealed dependency of our criteria on the various aspects of dataset such as the dimensionality, the description length approach proposed here provides a useful guidance to determine the number of clusters in a principled manner when underlying properties of the data are unknown and only inferred from observation of data.


title: Revisiting universality of the liquid-gas critical point in 2D

abstract: Critical point of liquid-gas (LG) transition does not conform with the paradigm of spontaneous symmetry breaking because there is no broken symmetry in both phases. This stimulated the ongoing debate about the nature of the universality class of the transition -- lasting since the creation of the theory of scaling. We revisit the conjecture that the LG criticality is that of the Ising model. Large scale Monte Carlo simulations of the LG criticality in 2D free space in combination with the numerical flowgram method give the critical indices agreeing with the Onsager values within the error of 1\%. The related problem about the role of higher order odd terms in the (real) $\varphi^4$ field model is addressed too. The scaling dimension of the $\varphi^5$ term at criticality is shown to be the same as of the linear one $\varphi$. We suggest that the role of all higher order odd terms is simply in generating the linear field operator at the criticality.


title: The Schubert normal form of a 3-bridge link and the 3-bridge link group

abstract: We introduce the Schubert form a $3$-bridge link diagram, as a generalization of the Schubert normal form of a $3$-bridge link. It consists of a set of six positive integers, written as $\left( p/n,q/m,s/l\right) $, with some conditions and it is based on the concept of $3$-butterfly. Using the Schubert normal form of a $3$-bridge link diagram, we give two presentations of the 3-bridge link group. These presentations are given by concrete formulas that depend on the integers $\left\{ p,n,q,m,s,l\right\} .$ The construction is a generalization of the form the link group presentation of the $2$-bridge link $p/q$ depends on the integers $p$ and $q$.


title: Reproducible experiments on dynamic resource allocation in cloud data centers

abstract: In Wolke et al. [1] we compare the efficiency of different resource allocation strategies experimentally. We focused on dynamic environments where virtual machines need to be allocated and deallocated to servers over time. In this companion paper, we describe the simulation framework and how to run simulations to replicate experiments or run new experiments within the framework.


title: Tree tribes and lower bounds for switching lemmas

abstract: We show tight upper and lower bounds for switching lemmas obtained by the action of random $p$-restrictions on boolean functions that can be expressed as decision trees in which every vertex is at a distance of at most $t$ from some leaf, also called $t$-clipped decision trees. More specifically, we show the following: $\bullet$ If a boolean function $f$ can be expressed as a $t$-clipped decision tree, then under the action of a random $p$-restriction $\rho$, the probability that the smallest depth decision tree for $f|_{\rho}$ has depth greater than $d$ is upper bounded by $(4p2^{t})^{d}$. $\bullet$ For every $t$, there exists a function $g_{t}$ that can be expressed as a $t$-clipped decision tree, such that under the action of a random $p$-restriction $\rho$, the probability that the smallest depth decision tree for $g_{t}|_{\rho}$ has depth greater than $d$ is lower bounded by $(c_{0}p2^{t})^{d}$, for $0\leq p\leq c_{p}2^{-t}$ and $0\leq d\leq c_{d}\frac{\log n}{2^{t}\log t}$, where $c_{0},c_{p},c_{d}$ are universal constants.


title: Digital Cellular Solid Pressure Vessels: A Novel Approach for Human Habitation in Space

abstract: It is widely assumed that human exploration beyond Earth's orbit will require vehicles capable of providing long-duration habitats that simulate an Earthlike environment: consistent artificial gravity, breathable atmosphere, and sufficient living space- while requiring the minimum possible launch mass. This paper examines how the qualities of digital cellular solids - high-performance, repairability, reconfigurability, tunable mechanical response - allow the accomplishment of long-duration habitat objectives at a fraction of the mass required for traditional structural technologies. To illustrate the impact digital cellular solids could make as a replacement to conventional habitat subsystems, we compare recent proposed deep space habitat structural systems with a digital cellular solids pressure vessel design that consists of a carbon fiber reinforced polymer (CFRP) digital cellular solid cylindrical framework that is lined with an ultra-high molecular weight polyethylene (UHMWPE) skin. We use the analytical treatment of a linear specific modulus scaling cellular solid to find the minimum mass pressure vessel for a structure and find that, for equivalent habitable volume and appropriate safety factors, the use of digital cellular solids provides clear methods for producing structures that are not only repairable and reconfigurable, but also higher performance than their conventionally-manufactured counterparts.


title: Aggregated knowledge from a small number of debates outperforms the wisdom of large crowds

abstract: The aggregation of many independent estimates can outperform the most accurate individual judgment. This centenarian finding, popularly known as the wisdom of crowds, has been applied to problems ranging from the diagnosis of cancer to financial forecasting. It is widely believed that social influence undermines collective wisdom by reducing the diversity of opinions within the crowd. Here, we show that if a large crowd is structured in small independent groups, deliberation and social influence within groups improve the crowd's collective accuracy. We asked a live crowd (N=5180) to respond to general-knowledge questions (e.g., what is the height of the Eiffel Tower?). Participants first answered individually, then deliberated and made consensus decisions in groups of five, and finally provided revised individual estimates. We found that averaging consensus decisions was substantially more accurate than aggregating the initial independent opinions. Remarkably, combining as few as four consensus choices outperformed the wisdom of thousands of individuals.


title: The Malgrange Form and Fredholm Determinants

abstract: We consider the factorization problem of matrix symbols relative to a closed contour, i.e., a Riemann-Hilbert problem, where the symbol depends analytically on parameters. We show how to define a function $\tau$ which is locally analytic on the space of deformations and that is expressed as a Fredholm determinant of an operator of "integrable" type in the sense of Its-Izergin-Korepin-Slavnov. The construction is not unique and the non-uniqueness highlights the fact that the tau function is really the section of a line bundle.


title: Motion of Massive Particles in Rindler Space and the Problem of Fall at the Centre

abstract: The motion of a massive particle in Rindler space has been studied and obtained the geodesics of motion. The orbits in Rindler space are found to be quite different from that of Schwarzschild case. The paths are not like the Perihelion Precession type. Further we have set up the non-relativistic Schrodinger equation for the particle in the quantum mechanical scenario in presence of background constant gravitational field and investigated the problem of fall of the particle at the center. This problem is also treated classically. Unlike the conventional scenario, here the fall occurs at the surface of a sphere of unit radius.


title: Provably Optimal Algorithms for Generalized Linear Contextual Bandits

abstract: Contextual bandits are widely used in Internet services from news recommendation to advertising, and to Web search. Generalized linear models (logistical regression in particular) have demonstrated stronger performance than linear models in many applications where rewards are binary. However, most theoretical analyses on contextual bandits so far are on linear bandits. In this work, we propose an upper confidence bound based algorithm for generalized linear contextual bandits, which achieves an $\tilde{O}(\sqrt{dT})$ regret over $T$ rounds with $d$ dimensional feature vectors. This regret matches the minimax lower bound, up to logarithmic terms, and improves on the best previous result by a $\sqrt{d}$ factor, assuming the number of arms is fixed. A key component in our analysis is to establish a new, sharp finite-sample confidence bound for maximum-likelihood estimates in generalized linear models, which may be of independent interest. We also analyze a simpler upper confidence bound algorithm, which is useful in practice, and prove it to have optimal regret for certain cases.


title: Hadron Resonance Gas Model with Induced Surface Tension

abstract: Here we present a physically transparent generalization of the multicomponent Van der Waals equation of state in the grand canonical ensemble. For the one-component case the third and fourth virial coefficients are calculated analytically. It is shown that an adjustment of a single model parameter allows us to reproduce the third and fourth virial coefficients of the gas of hard spheres with small deviations from their exact values. A thorough comparison of the compressibility factor and speed of sound of the developed model with the one and two component Carnahan-Starling equation of state is made. It is shown that the model with the induced surface tension is able to reproduce the results of the Carnahan-Starling equation of state up to the packing fractions 0.2-0.22 at which the usual Van der Waals equation of state is inapplicable. At higher packing fractions the developed equation of state is softer than the gas of hard spheres and, hence, it breaks causality in the domain where the hadronic description is expected to be inapplicable. Using this equation of state we develop an entirely new hadron resonance gas model and apply it to a description of the hadron yield ratios measured at AGS, SPS, RHIC and ALICE energies of nuclear collisions. The achieved quality of the fit per degree of freedom is about 1.08. We confirm that the strangeness enhancement factor has a peak at low AGS energies, while at and above the highest SPS energy of collisions the chemical equilibrium of strangeness is observed. We argue that the chemical equilibrium of strangeness, i.e. $\gamma_s \simeq 1$, observed above the center of mass collision energy 4.3 GeV may be related to the hadronization of quark gluon bags which have the Hagedorn mass spectrum, and, hence, it may be a new signal for the onset of deconfinement.


title: SceneSeer: 3D Scene Design with Natural Language

abstract: Designing 3D scenes is currently a creative task that requires significant expertise and effort in using complex 3D design interfaces. This effortful design process starts in stark contrast to the easiness with which people can use language to describe real and imaginary environments. We present SceneSeer: an interactive text to 3D scene generation system that allows a user to design 3D scenes using natural language. A user provides input text from which we extract explicit constraints on the objects that should appear in the scene. Given these explicit constraints, the system then uses a spatial knowledge base learned from an existing database of 3D scenes and 3D object models to infer an arrangement of the objects forming a natural scene matching the input description. Using textual commands the user can then iteratively refine the created scene by adding, removing, replacing, and manipulating objects. We evaluate the quality of 3D scenes generated by SceneSeer in a perceptual evaluation experiment where we compare against manually designed scenes and simpler baselines for 3D scene generation. We demonstrate how the generated scenes can be iteratively refined through simple natural language commands.


title: Strong thermal and electrostatic manipulation of the Casimir force in graphene multilayers

abstract: We show that graphene-dielectric multilayers give rise to an unusual tunability of the Casimir-Lifshitz forces, and allow to easily realize completely different regimes within the same structure. Concerning thermal effects, graphene-dielectric multilayers take advantage from the anomalous features predicted for isolated suspended graphene sheets, even though they are considerably affected by the presence of the dielectric substrate. They can also archive the anomalous non-monotonic thermal metallic behavior by increasing the graphene sheets density and their Fermi energy. In addition to a strong thermal modulation occurring at short separations, in a region where the force is orders of magnitude larger than the one occurring at large distances, the force can be also adjusted by varying the number of graphene layers as well as their Fermi energy levels, allowing for relevant force amplifications which can be tuned, very rapidly and in-situ, by simply applying an electric potential. Our predictions can be relevant for both Casimir experiments and micro/nano electromechanical systems and in new devices for technological applications.


title: Sloan Digital Sky Survey IV: Mapping the Milky Way, Nearby Galaxies, and the Distant Universe

abstract: We describe the Sloan Digital Sky Survey IV (SDSS-IV), a project encompassing three major spectroscopic programs. The Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2) is observing hundreds of thousands of Milky Way stars at high resolution and high signal-to-noise ratio in the near-infrared. The Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey is obtaining spatially-resolved spectroscopy for thousands of nearby galaxies (median redshift of z = 0.03). The extended Baryon Oscillation Spectroscopic Survey (eBOSS) is mapping the galaxy, quasar, and neutral gas distributions between redshifts z = 0.6 and 3.5 to constrain cosmology using baryon acoustic oscillations, redshift space distortions, and the shape of the power spectrum. Within eBOSS, we are conducting two major subprograms: the SPectroscopic IDentification of eROSITA Sources (SPIDERS), investigating X-ray AGN and galaxies in X-ray clusters, and the Time Domain Spectroscopic Survey (TDSS), obtaining spectra of variable sources. All programs use the 2.5-meter Sloan Foundation Telescope at Apache Point Observatory; observations there began in Summer 2014. APOGEE-2 also operates a second near-infrared spectrograph at the 2.5-meter du Pont Telescope at Las Campanas Observatory, with observations beginning in early 2017. Observations at both facilities are scheduled to continue through 2020. In keeping with previous SDSS policy, SDSS-IV provides regularly scheduled public data releases; the first one, Data Release 13, was made available in July 2016.


title: Verified Low-Level Programming Embedded in F*

abstract: We present Low*, a language for low-level programming and verification, and its application to high-assurance optimized cryptographic libraries. Low* is a shallow embedding of a small, sequential, well-behaved subset of C in F*, a dependently-typed variant of ML aimed at program verification. Departing from ML, Low* does not involve any garbage collection or implicit heap allocation; instead, it has a structured memory model \`a la CompCert, and it provides the control required for writing efficient low-level security-critical code. By virtue of typing, any Low* program is memory safe. In addition, the programmer can make full use of the verification power of F* to write high-level specifications and verify the functional correctness of Low* code using a combination of SMT automation and sophisticated manual proofs. At extraction time, specifications and proofs are erased, and the remaining code enjoys a predictable translation to C. We prove that this translation preserves semantics and side-channel resistance. We provide a new compiler back-end from Low* to C and, to evaluate our approach, we implement and verify various cryptographic algorithms, constructions, and tools for a total of about 28,000 lines of code, specification and proof. We show that our Low* code delivers performance competitive with existing (unverified) C cryptographic libraries, suggesting our approach may be applicable to larger-scale low-level software.


title: Probing light sterile neutrino signatures at reactor and Spallation Neutron Source neutrino experiments

abstract: We investigate the impact of a fourth sterile neutrino at reactor and Spallation Neutron Source neutrino detectors. Specifically, we explore the discovery potential of the TEXONO and COHERENT experiments to subleading sterile neutrino effects through the measurement of the coherent elastic neutrino-nucleus scattering event rate. Our dedicated $\chi^2$-sensitivity analysis employs realistic nuclear structure calculations adequate for high purity sub-keV threshold Germanium detectors.


title: A Monadic Framework for Relational Verification: Applied to Information Security, Program Equivalence, and Optimizations

abstract: Relational properties describe multiple runs of one or more programs. They characterize many useful notions of security, program refinement, and equivalence for programs with diverse computational effects, and they have received much attention in the recent literature. Rather than developing separate tools for special classes of effects and relational properties, we advocate using a general purpose proof assistant as a unifying framework for the relational verification of effectful programs. The essence of our approach is to model effectful computations using monads and to prove relational properties on their monadic representations, making the most of existing support for reasoning about pure programs. We apply this method in F* and evaluate it by encoding a variety of relational program analyses, including information flow control, program equivalence and refinement at higher order, correctness of program optimizations and game-based cryptographic security. By relying on SMT-based automation, unary weakest preconditions, user-defined effects, and monadic reification, we show that, compared to unary properties, verifying relational properties requires little additional effort from the F* programmer.


title: Fair prediction with disparate impact: A study of bias in recidivism prediction instruments

abstract: Recidivism prediction instruments (RPI's) provide decision makers with an assessment of the likelihood that a criminal defendant will reoffend at a future point in time. While such instruments are gaining increasing popularity across the country, their use is attracting tremendous controversy. Much of the controversy concerns potential discriminatory bias in the risk assessments that are produced. This paper discusses several fairness criteria that have recently been applied to assess the fairness of recidivism prediction instruments. We demonstrate that the criteria cannot all be simultaneously satisfied when recidivism prevalence differs across groups. We then show how disparate impact can arise when a recidivism prediction instrument fails to satisfy the criterion of error rate balance.


title: Rook placements and Jordan forms of upper-triangular nilpotent matrices

abstract: The set of n by n upper-triangular nilpotent matrices with entries in a finite field F_q has Jordan canonical forms indexed by partitions lambda of n. We present a combinatorial formula for computing the number F_\lambda(q) of matrices of Jordan type lambda as a weighted sum over standard Young tableaux. We also study a connection between these matrices and non-attacking rook placements, which leads to a refinement of the formula for F_\lambda(q).


title: On testing the simulation theory

abstract: Can the theory that reality is a simulation be tested? We investigate this question based on the assumption that if the system performing the simulation is finite (i.e. has limited resources), then to achieve low computational complexity, such a system would, as in a video game, render content (reality) only at the moment that information becomes available for observation by a player and not at the moment of detection by a machine (that would be part of the simulation and whose detection would also be part of the internal computation performed by the Virtual Reality server before rendering content to the player). Guided by this principle we describe conceptual wave/particle duality experiments aimed at testing the simulation theory.


title: Automorphisms of Products of Drinfeld Half Planes

abstract: The Drinfeld upper half-planes play the role of symmetric spaces in the $p$-adic analytic world. We find the automorphism group of a product of such spaces, where each may be defined over a different field. We deduce a rigidity theorem for quotients of such products by discrete and torsion free groups.


title: Achieving non-discrimination in prediction

abstract: Discrimination-aware classification is receiving an increasing attention in data science fields. The pre-process methods for constructing a discrimination-free classifier first remove discrimination from the training data, and then learn the classifier from the cleaned data. However, they lack a theoretical guarantee for the potential discrimination when the classifier is deployed for prediction. In this paper, we fill this gap by mathematically bounding the probability of the discrimination in prediction being within a given interval in terms of the training data and classifier. We adopt the causal model for modeling the data generation mechanism, and formally defining discrimination in population, in a dataset, and in prediction. We obtain two important theoretical results: (1) the discrimination in prediction can still exist even if the discrimination in the training data is completely removed; and (2) not all pre-process methods can ensure non-discrimination in prediction even though they can achieve non-discrimination in the modified training data. Based on the results, we develop a two-phase framework for constructing a discrimination-free classifier with a theoretical guarantee. The experiments demonstrate the theoretical results and show the effectiveness of our two-phase framework.


title: SceneSuggest: Context-driven 3D Scene Design

abstract: We present SceneSuggest: an interactive 3D scene design system providing context-driven suggestions for 3D model retrieval and placement. Using a point-and-click metaphor we specify regions in a scene in which to automatically place and orient relevant 3D models. Candidate models are ranked using a set of static support, position, and orientation priors learned from 3D scenes. We show that our suggestions enable rapid assembly of indoor scenes. We perform a user study comparing suggestions to manual search and selection, as well as to suggestions with no automatic orientation. We find that suggestions reduce total modeling time by 32%, that orientation priors reduce time spent re-orienting objects by 27%, and that context-driven suggestions reduce the number of text queries by 50%.


title: Optimal Investment and Pricing in the Presence of Defaults

abstract: We consider the optimal investment problem when the traded asset may default, causing a jump in its price. For an investor with constant absolute risk aversion, we compute indifference prices for defaultable bonds, as well as a price for dynamic protection against default. For the latter problem, our work complements Sircar & Zariphopoulou (2007), where it is implicitly assumed the investor is protected against default. We consider a factor model where the asset's instantaneous return, variance, correlation and default intensity are driven by a time-homogenous diffusion X taking values in an arbitrary region E. We identify the certainty equivalent with a semi-linear degenerate elliptic partial differential equation with quadratic growth in both function and gradient. Under a minimal integrability assumption on the market price of risk, we show the certainty equivalent is a classical solution. In particular, our results cover when X is a one-dimensional affine diffusion and when returns, variances and default intensities are also affine. Numerical examples highlight the relationship between the factor process and both the indifference price and default insurance. Lastly, we show the insurance protection price is not the default intensity under the dual optimal measure.


title: Quantum multi-parameter estimation with generalized balanced multi-mode NOON-like states

abstract: The simultaneous multi-parameter estimation problem using a class of multi-mode entangled states is investigated in this paper. Specifically, the problem of optical phase imaging is considered and the quantum probe is taken to be a balanced coherent superposition of components with an arbitrary quantum state in one mode and vacuum states in all the other modes, which is a generalization of the multi-mode NOON state. The analytical form for the quantum Cramer-Rao bound (QCRB) is presented, which shows the performance by providing a lower bound of the estimation uncertainty. It is shown that the NOON state has the worst performance among those in the class of the entangled states considered. We also analyze in detail four different scenarios, which are the NOON state, the entangled coherent state, the entangled squeezed coherent state, and the entangled squeezed vacuum state. From the comparison among these four states, we find that when the mean photon number is fixed, the squeezed vacuum state has the smallest QCRB, followed by the squeezed coherent state, entangled coherent state, and NOON state. We also illustrate that the balanced entangled state can perform better than a more generalized unbalanced form studied in previous works for certain scenarios. Finally, we give an experimental setup for producing a two-mode entangled state that can beat the NOON state in quantum metrology.


title: Ending the Anomaly: Achieving Low Latency and Airtime Fairness in WiFi

abstract: With more devices connected, delays and jitter at the WiFi hop become more prevalent, and correct functioning during network congestion becomes more important. However, two important performance issues prevent modern WiFi from reaching its potential: Increased latency under load caused by excessive queueing (i.e. bufferbloat) and the 802.11 performance anomaly. To remedy these issues, we present a novel two-part solution: We design a new queueing scheme that eliminates bufferbloat in the wireless setting. Leveraging this queueing scheme, we then design an airtime fairness scheduler that operates at the access point and doesn't require any changes to clients. We evaluate our solution using both a theoretical model and experiments in a testbed environment, formulating a suitable analytical model in the process. We show that our solution achieves an order of magnitude reduction in latency under load, large improvements in multi-station throughput, and nearly perfect airtime fairness for both TCP and downstream UDP traffic. Further experiments with application traffic confirm that the solution provides significant performance gains for real-world traffic.We develop a production quality implementation of our solution in the Linux kernel, the platform powering most access points outside of the managed enterprise setting. The implementation has been accepted into the mainline kernel distribution, making it available for deployment on billions of devices running Linux today.


title: Groups with one or two super-Brauer character theories

abstract: A super-Brauer character theory of a group $G$ and a prime $p$ is a pair consisting of a partition of the irreducible $p$-Brauer characters and a partition of the $p$-regular elements of $G$ that satisfy certain properties. We classify the groups and primes that have exactly one super-Brauer character theory. We discuss the groups with exactly two super-Brauer character theories.


title: On the Power of Learning from $k$-Wise Queries

abstract: Several well-studied models of access to data samples, including statistical queries, local differential privacy and low-communication algorithms rely on queries that provide information about a function of a single sample. (For example, a statistical query (SQ) gives an estimate of $Ex_{x \sim D}[q(x)]$ for any choice of the query function $q$ mapping $X$ to the reals, where $D$ is an unknown data distribution over $X$.) Yet some data analysis algorithms rely on properties of functions that depend on multiple samples. Such algorithms would be naturally implemented using $k$-wise queries each of which is specified by a function $q$ mapping $X^k$ to the reals. Hence it is natural to ask whether algorithms using $k$-wise queries can solve learning problems more efficiently and by how much. Blum, Kalai and Wasserman (2003) showed that for any weak PAC learning problem over a fixed distribution, the complexity of learning with $k$-wise SQs is smaller than the (unary) SQ complexity by a factor of at most $2^k$. We show that for more general problems over distributions the picture is substantially richer. For every $k$, the complexity of distribution-independent PAC learning with $k$-wise queries can be exponentially larger than learning with $(k+1)$-wise queries. We then give two approaches for simulating a $k$-wise query using unary queries. The first approach exploits the structure of the problem that needs to be solved. It generalizes and strengthens (exponentially) the results of Blum et al.. It allows us to derive strong lower bounds for learning DNF formulas and stochastic constraint satisfaction problems that hold against algorithms using $k$-wise queries. The second approach exploits the $k$-party communication complexity of the $k$-wise query function.


title: Conductance oscillations in graphene/nanoclusters hybrid material: towards large area single electron devices

abstract: Large assemblies of self-organized aluminum nanoclusters embedded in an oxide layer are formed on graphene templates and used to build tunnel-junction devices. Unexpectedly, single-electron-transport behavior with well-defined Coulomb oscillations is observed for a record junction area containing millions of metal islands. Such hybrid materials offer new prospects for single-electron electronics.


title: Magnetars

abstract: Magnetars are young and highly magnetized neutron stars which display a wide array of X-ray activity including short bursts, large outbursts, giant flares and quasi-periodic oscillations, often coupled with interesting timing behavior including enhanced spin-down, glitches and anti-glitches. The bulk of this activity is explained by the evolution and decay of an ultrastrong magnetic field, stressing and breaking the neutron star crust, which in turn drives twists of the external magnetosphere and powerful magnetospheric currents. The population of detected magnetars has grown to about 30 objects and shows unambiguous phenomenological connection with very highly magnetized radio pulsars. Recent progress in magnetar theory includes explanation of the hard X-ray component in the magnetar spectrum and development of surface heating models, explaining the sources' remarkable radiative output.


title: Deep Image Harmonization

abstract: Compositing is one of the most common operations in photo editing. To generate realistic composites, the appearances of foreground and background need to be adjusted to make them compatible. Previous approaches to harmonize composites have focused on learning statistical relationships between hand-crafted appearance features of the foreground and background, which is unreliable especially when the contents in the two layers are vastly different. In this work, we propose an end-to-end deep convolutional neural network for image harmonization, which can capture both the context and semantic information of the composite images during harmonization. We also introduce an efficient way to collect large-scale and high-quality training data that can facilitate the training process. Experiments on the synthesized dataset and real composite images show that the proposed network outperforms previous state-of-the-art methods.


title: Analytical calculation of the axis angle 2V from extinction measurements on the spindle stage

abstract: A concise derivation of the "Joel equations", which allow for the determination of the axis angle 2V from measurements of extinction directions on a spindle stage, is provided starting from the wave-equation. Only analytic methods and no geometric arguments referring to stereographic projections are invoked. For error free data, the resulting equations allow for a closed form solution. If angle data with measurement error are to be used, a maximum likelihood estimation methodology is proposed, which may be solved using e.g. iterative reweighting. The method is tested with and compared to published data.


title: The response of relativistic outflowing gas to the inner accretion disk of a black hole

abstract: Active galactic nucleus (AGN) feedback is the process by which supermassive black holes in the centres of galaxies may moderate the growth of their hosts. Gas outflows from supermassive black holes release huge quantities of energy into the interstellar medium, clearing the surrounding gas. The most extreme of these, the ultra-fast outflows (UFOs), are the subset of X-ray detected outflows with velocities higher than 10,000 km/s, believed to originate in relativistic disc winds, a few hundred gravitational radii from the black hole. The absorption features produced by these outflows are variable, but no clear link has been found between the behaviour of the X-ray continuum and the energy or equivalent width of the outflow features due to the long time-scales of quasar variability. Here, we present the detection of multiple absorption lines from an extreme ultra-fast gas flow in the X-ray spectrum of the active galactic nucleus IRAS 13224-3809, at 0.236+/-0.006 times the speed of light (71,000 km/s), where the absorption is strongly anti-correlated with the emission from the inner regions of the accretion disk. If the gas flow is identified as a genuine outflow then it is in the fastest 5 per cent of such winds, and its variability is hundreds of times faster than in other variable winds, allowing us to observe in hours what would take months in a quasar. We find signatures of the wind simultaneously in both low and high energy detectors, which are consistent with a single ionized outflow, linking the two phenomena. The detection of the wind responding to the emission from the inner disk demonstrates a connection between accretion processes occurring on very different scales, with the X-rays from within a few gravitational radii of the black hole ionizing the relativistically outflowing gas as the flux rises.


title: Experimental Synthetic Aperture Radar with Dynamic Metasurfaces

abstract: We investigate the use of a dynamic metasurface as the transmitting antenna for a synthetic aperture radar (SAR) imaging system. The dynamic metasurface consists of a one-dimensional microstrip waveguide with complementary electric resonator (cELC) elements patterned into the upper conductor. Integrated into each of the cELCs are two diodes that can be used to shift each cELC resonance out of band with an applied voltage. The aperture is designed to operate at K band frequencies (17.5 to 20.3 GHz), with a bandwidth of 2.8 GHz. We experimentally demonstrate imaging with a fabricated metasurface aperture using existing SAR modalities, showing image quality comparable to traditional antennas. The agility of this aperture allows it to operate in spotlight and stripmap SAR modes, as well as in a third modality inspired by computational imaging strategies. We describe its operation in detail, demonstrate high-quality imaging in both 2D and 3D, and examine various trade-offs governing the integration of dynamic metasurfaces in future SAR imaging platforms.


title: Physically unclonable function using initial waveform of ring oscillators on 65 nm CMOS technology

abstract: A silicon physically unclonable function (PUF) using ring oscillators (ROs) has the advantage of easy application in both an application specific integrated circuit (ASIC) and a field-programmable gate array (FPGA). Here, we provide a RO-PUF using the initial waveform of the ROs based on 65 nm CMOS technology. Compared with the conventional RO-PUF, the number of ROs is greatly reduced and the time needed to generate an ID is within a couple of system clocks.


title: Linking Magnon-Cavity Strong Coupling to Magnon-Polaritons through Effective Permeability

abstract: Strong coupling in cavity-magnon systems has shown great potential for use in spintronics and information processing technologies due to the low damping rates and long coherence times. Although such systems are conceptually similar to those coupled by magnon-polaritons (MPs), the link between magnon-cavity coupling and MPs has not been explicitly defined. In this work we establish such a connection by studying the frequency-wavevector dispersion of a strongly coupled magnon-cavity system, using a height-adjustable microwave cavity, and by modelling the observed behaviour through the system's effective permeability. A polariton gap between the upper and lower coupled modes of the magnon-cavity system is defined, and is seen to be dependent on the system's effective filling factor. This gap is equal to the MP polariton gap in the limit where filling factor = 1, corresponding to the removal of the microwave cavity. Thus, our work clarifies the connection between magnon-cavity and MP coupling, improving our understanding of magnon-photon interactions in coupled systems.


title: Discrete Wavelet Transform Based Algorithm for Recognition of QRS Complexes

abstract: This paper proposes the application of Discrete Wavelet Transform (DWT) to detect the QRS (ECG is characterized by a recurrent wave sequence of P, QRS and T-wave) of an electrocardiogram (ECG) signal. Wavelet Transform provides localization in both time and frequency. In preprocessing stage, DWT is used to remove the baseline wander in the ECG signal. The performance of the algorithm of QRS detection is evaluated against the standard MIT BIH (Massachusetts Institute of Technology, Beth Israel Hospital) Arrhythmia database. The average QRS complexes detection rate of 98.1 % is achieved.


title: A Global Compact Result for a Fractional Elliptic Problem with Critical Sobolev-Hardy Nonlinearities on ${\mathbb R}^N$

abstract: In this paper, we are concerned with the following type of elliptic problems: $$ (-\Delta)^{\alpha} u+a(x) u=\frac{|u|^{2^*_{s}-2}u}{|x|^s}+k(x)|u|^{q-2}u, u\,\in\,H^\alpha({\mathbb R}^N), $$ where $2<q< 2^*$, $0<\alpha<1$, $0<s<2\alpha$, $2^*_{s}=2(N-s)/(N-2\alpha)$ is the critical Sobolev-Hardy exponent, $2^*=2N/(N-2\alpha)$ is the critical Sobolev exponent, $a(x),k(x)\in C({\mathbb R}^N)$. Through a compactness analysis of the functional associated to the problem, we obtain the existence of positive solutions under certain assumptions on $a(x),k(x)$.


title: Branching fraction measurements of $\psi(3686) \to \gamma \chi_{cJ}$

abstract: Using a sample of 106 million $\psi(3686)$ decays, the branching fractions of $\psi(3686) \to \gamma \chi_{c0}, \psi(3686) \to \gamma \chi_{c1}$, and $\psi(3686) \to \gamma \chi_{c2}$ are measured with improved precision to be $(9.389 \pm 0.014 \pm 0.332)\,\%$, $(9.905 \pm 0.011 \pm 0.353)\,\% $, and $(9.621 \pm 0.013 \pm 0.272)\,\% $, respectively, where the first uncertainties are statistical and the second ones are systematic. The product branching fractions of $\psi(3686) \to \gamma \chi_{c1}, \chi_{c1} \to \gamma J/\psi$ and $\psi(3686) \to \gamma \chi_{c2}, \chi_{c2} \to \gamma J/\psi$ and the branching fractions of $\chi_{c1} \to \gamma J/\psi$ and $\chi_{c2} \to \gamma J/\psi$ are also presented.


title: Simple method for determining binding energies of fullerene negative ions

abstract: A robust potential wherein is embedded the crucial core polarization interaction is used in the Regge Pole methodology to calculate low energy electron elastic scattering total cross section (TCS) for the C60 fullerene in the electron impact energy range 0.02 through 10.0 eV. The energy position of the characteristic dramatically sharp resonance appearing at the second Ramsauer Townsend (RT) minimum of the TCS representing stable C60 fullerene negative ion formation agrees excellently with the measured electron affinity (EA) of C60 [Huang et al 2014 J. Chem. Phys. 140 224315]. The benchmarked potential and the Regge-pole method are then used to calculate electron elastic scattering TCSs for selected fullerenes, from C54 through C240. The TCSs are found to be characterized generally by RT minima, shape resonances (SRs) and dramatically sharp resonances representing long lived ground state fullerene negative ion formation. For the TCSs of C70, C76, C78, and C84 the agreement between the energy positions of the very sharp resonances, corresponding to the binding energies (BEs) of the resultant fullerene negative ions, and the measured EAs is outstanding. Additionally, we extract the BEs of the resultant fullerene negative ions from our calculated TCSs of the C86, C90 and C92 fullerenes with estimated EAs larger than 3.0 eV by the experiment [Boltalina et al, 1993 Rapid Commun. Mass Spectrom. 7 1009] as well as of other fullerenes, including C180 and C240. Most of the TCSs presented in this paper are the first and only. Our novel approach is general and should be applicable to other fullerenes as well and complex heavy atoms, such as the lanthanide atoms. We conclude with a remark on the catalytic properties of the fullerenes through their negative ions.


title: Fault Tolerant Thermal Control of Steam Turbine Shell Deflections

abstract: The metal-to-metal clearances of a steam turbine during full or part load operation are among the main drivers of efficiency. The requirement to add clearances is driven by a number of factors including the relative movements of the steam turbine shell and rotor during transient conditions such as startup and shutdown. This paper includes a description of a control algorithm to manage external heating blankets for the thermal control of the shell deflections during turbine shutdown. The proposed method is tolerant of changes in the heat loss characteristics of the system as well as simultaneous component failures.


title: Efficient Computation of Subspace Skyline over Categorical Domains

abstract: Platforms such as AirBnB, Zillow, Yelp, and related sites have transformed the way we search for accommodation, restaurants, etc. The underlying datasets in such applications have numerous attributes that are mostly Boolean or Categorical. Discovering the skyline of such datasets over a subset of attributes would identify entries that stand out while enabling numerous applications. There are only a few algorithms designed to compute the skyline over categorical attributes, yet are applicable only when the number of attributes is small. In this paper, we place the problem of skyline discovery over categorical attributes into perspective and design efficient algorithms for two cases. (i) In the absence of indices, we propose two algorithms, ST-S and ST-P, that exploits the categorical characteristics of the datasets, organizing tuples in a tree data structure, supporting efficient dominance tests over the candidate set. (ii) We then consider the existence of widely used precomputed sorted lists. After discussing several approaches, and studying their limitations, we propose TA-SKY, a novel threshold style algorithm that utilizes sorted lists. Moreover, we further optimize TA-SKY and explore its progressive nature, making it suitable for applications with strict interactive requirements. In addition to the extensive theoretical analysis of the proposed algorithms, we conduct a comprehensive experimental evaluation of the combination of real (including the entire AirBnB data collection) and synthetic datasets to study the practicality of the proposed algorithms. The results showcase the superior performance of our techniques, outperforming applicable approaches by orders of magnitude.


title: Construction of a blow-up solution for the Complex Ginzburg-Landau equation in some critical case

abstract: We construct a solution for the Complex Ginzburg-Landau equation in some critical case, which blows up in finite time $T$ only at one blow-up point. We also give a sharp description of its profile. The proof relies on the reduction of the problem to a finite dimensional one, and the use of index theory to conclude. The interpretation of the parameters of the finite dimension problem in terms of the blow-up point and time allows to prove the stability of the constructed solution.


title: Minimal models of quantum homotopy Lie algebras via the BV-formalism

abstract: Using the BV-formalism of mathematical physics an explicit construction for the minimal model of a quantum L-infinity-algebra is given as a formal super integral. The approach taken herein to these formal integrals is axiomatic; they can be approached using perturbation theory to obtain combinatorial formulae as shown in the appendix. Additionally, there exists a canonical differential graded Lie algebra morphism mapping formal functions on homology to formal functions on the whole space. An L-infinity-algebra morphism inverse to this differential graded Lie algebra morphism on the level of homology is constructed as a formal super integral.


title: Distributed Frequency Control with Operational Constraints, Part II: Network Power Balance

abstract: In Part I of this paper we propose a decentralized optimal frequency control of multi-area power system with operational constraints, where the tie-line powers remain unchanged in the steady state and the power mismatch is balanced within individual control areas. In Part II of the paper, we propose a distributed controller for optimal frequency control in the network power balance case, where the power mismatch is balanced over the whole system. With the proposed controller, the tie-line powers remain within the acceptable range at equilibrium, while the regulation capacity constraints are satisfied both at equilibrium and during transient. It is revealed that the closed-loop system with the proposed controller carries out primal-dual updates with saturation for solving an associated optimization problem. To cope with discontinuous dynamics of the closed-loop system, we deploy the invariance principle for nonpathological Lyapunov function to prove its asymptotic stability. Simulation results are provided to show the effectiveness of our controller.


title: Multi-Sensor Data Pattern Recognition for Multi-Target Localization: A Machine Learning Approach

abstract: Data-target pairing is an important step towards multi-target localization for the intelligent operation of unmanned systems. Target localization plays a crucial role in numerous applications, such as search, and rescue missions, traffic management and surveillance. The objective of this paper is to present an innovative target location learning approach, where numerous machine learning approaches, including K-means clustering and supported vector machines (SVM), are used to learn the data pattern across a list of spatially distributed sensors. To enable the accurate data association from different sensors for accurate target localization, appropriate data pre-processing is essential, which is then followed by the application of different machine learning algorithms to appropriately group data from different sensors for the accurate localization of multiple targets. Through simulation examples, the performance of these machine learning algorithms is quantified and compared.


title: Pulsating dynamics of thermal plumes and its implications for multiple eruption events in the Deccan Traps, India

abstract: In Earth's mantle gravity instabilities initiated by density inversion lead to upwelling of hot materials as plumes. This study focuses upon the problem of their ascent dynamics to provide an explanation of the periodic multiple eruption events in large igneous provinces and hotspots. We demonstrate from physical experiments that plumes can ascend in a continuous process to form a single large head trailing into a long slender tail, typically described in the literature only under specific physical conditions. Alternatively, they ascend in a pulsating fashion, and produce multiple in-axis heads of varying dimensions. Based on the Volume of Fluid (VOF) method, we performed computational fluid dynamics (CFD) simulations to constrain the thermo-mechanical conditions that decide the continuous versus pulsating dynamics. Our CFD simulations suggest the density ($\rho^*$) and the viscosity ($R$) ratios of the ambient to the plume materials and the influx rates ($Re$) are the prime factors in controlling the ascent dynamics. Conditions with large $R$ ($>$50) develop pulsating plumes, which are sustained preferentially under low $\rho^*$ and $Re$ conditions. Again, the increasing temperature difference between the plume and the ambient medium is found to promote pulsating behaviour. From these CFD simulations we show the thermal structures of a plume, and predict the peak thermal events near the surface that commence periodically as pulses with a time interval of 1.4-3 Ma. The pulsating plume model explains the multiple eruption events in the Deccan Traps in India, where the time scale of their intervals estimated from the available $^{40}Ar$/$^{39}Ar$ geochronological data shows an excellent match with our simulation results.


title: A comparison of observed and simulated absorption from HI, CIV, and SiIV around $z\approx2$ star-forming galaxies suggests redshift-space distortions are due to inflows

abstract: We study HI and metal-line absorption around $z\approx2$ star-forming galaxies by comparing an analysis of data from the Keck Baryonic Structure Survey to mock spectra generated from the EAGLE cosmological, hydrodynamical simulations. We extract sightlines from the simulations and compare the properties of the absorption by HI, CIV and SiIV around simulated and observed galaxies using pixel optical depths. We mimic the resolution, pixel size, and signal-to-noise ratio of the observations, as well as the distributions of impact parameters and galaxy redshift errors. We find that the EAGLE reference model is in excellent agreement with the observations. In particular, the simulation reproduces the high metal-line optical depths found at small galactocentric distances, the optical depth enhancements out to impact parameters of 2 proper Mpc, and the prominent redshift-space distortions which we find are due to peculiar velocities rather than redshift errors. The agreement is best for halo masses $\sim10^{12.0}$ M$_\odot$, for which the observed and simulated stellar masses also agree most closely. We examine the median ion mass-weighted radial gas velocities around the galaxies, and find that most of the gas is infalling, with the infall velocity depending on halo rather than stellar mass. From this we conclude that the observed redshift-space distortions are predominantly caused by infall rather than outflows.


title: Supervised Saliency Map Driven Segmentation of the Lesions in Dermoscopic Images

abstract: Lesion segmentation is the first step in most automatic melanoma recognition systems. Deficiencies and difficulties in dermoscopic images such as color inconstancy, hair occlusion, dark corners and color charts make lesion segmentation an intricate task. In order to detect the lesion in the presence of these problems, we propose a supervised saliency detection method tailored for dermoscopic images based on the discriminative regional feature integration (DRFI). DRFI method incorporates multi-level segmentation, regional contrast, property, background descriptors, and a random forest regressor to create saliency scores for each region in the image. In our improved saliency detection method, mDRFI, we have added some new features to regional property descriptors. Also, in order to achieve more robust regional background descriptors, a thresholding algorithm is proposed to obtain a new pseudo-background region. Findings reveal that mDRFI is superior to DRFI in detecting the lesion as the salient object in dermoscopic images. The proposed overall lesion segmentation framework uses detected saliency map to construct an initial mask of the lesion through thresholding and post-processing operations. The initial mask is then evolving in a level set framework to fit better on the lesion's boundaries. The results of evaluation tests on three public datasets show that our proposed segmentation method outperforms the other conventional state-of-the-art segmentation algorithms and its performance is comparable with most recent approaches that are based on deep convolutional neural networks.


title: Combinatorial models for Schubert polynomials

abstract: Schubert polynomials are a basis for the polynomial ring that represent Schubert classes for the flag manifold. In this paper, we introduce and develop several new combinatorial models for Schubert polynomials that relate them to other known bases including key polynomials and fundamental slide polynomials. We unify these and existing models by giving simple bijections between the combinatorial objects indexing each. In particular, we give a simple bijective proof that the balanced tableaux of Edelman and Greene enumerate reduced expressions and a direct combinatorial proof of Kohnert's algorithm for computing Schubert polynomials. Further, we generalize the insertion algorithm of Edelman and Greene to give a bijection between reduced expressions and pairs of tableaux of the same key diagram shape and use this to give a simple formula, directly in terms of reduced expressions, for the key polynomial expansion of a Schubert polynomial.


title: A Joint Identification Approach for Argumentative Writing Revisions

abstract: Prior work on revision identification typically uses a pipeline method: revision extraction is first conducted to identify the locations of revisions and revision classification is then conducted on the identified revisions. Such a setting propagates the errors of the revision extraction step to the revision classification step. This paper proposes an approach that identifies the revision location and the revision type jointly to solve the issue of error propagation. It utilizes a sequence representation of revisions and conducts sequence labeling for revision identification. A mutation-based approach is utilized to update identification sequences. Results demonstrate that our proposed approach yields better performance on both revision location extraction and revision type classification compared to a pipeline baseline.


title: Lagrangian Mean Curvature Flows and Moment maps

abstract: In this paper, we construct various examples of Lagrangian mean curvature flows in Calabi-Yau manifolds, using moment maps for actions of abelian Lie groups on them. The examples include Lagrangian self-shrinkers and translating solitons in the Euclidean spaces. Moreover, our method can be applied to construct examples of Lagrangian mean curvature flows in non-flat Calabi-Yau manifolds. In particular, we describe Lagrangian mean curvature flows in 4-dimensional Ricci-flat ALE spaces in detail and investigate their singularities.


title: Semi-analytical approximations to statistical moments of sigmoid and softmax mappings of normal variables

abstract: This note is concerned with accurate and computationally efficient approximations of moments of Gaussian random variables passed through sigmoid or softmax mappings. These approximations are semi-analytical (i.e. they involve the numerical adjustment of parametric forms) and highly accurate (they yield 5% error at most). We also highlight a few niche applications of these approximations, which arise in the context of, e.g., drift-diffusion models of decision making or non-parametric data clustering approaches. We provide these as examples of efficient alternatives to more tedious derivations that would be needed if one was to approach the underlying mathematical issues in a more formal way. We hope that this technical note will be helpful to modellers facing similar mathematical issues, although maybe stemming from different academic prospects.


title: Lectures on Structural Stability in Dynamics

abstract: These lectures present results and problems on the characterization of structurally stable dynamics. We will shed light those which do not seem to depend on the regularity class (holomorphic or differentiable). Furthermore, we will present some links between the problems of structural stability in dynamical systems and in singularity theory.


title: Accurate gradient computations at interfaces using finite element methods

abstract: New finite element methods are proposed for elliptic interface problems in one and two dimensions. The main motivation is not only to get an accurate solution but also an accurate first order derivative at the interface (from each side). The key in 1D is to use the idea from \cite{wheeler1974galerkin}. For 2D interface problems, the idea is to introduce a small tube near the interface and introduce the gradient as part of unknowns, which is similar to a mixed finite element method, except only at the interface. Thus the computational cost is just slightly higher than the standard finite element method. We present rigorous one dimensional analysis, which show second order convergence order for both of the solution and the gradient in 1D. For two dimensional problems, we present numerical results and observe second order convergence for the solution, and super-convergence for the gradient at the interface.


title: Extreme points and saturated polynomials

abstract: We consider the problem of characterizing the extreme points of the set of analytic functions f on the bidisk with positive real part and f(0)=1. If one restricts to those f whose Cayley transform is a rational inner function, one gets a more tractable problem. We construct families of such f that are extreme points and conjecture that these are all such extreme points. These extreme points are constructed from polynomials dubbed saturated, which roughly speaking means they have no zeros in the bidisk and as many zeros as possible on the boundary without having infinitely many zeros.


title: Active End-Effector Pose Selection for Tactile Object Recognition through Monte Carlo Tree Search

abstract: This paper considers the problem of active object recognition using touch only. The focus is on adaptively selecting a sequence of wrist poses that achieves accurate recognition by enclosure grasps. It seeks to minimize the number of touches and maximize recognition confidence. The actions are formulated as wrist poses relative to each other, making the algorithm independent of absolute workspace coordinates. The optimal sequence is approximated by Monte Carlo tree search. We demonstrate results in a physics engine and on a real robot. In the physics engine, most object instances were recognized in at most 16 grasps. On a real robot, our method recognized objects in 2--9 grasps and outperformed a greedy baseline.


title: Gram-CTC: Automatic Unit Selection and Target Decomposition for Sequence Labelling

abstract: Most existing sequence labelling models rely on a fixed decomposition of a target sequence into a sequence of basic units. These methods suffer from two major drawbacks: 1) the set of basic units is fixed, such as the set of words, characters or phonemes in speech recognition, and 2) the decomposition of target sequences is fixed. These drawbacks usually result in sub-optimal performance of modeling sequences. In this pa- per, we extend the popular CTC loss criterion to alleviate these limitations, and propose a new loss function called Gram-CTC. While preserving the advantages of CTC, Gram-CTC automatically learns the best set of basic units (grams), as well as the most suitable decomposition of tar- get sequences. Unlike CTC, Gram-CTC allows the model to output variable number of characters at each time step, which enables the model to capture longer term dependency and improves the computational efficiency. We demonstrate that the proposed Gram-CTC improves CTC in terms of both performance and efficiency on the large vocabulary speech recognition task at multiple scales of data, and that with Gram-CTC we can outperform the state-of-the-art on a standard speech benchmark.


title: Stability of Stationary Inverse Transport Equation in Diffusion Scaling

abstract: We consider the inverse problem of reconstructing the optical parameters for stationary radiative transfer equation (RTE) from velocity-averaged measurement. The RTE often contains multiple scales characterized by the magnitude of a dimensionless parameter---the Knudsen number ($K_n$). In the diffusive scaling ($K_n \ll 1$), the stationary RTE is well approximated by an elliptic equation in the forward setting. However, the inverse problem for the elliptic equation is acknowledged to be severely ill-posed as compared to the well-posedness of inverse transport equation, which raises the question of how uniqueness being lost as $K_n \rightarrow 0$. We tackle this problem by examining the stability of inverse problem with varying $K_n$. We show that, the discrepancy in two measurements is amplified in the reconstructed parameters at the order of $K_n^p~ (p = 1\text{ or} ~2)$, and as a result lead to ill-posedness in the zero limit of $K_n$. Our results apply to both continuous and discrete settings. Some numerical tests are performed in the end to validate these theoretical findings.