title: Stable representation homology and Koszul duality

abstract: This paper is a sequel to [BKR], where we studied the derived affine scheme DRep_n(A) of the classical representation scheme Rep_n(A) for an associative k-algebra A. In [BKR], we have constructed canonical trace maps Tr_n(A): HC(A) -> H[DRep_n(A)]^GL extending the usual characters of representations to higher cyclic homology. This raises a question whether a well known theorem of Procesi [P] holds in the derived setting: namely, is the algebra homomorphism Sym[Tr_n(A)]: Sym[HC(A)] -> H[DRep_n(A)]^GL defined by Tr_n(A) surjective ? In the present paper, we answer this question for augmented algebras. Given such an algebra, we construct a canonical dense DG subalgebra DRep_\infty(A)^Tr of the topological DG algebra DRep_\infty(A)^{GL_\infty}. It turns out that on passing to the inverse limit (as n -> \infty), the family of maps Sym[Tr_n(A)] "stabilizes" to an isomorphism Sym[\bar{HC}(A)] = H[DRep_\infty(A)^Tr]. The derived version of Procesi's theorem does therefore hold in the limit. However, for a fixed (finite) n, there exist homological obstructions to the surjectivity of Sym[Tr_n(A)], and we show on simple examples that these obstructions do not vanish in general. We compare our result with the classical theorem of Loday-Quillen and Tsygan on stable homology of matrix Lie algebras. We show that the relative Chevalley-Eilenberg complex C(gl_\infty(A), gl_\infty(k); k) equipped with the natural coalgebra structure is Koszul dual to the DG algebra DRep_\infty(A)^Tr. We also extend our main results to bigraded DG algebras, in which case we show that DRep_{\infty}(A)^Tr = DRep_{\infty}(A)^GL_{\infty}. As an application, we compute the (bigraded) Euler characteristics of DRep_\infty(A)^GL_{\infty} and \bar{HC}(A) and derive some interesting combinatorial identities.


title: Laser-Induced Magnetic Nanostructures with Tunable Topological Properties

abstract: We report the creation and real-space observation of magnetic structures with well-defined topological properties and a lateral size as low as about 150 nm. They are generated in a thin ferrimagnetic film by ultrashort single optical laser pulses. Thanks to their topological properties, such structures can be classified as Skyrmions of a particular type that does not require an externally applied magnetic field for stabilization. Besides Skyrmions, we are able to generate magnetic features with topological characteristics that can be tuned by changing the laser fluence. The stability of such features is accounted for by an analytical model based on the interplay between the exchange and the magnetic dipole-dipole interactions


title: Block Triangular Preconditioning for Stochastic Galerkin Method

abstract: In this paper we study fast iterative solvers for the large sparse linear systems resulting from the stochastic Galerkin discretization of stochastic partial differential equations. A block triangular preconditioner is introduced and applied to the Krylov subspace methods, including the generalized minimum residual method and the generalized preconditioned conjugate gradient method. This preconditioner utilizes the special structures of the stochastic Galerkin matrices to achieve high efficiency. Spectral bounds for the preconditioned matrix are provided for convergence analysis. The preconditioner system can be solved approximately by geometric multigrid V-cycle. Numerical results indicate that the block triangular preconditioner has better performance than the traditional block diagonal preconditioner for stochastic problems with large variance.


title: Trouble With The Curve: Improving MLB Pitch Classification

abstract: The PITCHf/x database has allowed the statistical analysis of of Major League Baseball (MLB) to flourish since its introduction in late 2006. Using PITCHf/x, pitches have been classified by hand, requiring considerable effort, or using neural network clustering and classification, which is often difficult to interpret. To address these issues, we use model-based clustering with a multivariate Gaussian mixture model and an appropriate adjustment factor as an alternative to current methods. Furthermore, we describe a new pitch classification algorithm based on our clustering approach to address the problems of pitch misclassification. We illustrate our methods for various pitchers from the PITCHf/x database that covers a wide variety of pitch types.


title: Asynchronous Gossip-Based Random Projection Algorithms Over Networks

abstract: We consider a fully distributed constrained convex optimization problem over a multi-agent (no central coordinator) network. We propose an asynchronous gossip-based random projection (GRP) algorithm that solves the distributed problem using only local communications and computations. We analyze the convergence properties of the algorithm for an uncoordinated diminishing stepsize and a constant stepsize. For a diminishing stepsize, we prove that the iterates of all agents converge to the same optimal point with probability 1. For a constant stepsize, we establish an error bound on the expected distance from the iterates of the algorithm to the optimal point. We also provide simulation results on a distributed robust model predictive control problem.


title: Complementarity of dark matter direct detection: the role of bolometric targets

abstract: We study how the combined observation of dark matter in various direct detection experiments can be used to determine the phenomenological properties of WIMP dark matter: mass, spin-dependent (SD) and spin-independent (SI) scattering cross section off nucleons. A convenient choice of target materials, including nuclei that couple to dark matter particles through a significantly different ratio of SD vs SI interactions, could break the degeneracies in the determination of those parameters that a single experiment cannot discriminate. In this work we investigate different targets that can be used as scintillating bolometers and could provide complementary information to germanium and xenon detectors. We observe that Al2O3 and LiF bolometers could allow a good reconstruction of the DM properties over regions of the parameter space with a SD scattering cross section as small as 10^(-5) pb and a SI cross section as small as 5x10^(-10) pb for a 50 GeV WIMP. In the case of a CaWO4 bolometer the area in which full complementarity is obtained is smaller but we show that it can be used to determine the WIMP mass and its SI cross section. For each target we study the required exposure and background.


title: Prospects for precision measurements in nuclear beta decay at the LHC era

abstract: Precision measurements in nuclear beta decay offer a sensitive window to search for new physics beyond the standard electroweak model and allow also the determination of the fundamental weak vector coupling in processes involving the lightest quarks. Searches for new physics are also a strong motivation for experiments carried out at the high energy frontier reached at the most powerful particle colliders. It is instructive to confront results from the low energy and the high energy frontiers in order to look for possible complementarities and orient new avenues for experiments at low energies. We review here the status of constraints on new physics obtained from nuclear and neutron decays and compare them to those from other semi-leptonic processes and from the LHC. We stress the requirements of new precision experiments in beta decay in order to impact the search for new physics at the light of current and projected LHC results. We describe recent experimental results and ongoing developments in nuclear and neutron beta decay, with emphasis on their planned goals to improve present limits on exotic weak couplings.


title: Note: interpreting iterative methods convergence with diffusion point of view

abstract: In this paper, we explain the convergence speed of different iteration schemes with the fluid diffusion view when solving a linear fixed point problem. This interpretation allows one to better understand why power iteration or Jacobi iteration may converge faster or slower than Gauss-Seidel iteration.


title: On Bayesian supremum norm contraction rates

abstract: Building on ideas from Castillo and Nickl [Ann. Statist. 41 (2013) 1999-2028], a method is provided to study nonparametric Bayesian posterior convergence rates when "strong" measures of distances, such as the sup-norm, are considered. In particular, we show that likelihood methods can achieve optimal minimax sup-norm rates in density estimation on the unit interval. The introduced methodology is used to prove that commonly used families of prior distributions on densities, namely log-density priors and dyadic random density histograms, can indeed achieve optimal sup-norm rates of convergence. New results are also derived in the Gaussian white noise model as a further illustration of the presented techniques.


title: Fractal properties of Bessel functions

abstract: A fractal oscillatority of solutions of second-order differential equations near infinity is measured by oscillatory and phase dimensions. The phase dimension is defined as a box dimension of the trajectory $(x,\dot{x})$ in $\mathbb{R}^2$ of a solution $x=x(t)$, assuming that $(x,\dot{x})$ is a spiral converging to the origin. In this work, we study the phase dimension of the class of second-order nonautonomous differential equations with oscillatory solutions including the Bessel equation. We prove that the phase dimension of Bessel functions is equal to 4/3, and that the corresponding trajectory is a wavy spiral, exhibiting an interesting behavior. The phase dimension of a generalization of the Bessel equation has been also computed.


title: Gravitation at Very Short Distances

abstract: In this note, it is discussed why attempts to renormalize conventional quantum gravity may be based on a false presumption, because, according to Verlinde's entropic view of gravity, the gravitational force between two elementary particles vanishes.


title: Novel distance indicator for gamma-ray bursts associated with supernovae

abstract: Context. It has been proposed that the temporal coincidence of a gamma-ray burst (GRB) and a type Ib/c supernova (SN) can be explained with the concept of induced gravitational collapse (IGC), induced by the matter ejected from an SN Ib/c accreting onto a neutron star (NS). The NS is expected to reach the critical mass necessary for it to collapse to a black hole (BH) and emit a GRB. We found a standard luminosity light curve behavior in the late-time X-ray emission of this subclass of GRBs. Aims. We test if this standard behavior in the luminosity found in this subclass of GRBs can become a redshift estimator of these sources. Methods. We selected a sample of GRBs that belong to this subclass of IGC GRBs associated to an SN (IGC GRB-SN sources). These sources have an isotropic energy $E_{iso} > 10^{52}$ erg and their cosmological redshifts are in the range of $z = 0.49-1.261$. We focused on the corresponding X-ray luminosity light curves. Results. We find that all GRBs of the sample with measured redshift present a standard luminosity late-time light curve in the $0.3-10$ keV rest-frame energy range. We used these results to estimate the GRB redshift of the sample without a measured redshift and found results consistent with other possible redshift indicators. Conclusions. The standard late-time X-ray luminosity light curve of all GRBs of the sample shows a common physical mechanism in this particular phase of the X-ray emission, possibly related to the creation of the NS from the SN process. This scaling law moreover represents strong evidence of very low or even absent beaming in this late phase of the X-ray afterglow emission process. This could be a fundamental tool for estimating the redshift of GRBs that belong to this subclass of events. We are currently expanding this subclass of GRBs to further verify the universal validity of this new redshift estimation method.


title: Strongly residual coordinates over A[x]

abstract: For a domain A of characteristic zero, a polynomial f over A[x] is called a strongly residual coordinate if f becomes a coordinate (over A) upon going modulo x, and f becomes a coordinate upon inverting x. We study the question of when a strongly residual coordinate is a coordinate, a question closely related to the Dolgachev-Weisfeiler conjecture. It is known that all strongly residual coordinates are coordinates for n=2 . We show that a large class of strongly residual coordinates that are generated by elementaries upon inverting x are in fact coordinates for arbitrary n, with a stronger result in the n=3 case. As an application, we show that all Venereau-type polynomials are 1-stable coordinates.


title: On the supernova remnants with flat radio spectra

abstract: A considerable fraction of Galactic supernova remnants (SNRs) characterize flat spectral indices ($\alpha<0.5$). There are several explanations of the flat radio spectra of SNRs in the present literature. The most of models involve a significant contribution of the second-order Fermi mechanism but some of them also discuss high compressions (>4), contribution of secondary electrons left over from the decay of charged pions, as well as the possibility of thermal contamination. In the case of expansion in high density environment, intrinsic thermal bremsstrahlung could theoretically shape the radio spectrum of an SNR and also account for observable curved -- "concave up" radio spectra of some Galactic SNRs. This model could also shed a light on the question of flat spectral indices determined in some Galactic SNRs. On the other hand, present knowledge of the radio continuum spectra (integrated flux densities at different frequencies) of SNRs prevent definite conclusions about the significance of proposed models so the question on flat spectral indices still remains open. New observations, especially at high radio continuum frequencies, are expected to solve these questions in the near future. Finally, as there is a significant connection between the majority of Galactic SNRs with flat integrated radio spectrum and their detection in $\gamma$-rays as well as detection of radiative recombination continua in their X-ray spectra, the analysis of high energy properties of these SNRs is very important.


title: On time and space double-slit experiments

abstract: Time double-slit interference experiments have been achieved and presented as complementary to spatial double-slit interference experiments, providing a further confirmation of the wave-particle duality. Numerical solutions of the free particle time dependent Schr\"odinger equation were presented as explanation of the experimental results, but have been objected to on the basis that the standard non relativistic quantum theory does not have the property of coherence in time. In this note the theoretical and experimental results are derived in a schematic but analytic solution of the TDSE with appropiate initial boundary conditions. The particular boundary conditions are justified by the experimental setups that actually result in having only a single electron at any given time in the double-slit arrangement; and consequently achieve the construction of double peak single electron wave packets. The progressive complementarity of "which-path" ("which-time") information and "space interference" ("oscillating time transient") pattern build up is also exhibited.


title: Tidal turbine array optimisation using the adjoint approach

abstract: Oceanic tides have the potential to yield a vast amount of renewable energy. Tidal stream generators are one of the key technologies for extracting and harnessing this potential. In order to extract an economically useful amount of power, hundreds of tidal turbines must typically be deployed in an array. This naturally leads to the question of how these turbines should be configured to extract the maximum possible power: the positioning and the individual tuning of the turbines could significantly influence the extracted power, and hence is of major economic interest. However, manual optimisation is difficult due to legal site constraints, nonlinear interactions of the turbine wakes, and the cubic dependence of the power on the flow speed. The novel contribution of this paper is the formulation of this problem as an optimisation problem constrained by a physical model, which is then solved using an efficient gradient-based optimisation algorithm. In each optimisation iteration, a two-dimensional finite element shallow water model predicts the flow and the performance of the current array configuration. The gradient of the power extracted with respect to the turbine positions and their tuning parameters is then computed in a fraction of the time taken for a flow solution by solving the associated adjoint equations. These equations propagate causality backwards through the computation, from the power extracted back to the turbine positions and the tuning parameters. This yields the gradient at a cost almost independent of the number of turbines, which is crucial for any practical application. The utility of the approach is demonstrated by optimising turbine arrays in four idealised scenarios and a more realistic case with up to 256 turbines in the Inner Sound of the Pentland Firth, Scotland.


title: Nonlocal Gravity: Damping of Linearized Gravitational Waves

abstract: In nonlocal general relativity, linearized gravitational waves are damped as they propagate from the source to the receiver in the Minkowski vacuum. Nonlocal gravity is a generalization of Einstein's theory of gravitation in which nonlocality is due to the gravitational memory of past events. That nonlocal gravity is dissipative is demonstrated in this paper within certain approximation schemes. The gravitational memory drag leads to the decay of the amplitude of gravitational waves given by the exponential damping factor exp (-t/\tau), where $\tau$ depends on the kernel of nonlocal gravity. The damping time $\tau$ is estimated for gravitational waves of current observational interest and is found to be of the order of, or longer than, the age of the universe.


title: The classification of compact simply connected biquotients in dimension 4 and 5

abstract: We classify all compact simply connected biquotients of dimension 4 and 5. In particular, all pairs of groups $(G,H)$ and embeddings $H\rightarrow G\times G$ giving rise to a particular biquotient are classified.


title: Cabinet of curiosities: the interesting geometry of the angle {\beta} = arccos((3{\phi} - 1)/4)

abstract: In this paper we present the construction of several aggregates of tetrahedra. Each construction is obtained by performing rotations on an initial set of tetrahedra that either (1) contains gaps between adjacent tetrahedra, or (2) exhibits an aperiodic nature. Following this rotation, gaps of the former case are "closed" (in the sense that faces of adjacent tetrahedra are brought into contact to form a "face junction") while translational and rotational symmetries are obtained in the latter case. In all cases, an angular displacement of {\beta} = arccos((3{\phi} - 1)/4) (or a closely related angle), where {\phi} is the golden ratio, is observed between faces of a junction. Additionally, the overall number of plane classes, defined as the number of distinct facial orientations in the collection of tetrahedra, is reduced following the transformation. Finally, we present several "curiosities" involving the structures discussed here with the goal of inspiring the reader's interest in constructions of this nature and their attending, interesting properties.


title: Some algebraic identities for the alpha-permanent

abstract: We show that the permanent of a matrix is a linear combination of determinants of block diagonal matrices which are simple functions of the original matrix. To prove this, we first show a more general identity involving \alpha-permanents: for arbitrary complex numbers \alpha and \beta, we show that the \alpha-permanent of any matrix can be expressed as a linear combination of \beta-permanents of related matrices. Some other identities for the \alpha-permanent of sums and products of matrices are shown, as well as a relationship between the \alpha-permanent and general immanants. We conclude with a discussion of the computational complexity of the \alpha-permanent and provide some numerical illustrations.


title: Minimal surfaces in finite volume hyperbolic 3-manifolds N and in MxS(1), M a finite area hyperbolic surface

abstract: We consider properly immersed finite topology minimal surfaces S in complete finite volume hyperbolic 3-manifolds N, and in M x S(1), where M is a complete hyperbolic surface of finite area. We prove S has finite total curvature equal to 2\pi times the Euler characteristic of S, and we describe the geometry of the ends of S.


title: Brinicles as a case of inverse chemical gardens

abstract: Brinicles are hollow tubes of ice from centimetres to metres in length that form under floating sea ice in the polar oceans when dense, cold brine drains downwards from sea ice into sea water close to its freezing point. When this extremely cold brine leaves the ice it freezes the water it comes into contact with; a hollow tube of ice --- a brinicle --- growing downwards around the plume of descending brine. We show that brinicles can be understood as a form of the self-assembled tubular precipitation structures termed chemical gardens, plant-like structures formed on placing together a soluble metal salt, often in the form of a seed crystal, and an aqueous solution of one of many anions, often silicate. On one hand, in the case of classical chemical gardens, an osmotic pressure difference across a semipermeable precipitation membrane that filters solutions by rejecting the solute leads to an inflow of water and to its rupture. The internal solution, generally being lighter than the external solution, flows up through the break, and as it does so a tube grows upwards by precipitation around the jet of internal solution. Such chemical-garden tubes can grow to many centimetres in length. In the case of brinicles, on the other hand, in floating sea ice we have porous ice in a mushy layer that filters out water, by freezing it, and allows concentrated brine through. Again there is an osmotic pressure difference leading to a continuing ingress of sea water in a siphon pump mechanism that is sustained as long as the ice continues to freeze. Since the brine that is pumped out is denser than the sea water, and descends rather rises, a brinicle is a downwards growing tube of ice; an inverse chemical garden.


title: Parameter estimation for compact binary coalescence signals with the first generation gravitational-wave detector network

abstract: Compact binary systems with neutron stars or black holes are one of the most promising sources for ground-based gravitational wave detectors. Gravitational radiation encodes rich information about source physics; thus parameter estimation and model selection are crucial analysis steps for any detection candidate events. Detailed models of the anticipated waveforms enable inference on several parameters, such as component masses, spins, sky location and distance that are essential for new astrophysical studies of these sources. However, accurate measurements of these parameters and discrimination of models describing the underlying physics are complicated by artifacts in the data, uncertainties in the waveform models and in the calibration of the detectors. Here we report such measurements on a selection of simulated signals added either in hardware or software to the data collected by the two LIGO instruments and the Virgo detector during their most recent joint science run, including a "blind injection" where the signal was not initially revealed to the collaboration. We exemplify the ability to extract information about the source physics on signals that cover the neutron star and black hole parameter space over the individual mass range 1 Msun - 25 Msun and the full range of spin parameters. The cases reported in this study provide a snap-shot of the status of parameter estimation in preparation for the operation of advanced detectors.


title: Local discrete velocity grids for deterministic rarefied flow simulations

abstract: Most of numerical methods for deterministic simulations of rarefied gas flows use the discrete velocity (or discrete ordinate) approximation. In this approach, the kinetic equation is approximated with a global velocity grid. The grid must be large and fine enough to capture all the distribution functions, which is very expensive for high speed flows (like in hypersonic aerodynamics). In this article, we propose to use instead different velocity grids that are local in time and space: these grids dynamically adapt to the width of the distribution functions. The advantages and drawbacks of the method are illustrated in several 1D test cases.


title: High-resolution x-ray analysis with multilayer gratings

abstract: Periodic multilayers are nowadays widely used to perform x-ray analysis in the soft x-ray range (photon energy lower than 1 keV). However, they do not permit to obtain high-resolution spectra like natural or synthetic crystals. Thus, multilayers cannot resolve interferences between close x-ray lines. It has been shown and demonstrated experimentally that patterning a grating profile within a multilayer structure leads to a diffractive optics with improved resolving power. We illustrate the use of a Mo/B4C multilayer grating in the Fe L and C K spectral ranges, around 700 eV and 280 eV respectively. First, in the Fe L range, the improved spectral resolution enables us to distinguish the Fe L\alpha and L\beta emissions (separated by 13 eV). In addition, using a sample made of a mix of LiF and an iron ore, we show that it is possible to easily resolve the F K and Fe L emissions. These examples demonstrate that an improved x-ray analysis can be obtained with multilayer gratings when there is the need to study samples having elements giving rise to close emission lines. Second, in the C K range, by comparing C K\alpha spectra from B4C and cellulose, we show that the shape of the emission band is sensitive to the chemical state of the carbon atom.


title: Sampling from Dirichlet process mixture models with unknown concentration parameter: Mixing issues in large data implementations

abstract: We consider the question of Markov chain Monte Carlo sampling from a general stick-breaking Dirichlet process mixture model, with concentration parameter alpha. This paper introduces a Gibbs sampling algorithm that combines the slice sampling approach of Walker (2007) and the retrospective sampling approach of Papaspiliopoulos and Roberts (2008). Our general algorithm is implemented as efficient open source C++ software, available as an R package, and is based on a blocking strategy similar to that suggested by Papaspiliopoulos (2008) and implemented by Yau et al (2011). We discuss the difficulties of achieving good mixing in MCMC samplers of this nature and investigate sensitivity to initialisation. We additionally consider the challenges when an additional layer of hierarchy is added such that joint inference is to be made on alpha. We introduce a new label switching move and compute the marginal model posterior to help to surmount these difficulties. Our work is illustrated using a profile regression (Molitor et al, 2010) application, where we demonstrate good mixing behaviour for both synthetic and real examples.


title: Hitting time theorems for random matrices

abstract: Starting from an n-by-n matrix of zeros, choose uniformly random zero entries and change them to ones, one-at-a-time, until the matrix becomes invertible. We show that with probability tending to one as n tends to infinity, this occurs at the very moment the last zero row or zero column disappears. We prove a related result for random symmetric Bernoulli matrices, and give quantitative bounds for some related problems. These results extend earlier work by Costello and Vu [arXiv:math/0606414].


title: Equivalence of two definitions of the effective mass of a polaron

abstract: Two definitions of the effective mass of a particle interacting with a quantum field, such as a polaron, are considered and shown to be equal in models similar to the Froehlich polaron model. These are: 1. the mass defined by the low momentum energy $E(P) \approx E(0)+P^2/2M$ of the translation invariant system constrained to have momentum $P$ and 2. the mass $M$ of a simple particle in an arbitrary slowly varying external potential, $V$, described by the nonrelativistic Schroedinger equation, whose ground state energy equals that of the combined particle/field system in a bound state in the same $V$.


title: Characterization of intra device mutual thermal coupling in multi finger SiGe:C HBTs

abstract: This paper studies the mutual coupling in trench isolated multi emitter bipolar transistors fabricated in a Si/SiGe:C HBT technology STMicroelectronics featuring fT and fmax of ~300GHz and ~400GHz, respectively. Thermal coupling parameters are extracted using three dimensional (3D) thermal TCAD simulations. The obtained parameters are implemented in a distributed transistor model that considers self-heating as well as thermal coupling between emitter fingers. Very good agreement is achieved between circuit simulations and DC measurements carried out on an in house designed test structure.


title: Arbitrarily Large Residual Finiteness Growth

abstract: The residual finiteness growth of a group quantifies how well approximated the group is by its finite quotients. In this paper, we construct groups with arbitrarily large residual finiteness growth. We also demonstrate a new relationship between residual finiteness growth and some decision problems in groups, which we apply to our new groups.


title: A convolution method for numerical solution of backward stochastic differential equations

abstract: We propose a new method for the numerical solution of backward stochastic differential equations (BSDEs) which finds its roots in Fourier analysis. The method consists of an Euler time discretization of the BSDE with certain conditional expectations expressed in terms of Fourier transforms and computed using the fast Fourier transform (FFT). The problem of error control is addressed and a local error analysis is provided. We consider the extension of the method to forward-backward stochastic differential equations (FBSDEs) and reflected FBSDEs. Numerical examples are considered from finance demonstrating the performance of the method.


title: Triboelectric charging of volcanic ash from the 2011 Gr\'{i}msv\"{o}tn eruption

abstract: The plume from the 2011 eruption of Gr\'{i}msv\"{o}tn was highly electrically charged, as shown by the considerable lightning activity measured by the UK Met Office's low-frequency lightning detection network. Previous measurements of volcanic plumes have shown that ash particles are electrically charged up to hundreds of km away from the vent, which indicates that the ash continues to charge in the plume [Harrison et al., Env. Res. Lett. 5 024004 (2010), Hatakeyama J. Met. Soc. Japan 27 372 (1949)]. In this paper we study triboelectric charging of different size fractions of a sample of volcanic ash experimentally. Consistently with previous work, we find that the particle size distribution is a determining factor in the charging. Specifically, our laboratory experiments demonstrate that the normalised span of the particle size distribution plays an important role in the magnitude of charging generated. The influence of the normalised span on plume charging suggests that all ash plumes are likely to be charged, with implications for remote sensing and plume lifetime through scavenging effects.


title: Capacity Considerations for Secondary Networks in TV White Space

abstract: The so-called `TV white spaces' (TVWS) - representing unused TV channels in any given location as the result of the transition to digital broadcasting - designated by U.S. Federal Communications Commission (FCC) for unlicensed use presents significant new opportunities within the context of emerging 4G networks for developing new wireless access technologies that meet the goals of the US National Broadband Plan (notably true broadband access for an increasing fraction of the population). There are multiple challenges in realizing this goal; the most fundamental being the fact that the available WS capacity is currently not accurately known, since it depends on a multiplicity of factors - including system parameters of existing incumbents (broadcasters), propagation characteristics of local terrain as well as FCC rules. In this paper, we explore the capacity of white space networks by developing a detailed model that includes all the major variables, and is cognizant of FCC regulations that provide constraints on incumbent protection. Real terrain information and propagation models for the primary broadcaster and adjacent channel interference from TV transmitters are included to estimate their impact on achievable WS capacity. The model is later used to explore various trade-offs between network capacity and system parameters and suggest possible amendments to FCC's incumbent protection rules in the favor of furthering white space capacity.


title: The Variable Optical Polarization of J0948+0022

abstract: The results of monitoring the optical polarization and the gamma-ray flux for the radio-loud narrow-line Seyfert 1 galaxy is reported. We have detected a weak correlation between the gamma-ray variability and the optical variability. However, these observations do not provide conclusive evidence that the emissions are due to a relativistic jet oriented close to the line-of-sight to the observer.


title: The MSSM Higgs sector at a high $M_{SUSY}$: reopening the low tan$\beta$ regime and heavy Higgs searches

abstract: One of the main implications of the LHC discovery of a Higgs boson with a mass $M_h \approx 126$ GeV is that the scale of supersymmetry-breaking in the Minimal Supersymmetric Standard Model (MSSM) might be rather high, $M_S \gg M_Z$. In this paper, we consider the high $M_S$ regime and study the spectrum of the extended Higgs sector of the MSSM, including the LHC constraints on the mass and the rates of the observed light $h$ state. In particular, we show that in a simplified model that approximates the important radiative corrections, the unknown scale $M_S$ (and some other leading SUSY parameters) can be traded against the measured value of $M_h$. One would be then essentially left with only two free parameters to describe the Higgs sector, tan$\beta$ and the pseudoscalar Higgs mass $M_A$, even at higher orders. The main phenomenological consequence of these high $M_S$ values is to reopen the low tan$\beta$ region, tan$\beta \lsim 3-5$, which was for a long time buried under the LEP constraint on the lightest $h$ mass when a low SUSY scale was assumed. We show that, in this case, the heavier MSSM neutral $H/A$ and charged $H^\pm$ states can be searched for in a variety of interesting final states such as decays into gauge and lighter Higgs bosons (in pairs on in mixed states) and decays into heavy top quarks. Examples of sensitivity on the $[tan\beta, M_A]$ parameter space at the LHC in these channels are given.


title: Unimodularity and preservation of volumes in nonholonomic mechanics

abstract: The equations of motion of a mechanical system subjected to nonholonomic linear constraints can be formulated in terms of a linear almost Poisson structure in a vector bundle. We study the existence of invariant measures for the system in terms of the unimodularity of this structure. In the presence of symmetries, our approach allows us to give necessary and sufficient conditions for the existence of an invariant volume, that unify and improve results existing in the literature. We present an algorithm to study the existence of a smooth invariant volume for nonholonomic mechanical systems with symmetry and we apply it to several concrete mechanical examples.


title: Effective Field Theory of Precision Electroweak Physics at One Loop

abstract: The one loop effects of two dimension-six operators on gauge boson self energies are computed within an effective field theory framework. These self energies are translated into effects on precision electroweak observables, and bounds are obtained on the operator coefficients. The effective field theory framework allows for the divergences that arise in the loop calculations to be properly handled, and for unambiguous bounds on the coefficients to be obtained. We find that the coefficients are only weakly bounded, in contrast to previous calculations that obtained much stronger bounds. We argue that the results of these previous calculations are specious.


title: Upgraded Approximation of Non-Binary Alphabets for Polar Code Construction

abstract: An algorithm is presented for approximating a single-user channel with a prime input alphabet size. The result is an upgraded version of the channel with a reduced output alphabet size. It is shown that this algorithm can be used to reduce the output alphabet size to the input alphabet size in most cases.


title: Probing the Pulsar Origin of the Anomalous Positron Fraction with AMS-02 and Atmospheric Cherenkov Telescopes

abstract: Recent observations by PAMELA, Fermi-LAT, and AMS-02 have conclusively indicated a rise in the cosmic-ray positron fraction above 10 GeV, a feature which is impossible to mimic under the paradigm of secondary positron production with self-consistent Galactic cosmic-ray propagation models. A leading explanation for the rising positron fraction is an additional source of electron-positron pairs, for example one or more mature, energetic, and relatively nearby pulsars. We point out that any one of two well-known nearby pulsars, Geminga and Monogem, can satisfactorily provide enough positrons to reproduce AMS-02 observations. A smoking-gun signature of this scenario is an anisotropy in the arrival direction of the cosmic-ray electrons and positrons, which may be detectable by existing, or future, telescopes. The predicted anisotropy level is, at present, consistent with limits from Fermi-LAT and AMS-02. We argue that the large collecting area of Atmospheric Cherenkov Telescopes (ACTs) makes them optimal tools for detecting such an anisotropy. Specifically, we show that much of the proton and gamma-ray background, which affects measurements of the cosmic-ray electron-positron spectrum with ACTs, may be controlled in the search for anisotropies. We conclude that observations using archival ACT data could already constrain or substantiate the pulsar origin of the positron anomaly, while upcoming instruments (such as the Cherenkov Telescope Array) will provide strong constraints on the source of the rising positron fraction.


title: Properties of Two-Dimensional Silicon grown on Graphene Substrate

abstract: The structure and electrical properties of a two-dimensional (2D) sheet of silicon on a graphene substrate are studied using first-principles calculations. A new corrugated rectangular structure of silicon is proposed to be the most energetically favorable structure. The shifting of the Fermi energy level indicates self-doping. Calculation of electron density shows a weak coupling between the silicon layer and graphene substrate. The 2D silicon sheet turns to be metallic and has a much higher value of transmission efficiency (TE) than the underlying graphene substrate.


title: Thermotropic nematic order upon nano-capillary filling

abstract: Optical birefringence and light absorption measurements reveal four regimes for the thermotropic behavior of a nematogen liquid (7CB) upon sequential filling of parallel-aligned capillaries of 12 nm diameter in a monolithic, mesoporous silica membrane. No molecular reorientation is observed for the first adsorbed monolayer. In the film-condensed state (up to 1 nm thickness) a weak, continuous paranematic-to-nematic (P-N) transition is found, which is shifted by 10 K below the discontinuous bulk transition at T_IN=305K. The capillary-condensed state exhibits a more pronounced, albeit still continuous P-N reordering, located 4 K below T_IN. This shift vanishes abruptly on complete filling of the capillaries, which we tentatively trace to a 10 MPa tensile pressure release associated with the disappearance of concave menisci in the confined liquid.

categories: cond-mat.soft cond-mat.mes-hall cond-mat.mtrl-sci physics.chem-ph


title: Is every matrix similar to a polynomial in a companion matrix?

abstract: Given a field $F$, an integer $n\geq 1$, and a matrix $A\in M_n(F)$, are there polynomials $f,g\in F[X]$, with $f$ monic of degree $n$, such that $A$ is similar to $g(C_f)$, where $C_f$ is the companion matrix of $f$? For infinite fields the answer is easily seen to positive, so we concentrate on finite fields. In this case we give an affirmative answer, provided $|F|\geq n-2$. Moreover, for any finite field $F$, with $|F|=m$, we construct a matrix $A\in M_{m+3}(F)$ that is not similar to any matrix of the form $g(C_f)$. Of use above, but also of independent interest, is a constructive procedure to determine the similarity type of any given matrix $g(C_f)$ purely in terms of $f$ and $g$, without resorting to polynomial roots in $F$ or in any extension thereof. This, in turn, yields an algorithm that, given $g$ and the invariant factors of any $A$, returns the elementary divisors of $g(A)$. It is a rational procedure, as opposed to the classical method that uses the Jordan decomposition of $A$ to find that of $g(A)$. Finally, extending prior results by the authors, we show that for an integrally closed ring $R$ with field of fractions $F$ and companion matrices $C,D$ the subalgebra $R< C,D>$ of $M_n(R)$ is a free $R$-module of rank $n+(n-m)(n-1)$, where $m$ is the degree of $\gcd (f,g)\in F[X]$, and a presentation for $R< C,D>$ is given in terms of $C$ and $D$. A counterexample is furnished to show that $R< C,D>$ need not be a free $R$-module if $R$ is not integrally closed. The preceding information is used to study $M_n(R)$, and others, as $R[X]$-modules.


title: Dynamics of a Stabilized Radion and Duality

abstract: We construct the effective theory of the graviscalar radion in the Randall-Sundrum scenario, taking into account effects arising from the stabilization of the extra dimension through the Goldberger-Wise mechanism. We explore the conditions under which the radion can remain light, and determine the corrections to its couplings to Standard Model (SM) states when the effects of stabilization are taken into account. We show that in the theories of interest for electroweak symmetry breaking that have a holographic dual, the presence of a light radion in the spectrum is not a robust prediction of the framework, but is in fact associated with mild tuning. We find that corrections to the form of the radion couplings to Standard Model particles arising from effects associated with brane stabilization are suppressed by the square of the ratio of the radion mass to the Kaluza-Klein scale. These corrections are small in theories where the radion is light, and are generally subleading, except in the case of couplings to the SM gluons and photon, when they can sometimes dominate. The AdS/CFT correspondence relates the radion in Randall-Sundrum models to the dilaton in theories where a strongly coupled conformal symmetry is spontaneously broken. We show that the discrepancies in the literature between the results for the dilaton and the radion can be traced to the omission of self-interaction terms that would otherwise dominate the potential for the Goldberger-Wise scalar near the infrared brane. In the dual picture, this corresponds to neglecting the corrections to the scaling behavior of the operator that breaks conformal symmetry when it grows large. With the inclusion of these self-interaction terms, we find good agreement between the results on the two sides of the correspondence.


title: The Round Complexity of Small Set Intersection

abstract: The set disjointness problem is one of the most fundamental and well-studied problems in communication complexity. In this problem Alice and Bob hold sets $S, T \subseteq [n]$, respectively, and the goal is to decide if $S \cap T = \emptyset$. Reductions from set disjointness are a canonical way of proving lower bounds in data stream algorithms, data structures, and distributed computation. In these applications, often the set sizes $|S|$ and $|T|$ are bounded by a value $k$ which is much smaller than $n$. This is referred to as small set disjointness. A major restriction in the above applications is the number of rounds that the protocol can make, which, e.g., translates to the number of passes in streaming applications. A fundamental question is thus in understanding the round complexity of the small set disjointness problem. For an essentially equivalent problem, called OR-Equality, Brody et. al showed that with $r$ rounds of communication, the randomized communication complexity is $\Omega(k \ilog^r k)$, where$\ilog^r k$ denotes the $r$-th iterated logarithm function. Unfortunately their result requires the error probability of the protocol to be $1/k^{\Theta(1)}$. Since na\"ive amplification of the success probability of a protocol from constant to $1-1/k^{\Theta(1)}$ blows up the communication by a $\Theta(\log k)$ factor, this destroys their improvements over the well-known lower bound of $\Omega(k)$ which holds for any number of rounds. They pose it as an open question to achieve the same $\Omega(k \ilog^r k)$ lower bound for protocols with constant error probability. We answer this open question by showing that the $r$-round randomized communication complexity of ${\sf OREQ}_{n,k}$, and thus also of small set disjointness, with {\it constant error probability} is $\Omega(k \ilog^r k)$, asymptotically matching known upper bounds for ${\sf OREQ}_{n,k}$ and small set disjointness.


title: On magnetic-field-induced dissipationless electric current in helicoidal graphene nanoribbons

abstract: We argue that twisted (helicoidal) graphene nanoribbons may support dissipationless electric current in the presence of static uniform magnetic field. The non-resistive charge transfer in this parity-odd system should be enhanced by thermal fluctuations.


title: Two-dimensional S-matrices from unitarity cuts

abstract: Using unitarity methods, we compute, for several massive two-dimensional models, the cut-constructible part of the one-loop 2->2 scattering S-matrices from the tree-level amplitudes. We apply our method to various integrable theories, finding evidence that for supersymmetric models the one-loop S-matrix is cut-constructible, while for models without supersymmetry (but with integrability) the missing rational terms are proportional to the tree-level S-matrix and therefore can be interpreted as a shift in the coupling. Finally, applying our procedure to the world-sheet theory for the light-cone gauge-fixed AdS_5 x S^5 superstring we reproduce, at one-loop in the near-BMN expansion, the S-matrix known from integrability techniques.


title: Intrinsic and extrinsic mirror symmetry breaking in anti-dot spin-wave waveguides

abstract: We theoretically study the spin-wave spectra in magnonic waveguides periodically patterned with square anti-dots in nanoscale with pinned magnetization at the edges. We show that the breaking of the mirror symmetry of the waveguide by the structural changes can result in a magnonic band gap closing. These intrinsic symmetry breaking can be compensated by properly chosen asymmetric external bias magnetic field, i.e. in an extrinsic way. As a result the magnonic gaps existing in the ideal symmetric structure can be recovered. The model used for the explanation also suggests that this idea could be generic both for exchange and dipolar interaction regimes of spin-waves and also for other types of waves, e.g., electrons in the graphene ribbons.


title: Effective Hamiltonian for the hybrid double quantum dot qubit

abstract: Quantum dot hybrid qubits formed from three electrons in double quantum dots represent a promising compromise between high speed and simple fabrication for solid state implementations of single qubit and two qubits quantum logic ports. We derive the Schrieffer-Wolff effective Hamiltonian that describes in a simple and intuitive way the qubit by combining a Hubbard-like model with a projector operator method. As a result, the Hubbard-like Hamiltonian is transformed in an equivalent expression in terms of the exchange coupling interactions between pairs of electrons. The effective Hamiltonian is exploited to derive the dynamical behaviour of the system and its eigenstates on the Bloch sphere to generate qubits operation for quantum logic ports. A realistic implementation in silicon and the coupling of the qubit with a detector are discussed.


title: The PdBI Arcsecond Whirlpool Survey (PAWS). I. A Cloud-Scale/Multi-Wavelength View of the Interstellar Medium in a Grand-Design Spiral Galaxy

abstract: The PdBI (Plateau de Bure Interferometer) Arcsecond Whirlpool Survey (PAWS) has mapped the molecular gas in the central ~9kpc of M51 in its 12CO(1-0) line emission at cloud-scale resolution of ~40pc using both IRAM telescopes. We utilize this dataset to quantitatively characterize the relation of molecular gas (or CO emission) to other tracers of the interstellar medium (ISM), star formation and stellar populations of varying ages. Using 2-dimensional maps, a polar cross-correlation technique and pixel-by-pixel diagrams, we find: (a) that (as expected) the distribution of the molecular gas can be linked to different components of the gravitational potential, (b) evidence for a physical link between CO line emission and radio continuum that seems not to be caused by massive stars, but rather depend on the gas density, (c) a close spatial relation between the PAH and molecular gas emission, but no predictive power of PAH emission for the molecular gas mass,(d) that the I-H color map is an excellent predictor of the distribution (and to a lesser degree the brightness) of CO emission, and (e) that the impact of massive (UV-intense) young star-forming regions on the bulk of the molecular gas in central ~9kpc can not be significant due to a complex spatial relation between molecular gas and star-forming regions that ranges from co-spatial to spatially offset to absent. The last point, in particular, highlights the importance of galactic environment -- and thus the underlying gravitational potential -- for the distribution of molecular gas and star formation.


title: Cosmic Chronometers in the R_h=ct Universe

abstract: The use of luminous red galaxies as cosmic chronometers provides us with an indispensable method of measuring the universal expansion rate H(z) in a model-independent way. Unlike many probes of the cosmological history, this approach does not rely on integrated quantities, such as the luminosity distance, and therefore does not require the pre-assumption of any particular model, which may bias subsequent interpretations of the data. We employ three statistical tools -- the Akaike, Kullback, and Bayes Information Criteria (AIC, KIC and BIC) -- to compare the LCDM model and the R_h=ct Universe with the currently available measurements of H(z), and show that the R_h=ct Universe is favored by these model selection criteria. The parameters in each model are individually optimized by maximum likelihood estimation. The R_h=ct Universe fits the data with a reduced chi^2_dof=0.745 for a Hubble constant H_0=63.2+/-2.5 km/s/Mpc, and H_0 is the sole parameter in this model. By comparison, the optimal LCDM model, which has three free parameters (including H_0=68.9+/-2.4 km/s/Mpc, Omega_m=0.32, and a dark-energy equation of state p_de=-rho_de), fits the H(z) data with a reduced chi^2_dof=0.777. With these chi^2_dof values, the AIC yields a likelihood of about 82 per cent that the distance--redshift relation of the R_h=ct Universe is closer to the correct cosmology, than is the case for LCDM. If the alternative BIC criterion is used, the respective Bayesian posterior probabilities are 91.2 per cent (R_h=ct) versus 8.8 per cent (LCDM). Using the concordance LCDM parameter values, rather than those obtained by fitting LCDM to the cosmic chronometer data, would further disfavor LCDM.


title: The N=4 Superconformal Bootstrap

abstract: We implement the conformal bootstrap for N=4 superconformal field theories in four dimensions. Consistency of the four-point function of the stress-energy tensor multiplet imposes significant upper bounds for the scaling dimensions of unprotected local operators as functions of the central charge of the theory. At the threshold of exclusion, a particular operator spectrum appears to be singled out by the bootstrap constraints. For large values of the central charge, this extremal spectrum is compatible with that of supergravity in AdS_5 x S^5. For finite central charge, we conjecture that the extremal spectrum is that of N=4 SYM at an S-duality invariant value of the complexified gauge coupling.


title: The Cosmophenomenology of Axionic Dark Radiation

abstract: Relativistic axions are good candidates for the dark radiation for which there are mounting observational hints. The primordial decays of heavy fields produce axions which are ultra-energetic compared to thermalised matter and inelastic axion-matter scattering can occur with $E_{CoM} \gg T_{\gamma}$, thus accessing many interesting processes which are otherwise kinematically forbidden in standard cosmology. Axion-photon scattering into quarks and leptons during BBN affects the light element abundances, and bounds on overproduction of $^4$He constrain a combination of the axion decay constant and the reheating temperature. For supersymmetric models, axion scattering into visible sector superpartners can give direct non-thermal production of dark matter at $T_{\gamma} \ll T_{freezeout}$. Most axions --- or any other dark radiation candidate from modulus decay --- still linger today as a Cosmic Axion Background with $E_{axion} \sim \mathcal{O}(100) eV$, and a flux of $\sim 10^6 cm^{-2} s^{-1}$.


title: General relativistic simulations of binary black hole-neutron stars: Precursor electromagnetic signals

abstract: We perform the first general relativistic force-free simulations of neutron star (NS) magnetospheres in orbit about spinning and non-spinning black holes. We find promising precursor electromagnetic emission: typical Poynting luminosities at, e.g., an orbital separation of 6.6 times the NS radius are L ~ 6 x 10^{42} erg/s for a 1.4 solar-mass NS with a 10^{13}G polar magnetic field. The Poynting flux peaks within a broad beam of ~40 degrees in the azimuthal direction and within ~60 degrees from the orbital plane, establishing a possible lighthouse effect. Our calculations, though preliminary, preview more detailed simulations of these systems that we plan to perform in the future.


title: Multiloop integrals in dimensional regularization made simple

abstract: Scattering amplitudes at loop level can be expressed in terms of Feynman integrals. The latter satisfy partial differential equations in the kinematical variables. We argue that a good choice of basis for (multi-)loop integrals can lead to significant simplifications of the differential equations, and propose criteria for finding an optimal basis. This builds on experience obtained in supersymmetric field theories that can be applied successfully to generic quantum field theory integrals. It involves studying leading singularities and explicit integral representations. When the differential equations are cast into canonical form, their solution becomes elementary. The class of functions involved is easily identified, and the solution can be written down to any desired order in epsilon within dimensional regularization. Results obtained in this way are particularly simple and compact. In this letter, we outline the general ideas of the method and apply them to a two-loop example.


title: Entropy production in classical Yang-Mills theory from Glasma initial conditions

abstract: We study the thermalization process in classical Yang-Mills (CYM) field theory starting from noisy glasma-like initial conditions by investigating the initial-value sensitivity of trajectories. Kunihiro et al. linked entropy generation to the Kolmogorov-Sinai entropy, which gives the entropy production rate in classical chaotic systems, calculated numerically for CYM fields starting from purely random initial field configurations. In contrast, we here study glasma-like initial conditions. For small random fluctuations we obtain qualitatively similar results while no entropy increase is observed when such fluctuations are absent. We analyze the intermediate time Lyapunov spectrum for several time windows and calculate the Kolmogorov-Sinai entropy. We find a large number of positive Lyapunov exponents at the early stages of time evolution. Also for later times their number is a sizeable fraction of the total number of degrees of freedom. The spectrum of positive Lyapunov exponents at first changes rapidly, but then stabilizes, indicating that the dynamics of the gauge fields approaches a steady state. Thus we conclude that also for glasma-like initial conditions a significant amount of entropy is produced by classical gluon field dynamics.

categories: hep-ph cond-mat.stat-mech nlin.CD nucl-th


title: The not-so-massive black hole in the microquasar GRS1915+105

abstract: We present a new dynamical study of the black hole X-ray transient GRS1915+105 making use of near-infrared spectroscopy obtained with X-shooter at the VLT. We detect a large number of donor star absorption features across a wide range of wavelengths spanning the H and K bands. Our 24 epochs covering a baseline of over 1 year permit us to determine a new binary ephemeris including a refined orbital period of P=33.85 +/- 0.16 d. The donor star radial velocity curves deliver a significantly improved determination of the donor semi-amplitude which is both accurate (K_2=126 +/- 1 km/s) and robust against choice of donor star template and spectral features used. We furthermore constrain the donor star's rotational broadening to vsini=21 +/-4 km/s, delivering a binary mass ratio of q=0.042 +/- 0.024. If we combine these new constraints with distance and inclination estimates derived from modelling the radio emission, a black hole mass of M_BH=10.1 +/- 0.6 M_sun is inferred, paired with an evolved mass donor of M_2=0.47 +/- 0.27 M_sun. Our analysis suggests a more typical black hole mass for GRS1915+105 rather than the unusually high values derived in the pioneering dynamical study by Greiner et al. (2001). Our data demonstrate that high-resolution infrared spectroscopy of obscured accreting binaries can deliver dynamical mass determinations with a precision on par with optical studies.


title: Heterotic Moduli Stabilisation

abstract: We perform a systematic analysis of moduli stabilisation for weakly coupled heterotic string theory compactified on manifolds which are Calabi-Yau up to alpha' effects. We review how to fix all geometric and bundle moduli in a supersymmetric way by fractional fluxes, the requirement of a holomorphic gauge bundle, D-terms, higher order perturbative contributions to W, non-perturbative and threshold effects. We then show that alpha' corrections to K lead to new stable Minkowski (or dS) vacua where the complex structure moduli Z and the dilaton are fixed supersymmetrically, while the fixing of the Kahler moduli at a lower scale leads to spontaneous SUSY breaking. The minimum lies at moderately large volumes of all geometric moduli, at a perturbative string coupling and at the right value of the GUT coupling. We also give a dynamical derivation of anisotropic compactifications which allow for gauge coupling unification around 10^16 GeV. The gravitino mass can be anywhere between the GUT and TeV scale depending on the fixing of the Z-moduli. In general, these are fixed by turning on background fluxes, leading to a gravitino mass around the GUT scale since the heterotic 3-form flux does not contain enough freedom to tune W to small values. Moreover accommodating the observed value of the cosmological constant (CC) is a challenge. Low-energy SUSY could instead be obtained in particular situations where the gauge bundle is holomorphic only at a point-like sub-locus of Z-moduli space, or where the number of Z-moduli is small (like orbifold models), since in these cases one may fix all moduli without turning on any quantised flux. However tuning the CC is even more of a challenge in these cases. Another option is to focus on non-complex manifolds since these allow for new geometric fluxes which can be used to tune W and the CC, even if their moduli space is presently only poorly understood.


title: A near-optimal approximation algorithm for Asymmetric TSP on embedded graphs

abstract: We present a near-optimal polynomial-time approximation algorithm for the asymmetric traveling salesman problem for graphs of bounded orientable or non-orientable genus. Our algorithm achieves an approximation factor of O(f(g)) on graphs with genus g, where f(n) is the best approximation factor achievable in polynomial time on arbitrary n-vertex graphs. In particular, the O(log(n)/loglog(n))-approximation algorithm for general graphs by Asadpour et al. [SODA 2010] immediately implies an O(log(g)/loglog(g))-approximation algorithm for genus-g graphs. Our result improves the O(sqrt(g)*log(g))-approximation algorithm of Oveis Gharan and Saberi [SODA 2011], which applies only to graphs with orientable genus g; ours is the first approximation algorithm for graphs with bounded non-orientable genus. Moreover, using recent progress on approximating the genus of a graph, our O(log(g) / loglog(g))-approximation can be implemented even without an embedding when the input graph has bounded degree. In contrast, the O(sqrt(g)*log(g))-approximation algorithm of Oveis Gharan and Saberi requires a genus-g embedding as part of the input. Finally, our techniques lead to a O(1)-approximation algorithm for ATSP on graphs of genus g, with running time 2^O(g)*n^O(1).


title: The Topology and Size of the Universe from CMB Temperature and Polarization Data

abstract: We analyze seven year and nine year WMAP temperature maps for signatures of three finite flat topologies M_0=T^3, M_1=T^2 x R^1, and M_2=S^1 x R^2. We use Monte-Carlo simulations with the Feldman-Cousins method to obtain confidence intervals for the size of the topologies considered. We analyze the V, W, and Q frequency bands along with the ILC map and find no significant difference in the results. The 95.5% confidence level lower bound on the size of the topology is 1.5L_0 for M_0, 1.4L_0 for M_1, and 1.1L_0 for M_2, where L_0 is the radius of the last scattering surface. Our results agree very well with the recently released results from the Planck temperature data. We show that the likelihood function is not Gaussian in the size, and therefore simulations are important for obtaining accurate bounds on the size. We then introduce the formalism for including polarization data in the analysis. The improvement that we find from WMAP polarization maps is small because of the high level of instrumental noise, but our forecast for Planck maps shows a much better improvement on the lower bound for L. For the M_0 topology we expect an improvement on the lower bound of L from 1.7L_0 to 1.9L_0 at 95.5% confidence level. Using both polarization and temperature data is important because it tests the hypothesis that deviations in the TT spectrum at small l originate in the primordial perturbation spectrum.


title: Tripartite quantum state violating the hidden influence constraints

abstract: The possibility to explain quantum correlations via (possibly) unknown causal influences propagating gradually and continuously at a finite speed v > c has attracted a lot of attention recently. In particular, it could be shown that this assumption leads to correlations that can be exploited for superluminal communication. This was achieved studying the set of possible correlations that are allowed within such a model and comparing them to correlations produced by local measurements on a four-party entangled quantum state. Here, we report on a quantum state that allows for the same conclusion involving only three parties.


title: Characterization of projectively flat Finsler manifolds of constant curvature with finite dimensional holonomy group

abstract: In this paper we prove that the holonomy group of a simply connected locally projectively flat Finsler manifold of constant curvature is a finite dimensional Lie group if and only if it is flat or it is Riemannian.


title: Dynamics of Large Fragments in the Tail of Active Asteroid P/2010 A2

abstract: We examine the motions of large fragments at the head of the dust tail of active asteroid P/2010 A2. In previous work we showed that these fragments were ejected from the primary nucleus in early 2009, either following a hypervelocity impact or by rotationally induced break-up. Here, we follow their positions through a series of Hubble Space Telescope images taken during the first half of 2010. The orbital evolution of each fragment allows us to constrain its velocity relative to the main nucleus after leaving its sphere of gravitational influence. We find that the fragments constituting a prominent X-shaped tail feature were emitted in a direction opposite to the motion of the asteroid and towards the south of its orbital plane. Derived emission velocities of these primary fragments range between 0.02 and 0.3 m/s, comparable to the ~0.08 m/s gravitational escape speed from the nucleus. Their sizes are on the order of decimeters or larger. We obtain the best fits to our data with ejection velocity vectors lying in a plane that includes the nucleus. This may suggest that the cause of the disruption of P/2010 A2 is rotational break-up.


title: On the S-Euclidean minimum of an ideal class

abstract: We show that the S-Euclidean minimum of an ideal class is a rational number, generalizing a result of Cerri. We also give some corollaries which explain the relationship of our results with Lenstra's notion of a norm-Euclidean ideal class and the conjecture of Barnes and Swinnerton-Dyer on quadratic forms. The proof is self-contained but uses ideas from ergodic theory and topological dynamics, particularly those of Berend.


title: Asymptotically inspired moment-closure approximation for adaptive networks

abstract: Adaptive social networks, in which nodes and network structure co-evolve, are often described using a mean-field system of equations for the density of node and link types. These equations constitute an open system due to dependence on higher order topological structures. We propose a new approach to moment closure based on the analytical description of the system in an asymptotic regime. We apply the proposed approach to two examples of adaptive networks: recruitment to a cause model and adaptive epidemic model. We show a good agreement between the improved mean-field prediction and simulations of the full network system.


title: Non-Gaussian Halo Bias Beyond the Squeezed Limit

abstract: Primordial non-Gaussianity, in particular the coupling of modes with widely different wavelengths, can have a strong impact on the large-scale clustering of tracers through a scale-dependent bias with respect to matter. We demonstrate that the standard derivation of this non-Gaussian scale-dependent bias is in general valid only in the extreme squeezed limit of the primordial bispectrum, i.e. for clustering over very large scales. We further show how the treatment can be generalized to describe the scale-dependent bias on smaller scales, without making any assumptions on the nature of tracers apart from a dependence on the small-scale fluctuations within a finite region. If the leading scale-dependent bias \Delta b \propto k^{\alpha}, then the first subleading term will scale as k^{\alpha+2}. This correction typically becomes relevant as one considers clustering over scales k >~ 0.01 h Mpc^{-1}.


title: Analytic Comparison between X-ray Fluorescence CT and K-edge CT

abstract: X-ray fluorescence computed tomography (XFCT) and K-edge computed tomography (CT) are two important modalities to quantify a distribution of gold nanoparticles (GNPs) in a small animal for preclinical studies. It is valuable to determine which modality is more efficient for a given application. In this paper, we report a theoretical analysis in terms of signal-to-noise ratio (SNR) for the two modalities, showing that there is a threshold for the GNPs concentration and XFCT has a better SNR than K-edge CT if GNPs concentration is less than this threshold. Numerical simulations are performed and two kinds of phantoms are used to represent multiple concentration levels and feature sizes. Experimental results illustrate that XFCT is superior to K-edge CT when contrast concentration is lower than 0.4% which coincides with the theoretical analysis.


title: Model-based Bayesian Reinforcement Learning for Dialogue Management

abstract: Reinforcement learning methods are increasingly used to optimise dialogue policies from experience. Most current techniques are model-free: they directly estimate the utility of various actions, without explicit model of the interaction dynamics. In this paper, we investigate an alternative strategy grounded in model-based Bayesian reinforcement learning. Bayesian inference is used to maintain a posterior distribution over the model parameters, reflecting the model uncertainty. This parameter distribution is gradually refined as more data is collected and simultaneously used to plan the agent's actions. Within this learning framework, we carried out experiments with two alternative formalisations of the transition model, one encoded with standard multinomial distributions, and one structured with probabilistic rules. We demonstrate the potential of our approach with empirical results on a user simulator constructed from Wizard-of-Oz data in a human-robot interaction scenario. The results illustrate in particular the benefits of capturing prior domain knowledge with high-level rules.


title: Gromov-Hausdorff collapsing of Calabi-Yau manifolds

abstract: This paper is a sequel to arXiv:1108.0967. We further study Gromov-Hausdorff collapsing limits of Ricci-flat K\"ahler metrics on abelian fibered Calabi-Yau manifolds. Firstly, we show that in the same setup as arXiv:1108.0967, if the dimension of the base manifold is one, the limit metric space is homeomorphic to the base manifold. Secondly, if the fibered Calabi-Yau manifolds are Lagrangian fibrations of holomorphic symplectic manifolds, the metrics on the regular parts of the limits are special K\"ahler metrics. By combining these two results, we extend arXiv:math/0008018 to any fibered projective K3 surface without any assumption on the type of singular fibers.


title: Optimal initiation of a GLWB in a variable annuity: no arbitrage approach

abstract: This paper offers a financial economic perspective on the optimal time (and age) at which the owner of a Variable Annuity (VA) policy with a Guaranteed Living Withdrawal Benefit (GLWB) rider should initiate guaranteed lifetime income payments. We abstract from utility, bequest and consumption preference issues by treating the VA as liquid and tradable. This allows us to use an American option pricing framework to derive a so-called optimal initiation region. Our main practical finding is that given current design parameters in which volatility (asset allocation) is restricted to less than 20%, while guaranteed payout rates (GPR) as well as bonus (roll-up) rates are less than 5%, GLWBs that are in-the-money should be turned on by the late 50s and certainly the early 60s. The exception to the rule is when a non-constant GPR is about to increase (soon) to a higher age band, in which case the optimal policy is to wait until the new GPR is hit and then initiate immediately. Also, to offer a different perspective, we invert the model and solve for the bonus (roll-up) rate that is required to justify delaying initiation at any age. We find that the required bonus is quite high and more than what is currently promised by existing products. Our methodology and results should be of interest to researchers as well as to the individuals that collectively have over \$1 USD trillion in aggregate invested in these products. We conclude by suggesting that much of the non-initiation at older age is irrational (which obviously benefits the insurance industry.)


title: Direct constraints on the top-Higgs coupling from the 8 TeV LHC data

abstract: The LHC experiments have analyzed the 7 and 8 TeV LHC data in the main Higgs production and decay modes. Current analyses only loosely constrain an anomalous top-Higgs coupling in a direct way. In order to strongly constrain this coupling, the Higgs-top associated production is reanalyzed. Thanks to the strong destructive interference in the t-channel for standard model couplings, this process can be very sensitive to both the magnitude and the sign of a non-standard top-Higgs coupling. We project the sensitivity to anomalous couplings to the integrated luminosity of 50 fb^{-1}, corresponding to the data collected by the ATLAS and CMS experiments in 7 and 8 TeV collisions, as of 2012. We show that the combination of di-photon and multi-lepton signatures, originating from different combinations of the top and Higgs decay modes, can be a potential probe to constrain a large portion of the negative top-Higgs coupling space presently allowed by the ATLAS and CMS global fits.


title: A scaling approach to Caffarelli-Kohn-Nirenberg inequality

abstract: We mainly consider the general Caffarelli-Kohn-Nirenberg inequality in the Euclidean and Riemannian setting. In both cases, our proof relies mostly on a new parameter s conveniently introduced, see (2.7).


title: Backbone structure of the Edwards-Anderson spin-glass model

abstract: We study the ground-state spatial heterogeneities of the Edwards-Anderson spin-glass model with both bimodal and Gaussian bond distributions. We characterize these heterogeneities by using a general definition of bond rigidity, which allows us to classify the bonds of the system into two sets, the backbone and its complement, with very different properties. This generalizes to continuous distributions of bonds the well known definition of a backbone for discrete bond distributions. By extensive numerical simulations we find that the topological structure of the backbone for a given lattice dimensionality is very similar for both discrete and continuous bond distributions. We then analyze how these heterogeneities influence the equilibrium properties at finite temperature and we discuss the possibility that a suitable backbone picture can be relevant to describe spin-glass phenomena.


title: Three years of Fermi GBM Earth Occultation Monitoring: Observations of Hard X-ray/Soft Gamma-Ray Sources

abstract: The Gamma ray Burst Monitor (GBM) on board Fermi Gamma-ray Space Telescope has been providing continuous data to the astronomical community since 2008 August 12. We will present the results of the analysis of the first three years of these continuous data using the Earth occultation technique to monitor a catalog of 209 sources. Although the occultation technique is in principle quite simple, in practice there are many complications including the dynamic instrument response, source confusion, and scattering in the Earth's atmosphere, which will be described. We detect 99 sources, including 40 low-mass X-ray binary/neutron star systems, 31 high-mass X-ray binary/neutron star systems, 12 black hole binaries, 12 active galaxies, 2 other sources, plus the Crab Nebula and the Sun. Nine of these sources are detected in the 100-300 keV band, including seven black-hole binaries, the active galaxy Cen A, and the Crab. The Crab and Cyg X-1 are also detected in the 300-500 keV band. GBM provides complementary data to other sky monitors below 100 keV and is the only all-sky monitor above 100 keV. In our fourth year of monitoring, we have already increased the number of transient sources detected and expect several of the weaker persistent sources to cross the detection threshold. I will briefly discuss these new sources and what to expect from our five year occultation catalog.


title: Concentration inequalities for non-Lipschitz functions with bounded derivatives of higher order

abstract: Building on the inequalities for homogeneous tetrahedral polynomials in independent Gaussian variables due to R. Lata{\l}a we provide a concentration inequality for non-necessarily Lipschitz functions $f\colon \R^n \to \R$ with bounded derivatives of higher orders, which hold when the underlying measure satisfies a family of Sobolev type inequalities $\|g- \E g\|_p \le C(p)\|\nabla g\|_p.$ Such Sobolev type inequalities hold, e.g., if the underlying measure satisfies the log-Sobolev inequality (in which case $C(p) \le C\sqrt{p}$) or the Poincar\'e inequality (then $C(p) \le Cp$). Our concentration estimates are expressed in terms of tensor-product norms of the derivatives of $f$. When the underlying measure is Gaussian and $f$ is a polynomial (non-necessarily tetrahedral or homogeneous), our estimates can be reversed (up to a constant depending only on the degree of the polynomial). We also show that for polynomial functions, analogous estimates hold for arbitrary random vectors with independent sub-Gaussian coordinates. We apply our inequalities to general additive functionals of random vectors (in particular linear eigenvalue statistics of random matrices) and the problem of counting cycles of fixed length in Erd\H{o}s-R{\'e}nyi random graphs, obtaining new estimates, optimal in a certain range of parameters.


title: Fuzzy Aggregates in Fuzzy Answer Set Programming

abstract: Fuzzy answer set programming is a declarative framework for representing and reasoning about knowledge in fuzzy environments. However, the unavailability of fuzzy aggregates in disjunctive fuzzy logic programs, DFLP, with fuzzy answer set semantics prohibits the natural and concise representation of many interesting problems. In this paper, we extend DFLP to allow arbitrary fuzzy aggregates. We define fuzzy answer set semantics for DFLP with arbitrary fuzzy aggregates including monotone, antimonotone, and nonmonotone fuzzy aggregates. We show that the proposed fuzzy answer set semantics subsumes both the original fuzzy answer set semantics of DFLP and the classical answer set semantics of classical disjunctive logic programs with classical aggregates, and consequently subsumes the classical answer set semantics of classical disjunctive logic programs. We show that the proposed fuzzy answer sets of DFLP with fuzzy aggregates are minimal fuzzy models and hence incomparable, which is an important property for nonmonotonic fuzzy reasoning.


title: Network Compression: Worst-Case Analysis

abstract: We study the problem of communicating a distributed correlated memoryless source over a memoryless network, from source nodes to destination nodes, under quadratic distortion constraints. We establish the following two complementary results: (a) for an arbitrary memoryless network, among all distributed memoryless sources of a given correlation, Gaussian sources are least compressible, that is, they admit the smallest set of achievable distortion tuples, and (b) for any memoryless source to be communicated over a memoryless additive-noise network, among all noise processes of a given correlation, Gaussian noise admits the smallest achievable set of distortion tuples. We establish these results constructively by showing how schemes for the corresponding Gaussian problems can be applied to achieve similar performance for (source or noise) distributions that are not necessarily Gaussian but have the same covariance.


title: On Double 3-Term Arithmetic Progressions

abstract: In this note we are interested in the problem of whether or not every increasing sequence of positive integers $x_1x_2x_3...$ with bounded gaps must contain a double 3-term arithmetic progression, i.e., three terms $x_i$, $x_j$, and $x_k$ such that $i + k = 2j$ and $x_i + x_k = 2x_j$. We consider a few variations of the problem, discuss some related properties of double arithmetic progressions, and present several results obtained by using RamseyScript, a high-level scripting language.


title: Linear magnetoresistance from Dirac-like fermions in graphite

abstract: We show that magnetoresistance of Bernal-stacked graphite (with the magnetic field ${\bf B}$ parallel to the c-axis and the current in the ab plane) scales linearly with the magnetic field over an interval of classically weak fields. The linearity is related to the presence of extremely light, Dirac-like carriers near the $H$ ($H^{\prime}$)- points of the Brillouin zone. The Hall resistivity in this interval also shows a non-analytic, $B\ln |B|$ behavior, and is dominated by holes.


title: Limits of local algorithms over sparse random graphs

abstract: Local algorithms on graphs are algorithms that run in parallel on the nodes of a graph to compute some global structural feature of the graph. Such algorithms use only local information available at nodes to determine local aspects of the global structure, while also potentially using some randomness. Recent research has shown that such algorithms show significant promise in computing structures like large independent sets in graphs locally. Indeed the promise led to a conjecture by Hatami, \Lovasz and Szegedy \cite{HatamiLovaszSzegedy} that local algorithms may be able to compute maximum independent sets in (sparse) random $d$-regular graphs. In this paper we refute this conjecture and show that every independent set produced by local algorithms is multiplicative factor $1/2+1/(2\sqrt{2})$ smaller than the largest, asymptotically as $d\rightarrow\infty$. Our result is based on an important clustering phenomena predicted first in the literature on spin glasses, and recently proved rigorously for a variety of constraint satisfaction problems on random graphs. Such properties suggest that the geometry of the solution space can be quite intricate. The specific clustering property, that we prove and apply in this paper shows that typically every two large independent sets in a random graph either have a significant intersection, or have a nearly empty intersection. As a result, large independent sets are clustered according to the proximity to each other. While the clustering property was postulated earlier as an obstruction for the success of local algorithms, such as for example, the Belief Propagation algorithm, our result is the first one where the clustering property is used to formally prove limits on local algorithms.


title: Distinguishing eigenforms modulo a prime ideal

abstract: Consider the Fourier expansions of two elements of a given space of modular forms. How many leading coefficients must agree in order to guarantee that the two expansions are the same? Sturm gave an upper bound for modular forms of a given weight and level. This was adapted by Ram Murty, Kohnen and Ghitza to the case of two eigenforms of the same level but having potentially different weights. We consider their expansions modulo a prime ideal, presenting a new bound. In the process of analysing this bound, we generalise a result of Bach and Sorenson, who provide a practical upper bound for the least prime in an arithmetic progression.


title: Fermi surface topology and de Hass-van Alphen orbits in PuIn$_{\rm 3}$ and PuSn$_{\rm 3}$ compounds

abstract: Since the recent discovery of plutonium-based superconductors such as PuCoIn$_{\rm 5}$, systematic studies of the electronic properties for plutonium compounds are providing insight into the itinerancy-localization crossover of Pu 5$f$ electrons. We are particularly interested in understanding the Fermi surface properties of PuIn$_{3}$ compound, which serves as the building block for the PuCoIn$_{\rm 5}$ superconductor. Motivated by the first observation of quantum oscillation and renewed interest in the de-Hass van-Alphen (dHvA) measurements on PuIn$_{\rm 3}$, we study the Fermi surface (FS) topology and the band dispersion in both paramagnetic and antiferromagnetic state of PuIn$_{\rm 3}$, based on density functional theory with generalized gradient approximation. We compare the results with its isostructural paramagnetic compound PuSn$_{\rm 3}$. We show the detailed Fermi surfaces of compounds PuIn$_{\rm 3}$ and PuSn$_{\rm 3}$ and conclude that the FS topology of an antiferromagnetic PuIn$_3$ agrees better with dHvA measurements. In addition, we show the magnetization of the antiferromagnetic order can alter the angle dependence and values of the effective mass for the dHvA orbits. We point out the accurate determination of the magnetic order orientation with respect to the crystal orientation is crucial to advance the understanding of the electronic structures of the PuIn$_{\rm 3}$ compound.


title: Detecting the structure of haplotypes, local ancestry and excessive local European ancestry in Mexicans

abstract: We present a two-layer hidden Markov model to detect structure of haplotypes for unrelated individuals. This allows modeling two scales of linkage disequilibrium (one within a group of haplotypes and one between groups), thereby taking advantage of rich haplotype information to infer local ancestry for admixed individuals. Our method outperforms competing state-of-art methods, particularly for regions of small ancestral track lengths. Applying our method to Mexican samples in HapMap3, we found five coding regions, ranging from $0.3 -1.3$ megabase (Mb) in lengths, that exhibit excessive European ancestry (average dosage > 1.6). A particular interesting region of 1.1Mb (with average dosage 1.95) locates on Chromosome 2p23 that harbors two genes, PXDN and MYT1L, both of which are associated with autism and schizophrenia. In light of the low prevalence of autism in Hispanics, this region warrants special attention. We confirmed our findings using Mexican samples from the 1000 genomes project. A software package implementing methods described in the paper is freely available at \url{http://bcm.edu/cnrc/mcmcmc}


title: Locality Optimization for Data Parallel Programs

abstract: Productivity languages such as NumPy and Matlab make it much easier to implement data-intensive numerical algorithms. However, these languages can be intolerably slow for programs that don't map well to their built-in primitives. In this paper, we discuss locality optimizations for our system Parakeet, a just-in-time compiler and runtime system for an array-oriented subset of Python. Parakeet dynamically compiles whole user functions to high performance multi-threaded native code. Parakeet makes extensive use of the classic data parallel operators Map, Reduce, and Scan. We introduce a new set of data parallel operators,TiledMap, TiledReduce, and TiledScan, that break up their computations into local pieces of bounded size so as better to make use of small fast memories. We introduce a novel tiling transformation to generate tiled operators automatically. Applying this transformation once tiles the program for cache, and applying it again enables tiling for registers. The sizes for cache tiles are left unspecified until runtime, when an autotuning search is performed. Finally, we evaluate our optimizations on benchmarks and show significant speedups on programs that exhibit data locality.


title: A Simulation and Modeling of Access Points with Definition Language

abstract: This submission has been withdrawn by arXiv administrators because it contains fictitious content and was submitted under a pseudonym, which is against arXiv policy.


title: Data enriched linear regression

abstract: We present a linear regression method for predictions on a small data set making use of a second possibly biased data set that may be much larger. Our method fits linear regressions to the two data sets while penalizing the difference between predictions made by those two models. The resulting algorithm is a shrinkage method similar to those used in small area estimation. We find a Stein-type finding for Gaussian responses: when the model has 5 or more coefficients and 10 or more error degrees of freedom, it becomes inadmissible to use only the small data set, no matter how large the bias is. We also present both plug-in and AICc-based methods to tune our penalty parameter. Most of our results use an $L_2$ penalty, but we obtain formulas for $L_1$ penalized estimates when the model is specialized to the location setting. Ordinary Stein shrinkage provides an inadmissibility result for only 3 or more coefficients, but we find that our shrinkage method typically produces much lower squared errors in as few as 5 or 10 dimensions when the bias is small and essentially equivalent squared errors when the bias is large.


title: Towards a Workload for Evolutionary Analytics

abstract: Emerging data analysis involves the ingestion and exploration of new data sets, application of complex functions, and frequent query revisions based on observing prior query answers. We call this new type of analysis evolutionary analytics and identify its properties. This type of analysis is not well represented by current benchmark workloads. In this paper, we present a workload and identify several metrics to test system support for evolutionary analytics. Along with our metrics, we present methodologies for running the workload that capture this analytical scenario.


title: Reconstruction of Signals from Magnitudes of Redundant Representations: The Complex Case

abstract: This paper is concerned with the question of reconstructing a vector in a finite-dimensional complex Hilbert space when only the magnitudes of the coefficients of the vector under a redundant linear map are known. We present new invertibility results as well an iterative algorithm that finds the least-square solution and is robust in the presence of noise. We analyze its numerical performance by comparing it to the Cramer-Rao lower bound.


title: Dark Matter and Pulsar Origins of the Rising Cosmic Ray Positron Fraction in Light of New Data From AMS

abstract: The rise of the cosmic ray positron fraction with energy, as first observed with high confidence by PAMELA, implies that a large flux of high energy positrons has been recently (or is being currently) injected into the local volume of the Milky Way. With the new and much more precise measurement of the positron fraction recently provided by AMS, we revisit the question of the origin of these high energy positrons. We find that while some dark matter models (annihilating directly to electrons or muons) no longer appear to be capable of accommodating these data, other models in which ~1-3 TeV dark matter particles annihilate to unstable intermediate states could still be responsible for the observed signal. Nearby pulsars also remain capable of explaining the observed positron fraction. Future measurements of the positron fraction by AMS (using a larger data set), combined with their anticipated measurements of various cosmic ray secondary-to-primary ratios, may enable us to further discriminate between these remaining scenarios.


title: Quadratic Hyperbolicity Preservers & Multiplier Sequences

abstract: It is known (see \cite[Br\"and\'en, Lemma 2.7]{B10}) that a necessary condition for $T:=\sum Q_k(x) D^k$ to be hyperbolicity preserving is that $Q_k(x)$ and $Q_{k-1}(x)$ have interlacing zeros. We characterize all quadratic linear operators, as a consequence we find several classes of $P_n$-multiplier sequence.


title: Information Measures of Complexity, Emergence, Self-organization, Homeostasis, and Autopoiesis

abstract: This chapter reviews measures of emergence, self-organization, complexity, homeostasis, and autopoiesis based on information theory. These measures are derived from proposed axioms and tested in two case studies: random Boolean networks and an Arctic lake ecosystem. Emergence is defined as the information a system or process produces. Self-organization is defined as the opposite of emergence, while complexity is defined as the balance between emergence and self-organization. Homeostasis reflects the stability of a system. Autopoiesis is defined as the ratio between the complexity of a system and the complexity of its environment. The proposed measures can be applied at different scales, which can be studied with multi-scale profiles.


title: Bayesian Lensing Shear Measurement

abstract: We derive an estimator of weak gravitational lensing shear from background galaxy images that avoids noise-induced biases through a rigorous Bayesian treatment of the measurement. The derived shear estimator disposes with the assignment of ellipticities to individual galaxies that is typical of previous approaches to galaxy lensing. Shear estimates from the mean of the Bayesian posterior are unbiased in the limit of large number of background galaxies, regardless of the noise level on individual galaxies. The Bayesian formalism requires a prior describing the (noiseless) distribution of the target galaxy population over some parameter space; this prior can be constructed from low-noise images of a subsample of the target population, attainable from long integrations of a fraction of the survey field. We find two ways to combine this exact treatment of noise with rigorous treatment of the effects of the instrumental point-spread function and sampling. The Bayesian model fitting (BMF) method assigns a likelihood of the pixel data to galaxy models (e.g. Sersic ellipses). The Bayesian Fourier domain (BFD) method compresses the pixel data to a small set of weighted moments calculated after PSF correction in Fourier space. A numerical test using a simplified model of a biased galaxy measurement process demonstrates that the Bayesian formalism recovers applied shears to $<1$ part in $10^3$ accuracy and provides accurate uncertainty estimates. BFD is the first shear measurement algorithm that is model-free and requires no approximations or ad hoc assumptions in correcting for the effects of PSF, noise, or sampling on the galaxy images. These algorithms are good candidates for attaining the part-per-thousand shear inference required for hemisphere-scale weak gravitational lensing surveys. (abridged)


title: The causal meaning of Fisher's average effect

abstract: In order to formulate the Fundamental Theorem of Natural Selection, Fisher defined the average excess and average effect of a gene substitution. Finding these notions to be somewhat opaque, some authors have recommended reformulating Fisher's ideas in terms of covariance and regression, which are classical concepts of statistics. We argue that Fisher intended his two averages to express a distinction between correlation and causation. On this view the average effect is a specific weighted average of the actual phenotypic changes that result from physically changing the allelic states of homologous genes. We show that the statistical and causal conceptions of the average effect, perceived as inconsistent by Falconer, can be reconciled if certain relationships between the genotype frequencies and non-additive residuals are conserved. There are certain theory-internal considerations favoring Fisher's original formulation in terms of causality; for example, the frequency-weighted mean of the average effects equaling zero at each locus becomes a derivable consequence rather than an arbitrary constraint. More broadly, Fisher's distinction between correlation and causation is of critical importance to gene-trait mapping studies and the foundations of evolutionary biology.


title: Potential Networks, Contagious Communities, and Understanding Social Network Structure

abstract: In this paper we study how the network of agents adopting a particular technology relates to the structure of the underlying network over which the technology adoption spreads. We develop a model and show that the network of agents adopting a particular technology may have characteristics that differ significantly from the social network of agents over which the technology spreads. For example, the network induced by a cascade may have a heavy-tailed degree distribution even if the original network does not. This provides evidence that online social networks created by technology adoption over an underlying social network may look fundamentally different from social networks and indicates that using data from many online social networks may mislead us if we try to use it to directly infer the structure of social networks. Our results provide an alternate explanation for certain properties repeatedly observed in data sets, for example: heavy-tailed degree distribution, network densification, shrinking diameter, and network community profile. These properties could be caused by a sort of `sampling bias' rather than by attributes of the underlying social structure. By generating networks using cascades over traditional network models that do not themselves contain these properties, we can nevertheless reliably produce networks that contain all these properties. An opportunity for interesting future research is developing new methods that correctly infer underlying network structure from data about a network that is generated via a cascade spread over the underlying network.


title: Trapping hot quasi-particles in a high-power superconducting electronic cooler

abstract: The performance of hybrid superconducting electronic coolers is usually limited by the accumulation of hot quasi-particles in the superconducting leads. This issue is all the more stringent in large-scale and high-power devices, as required by applications. Introducing a metallic drain connected to the superconducting electrodes via a fine-tuned tunnel barrier, we efficiently remove quasi-particles and obtain electronic cooling from 300 mK down to 130 mK with a 400 pW cooling power. A simple thermal model accounts for the experimental observations.