title: The transmission coefficient distribution of highly scattering sparse random media

abstract: We consider the distribution of the transmission coefficients, i.e. the singular values of the modal transmission matrix, for 2D random media with periodic boundary conditions composed of a large number of point-like non-absorbing scatterers. The scatterers are placed at random locations in the medium and have random refractive indices that are drawn from an arbitrary, known distribution. We construct a randomized model for the scattering matrix that retains scatterer dependent properties essential to reproduce the transmission coefficient distribution and analytically characterize the distribution of this matrix as a function of the refractive index distribution, the number of modes, and the number of scatterers. We show that the derived distribution agrees remarkably well with results obtained using a numerically rigorous spectrally accurate simulation. Analysis of the derived distribution provides the strongest principled justification yet of why we should expect perfect transmission in such random media regardless of the refractive index distribution of the constituent scatterers. The analysis suggests a sparsity condition under which random media will exhibit a perfect transmission-supporting universal transmission coefficient distribution in the deep medium limit.


title: Three analogues of Stern's diatomic sequence

abstract: We generalize Stern's diatomic sequence in three ways.


title: Bicriteria Optimization of Technological Parameters in Algorithm for Designing Magnetic Composites

abstract: Novel algorithm for designing values of technological parameters for production of Soft Magnetic Composites (SMC) has been created. These parameters are the following magnitudes: hardening temperature $T$ and compaction pressure $p$. They enable us to optimize of power losses and induction. The advantage of the presented algorithm consists in the bicriteria optimization. The crucial role in the presented algorithm play scaling and notion of pseudo-state equation. On the base of these items the mathematical models of the power losses and induction have been created. The models parameters have been calculated on the basis of the power losses characteristics and hysteresis loops. The created optimization system has been applied to specimens of Somaloy 500. Obtained output consists of finite set of feasible solutions. In order to select unique solution an example of additional criterion has been formulated.


title: Covariant Perturbations of f(R) Black Holes: The Weyl Terms

abstract: In this paper we revisit non-spherical perturbations of the Schwarzschild black hole in the context of $f(R)$ gravity. Previous studies were able to demonstrate the stability of the $f(R)$ Schwarzschild black hole against gravitational perturbations in both the even and odd parity sectors. In particular, it was seen that the Regge-Wheeler and Zerilli equations in $f(R)$ gravity obey the same equations as their General Relativity counterparts. More recently, the 1+1+2 semi-tetrad formalism has been used to derive a set of two wave equations: one for transverse, trace-free (tensor) perturbations and one for the additional scalar modes that characterise fourth-order theories of gravitation. The master variable governing tensor perturbations was shown to be a modified Regge-Wheeler tensor obeying the same equation as in General Relativity. However, it is well known that there is a non-uniqueness in the definition of the master variable. In this paper we derive a set of two perturbation variables and their concomitant wave equations that describe gravitational perturbations in a covariant and gauge invariant manner. These variables can be related to the Newman-Penrose (NP) Weyl scalars as well as the master variables from the 2+2 formalism.


title: GRChombo : Numerical Relativity with Adaptive Mesh Refinement

abstract: In this work, we introduce GRChombo: a new numerical relativity code which incorporates full adaptive mesh refinement (AMR) using block structured Berger-Rigoutsos grid generation. The code supports non-trivial "many-boxes-in-many-boxes" mesh hierarchies and massive parallelism through the Message Passing Interface (MPI). GRChombo evolves the Einstein equation using the standard BSSN formalism, with an option to turn on CCZ4 constraint damping if required. The AMR capability permits the study of a range of new physics which has previously been computationally infeasible in a full 3+1 setting, whilst also significantly simplifying the process of setting up the mesh for these problems. We show that GRChombo can stably and accurately evolve standard spacetimes such as binary black hole mergers and scalar collapses into black holes, demonstrate the performance characteristics of our code, and discuss various physics problems which stand to benefit from the AMR technique.


title: Microscopic analysis of fusion hindrance in heavy systems

abstract: Background: Heavy-ion fusion reactions involving heavy nuclei at energies around the Coulomb barrier exhibit fusion hindrance, where the probability of compound nucleus formation is strongly hindered compared with that in light- and medium-mass systems. The origin of this fusion hindrance has not been well understood from a microscopic point of view. Purpose: Analyze the fusion dynamics in heavy systems by a microscopic reaction model and understand the origin of the fusion hindrance. Method: We employ the time-dependent Hartree-Fock (TDHF) theory. We extract nucleus--nucleus potential and energy dissipation by the method combining TDHF dynamics of the entrance channel of fusion reactions with one-dimensional Newton equation including a dissipation term. Then, we analyze the origin of the fusion hindrance using the properties of the extracted potential and energy dissipation. Results: Extracted potentials show monotonic increase as the relative distance of two nuclei decreases, which induces the disappearance of an ordinary barrier structure of the potential. This is different from those in light- and medium-mass systems and from density-constraint TDHF calculations. Extracted friction coefficients show sizable energy dependence and universal value of their magnitude, which are rather similar to those in light- and medium-mass systems. Using these properties, we analyze the origin of the fusion hindrance and find that contribution of the increase in potential to the extra-push energy is larger than that of the accumulated dissipation energy in most systems studied in this article. Conclusions: By the analysis of the origin of the fusion hindrance, we conclude that, as the system becomes heavier, the dynamical increase in potential at small relative distances plays a more important role than the dissipation during the fusion reaction for understanding the origin of the fusion hindrance.


title: A mathematical motivation for complex-valued convolutional networks

abstract: A complex-valued convolutional network (convnet) implements the repeated application of the following composition of three operations, recursively applying the composition to an input vector of nonnegative real numbers: (1) convolution with complex-valued vectors followed by (2) taking the absolute value of every entry of the resulting vectors followed by (3) local averaging. For processing real-valued random vectors, complex-valued convnets can be viewed as "data-driven multiscale windowed power spectra," "data-driven multiscale windowed absolute spectra," "data-driven multiwavelet absolute values," or (in their most general configuration) "data-driven nonlinear multiwavelet packets." Indeed, complex-valued convnets can calculate multiscale windowed spectra when the convnet filters are windowed complex-valued exponentials. Standard real-valued convnets, using rectified linear units (ReLUs), sigmoidal (for example, logistic or tanh) nonlinearities, max. pooling, etc., do not obviously exhibit the same exact correspondence with data-driven wavelets (whereas for complex-valued convnets, the correspondence is much more than just a vague analogy). Courtesy of the exact correspondence, the remarkably rich and rigorous body of mathematical analysis for wavelets applies directly to (complex-valued) convnets.


title: Possibility of higher dimensional anisotropic compact star

abstract: We have provided here a new class of interior solutions for anisotropic stars admitting conformal motion in higher dimensional noncommutative spacetime. The Einstein fields equations are solved by choosing a particular density distribution function of Lorentzian type \cite{Nozari} under noncommutative geometry. Several cases with dimensions $4D$ and higher, e.g. $5D$, $6D$ and $11D$ have been discussed separately. An overall observation is that the model parameters, such as density, radial pressure, transverse pressure, anisotropy all are well behaved and represent a compact star with radius $4.17$ km. However, emphasis has been given on the acceptability of the model from physical point of view. As a consequence it is observed that higher dimensions, i.e. beyond $4D$ spacetime, exhibit several interesting yet bizarre features which are not at all untenable for a compact stellar model of strange quark type and thus dictates a possibility of its extra dimensional existence.


title: Fano resonance and spectrally modified photoluminescence enhancement in monolayer MoS2 integrated with plasmonic nanoantenna array

abstract: The manipulation of light-matter interactions in two-dimensional atomically thin crystals is critical for obtaining new optoelectronic functionalities in these strongly confined materials. Here, by integrating chemically grown monolayers of MoS2 with a silver-bowtie nanoantenna array supporting narrow surface-lattice plasmonic resonances, a unique two-dimensional optical system has been achieved. The enhanced exciton-plasmon coupling enables profound changes in the emission and excitation processes leading to spectrally tunable, large photoluminescence enhancement as well as surface-enhanced Raman scattering at room temperature. Furthermore, at low temperatures, due to the decreased damping of MoS2 excitons interacting with the plasmonic resonances of the bowtie array, stronger exciton-plasmon coupling is achieved resulting in a Fano lineshape in the reflection spectrum. The Fano lineshape, which is due to the interference between the pathways involving the excitation of the exciton and plasmon, can be tuned by altering the coupling strengths between the two systems via changing the design of the bowties lattice. The ability to manipulate the optical properties of two-dimensional systems with tunable plasmonic resonators offers a new platform for the design of novel optical devices with precisely tailored responses.


title: A Density Tur\'an Theorem

abstract: Let $F$ be a graph which contains an edge whose deletion reduces its chromatic number. For such a graph $F,$ a classical result of Simonovits from 1966 shows that every graph on $n\ge n_0(F)$ vertices with more than $\frac{\chi(F)-2}{\chi(F)-1}\cdot \frac{n^2}{2}$ edges contains a copy of $F$. In this paper we derive a similar theorem for multipartite graphs. For a graph $H$ and an integer $\ell \geq v(H)$, let $d_{\ell}(H)$ be the minimum real number such that every $\ell$-partite graph whose edge density between any two parts is greater than $d_{\ell}(H)$ contains a copy of $H$. Our main contribution is to show that $d_{\ell}(H)=\frac{\chi(H)-2}{\chi(H)-1}$ for $\ell \ge \ell_0(H)$ sufficiently large if and only if $H$ admits a vertex-colouring with $\chi(H)-1$ colours such that all colour classes but one are independent sets, and the exceptional class induces just a matching. When $H$ is a clique, this recovers a result of Pfender [Complete subgraphs in multipartite graphs, Combinatorica 32 (2012), 483--495]. We also consider several extensions of Pfender's result.


title: Magnetic susceptibility, magnetization, magnetic moment and characterization of Carancas meteorite

abstract: On September, 15th, 2007, in the community of Carancas (Puno, Peru) a stony meteorite formed a crater explosive type with a mean diameter of 13.5 m. some samples meteorite fragments were collected. The petrologic analysis performed corresponds to a meteorite ordinary chondrite H 4-5. In this paper we have analyzed the magnetic properties of a meteorite fragment with a proton magnetometer. Also in order to have a complete characterization of the Carancas meteorite and its crater, from several papers, articles and reports, we have made a compilation of the most important characteristics and properties of this meteorite.


title: The pseudoarc is a co-existentially closed continuum

abstract: Answering a question of P. Bankston, we show that the pseudoarc is a co-existentially closed continuum. We also show that $C(X)$, for $X$ a nondegenerate continuum, can never have quantifier elimination, answering a question of the the first and third named authors and Farah and Kirchberg.


title: Inert Scalar Doublet Asymmetry as Origin of Dark Matter

abstract: In the inert scalar doublet framework, we analyze what would be the effect of a $B-L$ asymmetry that could have been produced in the Universe thermal bath at high temperature. We show that, unless the "$\lambda_5$" scalar interaction is tiny, this asymmetry is automatically reprocessed in part into an inert scalar asymmetry that could be at the origin of dark matter today. Along this scenario, the inert mass scale lies in the few-TeV range and direct detection constraints require that the inert scalar particles decay into a lighter dark matter particle which, as the inert doublet, is odd under a $Z_2$ symmetry.


title: A mathematically assisted reconstruction of the initial focus of the yellow fever outbreak in Buenos Aires (1871)

abstract: We discuss the historic mortality record corresponding to the initial focus of the yellow fever epidemic outbreak registered in Buenos Aires during the year 1871 as compared to simulations of a stochastic population dynamics model. This model incorporates the biology of the urban vector of yellow fever, the mosquito Aedes aegypti, the stages of the disease in the human being as well as the spatial extension of the epidemic outbreak. After introducing the historical context and the restrictions it puts on initial conditions and ecological parameters, we discuss the general features of the simulation and the dependence on initial conditions and available sites for breeding the vector. We discuss the sensitivity, to the free parameters, of statistical estimators such as: final death toll, day of the year when the outbreak reached half the total mortality and the normalized daily mortality, showing some striking regularities. The model is precise and accurate enough to discuss the truthfulness of the presently accepted historic discussions of the epidemic causes, showing that there are more likely scenarios for the historic facts.


title: Extremal quantum states and their Majorana constellations

abstract: The characterization of quantum polarization of light requires knowledge of all the moments of the Stokes variables, which are appropriately encoded in the multipole expansion of the density matrix. We look into the cumulative distribution of those multipoles and work out the corresponding extremal pure states. We find that SU(2) coherent states are maximal to any order whereas the converse case of minimal states (which can be seen as the most quantum ones) is investigated for a diverse range of the number of photons. Taking advantage of the Majorana representation, we recast the problem as that of distributing a number of points uniformly over the surface of the Poincare sphere.


title: On the Nature of Rapidly Rotating Single Evolved Stars

abstract: We present an analysis of the nature of the rapidly rotating, apparently single giant based on rotational and radial velocity measurements carried out by the CORAVEL spectrometers. From the analyzed sample, composed of 2010 spectroscopic, apparently single, evolved stars of luminosity classes IV, III, II, and Ib with spectral types G and K, we classified 30 stars that presented unusual, moderate to rapid rotation. This work reports, for the first time, the presence of these abnormal rotators among subgiant, bright giant, and Ib supergiant stars. To date, this class of stars was reported only among giant stars of luminosity class III. Most of these abnormal rotators present an IRAS infrared excess, which, in principle, can be related to dust around these stars.


title: GAC, savings, and unbounded inputs

abstract: Let a control system and a target be given on an open subset of an Euclidean space. The existence of a Control Lyapunov Function - namely a positive definite, semiconcave, solution of the Hamilton-Jacobi inequality corresponding to the control vector field -- guarantees Global Asymptotic Controllability (GAC). In this case, however, minimization is not an issue. Instead, if a Lagrangean with non-negative values is considered as well, an optimal control problem can be defined in relation to the corresponding integral functional. In the first part of the present paper we show that the existence of a Minimum Restraint Function -- a solution of a strict Hamilton-Jacobi inequality involving the Lagrangian and a non-negative "savings multiplier" -- provides not only global asymptotic controllability but also savings, namely a state-dependent upper bound for the infima. This extends a former result, where the control set was assumed to be compact. Here we allow unbounded controls and replace inputs' values' compactness with a quite mild hypothesis concerning the dependence of the data on inputs: such condition is met, for instance, by control vector fields that are compositions of Lipschitz maps with polynomials and exponentials of the control variable. In the second part of the paper we focus on the case when the dynamics is a polynomial in the control variable. Through some analysis of convexity properties of vector-valued polynomials' ranges, we prove some simplified versions of the main result, in terms of either affine representability or reduction to weak subsystems for the original dynamics.


title: Two-User Erasure Interference Channels with Local Delayed CSIT

abstract: We study the capacity region of two-user erasure interference channels with local delayed channel state information at the transmitters. In our model, transmitters have local mismatched outdated knowledge of the channel gains. We propose a transmission strategy that only relies on the delayed knowledge of the {\it outgoing} links at each transmitter and achieves the outer-bound for the scenario in which transmitters learn the entire channel state with delay. Our result reveals the subset of the channel state information that affects the capacity region the most. We also identify cases in which local delayed knowledge of the channel state does not provide any gain over the zero knowledge assumption. To do so, we revisit a long-known intuition about interference channels that as long as the marginal distributions at the receivers are conserved, the capacity remains the same. We take this intuition and impose a certain spatial correlation among channel gains such that the marginal distributions remain unchanged. Then we provide an outer-bound on the capacity region of the channel with correlation that matches the capacity region when transmitters do not have access to channel state information.


title: Low-energy constants and condensates from ALEPH hadronic $\tau$ decay data

abstract: We determine the NLO chiral low-energy constant $L_{10}^r$, and combinations $C_{12}^r\pm C_{61}^r+C_{80}^r$, $C_{13}^r-C_{62}^r+C_{81}^r$, $C_{61}^r$, and $C_{87}^r$, of the NNLO chiral low-energy constants incorporating recently revised ALEPH results for the non-strange vector ($V$) and axial-vector ($A$) hadronic $\tau$ decay distributions and recently updated RBC/UKQCD lattice data for the non-strange $V-A$ two-point function. In the $\bar{\rm MS}$ scheme, at renormalization scale $\mu=770\ {MeV}$, we find $L_{10}^r=-0.00350(17)$, $C_{12}^r+C_{61}^r+C_{80}^r=0.00237(16)\ {GeV}^{-2}$, $C_{12}^r-C_{61}^r+C_{80}^r=-0.00056(15)\ {GeV}^{-2}$, $C_{13}^r-C_{62}^r+C_{81}^r=0.00046(9)\ {GeV}^{-2}$, $C_{61}^r=0.00146(15)\ {GeV}^{-2}$, and $C_{87}^r=0.00510(22)\ {GeV}^{-2}$. With errors here at or below the level expected for contributions of yet higher order in the chiral expansion, the analysis exhausts the possibilities of what can be meaningfully achieved in an NNLO analysis. We also consider the dimension six and dimension eight coefficients in the operator product expansion in the $V-A$ channel.


title: Exponential Family Models from Bayes' Theorem under Expectation Constraints

abstract: It is shown that a consistent application of Bayesian updating from a prior probability density to a posterior using evidence in the form of expectation constraints leads to exactly the same results as the application of the maximum entropy principle, namely a posterior belonging to the exponential family. The Bayesian updating procedure presented in this work is not expressed as a variational principle, and does not involve the concept of entropy. Therefore it conceptually constitutes a complete alternative to entropic methods of inference.


title: Design and implementation of an Android application (MobilitApp+) to analyze the mobility patterns of citizens in the Metropolitan Region of Barcelona

abstract: In our project we have designed an Android application to obtain mobility data of the citizens in the metropolitan area of Barcelona. Our implementation synchronously obtains in background on the one hand, periodic location updates and, on the other hand, the type of activity citizens are doing. At the end of the day, all this data is processed and sent to a server where are stored to obtain mobility patterns from citizens that could help to improve the current transportation infrastructure. MobilitApp is fully functional and stable although the results can be improved in some situations. In future releases we will implement machine learning technics to obtain significant improvements, especially in the activity recognition modules.


title: On the Variability Estimation of Lognormal Distribution Based on Sample Harmonic and Arithmetic Means

abstract: For the lognormal distribution, an unbiased estimator of the squared coefficient of variation is derived from the relative ratio of sample arithmetic to harmonic means. Analytical proofs and simulation results are presented.


title: Spectral solution of urn models for interacting particle systems

abstract: Using generating function methods for diagonalizing the transition matrix in 2-Urn models, we provide a complete classification into solvable and unsolvable subclasses, with further division of the solvable models into the Martingale and non-Martingale subcategories, and prove that the stationary distribution is a Gaussian function in the latter. We also give a natural condition related to the symmetry of the random walk in which the non-Martingale Urn models lead to an increase in entropy from Gaussian states. The condition also shows that universal symmetry in the macro-state is equivalent to increasing entropy. Certain models of social opinion dynamics, treated as Urn models, do not increase in entropy, unlike isolated mechanical systems.


title: Thermodynamic universality of quantum Carnot engines

abstract: The Carnot statement of the second law of thermodynamics poses an upper limit on the efficiency of all heat engines. Recently, it has been studied whether generic quantum features such as coherence and quantum entanglement could allow for quantum devices with efficiencies larger than the Carnot efficiency. The present study shows that this is not permitted by the laws of thermodynamics. In particular, we will show that rather the definition of heat has to be modified to account for the thermodynamic cost for maintaining coherence and entanglement. Our theoretical findings are numerically illustrated for an experimentally relevant example from optomechanics.


title: `Holographic' treatment of surface disorder on a topological insulator

abstract: The effect of surface disorder on electronic systems is particularly interesting for topological phases with surface and edge states. Using exact diagonalization, it has been demonstrated that the surface states of a 3D topological insulator survive strong surface disorder, and simply get pushed to a clean part of the bulk. Here we explore a new method which analytically eliminates the clean bulk, and reduces a $D$-dimensional problem to a Hamiltonian-diagonalization problem within the $(D-1)$-dimensional disordered surface. This dramatic reduction in complexity allows the analysis of significantly bigger systems than is possible with exact diagonalization. We use our method to analyze a 2D topological spin-Hall insulator with non-magnetic and magnetic edge impurities, and we calculate the probability density (or local density of states) of the zero-energy eigenstates as a function of edge-parallel momentum and layer index. Our analysis reveals that the system size needed to reach behavior in the thermodynamic limit increases with disorder. We also compute the edge conductance as a function of disorder strength, and chart a lower bound for the length scale marking the crossover to the thermodynamic limit.


title: Classical to quantum correspondence in dissipative directed transport

abstract: We compare the quantum and classical properties of the (Quantum) Isoperiodic Stable Structures -- (Q)ISSs -- which organize the parameter space of a paradigmatic dissipative ratchet model, i.e. the dissipative modified kicked rotator. We study the spectral behavior of the corresponding classical Perron-Frobenius operators with thermal noise and the quantum superoperators without it for small $\hbar_{\rm eff}$ values. We find a remarkable similarity between the classical and quantum spectra. This finding significantly extends previous results -- obtained for the mean currents and asymptotic distributions only -- and on the other hand unveils a classical to quantum correspondence mechanism where the classical noise is qualitatively different from the quantum one. This is crucial not only for simple attractors but also for chaotic ones, where just analyzing the asymptotic distribution reveals insufficient. Moreover, we provide with a detailed characterization of relevant eigenvectors by means of the corresponding Weyl-Wigner distributions, in order to better identify similarities and differences. Finally, this model being generic, it allows us to conjecture that this classical to quantum correspondence mechanism is a universal feature of dissipative systems.


title: Electron-beam evaporated cobalt films on molecular beam epitaxy prepared GaAs(001)

abstract: We have deposited Co films on the GaAs(001) surface by using an e-beam evaporation method. The thicknesses of the Co films are measured by using x-ray reflectivity and Rutherford backscattering. The magnetic properties of the films have been measured using superconducting quantum interference device. The magnetization of the films was found to decrease with increasing film thickness. The slight degradation of magnetic properties is attributed to increasing roughness on the Co surface or the Co/GaAs interface during the Co deposition.


title: Robust generation of elementary flux modes

abstract: Elementary flux modes (EFMs) are vectors defined from a metabolic reaction network, giving the connections between substrates and products. EFMs-based metabolic flux analysis (MFA) estimates the flux over each EFM from external flux measurements through least-squares data fitting. In previous work we presented an optimization method of column generation type that facilitates EFMs-based MFA when the metabolic reaction network is so large that enumerating all EFMs is prohibitive. In this work we extend this model by including errors on measurements in a robust optimization framework. In the robust optimization problem, the least-squares data fitting is minimized subject to the error on each metabolite being as unfavourable as it can be, within a given interval. In general, inclusion of robustness may make the optimization problem significantly harder. However, we show that in our case the robust problem can be stated as a convex quadratic programming problem, i.e., of the same form as the original non-robust problem. Additionally, we demonstrate that the column-generation technique of the non-robust problem can be extended also to the robust problem. Furthermore, the option to indicate intervals on metabolites that are not measured is introduced in this column generation framework. The effect of including robustness in the model is evaluated in a case-study, which indicated that the solutions of our non-robust problems are in fact near-optimal also when robustness is considered. On the other hand, the addition of intervals on unmeasured metabolites resulted in a change of optimal solution. Implying that the inclusion of intervals on unmeasured metabolites is more important than the direct consideration of measurement errors, this despite up to 20% errors.


title: Efficient Service Broker Policy For Large-Scale Cloud Environments

abstract: Algorithms, policies, and methodologies are necessary to achieve high user satisfaction and practical utilization in cloud computing by ensuring the efficient and fair allocation of every computing resource. Whenever a new job arrives in cloud environments, the service broker is responsible for selecting the data center that will execute that job. Selecting data centers serves an important function in enhancing the performance of a cloud environment. This study proposes a new service broker policy for large-scale cloud applications based on the round-robin algorithm. The proposed policy is implemented and evaluated using a CloudAnalyst simulator. It is then compared with three existing policies in terms of overall average response time by using different virtual machine load balancing algorithms. Simulation results show that the proposed policy improves the overall average response time relative to that of the other policies.


title: Skew polynomial rings over abelian and idempotent reflexive rings

abstract: Let $R$ be a ring and $\sigma$ an endomorphism of $R$. In this note, we study skew polynomial rings and skew power series rings over idempotent reflexive rings and abelian rings. Also, we introduce the concept of right (resp., left) $\sigma$-idempotent reflexive rings which generalizes right (resp., left) idempotent reflexive rings and $\sigma$-abelian rings. Certain results are obtained as corollaries from our results.


title: On the zone of a circle in an arrangement of lines

abstract: Let $\mathcal L$ be a set of $n$ lines in the plane, and let $C$ be a convex curve in the plane, like a circle or a parabola. The "zone" of $C$ in $\mathcal L$, denoted $\mathcal Z(C,\mathcal L)$, is defined as the set of all cells in the arrangement $\mathcal A(\mathcal L)$ that are intersected by $C$. Edelsbrunner et al. (1992) showed that the complexity (total number of edges or vertices) of $\mathcal Z(C,\mathcal L)$ is at most $O(n\alpha(n))$, where $\alpha$ is the inverse Ackermann function. They did this by translating the sequence of edges of $\mathcal Z(C,\mathcal L)$ into a sequence $S$ that avoids the subsequence $ababa$. Whether the worst-case complexity of $\mathcal Z(C,\mathcal L)$ is only linear is a longstanding open problem. Since the relaxation of the problem to pseudolines does have a $\Theta(n\alpha(n))$ bound, any proof of $O(n)$ for the case of straight lines must necessarily use geometric arguments. In this paper we present some such geometric arguments. We show that, if $C$ is a circle, then certain configurations of straight-line segments with endpoints on $C$ are impossible. In particular, we show that there exists a Hart-Sharir sequence that cannot appear as a subsequence of $S$. The Hart-Sharir sequences are essentially the only known way to construct $ababa$-free sequences of superlinear length. Hence, if it could be shown that every family of $ababa$-free sequences of superlinear-length eventually contains all Hart-Sharir sequences, it would follow that the complexity of $\mathcal Z(C,\mathcal L)$ is $O(n)$ whenever $C$ is a circle.


title: Towards the Design and Implementation of Aspect-Oriented Programming for Spreadsheets

abstract: A spreadsheet usually starts as a simple and single-user software artifact, but, as frequent as in other software systems, quickly evolves into a complex system developed by many actors. Often, different users work on different aspects of the same spreadsheet: while a secretary may be only involved in adding plain data to the spreadsheet, an accountant may define new business rules, while an engineer may need to adapt the spreadsheet content so it can be used by other software systems. Unfortunately, spreadsheet systems do not offer modular mechanisms, and as a consequence, some of the previous tasks may be defined by adding intrusive "code" to the spreadsheet. In this paper we go through the design and implementation of an aspect-oriented language for spreadsheets so that users can work on different aspects of a spreadsheet in a modular way. For example, aspects can be defined in order to introduce new business rules to an existing spreadsheet, or to manipulate the spreadsheet data to be ported to another system. Aspects are defined as aspect-oriented program specifications that are dynamically woven into the underlying spreadsheet by an aspect weaver. In this aspect-oriented style of spreadsheet development, different users develop, or reuse, aspects without adding intrusive code to the original spreadsheet. Such code is added/executed by the spreadsheet weaving mechanism proposed in this paper.


title: Curvilinear integral theorems for monogenic functions in commutative associative algebras

abstract: We consider an arbitrary finite-dimensional commutative associative algebra, $\mathbb{A}_n^m$, with unit over the field of complex number with $m$ idempotents. Let $e_1=1,e_2,e_3$ be elements of $\mathbb{A}_n^m$ which are linearly independent over the field of real numbers. We consider monogenic (i.e. continuous and differentiable in the sense of Gateaux) functions of the variable $xe_1+ye_2+ze_3$, where $x,y,z$ are real. For mentioned monogenic function we prove curvilinear analogues of the Cauchy integral theorem, the Morera theorem and the Cauchy integral formula.


title: Faster 64-bit universal hashing using carry-less multiplications

abstract: Intel and AMD support the Carry-less Multiplication (CLMUL) instruction set in their x64 processors. We use CLMUL to implement an almost universal 64-bit hash family (CLHASH). We compare this new family with what might be the fastest almost universal family on x64 processors (VHASH). We find that CLHASH is at least 60% faster. We also compare CLHASH with a popular hash function designed for speed (Google's CityHash). We find that CLHASH is 40% faster than CityHash on inputs larger than 64 bytes and just as fast otherwise.


title: Self-Consistent Projection Operator Theory in Nonlinear Quantum Optical Systems: A case study on Degenerate Optical Parametric Oscillators

abstract: Nonlinear quantum optical systems are of paramount relevance for modern quantum technologies, as well as for the study of dissipative phase transitions. Their nonlinear nature makes their theoretical study very challenging and hence they have always served as great motivation to develop new techniques for the analysis of open quantum systems. In this article we apply the recently developed self-consistent projection operator theory to the degenerate optical parametric oscillator to exemplify its general applicability to quantum optical systems. We show that this theory provides an efficient method to calculate the full quantum state of each mode with high degree of accuracy, even at the critical point. It is equally successful in describing both the stationary limit and the dynamics, including regions of the parameter space where the numerical integration of the full problem is significantly less efficient. We further develop a Gaussian approach consistent with our theory, which yields sensibly better results than the previous Gaussian methods developed for this system, most notably standard linearization techniques.


title: Multigrid with rough coefficients and Multiresolution operator decomposition from Hierarchical Information Games

abstract: We introduce a near-linear complexity (geometric and meshless/algebraic) multigrid/multiresolution method for PDEs with rough ($L^\infty$) coefficients with rigorous a-priori accuracy and performance estimates. The method is discovered through a decision/game theory formulation of the problems of (1) identifying restriction and interpolation operators (2) recovering a signal from incomplete measurements based on norm constraints on its image under a linear operator (3) gambling on the value of the solution of the PDE based on a hierarchy of nested measurements of its solution or source term. The resulting elementary gambles form a hierarchy of (deterministic) basis functions of $H^1_0(\Omega)$ (gamblets) that (1) are orthogonal across subscales/subbands with respect to the scalar product induced by the energy norm of the PDE (2) enable sparse compression of the solution space in $H^1_0(\Omega)$ (3) induce an orthogonal multiresolution operator decomposition. The operating diagram of the multigrid method is that of an inverted pyramid in which gamblets are computed locally (by virtue of their exponential decay), hierarchically (from fine to coarse scales) and the PDE is decomposed into a hierarchy of independent linear systems with uniformly bounded condition numbers. The resulting algorithm is parallelizable both in space (via localization) and in bandwith/subscale (subscales can be computed independently from each other). Although the method is deterministic it has a natural Bayesian interpretation under the measure of probability emerging (as a mixed strategy) from the information game formulation and multiresolution approximations form a martingale with respect to the filtration induced by the hierarchy of nested measurements.


title: Algorithms and Properties for Positive Symmetrizable Matrices

abstract: Matrices are the most common representations of graphs. They are also used for the representation of algebras and cluster algebras. This paper shows some properties of matrices in order to facilitate the understanding and locating symmetrizable matrices with specific characteristics, called positive quasi-Cartan companion matrices. Here, symmetrizable matrix are those which are symmetric when multiplied by a diagonal matrix with positive entries called symmetrizer matrix. Four algorithms are developed: one to decide whether there is a symmetrizer matrix; second to find such symmetrizer matrix; another to decide whether the matrix is positive or not; and the last to find a positive quasi-Cartan companion matrix, if there exists. The third algorithm is used to prove that the problem to decide if a matrix has a positive quasi-Cartan companion is NP.


title: HATS-13b and HATS-14b: two transiting hot Jupiters from the HATSouth survey

abstract: We report the discovery of HATS-13b and HATS-14b, two hot-Jupiter transiting planets discovered by the HATSouth survey. The host stars are quite similar to each other (HATS-13: V = 13.9 mag, M* = 0.96 Msun, R* = 0.89 Rsun, Teff = 5500 K, [Fe/H] = 0.05; HATS-14: V = 13.8 mag, M* = 0.97 Msun, R* = 0.93 Rsun, Teff = 5350 K, [Fe/H] = 0.33) and both the planets orbit around them with a period of roughly 3 days and a separation of roughly 0.04 au. However, even though they are irradiated in a similar way, the physical characteristics of the two planets are very different. HATS-13b, with a mass of Mp = 0.543 MJ and a radius of Rp = 1.212 RJ, appears as an inflated planet, while HATS-14b, having a mass of Mp = 1.071 MJ and a radius of Rp = 1.039 RJ, is only slightly larger in radius than Jupiter.


title: Analytic results for the Casimir free energy between ferromagnetic metals

abstract: We derive perturbation analytic expressions for the Casimir free energy and entropy between two dissimilar ferromagnetic plates which are applicale at arbitrarily low temperature. The dielectric properties of metals are described using either the nondissipative plasma model or the Drude model taking into account the dissipation of free charge carriers. Both cases of constant and frequency-dependent magnetic permeability are considered. It is shown that for ferromagnetic metals described by the plasma model the Casimir entropy goes to zero when the temperature vanishes, i.e., the Nernst heat theorem is satisfied. For ferromagnetic metals with perfect crystal lattices described by the Drude model the Casimir entropy goes to a nonzero constant depending on the parameters of a system with vanishing temperature, i.e., the Nernst heat theorem is violated. This constant can be positive which is quite different from the earlier investigated case of two nonmagnetic metals.


title: Are narrow unflavoured mesons a signature of new physics?

abstract: New data of narrow low mass unflavoured mesonic structures are presented. A table of these exotic masses is obtained adding previously published data. The mass sequence shows a significant coupling of some of these clusters with stable hadrons: pion, nucleon, and deuteron. Indeed this coupling allows to reproduce rather well the masses of exotic narrow baryons and dibaryons. A discussion is presented to suggest a possible interpretation of these exotic hadronic structures.


title: High frequency nano-optomechanical disk resonators in liquids

abstract: Vibrating nano- and micromechanical resonators have been the subject of research aiming at ultrasensitive mass sensors for mass spectrometry, chemical analysis and biomedical diagnosis. Unfortunately, their merits diminish dramatically in liquids due to dissipative mechanisms like viscosity and acoustic losses. A push towards faster and lighter miniaturized nanodevices would enable improved performances, provided dissipation was controlled and novel techniques were available to efficiently drive and read-out their minute displacement. Here we report on a nano-optomechanical approach to this problem using miniature semiconductor disks. These devices combine mechanical motion at high frequency above the GHz, ultra-low mass of a few picograms, and moderate dissipation in liquids. We show that high-sensitivity optical measurements allow to direct resolve their thermally driven Brownian vibrations, even in the most dissipative liquids. Thanks to this novel technique, we experimentally, numerically and analytically investigate the interaction of these resonators with arbitrary liquids. Nano-optomechanical disks emerge as probes of rheological information of unprecedented sensitivity and speed, opening applications in sensing and fundamental science.


title: Characterizations of Symmetrized Polydisc

abstract: Let $\Gamma_n$, $n \geq 2$, denote the symmetrized polydisc in $\mathbb{C}^n$, and $\Gamma_1$ be the closed unit disc in $\mathbb{C}$. We provide some characterizations of elements in $\Gamma_n$. In particular, an element $(s_1, \ldots, s_{n-1}, p) \in \mathbb{C}^n$ is in $\Gamma_n$ if and only if $s_j = \beta_j + \overline{\beta_{n-j}} p$, $j = 1, \ldots, n-1$, for some $(\beta_1, \ldots, \beta_{n-1}) \in \Gamma_{n-1}$, and $|p| \leq 1$.


title: Bohr-like black holes

abstract: The idea that black holes (BHs) result in highly excited states representing both the "hydrogen atom" and the "quasi-thermal emission" in quantum gravity is today an intuitive but general conviction. In this paper it will be shown that such an intuitive picture is more than a picture. In fact, we will discuss a model of quantum BH somewhat similar to the historical semi-classical model of the structure of a hydrogen atom introduced by Bohr in 1913. The model is completely consistent with existing results in the literature, starting from the celebrated result of Bekenstein on the area quantization.


title: Water Formation During the Epoch of First Metal Enrichment

abstract: We demonstrate that high abundances of water vapor could have existed in extremely low metallicity ($10^{-3}$ solar) partially shielded gas, during the epoch of first metal enrichment of the interstellar medium of galaxies at high redshifts.


title: A coordinated X-ray and Optical Campaign of the Nearest Massive Eclipsing Binary, $\delta$ Orionis Aa: IV. A multiwavelength, non-LTE spectroscopic analysis

abstract: Eclipsing systems of massive stars allow one to explore the properties of their components in great detail. We perform a multi-wavelength, non-LTE analysis of the three components of the massive multiple system $\delta$ Ori A, focusing on the fundamental stellar properties, stellar winds, and X-ray characteristics of the system. The primary's distance-independent parameters turn out to be characteristic for its spectral type (O9.5 II), but usage of the ${\rm \it Hipparcos}$ parallax yields surprisingly low values for the mass, radius, and luminosity. Consistent values follow only if $\delta$ Ori lies at about twice the ${\rm \it Hipparcos}$ distance, in the vicinity of the $\sigma$-Orionis cluster. The primary and tertiary dominate the spectrum and leave the secondary only marginally detectable. We estimate the V-band magnitude difference between primary and secondary to be $\Delta V \approx 2.\!\!^{\rm m}8$. The inferred parameters suggest the secondary is an early B-type dwarf ($\approx$ B1 V), while the tertiary is an early B-type subgiant ($\approx$ B0 IV). We find evidence for rapid turbulent velocities ($\sim 200$ km ${\rm s}^{-1}$) and wind inhomogeneities, partially optically thick, in the primary's wind. The bulk of the X-ray emission likely emerges from the primary's stellar wind ($\log L_{\text{X}} / L_{\text{Bol}} \approx -6.85$), initiating close to the stellar surface at $R_0 \sim 1.1\,R_*$. Accounting for clumping, the mass-loss rate of the primary is found to be $\log \dot{M} \approx -6.4\,[M_\odot\,{\rm yr}^{-1}]$, which agrees with hydrodynamic predictions, and provides a consistent picture along the X-ray, UV, optical and radio spectral domains.


title: Addressing the LHC flavour anomalies with horizontal gauge symmetries

abstract: We study the impact of an additional $U(1)'$ gauge symmetry with flavour-dependent charges for quarks and leptons on the LHC flavour anomalies observed in $B \to K^* \mu^+\mu^-$, $R(K) = B \to K \mu^+\mu^-/B \to K e^+e^-$, and $h \to \mu\tau$. In its minimal version with two scalar doublets, the resulting model naturally explains the deviations from the Standard Model observed in $B \to K^* \mu^+\mu^-$ and $R(K)$. The CMS access in $h \to \mu\tau$ can be explained by introducing a third scalar doublet, which gives rise to a prediction for $\tau \to 3\mu$. We investigate constraints from flavour observables and direct LHC searches for $pp \to Z' \to \mu^+\mu^-$. Our model successfully generates the measured fermion-mixing matrices and does not require vector-like fermions, unlike previous attempts to explain these anomalies.


title: Light stops, blind spots, and isospin violation in the MSSM

abstract: In the framework of the MSSM, we examine several simplified models where only a few superpartners are light. This allows us to study WIMP--nucleus scattering in terms of a handful of MSSM parameters and thereby scrutinize their impact on dark matter direct-detection experiments. Focusing on spin-independent WIMP--nucleon scattering, we derive simplified, analytic expressions for the Wilson coefficients associated with Higgs and squark exchange. We utilize these results to study the complementarity of constraints due to direct-detection, flavor, and collider experiments. We also identify parameter configurations that produce (almost) vanishing cross sections. In the proximity of these so-called blind spots, we find that the amount of isospin violation may be much larger than typically expected in the MSSM. This feature is a generic property of parameter regions where cross sections are suppressed, and highlights the importance of a careful analysis of the nucleon matrix elements and the associated hadronic uncertainties. This becomes especially relevant once the increased sensitivity of future direct-detection experiments corners the MSSM into these regions of parameter space.


title: A Dynamical Potential-Density Pair for Star Clusters With Nearly Isothermal Interiors

abstract: We present a potential-density pair designed to model nearly isothermal star clusters (and similar self-gravitating systems) with a central core and an outer turnover radius, beyond which density falls off as $r^{-4}$. In the intermediate zone, the profile is similar to that of an isothermal sphere (density $\rho \propto r^{-2}$), somewhat less steep than the King 62 profile, and with the advantage that many dynamical quantities can be written in a simple closed form. We derive new analytic expressions for the cluster binding energy and velocity dispersion, and apply these to create toy models for cluster core collapse and evaporation. We fit our projected surface brightness profiles to observed globular and open clusters, and find that the quality of the fit is generally at least as good as that for the surface brightness profiles of King 62. This model can be used for convenient computation of the dynamics and evolution of globular and nuclear star clusters.


title: Probing the Small Scale Matter Power Spectrum through Dark Matter Annihilation in the Early Universe

abstract: Recent observations of the cosmic microwave background (CMB) anisotropies and the distribution of galaxies, galaxy clusters, and the Lyman Alpha forest have constrained the shape of the power spectrum of matter fluctuations on large scales k < few h/Mpc. We explore a new technique to constrain the matter power spectrum on smaller scales, assuming the dark matter is a Weakly Interacting Massive Particle (WIMP) that annihilates at early epochs. Energy released by dark matter annihilation can modify the spectrum of CMB temperature fluctuations and thus CMB experiments such as Planck have been able to constrain the quantity f <sigma v> /m < 1/88 picobarn c / GeV, where f is the fraction of energy absorbed by gas, <sigma v> is the annihilation rate assumed constant, and m is the particle mass. We assume the standard scale-invariant primordial matter power spectrum of P_prim(k) ~ k^{n_s} at large scales k < k_p, while we adopt the modified power law of P_prim(k) ~ k_p^{n_s} (k/k_p)^{m_s} at small scales. We then aim at deriving constraints on m_s. For m_s > n_s, the excess small-scale power results in a much larger number of nonlinear small mass halos, particularly at high redshifts. Dark matter annihilation in these halos releases sufficient energy to partially ionize the gas, and consequently modify the spectrum of CMB fluctuations. We show that the recent Planck data can already be used to constrain the power spectrum on small scales. For a simple model with an NFW profile with halo concentration parameter c_200 = 5 and f <sigma v> / m = 1/100 picobarn c / GeV, we can limit the mass variance sigma_{max} < 100 at the 95% confidence level, corresponding to a power law index m_s < 1.43 (1.63) for k_p = 100 (1000) h/Mpc. Our results are also relevant to theories that feature a running spectral index.


title: On Holographic Insulators and Supersolids

abstract: We obtain holographic realizations for systems that have strong similarities to Mott insulators and supersolids, after examining the ground states of Einstein-Maxwell-scalar systems. The real part of the AC conductivity has a hard gap and a discrete spectrum only. We add momentum dissipation to resolve the delta function in the conductivity due to translational invariance. We develop tools to directly calculate the Drude weight for a large class of solutions and to support our claims. Numerical RG flows are also constructed to verify that such saddle points are IR fixed points of asymptotically AdS_4 geometries.


title: Polarization alignments of radio quasars in JVAS/CLASS surveys

abstract: We test the hypothesis that the polarization vectors of flat-spectrum radio sources (FSRS) in the JVAS/CLASS 8.4-GHz surveys are randomly oriented on the sky. The sample with robust polarization measurements is made of $4155$ objects and redshift information is known for $1531$ of them. We performed two statistical analyses: one in two dimensions and the other in three dimensions when distance is available. We find significant large-scale alignments of polarization vectors for samples containing only quasars (QSO) among the varieties of FSRS's. While these correlations prove difficult to explain either by a physical effect or by biases in the dataset, the fact that the QSO's which have significantly aligned polarization vectors are found in regions of the sky where optical polarization alignments were previously found is striking.


title: The low-frequency radio catalog of flat spectrum sources

abstract: A well known property of the gamma-ray sources detected by COS-B in the 1970s, by the Compton Gamma-ray Observatory in the 1990s and recently by the Fermi observations is the presence of radio counterparts, in particular for those associated to extragalactic objects. This observational evidence is the basis of the radio-gamma-ray connection established for the class of active galactic nuclei known as blazars. In particular, the main spectral property of the radio counterparts associated with gamma-ray blazars is that they show a flat spectrum in the GHz frequency range. Our recent analysis dedicated to search blazar-like candidates as potential counterparts for the unidentified gamma-ray sources (UGSs) allowed us to extend the radio-gamma-ray connection in the MHz regime. We also showed that below 1 GHz blazars maintain flat radio spectra. Thus on the basis of these new results, we assembled a low-frequency radio catalog of flat spectrum sources built by combining the radio observations of the Westerbork Northern Sky Survey (WENSS) and of the Westerbork in the southern hemisphere (WISH) catalog with those of the NRAO Very Large Array Sky survey (NVSS). This could be used in the future to search for new, unknown blazar-like counterparts of the gamma-ray sources. First we found NVSS counterparts of WSRT radio sources and then we selected flat spectrum radio sources according to a new spectral criterion specifically defined for radio observations performed below 1 GHz. We also described the main properties of the catalog listing 28358 radio sources and their logN-logS distributions. Finally a comparison with with the Green Bank 6-cm radio source catalog has been performed to investigate the spectral shape of the low-frequency flat spectrum radio sources at higher frequencies.


title: Anisotropic Propagator for the Goldstone Modes in Color-flavor Locked Phase in the Presence of a Magnetic Field

abstract: We consider the phase diagram of QCD at very high baryon density and at zero temperature in the presence of a strong magnetic field. The state of matter at such high densities and low temperatures is believed to be a phase known as the color-flavor locked phase which breaks color and electromagnetic gauge invariance leaving a linear combination of them unbroken. Of the 9 quarks (three flavors and three colors), five are neutral under this unbroken generator and four are oppositely charged. In the presence of a magnetic field corresponding to the unbroken generator however, the properties of the condensate changes and a new phase known as the magnetic color flavor locked (MCFL)phase is realized. This phase breaks some of the color-flavor symmetry of the Lagrangian spontaneously, giving rise to 6 Goldstone modes, 5 of which are pseudo Goldstone modes. These Goldstone modes are composed of excitations that correspond to both neutral quarks and charged quarks. Hence it is natural to expect that the propagators of these Goldstone modes get affected in the presence of a magnetic field and their speed becomes considerably anisotropic. Although this anisotropy is self-evident from symmetry arguments, it has not been quantified yet. We calculate this anisotropy in the speed of the Goldstone modes using an NJL model type of interaction between the quarks and comment on the impact of such anisotropic modes on the transport properties of the MCFL phase.


title: Self-induced flavor instabilities of a dense neutrino stream in a two-dimensional model

abstract: We consider a simplifed model for self-induced flavor conversions of a dense neutrino gas in two dimensions, showing new solutions that spontaneously break the spatial symmetries of the initial conditions. As a result of the symmetry breaking induced by the neutrino-neutrino interactions, the coherent behavior of the neutrino gas becomes unstable. This instability produces large spatial variations in the flavor content of the ensemble. Furthermore, it also leads to the creation of domains of different net lepton number flux. The transition of the neutrino gas from a coherent to incoherent behavior shows an intriguing analogy with a streaming flow changing from laminar to turbulent regime. These finding would be relevant for the self-induced conversions of neutrinos streaming-off a supernova core.


title: Gromov compactness for squiggly strip shrinking in pseudoholomorphic quilts

abstract: We establish a Gromov compactness theorem for strip shrinking in pseudoholomorphic quilts when composition of Lagrangian correspondences is immersed. In particular, we show that figure eight bubbling occurs in the limit, argue that this is a codimension-$0$ effect, and predict its algebraic consequences -- geometric composition extends to a curved $A_\infty$-bifunctor, in particular the associated Floer complexes are isomorphic after a figure eight correction of the bounding cochain. An appendix with Felix Schm\"{a}schke provides examples of nontrivial figure eight bubbles.


title: The angle-averaged squeezed limit of nonlinear matter N-point functions

abstract: We show that in a certain, angle-averaged squeezed limit, the $N$-point function of matter is related to the response of the matter power spectrum to a long-wavelength density perturbation, $P^{-1}d^nP(k|\delta_L)/d\delta_L^n|_{\delta_L=0}$, with $n=N-2$. By performing N-body simulations with a homogeneous overdensity superimposed on a flat Friedmann-Robertson-Lema\^itre-Walker (FRLW) universe using the \emph{separate universe} approach, we obtain measurements of the nonlinear matter power spectrum response up to $n=3$, which is equivalent to measuring the fully nonlinear matter $3-$ to $5-$point function in this squeezed limit. The sub-percent to few percent accuracy of those measurements is unprecedented. We then test the hypothesis that nonlinear $N$-point functions at a given time are a function of the linear power spectrum at that time, which is predicted by standard perturbation theory (SPT) and its variants that are based on the ideal pressureless fluid equations. Specifically, we compare the responses computed from the separate universe simulations and simulations with a rescaled initial (linear) power spectrum amplitude. We find discrepancies of 10\% at $k\simeq 0.2 - 0.5 \,h\,{\rm Mpc}^{-1}$ for $5-$ to $3-$point functions at $z=0$. The discrepancy occurs at higher wavenumbers at $z=2$. Thus, SPT and its variants, carried out to arbitrarily high order, are guaranteed to fail to describe matter $N$-point functions ($N>2$) around that scale.


title: Estimating the Mean Number of K-Means Clusters to Form

abstract: Utilizing the sample size of a dataset, the random cluster model is employed in order to derive an estimate of the mean number of K-Means clusters to form during classification of a dataset.


title: The SWIFT AGN and Cluster Survey I: Number Counts of AGN and Galaxy Clusters

abstract: The Swift AGN and Cluster Survey (SACS) uses 125 deg^2 of Swift XRT serendipitous fields with variable depths surrounding gamma-ray bursts to provide a medium depth (4e-15 erg/s/cm^2) and area survey filling the gap between deep, narrow Chandra/XMM-Newton surveys and wide, shallow ROSAT surveys. Here we present a catalog of 22,563 point sources and 442 extended sources and examine the number counts of the AGN and galaxy cluster populations. SACS provides excellent constraints on the AGN number counts at the bright end with negligible uncertainties due to cosmic variance, and these constraints are consistent with previous measurements. We use Wise mid-infrared (MIR) colors to classify the sources. For AGN we can roughly separate the point sources into MIR-red and MIR-blue AGN, finding roughly equal numbers of each type in the soft X-ray band (0.5-2 keV), but fewer MIR-blue sources in the hard X-ray band (2-8 keV). The cluster number counts, with 5% uncertainties from cosmic variance, are also consistent with previous surveys but span a much larger continuous flux range. Deep optical or IR follow-up observations of this cluster sample will significantly increase the number of higher redshift (z > 0.5) X-ray-selected clusters.


title: On robust solutions to uncertain linear complementarity problems and their variants

abstract: A popular approach for addressing uncertainty in variational inequality problems is by solving the expected residual minimization (ERM) problem. This avenue necessitates distributional information associated with the uncertainty and requires minimizing nonconvex expectation-valued functions. We consider a distinctly different approach in the context of uncertain linear complementarity problems with a view towards obtaining robust solutions. Specifically, we define a robust solution to a complementarity problem as one that minimizes the worst-case of the gap function. In what we believe is amongst the first efforts to comprehensively address such problems in a distribution-free environment, we show that under specified assumptions on the uncertainty sets, the robust solutions to uncertain monotone linear complementarity problem can be tractably obtained through the solution of a single convex program. We also define uncertainty sets that ensure that robust solutions to non-monotone generalizations can also be obtained by solving convex programs. More generally, robust counterparts of uncertain non-monotone LCPs are proven to be low-dimensional nonconvex quadratically constrained quadratic programs. We show that these problems may be globally resolved by customizing an existing branching scheme. We further extend the tractability results to include uncertain affine variational inequality problems defined over uncertain polyhedral sets as well as to hierarchical regimes captured by mathematical programs with uncertain complementarity constraints. Preliminary numerics on uncertain linear complementarity and traffic equilibrium problems suggest that the presented avenues hold promise.


title: Properties of simple sets in digital spaces. Contractions of simple sets preserving the homotopy type of a digital space

abstract: A point of a digital space is called simple if it can be deleted from the space without altering topology. This paper introduces the notion simple set of points of a digital space. The definition is based on contractible spaces and contractible transformations. A set of points in a digital space is called simple if it can be contracted to a point without changing topology of the space. It is shown that contracting a simple set of points does not change the homotopy type of a digital space, and the number of points in a digital space without simple points can be reduces by contracting simple sets. Using the process of contracting, we can substantially compress a digital space while preserving the topology. The paper proposes a method for thinning a digital space which shows that this approach can contribute to computer science such as medical imaging, computer graphics and pattern analysis.


title: Parallel Statistical Multi-resolution Estimation

abstract: We discuss several strategies to implement Dykstra's projection algorithm on NVIDIA's compute unified device architecture (CUDA). Dykstra's algorithm is the central step in and the computationally most expensive part of statistical multi-resolution methods. It projects a given vector onto the intersection of convex sets. Compared with a CPU implementation our CUDA implementation is one order of magnitude faster. For a further speed up and to reduce memory consumption we have developed a new variant, which we call incomplete Dykstra's algorithm. Implemented in CUDA it is one order of magnitude faster than the CUDA implementation of the standard Dykstra algorithm. As sample application we discuss using the incomplete Dykstra's algorithm as preprocessor for the recently developed super-resolution optical fluctuation imaging (SOFI) method (Dertinger et al. 2009). We show that statistical multi-resolution estimation can enhance the resolution improvement of the plain SOFI algorithm just as the Fourier-reweighting of SOFI. The results are compared in terms of their power spectrum and their Fourier ring correlation (Saxton and Baumeister 1982). The Fourier ring correlation indicates that the resolution for typical second order SOFI images can be improved by about 30 per cent. Our results show that a careful parallelization of Dykstra's algorithm enables its use in large-scale statistical multi-resolution analyses.


title: Aberration Corrected Emittance Exchange

abstract: Full exploitation of emittance exchange (EEX) requires aberration-free performance of a complex imaging system including active radio-frequency (RF) elements which can add temporal distortions. We investigate the performance of an EEX line where the exchange occurs between two dimensions with normalized emittances which differ by multiple orders of magnitude. The transverse emittance is exchanged into the longitudinal dimension using a double dog-leg emittance exchange setup with a five cell RF deflector cavity. Aberration correction is performed on the four most dominant aberrations. These include temporal aberrations that are corrected with higher order magnetic optical elements located where longitudinal and transverse emittance are coupled. We demonstrate aberration-free performance of an EEX line with emittances differing by four orders of magnitude, \textit{i.e.} an initial transverse emittance of 1~pm-rad is exchanged with a longitudinal emittance of 10~nm-rad.


title: Millennial-scale stable oscillations between sea ice and convective deep water formation

abstract: During the last ice age there were several quasi-periodic abrupt warming events. The climatic effects of the so-called Dansgaard-Oeschger (DO) events were felt globally, although the North Atlantic experienced the largest and most abrupt temperature anomalies. Similar but weaker oscillations also took place during the interglacial period. This paper proposes an auto-oscillatory mechanism between sea ice and convective deep water formation in the north Atlantic as the source of the persistent cycles. A simple dynamical model is constructed by coupling and slightly modifying two existing models of ocean circulation and sea ice. The model exhibits mixed mode oscillations, consisting of decadal scale small amplitude oscillations, and a large amplitude relaxation fluctuation. The decadal oscillations occur due to the insulating effect of sea ice and leads to periodic ventilation of heat from the polar ocean. Gradually an instability builds up in the polar column and results in an abrupt initiation of convection and polar warming. The unstable convective state relaxes back to the small amplitude oscillations from where the process repeats in a self-sustained manner. Freshwater pulses mimicking Heinrich events cause the oscillations to be grouped into packets of progressively weaker fluctuations, as also observed in the proxy records. Modulation of this stable oscillation mechanism by freshwater and insolation variations could account for the distribution and pacing of DO and Bond events. Physical aspects of the system such as sea ice extent and oceanic advective flow rates could determine the characteristic 1,500 year timescale of DO events.


title: Realistic estimate of valence transversity distributions from inclusive dihadron production

abstract: We present an updated extraction of the transversity parton distribution based on the analysis of pion-pair production in deep-inelastic scattering off transversely polarized targets in collinear factorization. Data for proton and deuteron targets make it possible to perform a flavor separation of the valence components of the transversity distribution, using di-hadron fragmentation functions taken from the semi-inclusive production of two pion pairs in back-to-back jets in e+e- annihilation. The e+e- data from Belle have been reanalyzed using the replica method and a more realistic estimate of the uncertainties on the chiral-odd interference fragmentation function has been obtained. Then, the transversity distribution has been extracted by using the most recent and more precise COMPASS data for deep-inelastic scattering off proton targets. Our results represent the most accurate estimate of the uncertainties on the valence components of the transversity distribution currently available.


title: An introduction to the New SI

abstract: Plans are underway to redefine the International System of Units (SI) around 2018. The New SI specifies the values of certain physical constants to define units. This article explains the New SI in a way that could be used to present it to high-school physics classes.


title: Pseudo prolate spheroidal functions

abstract: Let $D_{T}$ and $B_{\Omega }$ denote the operators which cut the time content outside $T$ and the frequency content outside $\Omega $, respectively. The prolate spheroidal functions are the eigenfunctions of the operator $P_{T,\Omega }=D_{T}B_{\Omega }D_{T}$. With the aim of formulating in precise mathematical terms the notion of Nyquist rate, Landau and Pollack have shown that, asymptotically, the number of such functions with eigenvalue close to one is $\approx \frac{\left\vert T\right\vert \left\vert \Omega \right\vert }{2\pi }$. We have recently revisited this problem with a new approach: instead of counting the number of eigenfunctions with eigenvalue close to one, we count the maximum number of orthogonal $\epsilon$-pseudoeigenfunctions with $\epsilon $-pseudoeigenvalue one. Precisely, we count how many orthogonal functions have a maximum of energy $\epsilon $ outside the domain $T\times \Omega $, in the sense that $\left\Vert P_{T,\Omega }f-f\right\Vert ^{2}\leq \epsilon $. We have recently discovered that the sharp asymptotic number is $\approx (1-\epsilon )^{-1}\frac{\left\vert T\right\vert \left\vert \Omega \right\vert }{2\pi }$. The proof involves an explicit construction of the pseudoeigenfunctions of $P_{T,\Omega }$. When $T$ and $\Omega$ are intervals we call them pseudo prolate spheroidal functions. In this paper we explain how they are constructed.


title: Boron Nitride Nanosheets as Improved and Reusable Substrates for Gold Nanoparticles Enabled Surface Enhanced Raman Spectroscopy

abstract: Atomically thin boron nitride (BN) nanosheets have been found an excellent substrate for noble metal particles enabled surface enhanced Raman spectroscopy (SERS), thanks to their good adsorption of aromatic molecules, high thermal stability and weak Raman scattering. Faceted gold (Au) nanoparticles have been synthesized on BN nanosheets by a simple but controllable and reproducible sputtering and annealing method. The size and density of the Au particles can be controlled by sputtering time, current and annealing temperature etc. Under the same sputtering and annealing conditions, the Au particles on BN of different thicknesses show various sizes because the surface diffusion coefficients of Au depends on the thickness of BN. Intriguingly, decorated with similar morphology and distribution of Au particles, BN nanosheets exhibit better Raman enhancements than silicon substrate as well as bulk BN crystals. Additionally, BN nanosheets show no noticeable SERS signal and hence cause no interference to the Raman signal of analyte. The Au/BN substrates can be reused by heating in air to remove adsorbed analyte without loss of SERS enhancement.


title: A Low-Resolution Spectroscopic Exploration of Puzzling OGLE Variable Stars

abstract: We present the results of a spectroscopic follow-up of various puzzling variable objects detected in the OGLE-III Galactic disk and bulge fields. The sample includes mainly short-period multi-mode pulsating stars that could not have been unambiguously classified as either delta Sct or beta Cep type stars based on photometric data only, also stars with irregular fluctuations mimicking cataclysmic variables and stars with dusty shells, and periodic variables displaying brightenings in their light curves that last for more than half of the period. The obtained low-resolution spectra show that all observed short-period pulsators are of delta Sct type, the stars with irregular fluctuations are young stellar objects, and the objects with regular brightenings are A type stars or very likely Ap stars with strong magnetic field responsible for the presence of bright caps around magnetic poles on their surface. We also took spectra of objects designated OGLE-GD-DSCT-0058 and OGLE-GD-CEP-0013. An estimated effective temperature of 33,000 K in OGLE-GD-DSCT-0058 indicates that it cannot be a delta Sct type variable. This very short-period (0.01962 d) high-amplitude (0.24 mag in the I-band) object remains a mystery. It may represent a new class of variable stars. The spectrum of OGLE-GD-CEP-0013 confirms that this is a classical Cepheid despite a peculiar shape of its light curve. The presented results will help in proper classification of variable objects in the OGLE Galactic Variability Survey.


title: Giant Triplet Proximity Effect in $\pi$-biased Josephson Junctions with Spin-Orbit Coupling

abstract: In diffusive Josephson junctions the phase-difference $\phi$ between the superconductors strongly influences the spectroscopic features of the layer separating them. The observation of a uniform minigap and its phase modulation were only recently experimentally reported, demonstrating agreement with theoretical predictions up to now - a vanishing minigap at $\phi=\pi$. Remarkably, we find that in the presence of intrinsic spin-orbit coupling a giant proximity effect due to spin-triplet Cooper pairs can develop at $\phi=\pi$, in complete contrast to the suppressed proximity effect without spin-orbit coupling. We here report a combined numerical and analytical study of this effect, proving its presence solely based on symmetry arguments, which makes it independent of the specific parameters used in experiments. We show that the spectroscopic signature of the triplets is present throughout the entire ferromagnetic layer. Our finding offers a new way to artificially create, control and isolate spin-triplet superconductivity.


title: Connecting Angular Momentum and Galactic Dynamics: The complex Interplay between Spin, Mass, and Morphology

abstract: The evolution and distribution of the angular momentum of dark matter (DM) halos have been discussed in several studies over the past decades. In particular, the idea arose that angular momentum conservation should allow to infer the total angular momentum of the entire DM halo from measuring the angular momentum of the baryonic component, which is populating the center of the halo, especially for disk galaxies. To test this idea and to understand the connection between the angular momentum of the DM halo and its galaxy, we use the Magneticum simulations. We successfully produce populations of spheroidal and disk galaxies self-consistently. Thus, we are able to study the dependence of galactic properties on their morphology. We find that (1) the specific angular momentum of stars in disk and spheroidal galaxies as a function of their stellar mass compares well with observational results; (2) the specific angular momentum of the stars in disk galaxies is slightly smaller compared to the specific angular momentum of the cold gas, in good agreement with observations; (3) simulations including the baryonic component show a dichotomy in the specific stellar angular momentum distribution when splitting the galaxies according to their morphological type (this dichotomy can also be seen in the spin parameter, where disk galaxies populate halos with slightly larger spin compared to spheroidal galaxies); (4) disk galaxies preferentially populate halos in which the angular momentum vector of the DM component in the central part shows a better alignment to the angular momentum vector of the entire halo; and (5) the specific angular momentum of the cold gas in disk galaxies is approximately 40 percent smaller than the specific angular momentum of the total DM halo and shows a significant scatter.


title: Axion Like Particles and the Inverse Seesaw Mechanism

abstract: Light pseudoscalars known as axion like particles (ALPs) may be behind physical phenomena like the Universe transparency to ultra-energetic photons, the soft $\gamma$-ray excess from the Coma cluster, and the 3.5 keV line. We explore the connection of these particles with the inverse seesaw (ISS) mechanism for neutrino mass generation. We propose a very restrictive setting where the scalar field hosting the ALP is also responsible for generating the ISS mass scales through its vacuum expectation value on gravity induced nonrenormalizable operators. A discrete gauge symmetry protects the theory from the appearance of overly strong gravitational effects and discrete anomaly cancellation imposes strong constraints on the order of the group. The anomalous U$(1)$ symmetry leading to the ALP is an extended lepton number and the protective discrete symmetry can be always chosen as a subgroup of a combination of the lepton number and the baryon number.


title: Structural properties of artificial halos in non-standard dark matter simulations

abstract: Artificial fragmentation of the matter density field causes the formation of spurious groups of particles in N-body simulations of non-standard Dark Matter (DM) models which are characterized by a small scale cut-off in the linear matter power spectrum. These spurious halos alter the prediction of the mass function in a range of masses where differences among DM models are most relevant to observational tests. Using a suite of high resolution simulations we show that the contamination of artificial groups of particles significantly affect the statistics of halo spin, shape and virial state parameters. We find that spurious halos have systematically larger spin values, are highly elliptical or prolate and significantly deviate from virial equilibrium. These characteristics allow us to detect the presence of spurious halos even in non-standard DM models for which the low-mass end of the mass function remains well behaved. We show that selecting halos near the virial equilibrium provides a simple and effective method to remove the bulk of spurious halos from numerical halo catalogs and consistently recover the halo mass function at low masses.


title: Extracting Periodic Transit Signals from Noisy Light Curves using Fourier Series

abstract: We present a simple and powerful method for extracting transit signals associated with a known transiting planet from noisy light curves. Assuming the orbital period of the planet is known and the signal is periodic, we illustrate that systematic noise can be removed in Fourier space at all frequencies, by only using data within a fixed time frame with a width equal to an integer number of orbital periods. This results in a reconstruction of the full transit signal which on average is unbiased, despite that no prior knowledge of either the noise or the transit signal itself is used in the analysis. The method has therefore clear advantages over standard phase folding, which normally requires external input such as nearby stars or noise models for removing systematic components. In addition, we can extract the full orbital transit signal ($360$ degrees) simultaneously, and \emph{Kepler} like data can be analyzed in just a few seconds. We illustrate the performance of our method by applying it to a dataset composed of light curves from \emph{Kepler} with a fake injected signal emulating a planet with rings. For extracting periodic transit signals, our presented method is in general the optimal and least biased estimator and could therefore lead the way toward the first detections of, e.g., planet rings and exo-trojan asteroids.


title: Measurement of very forward neutron energy spectra for 7 TeV proton--proton collisions at the Large Hadron Collider

abstract: The Large Hadron Collider forward (LHCf) experiment is designed to use the LHC to verify the hadronic-interaction models used in cosmic-ray physics. Forward baryon production is one of the crucial points to understand the development of cosmic-ray showers. We report the neutron-energy spectra for LHC $\sqrt{s}$ = 7 TeV proton--proton collisions with the pseudo-rapidity $\eta$ ranging from 8.81 to 8.99, from 8.99 to 9.22, and from 10.76 to infinity. The measured energy spectra obtained from the two independent calorimeters of Arm1 and Arm2 show the same characteristic feature before unfolding the difference in the detector responses. We unfolded the measured spectra by using the multidimensional unfolding method based on Bayesian theory, and the unfolded spectra were compared with current hadronic-interaction models. The QGSJET II-03 model predicts a high neutron production rate at the highest pseudo-rapidity range similar to our results and the DPMJET 3.04 model describes our results well at the lower pseudo-rapidity ranges. However no model perfectly explains the experimental results in the whole pseudo-rapidity range. The experimental data indicate the most abundant neutron production rate relative to the photon production, which does not agree with predictions of the models.


title: Diverse Landmark Sampling from Determinantal Point Processes for Scalable Manifold Learning

abstract: High computational costs of manifold learning prohibit its application for large point sets. A common strategy to overcome this problem is to perform dimensionality reduction on selected landmarks and to successively embed the entire dataset with the Nystr\"om method. The two main challenges that arise are: (i) the landmarks selected in non-Euclidean geometries must result in a low reconstruction error, (ii) the graph constructed from sparsely sampled landmarks must approximate the manifold well. We propose the sampling of landmarks from determinantal distributions on non-Euclidean spaces. Since current determinantal sampling algorithms have the same complexity as those for manifold learning, we present an efficient approximation running in linear time. Further, we recover the local geometry after the sparsification by assigning each landmark a local covariance matrix, estimated from the original point set. The resulting neighborhood selection based on the Bhattacharyya distance improves the embedding of sparsely sampled manifolds. Our experiments show a significant performance improvement compared to state-of-the-art landmark selection techniques.


title: Hypersurfaces with constant principal curvatures in $\mathbb{S}^{n}\times\mathbb{R}$ and $\mathbb{H}^{n}\times\mathbb{R}$

abstract: In this paper, we classify the hypersurfaces in $\mathbb{S}^{n}\times \mathbb{R}$ and $\mathbb{H}^{n}\times\mathbb{R}$, $n\neq 3$, with $g$ distinct constant principal curvatures, $g\in\{1,2,3\}$, where $\mathbb{S}^{n}$ and $\mathbb{H}^{n}$ denote the sphere and hyperbolic space of dimension $n$, respectively. We prove that such hypersurfaces are isoparametric in those spaces. Furthermore, we find a necessary and sufficient condition for an isoparametric hypersurface in $\mathbb{S}^{n}\times \mathbb{R}\subset \mathbb{R}^{n+2}$ and $\mathbb{H}^{n}\times\mathbb{R}\subset \mathbb{L}^{n+2}$ with flat normal bundle, having constant principal curvatures.


title: Fall-off of eigenfunctions for non-local Schr\"odinger operators with decaying potentials

abstract: We study the spatial decay of eigenfunctions of non-local Schr\"odinger operators whose kinetic terms are generators of symmetric jump-paring L\'evy processes with Kato-class potentials decaying at infinity. This class of processes has the property that the intensity of single large jumps dominates the intensity of all multiple large jumps. We find that the decay rates of eigenfunctions depend on the process via specific preference rates in particular jump scenarios, and depend on the potential through the distance of the corresponding eigenvalue from the edge of the continuous spectrum. We prove that the conditions of the jump-paring class imply that for all eigenvalues the corresponding positive eigenfunctions decay at most as rapidly as the L\'evy intensity. This condition is sharp in the sense that if the jump-paring property fails to hold, then eigenfunction decay becomes slower than the decay of the L\'evy intensity. We furthermore prove that under reasonable conditions the L\'evy intensity also governs the upper bounds of eigenfunctions, and ground states are comparable with it, i.e., two-sided bounds hold. As an interesting consequence, we identify a sharp regime change in the decay of eigenfunctions as the L\'evy intensity is varied from sub-exponential to exponential order, and dependent on the location of the eigenvalue, in the sense that through the transition L\'evy intensity-driven decay becomes slower than the rate of decay of L\'evy intensity. Our approach is based on path integration and probabilistic potential theory techniques, and all results are also illustrated by specific examples.


title: Safe Leads and Lead Changes in Competitive Team Sports

abstract: We investigate the time evolution of lead changes within individual games of competitive team sports. Exploiting ideas from the theory of random walks, the number of lead changes within a single game follows a Gaussian distribution. We show that the probability that the last lead change and the time of the largest lead size are governed by the same arcsine law, a bimodal distribution that diverges at the start and at the end of the game. We also determine the probability that a given lead is "safe" as a function of its size $L$ and game time $t$. Our predictions generally agree with comprehensive data on more than 1.25 million scoring events in roughly 40,000 games across four professional or semi-professional team sports, and are more accurate than popular heuristics currently used in sports analytics.


title: Charge transport through a semiconductor quantum dot-ring nanostructure

abstract: Transport properties of a gated nanostructure depend crucially on the coupling of its states to the states of electrodes. In the case of a single quantum dot the coupling, for a given quantum state, is constant or can be slightly modified by additional gating. In this paper we consider a concentric dot-ring nanostructure (DRN) and show that its transport properties can be drastically modified due to the unique geometry. We calculate the dc current through a DRN in the Coulomb blockade regime and show that it can efficiently work as a single electron transistor or a current rectifier. In both cases the transport characteristics strongly depends on the details of the confinement potential. The calculations are carried out for low and high bias regime, the latter being especially interesting in the context of current rectification due to fast relaxation processes.


title: Cross platform app: a comparative study

abstract: The use of mobile applications is now so common that users now expect companies whose services which they consume already have an application to provide these services or a mobile version of your site, but this is not always simple to do or cheap. Thus, the hybrid development has emerged as a potential alternative to this need. The evolution of this new paradigm has taken the attention of researchers and companies as viable alternative to the mobile development. This paper shows how hybrid development can be an alternative for companies provide their services with a low investment and still offer a great service to their clients.


title: Is language evolution grinding to a halt? The scaling of lexical turbulence in English fiction suggests it is not

abstract: Of basic interest is the quantification of the long term growth of a language's lexicon as it develops to more completely cover both a culture's communication requirements and knowledge space. Here, we explore the usage dynamics of words in the English language as reflected by the Google Books 2012 English Fiction corpus. We critique an earlier method that found decreasing birth and increasing death rates of words over the second half of the 20th Century, showing death rates to be strongly affected by the imposed time cutoff of the arbitrary present and not increasing dramatically. We provide a robust, principled approach to examining lexical evolution by tracking the volume of word flux across various relative frequency thresholds. We show that while the overall statistical structure of the English language remains stable over time in terms of its raw Zipf distribution, we find evidence of an enduring `lexical turbulence': The flux of words across frequency thresholds from decade to decade scales superlinearly with word rank and exhibits a scaling break we connect to that of Zipf's law. To better understand the changing lexicon, we examine the contributions to the Jensen-Shannon divergence of individual words crossing frequency thresholds. We also find indications that scholarly works about fiction are strongly represented in the 2012 English Fiction corpus, and suggest that a future revision of the corpus should attempt to separate critical works from fiction itself.

categories: cs.CL cs.IT math.IT physics.soc-ph stat.AP


title: The Effects of QCD Equation of State on the Relic Density of WIMP Dark Matter

abstract: Weakly Interactive Massive Particles (WIMPs) are the most widely studied candidate particles forming the cold dark matter (CDM) whose existence can be inferred from a wealth of astrophysical and cosmological observations. In the framework of the minimal cosmological model detailed measurements on the cosmic microwave background by the PLANCK collaboration fix the scaled CDM relic density to $\Omega_{c}h^2=0.1193\pm0.0014$, with an error of less than 1.5%. In order to fully exploit this observational precision, theoretical calculations should have a comparable or smaller error. In this paper we use recent lattice QCD calculations to improve the description of the thermal plasma. This affects the predicted relic density of "thermal WIMPs", which once were in chemical equilibrium with Standard Model particles. For WIMP masses between 3 and 15 GeV, where QCD effects are most important, our predictions differ from earlier results by up to 9% (12%) for pure S-wave (P-wave) annihilation. We use these results to compute the thermally averaged WIMP annihilation cross section that reproduces the correct CDM relic density, for WIMP masses between 0.1 GeV and 10 TeV.


title: Appearance-based indoor localization: A comparison of patch descriptor performance

abstract: Vision is one of the most important of the senses, and humans use it extensively during navigation. We evaluated different types of image and video frame descriptors that could be used to determine distinctive visual landmarks for localizing a person based on what is seen by a camera that they carry. To do this, we created a database containing over 3 km of video-sequences with ground-truth in the form of distance travelled along different corridors. Using this database, the accuracy of localization - both in terms of knowing which route a user is on - and in terms of position along a certain route, can be evaluated. For each type of descriptor, we also tested different techniques to encode visual structure and to search between journeys to estimate a user's position. The techniques include single-frame descriptors, those using sequences of frames, and both colour and achromatic descriptors. We found that single-frame indexing worked better within this particular dataset. This might be because the motion of the person holding the camera makes the video too dependent on individual steps and motions of one particular journey. Our results suggest that appearance-based information could be an additional source of navigational data indoors, augmenting that provided by, say, radio signal strength indicators (RSSIs). Such visual information could be collected by crowdsourcing low-resolution video feeds, allowing journeys made by different users to be associated with each other, and location to be inferred without requiring explicit mapping. This offers a complementary approach to methods based on simultaneous localization and mapping (SLAM) algorithms.


title: Bi-cross-validation for factor analysis

abstract: Factor analysis is over a century old, but it is still problematic to choose the number of factors for a given data set. The scree test is popular but subjective. The best performing objective methods are recommended on the basis of simulations. We introduce a method based on bi-cross-validation, using randomly held-out submatrices of the data to choose the number of factors. We find it performs better than the leading methods of parallel analysis (PA) and Kaiser's rule. Our performance criterion is based on recovery of the underlying factor-loading (signal) matrix rather than identifying the true number of factors. Like previous comparisons, our work is simulation based. Recent advances in random matrix theory provide principled choices for the number of factors when the noise is homoscedastic, but not for the heteroscedastic case. The simulations we choose are designed using guidance from random matrix theory. In particular, we include factors too small to detect, factors large enough to detect but not large enough to improve the estimate, and two classes of factors large enough to be useful. Much of the advantage of bi-cross-validation comes from cases with factors large enough to detect but too small to be well estimated. We also find that a form of early stopping regularization improves the recovery of the signal matrix.


title: Understanding The Effects Of Stellar Multiplicity On The Derived Planet Radii From Transit Surveys: Implications for Kepler, K2, and TESS

abstract: We present a study on the effect of undetected stellar companions on the derived planetary radii for the Kepler Objects of Interest (KOIs). The current production of the KOI list assumes that the each KOI is a single star. Not accounting for stellar multiplicity statistically biases the planets towards smaller radii. The bias towards smaller radii depends on the properties of the companion stars and whether the planets orbit the primary or the companion stars. Defining a planetary radius correction factor $X_R$, we find that if the KOIs are assumed to be single, then, {\it on average}, the planetary radii may be underestimated by a factor of $\langle X_R \rangle \approx 1.5$. If typical radial velocity and high resolution imaging observations are performed and no companions are detected, this factor reduces to $\langle X_R \rangle \approx 1.2$. The correction factor $\langle X_R \rangle$ is dependent upon the primary star properties and ranges from $\langle X_R \rangle \approx 1.6$ for A and F stars to $\langle X_R \rangle \approx 1.2$ for K and M stars. For missions like K2 and TESS where the stars may be closer than the stars in the Kepler target sample, observational vetting (primary imaging) reduces the radius correction factor to $\langle X_R \rangle \approx 1.1$. Finally, we show that if the stellar multiplicity rates are not accounted for correctly, occurrence rate calculations for Earth-sized planets may overestimate the frequency of small planets by as much as $15-20$\%.


title: Switching to Learn

abstract: A network of agents attempt to learn some unknown state of the world drawn by nature from a finite set. Agents observe private signals conditioned on the true state, and form beliefs about the unknown state accordingly. Each agent may face an identification problem in the sense that she cannot distinguish the truth in isolation. However, by communicating with each other, agents are able to benefit from side observations to learn the truth collectively. Unlike many distributed algorithms which rely on all-time communication protocols, we propose an efficient method by switching between Bayesian and non-Bayesian regimes. In this model, agents exchange information only when their private signals are not informative enough; thence, by switching between the two regimes, agents efficiently learn the truth using only a few rounds of communications. The proposed algorithm preserves learnability while incurring a lower communication cost. We also verify our theoretical findings by simulation examples.


title: Achievements and open issues in the determination of polarized parton distribution functions

abstract: I review the current status of the determination of helicity-dependent, or polarized, parton distribution functions from a comprehensive analysis of experimental data in perturbative quantum chromodynamics. I illustrate the latest achievements driven by new measurements in polarized proton-proton collisions at the Relativistic Heavy Ion Collider, namely the first evidence of a sizable polarized light sea quark asymmetry and of a positive polarized gluon distribution in the proton. I discuss which are the open issues in the determination of polarized distributions, and how these may be addressed in the future by ongoing, planned and proposed experimental programs.