title: Introduction of Curvilinear Coordinates into Numerical Analysis

abstract: Introduction of curvilinear coordinates might be very convenient in many cases. Theoretically, tensor analysis would be most suited. However, tensor notation can't be used in numerical procedure. For example, the strict discrimination of upper and lower suffices is impossible in the present computer languages. In the present paper, a rather elementary approach more suited to write codes of computer programming is adapted. If we introduce curvilinear coordinates, we would be able to treat in more natural way to handle fixed discontinuity, moving discontinuity and curved boundary. We could generate fine mesh in the neighborhood of fixed and moving discontinuities. Furthermore, if we introduce curvilinear coordinates, we can transform a non-square region into square one and can use a regular mesh. Usually in numerical calculation, a curved boundary is approximated by a jagged or non-smooth boundary. An introduction of curvilinear coordinates could solve many problems in the numerical analyses.


title: Two transitional type~Ia supernovae located in the Fornax cluster member NGC 1404: SN 2007on and SN 2011iv

abstract: We present an analysis of ultraviolet (UV) to near-infrared observations of the fast-declining Type Ia supernovae (SNe Ia) 2007on and 2011iv, hosted by the Fornax cluster member NGC 1404. The B-band light curves of SN 2007on and SN 2011iv are characterised by dm_15(B) decline-rate values of 1.96 mag and 1.77 mag, respectively. Although they have similar decline rates, their peak B- and H-band magnitudes differ by ~0.60 mag and ~0.35 mag, respectively. After correcting for the luminosity vs. decline rate and the luminosity vs. colour relations, the peak B-band and H-band light curves provide distances that differ by ~14% and ~9%, respectively. These findings serve as a cautionary tale for the use of transitional SNe Ia located in early-type hosts in the quest to measure cosmological parameters. Interestingly, even though SN 2011iv is brighter and bluer at early times, by three weeks past maximum and extending over several months, its B-V colour is 0.12 mag redder than that of SN 2007on. To reconcile this unusual behaviour, we turn to guidance from a suite of spherical one-dimensional Chandrasekhar-mass delayed-detonation explosion models. In this context, 56Ni production depends on both the so-called transition density and the central density of the progenitor white dwarf. To first order, the transition density drives the luminosity-width relation, while the central density is an important second-order parameter. Within this context, the differences in the B-V color evolution along the Lira regime suggests the progenitor of SN~2011iv had a higher central density than SN~2007on.


title: Gasoline2: A Modern SPH Code

abstract: The methods in the Gasoline2 Smoothed Particle Hydrodynamics (SPH) code are described and tested. Gasoline2 is the most recent version of the Gasoline code for parallel hydrodynamics and gravity with identical hydrodynamics to the Changa code. As with other Modern SPH codes, we prevent sharp jumps in time steps, use upgraded kernels and larger neighbour numbers and employ local viscosity limiters. Unique features in Gasoline2 include its Geometric-Density-Average Force expression, explicit Turbulent Diffusion terms and Gradient-Based shock detection to limit artificial viscosity. This last feature allows Gasoline2 to completely avoid artificial viscosity in non-shocking compressive flows. We present a suite of tests demonstrating the value of these features with the same code configuration and parameter choices used for production simulations.


title: On Universality of Holographic Results for (2+1)-Dimensional CFTs on Curved Spacetimes

abstract: The behavior of holographic CFTs is constrained by the existence of a bulk dual geometry. For example, in (2+1)-dimensional holographic CFTs living on a static spacetime with compact spatial slices, the vacuum energy must be nonpositive, certain averaged energy densities must be nonpositive, and the spectrum of scalar operators is bounded from below by the Ricci scalar of the CFT geometry. Are these results special to holographic CFTs? Here we show that for perturbations about appropriate backgrounds, they are in fact universal to all CFTs, as they follow from the universal behavior of two- and three-point correlators of known operators. In the case of vacuum energy, we extend away from the perturbative regime and make global statements about its negativity properties on the space of spatial geometries. Finally, we comment on the implications for dynamics which are dissipative and driven by such a vacuum energy and we remark on similar results for the behavior of the Euclidean partition function on deformations of flat space or the round sphere.


title: Active Galactic Nuclei outflows in galaxy discs

abstract: Galactic outflows, driven by active galactic nuclei (AGN), play a crucial role in galaxy formation and in the self-regulated growth of supermassive black holes (BHs). AGN feedback couples to and affects gas, rather than stars, and in many, if not most, gas-rich galaxies cold gas is rotationally supported and settles in a disc. We present a 2D analytical model for AGN-driven outflows in a gaseous disc and demonstrate the main improvements, compared to existing 1D solutions. We find significant differences for the outflow dynamics and wind efficiency. The outflow is energy-driven due to inefficient cooling up to a certain AGN luminosity ($\sim 10^{43}$erg/s in our fiducial model), above which the outflow remains momentum-driven in the disc up to galactic scales. We reproduce results of 3D simulations that gas is preferentially ejected perpendicular to the disc and find that the fraction of ejected interstellar medium is lower than in 1D models. The recovery time of gas in the disc, defined as the freefall time from the radius to which the AGN pushes the ISM at most, is remarkably short, of the order 1Myr. This indicates that AGN-driven winds cannot suppress BH growth for long. Without the inclusion of supernova feedback, we find a scaling of the black hole mass with the halo velocity dispersion of $M_\mathrm{BH} \propto \sigma ^{4.8}$.


title: ALMA observations of dust polarization and molecular line emission from the Class 0 protostellar source Serpens SMM1

abstract: We present high angular resolution dust polarization and molecular line observations carried out with the Atacama Large Millimeter/submillimeter Array (ALMA) toward the Class 0 protostar Serpens SMM1. By complementing these observations with new polarization observations from the Submillimeter Array (SMA) and archival data from the Combined Array for Research in Millimeter-wave Astronomy (CARMA) and the James Clerk Maxwell Telescopes (JCMT), we can compare the magnetic field orientations at different spatial scales. We find major changes in the magnetic field orientation between large (~0.1 pc) scales -- where the magnetic field is oriented E-W, perpendicular to the major axis of the dusty filament where SMM1 is embedded -- and the intermediate and small scales probed by CARMA (~1000 AU resolution), the SMA (~350 AU resolution), and ALMA (~140 AU resolution). The ALMA maps reveal that the redshifted lobe of the bipolar outflow is shaping the magnetic field in SMM1 on the southeast side of the source; however, on the northwestern side and elsewhere in the source, low velocity shocks may be causing the observed chaotic magnetic field pattern. High-spatial-resolution continuum and spectral-line observations also reveal a tight (~130 AU) protobinary system in SMM1-b, the eastern component of which is launching an extremely high-velocity, one-sided jet visible in both CO(2-1) and SiO(5-4); however, that jet does not appear to be shaping the magnetic field. These observations show that with the sensitivity and resolution of ALMA, we can now begin to understand the role that feedback (e.g., from protostellar outflows) plays in shaping the magnetic field in very young, star-forming sources like SMM1.


title: On the Progenitor of the Type Ic Supernova 2012fh

abstract: Little is observationally known about the progenitors of Type Ibc supernovae (SNe) or the typical activity of SNe progenitors in their final years. Here, we analyze deep Large Binocular Telescope imaging data spanning the 4 years before and after the Type Ic SN 2012fh using difference imaging. We place 1$\sigma$ upper limits on the detection of the progenitor star at $M_R>-4.0$, $M_V>-3.8$, $M_B>-3.1$, and $M_U>-3.8$ mag. These limits are the tightest placed on a Type Ic SNe and they largely rule out single star evolutionary models in favor of a binary channel as the origin of this SN. We also constrain the activity of the progenitor to be small on an absolute scale, with the RMS $UBVR$ optical variability $<2500L_\odot$ and long-term dimming or brightening trends $<1000L_\odot/\text{year}$ in all four bands.


title: Collapsed Dark Matter Structures

abstract: The distributions of dark matter and baryons in the Universe are known to be very different: the dark matter resides in extended halos, while a significant fraction of the baryons have radiated away much of their initial energy and fallen deep into the potential wells. This difference in morphology leads to the widely held conclusion that dark matter cannot cool and collapse on any scale. We revisit this assumption, and show that a simple model where dark matter is charged under a "dark electromagnetism" can allow dark matter to form gravitationally collapsed objects with characteristic mass scales much smaller than that of a Milky Way-type galaxy. Though the majority of the dark matter in spiral galaxies would remain in the halo, such a model opens the possibility that galaxies and their associated dark matter play host to a significant number of collapsed substructures. The observational signatures of such structures are not well explored, but potentially interesting.


title: Exploring extra dimensions through inflationary tensor modes

abstract: Predictions of inflationary schemes can be influenced by the presence of extra dimensions. This could be of particular relevance for the spectrum of gravitational waves in models where the extra dimensions provide a brane-world solution to the hierarchy problem. Apart from models of large as well as exponentially warped extra dimensions, we analyze the size of tensor modes in the Linear Dilaton scheme recently revived in the discussion of the "clockwork mechanism". The results are model dependent, significantly enhanced tensor modes on one side and a suppression on the other. In some cases we are led to a scheme of "remote inflation", where the expansion is driven by energies at a hidden brane. In all cases where tensor modes are enhanced, the requirement of perturbativity of gravity leads to a stringent upper limit on the allowed Hubble rate during inflation.


title: Superconductivity in three-dimensional spin-orbit coupled semimetals

abstract: Motivated by the experimental detection of superconductivity in the low-carrier density half-Heusler compound YPtBi, we study the pairing instabilities of three-dimensional strongly spin-orbit coupled semimetals with a quadratic band touching point. In these semimetals the electronic structure at the Fermi energy is described by spin j=3/2 quasiparticles, which are fundamentally different from those in ordinary metals with spin j=1/2. We develop a general approach to analyzing pairing instabilities in j=3/2 materials by decomposing the pair scattering interaction into irreducible channels, projecting them to the Fermi surface and deriving the corresponding Eliashberg theory. Applying our method to a generic density-density interaction in YPtBi we establish the following results: (i) The pairing strength in each channel uniquely encodes the j=3/2 nature of the Fermi surface band structure--a manifestation of the fundamental difference with ordinary metals. In particular, this implies that Anderson's theorem, which addresses the effect of spin-orbit coupling and disorder on pairing states of spin-1/2 electrons, cannot be applied in this case. (ii) The leading pairing instabilities are different for electron and hole doping. This implies that superconductivity depends on carrier type. (iii) In the case of hole doping--relevant to YPtBi, we find two odd-parity channels in close competition with s-wave pairing. One of these two channels is a multicomponent pairing channel, allowing for the possibility of time-reversal symmetry breaking. (iv) In the case of Coulomb interactions mediated by the long-ranged electric polarization of optical phonon modes, a significant coupling strength is generated in spite of the extremely low density of carriers. Furthermore, non-linear response and Fermi liquid corrections can favor non-s-wave pairing and potentially account for the experimentally-observed Tc.


title: Black Holes on FIRE: Stellar Feedback Limits Early Feeding of Galactic Nuclei

abstract: We introduce massive black holes (BHs) in the Feedback In Realistic Environments project and perform high-resolution cosmological hydrodynamic simulations of quasar-mass halos ($M_{\rm halo}(z=2) \approx 10^{12.5}\,\rm{M}_{\odot}$) down to $z=1$. These simulations model stellar feedback by supernovae, stellar winds, and radiation, and BH growth using a gravitational torque-based prescription tied to resolved properties of galactic nuclei. We do not include BH feedback. We show that early BH growth occurs through short ($\lesssim 1\,$Myr) accretion episodes that can reach or even exceed the Eddington rate. In this regime, BH growth is limited by bursty stellar feedback continuously evacuating gas from galactic nuclei, and BHs remain under-massive relative to the local $M_{\rm BH}$-$M_{\rm bulge}$ relation. BH growth is more efficient at later times, when the nuclear stellar potential retains a significant gas reservoir, star formation becomes less bursty, and galaxies settle into a more ordered state, with BHs rapidly converging onto the scaling relation when the host reaches $M_{\rm bulge} \sim 10^{10}\,\rm{M}_{\odot}$. Our results are not sensitive to the details of the accretion model so long as BH growth is tied to the gas content within $\sim 100\,$pc of the BH. Our simulations imply that bursty stellar feedback has strong implications for BH and AGN demographics, especially in the early Universe and for low-mass galaxies.


title: The first all-sky view of the Milky Way stellar halo with Gaia+2MASS RR Lyrae

abstract: We exploit the first \gaia data release to study the properties of the Galactic stellar halo as traced by RR Lyrae. We demonstrate that it is possible to select a pure sample of RR Lyrae using only photometric information available in the Gaia+2MASS catalogue. The final sample contains about 21600 RR Lyrae covering an unprecedented fraction ($\sim60\%$) of the volume of the Galactic inner halo ($\text{R}<28$ kpc). We study the morphology of the stellar halo by analysing the RR Lyrae distribution with parametric and non-parametric techniques. Taking advantage of the uniform all-sky coverage, we test halo models more sophisticated than usually considered in the literature, such as those with varying flattening, tilt and/or offset of the halo with respect to the Galactic disc. A consistent picture emerges: the inner halo is well reproduced by a smooth distribution of stars settled on triaxial ellipsoids. The minor axis is perpendicular to the Milky Ways disc, while the major axis is misaligned by $\sim20^{\circ}$ from the Galactic Y axis. The elongation along the major axis is mild ($\text{p}=1.27$), and the vertical flattening is shown to evolve from a squashed state with $\text{q}\approx0.57$ in the centre to a more spherical $\text{q}\approx0.75$ at the outer edge of our dataset. The density slope is well approximated by a single power-law with exponent $\alpha=-2.96$. Within the range probed, we see no significant evidence for a change of the radial density slope, out of the plane tilt or an offset of the halo with respect to the Galaxy's centre.


title: Screened moments in a Kondo insulator

abstract: Long known to have thermodynamic properties at odds with its insulating electrical transport, SmB6 has been the subject of great debate as it is unclear whether its unusual properties are related to the bulk or novel metallic surface states. We have observed a bulk moment-screening effect in nominally pure and Gd-doped SmB6 via heat capacity, magnetization, and resistivity measurements, and show this new Kondo-impurity like effect provides an unexpected but intuitive explanation for metal-like phenomena stemming from the strongly interacting host system. This affords a coherent understanding for decades of mysteries in strongly-correlated insulators, reveals the expanded utility of techniques previously only utilized for metals, and presents the novel effect of even highly-dilute impurities in strongly correlated insulators.


title: Impact of Resonance on Thermal Targets for Invisible Dark Photon Searches

abstract: Dark photons in the MeV to GeV mass range are important targets for experimental searches. We consider the case where dark photons $A'$ decay invisibly to hidden dark matter $X$ through $A' \to XX$. For generic masses, proposed accelerator searches are projected to probe the thermal target region of parameter space, where the $X$ particles annihilate through $XX \to A' \to \text{SM}$ in the early universe and freeze out with the correct relic density. However, if $m_{A'} \approx 2m_X$, dark matter annihilation is resonantly enhanced, shifting the thermal target region to weaker couplings. For $\sim 10\%$ degeneracies, we find that the annihilation cross section is generically enhanced by four (two) orders of magnitude for scalar (pseudo-Dirac) dark matter. For such moderate degeneracies, the thermal target region drops to weak couplings beyond the reach of all proposed accelerator experiments in the scalar case and becomes extremely challenging in the pseudo-Dirac case. Proposed direct detection experiments can probe moderate degeneracies in the scalar case. For greater degeneracies, the effect of the resonance can be even more significant, and both scalar and pseudo-Dirac cases are beyond the reach of all proposed accelerator and direct detection experiments. For scalar dark matter, we find an absolute minimum that sets the ultimate experimental sensitivity required to probe the entire thermal target parameter space, but for pseudo-Dirac fermions, we find no such thermal target floor.


title: Thermal transport in a two-dimensional $\mathbb{Z}_2$ spin liquid

abstract: We study the dynamical thermal conductivity of the two-dimensional Kitaev spin-model on the honeycomb lattice. We find a strongly temperature dependent low-frequency spectral intensity as a direct consequence of fractionalization of spins into mobile Majorana matter and a static $\mathbb{Z}_{2}$ gauge field. The latter acts as an emergent thermally activated disorder, leading to the appearance of a pseudogap which partially closes in the thermodynamic limit, indicating a dissipative heat conductor. Our analysis is based on complementary calculations of the current correlation function, comprising exact diagonalization by means of a complete summation over all gauge sectors, as well as a phenomenological mean-field treatment of thermal gauge fluctuations, valid at intermediate and high temperatures. The results will also be contrasted against the conductivity discarding gauge fluctuations.


title: Note on gauge and gravitational anomalies of discrete $Z_N$ symmetries

abstract: In this note, we discuss the consistency conditions which a discrete $Z_N$ symmetry should satisfy in order that it is not violated by gauge and gravitational instantons. As examples, we enlist all the $Z_N$ ${\cal R}$-symmetries as well as non-${\cal R}$ $Z_N$ symmetries (N=2,3,4) in the minimally supersymmetric standard model (MSSM) that are free from gauge and gravitational anomalies. We show that there exists non-anomalous discrete symmetries that forbid Baryon number violation up to dimension 6 level (in superspace). We also observe that there exists no non-anomalous $Z_3$ ${\cal R}$-symmetry in the MSSM. Furthermore, we point out that in a theory with one Majorana spin 3/2 gravitino, a large class of $Z_4$ ${\cal R}$-symmetries are violated in the presence of Eguchi-Hanson (EH) gravitational instanton. This is also in general true for higher $Z_N$ ${\cal R}$-symmetries. We also notice that in 4 dimensional ${\cal N}=1$ supergravity, the global $U(1)$ ${\cal R}$-symmetry is always violated by the EH instanton irrespective of the matter content of the theory.


title: Higher-Spin Witten Effect and Two-Dimensional Fracton Phases

abstract: We study the role of "$\theta$ terms" in the action for three-dimensional $U(1)$ symmetric tensor gauge theories, describing quantum phases of matter hosting gapless higher-spin gauge modes and gapped subdimensional particle excitations, such as fractons. In Maxwell theory, the $\theta$ term is a total derivative which has no effect on the gapless photon, but has two important, closely related consequences: attaching electric charge to magnetic monopoles (the Witten effect) and leading to a Chern-Simons theory on the boundary. We will find that a similar story holds in the higher-spin $U(1)$ gauge theories. These theories admit generalized $\theta$ terms which have no effect on the gapless gauge mode, but which bind together electric and magnetic charges (both of which are generally subdimensional) in specific combinations, in a higher-spin manifestation of the Witten effect. We derive the corresponding Witten quantization condition. We find that, as in Maxwell theory, imposing time-reversal invariance restricts $\theta$ to certain discrete values. We also find that these new $\theta$ terms imply a non-trivial boundary structure. The boundaries host fracton excitations coupled to a tensor $U(1)$ gauge field with a Chern-Simons-like action, in both chiral and non-chiral varieties. These boundary theories open a door to the study of $U(1)$ fracton phases described by tensor Chern-Simons theories, not only on boundaries of three-dimensional systems, but also in strictly two spatial dimensions. We explicitly work through three examples of bulk and boundary theories, the principles of which can be readily extended to arbitrary higher-spin theories.


title: Hard processes in small systems

abstract: These proceedings from the Quark Matter 2017 conference give an overview of the latest experimental results on hard processes in small collision systems at RHIC and the LHC, discuss their implications, and consider several prospects for future measurements.


title: Time Evolution of Complexity in Abelian Gauge Theories - And Playing Quantum Othello Game -

abstract: Quantum complexity is conjectured to probe inside of black hole horizons (or wormhole) via gauge gravity correspondence. In order to have a better understanding of this correspondence, we study time evolutions of complexities for generic Abelian pure gauge theories. For this purpose, we discretize $U(1)$ gauge group as $\mathbf{Z}_N$ and also continuum spacetime as lattice spacetime, and this enables us to define a universal gate set for these gauge theories, and evaluate time evolutions of the complexities explicitly. We find that for a generic class of diagonal Hamiltonians to achieve a large complexity $\sim \exp(\mbox{entropy})$, which is one of the conjectured criteria necessary to have a dual black hole, the Abelian gauge theory needs to be maximally nonlocal.


title: A tale of seven narrow spikes and a long trough: constraining the timing of the percolation of HII bubbles at the tail-end of reionization with ULAS J1120+0641

abstract: High-signal to noise observations of the Ly$\alpha$ forest transmissivity in the z = 7.085 QSO ULAS J1120+0641 show seven narrow transmission spikes followed by a long 240 cMpc/h trough. Here we use radiative transfer simulations of cosmic reionization previously calibrated to match a wider range of Ly$\alpha$ forest data to show that the occurrence of seven transmission spikes in the narrow redshift range z = 5.85 - 6.1 is very sensitive to the exact timing of reionization. Occurrence of the spikes requires the most under dense regions of the IGM to be already fully ionised. The rapid onset of a long trough at z = 6.12 requires a strong decrease of the photo-ionisation rate at z$\sim$6.1 in this line-of-sight, consistent with the end of percolation at this redshift. The narrow range of reionisation histories that we previously found to be consistent with a wider range of Ly$\alpha$ forest data have a reasonable probability of showing seven spikes and the mock absorption spectra provide an excellent match to the spikes and the trough in the observed spectrum of ULAS J1120+0641. Despite the large overall opacity of Ly$\alpha$ at z > 5.8, larger samples of high signal-to-noise observations of rare transmission spikes should therefore provide important further insights into the exact timing of the percolation of HII bubbles at the tail-end of reionization


title: Locating the intense interstellar scattering towards the inner Galaxy

abstract: We use VLBA+VLA observations to measure the sizes of the scatter-broadened images of 6 of the most heavily scattered known pulsars: 3 within the Galactic Centre (GC) and 3 elsewhere in the inner Galactic plane. By combining the measured sizes with temporal pulse broadening data from the literature and using the thin-screen approximation, we locate the scattering medium along the line of sight to these 6 pulsars. At least two scattering screens are needed to explain the observations of the GC sample. We show that the screen inferred by previous observations of SGR J1745-2900 and Sgr A*, which must be located far from the GC, falls off in strength on scales < 0.2 degree. A second scattering component closer to (< 2 kpc) or even (tentatively) within (< 700 pc) the GC produces most or all of the temporal broadening observed in the other GC pulsars. Outside the GC, the scattering locations for all three pulsars are ~2 kpc from Earth, consistent with the distance of the Carina-Sagittarius or Scutum spiral arm. For each object the 3D scattering origin coincides with a known HII region (and in one case also a supernova remnant), suggesting that such objects preferentially cause the intense interstellar scattering seen towards the Galactic plane. We show that the HII regions should contribute > 25% of the total dispersion measure (DM) towards these pulsars, and calculate reduced DM distances. Those distances for other pulsars lying behind HII regions may be similarly overestimated.


title: Hahn polynomials on polyhedra and quantum integrability

abstract: Orthogonal polynomials with respect to the hypergeometric distribution on lattices in polyhedral domains in ${\mathbb R}^d$, which include hexagons in ${\mathbb R}^2$ and truncated tetrahedrons in ${\mathbb R}^3$, are defined and studied. The polynomials are given explicitly in terms of the classical one-dimensional Hahn polynomials. They are also characterized as common eigenfunctions of a family of commuting partial difference operators. These operators provide symmetries for a system that can be regarded as a discrete extension of the generic quantum superintegrable system on the $d$-sphere. Moreover, the discrete system is proved to possess all essential properties of the continuous system. In particular, the symmetry operators for the discrete Hamiltonian define a representation of the Kohno-Drinfeld Lie algebra on the space of orthogonal polynomials, and an explicit set of $2d-1$ generators for the symmetry algebra is constructed. Furthermore, other discrete quantum superintegrable systems, which extend the quantum harmonic oscillator, are obtained by considering appropriate limits of the parameters.


title: First results from SAM-FP: Fabry-Perot observations with ground-layer adaptive optics - the structure and kinematics of the core of 30 Doradus

abstract: The aim of this paper is to present the first data set obtained with SOAR Adaptive Module-Fabry-Parot (SAM-FP), a Fabry-Perot instrument mounted inside the SOAR telescope Adaptive-Optics Module. This is the only existing imaging Fabry-Perot interferometer using laser-assisted ground-layer adaptive optics. SAM-FP was used to observe the ionized gas, traced by Halpha, in the centre of the 30 Doradus starburst (the Tarantula Nebula) in the Large Magellanic Cloud, with high spatial (~0.6" or 0.15 pc) and spectral (R=11200) resolution. Radial velocity, velocity dispersion and monochromatic maps were derived. The region displays a mix of narrow, sigma ~ 20 km/s profiles and multiple broader profiles with sigma ~ 70-80 km/s, indicating the complex nature of the nebula kinematics. A comparison with previously obtained VLT/FLAMES spectroscopy demonstrates that the data agree well in the regions of overlap, but the Fabry-Perot data are superior in spatial coverage. A preliminary analysis of the observations finds a new expanding bubble south of R136, with a projected radius of r=5.6 pc and an expansion velocity of 29 +/- 4 km/s. In addition, the first-time detailed kinematic maps derived here for several complexes and filaments of 30 Doradus allow identification of kinematically independent structures. These data exemplify the power of the combination of a high-order Fabry-Perot with a wide-field imager (3' x 3' GLAO-corrected field of view) for high-resolution spatial and spectral studies. In particular, SAM-FP data cubes are highly advantageous over multifibre or long-slit data sets for nebula structure studies and to search for small-scale bubbles, given their greatly improved spatial coverage. For reference, this paper also presents two appendices with detailed descriptions of the usage of Fabry-Perot devices, including formulae and explanations for understanding Fabry-Perot observations.


title: Limit linear series and the Amini-Baker construction

abstract: We draw comparisons between the author's recent construction of limit linear series for curves not of compact type and the Amini-Baker theory of limit linear series on metrized complexes, as well as the related theories of divisors on discrete graphs and on metric graphs. From these we conclude that the author's theory (like the others) satisfies the Riemann and Clifford inequalities. Motivated by our comparisons, we also develop negative results on Brill-Noether generality for certain families of metric graphs. Companion work of He develops our comparisons further and uses them to prove new results on smoothability of Amini-Baker limit linear series and of divisors on metric graphs.


title: Dyonic black holes at arbitrary locations

abstract: We construct and study stationary, asymptotically flat multicenter solutions describing regular black holes with non-Abelian hair (colored magnetic-monopole and dyon fields) in two models of $\mathcal{N}=2,d=4$ Super-Einstein-Yang-Mills theories: the quadratic model $\overline{\mathbb{CP}}^{3}$ and the cubic model ST$[2,6]$, which can be embedded in 10-dimensional Heterotic Supergravity. These solutions are based on the multicenter dyon recently discovered by one of us, which solves the SU$(2)$ Bogomol'nyi and dyon equations on $\mathbb{E}^{3}$. In contrast to the well-known Abelian multicenter solutions, the relative positions of the non-Abelian black hole centers are unconstrained. We study necessary conditions on the parameters of the solutions that ensure the regularity of the metric. In the case of the $\overline{\mathbb{CP}}^{3}$ model we show that it is enough to require the positivity of the "masses" of the individual black holes, the finiteness of each of their entropies and their superadditivity. In the case of the $ST[2,6]$ model we have not been able to show that analogous conditions are sufficient, but we give an explicit example of a regular solution describing thousands of non-Abelian dyonic black holes in equilibrium at arbitrary relative positions. We also construct non-Abelian solutions that interpolate smoothly between just two aDS$_{2}\times$S$^{2}$ vacua with different radii (dumbbell solutions).


title: Extended Main Sequence Turnoffs in Intermediate-Age Star Clusters: Stellar Rotation diminishes, but does not eliminate, Age Spreads

abstract: Extended main sequence turn-off (eMSTO) regions are a common feature in color-magnitude diagrams of young and intermediate-age star clusters in the Magellanic Clouds. The nature of eMSTOs remains debated in the literature. The currently most popular scenarios are extended star formation activity and ranges of stellar rotation rates. Here we study details of differences in MSTO morphology expected from spreads in age versus spreads in rotation rates, using Monte Carlo simulations with the Geneva SYCLIST isochrone models that include the effects of stellar rotation. We confirm a recent finding of Niederhofer et al. that a distribution of stellar rotation velocities yields an MSTO extent that is proportional to the cluster age, as observed. However, we find that stellar rotation yields MSTO crosscut widths that are generally smaller than observed ones at a given age. We compare the simulations with high-quality Hubble Space Telescope data of NGC 1987 and NGC 2249, which are the two only relatively massive star clusters with an age of ~1 Gyr for which such data is available. We find that the distribution of stars across the eMSTOs of these clusters cannot be explained solely by a distribution of stellar rotation velocities, unless the orientations of rapidly rotating stars are heavily biased towards an equator-on configuration. Under the assumption of random viewing angles, stellar rotation can account for ~60% and ~40% of the observed FWHM widths of the eMSTOs of NGC 1987 and NGC 2249, respectively. In contrast, a combination of distributions of stellar rotation velocities and stellar ages fits the observed eMSTO morphologies very well.


title: Reduced Electron Exposure for Energy-Dispersive Spectroscopy using Dynamic Sampling

abstract: Analytical electron microscopy and spectroscopy of biological specimens, polymers, and other beam sensitive materials has been a challenging area due to irradiation damage. There is a pressing need to develop novel imaging and spectroscopic imaging methods that will minimize such sample damage as well as reduce the data acquisition time. The latter is useful for high-throughput analysis of materials structure and chemistry. In this work, we present a novel machine learning based method for dynamic sparse sampling of EDS data using a scanning electron microscope. Our method, based on the supervised learning approach for dynamic sampling algorithm and neural networks based classification of EDS data, allows a dramatic reduction in the total sampling of up to 90%, while maintaining the fidelity of the reconstructed elemental maps and spectroscopic data. We believe this approach will enable imaging and elemental mapping of materials that would otherwise be inaccessible to these analysis techniques.


title: High Visibility in Two-Color Above-Threshold Photoemission from Tungsten Nanotips in a Coherent Control Scheme

abstract: In this article we present coherent control of above-threshold photoemission from a tungsten nanotip achieving nearly perfect modulation. Depending on the pulse delay between fundamental (1560 nm) and second harmonic (780 nm) pulses of a femtosecond fiber laser at the nanotip, electron emission is significantly enhanced or depressed during temporal overlap. Electron emission is studied as a function of pulse delay, optical near-field intensities, DC bias field, and final photoelectron energy. Under optimized conditions modulation amplitudes of the electron emission of 97.5% are achieved. Experimental observations are discussed in the framework of quantum- pathway interference supported by local density of states (LDOS) simulations.


title: Simulation of Drop Impact on a Hot Wall using SPH Method with Peng-Robinson Equation of State

abstract: This study presents a smoothed particle hydrodynamics (SPH) method with Peng-Robinson equation of state for simulating drop vaporization and drop impact on a hot surface. The conservation equations of momentum and energy and Peng-Robinson equation of state are applied to describe both the liquid and gas phases. The governing equations are solved numerically by the SPH method. The phase change between the liquid and gas phases are simulated directly without using any phase change models. The numerical method is validated by comparing numerical results with analytical solutions for the vaporization of n-heptane drops at different temperatures. Using the SPH method, the processes of n-heptane drops impacting on a solid wall with different temperatures are studied numerically. The results show that the size of the film formed by drop impact decreases when temperature increases. When the temperature is high enough, the drop will rebound.


title: Positive curvature and scalar field tunneling in the landscape

abstract: We present a model of vacuum tunneling through a classically forbidden region where a scalar field changes its value simultaneously over the entire volume of a (meta)stable ancestor vacuum with spherical curvature. The tunneling leaves the geometry unchanged but rearranges the energetic contributions of the matter sources, leading to an inflating solution with residual positive curvature. We show that there exists a parametric regime where this solution is self-consistent and dominates the overall tunneling rate. We conclude that an experimental detection of positive curvature, while not necessarily likely, therefore does not rule out the possibility that our present observer patch originated from semiclassical vacuum tunneling in a string or field theoretic landscape.


title: Big Bang Nucleosynthesis with Stable $^8$Be and the Primordial Lithium Problem

abstract: A change in the fundamental constants of nature or plasma effects in the early universe could stabilize $^8$Be against decay into two $^4$He nuclei. Coc et al. examined this effect on big bang nucleosynthesis as a function of $B_8$, the mass difference between two $^4$He nuclei and a single $^8$Be nucleus, and found no effects for $B_8 \le 100$ keV. Here we examine stable $^8$Be with larger $B_8$ and also allow for a variation in the rate for $^4$He + $^4$He $\longrightarrow$ $^8$Be to determine the threshold for interesting effects. We find no change to standard big bang nucleosynthesis for $B_8 < 1$ MeV. For $B_8 \gtrsim 1$ MeV and a sufficiently large reaction rate, a significant fraction of $^4$He is burned into $^8$Be, which fissions back into $^4$He when $B_8$ assumes its present-day value, leaving the primordial $^4$He abundance unchanged. However, this sequestration of $^4$He results in a decrease in the primordial $^7$Li abundance. Primordial abundances of $^7$Li consistent with observationally-inferred values can be obtained for reaction rates similar to those calculated for the present-day (unbound $^8$Be) case. Even for the largest binding energies and largest reaction rates examined here, only a small fraction of $^8$Be is burned into heavier elements, consistent with earlier studies. There is no change in the predicted deuterium abundance for any model we examined.


title: Supergravity Model of Inflation and Explaining IceCube HESE Data via PeV Dark Matter Decay

abstract: We construct an unified model of inflation and PeV dark matter with an appropriate choice of no-scale K\"ahler potential, superpotential and gauge kinetic function in terms of MSSM fields and hidden sector Polonyi field. The model is consistent with the CMB observations and can explain the PeV neutrino flux observed at IceCube HESE. A Starobinsky like Higgs-sneutrino plateau inflation is obtained from the $D$-term SUGRA potential while $F$-term being subdominant during inflation. To get PeV dark matter, SUSY breaking at PeV scale is achieved through Polonyi field. This sets the scale for soft SUSY breaking parameters $m_0, m_{1/2}, A_0 $ at the GUT scale in terms of the parameters of the model. The low energy particles spectrum is obtained by running the RGEs. We show that the $\sim$125 GeV higgs and the gauge coupling unification can be obtained in this model. The $6$ PeV bino-type dark matter is a subdominant fraction ($\sim 11\%$) of the relic density and its decay gives the PeV scale neutrino flux observed at IceCube by appropriately choosing the couplings of the $R$-parity violating operators. Also we find that there is a degeneracy in scalar field parameters $(\gamma, \beta)$ and coupling $\zeta$ value in producing the correct amplitude of CMB power spectrum. However the value of parameter $\tan(\beta)=1.8$, which is tightly fixed from the requirement of PeV scale SUSY breaking, removes the degeneracy in the values of the scalar field parameters to provide a unique solution for inflation. In this way it brings the explanation for dark matter, PeV neutrinos and inflation within the same framework.


title: Estimating the unseen from multiple populations

abstract: Given samples from a distribution, how many new elements should we expect to find if we continue sampling this distribution? This is an important and actively studied problem, with many applications ranging from unseen species estimation to genomics. We generalize this extrapolation and related unseen estimation problems to the multiple population setting, where population $j$ has an unknown distribution $D_j$ from which we observe $n_j$ samples. We derive an optimal estimator for the total number of elements we expect to find among new samples across the populations. Surprisingly, we prove that our estimator's accuracy is independent of the number of populations. We also develop an efficient optimization algorithm to solve the more general problem of estimating multi-population frequency distributions. We validate our methods and theory through extensive experiments. Finally, on a real dataset of human genomes across multiple ancestries, we demonstrate how our approach for unseen estimation can enable cohort designs that can discover interesting mutations with greater efficiency.


title: Capacity, Fidelity, and Noise Tolerance of Associative Spatial-Temporal Memories Based on Memristive Neuromorphic Network

abstract: We have calculated the key characteristics of associative (content-addressable) spatial-temporal memories based on neuromorphic networks with restricted connectivity - "CrossNets". Such networks may be naturally implemented in nanoelectronic hardware using hybrid CMOS/memristor circuits, which may feature extremely high energy efficiency, approaching that of biological cortical circuits, at much higher operation speed. Our numerical simulations, in some cases confirmed by analytical calculations, have shown that the characteristics depend substantially on the method of information recording into the memory. Of the four methods we have explored, two look especially promising - one based on the quadratic programming, and the other one being a specific discrete version of the gradient descent. The latter method provides a slightly lower memory capacity (at the same fidelity) then the former one, but it allows local recording, which may be more readily implemented in nanoelectronic hardware. Most importantly, at the synchronous retrieval, both methods provide a capacity higher than that of the well-known Ternary Content-Addressable Memories with the same number of nonvolatile memory cells (e.g., memristors), though the input noise immunity of the CrossNet memories is somewhat lower.


title: Buffer Size for Routing Limited-Rate Adversarial Traffic

abstract: We consider the slight variation of the adversarial queuing theory model, in which an adversary injects packets with routes into the network subject to the following constraint: For any link $e$, the total number of packets injected in any time window $[t,t')$ and whose route contains $e$, is at most $\rho(t'-t)+\sigma$, where $\rho$ and $\sigma$ are non-negative parameters. Informally, $\rho$ bounds the long-term rate of injections and $\sigma$ bounds the "burstiness" of injection: $\sigma=0$ means that the injection is as smooth as it can be. It is known that greedy scheduling of the packets (under which a link is not idle if there is any packet ready to be sent over it) may result in $\Omega(n)$ buffer size even on an $n$-line network and very smooth injections ($\sigma=0$). In this paper we propose a simple non-greedy scheduling policy and show that, in a tree where all packets are destined at the root, no buffer needs to be larger than $\sigma+2\rho$ to ensure that no overflows occur, which is optimal in our model. The rule of our algorithm is to forward a packet only if its next buffer is completely empty. The policy is centralized: in a single step, a long "train" of packets may progress together. We show that in some sense central coordination is required, by presenting an injection pattern with $\sigma=0$ for the $n$-node line that results in $\Omega(n)$ packets in a buffer if local control is used, even for the more sophisticated "downhill" algorithm, which forwards a packet only if its next buffer is less occupied than its current one.


title: Generalizing spin and pseudospin symmetries for relativistic spin 1/2 fermions

abstract: We propose a generalization of pseudospin and spin symmetries, the SU(2) symmetries of Dirac equation with scalar and vector mean-field potentials originally found independently in the 70's by Smith and Tassie, and Bell and Ruegg. As relativistic symmetries, they have been extensively researched and applied to several physical systems for the last 18 years. The main feature of these symmetries is the suppression of the spin-orbit coupling either in the upper or lower components of the Dirac spinor, thereby turning the respective second-order equations into Schr\"odinger-like equations, i.e, without a matrix structure. In this paper we use the original formalism of Bell and Ruegg to derive general requirements for the Lorentz structures of potentials in order to have these SU(2) symmetries in the Dirac equation, again allowing for the suppression of the matrix structure of the second-order equation of either the upper or lower components of the Dirac spinor. Furthermore, we derive equivalent conditions for spin and pseudospin symmetries with 2- and 1-dimensional potentials and list some possible candidates for 3, 2, and 1 dimensions. We suggest applications for physical systems in three and two dimensions, namely electrons in graphene.


title: Gradient Coding from Cyclic MDS Codes and Expander Graphs

abstract: Gradient coding is a technique for straggler mitigation in distributed learning. In this paper we design novel gradient codes using tools from classical coding theory, namely, cyclic MDS codes, which compare favorably with existing solutions, both in the applicable range of parameters and in the complexity of the involved algorithms. Second, we introduce an approximate variant of the gradient coding problem, in which we settle for approximate gradient computation instead of the exact one. This approach enables graceful degradation, i.e., the $\ell_2$ error of the approximate gradient is a decreasing function of the number of stragglers. Our main result is that normalized adjacency matrices of expander graphs yield excellent approximate gradient codes, which enable significantly less computation compared to exact gradient coding, and guarantee faster convergence than trivial solutions under standard assumptions. We experimentally test our approach on Amazon EC2, and show that the generalization error of approximate gradient coding is very close to the full gradient while requiring significantly less computation from the workers.


title: Intuitionistic modal logic based on neighborhood semantics without superset axiom

abstract: In this paper we investigate certain systems of propositional intuitionistic modal logic defined semantically in terms of neighborhood structures. We discuss various restrictions imposed on those frames but our constant approach is to discard superset axiom. Such assumption allows us to think about specific modalities $\Delta, \nabla$ and new functor $\rightsquigarrow$ depending on the notion of maximal neighborhood. We show how it is possible to treat our models as bi-relational ones. We prove soundness and completeness of proposed axiomatization by means of canonical model. Moreover, we show that finite model property holds. Then we describe properties of bounded morphism, behavioral equivalence, bisimulation and n-bisimulation. Finally, we discuss further researches (like some interesting classical cases and particular topological issues).


title: Ideal classes and Cappell-Shaneson homotopy 4-spheres

abstract: Gompf proposed a conjecture on Cappell-Shaneson matrices whose affirmative answer implies that all Cappell-Shaneson homotopy 4-spheres are diffeomorphic to the standard 4-sphere. We study Gompf conjecture on Cappell-Shaneson matrices using various algebraic number theoretic techniques. We find a hidden symmetry between trace $n$ Cappell-Shaneson matrices and trace $5-n$ Cappell-Shaneson matrices which was suggested by Gompf experimentally. Using this symmetry, we prove that Gompf conjecture for the trace $n$ case is equivalent to the trace $5-n$ case. We confirm Gompf conjecture for the special cases that $-64\leq \textrm{trace}\leq 69$ and corresponding Cappell-Shaneson homotopy 4-spheres are diffeomorphic to the standard 4-sphere. We also give a new infinite family of Cappell-Shaneson spheres which are diffeomorphic to the standard 4-sphere.


title: Two Non-Commutative Binomial Theorems

abstract: We derive two formulae for (A + B)^n, where A and B are elements in a non-commutative, associative algebra with identity.


title: Duality for finite Gelfand pairs

abstract: Let $\mathrm{G}$ be a split reductive group, $K$ be a non-Archimedean local field, and $O$ be its ring of integers. Satake isomorphism identifies the algebra of compactly supported invariants $\mathbb{C}_c[\mathrm{G}(K)/\mathrm{G}(O))]^{\mathrm{G}(O)}$ with a complexification of the algebra of characters of finite-dimensional representations $\mathcal{O}(\mathrm{G}^L(\mathbb{C}))^{\mathrm{G}^L(\mathbb{C})}$ of the Langlands dual group. In this note we report on the results of the study of analogues of such an isomorphism for finite groups. In our setup we replaced Gelfand pair $\mathrm{G}(O)\subset \mathrm{G}(K)$ by a finite pair $H\subset G$. It is convenient to rewrite the character side of the isomorphism as $\mathcal{O}(\mathrm{G}^L(\mathbb{C}))^{\mathrm{G}^L(\mathbb{C})}=\mathcal{O}((\mathrm{G}^L(\mathbb{C})\times \mathrm{G}^L(\mathbb{C}))/\mathrm{G}^L(\mathbb{C}))^{\mathrm{G}^L(\mathbb{C})}$. We replace diagonal Gelfand pair $\mathrm{G}^L(\mathbb{C})\subset \mathrm{G}^L(\mathbb{C})\times \mathrm{G}^L(\mathbb{C})$ by a dual finite pair $\check{H}\subset \check{G}$ and use Satake isomorphism as a defining property of the duality. In this text we make a preliminary study of such duality and compute a number of nontrivial examples of dual pairs $(H,G)$ and $(\check{H}, \check{G})$. We discuss a possible relation of our constructions to String Topology.


title: Simplicial Structures Over the 3-Sphere and Generalized Higher Order Hochschild Homology

abstract: In this paper we investigate the simplicial structure of a chain complex associated to the higher order Hochschild homology over the $3$-sphere. We also introduce the tertiary Hochschild homology corresponding to a quintuple $(A,B,C,\varepsilon,\theta)$, which becomes natural after we organize the elements in a convenient manner. We establish these results by way of a bar-like resolution in the context of simplicial modules. Finally, we generalize the higher order Hochschild homology over a trio of simplicial sets, which also grants natural geometric realizations.


title: Discordance between genomic divergence and phenotypic variation in a rapidly evolving avian genus (Motacilla)

abstract: Generally, genotypes and phenotypes are expected to be spatially congruent, however, in widespread species complexes with few barriers to dispersal, multiple contact zones, and limited reproductive isolation, discordance between phenotypes and phylogeographic groups is more probable. Wagtails (Aves: Motacilla) are a genus of birds with striking plumage pattern variation across Eurasia. Up to 13 subspecies are recognized within a single species, yet previous studies using mitochondrial DNA have supported phylogeographic groups that are inconsistent with subspecies plumage characteristics. In this study, we investigate the link between phenotypes and genotype by comparing populations thought to be at different stages along the speciation continuum. We take a phylogeographic approach by estimating population structure, testing for isolation by distance, conducting demographic modeling, and estimating the first time-calibrated species tree for the genus. Our study provides strong evidence for species-level patterns of differentiation in wagtails, however population-level differentiation is less pronounced. We find evidence that three of four widespread Eurasian species exhibit an east-west divide that contradicts both subspecies taxonomy and phenotypic variation. Both the geographic location of this divide and time estimates from demographic models are overlapping in two sympatric species, indicating that coincident Pleistocene events shaped their histories.


title: Mechanics Automatically Recognized via Interactive Observation: Jumping

abstract: Jumping has been an important mechanic since its introduction in Donkey Kong. It has taken a variety of forms and shown up in numerous games, with each jump having a different feel. In this paper, we use a modified Nintendo Entertainment System (NES) emulator to semi-automatically run experiments on a large subset (30%) of NES platform games. We use these experiments to build models of jumps from different developers, series, and games across the history of the console. We then examine these models to gain insights into different forms of jumping and their associated feel.


title: Bosonic Tensor Models at Large $N$ and Small $\epsilon$

abstract: We study the spectrum of the large $N$ quantum field theory of bosonic rank-$3$ tensors, whose quartic interactions are such that the perturbative expansion is dominated by the melonic diagrams. We use the Schwinger-Dyson equations to determine the scaling dimensions of the bilinear operators of arbitrary spin. Using the fact that the theory is renormalizable in $d=4$, we compare some of these results with the $4-\epsilon$ expansion, finding perfect agreement. This helps elucidate why the dimension of operator $\phi^{abc}\phi^{abc}$ is complex for $d<4$: the large $N$ fixed point in $d=4-\epsilon$ has complex values of the couplings for some of the $O(N)^3$ invariant operators. We show that a similar phenomenon holds in the $O(N)^2$ symmetric theory of a matrix field $\phi^{ab}$, where the double-trace operator has a complex coupling in $4-\epsilon$ dimensions. We also study the spectra of bosonic theories of rank $q-1$ tensors with $\phi^q$ interactions. In dimensions $d>1.93$ there is a critical value of $q$, above which we have not found any complex scaling dimensions. The critical value is a decreasing function of $d$, and it becomes $6$ in $d\approx 2.97$. This raises a possibility that the large $N$ theory of rank-$5$ tensors with sextic potential has an IR fixed point which is free of perturbative instabilities for $2.97<d<3$. This theory may be studied using renormalized perturbation theory in $d=3-\epsilon$.


title: Magnetised Polish doughnuts revisited

abstract: We discuss a procedure to build new sequences of magnetised, equilibrium tori around Kerr black holes which combines two approaches previously considered in the literature. For simplicity we assume that the test-fluid approximation holds, and hence we neglect the self-gravity of the fluid. The models are built assuming a particular form of the angular momentum distribution from which the location and morphology of equipotential surfaces can be computed. This ansatz includes, in particular, the constant angular momentum case originally employed in the construction of thick tori - or Polish doughnuts - and it has already been used to build equilibrium sequences of purely hydrodynamical models. We discuss the properties of the new models and their dependence on the initial parameters. These new sequences can be used as initial data for magnetohydrodynamical evolutions in general relativity.


title: Effective Interactions in a Graphene Layer Induced by the Proximity to a Ferromagnet

abstract: The proximity-induced couplings in graphene due to the vicinity of a ferromagnetic insulator are analyzed. We combine general symmetry principles and simple tight-binding descriptions to consider different orientations of the magnetization. We find that, in addition to a simple exchange field, a number of other terms arise. Some of these terms act as magnetic orbital couplings, and others are proximity-induced spin-orbit interactions. The couplings are of similar order of magnitude, and depend on the orientation of the magnetization. A variety of phases, and anomalous Hall effect regimes, are possible.


title: Cognitive Biases in Software Engineering: A Systematic Mapping Study

abstract: One source of software project challenges and failures is the systematic errors introduced by human cognitive biases. Although extensively explored in cognitive psychology, investigations concerning cognitive biases have only recently gained popularity in software engineering (SE) research. This paper therefore systematically maps, aggregates and synthesizes the literature on cognitive biases in software engineering to generate a comprehensive body of knowledge, understand state of the art research and provide guidelines for future research and practise. Focusing on bias antecedents, effects and mitigation techniques, we identified 65 articles, which investigate 37 cognitive biases, published between 1990 and 2016. Despite strong and increasing interest, the results reveal a scarcity of research on mitigation techniques and poor theoretical foundations in understanding and interpreting cognitive biases. Although bias-related research has generated many new insights in the software engineering community, specific bias mitigation techniques are still needed for software professionals to overcome the deleterious effects of cognitive biases on their work.


title: Lyapunov Conditions for Differentiability of Markov Chain Expectations: the Absolutely Continuous Case

abstract: We consider a family of Markov chains whose transition dynamics are affected by model parameters. Understanding the parametric dependence of (complex) performance measures of such Markov chains is often of significant interest. The derivatives of the performance measures w.r.t. the parameters play important roles, for example, in numerical optimization of the performance measures, and quantification of the uncertainties in the performance measures when there are uncertainties in the parameters from the statistical estimation procedures. In this paper, we establish conditions that guarantee the differentiability of various types of intractable performance measures---such as the stationary and random horizon discounted performance measures---of general state space Markov chains and provide probabilistic representations for the derivatives.


title: Particle Simulation of Fractional Diffusion Equations

abstract: This work explores different particle-based approaches to the simulation of one-dimensional fractional subdiffusion equations in unbounded domains. We rely on smooth particle approximations, and consider four methods for estimating the fractional diffusion term. The first method is based on direct differentiation of the particle representation, it follows the Riesz definition of the fractional derivative and results in a non-conservative scheme. The other three methods follow the particle strength exchange (PSE) methodology and are by construction conservative, in the sense that the total particle strength is time invariant. The first PSE algorithm is based on using direct differentiation to estimate the fractional diffusion flux, and exploiting the resulting estimates in an integral representation of the divergence operator. Meanwhile, the second one relies on the regularized Riesz representation of the fractional diffusion term to derive a suitable interaction formula acting directly on the particle representation of the diffusing field. A third PSE construction is considered that exploits the Green's function of the fractional diffusion equation. The performance of all four approaches is assessed for the case of a one-dimensional diffusion equation with constant diffusivity. This enables us to take advantage of known analytical solutions, and consequently conduct a detailed analysis of the performance of the methods. This includes a quantitative study of the various sources of error, namely filtering, quadrature, domain truncation, and time integration, as well as a space and time self-convergence analysis. These analyses are conducted for different values of the order of the fractional derivatives, and computational experiences are used to gain insight that can be used for generalization of the present constructions.


title: Independence, Conditionality and Structure of Dempster-Shafer Belief Functions

abstract: Several approaches of structuring (factorization, decomposition) of Dempster-Shafer joint belief functions from literature are reviewed with special emphasis on their capability to capture independence from the point of view of the claim that belief functions generalize bayes notion of probability. It is demonstrated that Zhu and Lee's {Zhu:93} logical networks and Smets' {Smets:93} directed acyclic graphs are unable to capture statistical dependence/independence of bayesian networks {Pearl:88}. On the other hand, though Shenoy and Shafer's hypergraphs can explicitly represent bayesian network factorization of bayesian belief functions, they disclaim any need for representation of independence of variables in belief functions. Cano et al. {Cano:93} reject the hypergraph representation of Shenoy and Shafer just on grounds of missing representation of variable independence, but in their frameworks some belief functions factorizable in Shenoy/Shafer framework cannot be factored. The approach in {Klopotek:93f} on the other hand combines the merits of both Cano et al. and of Shenoy/Shafer approach in that for Shenoy/Shafer approach no simpler factorization than that in {Klopotek:93f} approach exists and on the other hand all independences among variables captured in Cano et al. framework and many more are captured in {Klopotek:93f} approach.%


title: The Discrete-Time Geometric Maximum Principle

abstract: We establish a variety of results extending the well-known Pontryagin maximum principle of optimal control to discrete-time optimal control problems posed on smooth manifolds. These results are organized around a new theorem on critical and approximate critical points for discrete-time geometric control systems. We show that this theorem can be used to derive Lie group variational integrators in Hamiltonian form; to establish a maximum principle for control problems in the absence of state constraints; and to provide sufficient conditions for exact penalization techniques in the presence of state or mixed constraints. Exact penalization techniques are used to study sensitivity of the optimal value function to constraint perturbations and to prove necessary optimality conditions, including in the form of a maximum principle, for discrete-time geometric control problems with state or mixed constraints.


title: Cryogenic LED pixel-to-frequency mapper for kinetic inductance detector arrays

abstract: We present a cryogenic wafer mapper based on light emitting diodes (LEDs) for spatial mapping of a large microwave kinetic inductance detector (MKID) array. In this scheme, an array of LEDs, addressed by DC wires and collimated through horns onto the detectors, is mounted in front of the detector wafer. By illuminating each LED individually and sweeping the frequency response of all the resonators, we can unambiguously correspond a detector pixel to its measured resonance frequency. We have demonstrated mapping a 76.2 mm 90-pixel MKID array using a mapper containing 126 LEDs with 16 DC bias wires. With the frequency to pixel-position correspondence data obtained by the LED mapper, we have found a radially position-dependent frequency non-uniformity < 1.6% over the 76.2 mm wafer. Our LED wafer mapper has no moving parts and is easy to implement. It may find broad applications in superconducting detector and quantum computing/information experiments.


title: Filamentary superconductivity in semiconducting policrystalline ZrSe2 compound with Zr vacancies

abstract: ZrSe2 is a band semiconductor studied long time ago. It has interesting electronic properties, and because its layers structure can be intercalated with different atoms to change some of the physical properties. In this investigation we found that Zr deficiencies alter the semiconducting behavior and the compound can be turned into a superconductor. In this paper we report our studies related to this discovery. The decreasing of the number of Zr atoms in small proportion according to the formula ZrxSe2, where x is varied from about 8.1 to 8.6 K, changing the semiconducting behavior to a superconductor with transition temperatures ranging between 7.8 to 8.5 K, it depending of the deficiencies. Outside of those ranges the compound behaves as semiconducting with the properties already known. In our experiments we found that this new superconductor has only a very small fraction of superconducting material determined by magnetic measurements with applied magnetic field of 10 Oe. Our conclusions is that superconductivity is filamentary. However, in one studied sample the fraction was about 10.2 %, whereas in others is only about 1 % or less. We determined the superconducting characteristics; the critical fields that indicate a type two superonductor with Ginzburg-Landau ? parameter of the order about 2.7. The synthesis procedure is quite normal fol- lowing the conventional solid state reaction. In this paper are included, the electronic characteristics, transition temperature, and evolution with temperature of the critical fields.


title: Heavy traffic analysis of a polling model with retrials and glue periods

abstract: We present a heavy traffic analysis of a single-server polling model, with the special features of retrials and glue periods. The combination of these features in a polling model typically occurs in certain optical networking models, and in models where customers have a reservation period just before their service period. Just before the server arrives at a station there is some deterministic glue period. Customers (both new arrivals and retrials) arriving at the station during this glue period will be served during the visit of the server. Customers arriving in any other period leave immediately and will retry after an exponentially distributed time. As this model defies a closed-form expression for the queue length distributions, our main focus is on their heavy-traffic asymptotics, both at embedded time points (beginnings of glue periods, visit periods and switch periods) and at arbitrary time points. We obtain closed-form expressions for the limiting scaled joint queue length distribution in heavy traffic and use these to accurately approximate the mean number of customers in the system under different loads.


title: Foresight: Recommending Visual Insights

abstract: Current tools for exploratory data analysis (EDA) require users to manually select data attributes, statistical computations and visual encodings. This can be daunting for large-scale, complex data. We introduce Foresight, a system that helps the user rapidly discover visual insights from large high-dimensional datasets. Formally, an "insight" is a strong manifestation of a statistical property of the data, e.g., high correlation between two attributes, high skewness or concentration about the mean of a single attribute, a strong clustering of values, and so on. For each insight type, Foresight initially presents visualizations of the top k instances in the data, based on an appropriate ranking metric. The user can then look at "nearby" insights by issuing "insight queries" containing constraints on insight strengths and data attributes. Thus the user can directly explore the space of insights, rather than the space of data dimensions and visual encodings as in other visual recommender systems. Foresight also provides "global" views of insight space to help orient the user and ensure a thorough exploration process. Furthermore, Foresight facilitates interactive exploration of large datasets through fast, approximate sketching.


title: Transient nature of J195754+353513

abstract: We have searched for transient and/or variable radio sources in the field of Galactic micro-quasar Cygnus X-1 in 1.4 GHz (L band) using data from the Karl G. Jansky Very Large Array (JVLA). We used twenty years of data between 1983 and 2003. We found a source J195754+353513 showing transient behavior. The source was also mentioned earlier in NVSS and WENSS catalog but its transient nature was not reported earlier. The source is located 23.8 arcminutes far from Cygnus X-1. It is detected many times during the span of our study and it varied between less than 0.3 mJy to 201 mJy. J195754+353513 also showed high intra-day variability. In one occasion, the source rose from ~20 mJy to ~180 mJy 700 second. For limited number of cases circularly polarized emission could be detected from the source with V/I vary between 0.15 to 0.25. 2MASS J19575420+3535152 may be the near-infrared counterpart of the source. We compared the properties of the source with other Galactic transient sources having similar properties. The nature of the source is still unknown. We speculate on its possible nature.

categories: astro-ph.HE astro-ph.GA astro-ph.IM astro-ph.SR


title: The SAMI Galaxy Survey: Disk-halo interactions in radio-selected star-forming galaxies

abstract: In this paper, we compare the radio emission at 1.4 GHz with optical outflow signatures of edge-on galaxies. We report observations of six edge-on star-forming galaxies in the Sydney-AAO Multi-object Integral-field spectrograph (SAMI) Galaxy Survey with 1.4 GHz luminosities $> 1\times10^{21}$ W Hz$^{-1}$. Extended minor axis optical emission is detected with enhanced \nii/H$\alpha$ line ratios and velocity dispersions consistent with galactic winds in three of six galaxies. These galaxies may host outflows driven by a combination of thermal and cosmic ray processes. We find that galaxies with the strongest wind signatures have extended radio morphologies. Our results form a baseline for understanding the driving mechanisms of galactic winds.


title: Flux sensitivity requirements for the detection of Lyman continuum radiation drop-ins from star-forming galaxies below redshifts of 3

abstract: Flux estimates for ionizing radiation escaping from star-forming galaxies with characteristic UV luminosities ($L^{*}_{1500(1+z)}$) derived from GALEX and the VIMOS-VLT Deep Survey, are presented as a function of redshift and assumed escape fraction. These estimates offer guidance to the design of instrumentation and observing strategies, be they spectroscopic or photometric, attempting to detect LyC escaping star-forming galaxies for redshifts $ z < 3$. Examples are given that relate the absolution escape fraction ($f^{e}_{LyC}$) of ionizing photons, integrated over the entire extreme UV (EUV) bandpass, to the relative escape fraction ($f^{e}_{900}$) observed just shortward of the ionization edge at 911.8 \AA\ as a function of HI, HeI, and HeII column densities. We find that for $\log{N_{HI}(cm^{-2})} \gtrsim 17.0$ $f^{e}_{LyC}$ is significantly greater than $f^{e}_{900}$. Detection of LyC "drop-ins" in the rest-frame EUV will provide enhanced fidelity to determinations of the integrated fraction of ionizing photons $f^{e}_{LyC}$ that escape star-forming galaxies and contribute to the meta-galactic ionizing background.


title: Identification and Interpretation of Belief Structure in Dempster-Shafer Theory

abstract: Mathematical Theory of Evidence called also Dempster-Shafer Theory (DST) is known as a foundation for reasoning when knowledge is expressed at various levels of detail. Though much research effort has been committed to this theory since its foundation, many questions remain open. One of the most important open questions seems to be the relationship between frequencies and the Mathematical Theory of Evidence. The theory is blamed to leave frequencies outside (or aside of) its framework. The seriousness of this accusation is obvious: (1) no experiment may be run to compare the performance of DST-based models of real world processes against real world data, (2) data may not serve as foundation for construction of an appropriate belief model. In this paper we develop a frequentist interpretation of the DST bringing to fall the above argument against DST. An immediate consequence of it is the possibility to develop algorithms acquiring automatically DST belief models from data. We propose three such algorithms for various classes of belief model structures: for tree structured belief networks, for poly-tree belief networks and for general type belief networks.


title: Relevance of the speed and direction of pulling in simple modular proteins

abstract: A theoretical analysis of the unfolding pathway of simple modular proteins in length- controlled pulling experiments is put forward. Within this framework, we predict the first module to unfold in a chain of identical units, emphasizing the ranges of pulling speeds in which we expect our theory to hold. These theoretical predictions are checked by means of steered molecular dynamics of a simple construct, specifically a chain composed of two coiled-coils motives, where anisotropic features are revealed. These simulations also allow us to give an estimate for the range of pulling velocities in which our theoretical approach is valid.


title: Almost complex structures on spheres

abstract: In this paper we review the well-known fact that the only spheres admitting an almost complex structure are S^2 and S^6. The proof described here uses characteristic classes and the Bott periodicity theorem in topological K-theory. This paper originates from the talk "Almost Complex Structures on Spheres" given by the second author at the MAM1 workshop "(Non)-existence of complex structures on S^6", held in Marburg from March 27th to March 30th, 2017. It is a review paper, and as such no result is intended to be original. We tried to produce a clear, motivated and as much as possible self-contained exposition.


title: Braid Group Representations from Braiding Gapped Boundaries of Dijkgraaf-Witten Theories

abstract: We study representations of the braid groups from braiding gapped boundaries of Dijkgraaf-Witten theories and their twisted generalizations, which are (twisted) quantum doubled topological orders in two spatial dimensions. We show that the resulting braid (pure braid) representations are all monomial with respect to some specific bases, hence all such representation images of the braid groups are finite groups. We give explicit formulas for the monomial matrices and the ground state degeneracy of the Kitaev models that are Hamiltonian realizations of Dijkgraaf-Witten theories. Our results imply that braiding gapped boundaries alone cannot provide universal gate sets for topological quantum computing with gapped boundaries.


title: On the phantom barrier crossing and the bounds on the speed of sound in non-minimal derivative coupling theories

abstract: In this paper we investigate the so called "phantom barrier crossing" issue in a cosmological model based in the scalar-tensor theory with non-minimal derivative coupling to the Einstein's tensor. Special attention will be paid to the physical bounds on the squared sound speed. The numeric results are geometrically illustrated by means of a qualitative procedure of analysis that is based on the mapping of the orbits in the phase plane onto the surfaces that represent physical quantities in the extended phase space, that is: the phase plane complemented with an additional dimension relative to the given physical parameter. We find that the cosmological model based in the non-minimal derivative coupling theory -- this includes both the quintessence and the pure derivative coupling cases -- has serious causality problems related with superluminal propagation of the scalar and tensor perturbations. Even more disturbing is the finding that, despite that the underlying theory is free of the Ostrogradsky instability, the corresponding cosmological model is plagued by the Laplacian (classical) instability related with negative squared sound speed. This instability leads to an uncontrollable growth of the energy density of the perturbations that is inversely proportional to their wavelength. We show that independent of the self-interaction potential, for the positive coupling the tensor perturbations propagate superluminally, while for the negative coupling a Laplacian instability arises. This latter instability invalidates the possibility for the model to describe the primordial inflation.


title: A Formal Framework to Characterize Interpretability of Procedures

abstract: We provide a novel notion of what it means to be interpretable, looking past the usual association with human understanding. Our key insight is that interpretability is not an absolute concept and so we define it relative to a target model, which may or may not be a human. We define a framework that allows for comparing interpretable procedures by linking it to important practical aspects such as accuracy and robustness. We characterize many of the current state-of-the-art interpretable methods in our framework portraying its general applicability.


title: Quantum-kinetic perspective on photovoltaic device operation in nanostructure-based solar cells

abstract: The implementation of a wide range of novel concepts for next-generation high-efficiency solar cells is based on nanostructures with configuration-tunable optoelectronic properties. On the other hand, effective nano-optical light-trapping concepts enable the use of ultra-thin absorber architectures. In both cases, the local density of electronic and optical states deviates strongly from that in a homogeneous bulk material. At the same time, non-local and coherent phenomena like tunneling or ballistic transport become increasingly relevant. As a consequence, the semi-classical, diffusive bulk picture conventionally assumed may no longer be appropriate to describe the physical processes of generation, transport, and recombination governing the photovoltaic operation of such devices. In this review, we provide a quantum-kinetic perspective on photovoltaic device operation that reaches beyond the limits of the standard simulation models for bulk solar cells. Deviations from bulk physics are assessed in ultra-thin film and nanostructure-based solar cell architectures by comparing the predictions of the semi-classical models for key physical quantities such as absorption coefficients, emission spectra, generation and recombination rates as well as potentials, densities and currents with the corresponding properties as given by a more fundamental description based on non-equilibrium quantum statistical mechanics. This advanced approach, while paving the way to a comprehensive quantum theory of photovoltaics, bridges simulations at microscopic material and macroscopic device levels by providing the charge carrier dynamics at the mesoscale.


title: Additive non-approximability of chromatic number in proper minor-closed classes

abstract: Robin Thomas asked whether for every proper minor-closed class C, there exists a polynomial-time algorithm approximating the chromatic number of graphs from C up to a constant additive error independent on the class C. We show this is not the case: unless P=NP, for every integer k>=1, there is no polynomial-time algorithm to color a K_{4k+1}-minor-free graph G using at most chi(G)+k-1 colors. More generally, for every k>=1 and 1<=\beta<=4/3, there is no polynomial-time algorithm to color a K_{4k+1}-minor-free graph G using less than beta.chi(G)+(4-3beta)k colors. As far as we know, this is the first non-trivial non-approximability result regarding the chromatic number in proper minor-closed classes. We also give somewhat weaker non-approximability bound for K_{4k+1}-minor-free graphs with no cliques of size 4. On the positive side, we present additive approximation algorithm whose error depends on the apex number of the forbidden minor, and an algorithm with additive error 6 under the additional assumption that the graph has no 4-cycles.


title: Efficient eigenvalue determination for arbitrary Pauli products based on generalized spin-spin interactions

abstract: Effective spin-spin interactions between N+1 qubits enable the determination of the eigenvalue of an arbitrary Pauli product of dimension N with a constant, small number of multi-qubit gates that is independent of N and encodes the eigenvalue in the measurement basis states of an extra ancilla qubit. Such interactions are available whenever qubits can be coupled to a shared harmonic oscillator, a situation that can be realized in several physical qubit implementations. For example, suitable interactions have already been realized for up to 14 qubits in ion traps. It should be possible to implement stabilizer codes for quantum error correction with a constant number of multi-qubit gates, in contrast to typical constructions using a number of two-qubit gates that increases as a function of N. The special case of finding the parity of N qubits only requires a small number of operations that is independent of N. This compares favorably to algorithms for computing the parity on conventional machines, which implies a genuine quantum advantage.


title: Constraints on Higgs effective couplings in $H\nu \bar{\nu}$ production of CLIC at 380 GeV

abstract: The potential of the $e^+e^-\to\nu \bar{\nu} H$ process in the first stage of CLIC considering center of mass energy of 380 GeV and assuming the baseline integrated luminosity of 500 fb$^{-1}$ is examined to probe CP conserving dimension-six operators in a model-independent Standard Model effective field theory framework. In the analysis, a detailed fast simulation on $e^+e^-\to\nu \bar{\nu} H$ signal processes and dominant backgrounds are performed including parton showering with PYTHIA and detector simulation based on ILD type detector with DELPHES in MadGraph. The obtained best limits on the $\bar{c}_{HB} $, $\bar{c}_{W}=- \bar{c}_{B}$ and $\bar{c}_{HW} $ are $[-4.82;1.53]\times10^{-2}$, $[-5.11;4.13]\times10^{-3}$ and $[-6.58;5.55]\times10^{-3}$, respectively.


title: Unsupervised Body Part Regression via Spatially Self-ordering Convolutional Neural Networks

abstract: Automatic body part recognition for CT slices can benefit various medical image applications. Recent deep learning methods demonstrate promising performance, with the requirement of large amounts of labeled images for training. The intrinsic structural or superior-inferior slice ordering information in CT volumes is not fully exploited. In this paper, we propose a convolutional neural network (CNN) based Unsupervised Body part Regression (UBR) algorithm to address this problem. A novel unsupervised learning method and two inter-sample CNN loss functions are presented. Distinct from previous work, UBR builds a coordinate system for the human body and outputs a continuous score for each axial slice, representing the normalized position of the body part in the slice. The training process of UBR resembles a self-organization process: slice scores are learned from inter-slice relationships. The training samples are unlabeled CT volumes that are abundant, thus no extra annotation effort is needed. UBR is simple, fast, and accurate. Quantitative and qualitative experiments validate its effectiveness. In addition, we show two applications of UBR in network initialization and anomaly detection.


title: A sharp Dirac-Erd\H{o}s type bound for large graphs

abstract: Let $k \geq 3$ be an integer, $h_{k}(G)$ be the number of vertices of degree at least $2k$ in a graph $G$, and $\ell_{k}(G)$ be the number of vertices of degree at most $2k-2$ in $G$. Dirac and Erd\H{o}s proved in 1963 that if $h_{k}(G) - \ell_{k}(G) \geq k^{2} + 2k - 4$, then $G$ contains $k$ vertex-disjoint cycles. For each $k\geq 2$, they also showed an infinite sequence of graphs $G_k(n)$ with $h_{k}(G_k(n)) - \ell_{k}(G_k(n)) = 2k-1$ such that $G_k(n)$ does not have $k$ disjoint cycles. Recently, the authors proved that, for $k \geq 2$, a bound of $3k$ is sufficient to guarantee the existence of $k$ disjoint cycles and presented for every $k$ a graph $G_0(k)$ with $h_{k}(G_0(k)) - \ell_{k}(G_0(k))=3k-1$ and no $k$ disjoint cycles. The goal of this paper is to refine and sharpen this result: We show that the Dirac-Erd\H{o}s construction is optimal in the sense that for every $k \geq 2$, there are only finitely many graphs $G$ with $h_{k}(G) - \ell_{k}(G) \geq 2k$ but no $k$ disjoint cycles. In particular, every graph $G$ with $|V(G)| \geq 19k$ and $h_{k}(G) - \ell_{k}(G) \geq 2k$ contains $k$ disjoint cycles.


title: Collective phases of identical particles interfering on linear multiports

abstract: We introduce collective geometric phases of bosons and fermions interfering on a linear unitary multiport, where each phase depends on the internal states of identical particles (i.e., not affected by the multiport) and corresponds to a cycle of the symmetric group. We show that quantum interference of $N$ particles in generic pure internal states, i.e., with no pair being orthogonal, is governed by $(N-1)(N-2)/2$ independent triad phases (each involving only three particles). The deterministic distinguishability, preventing quantum interference with two or three particles, allows for the genuine $(N\ge 4)$-particle phase (interference) on a multiport: setting each particle to be deterministically distinguishable from all others except two by their internal states allows for a novel (circle-dance) interference of $N\ge 4$ particles governed by a collective $N$-particle phase, while simultaneously preventing the $R$-particle interference for $3\le R\le N-1$. The genuine $N$-particle interference manifests the $N$th order quantum correlations between identical particles at a multiport output, it does not appear in the marginal probability for a subset of the particles, e.g., it cannot be detected if at least one of the particles is lost. This means that the collective phases are not detectable by the usual "quantumness" criteria based on the second-order quantum correlations. The results can be useful for quantum computation, quantum information, and other quantum technologies with single photons. \end{abstract}


title: The Generalized Nagell-Ljunggren Problem: Powers with Repetitive Representations

abstract: We consider a natural generalization of the Nagell-Ljunggren equation to the case where the qth power of an integer y, for q >= 2, has a base-b representation that consists of a length-l block of digits repeated n times, where n >= 2. Assuming the abc conjecture of Masser and Oesterl\'e, we completely characterize those triples (q, n, l) for which there are infinitely many solutions b. In all cases predicted by the abc conjecture, we are able (without any assumptions) to prove there are indeed infinitely many solutions.


title: Quantum dots with split enhancement gate tunnel barrier control

abstract: We introduce a silicon metal-oxide-semiconductor quantum dot architecture based on a single polysilicon gate stack. The elementary structure consists of two enhancement gates separated spatially by a gap, one gate forming a reservoir and the other a quantum dot. We demonstrate, in three devices based on two different versions of this elementary structure, that a wide range of tunnel rates is attainable while maintaining single-electron occupation. A characteristic change in slope of the charge transitions as a function of the reservoir gate voltage, attributed to screening from charges in the reservoir, is observed in all devices, and is expected to play a role in the sizable tuning orthogonality of the split enhancement gate structure. The all-silicon process is expected to minimize strain gradients from electrode thermal mismatch, while the single gate layer should avoid issues related to overlayers (e.g., additional dielectric charge noise) and help improve yield. Finally, reservoir gate control of the tunnel barrier has implications for initialization, manipulation and readout schemes in multi-quantum dot architectures.


title: Character bounds for finite groups of Lie type

abstract: We establish new bounds on character values and character ratios for finite groups $G$ of Lie type, which are considerably stronger than previously known bounds, and which are best possible in many cases. These bounds have the form $|\chi(g)| \le \chi(1)^{\alpha_g}$, and give rise to a variety of applications, for example to covering numbers and mixing times of random walks on such groups. In particular we deduce that, if $G$ is a classical group in dimension $n$, then, under some conditions on $G$ and $g \in G$, the mixing time of the random walk on $G$ with the conjugacy class of $g$ as a generating set is (up to a small multiplicative constant) $n/s$, where $s$ is the support of $g$.


title: ClustGeo: an R package for hierarchical clustering with spatial constraints

abstract: In this paper, we propose a Ward-like hierarchical clustering algorithm including spatial/geographical constraints. Two dissimilarity matrices $D_0$ and $D_1$ are inputted, along with a mixing parameter $\alpha \in [0,1]$. The dissimilarities can be non-Euclidean and the weights of the observations can be non-uniform. The first matrix gives the dissimilarities in the "feature space" and the second matrix gives the dissimilarities in the "constraint space". The criterion minimized at each stage is a convex combination of the homogeneity criterion calculated with $D_0$ and the homogeneity criterion calculated with $D_1$. The idea is then to determine a value of $\alpha$ which increases the spatial contiguity without deteriorating too much the quality of the solution based on the variables of interest i.e. those of the feature space. This procedure is illustrated on a real dataset using the R package ClustGeo.


title: Tidal stripping and the structure of dwarf galaxies in the Local Group

abstract: The shallow faint-end slope of the galaxy mass function is usually reproduced in $\Lambda$CDM galaxy formation models by assuming that the fraction of baryons that turns into stars drops steeply with decreasing halo mass and essentially vanishes in haloes with maximum circular velocities $V_{\rm max}<20$-$30$ km/s. Dark matter-dominated dwarfs should therefore have characteristic velocities of about that value, unless they are small enough to probe only the rising part of the halo circular velocity curve (i.e., half-mass radii, $r_{1/2}\ll 1$ kpc). Many dwarfs have properties in disagreement with this prediction: they are large enough to probe their halo $V_{\rm max}$ but their characteristic velocities are well below $20$ km/s. These `cold faint giants' (an extreme example is the recently discovered Crater 2 Milky Way satellite) can only be reconciled with our $\Lambda$CDM models if they are the remnants of once massive objects heavily affected by tidal stripping. We examine this possibility using the APOSTLE cosmological hydrodynamical simulations of the Local Group. Assuming that low velocity dispersion satellites have been affected by stripping, we infer their progenitor masses, radii, and velocity dispersions, and find them in remarkable agreement with those of isolated dwarfs. Tidal stripping also explains the large scatter in the mass discrepancy-acceleration relation in the dwarf galaxy regime: tides remove preferentially dark matter from satellite galaxies, lowering their accelerations below the $a_{\rm min}\sim 10^{-11} m/s^2$ minimum expected for isolated dwarfs. In many cases, the resulting velocity dispersions are inconsistent with the predictions from Modified Newtonian Dynamics, a result that poses a possibly insurmountable challenge to that scenario.


title: A Topologist's View of Kinematic Maps and Manipulation Complexity

abstract: In this paper we combine a survey of the most important topological properties of kinematic maps that appear in robotics, with the exposition of some basic results regarding the topological complexity of a map. In particular, we discuss mechanical devices that consist of rigid parts connected by joints and show how the geometry of the joints determines the forward kinematic map that relates the configuration of joints with the pose of the end-effector of the device. We explain how to compute the dimension of the joint space and describe topological obstructions for a kinematic map to be a fibration or to admit a continuous section. In the second part of the paper we define the complexity of a continuous map and show how the concept can be viewed as a measure of the difficulty to find a robust manipulation plan for a given mechanical device. We also derive some basic estimates for the complexity and relate it to the degree of instability of a manipulation plan.


title: kIP: a Measured Approach to IPv6 Address Anonymization

abstract: Privacy-minded Internet service operators anonymize IPv6 addresses by truncating them to a fixed length, perhaps due to long-standing use of this technique with IPv4 and a belief that it's "good enough." We claim that simple anonymization by truncation is suspect since it does not entail privacy guarantees nor does it take into account some common address assignment practices observed today. To investigate, with standard activity logs as input, we develop a counting method to determine a lower bound on the number of active IPv6 addresses that are simultaneously assigned, such as those of clients that access World-Wide Web services. In many instances, we find that these empirical measurements offer no evidence that truncating IPv6 addresses to a fixed number of bits, e.g., 48 in common practice, protects individuals' privacy. To remedy this problem, we propose kIP anonymization, an aggregation method that ensures a certain level of address privacy. Our method adaptively determines variable truncation lengths using parameter k, the desired number of active (rather than merely potential) addresses, e.g., 32 or 256, that can not be distinguished from each other once anonymized. We describe our implementation and present first results of its application to millions of real IPv6 client addresses active over a week's time, demonstrating both feasibility at large scale and ability to automatically adapt to each network's address assignment practice and synthesize a set of anonymous aggregates (prefixes), each of which is guaranteed to cover (contain) at least k of the active addresses. Each address is anonymized by truncating it to the length of its longest matching prefix in that set.


title: Markov Numbers, Mather's $\beta$ function and stable norm

abstract: V. Fock [7] introduced an interesting function $\psi(x)$, $x \in {\mathbb R}$ related to Markov numbers. We explain its relation to Federer-Gromov's stable norm and Mather's $\beta$-function, and use this to study its properties. We prove that $\psi$ and its natural generalisations are differentiable at every irrational $x$ and non-differentiable otherwise, by exploiting the relation with length of closed geodesics on the punctured or one-hole tori with the hyperbolic metric and the results by Bangert [3] and McShane- Rivin [19].


title: Autoencoder-augmented Neuroevolution for Visual Doom Playing

abstract: Neuroevolution has proven effective at many reinforcement learning tasks, but does not seem to scale well to high-dimensional controller representations, which are needed for tasks where the input is raw pixel data. We propose a novel method where we train an autoencoder to create a comparatively low-dimensional representation of the environment observation, and then use CMA-ES to train neural network controllers acting on this input data. As the behavior of the agent changes the nature of the input data, the autoencoder training progresses throughout evolution. We test this method in the VizDoom environment built on the classic FPS Doom, where it performs well on a health-pack gathering task.


title: Negative Sampling Improves Hypernymy Extraction Based on Projection Learning

abstract: We present a new approach to extraction of hypernyms based on projection learning and word embeddings. In contrast to classification-based approaches, projection-based methods require no candidate hyponym-hypernym pairs. While it is natural to use both positive and negative training examples in supervised relation extraction, the impact of negative examples on hypernym prediction was not studied so far. In this paper, we show that explicit negative examples used for regularization of the model significantly improve performance compared to the state-of-the-art approach of Fu et al. (2014) on three datasets from different languages.