title: High Dimensional Data Modeling Techniques for Detection of Chemical Plumes and Anomalies in Hyperspectral Images and Movies

abstract: We briefly review recent progress in techniques for modeling and analyzing hyperspectral images and movies, in particular for detecting plumes of both known and unknown chemicals. For detecting chemicals of known spectrum, we extend the technique of using a single subspace for modeling the background to a "mixture of subspaces" model to tackle more complicated background. Furthermore, we use partial least squares regression on a resampled training set to boost performance. For the detection of unknown chemicals we view the problem as an anomaly detection problem, and use novel estimators with low-sampled complexity for intrinsically low-dimensional data in high-dimensions that enable us to model the "normal" spectra and detect anomalies. We apply these algorithms to benchmark data sets made available by the Automated Target Detection program co-funded by NSF, DTRA and NGA, and compare, when applicable, to current state-of-the-art algorithms, with favorable results.


title: Void alignment and density profile applied to measuring cosmological parameters

abstract: We study the orientation and density profiles of the cosmological voids with SDSS10 data. Using voids to test Alcock-Paczynski effect has been proposed and tested in both simulations and actual SDSS data. Previous observations imply that there exist an empirical stretching factor which plays an important role in the voids' orientation. Simulations indicate that this empirical stretching factor is caused by the void galaxies' peculiar velocities. Recently Hamaus et al. found that voids' density profiles are universal and their average velocities satisfy linear theory very well. In this article we first confirm that the stretching effect exists using independent analysis. We then apply the universal density profile to measure the cosmological parameters. We find that the void density profile can be a tool to measure the cosmological parameters.


title: Ab initio investigation of the adsorption of zoledronic acid molecule on hydroxyapatite (001) surface: an atomistic insight of bone protection

abstract: We report a computational study of the adsorption of zoledronic acid molecule on hydroxyapatite (001) surface within ab initio density functional theory. The systematic study has been performed, from hydroxyapatite bulk and surface, and zoledronic acid molecule to the adsorption of the molecule on the surface. The optimized bond lengths and bond angles were obtained and analyzed, giving an evidence of structural similarity between subjects under study. The formation energies of hydroxyapatite (001) surfaces with two kinds of terminations were computed as about 1.2 and 1.5 J/m^2 with detailed atomistic structural information. We determined the adsorption energies of zoledronic acid molecule on the surfaces, which are -260 kJ/mol at 0.25 ML and -400 kJ/mol at 0.5 ML. An atomistic insight of strong binding affinity of zoledronic acid to the hydroxyapatite surface was given and discussed.


title: PT phase transition in a (2+1)-d relativistic system

abstract: We study a massless Dirac particle with PT symmetric non-Hermitian Rashba interaction in the background of Dirac oscillator potential to show the PT phase transition in a (2+1) dimensional relativistic system analytically. PT phase transition occurs when strength of the (i) imaginary Rashba interaction or (ii) transverse magnetic field exceed their respective critical values. Small mass gap in the spectrum, consistent with other approaches is generated as long as the system is in the unbroken phase. Relativistic Landau levels are constructed explicitly for such a system.


title: Cosmological Hints of Modified Gravity ?

abstract: The recent measurements of Cosmic Microwave Background temperature and polarization anisotropies made by the Planck satellite have provided impressive confirmation of the $\Lambda$CDM cosmological model. However interesting hints of slight deviations from $\Lambda$CDM have been found, including a $95 \%$ c.l. preference for a "modified gravity" structure formation scenario. In this paper we confirm the preference for a modified gravity scenario from Planck 2015 data, find that modified gravity solves the so-called $A_{lens}$ anomaly in the CMB angular spectrum, and constrains the amplitude of matter density fluctuations to $\sigma_8=0.815_{-0.048}^{+0.032}$, in better agreement with weak lensing constraints. Moreover, we find a lower value for the reionization optical depth of $\tau=0.059\pm0.020$ (to be compared with the value of $\tau= 0.079 \pm 0.017$ obtained in the standard scenario), more consistent with recent optical and UV data. We check the stability of this result by considering possible degeneracies with other parameters, including the neutrino effective number, the running of the spectral index and the amount of primordial helium. The indication for modified gravity is still present at about $95\%$ c.l., and could become more significant if lower values of $\tau$ were to be further confirmed by future cosmological and astrophysical data.


title: A quasi-exactly solvable model: two charges in a magnetic field, subject to a non-Coulomb mutual interaction

abstract: We extend the class of QM problems which permit for quasi-exact solutions. Specifically, we consider planar motion of two interacting charges in a constant uniform magnetic field. While Turbiner and Escobar-Ruiz (2013) addressed the case of the Coulomb interaction between the particles, we explore three other potentials. We do this by reducing the appropriate Hamiltonians to the second-order polynomials in the generators of the representation of $SL(2,C)$ group in the differential form. This allows us to perform partial diagonalisation of the Hamiltonian, and to reduce the search for the first few energies and the corresponding wave functions to an algebraic procedure.


title: An Orphan No Longer? Detection of the Southern Orphan Stream and a Candidate Progenitor

abstract: Using a shallow, two-color survey carried out with the Dark Energy Camera, we detect the southern, possibly trailing arm of the Orphan Stream. The stream is reliably detected to a declination of $-38^\circ$, bringing the total known length of the Orphan stream to $108^\circ$. We find a slight offset or "S" shape in the stream at $\delta \simeq -14^\circ$ that would be consistent with the transition from leading to trailing arms. This coincides with a moderate concentration of $137 \pm 25$ stars (to $g = 21.6$) that we consider a possible remnant of the Orphan progenitor. The position of this feature is in agreement with previous predictions.


title: Compositional evolution during rocky protoplanet accretion

abstract: The Earth appears non-chondritic in its abundances of refractory lithophile elements, posing a significant problem for our understanding of its formation and evolution. It has been suggested that this non-chondritic composition may be explained by collisional erosion of differentiated planetesimals of originally chondritic composition. In this work, we present N-body simulations of terrestrial planet formation that track the growth of planetary embryos from planetesimals. We simulate evolution through the runaway and oligarchic growth phases under the Grand Tack model and in the absence of giant planets. These simulations include a state-of-the-art collision model which allows multiple collision outcomes, such as accretion, erosion, and bouncing events, that enables tracking of the evolving core mass fraction of accreting planetesimals. We show that the embryos grown during this intermediate stage of planet formation exhibit a range of core mass fractions, and that with significant dynamical excitation, enough mantle can be stripped from growing embryos to account for the Earth's non-chondritic Fe/Mg ratio. We also find that there is a large diversity in the composition of remnant planetesimals, with both iron-rich and silicate-rich fragments produced via collisions.


title: Profile Stochasticity in PSR J1909-3744

abstract: We extend the recently introduced Bayesian framework `Generative Pulsar Timing Analysis' to incorporate both pulse jitter (high frequency variation in the arrival time of the pulse) and epoch to epoch stochasticity in the shape of the pulse profile. This framework allows for a full timing analysis to be performed on the folded profile data, rather than the site arrival times as is typical in most timing studies. We apply this extended framework both to simulations, and to an 11 yr, 10 cm data set for PSR J1909$-$3744. Using simulations, we show that temporal profile variation can induce timing noise in the residuals that when performing a standard timing analysis is highly covariant with the signal expected from a gravitational wave (GW) background. When working in the profile domain, these variations are de-correlated from the expected GW signal, resulting in significant improvement in the obtained upper limits. Using the PSR J1909$-$3744 data set from the Parkes Pulsar Timing Array project, we find significant evidence for systematic high-frequency profile variation resulting from non-Gaussian noise in the oldest observing system, but no evidence for either detectable pulse jitter, or low-frequency profile shape variation. Using our profile domain framework we therefore obtain upper limits on a red noise process with a spectral index of $\gamma = 13/3$ of $1\times10^{-15}$, consistent with previously published limits.


title: Density Jumps Near the Virial Radius of Galaxy Clusters

abstract: Recent simulations have indicated that the dark matter halos of galaxy clusters should feature steep density jumps near the virial radius. Since the member galaxies are expected to follow similar collisionless dynamics as the dark matter, the galaxy density profile should show such a feature as well. We examine the potential of current datasets to test this prediction by selecting cluster members for a sample of 56 low-redshift (0.1<z<0.3) galaxy clusters, constructing their projected number density profiles, and fitting them with two profiles, one with a steep density jump and one without. Additionally, we investigate the presence of a jump using a non-parametric spline approach. We find that some of these clusters show strong evidence for a model with a density jump. We discuss avenues for further analysis of the density jump with future datasets.


title: Spectro-photometric characterization of high proper motion sources from WISE

abstract: The census of the solar neighborhood is almost complete for stars and becoming more complete in the brown dwarf regime. Spectroscopic, photometric and kinematic characterization of nearby objects helps us to understand the local mass function, the binary fraction, and provides new targets for sensitive planet searches. We aim to derive spectral types and spectro-photometric distances of a sample of new high proper motion sources found with the WISE satellite, and obtain parallaxes for those objects that fall within the area observed by the Vista Variables in the V\'ia L\'actea survey (VVV). We used low resolution spectroscopy and template fitting to derive spectral types, multiwavelength photometry to characterize the companion candidates and obtain photometric distances. Multi-epoch imaging from the VVV survey was used to measure the parallaxes and proper motions for three sources. We confirm a new T2 brown dwarf within $\sim$15 pc. We derived optical spectral types for twenty four sources, mostly M dwarfs within 50 pc. We addressed the wide binary nature of sixteen objects found by the WISE mission and previously known high proper motion sources. Six of these are probably members of wide binaries, two of those are new, and present evidence against the physical binary nature of two candidate binary stars found in the literature, and eight that we selected as possible binary systems. We discuss a likely microlensing event produced by a nearby low mass star and a galaxy, that is to occur in the following five years.


title: Skyrmions in Chiral Magnets with Rashba and Dresselhaus Spin-Orbit Coupling

abstract: Skyrmions are topological spin textures of interest for fundamental science and applications. Previous theoretical studies have focused on systems with broken bulk inversion symmetry, where skyrmions are stabilized by easy-axis anisotropy. We investigate here systems that break surface inversion symmetry, in addition to possible broken bulk inversion. This leads to two distinct Dzyaloshinskii-Moriya (DM) terms with strengths $D_\perp$, arising from Rashba spin-orbit coupling (SOC), and $D_\parallel$ from Dresselhaus SOC. We show that skyrmions become progressively more stable with increasing $D_\perp/D_\parallel$, extending into the regime of easy-plane anisotropy. We find that the spin texture and topological charge density of skyrmions develops nontrivial spatial structure, with quantized topological charge in a unit cell given by a Chern number. Our results give a design principle for tuning Rashba SOC and magnetic anisotropy to stabilize skyrmions in thin films, surfaces, interfaces and bulk magnetic materials that break mirror symmetry.


title: Holographic Consequences of a No Transmission Principle

abstract: Two quantum field theories whose Hilbert spaces do not overlap cannot transmit a signal to one another. From this simple principle, we deduce some highly nontrivial consequences for holographic quantum gravity. These include: (1) certain cosmological bounces are forbidden, (2) generic singularities inside black holes cannot be resolved, and (3) traversable wormholes do not exist. At the classical level, this principle rules out certain types of naked singularities and suggests that new singularity theorems should exist.


title: Bayesian model selection on linear mixed-effects models for comparisons between multiple treatments and a control

abstract: We propose a novel Bayesian model selection technique on linear mixed-effects models to compare multiple treatments with a control. A fully Bayesian approach is implemented to estimate the marginal inclusion probabilities that provide a direct measure of the difference between treatments and the control, along with the model-averaged posterior distributions. Default priors are proposed for model selection incorporating domain knowledge and a component-wise Gibbs sampler is developed for efficient posterior computation. We demonstrate the proposed method based on simulated data and an experimental dataset from a longitudinal study of mouse lifespan and weight trajectories.


title: Bright vigorous winds as signposts of supermassive black hole birth

abstract: The formation of supermassive black holes is still an outstanding question. In the quasi-star scenario, black hole seeds experience an initial super-Eddington growth, that in less than a million years may leave a $10^4-10^5$ M$_{\odot}$ black hole at the centre of a protogalaxy at $z \sim 20-10$. Super-Eddington accretion, however, may be accompanied by vigorous mass loss that can limit the amount of mass that reaches the black hole. In this paper, we critically assess the impact of radiative driven winds, launched from the surface of the massive envelopes from which the black hole accretes. Solving the full wind equations coupled with the hydrostatic structure of the envelope, we find mass outflows with rates between a few tens and $10^4$ M$_{\odot}$ yr$^{-1}$, mainly powered by advection luminosity within the outflow. We therefore confirm the claim by Dotan, Rossi & Shaviv (2011) that mass losses can severely affect the black hole seed early growth within a quasi-star. In particular, seeds with mass $>10^4$ M$_{\odot}$ can only form within mass reservoirs $ \gtrsim 10^7$ M$_{\odot}$, unless they are refilled at huge rates ($ \gtrsim 100$ M$_{\odot}$ yr$^{-1}$). This may imply that only very massive halos ($>10^9$ M$_{\odot}$) at those redshifts can harbour massive seeds. Contrary to previous claims, these winds are expected to be relatively bright ($10^{44}-10^{47}$ erg s$^{-1}$), blue ($T_{\rm eff} \sim 8000$ K) objects, that while eluding the Hubble Space Telescope, could be observed by the James Webb Space Telescope.


title: Excitation Mechanisms for HCN (1-0) and HCO+ (1-0) in Galaxies from the Great Observatories All-sky LIRG Survey

abstract: We present new IRAM 30m spectroscopic observations of the $\sim88$ GHz band, including emission from the CCH (n=1-0) multiplet, HCN (1-0), HCO+ (1-0), and HNC (1-0), for a sample of 58 local luminous and ultraluminous infrared galaxies from the Great Observatories All-sky LIRG Survey (GOALS). By combining our new IRAM data with literature data and Spitzer/IRS spectroscopy, we study the correspondence between these putative tracers of dense gas and the relative contribution of active galactic nuclei (AGN) and star formation to the mid-infrared luminosity of each system. We find the HCN (1-0) emission to be enhanced in AGN-dominated systems ($\langle$L'$_{HCN (1-0)}$/L'$_{HCO^+ (1-0)}\rangle=1.84$), compared to composite and starburst-dominated systems ($\langle$L'$_{HCN (1-0)}$/L'$_{HCO^+ (1-0)}\rangle=1.14$, and 0.88, respectively). However, some composite and starburst systems have L'$_{HCN (1-0)}$/L'$_{HCO^+ (1-0)}$ ratios comparable to those of AGN, indicating that enhanced HCN emission is not uniquely associated with energetically dominant AGN. After removing AGN-dominated systems from the sample, we find a linear relationship (within the uncertainties) between $\log_{10}$(L'$_{HCN (1-0)}$) and $\log_{10}$(L$_{IR}$), consistent with most previous findings. L'$_{HCN (1-0)}$/L$_{IR}$, typically interpreted as the dense gas depletion time, appears to have no systematic trend with L$_{IR}$ for our sample of luminous and ultraluminous infrared galaxies, and has significant scatter. The galaxy-integrated HCN (1-0) and HCO+ (1-0) emission do not appear to have a simple interpretation, in terms of the AGN dominance or the star formation rate, and are likely determined by multiple processes, including density and radiative effects.


title: Description of the Odin Event Extraction Framework and Rule Language

abstract: This document describes the Odin framework, which is a domain-independent platform for developing rule-based event extraction models. Odin aims to be powerful (the rule language allows the modeling of complex syntactic structures) and robust (to recover from syntactic parsing errors, syntactic patterns can be freely mixed with surface, token-based patterns), while remaining simple (some domain grammars can be up and running in minutes), and fast (Odin processes over 100 sentences/second in a real-world domain with over 200 rules). Here we include a thorough definition of the Odin rule language, together with a description of the Odin API in the Scala language, which allows one to apply these rules to arbitrary texts.


title: Astrometric Confirmation and Preliminary Orbital Parameters of the Young Exoplanet 51 Eridani b with the Gemini Planet Imager

abstract: We present new Gemini Planet Imager observations of the young exoplanet 51 Eridani b which provide further evidence that the companion is physically associated with 51 Eridani. Combining this new astrometric measurement with those reported in the literature, we significantly reduce the posterior probability that 51 Eridani b is an unbound foreground or background T-dwarf in a chance alignment with 51 Eridani to $2\times10^{-7}$, an order of magnitude lower than previously reported. If 51 Eridani b is indeed a bound object, then we have detected orbital motion of the planet between the discovery epoch and the latest epoch. By implementing a computationally efficient Monte Carlo technique, preliminary constraints are placed on the orbital parameters of the system. The current set of astrometric measurements suggest an orbital semimajor axis of $14^{+7}_{-3}$ AU, corresponding to a period of $41^{+35}_{-12}$ years (assuming a mass of $1.75$ $M_{\odot}$ for the central star), and an inclination of $138^{+15}_{-13}$ deg. The remaining orbital elements are only marginally constrained by the current measurements. These preliminary values suggest an orbit which does not share the same inclination as the orbit of the distant M-dwarf binary, GJ 3305, which is a wide physically bound companion to 51 Eridani.


title: A highly-ionized region surrounding SN Refsdal revealed by MUSE

abstract: Supernova (SN) Refsdal is the first multiply-imaged, highly-magnified, and spatially-resolved SN ever observed. The SN exploded in a highly-magnified spiral galaxy at z=1.49 behind the Frontier Fields Cluster MACS1149, and provides a unique opportunity to study the environment of SNe at high z. We exploit the time delay between multiple images to determine the properties of the SN and its environment, before, during, and after the SN exploded. We use the integral-field spectrograph MUSE on the VLT to simultaneously target all observed and model-predicted positions of SN Refsdal. We find MgII emission at all positions of SN Refsdal, accompanied by weak FeII* emission at two positions. The measured ratios of [OII] to MgII emission of 10-20 indicate a high degree of ionization with low metallicity. Because the same high degree of ionization is found in all images, and our spatial resolution is too coarse to resolve the region of influence of SN Refsdal, we conclude that this high degree of ionization has been produced by previous SNe or a young and hot stellar population. We find no variability of the [OII] line over a period of 57 days. This suggests that there is no variation in the [OII] luminosity of the SN over this period, or that the SN has a small contribution to the integrated [OII] emission over the scale resolved by our observations.


title: Correspondences between WZNW models and CFTs with W-algebra symmetry

abstract: We study theories with W-algebra symmetries and their relation to WZNW models on (super-)groups. Correlation functions of the WZNW models are expressed in terms of correlators of CFTs with W-algebra symmetry. The symmetries of the theories involved in these correspondences are related by the Drinfeld-Sokolov reduction of Lie algebras to W-algebras. The W-algebras considered in this paper are the Bershadsky-Polyakov algebra for sl(3) and the quasi-superconformal algebra for generic sl(N|M). The quantum W-algebras obtained from affine sl(N) are constructed using embeddings of sl(2) into sl(N), and these can in turn be characterized by partitions of N. The above cases correspond to \underline{N+2} = \underline{2} + N \underline{1} and its supergroup extension. Finally, sl(2N) and the correspondence corresponding to \underline{2N} = N \underline{2} is also analyzed.


title: Is E=mc^2 an exclusively relativistic result?

abstract: The mass-energy formula E=mc^2 is thought to be derived by Einstein from special relativity. The present study shows that since the formula has also been derived from classical physics by Einstein, it is not an exclusively relativistic result. The formula is implied by the classical electromagnetic momentum P=E/c and the Newtonian definition of momentum P=mv. Like momentum P=mv, E=mc^2 applies to both classical physics and special relativity, if relativistic mass is used in the equation. The derivation by Einstein in 1905 is logically flawed as a relativistic proof and the truly relativistic formula should be E=(E0)/(1-v^2/c^2)^(1/2) derived by Laue in 1911 and Klein in 1918. If the energy measured in one reference frame is E0, it is E=(E0)/(1-v^2/c^2)^(1/2) in a reference frame moving at velocity v relative to the first frame.


title: Effects of the sources of reionization on 21-cm redshift-space distortions

abstract: The observed 21-cm signal from the epoch of reionization will be distorted along the line-of-sight by the peculiar velocities of matter particles. These redshift-space distortions will affect the contrast in the signal and will also make it anisotropic. This anisotropy contains information about the cross-correlation between the matter density field and the neutral hydrogen field, and could thus potentially be used to extract information about the sources of reionization. In this paper, we study a collection of simulated reionization scenarios assuming different models for the sources of reionization. We show that the 21-cm anisotropy is best measured by the quadrupole moment of the power spectrum. We find that, unless the properties of the reionization sources are extreme in some way, the quadrupole moment evolves very predictably as a function of global neutral fraction. This predictability implies that redshift-space distortions are not a very sensitive tool for distinguishing between reionization sources. However, the quadrupole moment can be used as a model-independent probe for constraining the reionization history. We show that such measurements can be done to some extent by first-generation instruments such as LOFAR, while the SKA should be able to measure the reionization history using the quadrupole moment of the power spectrum to great accuracy.


title: Coupled attitude and trajectory optimization for a launcher tilting maneuver

abstract: We study the minimum time control problem of the launchers. The optimal trajectories of the problem may contain singular arcs, and thus chattering arcs. The motion of the launcher is described by its attitude kinematics and dynamics and also by its trajectory dynamics. Based on the nature of the system, we implement an efficient indirect numerical method, combined with numerical continuation, to compute numerically the optimal solutions of the problem. Numerical experiments show the efficiency and the robustness of the proposed method.


title: Grain Growth in the Circumstellar Disks of the Young Stars CY Tau and DoAr 25

abstract: We present new results from the Disks@EVLA program for two young stars: CY Tau and DoAr 25. We trace continuum emission arising from their circusmtellar disks from spatially resolved observations, down to tens of AU scales, at {\lambda} = 0.9, 2.8, 8.0, and 9.8 mm for DoAr25 and at {\lambda} = 1.3, 2.8, and 7.1 mm for CY Tau. Additionally, we constrain the amount of emission whose origin is different from thermal dust emission from 5 cm observations. Directly from interferometric data, we find that observations at 7 mm and 1 cm trace emission from a compact disk while millimeter-wave observations trace an extended disk structure. From a physical disk model, where we characterize the disk structure of CY Tau and DoAr 25 at wavelengths shorter than 5 cm, we find that (1) dust continuum emission is optically thin at the observed wavelengths and over the spatial scales studied, (2) a constant value of the dust opacity is not warranted by our observations, and (3) a high-significance radial gradient of the dust opacity spectral index, {\beta}, is consistent with the observed dust emission in both disks, with low-{\beta} in the inner disk and high-{\beta} in the outer disk. Assuming that changes in dust properties arise solely due to changes in the maximum particle size (amax), we constrain radial variations of amax in both disks, from cm-sized particles in the inner disk (R < 40 AU) to millimeter sizes in the outer disk (R > 80 AU). These observational constraints agree with theoretical predictions of the radial-drift barrier, however, fragmentation of dust grains could explain our amax(R) constraints if these disks have lower turbulence and/or if dust can survive high-velocity collisions.


title: Scale-dependent gravitational waves from a rolling axion

abstract: We consider a model in which a pseudo-scalar field $\sigma$ rolls for some e-folds during inflation, sourcing one helicity of a gauge field. These fields are only gravitationally coupled to the inflaton, and therefore produce scalar and tensor primordial perturbations only through gravitational interactions. These sourced signals are localized on modes that exit the horizon while the roll of $\sigma$ is significant. We focus our study on cases in which the model can simultaneously produce (i) a large gravitational wave signal, resulting in observable B-modes of the CMB polarizations, and (ii) sufficiently small scalar perturbations, so to be in agreement with the current limits from temperature anisotropies. Different choice of parameters can instead lead to a localized and visible departure from gaussianity in the scalar sector, either at CMB or LSS scales.


title: $J_{E_T}^{\rm II}$: A Two-prong Jet Finding Algorithm

abstract: We propose a new global jet-finding algorithm for reconstructing two-prong objects like hadronic weak gauge bosons at a hadron collider. The selection of particles in a two-prong jet is required to maximize a $J_{E_T}^{\rm II}$ function, which contains a modified second Fox-Wolfram moment and prefers a two-prong structure for a fixed jet mass. Compared to the traditional jet-substructure method, our algorithm can provide a similar or better performance for identifying boosted weak gauge bosons that are produced from either Standard Model processes or heavy di-boson resonance decays.


title: Model atmospheres of irradiated exoplanets: The influence of stellar parameters, metallicity, and the C/O ratio

abstract: Many parameters constraining the spectral appearance of exoplanets are still poorly understood. We therefore study the properties of irradiated exoplanet atmospheres over a wide parameter range including metallicity, C/O ratio and host spectral type. We calculate a grid of 1-d radiative-convective atmospheres and emission spectra. We perform the calculations with our new Pressure-Temperature Iterator and Spectral Emission Calculator for Planetary Atmospheres (PETIT) code, assuming chemical equilibrium. The atmospheric structures and spectra are made available online. We find that atmospheres of planets with C/O ratios $\sim$ 1 and $T_{\rm eff}$ $\gtrsim$ 1500 K can exhibit inversions due to heating by the alkalis because the main coolants CH$_4$, H$_2$O and HCN are depleted. Therefore, temperature inversions possibly occur without the presence of additional absorbers like TiO and VO. At low temperatures we find that the pressure level of the photosphere strongly influences whether the atmospheric opacity is dominated by either water (for low C/O) or methane (for high C/O), or both (regardless of the C/O). For hot, carbon-rich objects this pressure level governs whether the atmosphere is dominated by methane or HCN. Further we find that host stars of late spectral type lead to planetary atmospheres which have shallower, more isothermal temperature profiles. In agreement with prior work we find that for planets with $T_{\rm eff}$ $<$ 1750 K the transition between water or methane dominated spectra occurs at C/O $\sim$ 0.7, instead of $\sim$ 1, because condensation preferentially removes oxygen.


title: Tatooine Nurseries: Structure and Evolution of Circumbinary Protoplanetary Disks

abstract: Recent discoveries of circumbinary planets by Kepler mission provide motivation for understanding their birthplaces - protoplanetary disks around stellar binaries with separations <1 AU. We explore properties and evolution of such circumbinary disks focusing on modification of their structure caused by tidal coupling to the binary. We develop a set of analytical scaling relations describing viscous evolution of the disk properties, which are verified and calibrated using 1D numerical calculations with realistic inputs. Injection of angular momentum by the central binary suppresses mass accretion onto the binary and causes radial distribution of the viscous angular momentum flux F_J to be different from that in a standard accretion disk around a single star with no torque at the center. Disks with no mass accretion at the center develop F_J profile which is flat in radius. Radial profiles of temperature and surface density are also quite different from those in disks around single stars. Damping of the density waves driven by the binary and viscous dissipation dominate heating of the inner disk (within 1-2 AU), pushing the iceline beyond 3-5 AU, depending on disk mass and age. Irradiation by the binary governs disk thermodynamics beyond ~10 AU. However, self-shadowing by the hot inner disk may render central illumination irrelevant out to ~20 AU. Spectral energy distribution of a circumbinary disk exhibits a distinctive bump around 10 micron, which may facilitate identification of such disks around unresolved binaries. Efficient tidal coupling to the disk drives orbital inspiral of the binary and may cause low-mass and compact binaries to merge into a single star within the disk lifetime. We generally find that circumbinary disks present favorable sites for planet formation (despite wider zone of volatile depletion), in agreement with the statistics of Kepler circumbinary planets.


title: Dark Photons from the Center of the Earth: Smoking-Gun Signals of Dark Matter

abstract: Dark matter may be charged under dark electromagnetism with a dark photon that kinetically mixes with the Standard Model photon. In this framework, dark matter will collect at the center of the Earth and annihilate into dark photons, which may reach the surface of the Earth and decay into observable particles. We determine the resulting signal rates, including Sommerfeld enhancements, which play an important role in bringing the Earth's dark matter population to their maximal, equilibrium value. For dark matter masses $m_X \sim$ 100 GeV - 10 TeV, dark photon masses $m_{A'} \sim$ MeV - GeV, and kinetic mixing parameters $\varepsilon \sim 10^{-10} - 10^{-8}$, the resulting electrons, muons, photons, and hadrons that point back to the center of the Earth are a smoking-gun signal of dark matter that may be detected by a variety of experiments, including neutrino telescopes, such as IceCube, and space-based cosmic ray detectors, such as Fermi-LAT and AMS. We determine the signal rates and characteristics, and show that large and striking signals---such as parallel muon tracks---are possible in regions of the $(m_{A'}, \varepsilon)$ plane that are not probed by direct detection, accelerator experiments, or astrophysical observations.


title: A brief note on the Karhunen-Lo\`eve expansion

abstract: We provide a detailed derivation of the Karhunen-Lo\`eve expansion of a stochastic process. We also discuss briefly Gaussian processes, and provide a simple numerical study for the purpose of illustration.


title: On the overlap distribution of branching random walks

abstract: In this paper, we study the overlap distribution and Gibbs measure of the Branching Random Walk with Gaussian increments on a binary tree. We first prove that the Branching Random Walk is 1 step Replica Symmetry Breaking and give a precise form for its overlap distribution, verifying a prediction of Derrida and Spohn. We then prove that the Gibbs measure of this system satisfies the Ghirlanda-Guerra identities. As a consequence, the limiting Gibbs measure has Poisson-Dirichlet statistics. The main technical result is a proof that the overlap distribution for the Branching Random Walk is supported on the set $\{0,1\}$.


title: Standard Bases in mixed Power Series and Polynomial Rings over Rings

abstract: In this paper we study standard bases for submodules of a mixed power series and polynomial ring $R[[t_1,\ldots,t_m]][x_1,\ldots,x_n]^s$ respectively of their localization with respect to a $t$-local monomial ordering for a certain class of noetherian rings $R$. The main steps are to prove the existence of a division with remainder generalizing and combining the division theorems of Grauert--Hironaka and Mora and to generalize the Buchberger criterion. Everything else then translates naturally. Setting either $m=0$ or $n=0$ we get standard bases for polynomial rings respectively for power series rings over $R$ as a special case.


title: VLBI for Gravity Probe B: The Guide Star IM Pegasi

abstract: We review the radio very long baseline interferometry (VLBI) observations of the guide star, IM Peg, and three compact extragalactic reference sources, 3C 454.3, B2250+194, and B2252+172, made in support of the NASA/Stanford gyroscope relativity mission, GP-B. The main goal of the observations was the determination of the proper motion of IM Peg relative to the distant universe. VLBI observations made between 1997 and 2005 yield a proper motion of IM Peg of -20.83 $\pm$ 0.09 mas yr$^{-1}$ in RA and -27.27 $\pm$ 0.09 mas yr$^{-1}$ in dec, in a celestial reference frame of extragalactic radio galaxies and quasars virtually identical to the International Celestial Reference Frame 2 (ICRF2). They also yield a parallax for IM Peg of 10.37 $\pm$ 0.07 mas, corresponding to a distance of 96.4 $\pm$ 0.7 pc. The uncertainties are standard errors with statistical and estimated systematic contributions added in quadrature. These results met the pre-launch requirements of the GP-B mission to not discernibly degrade the estimates of the geodetic and frame-dragging effects. The paper also reports on a 1$\sigma$ upper limit on the magnitude of the components of the proper motion of the 'core' of 3C 454.3 relative to the ICRF2 of 46 and 56 $\mu$as yr$^{-1}$ in RA and dec, respectively, and presents densely sampled ellipses of the parallax and the orbit of the giant of the binary system. It further gives a sequence of images of the flickering radio emission relative to the disk of the giant. For a 'movie of a star,' see http://www.yorku.ca/bartel/impeg.mpg .


title: Finite-Type Invariants of Order One for Framed Virtual Knots

abstract: A. Henrich proved the existence of the universal finite-type invariant of order one for virtual knots. We extend the construction and the methods of her paper to framed virtual knots. To do so, we introduce the notions of virtual strings and based matrices for framed flat virtual knots.


title: A practical method for the analysis of meteor spectra

abstract: The analysis of meteor spectra (photographic, CCD or video recording) is complicated by the fact that spectra obtained with objective gratings are curved and have a nonlinear dispersion. In this paper it is shown that with a simple image transformation the spectra can be linearized in such a way that individual spectra over the whole image plane are parallel and have a constant, linear dispersion. This simplifies the identification and measurement of meteor spectral lines. A practical method is given to determine the required image transformation.


title: Aharonov-Bohm oscillations caused by non-topological surface states in Dirac nanowires

abstract: One intriguing fingerprint of surface states in topological insulators is the Aharonov-Bohm effect in magnetoconductivity of nanowires. We show that surface states in nanowires of Dirac materials (bismuth, bismuth antimony, and lead tin chalcogenides) being in non-topological phase, exhibit the same effect as amendment to magnetoconductivity of the bulk states. We consider a simple model of a cylindrical nanowire, which is described by the 3D Dirac equation with a general $T$-invariant boundary condition. The boundary condition is determined by a single phenomenological parameter whose sign defines topological-like and non-topological surface states. The non-topological surface states emerge outside the gap. In longitudinal magnetic field $B$ they lead to Aharonov-Bohm amendment for the density of states and correspondingly for conductivity of the nanowire. The phase of these magnetooscillations increases with $B$ from $\pi$ to $2\pi$.


title: Graph games and the pizza problem

abstract: We propose a class of two person perfect information games based on weighted graphs. One of these games can be described in terms of a round pizza which is cut radially into pieces of varying size. The two players alternately take pieces subject to the following rule: Once the first piece has been chosen, all subsequent selections must be adjacent to the hole left by the previously taken pieces. Each player tries to get as much pizza as possible. The original pizza problem was to settle the conjecture that Player One can always get at least half of the pizza. The conjecture turned out to be false. Our main result is a complete solution of a somewhat simpler class of games, concatenations of stacks and two-ended stacks, and we provide a linear time algorithm for this. The algorithm and its output can be described without reference to games. It produces a certain kind of partition of a given finite sequence of real numbers. The conditions on the partition involve alternating sums of various segments of the given sequence. We do not know whether these partitions have applications outside of game theory. The algorithm leads to a quadratic time algorithm which gives the value and an optimal initial move for pizza games. We also provide some general theory concerning the semigroup of equivalence classes of graph games.


title: Cosmological nonlinear density and velocity power spectra including nonlinear vector and tensor modes

abstract: We present the leading order non-linear density and velocity power spectra in the complete form; previous studies have omitted the vector- and tensor-type perturbations simultaneously excited by the scalar-type perturbation in nonlinear order. These additional contributions are comparable to the scalar-type purely relativistic perturbations, and thus negligible in the current paradigm of concordance cosmology: i.e., concerning density and velocity perturbations of the pressureless matter in perturbation regime well inside of matter-dominated epoch, we show that pure Einstein's gravity contributions appearing from the third order are entirely negligible (five orders of magnitude smaller than the Newtonian contributions) in all scales. We thus prove that Newtonian perturbation theory is quite reliable in calculating the amplitude of matter fluctuations even in the precision era of cosmology. Therefore, besides the ones imprinted as the initial condition generated in the earlier phase, the other relativistic effect relevant for interpreting observational data must be the projection effect that occurs when mapping galaxies on to the observed coordinate.


title: Bayesian Nonparametric Graph Clustering

abstract: We present clustering methods for multivariate data exploiting the underlying geometry of the graphical structure between variables. As opposed to standard approaches that assume known graph structures, we first estimate the edge structure of the unknown graph using Bayesian neighborhood selection approaches, wherein we account for the uncertainty of graphical structure learning through model-averaged estimates of the suitable parameters. Subsequently, we develop a nonparametric graph clustering model on the lower dimensional projections of the graph based on Laplacian embeddings using Dirichlet process mixture models. In contrast to standard algorithmic approaches, this fully probabilistic approach allows incorporation of uncertainty in estimation and inference for both graph structure learning and clustering. More importantly, we formalize the arguments for Laplacian embeddings as suitable projections for graph clustering by providing theoretical support for the consistency of the eigenspace of the estimated graph Laplacians. We develop fast computational algorithms that allow our methods to scale to large number of nodes. Through extensive simulations we compare our clustering performance with standard clustering methods. We apply our methods to a novel pan-cancer proteomic data set, and evaluate protein networks and clusters across multiple different cancer types.


title: Torsion Wave Solutions in Yang-Mielke Theory of Gravity

abstract: The approach of metric-affine gravity initially distinguishes it from Einstein's general relativity. Using an independent affine connection produces a theory with 10+64 unknowns. We write down the Yang-Mills action for the affine connection and produce the Yang-Mills equation and the so called complementary Yang-Mills equation by independently varying with respect to the connection and the metric respectively. We call this theory the Yang-Mielke theory of gravity. We construct explicit spacetimes with pp-metric and purely axial torsion and show that they represent a solution of Yang-Mills theory. Finally we compare these spacetimes to existing solutions of metric-affine gravity and present future research possibilities.


title: Optical Control of Fluorescence through Plasmonic Eigenmode Extinction

abstract: We introduce the concept of optical control of the fluorescence yield of CdSe quantum dots through plasmon-induced structural changes in random semicontinuous nanostructured gold films. We demonstrate that the wavelength- and polarization dependent coupling between quantum dots and the semicontinuous films, and thus the fluorescent emission spectrum, can be controlled and significantly increased through the optical extinction of a selective band of eigenmodes in the films. This optical method of effecting controlled changes in the metal nanostructure allows for versatile functionality in a single sample and opens a pathway to in situ control over the fluorescence spectrum.


title: Design Aspects of Short Range Millimeter Wave Networks: A MAC Layer Perspective

abstract: Increased density of wireless devices, ever growing demands for extremely high data rate, and spectrum scarcity at microwave bands make the millimeter wave (mmWave) frequencies an important player in future wireless networks. However, mmWave communication systems exhibit severe attenuation, blockage, deafness, and may need microwave networks for coordination and fall-back support. To compensate for high attenuation, mmWave systems exploit highly directional operation, which in turn substantially reduces the interference footprint. The significant differences between mmWave networks and legacy communication technologies challenge the classical design approaches, especially at the medium access control (MAC) layer, which has received comparatively less attention than PHY and propagation issues in the literature so far. In this paper, the MAC layer design aspects of short range mmWave networks are discussed. In particular, we explain why current mmWave standards fail to fully exploit the potential advantages of short range mmWave technology, and argue for the necessity of new collision-aware hybrid resource allocation frameworks with on-demand control messages, the advantages of a collision notification message, and the potential of multihop communication to provide reliable mmWave connections.


title: The Dehn functions of Stallings-Bieri groups

abstract: We show that the Stallings-Bieri groups, along with certain other Bestvina-Brady groups, have quadratic Dehn function.


title: Surface Chemical Reactivity of Ultrathin Pd(111) Films on Ru(0001): Importance of Orbital Symmetry in the Application of the d-Band Model

abstract: The chemical bonding of adsorbate molecules on transition-metal surfaces is strongly influenced by the hybridization between the molecular orbitals and the metal d-band. The strength of this interaction is often correlated with the location of the metal d-band center relative to the Fermi level. Here, we exploit finite size effects in the electronic structure of ultrathin Pd(111) films grown on Ru(0001) to tune their reactivity by changing the film thickness one atom layer at a time, while keeping all other variables unchanged. Interestingly, while bulk Pd(111) is reactive towards oxygen, Pd(111) films below five monolayers are surprisingly inert. This observation is fully in line with the d-band model prediction when applied to the orbitals involved in the bonding. The shift of the d-band center with film thickness is primarily attributed to shifts in the partial density of states associated with the 4dxz and 4dyz orbitals. This study gives an in-depth look into the orbital specific contributions to the surface chemical reactivity, providing new insights that could be useful in surface catalysis.


title: Charged-particle multiplicities in proton-proton collisions at $\sqrt{s}$ = 0.9 to 8 TeV

abstract: A detailed study of pseudorapidity densities and multiplicity distributions of primary charged particles produced in proton-proton collisions, at $\sqrt{s} =$ 0.9, 2.36, 2.76, 7 and 8 TeV, in the pseudorapidity range $|\eta|<2$, was carried out using the ALICE detector. Measurements were obtained for three event classes: inelastic, non-single diffractive and events with at least one charged particle in the pseudorapidity interval $|\eta|<1$. The use of an improved track-counting algorithm combined with ALICE's measurements of diffractive processes allows a higher precision compared to our previous publications. A KNO scaling study was performed in the pseudorapidity intervals $|\eta|<$ 0.5, 1.0 and 1.5. The data are compared to other experimental results and to models as implemented in Monte Carlo event generators PHOJET and recent tunes of PYTHIA6, PYTHIA8 and EPOS.


title: Accurate calculation of Stokes drag for point-particle tracking in two-way coupled flows

abstract: In this work, we propose and test a method for calculating Stokes drag applicable to particle-laden fluid flows where two-way momentum coupling is important. In the point-particle formulation, particle dynamics are coupled to fluid dynamics via a source term that appears in the respective momentum equations. When the particle Reynolds number is small and the particle diameter is smaller than the fluid scales, it is common to approximate the momentum coupling source term as the Stokes drag. The Stokes drag force depends on the difference between the undisturbed fluid velocity evaluated at the particle location, and the particle velocity. However, owing to two-way coupling, the fluid velocity is modified in the neighborhood of a particle, relative to its undisturbed value. This causes the computed Stokes drag force to be underestimated in two-way coupled point-particle simulations. We develop estimates for the drag force error as function of the particle size relative to the grid size. We then develop a correction method that estimates the undisturbed fluid velocity from the computed disturbed velocity field. The correction scheme is tested for a particle settling in an otherwise quiescent fluid and is found to reduce the error in computed settling velocity by an order of magnitude compared with common interpolation schemes.


title: On Optimizing Human-Machine Task Assignments

abstract: When crowdsourcing systems are used in combination with machine inference systems in the real world, they benefit the most when the machine system is deeply integrated with the crowd workers. However, if researchers wish to integrate the crowd with "off-the-shelf" machine classifiers, this deep integration is not always possible. This work explores two strategies to increase accuracy and decrease cost under this setting. First, we show that reordering tasks presented to the human can create a significant accuracy improvement. Further, we show that greedily choosing parameters to maximize machine accuracy is sub-optimal, and joint optimization of the combined system improves performance.


title: Possible Effects of Fierz Transformations on Vacua of Some Four-Fermion Interaction Models

abstract: A theoretical research on possible effects of the Fierz transformations on the ground states (vacua) of some 2-flavor and $N_c$-color four-fermion (quark) interaction models has been systematically conducted. It has been shown that, based on the known criterions of the interplay between the antiquark-quark ($\bar{q}$-$q$) and diquark ($q$-$q$) condensates, in 4D space-time, for the given $\bar{q}$-$q$ channel couplings with chiral symmetry and from the heavy gluon exchange, the effects of the Fierz transformations are not enough to change the feature that the models' vacua could only be in the pure $\bar{q}$-$q$ condensate phases. However, for a given pure scalar $q$-$q$ channel coupling with the strength $H_S$, the Fierz transformations will lead to the nontrivial effect that the model's vacuum could be in the expected $q$-$q$ condensate phase only if $N_c<9$ and $H_S$ is small, and as the increase of $N_c$ and $H_S$, the vacuum will get first in a coexistence phase with $q$-$q$ and $\bar{q}$-$q$ condensates then in a pure $\bar{q}$-$q$ condensate phase until $N_c\rightarrow\infty$. Similar conclusions are also drawn from relevant four-fermion interaction models in 2D and 3D space-time. The general significance of the research is indicated.


title: Asymptotic properties of infinite directed unions of local quadratic transforms

abstract: We consider infinite sequences {R_n} of successive local quadratic transforms of a regular local ring. Let S denote the directed union of the sequence of regular local rings R_n. We previously showed the existence of a unique limit point V of the family of order valuation rings of the sequence. In this paper, we examine asymptotic properties of this family of order valuations. We link this asymptotic behavior to ring-theoretic properties of S, namely whether S is archimedean and whether S is completely integrally closed. We construct examples of such S that are archimedean and completely integrally closed but not valuation domains. We give an explicit description of V, where the description depends on whether S is archimedean or non-archimedean.


title: Benchmark Test of Differential Emission Measure Codes and Multi-Thermal Energies in Solar Active Regions

abstract: We compare the ability of 11 Differential Emission Measure (DEM) forward-fitting and inversion methods to constrain the properties of active regions and solar flares by simulating synthetic data using the instrumental response functions of SDO/AIA, SDO/EVE, RHESSI, and GOES/XRS. The codes include the single-Gaussian DEM, a bi-Gaussian DEM, a fixed-Gaussian DEM, a linear spline DEM, the spatial synthesis DEM, the Monte-Carlo Markov chain DEM, the regularized DEM inversion, the Hinode/XRT method, a polynomial spline DEM, an EVE+GOES, and an EVE+RHESSI method. Averaging the results from all 11 DEM methods, we find the following accuracies in the inversion of physical parameters: the EM-weighted temperature $T_w^{fit}/T_w^{sim}=0.9\pm0.1$, the peak emission measure $EM_p^{fit}/EM_p^{sim}=0.6\pm0.2$, the total emission measure $EM_t^{fit}/EM_t^{sim}=0.8\pm0.3$, and the multi-thermal energies $E_{th}^{fit}/EM_{th}^{sim}=1.2\pm0.4$. We find that the AIA spatial synthesis, the EVE+GOES, and the EVE+RHESSI method yield the most accurate results.


title: Quantum diffusion of a relativistic particle in a time-dependent random potential

abstract: We present a rigorous study of quantum diffusion of a relativistic particle subjected to a time-dependent random potential with $\delta$ correlation in time. We find that in the asymptotic time limit the particle wave packet spreads ballistically in contrast with the nonrelativistic case, which in the same situation exhibits superballistic diffusion. The relativistic suppression of wave packet diffusion is discussed in connection with statistical conservation laws that follow from relativistic dynamics.


title: End-point estimates for singular integrals with non-smooth kernels on product spaces

abstract: The main aim of this article is to establish boundedness of singular integrals with non-smooth kernels on product spaces. Let $L_1$ and $L_2$ be non-negative self-adjoint operators on $L^2(\mathbb{R}^{n_1})$ and $L^2(\mathbb{R}^{n_2})$, respectively, whose heat kernels satisfy Gaussian upper bounds. First, we obtain an atomic decomposition for functions in $H^1_{L_1,L_2}(\mathbb{R}^{n_1}\times\mathbb{R}^{n_2})$ where the Hardy space $H^1_{L_1,L_2}(\mathbb{R}^{n_1}\times\mathbb{R}^{n_2})$ associated with $L_1$ and $L_2$ is defined by square function norms, then prove an interpolation property for this space. Next, we establish sufficient conditions for certain singular integral operators to be bounded on the Hardy space $H^1_{L_1,L_2}(\mathbb{R}^{n_1}\times\mathbb{R}^{n_2})$ when the associated kernels of these singular integrals only satisfy regularity conditions significantly weaker than those of the standard Calder\'on--Zygmund kernels. As applications, we obtain endpoint estimates of the double Riesz transforms associated to Schr\"dingier operators and a Marcinkiewicz-type spectral multiplier theorem for non-negative self-adjoint operators on product spaces.


title: Critical Phenomena of Dynamical Delocalization in Quantum Anderson Map

abstract: Using a quantum map version of one-dimensional Anderson model, the localization-delocalization transition of quantum diffusion induced by coherent dynamical perturbation is investigated in comparison with quantum standard map. Existence of critical phenomena, which depends on the number of frequency component $M$, is demonstrated. Diffusion exponents agree with theoretical prediction for the transition, but the critical exponent of the localization length deviates from it with increase in the $M$. The critical power $\epsilon_c$ of the normalized perturbation at the transition point remarkably decreases as $\epsilon_c \sim (M-1)^{-1}$.


title: Spectral Gap and Edge Excitations of $d$-dimensional PVBS models on half-spaces

abstract: We analyze a class of quantum spin models defined on half-spaces in the $d$-dimensional hypercubic lattice bounded by a hyperplane with inward unit normal vector $m\in\mathbb{R}^d$. The family of models was previously introduced as the single species Product Vacua with Boundary States (PVBS) model, which is a spin-$1/2$ model with a XXZ-type nearest neighbor interactions depending on parameters $\lambda_j\in (0,\infty)$, one for each coordinate direction. For any given values of the parameters, we prove an upper bound for the spectral gap above the unique ground state of these models, which vanishes for exactly one direction of the normal vector $m$. For all other choices of $m$ we derive a positive lower bound of the spectral gap, except for the case $\lambda_1 =\cdots =\lambda_d=1$, which is known to have gapless excitations in the bulk.


title: Flavor tagging TeV jets for BSM and QCD

abstract: We present a new scheme for tagging high-$p_{T}$ bottom and charm jets using energetic muons. Contemporary track-based $b$ tags lose their ability to reject light jet background as jet $p_{T}\rightarrow\mathcal{O}(\mathrm{TeV})$, where the massive boost exposes fundamental limits in tracking resolution. For our "$\mu_{x}$" tag, the signal efficiency and light jet rejection is robust versus $p_{T}$. In the tested regime (jet $p_{T}\in[\mathrm{0.5,2.1}]$ TeV), $\mu_{x}$ tags $\sim14\%$ of $b$ jets, $\sim6.5\%$ of $c$ jets and $\sim0.65\%$ of light jets. Since $\mu_{x}$ tagging should be immediately useful in a searches for heavy resonances, we test it with a typical dijet search --- a heavy, leptophobic $Z^{\prime}$.


title: Formal Support for Standardizing Protocols with State

abstract: Many cryptographic protocols are designed to achieve their goals using only messages passed over an open network. Numerous tools, based on well-understood foundations, exist for the design and analysis of protocols that rely purely on message passing. However, these tools encounter difficulties when faced with protocols that rely on non-local, mutable state to coordinate several local sessions. We adapt one of these tools, {\cpsa}, to provide automated support for reasoning about state. We use Ryan's Envelope Protocol as an example to demonstrate how the message-passing reasoning can be integrated with state reasoning to yield interesting and powerful results. Keywords: protocol analysis tools, stateful protocols, TPM, PKCS#11.


title: Linear-time Learning on Distributions with Approximate Kernel Embeddings

abstract: Many interesting machine learning problems are best posed by considering instances that are distributions, or sample sets drawn from distributions. Previous work devoted to machine learning tasks with distributional inputs has done so through pairwise kernel evaluations between pdfs (or sample sets). While such an approach is fine for smaller datasets, the computation of an $N \times N$ Gram matrix is prohibitive in large datasets. Recent scalable estimators that work over pdfs have done so only with kernels that use Euclidean metrics, like the $L_2$ distance. However, there are a myriad of other useful metrics available, such as total variation, Hellinger distance, and the Jensen-Shannon divergence. This work develops the first random features for pdfs whose dot product approximates kernels using these non-Euclidean metrics, allowing estimators using such kernels to scale to large datasets by working in a primal space, without computing large Gram matrices. We provide an analysis of the approximation error in using our proposed random features and show empirically the quality of our approximation both in estimating a Gram matrix and in solving learning tasks in real-world and synthetic data.


title: Collinear Twist-3 Approach to Transverse Single-Spin Asymmetry in Proton-Proton Collision

abstract: We present our recent analysis on the transverse single-spin asymmetry (SSA) in inclusive pion and direct-photon production in pp collisions for RHIC kinematics. The analysis includes the contributions from twist-3 quark-gluon-quark correlations in the proton and twist-3 fragmentation effects for the pion. Some of the functions appearing in the formula for the SSA, such as the soft-gluon-pole Qiu-Sterman function, the nucleon transversity, and the Collins function were fixed consistently with the SSA data in semi-inclusive DIS and in e+e- annihilation, so that our analysis is free from the sign mismatch problem.


title: Enceladus's measured physical libration requires a global subsurface ocean

abstract: Several planetary satellites apparently have subsurface seas that are of great interest for, among other reasons, their possible habitability. The geologically diverse Saturnian satellite Enceladus vigorously vents liquid water and vapor from fractures within a south polar depression and thus must have a liquid reservoir or active melting. However, the extent and location of any subsurface liquid region is not directly observable. We use measurements of control points across the surface of Enceladus accumulated over seven years of spacecraft observations to determine the satellite's precise rotation state, finding a forced physical libration of 0.120 $\pm$ 0.014{\deg} (2{\sigma}). This value is too large to be consistent with Enceladus's core being rigidly connected to its surface, and thus implies the presence of a global ocean rather than a localized polar sea. The maintenance of a global ocean within Enceladus is problematic according to many thermal models and so may constrain satellite properties or require a surprisingly dissipative Saturn.


title: Rapid Doubling of the Critical Current of YBa$_2$Cu$_3$O$_{7-\delta}$ Coated Conductors for Viable High-Speed Industrial Processing

abstract: We demonstrate that 3.5-MeV oxygen irradiation can markedly enhance the in-field critical current of commercial 2nd generation superconducting tapes with an exposure time of just one second per 0.8 cm2. The speed demonstrated here is now at the level required for an industrial reel-to-reel post-processing. The irradiation is made on production line samples through the protective silver coating and does not require any modification of the growth process. From TEM imaging, we identify small clusters as the main source of increased vortex pinning.


title: Induced fermionic current at finite temperature

abstract: Here we analyze the finite temperature expectation values of the charge and current densities for a massive fermionic quantum field with nonzero chemical potential, $\mu$, induced by a magnetic flux running along the axis of an idealized cosmic string. These densities are decomposed into the vacuum expectation values and contributions coming from the particles and antiparticles. Specifically the charge density is an even periodic function of the magnetic flux with the period equal to the quantum flux and an odd function of the chemical potential. The only nonzero component of the current density corresponds to the azimuthal current and it is an odd periodic function of the magnetic flux and an even function of the chemical potential. Both analyzed are developed for the cases where $|\mu |$ is smaller than the mass of the field quanta, $m$.


title: The Hausdorff dimension of quasi-circles: a result of Ruelle and Bowen

abstract: We provide an expanded and clarified proof of the famous result of Bowen and Ruelle giving an asymptotic formula for the Hausdorff dimension of quasi-circles corresponding to the Julia sets of $f(z)=z^2+c$ for small $c$. The proof does not contain new material but has been rewritten to make it more accessible to MSc or PhD students with an interest in dimension theory.


title: Dehn surgeries and rational homology balls

abstract: We consider the question of which Dehn surgeries along a given knot bound rational homology balls. We use Ozsv\'ath and Szab\'o's correction terms in Heegaard Floer homology to obtain general constraints on the surgery coefficients. We then turn our attention to the case of integral surgeries, with particular emphasis on positive torus knots. Finally, combining these results with a lattice-theoretic obstruction based on Donaldon's theorem, we classify which integral surgeries along torus knots of the form $T_{kq\pm 1,q}$ bound rational homology balls.


title: A robust DPG method for singularly perturbed reaction-diffusion problems

abstract: We present and analyze a discontinuous Petrov-Galerkin method with optimal test functions for a reaction-dominated diffusion problem in two and three space dimensions. We start with an ultra-weak formulation that comprises parameters $\alpha$, $\beta$ to allow for general $\varepsilon$-dependent weightings of three field variables ($\varepsilon$ being the small diffusion parameter). Specific values of $\alpha$ and $\beta$ imply robustness of the method, that is, a quasi-optimal error estimate with a constant that is independent of $\varepsilon$. Moreover, these values lead to a norm for the field variables that is known to be balanced in $\varepsilon$ for model problems with typical boundary layers. Several numerical examples underline our theoretical estimates and reveal stability of approximations even for very small $\varepsilon$.


title: Surveying the TeV sky with HAWC

abstract: The High altitude Water Cherenkov (HAWC) Observatory has been completed and began full operation in early 2015. Located at an elevation of 4,100 m near the Sierra Negra volcano in the state of Puebla, Mexico, HAWC consists of 300 water tanks instrumented with 4 PMTs each. The array is optimized for detecting air showers produced by gamma rays with energies between 100 GeV and 100 TeV and can also be used to measure charged cosmic rays. A wide instantaneous field of view of ~2 steradians and a duty cycle >95% allow HAWC to survey two-thirds of the sky every day. These unique capabilities make it possible to monitor variable gamma-ray fluxes and search for gamma-ray bursts and other transient events, providing new insights into particle acceleration in galactic and extra-galactic sources. In this contribution, we will present first results from more than one year of observations with a partial array configuration. We will discuss how HAWC can map the gamma-ray sky as well as probe other physics including cosmic ray anisotropies and the search for signatures of dark matter annihilation.


title: Extending cosmological tests of General Relativity with the Square Kilometre Array

abstract: Tests of general relativity (GR) are still in their infancy on cosmological scales, but forthcoming experiments promise to greatly improve their precision over a wide range of distance scales and redshifts. One such experiment, the Square Kilometre Array (SKA), will carry out several wide and deep surveys of resolved and unresolved neutral hydrogen (HI) 21cm line-emitting galaxies, mapping a significant fraction of the sky from $0 \le z \lesssim 6$. I present forecasts for the ability of a suite of possible SKA HI surveys to detect deviations from GR by reconstructing the cosmic expansion and growth history. SKA Phase 1 intensity mapping surveys can achieve sub-1\% measurements of $f\sigma_8$ out to $z\approx 1$, with an SKA1-MID Band 2 survey out to $z \lesssim 0.6$ able to surpass contemporary spectroscopic galaxy surveys such as DESI and Euclid in terms of constraints on modified gravity parameters if challenges such as foreground contamination can be tackled effectively. A more futuristic Phase 2 HI survey of $\sim10^9$ spectroscopic galaxy redshifts would be capable of detecting a $\sim 2\%$ modification of the Poisson equation out to $z\approx 2$.


title: Dynamical transition between two mesons and a tetraquark

abstract: We consider a system composed of two identical light quarks ($qq$) and two identical antiquarks ($\bar Q\bar Q$) that can be linked either as two mesons or as a tetraquark, incorporating quantum correlations between identical particles and an effective many-body potential between particles. We perform a 3-D Monte Carlo simulation of the system, considering the configurations allowed to form: i) Only two mesons, ii) Only tetraquark and iii) two mesons and tetraquark . We characterize each case and determine whether it is energetically more favorable to form a tetraquark or two mesons, as a function of the interparticle separation distance which, for a fixed number of particles, can be identified as a particle density. We determine how the two mesons, which dominate the low density regime, mixes with a tetraquark state as the density increases. Properties like the mean square radius and the two-particle correlation function are found to reflect such transition, and we provide a parameterization of the diquark correlation function in the isolated case. We track the dynamical flipping among configurations to determine the recombination probability, exhibiting the importance of the tetraquark state. We analize the four-body potential evolution and show that its linear behavior is preserved, although the slope can reflect the presence of a mixed state. Results are shown for several light-quarks to heavy-antiquarks mass ratios whenever they are found to be relevant.


title: Reconsidering No-Go Theorems from a Practical Perspective

abstract: I argue that our judgements regarding the locally causal models which are compatible with a given quantum no-go theorem implicitly depend, in part, on the context of inquiry. It follows from this that certain no-go theorems, which are particularly striking in the traditional foundational context, have no force when the context switches to a discussion of the physical systems we are capable of building with the aim of classically reproducing quantum statistics. I close with a general discussion of the possible implications of this for our understanding of the limits of classical description, and for our understanding of the fundamental aim of physical investigation.


title: Moment estimates implied by modified log-Sobolev inequalities

abstract: We study a class of logarithmic Sobolev inequalities with a general form of the energy functional. The class generalizes various examples of modified logarithmic Sobolev inequalities considered previously in the literature. Refining a method of Aida and Stroock for the classical logarithmic Sobolev inequality, we prove that if a measure on $\mathbb{R}^n$ satisfies a modified logarithmic Sobolev inequality then it satisfies a family of $L^p$-Sobolev-type inequalities with non-Euclidean norms of gradients (and dimension-independent constants). The latter are shown to yield various concentration-type estimates for deviations of smooth (not necessarily Lipschitz) functions and measures of enlargements of sets corresponding to non-Euclidean norms. We also prove a two-level concentration result for functions of bounded Hessian and measures satisfying the classical logarithmic Sobolev inequality.


title: Detecting Sparse Mixtures: Rate of Decay of Error Probability

abstract: We study the rate of decay of the probability of error for distinguishing between a sparse signal with noise, modeled as a sparse mixture, from pure noise. This problem has many applications in signal processing, evolutionary biology, bioinformatics, astrophysics and feature selection for machine learning. We let the mixture probability tend to zero as the number of observations tends to infinity and derive oracle rates at which the error probability can be driven to zero for a general class of signal and noise distributions via the likelihood ratio test. In contrast to the problem of detection of non-sparse signals, we see the log-probability of error decays sublinearly rather than linearly and is characterized through the $\chi^2$-divergence rather than the Kullback-Leibler divergence for "weak" signals and can be independent of divergence for "strong" signals. Our contribution is the first characterization of the rate of decay of the error probability for this problem for both the false alarm and miss probabilities.


title: Majority Digraphs

abstract: A majority digraph is a finite simple digraph $G=(V,\to)$ such that there exist finite sets $A_v$ for the vertices $v\in V$ with the following property: $u\to v$ if and only if "more than half of the $A_u$ are $A_v$". That is, $u\to v$ if and only if $ |A_u \cap A_v | > \frac{1}{2} \cdot |A_u|$. We characterize the majority digraphs as the digraphs with the property that every directed cycle has a reversal. If we change $\frac{1}{2}$ to any real number $\alpha\in (0,1)$, we obtain the same class of digraphs. We apply the characterization result to obtain a result on the logic of assertions "most $X$ are $Y$" and the standard connectives of propositional logic.


title: Additional strange resonances from Lattice QCD

abstract: Recent Lattice QCD (LQCD) studies suggest that there are missing resonances in the strange sector of the Hadron Resonance Gas (HRG) model. By adopting the continuous Hagedorn mass spectrum, we present how different medium compositions influence the HRG predictions of conserved charge fluctuations. It is shown that missing strange resonances may be partially accounted for by applying the Hagedorn mass spectrum extracted from experimentally established hadrons. On the other hand, the strange-baryonic spectra, extracted from LQCD results for fluctuations, are found to be consistent with the unconfirmed states in the Particle Data Group (PDG) database, whilst the strange-mesonic spectrum points towards yet undiscovered states in the intermediate mass region.


title: Design and control of a soft, shape-changing, crawling robot

abstract: Soft materials have many important roles in animal locomotion and object manipulation. In robotic applications soft materials can store and release energy, absorb impacts, increase compliance and increase the range of possible shape profiles using minimal actuators. The shape changing ability is also a potential tool to manipulate friction forces caused by contact with the environment. These advantages are accompanied by challenges of soft material actuation and the need to exploit frictional interactions to generate locomotion. Accordingly, the design of soft robots involves exploitation of continuum properties of soft materials for manipulating frictional interactions that result in robot locomotion. The research presents design and control of a soft body robot that uses its shape change capability for locomotion. The bioinspired (caterpillar) modular robot design is a soft monolithic body which interacts with the environment at discrete contact points (caterpillar prolegs). The deformable body is actuated by muscle-like shape memory alloy coils and the discrete contact points manipulate friction in a binary manner. This novel virtual grip mechanism combines two materials with different coefficients of frictions (sticky-slippery) to control the robot-environment friction interactions. The research also introduces a novel control concept that discretizes the robot-environment-friction interaction into binary states. This facilitates formulation of a control framework that is independent of the specific actuator or soft material properties and can be applied to multi-limbed soft robots. The transitions between individual robot states are assigned a reward that allow optimized state transition control sequences to be calculated. This conceptual framework is extremely versatile and we show how it can be applied to situations in which the robot loses limb function.


title: Investigating transfer of learning in advanced quantum mechanics

abstract: Research suggests that students often have difficulty transferring their learning from one context to another. We examine upper-level undergraduate and graduate students' facility with questions about the interference pattern in the double-slit experiment (DSE) with single photons and polarizers of various orientations placed in front of one or both slits. Before answering these types of questions, students had worked through a tutorial on the Mach-Zehnder Interferometer (MZI) in which they learned about interference of single photons when polarizers of various orientations are placed in the two paths of the MZI. After working on the MZI tutorial, students were asked similar questions in the isomorphic context of the DSE. We discuss the extent to which they were able to transfer what they learned in the context of the MZI to analogous problems in the isomorphic context of the DSE.


title: Probing the Kosterlitz-Thouless transition in 1D Heisenberg antiferromagnet based on topological properties of its ground state

abstract: A Kosterlitz-Thouless phase transition in the ground state of an antiferromagnetic spin-$\frac{1}{2}$ Heisenberg chain with nearest and next-nearest-neighbor interactions is re-investigated from a different perspective: An unequivocal correspondence is found between components of the scalar product $\langle \Psi_0| \Psi_0 \rangle$ and geometrical objects. One can classify these objects according to whether any two of them can be transformed into each other in a continuous way (belong to the same homotopy class). A finite size scaling of the "connection term`` $\langle \Psi_0|\partial \Psi_0 \rangle$ with respect to chain length (16, 18, 20, 22, 24 spins) for each homotopy class of above mentioned objects leads to the critical value of $\lambda$ with rather high accuracy.


title: Strain-controlled switch between ferromagnetism and antiferromagnetism in 1T-CrX2 (X = Se, Te) monolayers

abstract: We report on the strain-induced switch between ferromagnetic (FM) and antiferromagnetic (AFM) orderings in 1T-CrX2 (X = Se, Te) monolayers based on the first-principles calculations. The CrSe2 and CrTe2 monolayers without strains are found to be AFM and FM, respectively. Under the biaxial tensile strain, the CrSe2 monolayer tends to be FM when the strain is larger than 2%. The FM state is further stabilized when the strain is increased. Moreover, the CrSe2 monolayer changes to be half-metallic when the tensile strain is larger than 10%. While for the CrTe2 monolayer, the critical strain at which the transition between the FM and AFM states occurs is compressive, of -1%. Relatively small tensile strains of 4% and 2%, respectively, can enhance the Curie temperature of CrSe2 and CrTe2 monolayers above the room temperature. The strain-induced switch between the FM and AFM states in CrSe2 (CrTe2) monolayer can be understood by the competition between the AFM Cr-Cr direct exchange and FM Cr-Se(Te)-Cr superexchange interactions. The tunable and attractive magnetic and electronic properties controlled by the flexible strain are desirable for the future nanoelectronic applications.


title: Regular irreducible characters of a hyperspecial compact group and Weil representations over finite fields

abstract: A method to construct irreducible unitary representations of a hyperspecial compact subgroup of a reductive group over p-adic field with odd p is presented. Our method is based upon Cliffods theory and Weil representations over finite fields. It works under an assumption of the triviality of certain Schur multipliers defined for an algebraic group over a finite field. The assumption of the triviality has good evidences in the case of general linear groups and highly probable in regular cases in general. We will give several examples of classical groups where the Schur multipliers are actually trivial.


title: Measure preserving actions of affine semigroups and {x+y,xy} patterns

abstract: Ergodic and combinatorial results obtained in [10] involved measure preserving actions of the affine group ${\mathcal A}_K$ of a countable field $K$. In this paper we develop a new approach based on ultrafilter limits which allows one to refine and extend the results obtained in [10] to a more general situation involving the measure preserving actions of the non-amenable affine semigroups of a large class of integral domains. (The results in [10] heavily depend on the amenability of the affine group of a field). Among other things, we obtain, as a corollary of an ultrafilter ergodic theorem, the following result: Let $K$ be a number field and let ${\mathcal O}_K$ be the ring of integers of $K$. For any finite partition $K=C_1\cup\cdots\cup C_r$ there exists $i\in\{1,\dots,r\}$ and many $x\in K$ and $y\in{\mathcal O}_K$ such that $\{x+y,xy\}\subset C_i$.


title: Extended use of superconducting magnets for bio-medical development

abstract: Magnetic fields interact with biological cells affecting them in variety of ways which are usually hard to predict. Among them, it was observed that strong fields can align dividing cells in a preferred direction. It was also demonstrated that dividing cancer cells are effectively destroyed by applying electric fields in vivo with a success rate dependent on the cell-to-field orientation. Based on these facts, the present note aims to suggest the use of magnetic and electric fields for improved cancer treatment. Several possibilities of generating the electric fields inside the magnetic field volume are reviewed, main tentative approaches are described and discussed. Most if not all of them require special magnet configuration research which can be based on existing magnet systems in operation or in development.


title: Fused RSOS Lattice Models as Higher-Level Nonunitary Minimal Cosets

abstract: We consider the Forrester-Baxter RSOS lattice models with crossing parameter $\lambda=(m'\!-\!m)\pi/m'$ in Regime~III. In the continuum scaling limit, these models are described by the minimal models ${\cal M}(m,m')$. We conjecture that, for $\lambda<\pi/n$, the $n\times n$ fused RSOS models with $n\ge 2$ are described by the higher-level coset $(A^{(1)}_1)_k\otimes (A^{(1)}_1)_n/(A^{(1)}_1)_{k+n}$ at fractional level $k=nM/(M'\!-\!M)-2$ with $(M,M')=\big(nm-(n\!-\!1)m',m'\big)$. To support this conjecture, we investigate the one-dimensional sums arising from Baxter's off-critical corner transfer matrices. In unitary cases ($m=m'\!-\!1$) it is known that, up to leading powers of $q$, these coincide with the branching functions $b_{r,s,\ell}^{m'\!-n,m'\!,n}(q)$. For general nonunitary cases ($m<m'\!-\!1$), we identify the ground state one-dimensional RSOS paths and relate them to the quantum numbers $(r,s,\ell)$ in the various sectors. For $n=1,2,3$, we obtain the local energy functions $H(a,b,c)$ in a suitable gauge and verify that the associated one-dimensional sums produce finitized forms that converge, as $N$ becomes large, to the fractional level branching functions $b_{r,s,\ell}^{M,M'\!,n}(q)$. Extending the work of Schilling, we also conjecture finitized bosonic branching functions $b_{r,s,\ell}^{M,M'\!,n;(N)}(q)$ for general $n$ and check that these agree with the one-dimensional sums for $n=1,2,3$ out to system sizes $N=14$. Lastly, the finitized Kac characters $\chi_{r,s,\ell}^{P,P'\!,n;(N)}(q)$ of the $n\times n$ fused logarithmic minimal models ${\cal LM}(p,p')$ are obtained by taking the {\em logarithmic limit\/} $m,m'\to\infty$ with $m/m'\to p/p'+$.


title: A Review of Feature Selection Methods Based on Mutual Information

abstract: In this work we present a review of the state of the art of information theoretic feature selection methods. The concepts of feature relevance, redundance and complementarity (synergy) are clearly defined, as well as Markov blanket. The problem of optimal feature selection is defined. A unifying theoretical framework is described, which can retrofit successful heuristic criteria, indicating the approximations made by each method. A number of open problems in the field are presented.


title: A Social Network Framework to Explore Healthcare Collaboration

abstract: A patient-centric approach to healthcare leads to an informal social network among medical professionals. This chapter presents a research framework to: identify the collaboration structure among physicians that is effective and efficient for patients, discover effective structural attributes of a collaboration network that evolves during the course of providing care, and explore the impact of socio-demographic characteristics of healthcare professionals, patients, and hospitals on collaboration structures, from the point of view of measurable outcomes such as cost and quality of care. The framework uses illustrative examples drawn from a data set of patients undergoing hip replacement surgery.


title: Symplectic Rigidity of Real Bidisc

abstract: Let $\mathbb{D}$ be the unit disc in $\mathbb{C}$, then $\mathbb{D}^n(r)$ is the complex or symplectic $n$-discs of radius $r$. Let $z_j = x_j+iy_j\in\mathbb{C}, j=1,2$ and $\mathbb{D}_{\mathbb{R}}^2=\{(z_1,z_2) : |x_1|^2+|x_2|^2<1,|y_1|^2+|y_2|^2<1\}$ be the real bidisc. In this paper we will prove the following two theorems: 1) If $T\in O(4)$ is an orthogonal transformation on $\mathbb{R}^4$, then $T(\mathbb{D}^2)$ is symplectomorphic to $\mathbb{D}^2$ w.r.t. the standard symplectic form on $\mathbb{R}^4$ if and only if $T$ is unitary or conjugate to unitary. 2) For $r\geq 1$ and $n\geq 2$, $\mathbb{D}_{\mathbb{R}}^2\times \mathbb{D}^{n-2}(r)$ and $\mathbb{D}^2\times \mathbb{D}^{n-2}(r)$ are not symplectomorphic w.r.t. the standard symplectic form on $\mathbb{C}^n$.


title: Surface tension and Laplace pressure in triangulated surface models for membranes without fixed boundary

abstract: A Monte Carlo (MC) study is performed to evaluate the surface tension $\gamma $ of spherical membranes that may be regarded as the models of the lipid layers. We use the canonical surface model defined on the self-avoiding triangulated lattices. The surface tension $\gamma $ is calculated by keeping the total surface area $A$ constant during the MC simulations. In the evaluation of $\gamma $, we use $A$ instead of the projected area $A_p$, which is unknown due to the fluctuation of the spherical surface without boundary. The pressure difference ${\it\Delta}p $ between the inner and the outer sides of the surface is also calculated by maintaining the enclosed volume constant. Using ${\it\Delta}p $ and the Laplace formula, we obtain the tension, which is considered to be equal to the frame tension $\tau$ conjugate to $A_p$, and check whether or not $\gamma $ is consistent with $\tau$. We find reasonable consistency between $\gamma$ and $\tau$ in the region of sufficiently large bending rigidity $\kappa$ or sufficiently large $A/N$. It is also found that $\tau$ becomes constant in the limit of $A/N\to \infty$ both in the tethered and fluid surfaces.


title: Pure states on Cuntz algebras arising from geometric progressions

abstract: Let ${\cal O}_n$ denote the Cuntz algebra for $n\geq 2$. We introduce an embedding $f$ of ${\cal O}_m$ into ${\cal O}_n$ arising from a geometric progression of Cuntz generaters of ${\cal O}_n$. By identifying ${\cal O}_m$ with $f({\cal O}_m)$, we extend Cuntz states on ${\cal O}_m$ to ${\cal O}_n$. We show (i) a necessary and sufficient condition of the uniqueness of the extension, (ii) the complete classification of all such extensions up to unitary equivalence of their GNS representations, and (iii) the decomposition formula of a mixing state into a convex hull of pure states. The complete set of invariants of all GNS representations by such pure states is given as a certain set of complex unit vectors.


title: Constructing Abstraction Hierarchies Using a Skill-Symbol Loop

abstract: We describe a framework for building abstraction hierarchies whereby an agent alternates skill- and representation-acquisition phases to construct a sequence of increasingly abstract Markov decision processes. Our formulation builds on recent results showing that the appropriate abstract representation of a problem is specified by the agent's skills. We describe how such a hierarchy can be used for fast planning, and illustrate the construction of an appropriate hierarchy for the Taxi domain.


title: mplot: An R Package for Graphical Model Stability and Variable Selection Procedures

abstract: The mplot package provides an easy to use implementation of model stability and variable inclusion plots (M\"uller and Welsh 2010; Murray, Heritier, and M\"uller 2013) as well as the adaptive fence (Jiang, Rao, Gu, and Nguyen 2008; Jiang, Nguyen, and Rao 2009) for linear and generalised linear models. We provide a number of innovations on the standard procedures and address many practical implementation issues including the addition of redundant variables, interactive visualisations and approximating logistic models with linear models. An option is provided that combines our bootstrap approach with glmnet for higher dimensional models. The plots and graphical user interface leverage state of the art web technologies to facilitate interaction with the results. The speed of implementation comes from the leaps package and cross-platform multicore support.


title: Brouwer's fixed-point theorem in real-cohesive homotopy type theory

abstract: We combine Homotopy Type Theory with axiomatic cohesion, expressing the latter internally with a version of "adjoint logic" in which the discretization and codiscretization modalities are characterized using a judgmental formalism of "crisp variables". This yields type theories that we call "spatial" and "cohesive", in which the types can be viewed as having independent topological and homotopical structure. These type theories can then be used to study formally the process by which topology gives rise to homotopy theory (the "fundamental $\infty$-groupoid" or "shape"), disentangling the "identifications" of Homotopy Type Theory from the "continuous paths" of topology. In a further refinement called "real-cohesion", the shape is determined by continuous maps from the real numbers, as in classical algebraic topology. This enables us to reproduce formally some of the classical applications of homotopy theory to topology. As an example, we prove Brouwer's fixed-point theorem.


title: Probing higher-order primordial non-Gaussianity with galaxy surveys

abstract: With a radio continuum galaxy survey by the Square Kilometre Array (SKA), a photometric galaxy survey by Euclid and their combination, we forecast future constraints on primordial non-Gaussianity. We focus on the potential impact of local-type higher-order nonlinear parameters on the parameter estimation and particularly the confirmation of the inflationary consistency inequality. Nonstandard inflationary models, such as multifield models, introduce the scale-dependent stochastic clustering of galaxies on large scales, which is a unique probe of mechanism for generating primordial density fluctuations. Our Fisher matrix analysis indicates that a deep and wide survey provided by SKA is more advantageous to constrain $\tau_{\rm NL}$, while Euclid has a strong constraining power for $f_{\rm NL}$ due to the redshift information, suggesting that the joint analysis between them is quite essential to break the degeneracy between the nonlinear parameters. The combination of the full SKA and Euclid will achieve the precision level needed to confirm the consistency inequality even for $f_{\rm NL}\approx 1.5$ and $\tau_{\rm NL}\approx 17$, though it is still hard for a single survey to confirm it when $f_{\rm NL}\lesssim 2.7$.


title: An extending result on spectral radius of bipartite graphs

abstract: Let $G$ denote a bipartite graph with $e$ edges without isolated vertices. It was known that the spectral radius of $G$ is at most the square root of $e$, and the upper bound is attained if and only if $G$ is a complete bipartite graph. Suppose that $G$ is not a complete bipartite graph, and $e-1$ and $e+1$ are not twin primes. We determine the maximal spectral radius of $G$. As a byproduct of our study, we obtain a spectral characterization of a pair $(e-1, e+1)$ of integers to be a pair of twin primes.


title: The S-matrix algebra of the AdS2 X S2 superstring

abstract: In this paper we find the Yangian algebra responsible for the integrability of the AdS2 X S2 X T^6 superstring in the planar limit. We demonstrate the symmetry of the corresponding exact S-matrix in the massive sector, including the presence of the secret symmetry. We give two alternative presentations of the Hopf algebra, along with related discussions on the issue of evaluation representations. We study the classical r-matrix, and re-discover the need for a secret symmetry also in this context. Finally, taking the simplifying zero-coupling limit of the S-matrix as a generating R-matrix for the Algebraic Bethe Ansatz, we obtain an effective model of free fermions on a periodic spin-chain. This limit should provide hints to the one-loop anomalous dimension of the mysterious superconformal quantum mechanics dual to the superstring theory in this geometry.


title: Fractional coverings, greedy coverings, and rectifier networks

abstract: A rectifier network is a directed acyclic graph with distinguished sources and sinks; it is said to compute a Boolean matrix $M$ that has a $1$ in the entry $(i,j)$ iff there is a path from the $j$th source to the $i$th sink. The smallest number of edges in a rectifier network computing $M$ is a classic complexity measure on matrices, which has been studied for more than half a century. We explore two well-known techniques that have hitherto found little to no applications in this theory. Both of them build on a basic fact that depth-$2$ rectifier networks are essentially weighted coverings of Boolean matrices with rectangles. We obtain new results by using fractional and greedy coverings (defined in the standard way). First, we show that all fractional coverings of the so-called full triangular matrix have cost at least $n\log n$. This provides (a fortiori) a new proof of the tight lower bound on its depth-$2$ complexity (the exact value has been known since 1965, but previous proofs are based on different arguments). Second, we show that the greedy heuristic is instrumental in tightening the upper bound on the depth-$2$ complexity of the Kneser-Sierpi\'nski (disjointness) matrix. The previous upper bound is $O(n^{1.28})$, and we improve it to $O(n^{1.17})$, while the best known lower bound is $\Omega(n^{1.16})$. Third, using fractional coverings, we obtain a form of direct product theorem that gives a lower bound on unbounded-depth complexity of Kronecker (tensor) products of matrices. In this case, the greedy heuristic shows (by an argument due to Lov\'asz) that our result is only a logarithmic factor away from the "full" direct product theorem. Our second and third results constitute progress on open problem 7.3 and resolve, up to a logarithmic factor, open problem 7.5 from a recent book by Jukna and Sergeev (in Foundations and Trends in Theoretical Computer Science (2013)).


title: 3-state Hamiltonians associated to solvable 33-vertex models

abstract: Using the nested coordinate Bethe ansatz, we study 33-vertex models, where only one global charge with degenerate eigenvalues exists and each site possesses three internal degrees of freedom. In the context of Markovian processes, they correspond to diffusing particles with two possible internal states which may be exchanged during the diffusion (transmutation). The first step of the nested coordinate Bethe ansatz is performed providing the eigenvalues in terms of rapidities. We give the constraints ensuring the consistency of the computations. These rapidities also satisfy Bethe equations involving $4\times 4$ R-matrices, solutions of the Yang--Baxter equation which implies new constraints on the models. We solve them allowing us to list all the solvable 33-vertex models.


title: Fractality in Persistence Decay and Domain Growth during Ferromagnetic Ordering: Dependence upon initial correlation

abstract: Dynamics of ordering in Ising model, following quench to zero temperature, have been studied via Glauber spin-flip Monte Carlo simulations in space dimensions $d=2$ and $3$. One of the primary objectives has been to understand phenomena associated with the persistent spins, viz., time decay in the number of unaffected spins, growth of the corresponding pattern and its fractal dimensionality, for varying correlation length in the initial configurations, prepared at different temperatures, at and above the critical value. It is observed that the fractal dimensionality and the exponent describing the power-law decay of persistence probability are strongly dependent upon the relative values of nonequilibrium domain size and the initial equilibrium correlation length. Via appropriate scaling analyses, these quantities have been estimated for quenches from infinite and critical temperatures. The above mentioned dependence is observed to be less pronounced in higher dimension. In addition to these findings for the local persistence, we present results for the global persistence as well. Further, important observations on the standard domain growth problem are reported. For the latter, a controversy in $d=3$, related to the value of the exponent for the power-law growth of the average domain size with time, has been resolved.


title: Ramanujan-type Congruences for $\ell$-Regular Partitions Modulo $3, 5, 11$ and $13$

abstract: Let $b_\ell(n)$ be the number of $\ell$-regular partitions of $n$. Recently, Hou et al established several infinite families of congruences for $b_\ell(n)$ modulo $m$, where $(\ell,m)=(3,3),(6,3),(5,5),(10,5)$ and $(7,7)$. In this paper, by the vanishing property given by Hou et al, we show an infinite family of congruence for $b_{11}(n)$ modulo $11$. Moreover, for $\ell= 3, 13$ and $25$, we obtain three infinite families of congruences for $b_{\ell}(n)$ modulo $3, 5$ and $13$ by the theory of Hecke eigenforms.