title: On a Hypergraph Approach to Multistage Group Testing Problems

abstract: Group testing is a well known search problem that consists in detecting up to $s$ defective elements of the set $[t]=\{1,\ldots,t\}$ by carrying out tests on properly chosen subsets of $[t]$. In classical group testing the goal is to find all defective elements by using the minimal possible number of tests. In this paper we consider multistage group testing. We propose a general idea how to use a hypergraph approach to searching defects. For the case $s=2$, we design an explicit construction, which makes use of $2\log_2t(1+o(1))$ tests in the worst case and consists of $4$ stages. For the general case $s>2$, we provide an explicit construction, which uses $(2s-1)\log_2t(1+o(1))$ tests and consists of $2s-1$ rounds.


title: On Multistage Learning a Hidden Hypergraph

abstract: Learning a hidden hypergraph is a natural generalization of the classical group testing problem that consists in detecting unknown hypergraph $H_{un}=H(V,E)$ by carrying out edge-detecting tests. In the given paper we focus our attention only on a specific family $F(t,s,\ell)$ of localized hypergraphs for which the total number of vertices $|V| = t$, the number of edges $|E|\le s$, $s\ll t$, and the cardinality of any edge $|e|\le\ell$, $\ell\ll t$. Our goal is to identify all edges of $H_{un}\in F(t,s,\ell)$ by using the minimal number of tests. We develop an adaptive algorithm that matches the information theory bound, i.e., the total number of tests of the algorithm in the worst case is at most $s\ell\log_2 t(1+o(1))$. We also discuss a probabilistic generalization of the problem.


title: Abelian cosmic string in the extended Starobinsky model of gravity

abstract: We analyze numerically the behaviour of the solutions corresponding to an Abelian cosmic string taking into account an extension of the Starobinsky model, where the action of general relativity is replaced by $f(R) = R - 2\Lambda + \eta R^2 + \rho R^m$, with $m > 2$. As an interesting result, we find that the angular deficit which characterizes the cosmic string decreases as the parameters $\eta$ and $\rho$ increase. We also find that the cosmic horizon due to the presence of a cosmological constant is affected in such a way that it can grows or shrinks, depending on the vacuum expectation value of the scalar field and on the value of the cosmological constant


title: Nonlinear perturbed integral equations related to nonlocal boundary value problems

abstract: By topological arguments, we prove new results on the existence, non-existence, localization and multiplicity of nontrivial solutions of a class of perturbed nonlinear integral equations. These type of integral equations arise, for example, when dealing with boundary value problems where nonlocal terms occur in the differential equation and/or in the boundary conditions. Some examples are given to illustrate the theoretical results.


title: Line defects in Graphene: How doping affects the electronic and mechanical properties

abstract: Graphene and carbon nanotubes have extraordinary mechanical and electronic properties. Intrinsic line defects such as local non-hexagonal reconstructions or grain boundaries, however, significantly reduce the tensile strength, but feature exciting electronic properties. Here, we address the properties of line defects in graphene from first-principles on the level of full-potential density-functional theory, and assess doping as one strategy to strengthen such materials. We carefully disentangle the global and local effect of doping by comparing results from the virtual crystal approximation with those from local substitution of chemical species, in order to gain a detailed understanding of the breaking and stabilization mechanisms. We find that n-type doping or local substitution with nitrogen increases the ultimate tensile strength significantly. In particular, it can stabilize the defects beyond the ultimate tensile strength of the pristine material. We therefore propose this as a key strategy to strengthen graphenic materials. Furthermore, we find that doping and/or applying external stress lead to tunable and technologically interesting metal/semi-conductor transitions.


title: Threshold Disjunctive Codes

abstract: Let $1 \le s < t$, $N \ge 1$ be integers and a complex electronic circuit of size $t$ is said to be an $s$-active, $\; s \ll t$, and can work as a system block if not more than $s$ elements of the circuit are defective. Otherwise, the circuit is said to be an $s$-defective and should be substituted for the similar $s$-active circuit. Suppose that there exists a possibility to check the $s$-activity of the circuit using $N$ non-adaptive group tests identified by a conventional disjunctive $s$-code $X$ of size~$t$ and length~$N$. As usually, we say that any group test yields the positive response if the group contains at least one defective element. In this case, there is no any interest to look for the defective elements. We need to decide on the number of the defective elements in the circuit without knowing the code~$X$. In addition, the decision has the minimal possible complexity because it is based on the simple comparison of a fixed threshold $T$, $0 \le T \le N - 1$, with the number of positive responses $p$, $0 \le p \le N$, obtained after carrying out $N$ non-adaptive tests prescribed by the disjunctive $s$-code~$X$. For the introduced group testing problem, a new class of the well-known disjunctive $s$-codes called the threshold disjunctive $s$-codes is defined. The aim of our paper is to discuss both some constructions of suboptimal threshold disjunctive $s$-codes and the best random coding bounds on the rate of threshold disjunctive $s$-codes.


title: A new version of Brakke's local regularity theorem

abstract: Consider an integral Brakke flow $(\mu_t)$, $t\in [0,T]$, inside some ball in Euclidean space. If $\mu_{0}$ has small height, its measure does not deviate too much from that of a plane and if $\mu_{T}$ is non-empty, then Brakke's local regularity theorem yields that $(\mu_t)$ is actually smooth and graphical inside a smaller ball for times $t\in (C,T-C)$ for some constant $C$. Here we extend this result to times $t\in (C,T)$. The main idea is to prove that a Brakke flow that is initially locally graphical with small gradient will remain graphical for some time. Moreover we use the new local regularity theorem to generalise White's regularity theorem to Brakke flows.


title: Scalable Anomaly Ranking of Attributed Neighborhoods

abstract: Given a graph with node attributes, what neighborhoods are anomalous? To answer this question, one needs a quality score that utilizes both structure and attributes. Popular existing measures either quantify the structure only and ignore the attributes (e.g., conductance), or only consider the connectedness of the nodes inside the neighborhood and ignore the cross-edges at the boundary (e.g., density). In this work we propose normality, a new quality measure for attributed neighborhoods. Normality utilizes structure and attributes together to quantify both internal consistency and external separability. It exhibits two key advantages over other measures: (1) It allows many boundary-edges as long as they can be "exonerated"; i.e., either (i) are expected under a null model, and/or (ii) the boundary nodes do not exhibit the subset of attributes shared by the neighborhood members. Existing measures, in contrast, penalize boundary edges irrespectively. (2) Normality can be efficiently maximized to automatically infer the shared attribute subspace (and respective weights) that characterize a neighborhood. This efficient optimization allows us to process graphs with millions of attributes. We capitalize on our measure to present a novel approach for Anomaly Mining of Entity Neighborhoods (AMEN). Experiments on real-world attributed graphs illustrate the effectiveness of our measure at anomaly detection, outperforming popular approaches including conductance, density, OddBall, and SODA. In addition to anomaly detection, our qualitative analysis demonstrates the utility of normality as a powerful tool to contrast the correlation between structure and attributes across different graphs.


title: Nonlinearities and Noise-Signal Relations in Electronic Heat Transport via Molecules

abstract: We examine the electronic heat transport phenomena in nanoscale junctions composed of organic molecules coupled to two metallic reservoirs of different temperatures. The electronic heat flux and its dynamical noise properties are calculated within the scattering (Landauer) formalism with the transmission probability determined by using non-equilibrium Green's functions (NEGF technique). The method based on Taylor series expansion is used to determine nonlinear corrections to the electronic heat flux and its noise power spectral density with up to the second order terms with respect to the temperature difference. Our results show only limited applicability of ballistic Fourier's law and fluctuation-dissipation theorem to heat transport in molecular systems. We derived and tested numerically several signal-signal, noise-signal, and noise-noise relations applicable to nanoscale heat flow carried by electrons at strongly non-equilibrium conditions (similar formulas are expected for phonons and photons). Importantly, the special treatment proposed by us may be extended to higher order terms in order to address a variety of problems related to nonlinear thermal and electro-thermal effects which may occur at nanoscale.


title: Azimuthal decorrelation of jets widely separated in rapidity in pp collisions at sqrt(s) = 7 TeV

abstract: The decorrelation in the azimuthal angle between the most forward and the most backward jets (Mueller-Navelet jets) is measured in data collected in pp collisions with the CMS detector at the LHC at sqrt(s) = 7 TeV. The measurement is presented in the form of distributions of azimuthal-angle differences, Delta(phi), between the Mueller-Navelet jets, the average cosines of (pi - Delta(phi)), 2(pi - Delta(phi)), and 3(pi - Delta(phi)), and ratios of these cosines. The jets are required to have transverse momenta, pT, in excess of 35 GeV and rapidities, |y|, of less than 4.7. The results are presented as a function of the rapidity separation, Delta(y), between the Mueller-Navelet jets, reaching Delta(y) up to 9.4 for the first time. The results are compared to predictions of various Monte Carlo event generators and to analytical predictions based on the DGLAP and BFKL parton evolution schemes.


title: Hiding Scalar Higgs Portal Dark Matter

abstract: We investigate a Higgs portal dark matter model by extending the Standard model (SM) with a complex singlet, $S=(s+ia )/\sqrt{2}$, where $a $ is a dark matter candidate and $S$ gets no vacuum expectation value but $s$ mixies with the SM Higgs via a trilinear interaction. We point out an interesting scenario, where only quartic coupling contributes to the dark matter relic abundance and there is no tree level contribution to the direct detection. Numerical analysis shows that the direction detection cross section, which arises at the one-loop level, is about $2\sim 5$ orders below the current LUX bound. Constraints from Higgs measurements as well as collider signatures of the model at the LHC are studied.


title: Long-Wavelength, Free-Free Spectral Energy Distributions from Porous Stellar Winds

abstract: The influence of macroclumps for free-free spectral energy distributions (SEDs) of ionized winds is considered. The goal is to emphasize distinctions between microclumping and macroclumping effects. Microclumping can alter SED slopes and flux levels if the volume filling factor of the clumps varies with radius; however, the modifications are independent of the clump geometry. To what extent does macroclumping alter SED slopes and flux levels? In addressing the question, two specific types of macroclump geometries are explored: shell fragments ("pancake"-shaped) and spherical clumps. Analytic and semi-analytic results are derived in the limiting case that clumps never obscure one another. Numerical calculations based on a porosity formalism is used when clumps do overlap. Under the assumptions of a constant expansion, isothermal, and fixed ionization wind, the fragment model leads to results that are essentially identical to the microclumping result. Mass-loss rate determinations are not affected by porosity effects for shell fragments. By contrast, spherical clumps can lead to a reduction in long-wavelength fluxes, but the reductions are only significant for extreme volume filling factors.


title: The Optical Chirality Flux as a Useful Far-Field Probe of Chiral Near Fields

abstract: To optimize the interaction between chiral matter and highly twisted light, quantities that can help characterize chiral electromagnetic fields near nanostructures are needed. Here, by analogy with Poynting's theorem, we formulate the time-averaged conservation law of optical chirality in lossy dispersive media and identify the optical chirality flux as an ideal far-field observable for characterizing chiral optical near fields. Bounded by the conservation law, we show that it provides precise information, unavailable from circular dichroism spectroscopy, on the magnitude and handedness of highly twisted fields near nanostructures.


title: Multi-Constituent Simulation of Thrombus Deposition

abstract: In this paper, we present a spatio-temporal mathematical model for simulating the formation and growth of a thrombus. Blood is treated as a multi-constituent mixture comprised of a linear fluid phase and a thrombus (solid) phase. The transport and reactions of 10 chemical and biological species are incorporated using a system of coupled convection-reaction-diffusion (CRD) equations to represent three processes in thrombus formation: initiation, propagation and stabilization. Computational fluid dynamic (CFD) simulations using the libraries of OpenFOAM were performed for two illustrative benchmark problems: in vivo thrombus growth in an injured blood vessel and in vitro thrombus deposition in micro-channels (1.5mm x 1.6mm x 0.1mm) with small crevices (125{\mu}m x 75{\mu}m and 125{\mu}m x 137{\mu}m). For both problems, the simulated thrombus deposition agreed very well with experimental observations, both spatially and temporally. Based on the success with these two benchmark problems, which have very different flow conditions and biological environments, we believe that the current model will provide useful insight into the genesis of thrombosis in blood-wetted devices, and provide a tool for the design of less thrombogenic devices.


title: Second order analysis of geometric functionals of Boolean models

abstract: This paper presents asymptotic covariance formulae and central limit theorems for geometric functionals, including volume, surface area, and all Minkowski functionals and translation invariant Minkowski tensors as prominent examples, of stationary Boolean models. Special focus is put on the anisotropic case. In the (anisotropic) example of aligned rectangles, we provide explicit analytic formulae and compare them with simulation results. We discuss which information about the grain distribution second moments add to the mean values.


title: Relief R-CNN : Utilizing Convolutional Features for Fast Object Detection

abstract: R-CNN style methods are sorts of the state-of-the-art object detection methods, which consist of region proposal generation and deep CNN classification. However, the proposal generation phase in this paradigm is usually time consuming, which would slow down the whole detection time in testing. This paper suggests that the value discrepancies among features in deep convolutional feature maps contain plenty of useful spatial information, and proposes a simple approach to extract the information for fast region proposal generation in testing. The proposed method, namely Relief R-CNN (R2-CNN), adopts a novel region proposal generator in a trained R-CNN style model. The new generator directly generates proposals from convolutional features by some simple rules, thus resulting in a much faster proposal generation speed and a lower demand of computation resources. Empirical studies show that R2-CNN could achieve the fastest detection speed with comparable accuracy among all the compared algorithms in testing.


title: The Generalized Quasilinear Approximation: Application to Zonal Jets

abstract: Quasilinear theory is often utilized to approximate the dynamics of fluids exhibiting significant interactions between mean flows and eddies. In this paper we present a generalization of quasilinear theory to include dynamic mode interactions on the large scales. This generalized quasilinear (GQL) approximation is achieved by separating the state variables into large and small zonal scales via a spectral filter rather than by a decomposition into a formal mean and fluctuations. Nonlinear interactions involving only small zonal scales are then removed. The approximation is conservative and allows for scattering of energy between small-scale modes via the large scale (through non-local spectral interactions). We evaluate GQL for the paradigmatic problems of the driving of large-scale jets on a spherical surface and on the beta-plane and show that it is accurate even for a small number of large-scale modes. As this approximation is formally linear in the small zonal scales it allows for the closure of the system and can be utilized in direct statistical simulation schemes that have proved an attractive alternative to direct numerical simulation for many geophysical and astrophysical problems.

categories: physics.flu-dyn astro-ph.EP astro-ph.SR physics.ao-ph


title: Generalized Rabi models: diagonalization in the spin subspace and differential operators of Dunkl type

abstract: A discrete parity $\mathbb{Z}_2$ symmetry of a two parameter extension of the quantum Rabi model which smoothly interpolates between the latter and the Jaynes-Cummings model, and of the two-photon and the two-mode quantum Rabi models enables their diagonalization in the spin subspace. A more general statement is that the respective sets of $2\times 2$ hermitian operators of the Fulton-Gouterman type and those diagonal in the spin subspace are unitary equivalent. The diagonalized representation makes it transparent that any question about integrability and solvability can be addressed only at the level of ordinary differential operators of Dunkl type. Braak's definition of integrability is shown (i) to contradict earlier numerical studies and (ii) to imply that any physically reasonable differential operator of Fulton-Gouterman type is integrable.


title: Optimal designs for comparing regression models with correlated observations

abstract: We consider the problem of efficient statistical inference for comparing two regression curves estimated from two samples of dependent measurements. Based on a representation of the best pair of linear unbiased estimators in continuous time models as a stochastic integral, an efficient pair of linear unbiased estimators with corresponding optimal designs for finite sample size is constructed. This pair minimises the width of the confidence band for the difference between the estimated curves. We thus extend results readily available in the literature to the case of correlated observations and provide an easily implementable and efficient solution. The advantages of using such pairs of estimators with corresponding optimal designs for the comparison of regression models are illustrated via numerical examples.


title: Perspectives and completely positive maps

abstract: We study the filtering of the perspective of a regular operator map of several variables through a completely positive linear map. By this method we are able to extend known operator inequalities of two variables to several variables; with applications in the theory of operator means of several variables. We also extend Lieb-Ruskai's convexity theorem from two to $ n+1 $ operator variables.


title: Breakdown of Elasticity in Amorphous Solids

abstract: What characterises a solid is its way to respond to external stresses. Ordered solids, such crystals, display an elastic regime followed by a plastic one, both well understood microscopically in terms of lattice distortion and dislocations. For amorphous solids the situation is instead less clear, and the microscopic understanding of the response to deformation and stress is a very active research topic. Several studies have revealed that even in the elastic regime the response is very jerky at low temperature, resembling very much the one of disordered magnetic materials. Here we show that in a very large class of amorphous solids this behaviour emerges by decreasing the temperature as a phase transition where standard elastic behaviour breaks down. At the transition all non-linear elastic modulii diverge and standard elasticity theory does not hold anymore. Below the transition the response to deformation becomes history and time-dependent.


title: Numerical simulations of oscillating tori in a pseudo-Newtonian potential

abstract: We have modeled hydrodynamical tori in the Klu{\'z}niak-Lee pseudo-Newtonian potential. The tori in equilibrium were perturbed with uniform sub-sonic velocity fields: vertical, radial and diagonal respectively, and allowed to evolve in time. We identify the eigenmodes corresponding to those of slender tori. The results of our simulations are relevant to the investigation of high-frequency quasi-periodic oscillations observed in stellar-mass black hole binaries.


title: On the warm absorber in AGN outflow

abstract: Warm absorber (WA) is an ionised gas present in the line of sight to the AGN central engine. The effect of the absorber is imprinted in the absorption lines observed in X-ray spectra of AGN. In this work, we model the WA in Seyfert 1 galaxy Mrk 509 using its recently published shape of broad band spectral energy distribution (SED) as a continuum illuminating the absorber. Using the photoionization code {\sc Titan}, recently we have shown that the absorption measure distribution (AMD) found for this object can be successfully modelled as a single slab of gas in total pressure (radiation+gas) equilibrium, contrary to the usual models of constant density multiple slabs. We discuss the transmitted spectrum that would be recorded by an observer after the radiation from the nucleus passes through the WA.


title: Optimal Bounds on Functions of Quantum States under Quantum Channels

abstract: Let $\rho_1, \rho_2$ be quantum states and $(\rho_1,\rho_2) \mapsto D(\rho_1, \rho_2)$ be a scalar function such as the trace norm, the fidelity, and the relative entropy, etc. We determine optimal bounds for $D(\rho_1, \Phi(\rho_2))$ for $\Phi \in \mathcal{S}$ for different class of functions $D(\cdot, \cdot)$, where $\mathcal{S}$ is the set of unitary quantum channels, the set of mixed unitary channels, the set of unital quantum channels, and the set of all quantum channels.


title: Topological nature of in-gap bound states in disordered large-gap monolayer transition metal dichalcogenides

abstract: We propose a physical model based on disordered (a hole punched inside a material) monolayer transition metal dichalcogenides (TMDs) to demonstrate a large-gap quantum valley Hall insulator. We find an emergence of bound states lying inside the bulk gap of the TMDs. They are strongly affected by spin-valley coupling, rest- and kinetic- mass terms and the hole size. In addition, in the whole range of the hole size, at least two in-gap bound states with opposite angular momentum, circulating around the edge of the hole, exist. Their topological insulator (TI) feature is analyzed by the Chern number, characterized by spacial distribution of their probabilities and confirmed by energy dispersion curves (Energy vs. angular momentum). It not only sheds light on overcoming low-temperature operating limitation of existing narrow-gap TIs, but also opens an opportunity to realize valley- and spin- qubits.


title: On rank one perturbations of Hamiltonian system with periodic coefficients

abstract: From a theory developed by C. Mehl, et al., a theory of the rank one perturbation of Hamiltonian systems with periodic coefficients is proposed. It is showed that the rank one perturbation of the fundamental solution of Hamiltonian system with periodic coefficients is solution of its rank one perturbation. Some results on the consequences of the strong stability of these types of systems on their rank one perturbation is proposed. Two numerical examples are given to illustrate this theory.


title: Derivation of Hodgkin-Huxley equations for a Na+ channel from a master equation for coupled activation and inactivation

abstract: The Na+ current in nerve and muscle membranes may be described in terms of the activation variable m(t) and the inactivation variable h(t), which are dependent on the transitions of S4 sensors of each of the Na+ channel domains DI to DIV. The time-dependence of the Na+ current and the rate equations satisfied by m(t) and h(t) may be derived from the solution to a master equation which describes the coupling between two or three activation sensors regulating the Na+ channel conductance and a two stage inactivation process. If the inactivation rate from the closed or open states increases as the S4 sensors activate, a more general form for the Hodgkin-Huxley expression for the open state probability may be derived where m(t) is dependent on both activation and inactivation processes. The voltage dependence of the rate functions for inactivation and recovery from inactivation are consistent with the empirically determined expressions, and exhibit saturation for both depolarized and hyperpolarized clamp potentials.


title: Resiliency and Robustness of Complex, Multi-Genre Networks

abstract: We explore the resiliency and robustness of systems while viewing them as complex, multi-genre networks. The term "complex, multi-genre networks" refers to networks that combine several distinct genres - networks of physical resources, communication networks, information networks, and social and cognitive networks. We show that this perspective is fruitful and adds to our understanding of fundamental challenges and tradeoffs in robustness and resiliency, as well as potential solutions to the challenges. Study of systems as multi-genre networks is relatively uncommon; instead, it is customary in research and engineering literature to focus on a view of a network comprised of homogeneous elements, (e.g., a network of communication devices, or a network of social beings). Yet, most if not all real-world networks are multi-genre - it is hard to find any real system of a significant complexity that does not include a combination of interconnected physical elements, communication devices and channels, data collections, and human users forming an integrated, inter-dependent whole. Most approaches to improving resiliency and robustness involve compromises, and the key challenge is to find a favorable compromise. Such compromises involve reducing or managing the complexity of the network: coupling, rigidity and dependency. We discuss several of these compromises, e.g., performance vs resiliency; resiliency to one type of disruption vs resiliency to another disruption type; and complexity vs resiliency.


title: Concept Generation in Language Evolution

abstract: This thesis investigates the generation of new concepts from combinations of existing concepts as a language evolves. We give a method for combining concepts, and will be investigating the utility of composite concepts in language evolution and thence the utility of concept generation.


title: Long Short-Term Memory-Networks for Machine Reading

abstract: In this paper we address the question of how to render sequence-level networks better at handling structured input. We propose a machine reading simulator which processes text incrementally from left to right and performs shallow reasoning with memory and attention. The reader extends the Long Short-Term Memory architecture with a memory network in place of a single memory cell. This enables adaptive memory usage during recurrence with neural attention, offering a way to weakly induce relations among tokens. The system is initially designed to process a single sequence but we also demonstrate how to integrate it with an encoder-decoder architecture. Experiments on language modeling, sentiment analysis, and natural language inference show that our model matches or outperforms the state of the art.


title: Cosmology in new gravitational scalar-tensor theories

abstract: We investigate the cosmological applications of new gravitational scalar-tensor theories, which are novel modifications of gravity possessing 2+2 propagating degrees of freedom, arising from a Lagrangian that includes the Ricci scalar and its first and second derivatives. Extracting the field equations we obtain an effective dark energy sector that consists of both extra scalar degrees of freedom, and we determine various observables. We analyze two specific models and we obtain a cosmological behavior in agreement with observations, i.e. transition from matter to dark energy era, with the onset of cosmic acceleration. Additionally, for a particular range of the model parameters, the equation-of-state parameter of the effective dark energy sector can exhibit the phantom-divide crossing. These features reveal the capabilities of these theories, since they arise solely from the novel, higher-derivative terms.


title: Position-dependent mass quantum Hamiltonians: General approach and duality

abstract: We analyze a general family of position-dependent mass quantum Hamiltonians which are not self-adjoint and include, as particular cases, some Hamiltonians obtained in phenomenological approaches to condensed matter physics. We build a general family of self-adjoint Hamiltonians which are quantum mechanically equivalent to the non self-adjoint proposed ones. Inspired in the probability density of the problem, we construct an ansatz for the solutions of the family of self-adjoint Hamiltonians. We use this ansatz to map the solutions of the time independent Schrodinger equations generated by the non self-adjoint Hamiltonians into the Hilbert space of the solutions of the respective dual self-adjoint Hamiltonians. This mapping depends on both the position-dependent mass and on a function of position satisfying a condition that assures the existence of a consistent continuity equation. We identify the non self-adjoint Hamiltonians here studied to a very general family of Hamiltonians proposed in a seminal article of Harrison [1] to describe varying band structures in different types of metals. Therefore, we have self-adjoint Hamiltonians that correspond to the non self-adjoint ones found in Harrison's article. We analyze three typical cases by choosing a physical position-dependent mass and a deformed harmonic oscillator potential . We completely solve the Schrodinger equations for the three cases; we also find and compare their respective energy levels.


title: Energy projection and modified Laughlin states

abstract: We develop a method to efficiently calculate trial wave functions for quantum Hall systems which involve projection onto the lowest Landau level. The method essentially replaces lowest Landau level projection by projection onto the $M$ lowest eigenstates of a suitably chosen hamiltonian acting within the lowest Landau level. The resulting "energy projection" is a controlled approximation to the exact lowest Landau level projection which improves with increasing $M$. It allows us to study projected trial wave functions for system sizes close to the maximal sizes that can be reached by exact diagonalization and can be straightforwardly applied in any geometry. As a first application and test case, we study a class of trial wave functions first proposed by Girvin and Jach, which are modifications of the Laughlin states involving a single real parameter. While these modified Laughlin states probably represent the same universality class exemplified by the Laughlin wave functions, we show by extensive numerical work for systems on the sphere and torus that they provide a significant improvement of the variational energy, overlap with the exact wave function and properties of the entanglement spectrum.


title: C^m Eigenfunctions of Perron-Frobenius Operators and a New Approach to Numerical Computation of Hausdorff Dimension

abstract: We develop a new approach to the computation of the Hausdorff dimension of the invariant set of an iterated function system or IFS. In the one dimensional case, our methods require only C^3 regularity of the maps in the IFS. The key idea, which has been known in varying degrees of generality for many years, is to associate to the IFS a parametrized family of positive, linear, Perron-Frobenius operators L_s. The operators L_s can typically be studied in many different Banach spaces. Here, unlike most of the literature, we study L_s in a Banach space of real-valued, C^k functions, k >= 2; and we note that L_s is not compact, but has a strictly positive eigenfunction v_s with positive eigenvalue lambda_s equal to the spectral radius of L_s. Under appropriate assumptions on the IFS, the Hausdorff dimension of the invariant set of the IFS is the value s=s_* for which lambda_s =1. This eigenvalue problem is then approximated by a collocation method using continuous piecewise linear functions (in one dimension) or bilinear functions (in two dimensions). Using the theory of positive linear operators and explicit a priori bounds on the derivatives of the strictly positive eigenfunction v_s, we give rigorous upper and lower bounds for the Hausdorff dimension s_*, and these bounds converge to s_* as the mesh size approaches zero.


title: A Label Semantics Approach to Linguistic Hedges

abstract: We introduce a model for the linguistic hedges `very' and `quite' within the label semantics framework, and combined with the prototype and conceptual spaces theories of concepts. The proposed model emerges naturally from the representational framework we use and as such, has a clear semantic grounding. We give generalisations of these hedge models and show that they can be composed with themselves and with other functions, going on to examine their behaviour in the limit of composition.


title: Adjustably Robust Optimal Power Flows with Demand Uncertainty via Path-based Flows

abstract: We study the optimal power flow problem with switching (or, equivalently, the line expansion problem) under demand uncertainty. Specifically, we consider the line-use variables at the first stage and the current- or power-flow at the second stage of two affinely adjustably robust formulations.


title: 50+ Metrics for Calendar Mining

abstract: In this report we propose 50+ metrics which can be measured by organizations in order to identify improvements in various areas such as meeting efficiency, capacity planning or leadership skills, just to new a few. The notion of calendar mining is introduced and support is provided for performing the measurement by a reference data model and queries for all metrics defined.


title: The z Distribution of Hydrogen Clouds and Masers with Kinematic Distances

abstract: Data on HII regions, molecular clouds, and methanol masers have been used to estimate the Sun's distance from the symmetry plane zo and the vertical disk scale height h. Kinematic distance estimates are available for all objects in these samples. The Local-arm (Orion-arm) objects are shown to affect noticeably the pattern of the z distribution. The deviations from the distribution symmetry are particularly pronounced for the sample of masers with measured trigonometric parallaxes, where the fraction of Local-arm masers is large. The situation with the sample of HII regions in the solar neighborhood is similar. We have concluded that it is better to exclude the Local arm from consideration. Based on the model of a self-gravitating isothermal disk, we have obtained the following estimates from objects located in the inner region of the Galaxy (R<= Ro): zo= -5.7+/-0.5 pc and h2=24.1+/-0.9 pc from the sample of 639 methanol masers, zo=-7.6+/-0.4 pc and h2=28.6+/-0.5 pc from 878 HII regions, zo=-10.1+/-0.5 pc and h2=28.2+/-0.6 pc from 538 giant molecular clouds.


title: On the probability that two elements of a finite semigroup have the same right matrix

abstract: In this paper we study the probability that two elements selected at random with replacement from a given finite semigroup act the same by right translation on the semigroup, that is, the chosen elements have the same right matrix.


title: On estimating causal controlled direct and mediator effects for count outcomes without assuming sequential ignorability

abstract: Causal mediation analysis is an important statistical method in social and medical studies, as it can provide insights about why an intervention works and inform the development of future interventions. Currently, most causal mediation methods focus on mediation effects defined on a mean scale. However, in health-risk studies, such as alcohol or risky sex, outcomes are typically count data and heavily skewed. Thus, mediation effects in these setting would be natural on a rate ratio scale, such as in Poisson and negative binomial regression methods. Existing methods also mainly rely on the assumption of no unmeasured confounding between mediator and outcome. To allow for potential confounders between the mediator and outcome, we define the direct and mediator effects on a new scale and propose a multiplicative structural mean model for mediation analysis with count outcomes. The estimator is compared with both Poisson and negative binomial regression methods assuming sequential ignorability using a simulation study and a real world example about an alcohol-related intervention study. Mediation analyses using the new methods confirm the study hypothesis that the intervention decreases drinking by decreasing individual's normative perceptions of alcohol use.


title: Octonion Algebra and Noise-Free Fully Homomorphic Encryption (FHE) Schemes

abstract: Brakerski showed that linearly decryptable fully homomorphic encryption (FHE) schemes cannot be secure in the chosen plaintext attack (CPA) model. In this paper, we show that linearly decryptable FHE schemes cannot be secure even in the ciphertext only security model. Then we consider the maximum security that a linearly decryptable FHE scheme could achieve. This paper designs fully homomorphic symmetric key encryption (FHE) schemes without bootstrapping (that is, noise-free FHE schemes). The proposed FHE schemes are based on quaternion/octonion algebra and Jordan algebra over finite rings Z_n and are secure in the weak ciphertext-only security model assuming the hardness of solving multivariate quadratic equation systems and solving univariate high degree polynomial equation systems in Z_n. It is up to our knowledge that this is the first noise-free FHE scheme that has ever been designed with a security proof (even in the weak ciphertext-only security model). It is argued that the weak ciphertext-only security model is sufficient for various applications such as privacy preserving computation in cloud. As an example, the proposed FHE schemes are used to construct obfuscated programs. This example could be further used to show that the scheme presented in this paper could be combined with existing FHE schemes with bootstrapping to obtain more efficient FHE schemes with bootstrapping in the fully CPA model. At the end of the paper, we point out the insecurity of several recently proposed noise-free FHE schemes.


title: Near-field integration of a SiN nanobeam and a SiO$_2$ microcavity for Heisenberg-limited displacement sensing

abstract: Placing a nanomechanical object in the evanescent near-field of a high-$Q$ optical microcavity gives access to strong gradient forces and quantum-noise-limited displacement readout, offering an attractive platform for precision sensing technology and basic quantum optics research. Robustly implementing this platform is challenging, however, as it requires separating optically smooth surfaces by $\lesssim\lambda/10$. Here we describe a fully-integrated evanescent opto-nanomechanical transducer based on a high-stress Si$_3$N$_4$ nanobeam monolithically suspended above a SiO$_2$ microdisk cavity. Employing a novel vertical integration technique based on planarized sacrificial layers, we achieve beam-disk gaps as little as 25 nm while maintaining mechanical $Q\times f>10^{12}$ Hz and intrinsic optical $Q\sim10^7$. The combined low loss, small gap, and parallel-plane geometry result in exceptionally efficient transduction, characterizing by radio-frequency flexural modes with vacuum optomechanical coupling rates of 100 kHz, single-photon cooperativities in excess of unity, and zero-point frequency (displacement) noise amplitudes of 10 kHz (fm)/$\surd$Hz. In conjunction with the high power handling capacity of SiO$_2$ and low extraneous substrate noise, the transducer operates particularly well as a sensor. Deploying it in a 4 K cryostat, we recently demonstrated a displacement imprecision 40 dB below that at the standard quantum limit (SQL) with an imprecision-back-action product $<5\cdot\hbar$. In this report we provide a comprehensive description of device design, fabrication, and characterization, with an emphasis on extending Heisenberg-limited readout to room temperature. Towards this end, we describe a room temperature experiment in which a displacement imprecision 30 dB below that at the SQL and an imprecision-back-action product $<75\cdot\hbar$ is achieved.

categories: physics.optics cond-mat.mes-hall physics.ins-det quant-ph


title: Quantum electrodynamics and the electron self-energy in a deformed space with a minimal length scale

abstract: The main motivation to study models in the presence of a minimal length is to obtain a quantum field theory free of the divergences. In this way, in this paper, we have constructed a new framework for quantum electrodynamics embedded in a minimal length scale background. New operators are introduced and the Green function method was used for the solution of the field equations, i.e., the Maxwell, Klein-Gordon and Dirac equations. We have analyzed specifically the scalar field and its one loop propagator. The mass of the scalar field regularized by the minimal length was obtained. The QED Lagrangian containing a minimal length was also constructed and the divergences were analyzed. The electron and photon propagators, and the electron self-energy at one loop as a function of the minimal length was also obtained.


title: Level of helium enhancement among M3's horizontal-branch stars

abstract: (Abridged) The color and luminosity distribution of horizontal branch (HB) stars in globular clusters are sensitive probes of their original helium abundances. However, recent analyses based on visual and near-ultraviolet (UV) CMDs have provided conflicting results. To clarify the situation, we address the optimum ranges of applicability (in terms of the $T_{eff}$ range covered by the HB stars) for visual and near-UV CMDs, as far as application of this HB $Y$ test goes. We consider both Stromgren and Hubble Space Telescope (HST) bandpasses. Our results indicate that the distributions of HB stars in the y vs (b-y) plane can be a reliable indicator of the He content in cool blue HB (BHB) stars, particularly when a differential comparison between blue and red HB stars is carried out in the range $T_{eff}<8300$K. Conversely, we demonstrate that CMDs using the F336W filter have a much less straightforward interpretation at the cool end of the BHB, as the distributions of HB stars in the F336W vs (F336W-F555W) plane are affected by a triple degeneracy effect. In other words, the position of an HB star is exactly the same, for a given chemical composition, for multiple combinations of the parameters $Y$, $M_{HB}$, and age along the HB evolutionary track. Other HST UV filters do not appear to be as severely affected by this degeneracy effect, to which visual bandpasses are also immune. We do find, on the other hand, that such near-UV CMDs can be extremely useful for the hottest stars along the cool BHB end. Based on a reanalysis of the distribution of HB stars in the y vs (b-y) plane, we find that the coolest among the BHB stars in M3 are likely enhanced in helium by $\Delta Y\approx 0.01$, compared with the red HB stars in the same cluster. Using near-UV HST photometry we find evidence of a progressive increase in $Y$ with increasing temperature, reaching $\Delta Y=0.02$ at $T_{eff}=10900$K.


title: BOLA: Near-Optimal Bitrate Adaptation for Online Videos

abstract: Modern video players employ complex algorithms to adapt the bitrate of the video that is shown to the user. Bitrate adaptation requires a tradeoff between reducing the probability that the video freezes (rebuffers) and enhancing the quality of the video. A bitrate that is too high leads to frequent rebuffering, while a bitrate that is too low leads to poor video quality. Video providers segment videos into short segments and encode each segment at multiple bitrates. The video player adaptively chooses the bitrate of each segment to download, possibly choosing different bitrates for successive segments. We formulate bitrate adaptation as a utility-maximization problem and devise an online control algorithm called BOLA that uses Lyapunov optimization to minimize rebuffering and maximize video quality. We prove that BOLA achieves a time-average utility that is within an additive term O(1/V) of the optimal value, for a control parameter V related to the video buffer size. Further, unlike prior work, BOLA does not require prediction of available network bandwidth. We empirically validate BOLA in a simulated network environment using a collection of network traces. We show that BOLA achieves near-optimal utility and in many cases significantly higher utility than current state-of-the-art algorithms. Our work has immediate impact on real-world video players and for the evolving DASH standard for video transmission. We also implemented an updated version of BOLA that is now part of the standard reference player dash.js and is used in production by several video providers such as Akamai, BBC, CBS, and Orange.


title: Empirical bayes formulation of the elastic net and mixed-norm models: application to the eeg inverse problem

abstract: The estimation of EEG generating sources constitutes an Inverse Problem (IP) in Neuroscience. This is an ill-posed problem, due to the non-uniqueness of the solution, and many kinds of prior information have been used to constrain it. A combination of smoothness (L2 norm-based) and sparseness (L1 norm-based) constraints is a flexible approach that have been pursued by important examples such as the Elastic Net (ENET) and mixed-norm (MXN) models. The former is used to find solutions with a small number of smooth non-zero patches, while the latter imposes sparseness and smoothness simultaneously along different dimensions of the spatio-temporal matrix solutions. Both models have been addressed within the penalized regression approach, where the regularization parameters are selected heuristically, leading usually to non-optimal solutions. The existing Bayesian formulation of ENET allows hyperparameter learning, but using computationally intensive Monte Carlo/Expectation Maximization methods. In this work we attempt to solve the EEG IP using a Bayesian framework for models based on mixtures of L1/L2 norms penalization functions (Laplace/Normal priors) such as ENET and MXN. We propose a Sparse Bayesian Learning algorithm based on combining the Empirical Bayes and the iterative coordinate descent procedures to estimate both the parameters and hyperparameters. Using simple but realistic simulations we found that our methods are able to recover complicated source setups more accurately and with a more robust variable selection than the ENET and LASSO solutions using classical algorithms. We also solve the EEG IP using data coming from a visual attention experiment, finding more interpretable neurophysiological patterns with our methods, as compared with other known methods such as LORETA, ENET and LASSO FUSION using the classical regularization approach.


title: A Robust UCB Scheme for Active Learning in Regression from Strategic Crowds

abstract: We study the problem of training an accurate linear regression model by procuring labels from multiple noisy crowd annotators, under a budget constraint. We propose a Bayesian model for linear regression in crowdsourcing and use variational inference for parameter estimation. To minimize the number of labels crowdsourced from the annotators, we adopt an active learning approach. In this specific context, we prove the equivalence of well-studied criteria of active learning like entropy minimization and expected error reduction. Interestingly, we observe that we can decouple the problems of identifying an optimal unlabeled instance and identifying an annotator to label it. We observe a useful connection between the multi-armed bandit framework and the annotator selection in active learning. Due to the nature of the distribution of the rewards on the arms, we use the Robust Upper Confidence Bound (UCB) scheme with truncated empirical mean estimator to solve the annotator selection problem. This yields provable guarantees on the regret. We further apply our model to the scenario where annotators are strategic and design suitable incentives to induce them to put in their best efforts.


title: Wgamma and Zgamma production at the LHC in NNLO QCD

abstract: We consider the production of Wgamma and Zgamma pairs at the LHC, and report on the fully differential computation of next-to-next-to-leading order (NNLO) corrections in QCD perturbation theory. The calculation includes leptonic vector-boson decays with the corresponding spin correlations, off-shell effects and final-state photon radiation. We present numerical results for pp collisions at 7 TeV, and compare them with available ATLAS data. In the case of Zgamma production, the impact of NNLO corrections is generally moderate, ranging from 8% to 17%, depending on the applied cuts. In the case of Wgamma production, the NNLO effects are more important, and range from 19% to 26%, thereby improving the agreement of the theoretical predictions with the data. As expected, a veto against jets significantly reduces the impact of QCD radiative corrections.


title: Device-independent two-party cryptography secure against sequential attacks

abstract: The goal of two-party cryptography is to enable two parties, Alice and Bob, to solve common tasks without the need for mutual trust. Examples of such tasks are private access to a database, and secure identification. Quantum communication enables security for all of these problems in the noisy-storage model by sending more signals than the adversary can store in a certain time frame. Here, we initiate the study of device-independent protocols for two-party cryptography in the noisy-storage model. Specifically, we present a relatively easy to implement protocol for a cryptographic building block known as weak string erasure and prove its security even if the devices used in the protocol are prepared by the dishonest party. Device-independent two-party cryptography is made challenging by the fact that Alice and Bob do not trust each other, which requires new techniques to establish security. We fully analyse the case of memoryless devices (for which sequential attacks are optimal) and the case of sequential attacks for arbitrary devices. The key ingredient of the proof, which might be of independent interest, is an explicit (and tight) relation between the violation of the Clauser-Horne-Shimony-Holt inequality observed by Alice and Bob and uncertainty generated by Alice against Bob who is forced to measure his system before finding out Alice's setting (guessing with postmeasurement information). In particular, we show that security is possible for arbitrarily small violation.


title: Homogenization of Fucik eigenvalues by optimal partition methods

abstract: Given a bounded domain $\Omega$ in $\mathbb{R}^N$, $N\geq 1$ we study the asymptotic behavior as $\varepsilon \to 0$ of the eigencurves of $$ -\Delta_p u_\varepsilon=\alpha_\varepsilon m(\tfrac{x}{\varepsilon})(u_\varepsilon^+ )^{p-1} - \beta_\varepsilon n(\tfrac{x}{\varepsilon})(u_\varepsilon^- )^{p-1} \quad \textrm{ in } \Omega $$ with Dirichlet boundary conditions, where $m$ and $n$ are bounded periodic weights. In this work we obtain accurate bounds of the convergence rates of these curves to some limit curves as $\varepsilon \to 0$.


title: On the Limited Sizes of Dusty Starbursting Regions at High Redshifts

abstract: Using the far-infrared data obtained by the Herschel Space Observatory, we study the relation between the infrared luminosity (L_IR) and the dust temperature (T) of dusty starbursting galaxies at high redshifts (high-z). We focus on the total infrared luminosity from the cold-dust component (L_IR^(cd)), whose emission can be described by a modified black body (MBB) of a single temperature (T_mbb). An object on the (L_IR^(cd), T_mbb) plane can be explained by the equivalent of the Stefan-Boltzmann law for a MBB with an effective radius of R_eff. We show that R_eff is a good measure of the combined size of the dusty starbursting regions (DSBRs) of the host galaxy. In at least one case where the individual DSBRs are well resolved through strong gravitational lensing, R_eff is consistent with the direct size measurement. We show that the observed L_IR-T relation is simply due to the limited R_eff (<~ 2 kpc). The small R_eff values also agree with the compact sizes of the DSBRs seen in the local universe. However, previous interferometric observations to resolve high-z dusty starbursting galaxies often quote much larger sizes. This inconsistency can be reconciled by the blending effect when considering that the current interferometry might still not be of sufficient resolution. From R_eff we infer the lower limits to the volume densities of the star formation rate ("minSFR3D") in the DSBRs, and find that the $L_{IR}$-$T$ relation outlines a boundary on the (L_IR^(cd), T) plane, below which is the "zone of avoidance" in terms of minSFR3D.


title: The Utility of Hedged Assertions in the Emergence of Shared Categorical Labels

abstract: We investigate the emergence of shared concepts in a community of language users using a multi-agent simulation. We extend results showing that negated assertions are of use in developing shared categories, to include assertions modified by linguistic hedges. Results show that using hedged assertions positively affects the emergence of shared categories in two distinct ways. Firstly, using contraction hedges like `very' gives better convergence over time. Secondly, using expansion hedges such as `quite' reduces concept overlap. However, both these improvements come at a cost of slower speed of development.


title: Privacy, Secrecy, and Storage with Multiple Noisy Measurements of Identifiers

abstract: The key-leakage-storage region is derived for a generalization of a classic two-terminal key agreement model. The additions to the model are that the encoder observes a hidden, or noisy, version of the identifier, and that the encoder and decoder can perform multiple measurements. To illustrate the behavior of the region, the theory is applied to binary identifiers and noise modeled via binary symmetric channels. In particular, the key-leakage-storage region is simplified by applying Mrs. Gerber's lemma twice in different directions to a Markov chain. The growth in the region as the number of measurements increases is quantified. The amount by which the privacy-leakage rate reduces for a hidden identifier as compared to a noise-free (visible) identifier at the encoder is also given. If the encoder incorrectly models the source as visible, it is shown that substantial secrecy leakage may occur and the reliability of the reconstructed key might decrease.

categories: cs.IT cs.CR cs.MM eess.SP math.IT math.PR


title: Frequency-dependent polarizability, plasmons, and screening in the 2D pseudospin-1 dice lattice

abstract: We calculate the dynamic polarizability under the random phase approximation for the dice lattice. This two-dimensional system gives rise to massless Dirac fermions with pseudospin-1 in the low-energy quantum excitation spectrum, providing a Dirac-cone plus flat-band dispersion. Due to the presence of the flat band, the polarizability shows key differences to that of graphene (the pseudospin-1/2 Dirac material). We find that the plasmon branch is pinched in to a single point, $\omega_p=q=\mu$, independent of the background dielectric constant. Finally, screening effects are discussed with regard to impurities.


title: Master Stability Islands for Amplitude Death in Networks of Delay-Coupled Oscillators

abstract: This paper presents a master stability function (MSF) approach for analyzing the stability of amplitude death (AD) in networks of delay-coupled oscillators. Unlike the familiar MSFs for instantaneously coupled networks, which typically have a single input encoding for the effects of the eigenvalues of the network Laplacian matrix, for delay-coupled networks we show that such MSFs generally require two additional inputs: the time delay and the coupling strength. To utilize the MSF for predicting the stability of AD of arbitrary networks for a chosen nonlinear system (node dynamics) and coupling function, we introduce the concept of master stability islands (MSIs), which are two-dimensional stability islands of the delay-coupling space together with a third dimension ("altitude") encoding for eigenvalues that result in stable AD. We compute the MSFs and show the corresponding MSIs for several common chaotic systems including the Rossler, the Lorenz, and Chen's system, and found that it is generally possible to achieve AD and that a nonzero time delay is necessary for the stabilization of the AD states.


title: Pixel Recurrent Neural Networks

abstract: Modeling the distribution of natural images is a landmark problem in unsupervised learning. This task requires an image model that is at once expressive, tractable and scalable. We present a deep neural network that sequentially predicts the pixels in an image along the two spatial dimensions. Our method models the discrete probability of the raw pixel values and encodes the complete set of dependencies in the image. Architectural novelties include fast two-dimensional recurrent layers and an effective use of residual connections in deep recurrent networks. We achieve log-likelihood scores on natural images that are considerably better than the previous state of the art. Our main results also provide benchmarks on the diverse ImageNet dataset. Samples generated from the model appear crisp, varied and globally coherent.


title: Symmetries of Analytic Curves

abstract: Analytic curves are classified w.r.t. their symmetry under a regular and separately analytic Lie group action on an analytic manifold. We show that an analytic curve is either exponential or splits into countably many analytic immersive curves, each of them discretely generated by the symmetry group (i.e., each such curve naturally decomposes into countably many symmetry free subcurves that are mutually and uniquely related by the Lie group action). We additionally extend the classification result to the analytic 1-submanifold case. Specifically, we show that an analytic 1-submanifold is either free or (exponential, i.e.) analytically diffeomorphic (via the exponential map) to the unit circle or an interval. The corresponding decomposition results in the free case are outlined in this paper, but proven in a separate one.


title: LHC probes of TeV-scale scalars in $\mathrm{SO}(10)$ grand unification

abstract: We investigate the possibility of TeV-scale scalars as low energy remnants arising in the non-supersymmetric $\mathrm{SO}(10)$ grand unification framework where the field content is minimal. We consider a scenario where the $\mathrm{SO}(10)$ gauge symmetry is broken into the gauge symmetry of the Standard Model (SM) through multiple stages of symmetry breaking, and a colored and hypercharged scalar $\chi$ picks a TeV-scale mass in the process. The last stage of the symmetry breaking occurs at the TeV scale where the left-right symmetry, i.e. $\mathrm{SU}(2)_L\otimes \mathrm{SU}(2)_R\otimes \mathrm{U}(1)_{B-L}\otimes \mathrm{SU}(3)_C$, is broken into that of the SM by a singlet scalar field $\mathcal{S}$ of mass $M_{\mathcal{S}}\sim 1$ TeV, which is a component of an $\mathrm{SU}(2)_R$-triplet scalar field, acquiring a TeV-scale vacuum expectation value. For the LHC phenomenology, we consider a scenario where $\mathcal{S}$ is produced via gluon-gluon fusion through loop interactions with $\chi$ and also decays to a pair of SM gauge bosons through $\chi$ in the loop. We find that the parameter space is heavily constrained from the latest LHC data. We use a multivariate analysis to estimate the LHC discovery reach of $\mathcal{S}$ into the diphoton channel.


title: Energy-Aware Optimization and Mechanism Design for Cellular Device-to-Device Local Area Networks

abstract: In a device-to-device (D2D) local area network (LAN), mobile users (MUs) must cooperate to download common real-time content from a wireless cellular network. However, sustaining such D2D LANs over cellular networks requires the introduction of mechanisms that will incentivize the MUs to cooperate. In this paper, the problem of energy-aware D2D LAN formation over cellular networks is studied. The problem is formulated using a game-theoretic framework in which each MU seeks to minimize its energy consumption while actively participating in the D2D LAN. To account for the selfish behavior of the MUs, a punishment and incentive protocol is proposed in order to ensure cooperation among MUs. Within this protocol, an estimation algorithm is proposed to simulate the process of D2D LAN formation and, then, adjust the mechanism parameters to maintain cooperation. Simulation results show that the proposed framework can improve energy efficiency up to 36% relative to the traditional multicast scenario.


title: Emerging Dimension Weights in a Conceptual Spaces Model of Concept Combination

abstract: We investigate the generation of new concepts from combinations of properties as an artificial language develops. To do so, we have developed a new framework for conjunctive concept combination. This framework gives a semantic grounding to the weighted sum approach to concept combination seen in the literature. We implement the framework in a multi-agent simulation of language evolution and show that shared combination weights emerge. The expected value and the variance of these weights across agents may be predicted from the distribution of elements in the conceptual space, as determined by the underlying environment, together with the rate at which agents adopt others' concepts. When this rate is smaller, the agents are able to converge to weights with lower variance. However, the time taken to converge to a steady state distribution of weights is longer.


title: FreeSASA: An open source C library for solvent accessible surface area calculations

abstract: Calculating solvent accessible surface areas (SASA) is a run-of-the-mill calculation in structural biology. Although there are many programs available for this calculation, there are no free-standing, open-source tools designed for easy tool-chain integration. FreeSASA is an open source C library for SASA calculations that provides both command-line and Python interfaces in addition to its C API. The library implements both Lee and Richards' and Shrake and Rupley's approximations, and is highly configurable to allow the user to control molecular parameters, accuracy and output granularity. It only depends on standard C libraries and should therefore be easy to compile and install on any platform. The source code is freely available from http://freesasa.github.io/. The library is well-documented, stable and efficient. The command-line interface can easily replace closed source legacy programs, with comparable or better accuracy and speed, and with some added functionality.


title: Counting roots of truncated hypergeometric series over finite fields

abstract: We consider natural polynomial truncations of hypergeometric power series defined over finite fields. For these truncations, we establish asymptotic upper bounds of order $O(p^{11/12})$ on the number of roots in the prime field $\mathbb{F}_p$. We discuss the correspondence to families of elliptic curves and K3 surfaces of certain such hypergeometric polynomials, for which sharp bounds are obtained in some cases. We include some computations to illustrate and supplement our results.


title: On the observation of nonclassical excitations in Bose-Einstein condensates

abstract: In the recent experimental and theoretical literature well-established nonclassicality criteria from the field of quantum optics have been directly applied to the case of excitations in matter-waves. Among these are violations of Cauchy-Schwarz inequalities, Glauber-Sudarshan P-nonclassicality, sub-Poissonian number-difference squeezing (also known as the two-mode variance) and the criterion of nonseparability. We review the strong connection of these criteria and their meaning in quantum optics, and point out differences in the interpretation between light and matter waves. We then calculate observables for a homogenous Bose-Einstein condensate undergoing an arbitrary modulation in the interaction parameter at finite initial temperature, within both the quantum theory as well as a classical reference. We conclude that to date in experiments relevant for analogue gravity, nonclassical effects have not conclusively been observed and conjecture that additional, noncommuting, observables have to be measured to this end.


title: Towards Fully Converged GW Calculations for Large Systems

abstract: Although the GW approximation is recognized as one of the most accurate theories for predicting materials excited states properties, scaling up conventional GW calculations for large systems remains a major challenge. We present a powerful and simple-to-implement method that can drastically accelerate fully converged GW calculations for large systems. We demonstrate the performance of this new method by calculating the quasiparticle band gap of MgO supercells. A speed-up factor of nearly two orders of magnitude is achieved for a system contaning 256 atoms (1024 velence electrons) with a negligibly small numerical error of $\pm 0.03$ eV.


title: The Spectrum in the Sachdev-Ye-Kitaev Model

abstract: The SYK model consists of $N\gg 1$ fermions in $0+1$ dimensions with a random, all-to-all quartic interaction. Recently, Kitaev has found that the SYK model is maximally chaotic and has proposed it as a model of holography. We solve the Schwinger-Dyson equation and compute the spectrum of two-particle states in SYK, finding both a continuous and discrete tower. The four-point function is expressed as a sum over the spectrum. The sum over the discrete tower is evaluated.


title: Dust Emissivity in the Star-Forming Filament OMC 2/3

abstract: We present new measurements of the dust emissivity index, beta, for the high-mass, star-forming OMC 2/3 filament. We combine 160-500 um data from Herschel with long-wavelength observations at 2 mm and fit the spectral energy distributions across a ~ 2 pc long, continuous section of OMC 2/3 at 15000 AU (0.08 pc) resolution. With these data, we measure beta and reconstruct simultaneously the filtered-out large-scale emission at 2 mm. We implement both variable and fixed values of beta, finding that beta = 1.7 - 1.8 provides the best fit across most of OMC 2/3. These beta values are consistent with a similar analysis carried out with filtered Herschel data. Thus, we show that beta values derived from spatial filtered emission maps agree well with those values from unfiltered data at the same resolution. Our results contradict the very low beta values (~ 0.9) previously measured in OMC 2/3 between 1.2 mm and 3.3 mm data, which we attribute to elevated fluxes in the 3.3 mm observations. Therefore, we find no evidence or rapid, extensive dust grain growth in OMC 2/3. Future studies with Herschel data and complementary ground-based long-wavelength data can apply our technique to obtain robust determinations of beta in nearby cold molecular clouds.


title: Dust grains from the heart of supernovae

abstract: Dust grains are classically thought to form in the winds of AGB stars. However, nowadays there is increasing evidence for dust formation in SNe. In order to establish the relative importance of these two classes of stellar sources of dust it is important to know what is the fraction of freshly formed dust in SN ejecta that is able to survive the passage of the reverse shock and be injected in the interstellar medium. With this aim, we have developed a new code, GRASH_Rev, that allows to follow the dynamics of dust grains in the shocked SN ejecta and to compute the time evolution of the mass, composition and size distribution of the grains. We consider four well studied SNe in the Milky Way and LMC: SN 1987a, Cas A, the Crab Nebula, and N49. For all the simulated models, we find good agreement with observations. Our study suggests that SN 1987A is too young for the reverse shock to have affected the dust mass. Conversely, in the other three SNe, the reverse shock has already destroyed between 10 and 40% of the initial dust mass. However, the largest dust mass destruction is predicted to occur between 10^3 and 10^5 yr after the explosions. Since the oldest SN in the sample has an estimated age of 4800 yr, current observations can only provide an upper limit to the mass of SN dust that will enrich the interstellar medium, the so-called effective dust yields. We find that only between 1 and 8% of the currently observed mass will survive. This is in good agreement with the values adopted in chemical evolution models which consider the effect of the SN reverse shock. We discuss the astrophysical implications of our results for dust enrichment in local galaxies and at high redshift.


title: On Axially Symmetric Space-Times Admitting Homothetic Vector Fields in Lyra's Geometry

abstract: This paper investigates axially symmetric space-times that admit a homothetic vector field based on Lyra's geometry. The cases when the displacement vector is a function of $t$ and when it is constant are studied. In the context of this geometry, we find and classify the solutions of the Einstein's field equations (EFE) for the space-time under consideration, which display a homothetic symmetry.


title: Tachyonic models of dark matter

abstract: We consider a spherically symmetric stationary problem in General Relativity, including a black hole, inflow of normal and tachyonic matter and outflow of tachyonic matter. Computations in a weak field limit show that the resulting concentration of matter around the black hole leads to gravitational effects equivalent to those associated with dark matter halo. In particular, the model reproduces asymptotically constant galactic rotation curves, if the tachyonic flows of the central supermassive black hole in the galaxy are considered as a main contribution.


title: Anisotropic String Cosmological Models in Heckmann-Suchuking Space-Time

abstract: In the present work we have searched the existence of the late time acceleration of the universe with string fluid as source of matter in anisotropic Heckmann-Suchking space-time by using 287 high red shift $(0.3 \leq z\leq 1.4)$ SN Ia data of observed absolute magnitude along with their possible error from Union 2.1 compilation. It is found that the best fit values for $(\Omega_{m})_{0}$, $(\Omega_{\Lambda})_{0}$, $(\Omega_{\sigma})_{0}$ and $(q)_{0}$ are 0.2820, 0.7177, 0.0002 $\&$ -0.5793 respectively. Several physical aspects and geometrical properties of the model are discussed in detail.


title: An Asymptotic Expansion for Perturbations in the Displacement Field due to the Presence of Thin Interfaces

abstract: We derive an asymptotic expansion for two-dimensional displacement field associated to thin elastic inhomogeneities having no uniform thickness. Our derivation is rigorous and based on layer potential techniques. We extend these techniques to determine a relationship between traction-displacement measurements and the shape of the object.


title: A family of semifields in characteristic 2

abstract: We construct and describe the basic properties of a family of semifields in characteristic $2.$ The construction relies on the properties of projective polynomials over finite fields. We start by associating non-associative products to each such polynomial. The resulting presemifields form the degenerate case of our family. They are isotopic to the Knuth semifields which are quadratic over left and right nucleus. The non-degenerate members of our family display a very different behaviour. Their left and right nucleus agrees with the center, the middle nucleus is quadratic over the center. None of those semifields is isotopic or Knuth equivalent to a commutative semifield. As a by-product we obtain the complete taxonomy of the characteristic $2$ semifields which are quadratic over the middle nucleus, bi-quadratic over the left and right nucleus and not isotopic to twisted fields. This includes {determining} when two such semifields are isotopic and the order of the autotopism group.


title: Composition operators on Orlicz-Sobolev spaces

abstract: The kernel of composition operator $C_T$ on Orlicz-Sobolev space is obtained. Using the kernel, a necessary and a sufficient condition for injectivity of composition operator $C_T$ has been established. Composition operators on Orlicz-Sobolev space with finite ascent as well as infinite ascent have been characterized.


title: Operator means of probability measures and generalized Karcher equations

abstract: In this article we consider means of positive bounded linear operators on a Hilbert space. We present a complete theory that provides a framework which extends the theory of the Karcher mean, its approximating matrix power means, and a large part of Kubo-Ando theory to arbitrary many variables, in fact, to the case of probability measures with bounded support on the cone of positive definite operators. This framework characterizes each operator mean extrinsically as unique solutions of generalized Karcher equations which are obtained by exchanging the matrix logarithm function in the Karcher equation to arbitrary operator monotone functions over the positive real half-line. If the underlying Hilbert space is finite dimensional, then these generalized Karcher equations are Riemannian gradients of convex combinations of strictly geodesically convex log-determinant divergence functions, hence these new means are the global minimizers of them, in analogue to the case of the Karcher mean as pointed out. Our framework is based on fundamental contraction results with respect to the Thompson metric, which provides us nonlinear contraction semigroups in the cone of positive definite operators that form a decreasing net approximating these operator means in the strong topology from above.


title: Enhancement of perovskite solar cells by plasmonic nanoparticles

abstract: Synthetic perovskites with photovoltaic properties open a new era in solar photovoltaics. Due to high optical absorption perovskite-based thin-film solar cells are usually considered as fully absorbing solar radiation on condition of ideal blooming. However, is it really so? The analysis of the literature data has shown that the absorbance of all photovoltaic pervoskites has the spectral hole at infrared frequencies where the solar radiation spectrum has a small local peak. This absorption dip results in the decrease of the optical efficiency of thin-film pervoskite solar cells by nearly 3% and close the ways of utilise them at this range for any other applications. In our work we show that to cure this shortage is possible complementing the basic structure by an inexpensive plasmonic array.


title: Non-Commutative Stochastic Independence and Cumulants

abstract: In a central lemma we characterize "generating functions" of certain functors on the category of algebraic non-commutative probability spaces. Special families of such generating functions correspond to "unital, associative universal products" on this category, which again define a notion of non-commutative stochastic independence. Using the central lemma, we can prove the existence of cumulants and of "cumulant Lie-algebras" for a wide class of independences. These include the five independences (tensor, free, Boolean, monotone, anti-monotone) appearing in N. Murakis classification, c-free independence of M. Bozejko and R. Speicher, the indented product of T. Hasebe and the bi-free independence of D. Voiculescu. We show that the non-commutative independence can be reconstructed from its cumulants and cumulant Lie algebras.


title: Multivariate information measures: a unification using M\"obius operators on subset lattices

abstract: Information related measures are useful tools for multi variable data analysis, as measures of dependence among variables, and as descriptions of order in biological and physical systems. Information related measures, like marginal entropies, mutual / interaction / multi-information, have been used in a number of fields including descriptions of systems complexity and biological data analysis. The mathematical relationships among these measures are therefore of significant interest. Relations between common information measures include the duality relations based on M\"obius inversion on lattices. These are the direct consequence of the symmetries of the lattices of the sets of variables (subsets ordered by inclusion). While the mathematical properties and relationships among these information-related measures are of significant interest, there has been, to our knowledge, no systematic examination of the full range of relationships and no unification of this diverse range of functions into a single formalism as we do here. In this paper we define operators on functions on these lattices based on the M\"obius inversion idea that map the functions into one another (M\"obius operators.) We show that these operators form a simple group isomorphic to the symmetric group S3. Relations among the set of functions on the lattice are transparently expressed in terms of the operator algebra, and, applied to the information measures, can be used to derive a wide range of relationships among measures. We describe a direct relation between sums of conditional log-likelihoods and previously defined dependency measures. The algebra is naturally generalized which yields more extensive relationships. This formalism provides a fundamental unification of information related measures, but isomorphism of all distributive lattices with the subset lattice implies broad potential application of these results.


title: Dark matter subhalos and unidentified sources in the Fermi 3FGL source catalog

abstract: If dark matter consists of weakly interacting massive particles (WIMPs), dark matter subhalos in the Milky Way could be detectable as gamma-ray point sources due to WIMP annihilation. In this work, we perform an updated study of the detectability of dark matter subhalos as gamma-ray sources with the Fermi Large Area Telescope (Fermi LAT). We use the results of the Via Lactea II simulation, scaled to the Planck 2015 cosmological parameters, to predict the local dark matter subhalo distribution. Under optimistic assumptions for the WIMP parameters --- a 40 GeV particle annihilating to $b\bar{b}$ with a thermal cross-section, as required to explain the Galactic center GeV excess --- we predict that at most $\sim 10$ subhalos might be present in the third Fermi LAT source catalog (3FGL). This is a smaller number than has been predicted by prior studies, and we discuss the origin of this difference. We also compare our predictions for the detectability of subhalos with the number of subhalo candidate sources in 3FGL, and derive upper limits on the WIMP annihilation cross-section as a function of the particle mass. If a dark matter interpretation could be excluded for all 3FGL sources, our constraints would be competitive with those found by indirect searches using other targets, such as known Milky Way satellite galaxies.

categories: astro-ph.HE astro-ph.CO astro-ph.GA hep-ph


title: Molecular Hydrogen Absorption from the Halo of a z ~ 0.4 Galaxy

abstract: Lyman- and Werner-band absorption of molecular hydrogen (H$_2$) is detected in $\sim$50% of low redshift ($z<1$) DLAs/sub-DLAs with $N$(H$_2$) > 10$^{14.4}$ cm$^{-2}$. However the true origin(s) of the H$_2$ bearing gas remain elusive. Here we report a new detection of an H$_{2}$ absorber at $z=$ 0.4298 in the HST/COS spectra of quasar PKS 2128-123. The total $N$(HI) of 10$^{19.50\pm0.15}$ cm$^{-2}$ classifies the absorber as a sub-DLA. H$_{2}$ absorption is detected up to the $J=3$ rotational level with a total $\log N$(H$_{2}$) = 16.36$\pm$0.08 corresponding to a molecular fraction of log $f$(H$_{2}$) = $-$2.84$\pm$0.17. The excitation temperature of $T_{ex}$ = 206$\pm$6K indicates the presence of cold gas. Using detailed ionization modelling we obtain a near-solar metallicity (i.e., [O/H]= $-$0.26$\pm$0.19) and a dust-to-gas ratio of $\log \kappa \sim -0.45$ for the H$_{2}$ absorbing gas. The host-galaxy of the sub-DLA is detected at an impact parameter of $\rho \sim$ 48 kpc with an inclination angle of $i \sim$ 48 degree and an azimuthal angle of $\Phi \sim$ 15 degree with respect to the QSO sightline. We show that co-rotating gas in an extended disk cannot explain the observed kinematics of Mg II absorption. Moreover, the inferred high metallicity is not consistent with the scenario of gas accretion. An outflow from the central region of the host-galaxy, on the other hand, would require a large opening angle (i.e., 2$\theta>$150 degree), much larger than the observed outflow opening angles in Seyfert galaxies, in order to intercept the QSO sightline. We thus favor a scenario in which the H$_2$ bearing gas is stemming from a dwarf-satellite galaxy, presumably via tidal and/or ram-pressure stripping. Detection of a dwarf galaxy candidate in the HST/WFPC2 image at an impact parameter of $\sim$12 kpc reinforces such an idea.


title: Properties of the Interstellar Medium in Star-Forming Galaxies at z~1.4 revealed with ALMA

abstract: We conducted observations of 12CO(J=5-4) and dust thermal continuum emission toward twenty star-forming galaxies on the main sequence at z~1.4 using ALMA to investigate the properties of the interstellar medium. The sample galaxies are chosen to trace the distributions of star-forming galaxies in diagrams of stellar mass-star formation rate and stellar mass-metallicity. We detected CO emission lines from eleven galaxies. The molecular gas mass is derived by adopting a metallicity-dependent CO-to-H2 conversion factor and assuming a CO(5-4)/CO(1-0) luminosity ratio of 0.23. Molecular gas masses and its fractions (molecular gas mass/(molecular gas mass + stellar mass)) for the detected galaxies are in the ranges of (3.9-12) x 10^{10} Msun and 0.25-0.94, respectively; these values are significantly larger than those in local spiral galaxies. The molecular gas mass fraction decreases with increasing stellar mass; the relation holds for four times lower stellar mass than that covered in previous studies, and that the molecular gas mass fraction decreases with increasing metallicity. Stacking analyses also show the same trends. The dust thermal emissions were clearly detected from two galaxies and marginally detected from five galaxies. Dust masses of the detected galaxies are (3.9-38) x 10^{7} Msun. We derived gas-to-dust ratios and found they are 3-4 times larger than those in local galaxies. The depletion times of molecular gas for the detected galaxies are (1.4-36) x 10^{8} yr while the results of the stacking analysis show ~3 x 10^{8} yr. The depletion time tends to decrease with increasing stellar mass and metallicity though the trend is not so significant, which contrasts with the trends in local galaxies.


title: On the R-Process Enrichment of Dwarf Spheroidal Galaxies

abstract: Recent observations of Reticulum II have uncovered an overabundance of r-process elements, compared to similar ultra-faint dwarf spheroidal galaxies (UFDs). Because the metallicity and star formation history of Reticulum II appear consistent with all known UFDs, the high r-process abundance of Reticulum II suggests enrichment through a single, rare event, such as a double neutron star (NS) merger. However, we note that this scenario is extremely unlikely, as binary stellar evolution models require significant supernova natal kicks to produce NS-NS or NS-black hole mergers, and these kicks would efficiently remove compact binary systems from the weak gravitational potentials of UFDs. We examine alternative mechanisms for the production of r-process elements in UFDs, including a novel mechanism wherein NSs in regions of high dark matter density implode after accumulating a black-hole-forming mass of dark matter. We find that r-process proto-material ejection by tidal forces, when a single neutron star implodes into a black hole, can occur at a rate matching the r-process abundance of both Reticulum II and the Milky Way. Remarkably, dark matter models which collapse a single neutron star in observed UFDs also solve the missing pulsar problem in the Milky Way Galactic center. We propose tests specific to dark matter r-process production which may uncover, or rule out, this model.

categories: astro-ph.HE astro-ph.CO astro-ph.GA astro-ph.SR hep-ph


title: Thin accretion disks are stabilized by a strong magnetic field

abstract: By studying three-dimensional, radiative, global simulations of sub-Eddington, geometrically thin black hole accretion flows we show that thin disks which are dominated by magnetic pressure are stable against thermal instability. Such disks are thicker than predicted by the standard model and show significant amount of dissipation inside the marginally stable orbit. Radiation released in this region, however, does not escape to infinity but is advected into the black hole. We find that the resulting accretion efficiency ($5.5\pm0.5\%$ for the simulated $0.8\dot M_{\rm Edd}$ disk) is very close to the predicted by the standard model ($5.7\%$).


title: Spectro-polarimetry confirms central powering of a Ly$\alpha$ nebula at z=3.09

abstract: We present a follow-up study to the imaging polarimetry performed by Hayes et. al. 2011 on LAB1 in the SSA22 protocluster region. Arguably the most well-known Lyman-$\alpha$ "blob", this radio-quiet emission-line nebula likely hosts a galaxy which is either undergoing significant star formation or hosts an AGN, or both. We obtain deep, spatially resolved spectro-polarimetry of the Ly$\alpha$ emission and detect integrated linear polarization of $9$-$13\%\pm2$-$3\%$ at a distance of approximately 15 kpc north and south of the peak of the Lyman-$\alpha$ surface brightness with polarization vectors lying tangential to the galactic central source. In these same regions, we also detect a wavelength dependence in the polarization which is low at the center of the Ly$\alpha$ line profile and rises substantially in the wings of the profile. These polarization signatures are easily explained by a weak out-flowing shell model. The spectral dependence of the polarization presented here provide a framework for future observations and interpretations of the southern portion of LAB1 in that any model for this system must be able to reproduce this particular spectral dependence. However, questions still remain for the northern-most spur of LAB1. In this region we detect total linear polarization of between $3$ and $20\%$ at the $5\%$ significance level. Simulations predict that polarization should increase with radius for a symmetric geometry. That the northern spur does not suggests either that this region is not symmetric (which is likely) and exhibits variations in columns density, or that it is kinematically distinct from the rest of LAB1 and powered by another mechanism altogether.


title: Tidal Disruption Event (TDE) Demographics

abstract: We survey the properties of stars destroyed in TDEs as a function of BH mass, stellar mass and evolutionary state, star formation history and redshift. For Mbh<10^7Msun, the typical TDE is due to a M*~0.3Msun M-dwarf, although the mass function is relatively flat for $M*<Msun. The contribution from older main sequence stars and sub-giants is small but not negligible. From Mbh~10^7.5-10^8.5Msun, the balance rapidly shifts to higher mass stars and a larger contribution from evolved stars, and is ultimately dominated by evolved stars at higher BH masses. The star formation history has little effect until the rates are dominated by evolved stars. TDE rates should decline very rapidly towards higher redshifts. The volumetric rate of TDEs is very high because the BH mass function diverges for low masses. However, any emission mechanism which is largely Eddington-limited for low BH masses suppresses this divergence in any observed sample and leads to TDE samples dominated by Mbh~10^6.0-10^7.5Msun BHs with roughly Eddington peak accretion rates. The typical fall back time is relatively long, with 16% having Tfb<10^(-1) years (37 days), and 84% having longer time scales. Many residual rate discrepancies can be explained if surveys are biased against TDEs with these longer Tfb, which seems very plausible if Tfb has any relation to the transient rise time. For almost any BH mass function, systematic searches for fainter, faster time scale TDEs in smaller galaxies, and longer time scale TDEs in more massive galaxies are likely to be rewarded.


title: The Gravity Dual of Renyi Entropy

abstract: A remarkable yet mysterious property of black holes is that their entropy is proportional to the horizon area. This area law inspired the holographic principle, which was later realized concretely in gauge/gravity duality. In this context, entanglement entropy is given by the area of a minimal surface in a dual spacetime. However, discussions of area laws have been constrained to entanglement entropy, whereas a full understanding of a quantum state requires Renyi entropies. Here we show that all Renyi entropies satisfy a similar area law in holographic theories and are given by the areas of dual cosmic branes. This geometric prescription is a one-parameter generalization of the minimal surface prescription for entanglement entropy. Applying this we provide the first holographic calculation of mutual Renyi information between two disks of arbitrary dimension. Our results provide a framework for efficiently studying Renyi entropies and understanding entanglement structures in strongly coupled systems and quantum gravity.


title: A Simple, Heuristic Derivation of our "No Backreaction" Results

abstract: We provide a simple discussion of our results on the backreaction effects of density inhomogeneities in cosmology, without mentioning one-parameter families or weak limits. Emphasis is placed on the manner in which "averaging" is done and the fact that one is solving Einstein's equation. The key assumptions and results that we rigorously derived within our original mathematical framework are thereby explained in a heuristic way.


title: Spending too much time at the Galactic bar: chaotic fanning of the Ophiuchus stream

abstract: The Ophiuchus stellar stream is peculiar: (1) its length is short given the age of its constituent stars, and (2) several probable member stars that lie close in both sky position and velocity have dispersions in these dimensions that far exceed those seen within the stream. The stream's proximity to the Galactic center suggests that the bar must have a significant influence on its dynamical history: The triaxiality and time-dependence of the bar may generate chaotic orbits in the vicinity of the stream that can greatly affect its morphology. We explore this hypothesis with models of stream formation along orbits consistent with Ophiuchus' properties in a Milky Way potential model that includes a rotating bar. We find that in all choices for the rotation parameters of the bar, orbits fit to the stream are strongly chaotic. Mock streams generated along these orbits qualitatively match the observed properties of the stream: because of chaos, stars stripped early generally form low-density, high-dispersion "fans" leaving only the most recently disrupted material detectable as a strong over-density. Our models predict that there should be more low-surface-brightness tidal debris than detected so far, likely with a complex phase-space morphology. The existence of or lack of these features around the Ophiuchus stream would provide an interesting constraint on the properties of the Milky Way bar and would help distinguish between formation scenarios for the stream. This is the first time that chaos has been used to explain the properties of a stellar stream and is the first demonstration of the dynamical importance of chaos in the Galactic halo. The existence of long, thin streams around the Milky Way---presumably formed along non- or weakly-chaotic orbits---may represent only a subset of the total population of disrupted satellites.


title: The stellar-to-halo mass relation of GAMA galaxies from 100 square degrees of KiDS weak lensing data

abstract: We study the stellar-to-halo mass relation of central galaxies in the range 9.7<log_10(M_*/h^-2 M_sun)<11.7 and z<0.4, obtained from a combined analysis of the Kilo Degree Survey (KiDS) and the Galaxy And Mass Assembly (GAMA) survey. We use ~100 deg^2 of KiDS data to study the lensing signal around galaxies for which spectroscopic redshifts and stellar masses were determined by GAMA. We show that lensing alone results in poor constraints on the stellar-to-halo mass relation due to a degeneracy between the satellite fraction and the halo mass, which is lifted when we simultaneously fit the stellar mass function. At M_sun>5x10^10 h^-2 M_sun, the stellar mass increases with halo mass as ~M_h^0.25. The ratio of dark matter to stellar mass has a minimum at a halo mass of 8x10^11 h^-1 M_sun with a value of M_h/M_*=56_-10^+16 [h]. We also use the GAMA group catalogue to select centrals and satellites in groups with five or more members, which trace regions in space where the local matter density is higher than average, and determine for the first time the stellar-to-halo mass relation in these denser environments. We find no significant differences compared to the relation from the full sample, which suggests that the stellar-to-halo mass relation does not vary strongly with local density. Furthermore, we find that the stellar-to-halo mass relation of central galaxies can also be obtained by modelling the lensing signal and stellar mass function of satellite galaxies only, which shows that the assumptions to model the satellite contribution in the halo model do not significantly bias the stellar-to-halo mass relation. Finally, we show that the combination of weak lensing with the stellar mass function can be used to test the purity of group catalogues.


title: Sample variance in weak lensing: how many simulations are required?

abstract: Constraining cosmology using weak gravitational lensing consists of comparing a measured feature vector of dimension $N_b$ with its simulated counterpart. An accurate estimate of the $N_b\times N_b$ feature covariance matrix $\mathbf{C}$ is essential to obtain accurate parameter confidence intervals. When $\mathbf{C}$ is measured from a set of simulations, an important question is how large this set should be. To answer this question, we construct different ensembles of $N_r$ realizations of the shear field, using a common randomization procedure that recycles the outputs from a smaller number $N_s\leq N_r$ of independent ray-tracing $N$--body simulations. We study parameter confidence intervals as a function of ($N_s,N_r$) in the range $1\leq N_s\leq 200$ and $1\leq N_r\lesssim 10^5$. Previous work has shown that Gaussian noise in the feature vectors (from which the covariance is estimated) lead, at quadratic order, to an $O(1/N_r)$ degradation of the parameter confidence intervals. Using a variety of lensing features measured in our simulations, including shear-shear power spectra and peak counts, we show that cubic and quartic covariance fluctuations lead to additional $O(1/N_r^2)$ error degradation that is not negligible when $N_r$ is only a factor of few larger than $N_b$. We study the large $N_r$ limit, and find that a single, 240Mpc$/h$ sized $512^3$-particle $N$--body simulation ($N_s=1$) can be repeatedly recycled to produce as many as $N_r={\rm few}\times10^4$ shear maps whose power spectra and high-significance peak counts can be treated as statistically independent. As a result, a small number of simulations ($N_s=1$ or $2$) is sufficient to forecast parameter confidence intervals at percent accuracy.


title: Systematic or Signal? How dark matter misalignments can bias strong lensing models of galaxy clusters

abstract: We explore how assuming that mass traces light in strong gravitational lensing models can lead to systematic errors in the predicted position of multiple images. Using a model based on the galaxy cluster MACSJ0416 (z = 0.397) from the Hubble Frontier Fields, we split each galactic halo into a baryonic and dark matter component. We then shift the dark matter halo such that it no longer aligns with the baryonic halo and investigate how this affects the resulting position of multiple images. We find for physically motivated misalignments in dark halo position, ellipticity, position angle and density profile, that multiple images can move on average by more than 0.2" with individual images moving greater than 1". We finally estimate the full error induced by assuming that light traces mass and find that this assumption leads to an expected RMS error of 0.5", almost the entire error budget observed in the Frontier Fields. Given the large potential contribution from the assumption that light traces mass to the error budget in mass reconstructions, we predict that it should be possible to make a first significant detection and characterisation of dark halo misalignments in the Hubble Frontier Fields with strong lensing. Finally, we find that it may be possible to detect ~1kpc offsets between dark matter and baryons, the smoking gun for self-interacting dark matter, should the correct alignment of multiple images be observed.