title: Investigating the Maximum Number of Real Solutions to the Power Flow Equations: Analysis of Lossless Four-Bus Systems

abstract: The power flow equations model the steady-state relationship between the power injections and voltage phasors in an electric power system. By separating the real and imaginary components of the voltage phasors, the power flow equations can be formulated as a system of quadratic polynomials. Only the real solutions to these polynomial equations are physically meaningful. This paper focuses on the maximum number of real solutions to the power flow equations. An upper bound on the number of real power flow solutions commonly used in the literature is the maximum number of complex solutions. There exist two- and three-bus systems for which all complex solutions are real. It is an open question whether this is also the case for larger systems. This paper investigates four-bus systems using techniques from numerical algebraic geometry and conjectures a negative answer to this question. In particular, this paper studies lossless, four-bus systems composed of PV buses connected by lines with arbitrary susceptances. Computing the Galois group, which is degenerate, enables conversion of the problem of counting the number of real solutions to the power flow equations into counting the number of positive roots of a univariate sextic polynomial. From this analysis, it is conjectured that the system has at most 16 real solutions, which is strictly less than the maximum number of complex solutions, namely 20. We also provide explicit parameter values where this system has 16 real solutions so that the conjectured upper bound is achievable.


title: Acylindrical hyperbolicity, non simplicity and SQ-universality of groups splitting over Z

abstract: We show, using acylindrical hyperbolicity, that a finitely generated group splitting over $\Z$ cannot be simple. We also obtain SQ-universality in most cases, for instance a balanced group (one where if two powers of an infinite order element are conjugate then they are equal or inverse) which is finitely generated and splits over $\Z$ must either be SQ-universal or it is one of exactly seven virtually abelian exceptions.


title: Reconstruction of the Higgs mass in events with Higgs bosons decaying into a pair of tau leptons using matrix element techniques

abstract: We present an algorithm for the reconstruction of the Higgs mass in events with Higgs bosons decaying into a pair of tau leptons. The algorithm is based on matrix element (ME) techniques and achieves a relative resolution on the Higgs boson mass of typically 15-20%. A previous version of the algorithm has been used in analyses of Higgs boson production performed by the CMS collaboration during LHC Run 1. The algorithm is described in detail and its performance on simulated events is assessed. The development of techniques to handle tau decays in the ME formalism represents an important result of this paper.


title: Simulation tests of galaxy cluster constraints on chameleon gravity

abstract: We use two new hydrodynamical simulations of $\Lambda$CDM and $f(R)$ gravity to test the methodology used by Wilcox et al. 2015 (W15) in constraining the effects of a fifth force on the profiles of clusters of galaxies. We construct realistic simulated stacked weak lensing and X-ray surface brightness cluster profiles from these cosmological simulations, and then use these data projected along various lines-of-sight to test the spherical symmetry of our stacking procedure. We also test the applicability of the NFW profile to model weak lensing profiles of clusters in $f(R)$ gravity. Finally, we test the validity of the analytical model developed in W15 against the simulated profiles. Overall, we find our methodology is robust and broadly agrees with these simulated data. We also apply our full Markov Chain Monte Carlo (MCMC) analysis from W15 to our simulated X-ray and lensing profiles, providing consistent constraints on the modified gravity parameters as obtained from the real cluster data, e.g. for our $\Lambda$CDM simulation we obtain $|f_{\rm{R}0}| < 8.3 \times 10^{-5}$ (95% CL), which is in good agreement with the W15 measurement of $|f_{\rm{R}0}| < 6 \times 10^{-5}$. Overall, these tests confirm the power of our methodology which can now be applied to larger cluster samples available with the next generation surveys.


title: Representations of Rota-Baxter algebras and regular singular decompositions

abstract: There is a Rota-Baxter algebra structure on the field $A=\mathbf{k}((t))$ with $ P$ being the projection map $A=\mathbf{k}[[t]]\oplus t^{-1}\mathbf{k}[t^{-1}]$ onto $ \mathbf{k}[[ t]]$. We study the representation theory and regular-singular decompositions of any finite dimensional $A$-vector space. The main result shows that the category of finite dimensional representations is semisimple and consists of exactly three isomorphism classes of irreducible representations which are all one-dimensional. As a consequence, the number of $GL_A(V)$-orbits in the set of all regular-singular decompositions of an $n$-dimensional $ A$-vector space $V$ is $(n+2)(n+1)/2$. We also use the result to compute the generalized class number, i.e., the number of the $GL_n(A)$-isomorphism classes of finitely generated $\mathbf{k}[[t]]$-submodules of $A^n$.


title: Suzaku Observation of the High-Inclination Binary EXO 0748-676 in the Hard State

abstract: Utilizing an archived Suzaku data acquired on 2007 December 25 for 46 ks, X-ray spectroscopic properties of the dipping and eclipsing low-mass X-ray binary EXO 0748$-$676 were studied. At an assumed distance of 7.1 kpc, the data gave a persistent unabsorbed luminosity of $3.4\times10^{36}$ erg cm$^{-2}$ s$^{-1}$ in 0.6 $-$ 55 keV. The source was in a relatively bright low/hard state, wherein the 0.6 $-$ 55 keV spectrum can be successfully explained by a "double-seed" Comptonization model, incorporating a common corona with an electron temperature of $\sim13$ keV. The seed photons are thought to be supplied from both the neutron star surface, and a cooler truncated disk. Compared to a sample of non-dipping low-mass X-ray binaries in the low/hard state, the spectrum is subject to stronger Comptonization, with a relatively larger Comptonizing $y$-parameter of $\sim1.4$ and a larger coronal optical depth of $\sim5$. This result, when attributed to the high inclination of EXO 0748$-$676, suggests that the Comptonizing corona may elongate along the disk plane, and give a longer path for the seed photons when viewed from edge-on inclinations.


title: Statistically validated network of portfolio overlaps and systemic risk

abstract: Common asset holding by financial institutions, namely portfolio overlap, is nowadays regarded as an important channel for financial contagion with the potential to trigger fire sales and thus severe losses at the systemic level. In this paper we propose a method to assess the statistical significance of the overlap between pairs of heterogeneously diversified portfolios, which then allows us to build a validated network of financial institutions where links indicate potential contagion channels due to realized portfolio overlaps. The method is implemented on a historical database of institutional holdings ranging from 1999 to the end of 2013, but can be in general applied to any bipartite network where the presence of similar sets of neighbors is of interest. We find that the proportion of validated network links (i.e., of statistically significant overlaps) increased steadily before the 2007-2008 global financial crisis and reached a maximum when the crisis occurred. We argue that the nature of this measure implies that systemic risk from fire sales liquidation was maximal at that time. After a sharp drop in 2008, systemic risk resumed its growth in 2009, with a notable acceleration in 2013, reaching levels not seen since 2007. We finally show that market trends tend to be amplified in the portfolios identified by the algorithm, such that it is possible to have an informative signal about financial institutions that are about to suffer (enjoy) the most significant losses (gains).


title: MSIQ: Joint Modeling of Multiple RNA-seq Samples for Accurate Isoform Quantification

abstract: Next-generation RNA sequencing (RNA-seq) technology has been widely used to assess full-length RNA isoform abundance in a high-throughput manner. RNA-seq data offer insight into gene expression levels and transcriptome structures, enabling us to better understand the regulation of gene expression and fundamental biological processes. Accurate isoform quantification from RNA-seq data is challenging due to the information loss in sequencing experiments. A recent accumulation of multiple RNA-seq data sets from the same tissue or cell type provides new opportunities to improve the accuracy of isoform quantification. However, existing statistical or computational methods for multiple RNA-seq samples either pool the samples into one sample or assign equal weights to the samples when estimating isoform abundance. These methods ignore the possible heterogeneity in the quality of different samples and could result in biased and unrobust estimates. In this article, we develop a method, which we call "joint modeling of multiple RNA-seq samples for accurate isoform quantification" (MSIQ), for more accurate and robust isoform quantification by integrating multiple RNA-seq samples under a Bayesian framework. Our method aims to (1) identify a consistent group of samples with homogeneous quality and (2) improve isoform quantification accuracy by jointly modeling multiple RNA-seq samples by allowing for higher weights on the consistent group. We show that MSIQ provides a consistent estimator of isoform abundance, and we demonstrate the accuracy and effectiveness of MSIQ compared with alternative methods through simulation studies on D. melanogaster genes. We justify MSIQ's advantages over existing approaches via application studies on real RNA-seq data from human embryonic stem cells, brain tissues, and the HepG2 immortalized cell line.


title: Riemannian geometry of the space of volume preserving immersions

abstract: Given a compact manifold $M$ and a Riemannian manifold $N$ of bounded geometry, we consider the manifold ${\rm Imm} (M,N)$ of immersions from $M$ to $N$ and its subset ${\rm Imm}_\mu (M,N)$ of those immersions with the property that the volume-form of the pull-back metric equals $\mu$. We first show that the non-minimal elements of ${\rm Imm}_\mu (M,N) $ form a splitting submanifold. On this submanifold we consider the Levi-Civita connection for various natural Sobolev metrics write down the geodesic equation and show local well-posedness in many cases. The question is a natural generalization of the corresponding well-posedness question for the group of volume-preserving diffeomorphisms, which is of great importance in fluid mechanics.


title: On the set of wild points of attracting surfaces in $\mathbb{R}^3$

abstract: Suppose that a closed surface $S \subseteq \mathbb{R}^3$ is an attractor, not necessarily global, for a discrete dynamical system. Assuming that its set of wild points $W$ is totally disconnected, we prove that (up to an ambient homeomorphism) it has to be contained in a straight line. Using this result and a modification of the classical construction of a wild sphere due to Antoine we show that there exist uncountably many different $2$--spheres in $\mathbb{R}^3$ none of which can be realized as an attractor for a homeomorphism.


title: Work Fluctuations in Bosonic Josephson Junctions

abstract: We calculate the first two moments and full probability distribution of the work performed on a system of bosonic particles in a two-mode Bose-Hubbard Hamiltonian when the self-interaction term is varied instantaneously or with a finite-time ramp. In the instantaneous case, we show how the irreversible work scales differently depending on whether the system is driven to the Josephson or Fock regime of the bosonic Josephson junction. In the finite-time case, we use optimal control techniques to substantially decrease the irreversible work to negligible values. Our analysis can be implemented in present-day experiments with ultracold atoms and we show how to relate the work statistics to that of the population imbalance of the two modes.


title: The Distributional Zeta-Function in Disordered Field Theory

abstract: In this paper we present a new mathematical rigorous technique for computing the average free energy of a disordered system with quenched randomness, using the replicas. The basic tool of this technique is a distributional zeta-function, a complex function whose derivative at the origin yields the average free energy of the system as the sum of two contributions: the first one is a series in which all the integer moments of the partition function of the model contribute; the second one, which can not be written as a series of the integer moments, can be made as small as desired. This result supports the use of integer moments of the partition function, computed via replicas, for expressing the average free energy of the system. One advantage of the proposed formalism is that it does not require the understanding of the properties of the permutation group when the number of replicas goes to zero. Moreover, the symmetry is broken using the saddle-point equations of the model. As an application for the distributional zeta-function technique, we obtain the average free energy of the disordered $\lambda\varphi^{4}$ model defined in a $d$-dimensional Euclidean space.


title: Partial and Total Dielectronic Recombination rate coefficients for W73+ to W56+

abstract: Dielectronic recombination (DR) is a key atomic process which affects the spectroscopic diagnostic modelling of tungsten, most of whose ionization stages will be found somewhere in the ITER fusion reactor: in the edge, divertor, or core plasma. Accurate DR data is sparse while complete DR coverage is unsophisticated (e.g. average-atom or Burgess General Formula) as illustrated by the large uncertainties which currently exist in the tungsten ionization balance. To this end, we present a series of partial final state-resolved and total DR rate coefficients for W73+ to W56+ Tungsten ions. This is part of a wider effort within The Tungsten Project to calculate accurate dielectronic recombination rate coefficients for the tungsten isonuclear sequence for use in collisional-radiative modelling of finite-density tokamak plasmas. The recombination rate coefficients have been calculated with autostructure using kappa-averaged relativistic wavefunctions in level resolution (intermediate coupling) and configuration resolution (configuration average). The results are available from OPEN-ADAS according to the adf09 and adf48 standard formats. Comparison with previous calculations of total DR rate coefficients for W63+ and W56+ yield agreement to within 20% and 10%, respectively, at peak temperature. It is also seen that the Juttner correction to the Maxwell distribution has a significant effect on the ionization balance of tungsten at the highest charge states, changing both the peak abundance temperatures and the ionization fractions of several ions.


title: Structural and Magnetic Characterization of Large Area, Free-Standing Thin Films of Magnetic Ion Intercalated Dichalcogenides Mn0.25TaS2 and Fe0.25TaS2

abstract: Free-standing thin films of magnetic ion intercalated transition metal dichalcogenides are produced using ultramicrotoming techniques. Films of thicknesses ranging from 30nm to 250nm were achieved and characterized using transmission electron diffraction and X-ray magnetic circular dichroism. Diffraction measurements visualize the long range crystallographic ordering of the intercalated ions, while the dichroism measurements directly assess the orbital contributions to the total magnetic moment. We thus verify the unquenched orbital moment in Fe0.25TaS2 and measure the fully quenched orbital contribution in Mn0.25TaS2. Such films can be used in a wide variety of ultrafast X-ray and electron techniques that benefit from transmission geometries, and allow measurements of ultrafast structural, electronic, and magnetization dynamics in space and time.


title: Empirical Evaluation of a Thread-Safe Dynamic Range Min-Max Tree using HTM

abstract: Succinct trees, such as wavelet trees and those based on, for instance, range Min-Max trees (RMMTs), are a family of practical data structures that store information close to their information-theoretic space lower bound. These structures are often static; meaning that once they are built, nodes cannot be added, deleted or modified. This read-only property simplifies concurrency. However, newer versions of these data structures allow for a fair degree of dynamism. Parallel programming using Hardware Transactional Memory(HTM), has been available in mainstream microprocessors since a few years ago. One limitation of HTM is still on the size of each transaction. This is why HTM's use, for the moment, is limited to operations that involve few memory addresses that need to be updated atomically, or where the level of concurrency is low. We provide the first available implementation of a concurrent, dynamic RMMT based on HTM, and we compare empirically how well HTM performs compared to a naive implementation using locks. We have shown that because of the formal properties of RMMTs, HTM is a good fit for adding concurrency to otherwise slow lock-based alternatives. We have also shown that HTM performs better than locks when the number of write operations increase, making it a practical structure to use in several write-intensive contexts. This is, as far as we know, the only practical implementation of RMMTs thoroughly tested using HTM.


title: Uncertainty Estimates for Theoretical Atomic and Molecular Data

abstract: Sources of uncertainty are reviewed for calculated atomic and molecular data that are important for plasma modeling: atomic and molecular structure and cross sections for electron-atom, electron-molecule, and heavy particle collisions. We concentrate on model uncertainties due to approximations to the fundamental many-body quantum mechanical equations and we aim to provide guidelines to estimate uncertainties as a routine part of computations of data for structure and scattering.


title: Rydberg polaritons in a thermal vapor

abstract: We present a pulsed four-wave mixing (FWM) scheme via a Rydberg state to create, store and retrieve collective Rydberg polaritons. The storage medium consists of a gas of thermal Rb atoms confined in a 220 {\mu}m thick cell, which are heated above room temperature. The experimental sequence consists of a pulsed excitation of Rydberg polaritons via the D1 line, a variable delay or storage time, and a final retrieval pulse via the D2 line. The lifetime of the Rydberg polaritons is around 1.2 ns, almost entirely limited by the excitation bandwidth and the corresponding motional dephasing of the atoms. The presented scheme combined with a tightly confined atomic ensemble is a good candidate for a deterministic single-photon source, as soon as strong interactions in terms of a Rydberg blockade are added.


title: The joint distribution of the Parisian ruin time and the number of claims until Parisian ruin in the classical risk model

abstract: In this paper we propose new iterative algorithm of calculating the joint distribution of the Parisian ruin time and the number of claims until Parisian ruin for the classical risk model. Examples are provided when the generic claim size is exponentially distributed.


title: Terahertz properties of Dirac fermions in HgTe films with optical doping

abstract: Terahertz properties of mercury telluride (HgTe) films with critical thickness are presented and discussed. The density of the charge carriers is controlled using contact-free optical doping by visible light. In the magneto-optical response of HgTe the contribution of two types of carriers (electrons and holes) can be identified. The density of the electrons can be modified by light illumination by more than one order of magnitude. As the hole density is roughly illuminationindependent, the terahertz response of the illuminated samples becomes purely electronic. In some cases, light illumination may switch the qualitative electrodynamic response from hole-like to the electron-like. The cyclotron mass of the electrons could be extracted from the data and shows a square root dependence upon the charge concentration in the broad range of parameters. This can be interpreted as a clear proof of a linear dispersion relations, i.e. Dirac-type charge carriers.


title: Shaken not stirred: Creating exotic angular momentum states by shaking an optical lattice

abstract: We propose a method to create higher orbital states of ultracold atoms in the Mott regime of an optical lattice. This is done by periodically modulating the position of the trap minima (known as shaking) and controlling the interference term of the lasers creating the lattice. These methods are combined with techniques of shortcuts to adiabaticity. As an example of this, we show specifically how to create an anti-ferromagnetic type ordering of angular momentum states of atoms. The specific pulse sequences are designed using Lewis-Riesenfeld invariants and a four-level model for each well. The results are compared with numerical simulations of the full Schroedinger equation.


title: Monoidal supercategories

abstract: In the literature, one finds several competing notions for the super (i.e., Z/2-graded) analog of a monoidal category. The goal of this paper is to clarify these definitions and the connections between them. We also discuss in detail the example of the odd Temperley-Lieb supercategory. In a forthcoming article, we will exploit the formalism developed here in order to define super analogs of the Kac-Moody 2-categories of Khovanov-Lauda and Rouquier.


title: Entanglement Entropy of the N=4 SYM spin chain

abstract: We present a detailed study of the Entanglement Entropy (EE) of excited states in all closed rank one subsectors of N=4 SYM, namely SU(2), SU(1|1) and SL(2). Exploiting the techniques of the Coordinate and the Algebraic Bethe Ansatz we obtain the EE for spin chains with up to seven magnons, at leading order in the coupling expansion but exact in the length of the spin chain and of the part of it that we cut. Focusing on the superconformal primary operator with two magnons in the BMN limit, we derive analytic and exact, in the coupling $\lambda'$, expressions for the Renyi and the EE. The interpolating functions for the Renyi and the EE monotonically increase as the coupling increases from the weak coupling $\lambda' \rightarrow 0$ regime to the strong coupling $\lambda' \rightarrow \infty$ regime. This results to a violation of a certain bound for the EE that is present at weak coupling and confirms the physical intuition that entanglement increases when the coupling increases.


title: Geometric Hypergraph Learning for Visual Tracking

abstract: Graph based representation is widely used in visual tracking field by finding correct correspondences between target parts in consecutive frames. However, most graph based trackers consider pairwise geometric relations between local parts. They do not make full use of the target's intrinsic structure, thereby making the representation easily disturbed by errors in pairwise affinities when large deformation and occlusion occur. In this paper, we propose a geometric hypergraph learning based tracking method, which fully exploits high-order geometric relations among multiple correspondences of parts in consecutive frames. Then visual tracking is formulated as the mode-seeking problem on the hypergraph in which vertices represent correspondence hypotheses and hyperedges describe high-order geometric relations. Besides, a confidence-aware sampling method is developed to select representative vertices and hyperedges to construct the geometric hypergraph for more robustness and scalability. The experiments are carried out on two challenging datasets (VOT2014 and Deform-SOT) to demonstrate that the proposed method performs favorable against other existing trackers.


title: Kinetic control of the coverage of oil droplets by DNA-functionalised colloids

abstract: We report a study of reversible adsorption of DNA-coated colloids on complementary functionalized oil droplets. We show that it is possible to control the surface coverage of oil droplets by colloidal particles, by exploiting the fact that during slow adsorption, compositional arrest takes place well before structural arrest occurs. As a consequence, we can prepare colloid-coated oil droplets with a `frozen' degree of loading, but with fully ergodic colloidal dynamics on the droplets. We illustrate the equilibrium nature of the adsorbed colloidal phase by exploring the quasi two-dimensional (2d) phase behaviour of the adsorbed colloids under the influence of depletion interactions. We present simulations of a simple model that illustrates the nature of the compositional arrest and the structural ergodicity.


title: Elastic vector and axial scattering of weakly interacting particles off nuclei

abstract: We analyze the elastic scattering of particles interacting with nuclei through vector and axial currents with overall couplings of the order of the Standard Model weak interaction, or smaller; the dominant contribution to the elastic cross section is identified as the coherent component and is therefore spin-independent. Differential and integrated cross sections are obtained for a wide range of incident particle masses and velocities and for nuclear targets with different masses; vector, axial and overall couplings of the incident particle and of the hadronic target to the massive exchanged boson are also kept general. This study naturally encompasses several kinds of possible dark matter components, including active and sterile neutrinos or neutralinos, and addresses the prospects for their direct detection through elastic scattering off nuclei.


title: Distributed Iterative Learning Control for a Team of Quadrotors

abstract: The goal of this work is to enable a team of quadrotors to learn how to accurately track a desired trajectory while holding a given formation. We solve this problem in a distributed manner, where each vehicle has only access to the information of its neighbors. The desired trajectory is only available to one (or few) vehicles. We present a distributed iterative learning control (ILC) approach where each vehicle learns from the experience of its own and its neighbors' previous task repetitions, and adapts its feedforward input to improve performance. Existing algorithms are extended in theory to make them more applicable to real-world experiments. In particular, we prove stability for any causal learning function with gains chosen according to a simple scalar condition. Previous proofs were restricted to a specific learning function that only depends on the tracking error derivative (D-type ILC). Our extension provides more degrees of freedom in the ILC design and, as a result, better performance can be achieved. We also show that stability is not affected by a linear dynamic coupling between neighbors. This allows us to use an additional consensus feedback controller to compensate for non-repetitive disturbances. Experiments with two quadrotors attest the effectiveness of the proposed distributed multi-agent ILC approach. This is the first work to show distributed ILC in experiment.


title: Mathematical modeling of T-cells activation dynamics for CD3 molecules recycling process

abstract: The down regulation of CD3 protein present on the surface of a T-cell triggers the process of activation. A rapid decrease of CD3 from the surface is required to reach the threshold that is necessary for the continuation of the activation process. By flow cytometry technique, it is well established that an amount of down regulated CD3 proteins recycle and may come on the surface of T-cell. This article is dedicated to study the impact of recycled CD3 protein that is recruited into the immune synapse at different rates. The population density of T-cells is distributed over surface protein (CD3) concentration under the framework of population balance models by using coupled system of PDEs. Such problems need non-conventional population balance modeling so that classical numerical technique can be applied. A recently developed efficient numerical method is validated against the analytical and numerical solutions.


title: Field Theory of the Eulerian Perfect Fluid

abstract: The Eulerian perfect-fluid theory is reformulated from its action principle in a pure field-theoretic manner. Conservation of the convective current is no longer imposed by Lin's constraints, but rather adopted as the central idea of the theory. Our formulation, for the first time, successfully reduces redundant degrees of freedom promoting one half of the Clebsch variables as the true dynamical fields. Interactions on these fields allow for the exchange of the convective current of quantities such as mass and charge, which are uniformly understood as the breaking of the underlying symmetry of the force-free fluid. The Clebsch fields play the essential role in the exchange of angular momentum with the force field producing vorticity.


title: Near field asymptotic behavior for the porous medium equation on the half-line

abstract: Kamin and V\'azquez proved in 1991 that solutions to the Cauchy-Dirichlet problem for the porous medium equation $u_t=(u^m)_{xx}$ on the half line with zero boundary data and nonnegative compactly supported integrable initial data behave for large times as a dipole type solution to the equation having the same first moment as the initial data, with an error which is $o(t^{-1/m})$. However, on sets of the form $0<x<g(t)$, with $g(t)=o(t^{1/(2m)})$ as $t\to\infty$, in the so called near field, the dipole solution is $o(t^{-1/m})$, and their result does not give neither the right rate of decay of the solution, nor a nontrivial asymptotic profile. In this paper we will show that the error is $o\big(t^{-(2m+1)/(2m^2)}(1+x)^{1/m}\big)$. This allows in particular to obtain a nontrivial asymptotic profile in the near field limit, which is a multiple of $x^{1/m}$, thus improving in this scale the results of Kamin and V\'azquez.


title: How to Combine a Billion Alphas

abstract: We give an explicit algorithm and source code for computing optimal weights for combining a large number N of alphas. This algorithm does not cost O(N^3) or even O(N^2) operations but is much cheaper, in fact, the number of required operations scales linearly with N. We discuss how in the absence of binary or quasi-binary clustering of alphas, which is not observed in practice, the optimization problem simplifies when N is large. Our algorithm does not require computing principal components or inverting large matrices, nor does it require iterations. The number of risk factors it employs, which typically is limited by the number of historical observations, can be sizably enlarged via using position data for the underlying tradables.


title: Efficient and powerful familywise error control in genome-wide association studies using generalized linear models

abstract: In genetic association studies, detecting phenotype-genotype association is a primary goal. We assume that the relationship between the data -phenotype, genetic markers and environmental covariates - can be modelled by a generalized linear model (GLM). The inclusion of environmental covariates makes it possible to account for important confounding factors, such as sex and population substructure. A multivariate score statistic, which under the complete null hypothesis of no phenotype-genotype association asymptotically has a multivariate normal distribution with a covariance matrix that can be estimated from the data, is used to test a large number of genetic markers for association with the phenotype. We stress the importance of controlling the familywise error rate (FWER), and use the asymptotic distribution of the multivariate score test statistic to find a local significance level for the individual test. Using real data (from one study on schizophrenia and bipolar disorder and one on maximal oxygen uptake) and constructed correlated structures, we show that our method is a powerful alternative to the popular Bonferroni and Sidak methods. For GLMs without environmental covariates, we show that our method is an efficient alternative to permutation methods for multiple testing. Further, we show that if environmental covariates and genetic markers are uncorrelated, the estimated covariance matrix of the score test statistic can be approximated by the estimated correlation matrix for just the genetic markers. As byproducts of our method, an effective number of independent tests can be defined, and FWER-adjusted $p$-values can be calculated as an alternative to using a local significance level.


title: Online Distributed Scheduling on a Fault-prone Parallel System

abstract: We consider a parallel system of $m$ identical machines prone to unpredictable crashes and restarts, trying to cope with the continuous arrival of tasks to be executed. Tasks have different computational requirements (i.e., processing time or size). The flow of tasks, their size, and the crash and restart of the machines are assumed to be controlled by an adversary. Then, we focus on the study of online distributed algorithms for the efficient scheduling of the tasks. We use competitive analysis, considering as efficiency metric the completed-load, i.e., the aggregated size of the completed tasks. We first present optimal completed-load competitiveness algorithms when the number of different task sizes that can be injected by the adversary is bounded. (It is known that, if it is not bounded, competitiveness is not achievable.) We first consider only two different task sizes, and then proceed to $k$ different ones, showing in both cases that the optimal completed-load competitiveness can be achieved. Then, we consider the possibility of having some form of resource augmentation, allowing the scheduling algorithm to run with a speedup $s \geq 1$. In this case, we show that the competitiveness of all work-conserving scheduling algorithms can be increased by using a large enough speedup.


title: Current of interacting particles inside a channel of exponential cavities: Application of a modified Fick--Jacobs equation

abstract: The Fick--Jacobs equation has been widely studied, because of its applications in the diffusion and transport of non-interacting particles in narrow channels. It is also known that a modified version of this equation can be used to describe the same system with particles interacting through a hard-core potential. In this work we present a system that can be exactly solved using the Fick--Jacobs equation. The exact results of the particle concentration profile along the channel $n$, the current, $J$, and the mobility, $\mu$, of particles as a function of an external force are contrasted with Monte Carlo simulations results of non-interacting particles. For interacting particles the behavior of $n$, $J$ and $\mu$, obtained from the modified Fick--Jacobs equation are in agreement with numerical simulations, where the hard-core interaction is taken into account. Even more, for interacting particles the modified Fick--Jacobs equation gives comparatively more accurate results of the current difference (when a force is applied in opposite directions) than the exact result for the non-interacting ones.


title: Uplink FBMC/OQAM-based Multiple Access Channel: Distortion Analysis under Strong Frequency Selectivity

abstract: This paper computes the distortion power at the receiver side of an FBMC/OQAM-based OFDMA uplink channel under strong frequency selectivity and/or user timing errors. More precisely, it provides a distortion expression that is valid for a wide class of prototype pulses (not necessarily perfect-reconstruction ones) when the number of subcarriers is sufficiently large. This result is a valuable instrument for analyzing how users interfere to one another and to justify, formally, the common choice of placing an empty guard band between adjacent users. Interestingly, the number of out-band subcarriers contaminated by each user only depends on the prototype pulses and not on the channel nor on the equalizer. To conclude, the distortion analysis presented in this paper, together with some simulation results for a realistic scenario, also provide convincing evidence that FBMC/OQAM-based OFDMA is superior to classic circular-prefix OFDMA in the case of asynchronous users.


title: Report on the Poblacion Flotante of Bogota (D.C.)

abstract: In this document we describe the size of the Poblacion Flotante of Bogota (D.C.). The Poblacion Flotante is composed by people who live outside Bogota (D.C.), but who rely on the city for performing their job. We estimate the Poblacion Flotante impact relying on a new data source provided by telecommunications operators in Colombia, which enables us to estimate how many people commute daily from every municipality of Colombia to a specific area of Bogota (D.C.). We estimate that the size of the Poblacion Flotante could represent a 5.4% increase of Bogota (D.C.)'s population. During weekdays, the commuters tend to visit the city center more.


title: Nanbu-Goto action and qubit theory in any signature and higher dimensions

abstract: We perform an extension of the relation between the Nambu-Goto action and qubit theory. Of course, the Cayley hyperdeterminant is the key mathematical tool in such generalization. Using the Wick rotation we find that in four dimensions such a relation can be established no only in (2+2)-dimensions but also in any signature. We generalize our result to a curved space-time of (2$^{2n}$+2$^{2n}$)-dimensions and (2$^{2n+1}$+2$^{2n+1}$)-dimensions.


title: Lower Bounds for Graph Exploration Using Local Policies

abstract: We give lower bounds for various natural node- and edge-based local strategies for exploring a graph. We consider this problem both in the setting of an arbitrary graph as well as the abstraction of a geometric exploration of a space by a robot, both of which have been extensively studied. We consider local exploration policies that use time-of-last- visit or alternatively least-frequently-visited local greedy strategies to select the next step in the exploration path. Both of these strategies were previously considered by Cooper et al. (2011) for a scenario in which counters for the last visit or visit frequency are attached to the edges. In this work we consider the case in which the counters are associated with the nodes, which for the case of dual graphs of geometric spaces could be argued to be intuitively more natural and likely more efficient. Surprisingly, these alternate strategies give worst-case superpolynomial/ exponential time for exploration, whereas the least-frequently visited strategy for edges has a polynomially bounded exploration time, as shown by Cooper et al. (2011).


title: Parameterized vertex deletion problems for hereditary graph classes with a block property

abstract: For a class of graphs $\mathcal{P}$, the Bounded $\mathcal{P}$-Block Vertex Deletion problem asks, given a graph $G$ on $n$ vertices and positive integers $k$ and $d$, whether there is a set $S$ of at most $k$ vertices such that each block of $G-S$ has at most $d$ vertices and is in $\mathcal{P}$. We show that when $\mathcal{P}$ satisfies a natural hereditary property and is recognizable in polynomial time, Bounded $\mathcal{P}$-Block Vertex Deletion can be solved in time $2^{O(k \log d)}n^{O(1)}$. When $\mathcal{P}$ contains all split graphs, we show that this running time is essentially optimal unless the Exponential Time Hypothesis fails. On the other hand, if $\mathcal{P}$ consists of only complete graphs, or only cycle graphs and $K_2$, then Bounded $\mathcal{P}$-Block Vertex Deletion admits a $c^{k}n^{O(1)}$-time algorithm for some constant $c$ independent of $d$. We also show that Bounded $\mathcal{P}$-Block Vertex Deletion admits a kernel with $O(k^2 d^7)$ vertices.


title: Existence of multi-bump solutions to biharmonic operator with critical exponential growth in $\mathbb{R}^4$

abstract: Using variational methods, we establish existence of multi-bump solutions for the following class of problems $$ \left\{ \begin{array}{l} \Delta^2 u +(\lambda V(x)+1)u = f(u), \quad \mbox{in} \quad \mathbb{R}^{4}, u \in H^{2}(\mathbb{R}^{4}), \end{array} \right. $$ where $\Delta^2$ is the biharmonic operator, $f$ is a continuous function with critical exponential growth and $V : \mathbb{R}^4 \rightarrow \mathbb{R}$ is a continuous function verifying some conditions.


title: Noisy Hypotheses in the Age of Discovery Science

abstract: We draw attention to one specific issue raised by Ioannidis (2005), that of very many hypotheses being tested in a given field of investigation. To better isolate the problem that arises in this (massive) multiple testing scenario, we consider a utopian setting where the hypotheses are tested with no additional bias. We show that, as the number of hypotheses being tested becomes much larger than the discoveries to be made, it becomes impossible to reliably identify true discoveries. This phenomenon, well-known to statisticians working in the field of multiple testing, puts in jeopardy any naive pursuit in (pure) discovery science.


title: Evaluating Generating Functions for Periodic Multiple Polylogarithms

abstract: The goal of the paper is to give a systematic way to numerically evaluate the generating function of a periodic multiple polylogarithm using a Chen-Fliess series with a rational generating series. The idea is to realize the corresponding Chen-Fliess series as a bilinear dynamical system. A standard form for such a realization is given. The method is also generalized to the case where the multiple polylogarithm has non-periodic components. This allows one, for instance, to numerically validate the Hoffman conjecture. Finally, a setting in terms of dendriform algebras is provided.


title: Problem optymalizacyjny de Finettiego dla proces\'ow L\'evy'ego

abstract: We give a review of the state of the art with regard to the dividend problem.


title: Quasinormal modes and the phase structure of strongly coupled matter

abstract: We investigate the poles of the retarded Green's functions of strongly coupled field theories exhibiting a variety of phase structures from a crossover up to different first order phase transitions. These theories are modeled by a dual gravitational description. The poles of the holographic Green's functions appear at the frequencies of the quasinormal modes of the dual black hole background. We focus on quantifying linearized level dynamical response of the system in the critical region of phase diagram. Generically non-hydrodynamic degrees of freedom are important for the low energy physics in the vicinity of a phase transition. For a model with linear confinement in the meson spectrum we find degeneracy of hydrodynamic and non-hydrodynamic modes close to the minimal black hole temperature, and we establish a region of temperatures with unstable non-hydrodynamic modes in a branch of black hole solutions.


title: Model Comparison of a Data-Driven and a Physical Model for Simulating HVAC Systems

abstract: Commercial buildings are responsible for a large fraction of energy consumption in developed countries, and therefore are targets of energy efficiency programs. Motivated by the large inherent thermal inertia of buildings, the power consumption can be flexibly scheduled without compromising occupant comfort. This temporal flexibility offers opportunities for the provision of frequency regulation to support grid stability. To realize energy savings and frequency regulation, it is of prime importance to identify a realistic model for the temperature dynamics of a building. We identify a low- dimensional data-driven model and a high-dimensional physics- based model for different spatial granularities and temporal seasons based on a case study of an entire floor of Sutardja Dai Hall, an office building on the University of California, Berkeley campus. A comparison of these contrasting models shows that, despite the higher forecasting accuracy of the physics-based model, both models perform almost equally well for energy efficient control. We conclude that the data-driven model is more amenable to controller design due to its low complexity, and could serve as a substitution for highly complex physics- based models with an insignificant loss of prediction accuracy for many applications. On the other hand, our physics-based approach is more suitable for modeling buildings with finer spatial granularities.


title: The Lepton Flavour Violating Higgs Decays at the HL-LHC and the ILC

abstract: Run-I results from the CMS collaboration show an excess of events in the decay $h \to \mu \tau_e$ with a local significances of 2.4$\sigma$. This could be the first hint of flavour violation in the Higgs sector. We summarise the bounds on the flavour violating Yukawa couplings from direct searches, low energy measurements and projected future experiments. We discuss the sensitivity of upcoming HL-LHC runs and future lepton colliders in measuring lepton-flavour violating couplings using an effective field theory framework. For the HL-LHC we find limits on BR($h\to \mu \tau$) and BR$(h\to e \tau) \lesssim \mathcal{O}(0.5)\%$ and on BR($h \to e \mu) \lesssim \mathcal{O}(0.02)\%$. For an ILC with center-of-mass energy of 1 TeV we expect BR($h \to e \tau$) and BR($h\to \mu \tau$) to be measurable down to $\mathcal{O}(0.2)\%$.


title: Katyusha: The First Direct Acceleration of Stochastic Gradient Methods

abstract: Nesterov's momentum trick is famously known for accelerating gradient descent, and has been proven useful in building fast iterative algorithms. However, in the stochastic setting, counterexamples exist and prevent Nesterov's momentum from providing similar acceleration, even if the underlying problem is convex and finite-sum. We introduce $\mathtt{Katyusha}$, a direct, primal-only stochastic gradient method to fix this issue. In convex finite-sum stochastic optimization, $\mathtt{Katyusha}$ has an optimal accelerated convergence rate, and enjoys an optimal parallel linear speedup in the mini-batch setting. The main ingredient is $\textit{Katyusha momentum}$, a novel "negative momentum" on top of Nesterov's momentum. It can be incorporated into a variance-reduction based algorithm and speed it up, both in terms of $\textit{sequential and parallel}$ performance. Since variance reduction has been successfully applied to a growing list of practical problems, our paper suggests that in each of such cases, one could potentially try to give Katyusha a hug.


title: The structure of combinatorial Markov processes

abstract: Every exchangeable Feller process taking values in a suitably nice combinatorial state space can be constructed by a system of iterated random Lipschitz functions. In discrete time, the construction proceeds by iterated application of independent, identically distributed functions, while in continuous time the random functions occur as the atoms of a time homogeneous Poisson point process. We further show that every exchangeable Feller process projects to a Feller process in an appropriate limit space, akin to the projection of partition-valued processes into the ranked-simplex and graph-valued processes into the space of graph limits. Together, our main theorems establish common structural features shared by all exchangeable combinatorial Feller processes, regardless of the dynamics or resident state space, thereby generalizing behaviors previously observed for exchangeable coalescent and fragmentation processes as well as other combinatorial stochastic processes. If, in addition, an exchangeable Feller process evolves on a state space satisfying the $n$-disjoint amalgamation property for all $n\geq1$, then its jump measure can be decomposed explicitly in the sense of L\'evy--It\^o--Khintchine.


title: Transferring Learned Microcalcification Group Detection from 2D Mammography to 3D Digital Breast Tomosynthesis Using a Hierarchical Model and Scope-based Normalization Features

abstract: A novel hierarchical model is introduced to solve a general problem of detecting groups of similar objects. Under this model, detection of groups is performed in hierarchically organized layers while each layer represents a scope for target objects. The processing of these layers involves sequential extraction of appearance features for an individual object, consistency measurement features for nearby objects, and finally the distribution features for all objects within the group. Using the concept of scope-based normalization, the extracted features not only enhance local contrast of an individual object, but also provide consistent characterization for all related objects. As an example, a microcalcification group detection system for 2D mammography was developed, and then the learned model was transferred to 3D digital breast tomosynthesis without any retraining or fine-tuning. The detection system demonstrated state-of-the-art performance and detected 96% of cancerous lesions at the rate of 1.2 false positives per volume as measured on an independent tomosynthesis test set.


title: Lie subalgebras of the matrix quantum pseudo differential operators

abstract: We give a complete description of the anti-involutions that preserve the principal gradation of the algebra of matrix quantum pseudodifferential operators and we describe the Lie subalgebras of its minus fixed points.


title: Solutions with a bounded support promote permanence of a distributed replicator equation

abstract: The now classical replicator equation describes a wide variety of biological phenomena, including those in theoretical genetics, evolutionary game theory, or in the theories of the origin of life. Among other questions, the permanence of the replicator equation is well studied in the local, well-mixed case. Inasmuch as the spatial heterogeneities are key to understanding the species coexistence at least in some cases, it is important to supplement the classical theory of the non-distributed replicator equation with a spatially explicit framework. One possible approach, motivated by the porous medium equation, is introduced. It is shown that the solutions to the spatially heterogeneous replicator equation may evolve to equilibrium states that have a bounded support, and, moreover, that these solutions are of paramount importance for the overall system permanence, which is shown to be a more commonplace phenomenon for the spatially explicit equation if compared with the local model.


title: Tutorial to SARAH

abstract: I give in this brief tutorial a short practical introduction to the Mathematica package SARAH. First, it is shown how an existing model file can be changed to implement a new model in SARAH. In the second part, masses, vertices and renormalisation group equations are calculated with SARAH. Finally, the main commands to generate model files and output for other tools are summarised.


title: Efficient Multi-Scale 3D CNN with Fully Connected CRF for Accurate Brain Lesion Segmentation

abstract: We propose a dual pathway, 11-layers deep, three-dimensional Convolutional Neural Network for the challenging task of brain lesion segmentation. The devised architecture is the result of an in-depth analysis of the limitations of current networks proposed for similar applications. To overcome the computational burden of processing 3D medical scans, we have devised an efficient and effective dense training scheme which joins the processing of adjacent image patches into one pass through the network while automatically adapting to the inherent class imbalance present in the data. Further, we analyze the development of deeper, thus more discriminative 3D CNNs. In order to incorporate both local and larger contextual information, we employ a dual pathway architecture that processes the input images at multiple scales simultaneously. For post-processing of the network's soft segmentation, we use a 3D fully connected Conditional Random Field which effectively removes false positives. Our pipeline is extensively evaluated on three challenging tasks of lesion segmentation in multi-channel MRI patient data with traumatic brain injuries, brain tumors, and ischemic stroke. We improve on the state-of-the-art for all three applications, with top ranking performance on the public benchmarks BRATS 2015 and ISLES 2015. Our method is computationally efficient, which allows its adoption in a variety of research and clinical settings. The source code of our implementation is made publicly available.


title: Solar models in light of new high metallicity measurements from solar wind data

abstract: We study the impact of new metallicity measurements, from solar wind data, on the solar model. The "solar modelling problem" refers to the persisting discrepancy between helioseismological observations and predictions of solar models computed implementing state-of-the-art photospheric abundances. We critically reassess the problem, in particular considering the new set of abundances of von Steiger \& Zurbuchen 2016, determined through the \textit{in situ} collection of solar wind samples from polar coronal holes. This new set of abundances indicates a solar metallicity $Z_{\odot} \geq 0.0196 \pm 0.0014$, significantly higher than the currently established value. The new values hint at an abundance of volatile elements (i.e. C, N, O, Ne) close to previous results of Grevesse \& Sauval 1998, whereas the abundance of refractory elements (i.e. Mg, Si, S, Fe) is considerably increased. Using the Linear Solar Model formalism, we determine the variation of helioseismological observables in response to the changes in elemental abundances, in order to explore the consistency of these new measurements with constraints from helioseismology. We find that, for observables that are particularly sensitive to the abundance of volatile elements, in particular the radius of the convective zone boundary (CZB) and the sound speed around the radius of CZB, improved agreement over previous models is obtained. Conversely, the high abundance of refractories correlates with a higher core temperature, resulting in an overproduction of neutrinos and a huge increase in the surface Helium abundance. We conclude that the "solar modelling problem" remains unsolved.


title: Plasma Relaxation and Topological Aspects in Electronmagnetohydrodynamics

abstract: Parker's formulation of isotopological plasma relaxation process toward minimum magnetics energy states in magnetohydrodynamics (MHD) is extended to electron MHD (EMHD). The lower bound on magnetic energy in EMHD is determined by both the magnetic field and the electron vorticity field topologies, and is shown to be reduced further in EMHD by an amount proportional to the sum of total electron-flow kinetic energy and total electron-flow enstrophy. The EMHD Beltrami condition becomes equivalent to the potential vorticity conservation equation in two-dimensional (2D) hydrodynamics, and the torsion coefficient and turns out to be proportional to potential vorticity. The winding pattern of the magnetic field lines appears to evolve therefore in the same way as "potential vorticity" lines in 2D hydrodynamics.


title: Document Neural Autoregressive Distribution Estimation

abstract: We present an approach based on feed-forward neural networks for learning the distribution of textual documents. This approach is inspired by the Neural Autoregressive Distribution Estimator(NADE) model, which has been shown to be a good estimator of the distribution of discrete-valued igh-dimensional vectors. In this paper, we present how NADE can successfully be adapted to the case of textual data, retaining from NADE the property that sampling or computing the probability of observations can be done exactly and efficiently. The approach can also be used to learn deep representations of documents that are competitive to those learned by the alternative topic modeling approaches. Finally, we describe how the approach can be combined with a regular neural network N-gram model and substantially improve its performance, by making its learned representation sensitive to the larger, document-specific context.


title: You Really Need A Good Ruler to Measure Caching Performance in Information-Centric Networks

abstract: Information-centric networks are an interesting new paradigm for distributing content on the Internet. They bring up many research challenges, such as addressing content by name, securing content, and wide-spread caching of content. Caching has caught a lot of attention in the research community, but a lot of the work suffers from a poor understanding of the different metrics with which caching performance can be measured. In this paper we not only present a comprehensive overview of different caching metrics that have been proposed for information-centric networks, but also propose the coupling factor as a new metric to capture the relation- ship between content popularity and network topology. As we show, many commonly used metrics have several failure modes which are largely ignored in literature. We identify these problems and propose remedies and new metrics to address these failures. Our work highlights the fundamental differences between information-centric caches and "traditional" cache networks and we demonstrate the need for a systematic understanding of the metrics for information- centric caching. We also discuss how experimental work should be done when evaluating networks of caches.


title: Electro-Optical Radiation of a Charged Particle in a Single Laser Beam Field

abstract: A peculiar radiation arising as a result of radiation interference of nonlinear oscillators excited by a monochromatic plane wave field of the incident particle is described. The radiation properties are determined by the fact that a phase of each oscillator radiation fields is synchronized by a wave field, while the radiation itself occurs due to the particle field influence on the oscillators. The consideration is performed for a thin film with negligible density effect. It is supposed that the contribution is given only by a long-wave part of the Weizsacker spectrum for which nonlinear polarization coefficients of medium are large.


title: Quantum Theory of Half-integer Spin Free Particles from the Perspective of the Majorana Equation

abstract: In this study, the Majorana equation for particles with arbitrary spin is solved for a half-integer spin free particle. The solution for the fundamental state, corresponding to the reference frame in which the particle is at rest, is compared with that obtained using the Dirac equation, especially as regards the approximation in the relativistic limit, in which the speed of the particle is close to that of light. Furthermore, the solutions that Majorana defines unphysical, proving that their occupation probability increases with the particle velocity, are taken into consideration. The anomalous behavior exhibited by these states also shows that for high-energy particles with small mass, transitions from a bradyonic state to a tachyonic state become possible.


title: Distributed Randomized Control for Demand Dispatch

abstract: The paper concerns design of control systems for Demand Dispatch to obtain ancillary services to the power grid by harnessing inherent flexibility in many loads. The role of "local intelligence" at the load has been advocated in prior work, randomized local controllers that manifest this intelligence are convenient for loads with a finite number of states. The present work introduces two new design techniques for these randomized controllers: (i) The Individual Perspective Design (IPD) is based on the solution to a one-dimensional family of Markov Decision Processes, whose objective function is formulated from the point of view of a single load. The family of dynamic programming equation appears complex, but it is shown that it is obtained through the solution of a single ordinary differential equation. (ii) The System Perspective Design (SPD) is motivated by a single objective of the grid operator: Passivity of any linearization of the aggregate input-output model. A solution is obtained that can again be computed through the solution of a single ordinary differential equation. Numerical results complement these theoretical results.


title: A computational high-throughput search for new ternary superalloys

abstract: In 2006, a novel cobalt-based superalloy was discovered [1] with mechanical properties better than some conventional nickel-based superalloys. As with conventional superalloys, its high performance arises from the precipitate-hardening effect of a coherent L1$_2$ phase, which is in two-phase equilibrium with the fcc matrix. Inspired by this unexpected discovery of an L1$_2$ ternary phase, we performed a first-principles search through 2224 ternary metallic systems for analogous precipitate-hardening phases of the form $X_{3}$[$A_{0.5}, B_{0.5}$], where $X$ = Ni, Co, or Fe, and [$A,B$] = Li, Be, Mg, Al, Si, Ca, Sc, Ti, V, Cr, Mn, Fe, Co, Ni, Cu, Zn Ga, Sr, Y, Zr, Nb, Mo, Tc, Ru, Rh, Pd, Ag, Cd, In, Sn, Sb, Hf, Ta, W, Re, Os, Ir, Pt, Au, Hg, or Tl. We found 102 systems that have a smaller decomposition energy and a lower formation enthalpy than the Co$_{3}$(Al, W) superalloy. They have a stable two-phase equilibrium with the host matrix within the concentration range $0<x<1$ ($X_{3}$[$A_{x}, B_{1-x}$]) and have a relative lattice mismatch with the host matrix of less than or equal to 5%. These new candidates, narrowed from 2224 systems, suggest possible experimental exploration for identifying new superalloys. Of these 102 systems, 37 are new; they have no reported phase diagrams in standard databases. Based on cost, experimental difficulty, and toxicity, we limit these 37 to a shorter list of six promising candidates of immediate interest. Our calculations are consistent with current experimental literature where data exists.


title: On the origin of the CDDR violation

abstract: The investigation of any violation in the Cosmic Distance Duality Relation (CDDR) is one of the most important sources of investigation for a new physic. In this paper we propose a new method to find the origin of a possible violation on the CDDR. Such violation is defined from the equation $\frac{d_L}{d_A(1+z)^2} = \eta$, with $\eta \neq 1$. We analyze the observational constraints from SNIa, BAO and $f_{gas}$ for the parameter $\eta$ using three different parameterizations. We create three data sets with the following combinations: (Set I) $d_L$ data from SNIa and $d_A$ data from BAO measurements, (Set II) $d_L$ data from $f_{gas}$ and $d_A$ from BAO measurements and (Set III) $d_L$ from SNIa and $d_A$ from $f_{gas}$ measurements. The sets have 18 points and a redshift difference between the observational pairs of $\Delta z \leq 0.08$. It is found that the difference between Sets I and II is up to five times higher than Sets I and III, what suggest that a violation in the CDDR come from a new physics related with the angular diameter distance.


title: Circuits and Hurwitz action in finite root systems

abstract: In a finite real reflection group, two factorizations of a Coxeter element into an arbitrary number of reflections are shown to lie in the same orbit under the Hurwitz action if and only if they use the same multiset of conjugacy classes. The proof makes use of a surprising lemma, derived from a classification of the minimal linear dependences (matroid circuits) in finite root systems: any set of roots forming a minimal linear dependence with positive coefficients has a disconnected graph of pairwise acuteness.


title: Coherent-state constellations and polar codes for thermal Gaussian channels

abstract: Optical communication channels are ultimately quantum-mechanical in nature, and we must therefore look beyond classical information theory to determine their communication capacity as well as to find efficient encoding and decoding schemes of the highest rates. Thermal channels, which arise from linear coupling of the field to a thermal environment, are of particular practical relevance; their classical capacity has been recently established, but their quantum capacity remains unknown. While the capacity sets the ultimate limit on reliable communication rates, it does not promise that such rates are achievable by practical means. Here we construct efficiently encodable codes for thermal channels which achieve the classical capacity and the so-called Gaussian coherent information for transmission of classical and quantum information, respectively. Our codes are based on combining polar codes with a discretization of the channel input into a finite "constellation" of coherent states. Encoding of classical information can be done using linear optics.


title: Polaritonic Rabi and Josephson Oscillations

abstract: The dynamics of coupled condensates is a wide-encompassing problem with relevance to superconductors, BECs in traps, superfluids, etc. Here, we provide a unified picture of this fundamental problem that includes i) detuning of the free energies, ii) different self-interaction strengths and iii) finite lifetime of the modes. At such, this is particularly relevant for the dynamics of polaritons, both for their internal dynamics between their light and matter constituents, as well as for the more conventional dynamics of two spatially separated condensates. Polaritons are short-lived, interact only through their material fraction and are easily detuned. At such, they bring several variations to their atomic counterpart. We show that the combination of these parameters results in important twists to the phenomenology of the Josephson effect, such as the behaviour of the relative phase (running or oscillating) or the occurrence of self-trapping. We undertake a comprehensive stability analysis of the fixed points on a normalized Bloch sphere, that allows us to provide a generalized criterion to identify the Rabi and Josephson regimes in presence of detuning and decay.


title: Anomalous diffusion and FRAP dynamics in the random comb model

abstract: We address the problem of diffusion on a comb whose teeth display a varying length. Specifically, the length $\ell$ of each tooth is drawn from a probability distribution displaying the large-$\ell$ behavior $P(\ell) \sim \ell^{-(1+\alpha)}$ ($\alpha>0$). Our method is based on the mean-field description provided by the well-tested CTRW approach for the random comb model, and the obtained analytical result for the diffusion coefficient is confirmed by numerical simulations. We subsequently incorporate retardation effects arising from binding/unbinding kinetics into our model and obtain a scaling law characterizing the corresponding change in the diffusion coefficient. Finally, our results for the diffusion coefficient are used as an input to compute concentration recovery curves mimicking FRAP experiments in comb-like geometries such as spiny dendrites. We show that such curves cannot be fitted perfectly by a model based on scaled Brownian motion, i.e., a standard diffusion equation with a time-dependent diffusion coefficient. However, differences between the exact curves and such fits are small, thereby providing justification for the practical use of models relying on scaled Brownian motion as a fitting procedure for recovery curves arising from particle diffusion in comb-like systems.


title: Producing Virtually Defect Free Nanoscale Ripples by Ion Bombardment of Rocked Solid Surfaces

abstract: Bombardment of a solid surface with a broad, obliquely-incident ion beam frequently produces nanoscale surface ripples. The primary obstacle that prevents the adoption of ion bombardment as a nano-fabrication tool is the high density of defects in the patterns that are typically formed. Our simulations indicate that ion bombardment can produce nearly defect free ripples on the surface of an elemental solid if the sample is concurrently and periodically rocked about an axis orthogonal to the surface normal and the incident beam direction. We also investigate the conditions necessary for rocking to produce highly ordered ripples and discuss how the results of our simulations can be reproduced experimentally.


title: First-principles study of ferroelectricity induced by p-d hybridization in ferrimagnetic NiFe2O4

abstract: We investigate the ferrimagnetism and ferroelectricity of bulk NiFe$_2$O$_4$ with tetragonal $P4_122$ ~symmetry by means of density functional calculations using generalized gradient approximation + Hubbard $U$ approach. Special attention is paid to finding the most energetically favorable configuration on magnetic ordering and further calculating the reliable spontaneous electric polarization. With the fully optimized crystalline structure of the most stable configuration, the spontaneous polarization is obtained to be 23 $\mu$C/cm$^2$ along the z direction, which originates from the hybridization between the 3d states of the Fe$^{3+}$ cation and the 2p states of oxygen induced by Jahn-Teller effect.


title: Gravitational wave emission by the high braking index pulsar PSR J1640-4631

abstract: Recently, a braking index for the pulsar PSR J1640-4631 has been measured. With a braking index of $n = 3.15 \pm 0.03$, this pulsar has the highest braking index ever measured. As it is well known, a pure magnetic dipole brake yields $n = 3$, whereas a pure gravitational wave brake yields $n = 5$. Therefore, each of these mechanisms alone can not account for the braking index found for PSR J1640-4631. Here we show that such a braking index can be accounted for if the spindown model combines magnetic dipole and gravitational wave brakes. Then, we briefly discuss the detectability of this pulsar by aLIGO and the planned Einstein Telescope. In particular, we show that the amplitude of the gravitational wave that comes from our model is around a factor four lower than the amplitude modeled exclusively by gravitational wave energy loss. Another interesting outcome of our modeling is that we are able to obtain the ellipticity from the braking index and other pulsar parameters.


title: BICEP2 / Keck Array VII: Matrix based E/B Separation applied to BICEP2 and the Keck Array

abstract: A linear polarization field on the sphere can be uniquely decomposed into an E-mode and a B-mode component. These two components are analytically defined in terms of spin-2 spherical harmonics. Maps that contain filtered modes on a partial sky can also be decomposed into E-mode and B-mode components. However, the lack of full sky information prevents orthogonally separating these components using spherical harmonics. In this paper, we present a technique for decomposing an incomplete map into E and B-mode components using E and B eigenmodes of the pixel covariance in the observed map. This method is found to orthogonally define E and B in the presence of both partial sky coverage and spatial filtering. This method has been applied to the BICEP2 and the Keck Array maps and results in reducing E to B leakage from LCDM E-modes to a level corresponding to a tensor-to-scalar ratio of $r<1\times10^{-4}$.


title: Connecting the new H.E.S.S. diffuse emission at the Galactic center with the Fermi GeV excess: a combination of millisecond pulsars and heavy dark matter?

abstract: The H.E.S.S. collaboration has reported a high-energy spherically symmetric diffuse gamma-ray emission in the inner 50 pc of the Milky Way, up to ~ 50 TeV. Here we propose a leptonic model which provides an alternative to the hadronic scenario presented by the H.E.S.S. collaboration, and connects the newly reported TeV emission to the Fermi-LAT Galactic center GeV excess. Our model relies on a combination of inverse Compton emission from a population of millisecond pulsars---which can account for the GeV excess---and a supermassive black hole-induced spike of heavy (~ 60 TeV) dark matter particles annihilating into electrons with a sub-thermal cross-section. With an up-to-date interstellar radiation field, as well as a standard magnetic field and diffusion set-up, our model accounts for the spectral morphology of the detected emission. Moreover, we show that the dark matter induced emission reproduces the spatial morphology of the H.E.S.S. signal above ~ 10 TeV, while we obtain a slightly more extended component from pulsars at lower energies, which could be used as a prediction for future H.E.S.S. observations.

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


title: Diphoton Resonance from a Warped Extra Dimension

abstract: We argue that extensions of the Standard Model (SM) with a warped extra dimension, which successfully address the hierarchy and flavor problems of elementary particle physics, can provide an elegant explanation of the 750 GeV diphoton excess recently reported by ATLAS and CMS. A gauge-singlet bulk scalar with ${\cal O}(1)$ couplings to fermions is identified as the new resonance $S$, and the vector-like Kaluza-Klein excitations of the SM quarks and leptons mediate its loop-induced couplings to photons and gluons. The electroweak gauge symmetry almost unambiguously dictates the bulk matter content and hence the hierarchies of the $S\to \gamma\gamma$, $WW$, $ZZ$, $Z\gamma$, $t\bar t$ and dijet decay rates. We find that the $S\to Z\gamma$ decay mode is strongly suppressed, such that $\mbox{Br}(S\to Z\gamma)/\mbox{Br}(S\to\gamma\gamma)<0.1$. The hierarchy problem for the new scalar boson is solved in analogy with the Higgs boson by localizing it near the infrared brane. The infinite sums over the Kaluza-Klein towers of fermion states are finite and can be calculated in closed form with a remarkably simple result. Reproducing the observed $pp\to S\to\gamma\gamma$ signal requires Kaluza-Klein masses in the multi-TeV range, consistent with bounds from flavor physics and electroweak precision observables. Useful side products of our analysis, which can be adapted to almost any model for the diphoton resonance, are the calculation of the gluon-fusion production cross section $\sigma(pp\to S)$ at NNLO in QCD, an exact expression for the inclusive $S\to gg$ decay rate at N$^3$LO, a study of the $S\to t\bar t h$ three-body decay and a phenomenological analysis of portal couplings connecting $S$ with the Higgs field.


title: The effects of short-lived radionuclides and porosity on the early thermo-mechanical evolution of planetesimals

abstract: The thermal history and internal structure of chondritic planetesimals, assembled before the giant impact phase of chaotic growth, potentially yield important implications for the final composition and evolution of terrestrial planets. These parameters critically depend on the internal balance of heating versus cooling, which is mostly determined by the presence of short-lived radionuclides (SLRs), such as aluminum-26 and iron-60, as well as the heat conductivity of the material. The heating by SLRs depends on their initial abundances, the formation time of the planetesimal and its size. It has been argued that the cooling history is determined by the porosity of the granular material, which undergoes dramatic changes via compaction processes and tends to decrease with time. In this study we assess the influence of these parameters on the thermo-mechanical evolution of young planetesimals with both 2D and 3D simulations. Using the code family I2ELVIS/I3ELVIS we have run numerous 2D and 3D numerical finite-difference fluid dynamic models with varying planetesimal radius, formation time and initial porosity. Our results indicate that powdery materials lowered the threshold for melting and convection in planetesimals, depending on the amount of SLRs present. A subset of planetesimals retained a powdery surface layer which lowered the thermal conductivity and hindered cooling. The effect of initial porosity was small, however, compared to those of planetesimal size and formation time, which dominated the thermo-mechanical evolution and were the primary factors for the onset of melting and differentiation. We comment on the implications of this work concerning the structure and evolution of these planetesimals, as well as their behavior as possible building blocks of terrestrial planets.


title: Stability mechanisms of high current transport in iron-chalcogenides superconducting films

abstract: The improvement in the fabrication techniques of iron-based superconductors have made these materials real competitors of high temperature superconductors and MgB$_2$. In particular, iron-chalcogenides have proved to be the most promising for the realization of high current carrying tapes. But their use on a large scale cannot be achieved without the understanding of the current stability mechanisms in these compounds. Indeed, we have recently observed the presence of flux flow instabilities features in Fe(Se,Te) thin films grown on CaF$_2$. Here we present the results of current-voltage characterizations at different temperatures and applied magnetic fields on Fe(Se,Te) microbridges grown on CaF$_2$. These results will be analyzed from the point of view of the most validated models with the aim to identify the nature of the flux flow instabilities features (i.e., thermal or electronic), in order to give a further advance to the high current carrying capability of iron-chalcogenide superconductors.


title: Galaxy Populations in the 26 most massive Galaxy Clusters in the South Pole Telescope SZE Survey

abstract: We present a study of the optical properties of the 26 most massive galaxy clusters selected within the SPT-SZ 2500 deg$^2$ survey. This Sunyaev-Zel'dovich effect selected sample spans a redshift range of 0.10 < z < 1.13. We measure the galaxy radial profile, the luminosity function (LF), and the halo occupation number (HON) using optical data with a typical depth of $m^*$ + 2. The stacked radial profiles are consistent with a NFW profile with a concentration of $2.84^{+0.40}_{-0.37}$ for the red sequence (RS) and $2.36^{+0.38}_{-0.35}$ for the total population. Stacking the data in multiple redshift bins shows a hint of redshift evolution in the concentration when both the total population is used, and when only RS galaxies are used (at 2.1$\sigma$ and 2.8$\sigma$, respectively). The stacked LF shows a faint end slope $\alpha = -1.06^{+0.04}_{-0.03}$ for the total and $\alpha = -0.80^{+0.04}_{-0.03}$ for the RS population. The redshift evolution of $m^*$ is found to be consistent with a passively evolving Composite Stellar Population (CSP) model. By adopting the CSP model predictions, we explore the redshift evolution of the schechter parameters $\alpha$ and $\phi^*$. We find $\alpha$ for the total population to be consistent with no evolution (0.3$\sigma$), while evidence of evolution for the red galaxies is mildly significant (1.1-2.1$\sigma$). The data show that the density $\phi^*$/E$^2$(z) decreases with redshift, in tension with the self-similar expectation at a 2.4$\sigma$ level for the total population. The measured HON-mass relation has a lower normalization than previous studies at low redshift. Finally, our data support HON redshift evolution at a 2.1$\sigma$ level, with clusters at higher redshift containing fewer galaxies per unit mass to $m^*$ + 3 than their low-z counterparts [abridged].


title: Anomalous vacuum energy and stability of a quantum liquid

abstract: We show that the vacuum (zero-point) energy of a low-temperature quantum liquid is a variable property which changes with the state of the system, in notable contrast to the static vacuum energy in solids commonly considered. We further show that this energy is inherently anomalous: it decreases with temperature and gives negative contribution to system's heat capacity. This effect operates in an equilibrium and macroscopic system, in marked contrast to small or out-of-equilibrium configurations discussed previously. We find that the negative contribution is over-compensated by the positive term from the excitation of longitudinal fluctuations and demonstrate how the overall positive heat capacity is related to the stability of a condensed phase at the microscopic level.


title: Scale-Independent Inflation and Hierarchy Generation

abstract: We discuss models involving two scalar fields coupled to classical gravity that satisfy the general criteria: (i) the theory has no mass input parameters, (ii) classical scale symmetry is broken only through $-\frac{1}{12}\varsigma \phi^2 R$ couplings where $\varsigma$ departs from the special conformal value of $1$; (iii) the Planck mass is dynamically generated by the vacuum expectations values (VEVs) of the scalars (iv) there is a stage of viable inflation associated with slow roll in the two--scalar potential; (v) the final vacuum has a small to vanishing cosmological constant and an hierarchically small ratio of the VEVs and the ratio of the scalar masses to the Planck scale. This assumes the paradigm of classical scale symmetry as a custodial symmetry of large hierarchies.


title: Bimodality of low-redshift circumgalactic O VI in non-equilibrium EAGLE zoom simulations

abstract: We introduce a series of 20 cosmological hydrodynamical simulations of Lstar (M_200 =10^11.7 - 10^12.3 Msol) and group-sized (M_200 = 10^12.7 - 10^13.3 Msol) haloes run with the model used for the EAGLE project, which additionally includes a non-equilibrium ionization and cooling module that follows 136 ions. The simulations reproduce the observed correlation, revealed by COS-Halos at z~0.2, between O VI column density at impact parameters b < 150 kpc and the specific star formation rate (sSFR=SFR/Mstar) of the central galaxy at z~0.2. We find that the column density of circumgalactic O VI is maximal in the haloes associated with Lstar galaxies, because their virial temperatures are close to the temperature at which the ionization fraction of O VI peaks (T~10^5.5 K). The higher virial temperature of group haloes (> 10^6 K) promotes oxygen to higher ionization states, suppressing the O VI column density. The observed NO VI-sSFR correlation therefore does not imply a causal link, but reflects the changing characteristic ionization state of oxygen as halo mass is increased. In spite of the mass-dependence of the oxygen ionization state, the most abundant circumgalactic oxygen ion in both Lstar and group haloes is O VII; O VI accounts for only 0.1% of the oxygen in group haloes and 0.9-1.3% with Lstar haloes. Nonetheless, the metals traced by O VI absorbers represent a fossil record of the feedback history of galaxies over a Hubble time; their characteristic epoch of ejection corresponds to z > 1 and much of the ejected metal mass resides beyond the virial radius of galaxies. For both Lstar and group galaxies, more of the oxygen produced and released by stars resides in the circumgalactic medium (within twice the virial radius) than in the stars and ISM of the galaxy.


title: A Network Approach to Unravel Correlated Ion Pair Dynamics in Protic Ionic Liquids. The Case of Triethylammonium Nitrate

abstract: The intermolecular interactions in the title compound are investigated using self-consistent charge density functional based tight binding molecular dynamics. Emphasis is put on the analysis of correlated motions of ion pairs using ideas of network theory. At equilibrium such correlations are not very pronounced on average. However, there exist sizeable local correlations for cases where two cations share the same anion via two NHO-hydrogen bonds. The effect of an external perturbation, which artificially introduces a sudden local heating of an NH-bond, is investigated using nonequilibrium molecular dynamics. Here, it is found that the average N-H bond vibrational relaxation time is about 5.3 ~ps. This energy redistribution is rather nonspecific with respect to the ion pairs and does not lead to long-range correlations spreading from the initially excited ion pair.


title: Green$'$s function approach to edge states in transition metal dichalcogenides

abstract: The semiconducting two-dimensional transition metal dichalcogenides MX$_{2}$ show an abundance of one-dimensional metallic edges and grain boundaries. Standard techniques for calculating edge states typically model nanoribbons, and require the use of supercells. In this paper we formulate a Green$'$s function technique for calculating edge states of (semi-)infinite two-dimensional systems with a single well-defined edge or grain boundary. We express Green$'$s functions in terms of Bloch matrices, constructed from the solutions of a quadratic eigenvalue equation. The technique can be applied to any localized basis representation of the Hamiltonian. Here we use it to calculate edge states of MX$_{2}$ monolayers by means of tight-binding models. Besides the basic zigzag and armchair edges, we study edges with a more general orientation, structurally modifed edges, and grain boundaries. A simple three-band model captures an important part of the edge electronic structures. An eleven-band model comprising all valence orbitals of the M and X atoms, is required to obtain all edge states with energies in the MX$_{2}$ band gap. Here states of odd symmetry with respect to a mirror plane through the layer of M atoms have a dangling-bond character, and tend to pin the Fermi level.


title: SO(d,1)-invariant Yang-Baxter operators and the dS/CFT correspondence

abstract: We propose a model for the dS/CFT correspondence. The model is constructed in terms of a "Yang-Baxter operator" $R$ for unitary representations of the deSitter group $SO(d,1)$. This $R$-operator is shown to satisfy the Yang-Baxter equation, unitarity, as well as certain analyticity relations, including in particular a crossing symmetry. With the aid of this operator we construct: a) A chiral (light-ray) conformal quantum field theory whose internal degrees of freedom transform under the given unitary representation of $SO(d,1)$. By analogy with the $O(N)$ non-linear sigma model, this chiral CFT can be viewed as propagating in a deSitter spacetime. b) A (non-unitary) Euclidean conformal quantum field theory on ${\mathbb R}^{d-1}$, where $SO(d,1)$ now acts by conformal transformations in (Euclidean) spacetime. These two theories can be viewed as dual to each other if we interpret ${\mathbb R}^{d-1}$ as conformal infinity of deSitter spacetime. Our constructions use semi-local generator fields defined in terms of $R$ and abstract methods from operator algebras.


title: Deflation as a Method of Variance Reduction for Estimating the Trace of a Matrix Inverse

abstract: Many fields require computing the trace of the inverse of a large, sparse matrix. The typical method used for such computations is the Hutchinson method which is a Monte Carlo (MC) averaging over matrix quadratures. To improve its convergence, several variance reductions techniques have been proposed. In this paper, we study the effects of deflating the near null singular value space. We make two main contributions. First, we analyze the variance of the Hutchinson method as a function of the deflated singular values and vectors. Although this provides good intuition in general, by assuming additionally that the singular vectors are random unitary matrices, we arrive at concise formulas for the deflated variance that include only the variance and mean of the singular values. We make the remarkable observation that deflation may increase variance for Hermitian matrices but not for non-Hermitian ones. This is a rare, if not unique, property where non-Hermitian matrices outperform Hermitian ones. The theory can be used as a model for predicting the benefits of deflation. Second, we use deflation in the context of a large scale application of "disconnected diagrams" in Lattice QCD. On lattices, Hierarchical Probing (HP) has previously provided an order of magnitude of variance reduction over MC by removing "error" from neighboring nodes of increasing distance in the lattice. Although deflation used directly on MC yields a limited improvement of 30% in our problem, when combined with HP they reduce variance by a factor of over 60 compared to MC. For this, we pre-computated 1000 smallest singular values of an ill-conditioned matrix of size 25 million. Using PRIMME and a domain-specific Algebraic Multigrid preconditioner, we perform one of the largest eigenvalue computations in Lattice QCD at a fraction of the cost of our trace computation.


title: Breadth versus depth: Interactions that stabilize particle assemblies to changes in density or temperature

abstract: We use inverse methods of statistical mechanics to explore trade-offs associated with designing interactions to stabilize self-assembled structures against changes in density or temperature. Specifically, we find isotropic,convex-repulsive pair potentials that maximize the density range for which a two-dimensional square lattice is the stable ground state subject to a constraint on the chemical potential advantage it exhibits over competing structures (i.e., 'depth' of the associated minimum on the chemical potential hypersurface). We formulate the design problem as a non-linear program, which we solve numerically. This allows us to efficiently find optimized interactions for a wide range of possible chemical potential constraints. We find that assemblies designed to exhibit a large chemical potential advantage at a specified density have a smaller overall range of densities for which they are stable. This trend can be understood by considering the separation-dependent features of the pair potential and its gradient required to enhance the stability of the target structure relative to competitors. Using molecular dynamics simulations, we further show that potentials designed with larger chemical potential advantages exhibit higher melting temperatures.


title: The LOFAR Solar Imaging Pipeline and the LOFAR Solar Data Center

abstract: LOFAR is a new and sensitive radio interferometer that can be used for dynamic high-resolution imaging spectroscopy at low radio frequencies from 10 to 90 and 110 to 250 MHz. Here we describe its usage for observations of the Sun and in particular of solar radio bursts. We also describe the processing, archiving and accessing of solar LOFAR data, which is accomplished via the LOFAR Solar Imaging Pipeline and the LOFAR Solar Data Center.


title: Flat line bundles and the Cappell-Miller torsion in Arakelov geometry

abstract: In this paper, we extend Deligne's functorial Riemann-Roch isomorphism for hermitian holomorphic line bundles on Riemann surfaces to the case of flat, not necessarily unitary connections. The Quillen metric and star-product of Gillet-Soule are replaced with complex valued logarithms. On the determinant of cohomology side, the idea goes back to Fay's holomorphic extension of determinants of Dolbeault laplacians, and it is shown here to be equivalent to the holomorphic Cappell-Miller torsion. On the Deligne pairing side, the logarithm is a refinement of the intersection connections considered in previous work. The construction naturally leads to an Arakelov theory for flat line bundles on arithmetic surfaces and produces arithmetic intersection numbers valued in ${\mathbb C}/\pi i {\mathbb Z}$. In this context we prove an arithmetic Riemann-Roch theorem. This realizes a program proposed by Cappell-Miller to show that the holomorphic torsion exhibits properties similar to those of the Quillen metric proved by Bismut, Gillet and Soule. Finally, we give examples that clarify the kind of invariants that the formalism captures; namely, periods of differential forms.


title: The VIMOS Ultra Deep Survey: Lyman Alpha Emission and Stellar Populations of Star-Forming Galaxies at 2<z<6

abstract: The extensive ground-based spectroscopy campaign from the VIMOS Ultra-Deep Survey (VUDS), and the deep multi-wavelength photometry in three very well observed extragalactic fields (ECDFS, COSMOS, VVDS), allow us to investigate physical properties of a large sample (~4000 galaxies) of spectroscopically confirmed faint (i_{AB}<~25 mag) star-forming galaxies, with and without Lyman alpha in emission, at z~2-6. The fraction of Lyman alpha emitters (LAEs; equivalent width (EW)=>20A) increases from ~10% at z~2 to ~40% at z~5-6, which is consistent with previous studies that employ higher Lyman alpha EW cut. This increase in the LAE fraction could be, in part, due to a decrease in the dust content of galaxies as redshift increases. When we compare best-fit SED estimated stellar parameters for LAEs and non-LAEs, we find that E(B-V) is smaller for LAEs at all redshifts and the difference in the median E(B-V) between LAEs and non-LAEs increases as redshift increases, from 0.05 at z~2 to 0.1 at z~3.5 to 0.2 at z~5-6. For the luminosities probed here (~L*), we find that star formation rates (SFRs) and stellar masses of galaxies, with and without Lyman alpha in emission, show small differences such that, LAEs have lower SFRs and stellar masses compared to non-LAEs. This result could be a direct consequence of the sample selection. Our sample of LAEs are selected based on their continuum magnitudes and they probe higher continuum luminosities compared to narrow-band/emission line selected LAEs. Based on our results, it is important to note that all LAEs are not universally similar and their properties are strongly dependent on the sample selection, and/or continuum luminosities.


title: A class of spherical, truncated, anisotropic models for application to globular clusters

abstract: Recently, a class of non-truncated radially-anisotropic models (the so-called $f^{(\nu)}$-models), originally constructed in the context of violent relaxation and modeling of elliptical galaxies, has been found to possess interesting qualities in relation to observed and simulated globular clusters. In view of new applications to globular clusters, we improve this class of models along two directions. To make them more suitable for the description of small stellar systems hosted by galaxies, we introduce a 'tidal' truncation (by means of a procedure that guarantees full continuity of the distribution function). The new $f_T^{(\nu)}$-models are shown to provide a better fit to the observed photometric and spectroscopic profiles for a sample of 13 globular clusters studied earlier by means of non-truncated models; interestingly, the best-fit models also perform better with respect to the radial-orbit instability. Then we design a flexible but simple two-component family of truncated models, to study the separate issues of mass segregation and of multiple populations. We do not aim at a fully realistic description of globular clusters, to compete with the description currently obtained by means of dedicated simulations. The goal here is to try to identify the simplest models, that is, those with the smallest number of free parameters, but still able to provide a reasonable description for clusters that are evidently beyond the reach of one-component models: with this tool we aim at identifying the key factors that characterize mass segregation or the presence of multiple populations. To reduce the relevant parameter space, we formulate a few physical arguments (based on recent observations and simulations). A first application to two well-studied globular clusters is briefly described and discussed.


title: Effective Field Theory of Dark Matter: a Global Analysis

abstract: We present global fits of an effective field theory description of real, and complex scalar dark matter candidates. We simultaneously take into account all possible dimension 6 operators consisting of dark matter bilinears and gauge invariant combinations of quark and gluon fields. We derive constraints on the free model parameters for both the real (five parameters) and complex (seven) scalar dark matter models obtained by combining Planck data on the cosmic microwave background, direct detection limits from LUX, and indirect detection limits from the Fermi Large Area Telescope. We find that for real scalars indirect dark matter searches disfavour a dark matter particle mass below 100 GeV. For the complex scalar dark matter particle current data have a limited impact due to the presence of operators that lead to p-wave annihilation, and also do not contribute to the spin-independent scattering cross- section. Although current data are not informative enough to strongly constrain the theory parameter space, we demonstrate the power of our formalism to reconstruct the theoretical parameters compatible with an actual dark matter detection, by assuming that the excess of gamma rays observed by the Fermi Large Area Telescope towards the Galactic centre is entirely due to dark matter annihilations. Please note that the excess can very well be due to astrophysical sources such as millisecond pulsars. We find that scalar dark matter interacting via effective field theory operators can in principle explain the Galactic centre excess, but that such interpretation is in strong tension with the non-detection of gamma rays from dwarf galaxies in the real scalar case. In the complex scalar case there is enough freedom to relieve the tension.


title: Diffeomorphism groups of compact convex sets

abstract: For a compact convex subset K with non-empty interior in a finite-dimensional vector space, let G be the group of all smooth diffeomorphisms of K which fix the boundary of K pointwise. We show that G is a C^0-regular infinite-dimensional Lie group. As a byproduct, we obtain results concerning solutions to ordinary differential equations on compact convex sets.