title: Vertex adjacencies in the set covering polyhedron

abstract: We describe the adjacency of vertices of the (unbounded version of the) set covering polyhedron, in a similar way to the description given by Chvatal for the stable set polytope. We find a sufficient condition for adjacency, and characterize it with similar conditions in the case where the underlying matrix is row circular. We apply our findings to show a new infinite family of minimally nonideal matrices.


title: Noncollinear magnetic phases and edge states in graphene quantum Hall bars

abstract: Application of a perpendicular magnetic field to charge neutral graphene is expected to result in a variety of broken symmetry phases, including antiferromagnetic, canted and ferromagnetic. All these phases open a gap in bulk but have very different edge states and non-collinear spin order, recently confirmed experimentally. Here we provide an integrated description of both edge and bulk for the various magnetic phases of graphene Hall bars making use of a non-collinear mean field Hubbard model. Our calculations show that, at the edges, the three types of magnetic order are either enhanced (zigzag) or suppressed (armchair). Interestingly, we find that preformed local moments in zigzag edges interact with the quantum Spin Hall like edge states of the ferromagnetic phase and can induce back-scattering.


title: Uniform-in-bandwidth consistency for nonparametric estimation of divergence measures

abstract: We propose nonparametric estimation of divergence measures between continuous distributions. Our approach is based on a plug-in kernel- type estimators of density functions. We give the uniform in bandwidth consistency for the proposal estimators. As a consequence, their asymp- totic 100% confidence intervals are also provided.


title: A brief history of long memory: Hurst, Mandelbrot and the road to ARFIMA

abstract: Long memory plays an important role in many fields by determining the behaviour and predictability of systems; for instance, climate, hydrology, finance, networks and DNA sequencing. In particular, it is important to test if a process is exhibiting long memory since that impacts the accuracy and confidence with which one may predict future events on the basis of a small amount of historical data. A major force in the development and study of long memory was the late Benoit B. Mandelbrot. Here we discuss the original motivation of the development of long memory and Mandelbrot's influence on this fascinating field. We will also elucidate the sometimes contrasting approaches to long memory in different scientific communities


title: Anisotropic shear viscosity of a strongly coupled non-Abelian plasma from magnetic branes

abstract: Recent estimates for the electromagnetic fields produced in the early stages of non-central ultra-relativistic heavy ion collisions indicate the presence of magnetic fields $B\sim \mathcal{O}(0.1-15\,m_\pi^2)$, where $m_\pi$ is the pion mass. It is then of special interest to study the effects of strong (Abelian) magnetic fields on the transport coefficients of strongly coupled non-Abelian plasmas, such as the quark-gluon plasma formed in heavy ion collisions. In this work we study the anisotropy in the shear viscosity induced by an external magnetic field in a strongly coupled $\mathcal{N} = 4$ SYM plasma. Due to the spatial anisotropy created by the magnetic field, the most general viscosity tensor of a magnetized plasma has 5 shear viscosity coefficients and 2 bulk viscosities. We use the holographic correspondence to evaluate two of the shear viscosities, $\eta_{\perp} \equiv \eta_{xyxy}$ (perpendicular to the magnetic field) and $\eta_{\parallel} \equiv \eta_{xzxz}=\eta_{yzyz}$ (parallel to the field). When $B\neq 0$ the shear viscosity perpendicular to the field saturates the viscosity bound $\eta_{\perp}/s = 1/(4\pi)$ while in the direction parallel to the field the bound is violated since $\eta_{\parallel}/s < 1/(4\pi)$. However, the violation of the bound in the case of strongly coupled SYM is minimal even for the largest value of $B$ that can be reached in heavy ion collisions.


title: Stationary Mixing Bandits

abstract: We study the bandit problem where arms are associated with stationary phi-mixing processes and where rewards are therefore dependent: the question that arises from this setting is that of recovering some independence by ignoring the value of some rewards. As we shall see, the bandit problem we tackle requires us to address the exploration/exploitation/independence trade-off. To do so, we provide a UCB strategy together with a general regret analysis for the case where the size of the independence blocks (the ignored rewards) is fixed and we go a step beyond by providing an algorithm that is able to compute the size of the independence blocks from the data. Finally, we give an analysis of our bandit problem in the restless case, i.e., in the situation where the time counters for all mixing processes simultaneously evolve.


title: Analytic Continuation of Black Hole Entropy in Loop Quantum Gravity

abstract: We define the analytic continuation of the number of black hole microstates in Loop Quantum Gravity to complex values of the Barbero-Immirzi parameter $\gamma$. This construction deeply relies on the link between black holes and Chern-Simons theory. Technically, the key point consists in writing the number of microstates as an integral in the complex plane of a holomorphic function, and to make use of complex analysis techniques to perform the analytic continuation. Then, we study the thermodynamical properties of the corresponding system (the black hole is viewed as a gas of indistinguishable punctures) in the framework of the grand canonical ensemble where the energy is defined \'a la Frodden-Gosh-Perez from the point of view of an observer located close to the horizon. The semi-classical limit occurs at the Unruh temperature $T_U$ associated to this local observer. When $\gamma=\pm i$, the entropy reproduces at the semi-classical limit the area law with quantum corrections. Furthermore, the quantum corrections are logarithmic provided that the chemical potential is fixed to the simple value $\mu=2T_U$.


title: Charge Storage in Cation Incorporated {\alpha}-MnO2

abstract: Electrochemical supercapacitors utilizing {\alpha}-MnO2 offer the possibility of both high power density and high energy density. Unfortunately, the mechanism of electrochemical charge storage in {\alpha}-MnO2 and the effect of operating conditions on the charge storage mechanism are generally not well understood. Here, we present the first detailed charge storage mechanism of {\alpha}-MnO2 and explain the capacity differences between {\alpha}- and {\beta}-MnO2 using a combined theoretical electrochemical and band structure analysis. We identify the importance of the band gap, work function, the point of zero charge, and the tunnel sizes of the electrode material, as well as the pH and stability window of the electrolyte in determining the viability of a given electrode material. The high capacity of {\alpha}-MnO2 results from cation induced charge-switching states in the band gap that overlap with the scanned potential allowed by the electrolyte. The charge-switching states originate from interstitial and substitutional cations (H+, Li+, Na+, and K+) incorporated into the material. Interstitial cations are found to induce charge-switching states by stabilizing Mn-O antibonding orbitals from the conduction band. Substitutional cations interact with O[2p] dangling bonds that are destabilized from the valence band by Mn vacancies to induce charge-switching states. We calculate the equilibrium electrochemical potentials at which these states are reduced and predict the effect of the electrochemical operating conditions on their contribution to charge storage. The mechanism and theoretical approach we report is general and can be used to computationally screen new materials for improved charge storage via ion incorporation.


title: Note on the self-duality of gauge fields in topologically nontrivial spacetime

abstract: We show the derivation of the self-duality relation of abelian higher-form gauge field strength in the topologically nontrivial spacetime background. The so-called Pasti-Sorokin-Tonin action for the self-dual abelian gauge field assumes that the spacetime topology is trivial to derive the self-duality relation using a gauge transformation of the action. In this paper we find a new gauge transformation of the same theory and show that this new gauge transformation enables us to derive the self-duality relation even in the topologically nontrivial spacetime background.


title: Monte Carlo simulation of MLC-shaped TrueBeam electron fields benchmarked against measurement

abstract: Modulated electron radiotherapy (MERT) and combined, modulated photon/electron radiotherapy (MPERT) have received increased research attention, having shown capacity for reduced low dose exposure to healthy tissue and comparable, if not improved, target coverage for a number of treatment sites. Accurate dose calculation tools are necessary for clinical treatment planning, and Monte Carlo (MC) is the gold standard for electron field simulation. With many clinics replacing older accelerators, MC source models of the new machines are needed for continued development, however, Varian has kept internal schematics of the TrueBeam confidential and electron phase-space sources have not been made available. TrueBeam electron fields are not substantially different from those generated by the Clinac 21EX, so we have modified the internal schematics of the Clinac 21EX to simulate TrueBeam electrons. BEAMnrc/DOSXYZnrc were used to simulate 5x5 and 20x20 cm$^2$ electron fields with MLC-shaped apertures. Secondary collimating jaws were set 0.5 cm beyond the MLC periphery and to 40x40 cm$^2$. Complete accelerator models agreed with diode measurement within 2%/2 mm at 12 and 20 MeV, and within 3%/3 mm at 6 MeV. Comparisons of measured depth and profile data showed dose-dependencies on jaw position; simulated energy fluences scored just above the phantom showed that, for small apertures, dose dependencies on jaw position are dominated by changes in the in-field energy fluence (not scattered by jaws or MLC) at 6 MeV, while at 20 MeV dose dependencies are dominated by the scattered component. Our models reproduce these jaw position dependencies, which is an asset as there is no consensus on the optimal position for jaws in modulated electron field delivery. Good agreement between simulation and measurement and flexibility in jaw position make our models appropriate for use in MERT and MPERT planning.


title: Initial Performance of the NEOWISE Reactivation Mission

abstract: NASA's Wide-field Infrared Survey Explorer (WISE) spacecraft has been brought out of hibernation and has resumed surveying the sky at 3.4 and 4.6 um. The scientific objectives of the NEOWISE reactivation mission are to detect, track, and characterize near-Earth asteroids and comets. The search for minor planets resumed on December 23, 2013, and the first new near-Earth object (NEO) was discovered six days later. As an infrared survey, NEOWISE detects asteroids based on their thermal emission and is equally sensitive to high and low albedo objects; consequently, NEOWISE-discovered NEOs tend to be large and dark. Over the course of its three-year mission, NEOWISE will determine radiometrically-derived diameters and albedos for approximately 2000 NEOs and tens of thousands of Main Belt asteroids. The 32 months of hibernation have had no significant effect on the mission's performance. Image quality, sensitivity, photometric and astrometric accuracy, completeness, and the rate of minor planet detections are all essentially unchanged from the prime mission's post-cryogenic phase.


title: The Quadratic Approximation for Quintessence with Arbitrary Initial Conditions

abstract: We examine quintessence models for dark energy in which the scalar field, $\phi$, evolves near the vicinity of a local maximum or minimum in the potential $V(\phi)$, so that $V(\phi)$ be approximated by a quadratic function of $\phi$ with no linear term. We generalize previous studies of this type by allowing the initial value of $d \phi/dt$ to be nonzero. We derive an analytic approximation for $w(a)$ and show that it is in excellent agreement with numerical simulations for a variety of scalar field potentials having local minima or maxima. We derive an upper bound on the present-day value of $w$ as a function of the other model parameters and present representative limits on these models from observational data. This work represents a final generalization of previous studies using linear or quadratic approximations for $V(\phi)$.


title: Updated cosmic-ray and radio constraints on light dark matter: Implications for the GeV gamma-ray excess at the Galactic center

abstract: The apparent gamma-ray excess in the Galactic center region and inner Galaxy has attracted considerable interest, notably because both its spectrum and radial distribution are consistent with an interpretation in terms of annihilating dark matter particles with a mass of about 10-40 GeV. We confront such an interpretation with an updated compilation of various indirect dark matter detection bounds, which we adapt to the specific form required by the observed signal. We find that cosmic-ray positron data strongly rule out dark matter annihilating to light leptons, or 'democratically' to all leptons, as an explanation of the signal. Cosmic-ray antiprotons, for which we present independent and significantly improved limits with respect to previous estimates, are already in considerable tension with DM annihilation to any combination of quark final states; the first set of AMS-02 data will thus be able to rule out or confirm the DM hypothesis with high confidence. For reasonable assumptions about the magnetic field in the Galactic center region, radio observations independently put very severe constraints on a DM interpretation of the excess, in particular for all leptonic annihilation channels.


title: Nonsmooth frameworks for an extended Budyko model

abstract: In latitude-dependent energy balance models, ice-free and ice-covered conditions form physical boundaries of the system. With carbon dioxide treated as a bifurcation parameter, the resulting bifurcation diagram is nonsmooth with curves of equilibria and boundaries forming corners at points of intersection. Over long time scales, atmospheric carbon dioxide varies dynamically and the nonsmooth diagram becomes a set of quasi-equilibria. {However, when introducing carbon dynamics, care must be taken with the physical boundaries and appropriate boundary motion specified. In this article, we extend an energy balance model to include slowly varying carbon dioxide and develop a nonsmooth framework based on physically relevant boundary dynamics. Within this framework, we prove existence and uniqueness of solutions, as well as invariance of the region of phase space bounded by ice-free and ice-covered states.


title: A note on the unit distance problem for planar configurations with Q-independent direction set

abstract: Let $T(n)$ denote the maximum number of unit distances that a set of $n$ points in the Euclidean plane $\mathbb{R}^2$ can determine with the additional condition that the distinct unit length directions determined by the configuration must be $\mathbb{Q}$-independent. This is related to the Erdos unit distance problem but with a simplifying additional assumption on the direction set which holds "generically". We show that $T(n+1)-T(n)$ is the Hamming weight of $n$, i.e., the number of nonzero binary coefficients in the binary expansion of $n$, and find a formula for $T(n)$ explicitly. In particular $T(n)$ is $\Theta(n log(n))$. Furthermore we describe a process to construct a set of $n$ points in the plane with $\mathbb{Q}$-independent unit length direction set which achieves exactly $T(n)$ unit distances. In the process of doing this, we show $T(n)$ is also the same as the maximum number of edges a subset of vertices of size $n$ determines in either the countably infinite lattice $\mathbb{Z}^{\infty}$ or the infinite hypercube graph $\{0,1\}^{\infty}$. The problem of determining T(n) can be viewed as either a type of packing or isoperimetric problem.


title: Categorical Probability Theory

abstract: We present a categorical viewpoint of probability measures by showing that a probability measure can be viewed as a weakly averaging affine measurable functional taking values in the unit interval which preserves limits. The probability measures on a space are the elements of a submonad of a double dualization monad on the category of measurable spaces into the unit interval, and this monad is naturally isomorphic to the Giry monad. We show this submonad is the codensity monad of a functor from the category of convex spaces to the category of measurable spaces. A theorem proving the integral operator acting on the space of measurable functions and the space of probability measures on the domain space of those functions is given using the strong monad structure of the Giry monad.


title: Robust estimation of multivariate location and scatter in the presence of cellwise and casewise contamination

abstract: Multivariate location and scatter matrix estimation is a cornerstone in multivariate data analysis. We consider this problem when the data may contain independent cellwise and casewise outliers. Flat data sets with a large number of variables and a relatively small number of cases are common place in modern statistical applications. In these cases global down-weighting of an entire case, as performed by traditional robust procedures, may lead to poor results. We highlight the need for a new generation of robust estimators that can efficiently deal with cellwise outliers and at the same time show good performance under casewise outliers.


title: Displacive-type ferroelectricity from magnetic correlations within spin-chain

abstract: Observation of ferroelectricity among non-d0 systems, which was believed for a long time an unrealistic concept, led to various proposals for the mechanisms to explain the same (i.e. magnetically induced ferroelectricity) during last decade. Here, we provide support for ferroelectricity of a displacive-type possibly involving magnetic ions due to short-range magnetic correlations within a spin-chain, through the demonstration of magnetoelectric coupling in a Haldane spin-chain compound Er2BaNiO5 well above its N\'eel temperature of (TN=) 32K. There is a distinct evidence for electric polarization setting in near 60 K around which there is an evidence for short-range magnetic correlations from other experimental methods. Raman studies also establish a softening of phonon modes in the same temperature (T) range and T-dependent x-ray diffraction (XRD) patterns also reveal lattice parameters anomalies. Density-functional theory based calculations establish a displacive component (similar to d0-ness) as the root-cause of ferroelectricity from (magnetic) NiO6 chain, thereby offering a new route to search for similar materials near room temperature to enable applications.


title: Mutations and short geodesics in hyperbolic 3-manifolds

abstract: In this paper, we explicitly construct large classes of incommensurable hyperbolic knot complements with the same volume and the same initial (complex) length spectrum. Furthermore, we show that these knot complements are the only knot complements in their respective commensurabiltiy classes by analyzing their cusp shapes. The knot complements in each class differ by a topological cut-and-paste operation known as mutation. Ruberman has shown that mutations of hyperelliptic surfaces inside hyperbolic 3-manifolds preserve volume. Here, we provide geometric and topological conditions under which such mutations also preserve the initial (complex) length spectrum. This work requires us to analyze when least area surfaces could intersect short geodesics in a hyperbolic 3-manifold.


title: Constraining the neutrino magnetic dipole moment from white dwarf pulsations

abstract: Pulsating white dwarf stars can be used as astrophysical laboratories to constrain the properties of weakly interacting particles. Comparing the cooling rates of these stars with the expected values from theoretical models allows us to search for additional sources of cooling due to the emission of axions, neutralinos, or neutrinos with magnetic dipole moment. In this work, we derive an upper bound to the neutrino magnetic dipole moment using an estimate of the rate of period change of the pulsating DB white dwarf star PG 1351+489. By comparing the theoretical rate of change of period expected for this star with the rate of change of period with time of PG 1351+489, we assess the possible existence of additional cooling by neutrinos with magnetic dipole moment. Our models suggest the existence of some additional cooling in this pulsating DB white dwarf, consistent with a non-zero magnetic dipole moment. Our upper limit for the neutrino magnetic dipole moment is somewhat less restrictive than, but still compatible with, other limits inferred from the white dwarf luminosity function or from the color-magnitude diagram of the Globular cluster M5. Further improvements of the measurement of the rate of period change of the dominant pulsation mode of PG 1351+489 will be necessary to confirm our bound.


title: Refinement Calculus of Reactive Systems

abstract: Refinement calculus is a powerful and expressive tool for reasoning about sequential programs in a compositional manner. In this paper we present an extension of refinement calculus for reactive systems. Refinement calculus is based on monotonic predicate transformers, which transform sets of post-states into sets of pre-states. To model reactive systems, we introduce monotonic property transformers, which transform sets of output traces into sets of input traces. We show how to model in this semantics refinement, sequential composition, demonic choice, and other semantic operations on reactive systems. We use primarily higher order logic to express our results, but we also show how property transformers can be defined using other formalisms more amenable to automation, such as linear temporal logic (suitable for specifications) and symbolic transition systems (suitable for implementations). Finally, we show how this framework generalizes previous work on relational interfaces so as to be able to express systems with infinite behaviors and liveness properties.


title: Generating non-classical states from spin coherent states via interaction with ancillary spins

abstract: The generation of non-classical states of large quantum systems has attracted much interest from a foundational perspective, but also because of the significant potential of such states in emerging quantum technologies. In this paper we consider the possibility of generating non-classical states of a system of spins by interaction with an ancillary system, starting from an easily prepared initial state . We extend previous results for an ancillary system comprising a single spin to bigger ancillary systems and the interaction strength is enhanced by a factor of the number of ancillary spins. Depending on initial conditions, we find -- by a combination of approximation and numerics -- that the system of spins can evolve to spin cat states, spin squeezed states or to multiple cat states. We also discuss some candidate systems for implementation of the Hamiltonian necessary to generate these non-classical states.


title: Preemptive Thread Block Scheduling with Online Structural Runtime Prediction for Concurrent GPGPU Kernels

abstract: Recent NVIDIA Graphics Processing Units (GPUs) can execute multiple kernels concurrently. On these GPUs, the thread block scheduler (TBS) uses the FIFO policy to schedule their thread blocks. We show that FIFO leaves performance to chance, resulting in significant loss of performance and fairness. To improve performance and fairness, we propose use of the preemptive Shortest Remaining Time First (SRTF) policy instead. Although SRTF requires an estimate of runtime of GPU kernels, we show that such an estimate of the runtime can be easily obtained using online profiling and exploiting a simple observation on GPU kernels' grid structure. Specifically, we propose a novel Structural Runtime Predictor. Using a simple Staircase model of GPU kernel execution, we show that the runtime of a kernel can be predicted by profiling only the first few thread blocks. We evaluate an online predictor based on this model on benchmarks from ERCBench, and find that it can estimate the actual runtime reasonably well after the execution of only a single thread block. Next, we design a thread block scheduler that is both concurrent kernel-aware and uses this predictor. We implement the SRTF policy and evaluate it on two-program workloads from ERCBench. SRTF improves STP by 1.18x and ANTT by 2.25x over FIFO. When compared to MPMax, a state-of-the-art resource allocation policy for concurrent kernels, SRTF improves STP by 1.16x and ANTT by 1.3x. To improve fairness, we also propose SRTF/Adaptive which controls resource usage of concurrently executing kernels to maximize fairness. SRTF/Adaptive improves STP by 1.12x, ANTT by 2.23x and Fairness by 2.95x compared to FIFO. Overall, our implementation of SRTF achieves system throughput to within 12.64% of Shortest Job First (SJF, an oracle optimal scheduling policy), bridging 49% of the gap between FIFO and SJF.


title: Exact fit of simple finite mixture models

abstract: How to forecast next year's portfolio-wide credit default rate based on last year's default observations and the current score distribution? A classical approach to this problem consists of fitting a mixture of the conditional score distributions observed last year to the current score distribution. This is a special (simple) case of a finite mixture model where the mixture components are fixed and only the weights of the components are estimated. The optimum weights provide a forecast of next year's portfolio-wide default rate. We point out that the maximum-likelihood (ML) approach to fitting the mixture distribution not only gives an optimum but even an exact fit if we allow the mixture components to vary but keep their density ratio fix. From this observation we can conclude that the standard default rate forecast based on last year's conditional default rates will always be located between last year's portfolio-wide default rate and the ML forecast for next year. As an application example, then cost quantification is discussed. We also discuss how the mixture model based estimation methods can be used to forecast total loss. This involves the reinterpretation of an individual classification problem as a collective quantification problem.


title: Boundary Harnack estimates in slit domains and applications to thin free boundary problems

abstract: We provide a higher order boundary Harnack inequality for harmonic functions in slit domains. As a corollary we obtain the $C^\infty$ regularity of the free boundary in the Signorini problem near non-degenerate points.


title: Stopping the SuperSpreader Epidemic, Part III: Prediction

abstract: In two previous papers, I introduced SuperSpreader (SS) epidemic models, offered some theoretical discussion of prevention issues, and fitted some models to data derived from published accounts of the ongoing MERS epidemic (concluding that a pandemic is likely). Continuing on this theme, here I discuss prediction: whether, in a disease outbreak driven by superspreader events, a rigorous decision point---meaning a declaration that a pandemic is imminent---can be defined. I show that all sources of prediction bias contribute to generating false negatives (i.e., discounting the chance of a pandemic when it is looming or has already started). Nevertheless, the statistical difficulties can be overcome by improved data gathering and use of known techniques that decrease bias. One peculiarity of the SS epidemic is that the prediction can sometimes be made long before the actual pandemic onset, generating lead time to alert the medical community and the public. Thus modeling is useful to overcome a false sense of security arising from the long "kindling times" characteristic of SS epidemics and certain political/psychological factors, as well as improve the public health response.


title: The moduli scheme of affine spherical varieties with a free weight monoid

abstract: We study Alexeev and Brion's moduli scheme $M_\Gamma$ of affine spherical varieties with weight monoid $\Gamma$ under the assumption that $\Gamma$ is free. We describe the tangent space to $M_\Gamma$ at its `most degenerate point' in terms of the combinatorial invariants of spherical varieties and deduce that the irreducible components of $M_\Gamma$, equipped with their reduced induced scheme structure, are affine spaces.


title: Global analysis of a model for capillary water waves in 2D

abstract: In this paper we prove a global regularity result for a quadratic quasilinear model associated to the water waves system with surface tension and no gravity in dimension two (the capillary waves system). The model we consider here retains most of the difficulties of the full capillary water waves system, including the delicate time-resonance structure and modified scattering. It is slightly simpler, however, at the technical level and our goal here is to present our method in this simplified situation. The full system is the subject of a forthcoming paper.


title: Geometry of the Infalling Causal Patch

abstract: The firewall paradox states that an observer falling into an old black hole must see a violation of unitarity, locality, or the equivalence principle. Motivated by this remarkable conflict, we analyze the causal structure of black hole spacetimes in order to determine whether all the necessary ingredients for the paradox fit within a single observer's causal patch. We particularly focus on the question of whether the interior partner modes of the outgoing Hawking quanta can, in principle, be measured by an infalling observer. Since the relevant modes are spread over the entire sphere, we answer a simple geometrical question: can any observer see an entire sphere behind the horizon? We find that for all static black holes in 3+1 and higher dimensions, with any value of the cosmological constant, no single observer can see both the early Hawking radiation and the interior modes. We present a detailed description of the causal patch geometry of the Schwarzschild black hole in 3+1 dimensions, where an infalling observer comes closest to being able to measure the relevant modes.


title: The rate of increase for recursion with quadratic non-linearity

abstract: We investigate in this short report the rate of increase of positive numerical recursion with quadratic non-linearity. More exactly, we intent to calculate the logarithmic index of its increasing. We present also the possible application in the theory of the Navier-Stokes equations.


title: Cognitive Surveillance: Why does it never appear among the AVSS Conferences topics?

abstract: Video Surveillance is a fast evolving field of research and development (R&D) driven by the urgent need for public security and safety (due to the growing threats of terrorism, vandalism, and anti-social behavior). Traditionally, surveillance systems are comprised of two components - video cameras distributed over the guarded area and human observer watching and analyzing the incoming video. Explosive growth of installed cameras and limited human operator's ability to process the delivered video content raise an urgent demand for developing surveillance systems with human like cognitive capabilities, that is - Cognitive surveillance systems. The growing interest in this issue is testified by the tens of workshops, symposiums and conferences held over the world each year. The IEEE International Conference on Advanced Video and Signal-Based Surveillance (AVSS) is certainly one of them. However, for unknown reasons, the term Cognitive Surveillance does never appear among its topics. As to me, the explanation for this is simple - the complexity and the indefinable nature of the term "Cognition". In this paper, I am trying to resolve the problem providing a novel definition of cognition equally suitable for biological as well as technological applications. I hope my humble efforts will be helpful.


title: Hybrid Epidemics - A Case Study on Computer Worm Conficker

abstract: Conficker is a computer worm that erupted on the Internet in 2008. It is unique in combining three different spreading strategies: local probing, neighbourhood probing, and global probing. We propose a mathematical model that combines three modes of spreading, local, neighbourhood and global to capture the worm's spreading behaviour. The parameters of the model are inferred directly from network data obtained during the first day of the Conifcker epidemic. The model is then used to explore the trade-off between spreading modes in determining the worm's effectiveness. Our results show that the Conficker epidemic is an example of a critically hybrid epidemic, in which the different modes of spreading in isolation do not lead to successful epidemics. Such hybrid spreading strategies may be used beneficially to provide the most effective strategies for promulgating information across a large population. When used maliciously, however, they can present a dangerous challenge to current internet security protocols.


title: Efficient Algorithms for de novo Assembly of Alternative Splicing Events from RNA-seq Data

abstract: In this thesis, we address the problem of identifying and quantifying variants (alternative splicing and genomic polymorphism) in RNA-seq data when no reference genome is available, without assembling the full transcripts. Based on the fundamental idea that each variant corresponds to a recognizable pattern, a bubble, in a de Bruijn graph constructed from the RNA-seq reads, we propose a general model for all variants in such graphs. We then introduce an exact method, called KisSplice, to extract alternative splicing events. Finally, we show that it enables to identify more correct events than general purpose transcriptome assemblers. In order to deal with ever-increasing volumes of NGS data, we put an extra effort to make KisSplice as scalable as possible. First, to improve its running time, we propose a new polynomial delay algorithm to enumerate bubbles. We show that it is several orders of magnitude faster than previous approaches. Then, to reduce its memory consumption, we propose a new compact way to build and represent a de Bruijn graph. We show that our approach uses 30% to 40% less memory than the state of the art, with an insignificant impact on the construction time. Additionally, we apply the same techniques developed to list bubbles in two classical problems: cycle enumeration and the K-shortest paths problem. We give the first optimal algorithm to list cycles in undirected graphs, improving over Johnson's algorithm. This is the first improvement to this problem in almost 40 years. We then consider a different parameterization of the classical K-shortest (simple) paths problem: instead of bounding the number of st-paths, we bound the weight of the st-paths. We present new algorithms with the same time complexities but using exponentially less memory than previous approaches.


title: A Statistical Reconstruction of the Planet Population Around Kepler Solar-Type Stars

abstract: Using the cumulative catalog of planets detected by the NASA Kepler mission, we reconstruct the intrinsic occurrence of Earth- to Neptune-size (1 - 4$R_{\oplus}$) planets and their distributions with radius and orbital period. We analyze 76,711 solar-type ($0.8<R_*/R_{\odot}<1.2 $) stars with 430 planets on 20-200~d orbits, excluding close-in planets that may have been affected by the proximity to the host star. Our analysis considers errors in planet radii and includes an "iterative simulation" technique that does not bin the data. We find a radius distribution that peaks at 2-2.8 Earth radii, with lower numbers of smaller and larger planets. These planets are uniformly distributed with logarithmic period, and the mean number of such planets per star is $0.46 \pm 0.03$. The occurrence is $\sim 0.66$ if planets interior to 20~d are included. We estimate the occurrence of Earth-size planets in the "habitable zone" (defined as 1-2$R_{\oplus}$, 0.99-1.7 AU for solar-like stars) as $6.4^{+3.4}_{-1.1} \%$. Our results largely agree with those of Petigura et al. (2013), although we find a higher occurrence of 2.8-4 Earth-radii planets. The reasons for this excess are the inclusion of errors in planet radius, updated Huber et al. (2014) stellar parameters, and also the exclusion of planets which may have been affected by proximity to the host star.


title: Matter creation in a nonsingular bouncing cosmology

abstract: We examine reheating in the two-field matter bounce cosmology. In this model, the universe evolves from a matter-dominated phase of contraction to an Ekpyrotic phase of contraction before the nonsingular bounce. The Ekpyrotic phase frees the model from unwanted anisotropies, but leaves the universe cold and empty of particles after the bounce. For this reason, we explore two particle production mechanisms which take place during the course of the cosmological evolution: Parker particle production where the matter field couples only to gravity and particle creation via interactions between the matter field and the bounce field. Although we show that both mechanisms can produce particles in this model, we find that Parker particle production is sufficient to reheat the universe to high temperatures. Thus there is {\it{a priori}} no need to add an interaction term to the Lagrangian of the model. Still, particle creation via interactions can contribute to the formation of matter and radiation, but only if the coupling between the fields is tuned to be large.


title: Very late photometry of SN 2011fe

abstract: The Type Ia supernova SN 2011fe is one of the closest supernovae of the past decades. Due to its proximity and low dust extinction, this object provides a very rare opportunity to study the extremely late time evolution ($>900$ days) of thermonuclear supernovae. In this Letter, we present our photometric data of SN 2011fe taken at an unprecedented late epoch of $\approx 930$ days with GMOS-N mounted on the Gemini North telescope ($g=23.43 \pm 0.28$, $r=24.14 \pm 0.14$, $i=23.91 \pm 0.18$, and $z=23.90 \pm 0.17$) to study the energy production and retention in the ejecta of SN 2011fe. Together with previous measurements by other groups, our result suggests that the optical supernova light curve can still be explained by the full thermalization of the decay positrons of Co56. This is in spite of theoretical predicted effects (e.g. infrared catastrophe, positron escape, and dust) that advocate a substantial energy redistribution and/or loss via various processes that result in a more rapid dimming at these very late epochs.


title: Propulsion of a domain wall in an antiferromagnet by magnons

abstract: We analyze the dynamics of a domain wall in an easy-axis antiferromagnet driven by circularly polarized magnons. Magnons pass through a stationary domain wall without reflection and thus exert no force on it. However, they reverse their spin upon transmission, thereby transferring two quanta of angular momentum to the domain wall and causing it to precess. A precessing domain wall partially reflects magnons back to the source. The reflection of spin waves creates a previously identified reactive force. We point out a second mechanism of propulsion, which we term redshift: magnons passing through a precessing domain wall lower their frequency by twice the angular velocity of the domain wall; the concomitant reduction of magnons' linear momentum indicates momentum transfer to the domain wall. We solve the equations of motion for spin waves in the background of a uniformly precessing domain wall with the aid of supersymmetric quantum mechanics and compute the net force and torque applied by magnons to the domain wall. Redshift is the dominant mechanism of propulsion at low spin-wave intensities; reflection dominates at higher intensities. We derive a set of coupled algebraic equations to determine the linear velocity and angular frequency of the domain wall in a steady state. The theory agrees well with numerical micromagnetic simulations.


title: Resistance of helical edges formed in a semiconductor heterostructure

abstract: Time-reversal symmetry prohibits elastic backscattering of electrons propagating within a helical edge of a two-dimensional topological insulator. However, small band gaps in these systems make them sensitive to doping disorder, which may lead to the formation of electron and hole puddles. Such a puddle -- a quantum dot -- tunnel-coupled to the edge may significantly enhance the inelastic backscattering rate, due to the long dwelling time of an electron in the dot. The added resistance is especially strong for dots carrying an odd number of electrons, due to the Kondo effect. For the same reason, the temperature dependence of the added resistance becomes rather weak. We present a detailed theory of the quantum dot effect on the helical edge resistance. It allows us to make specific predictions for possible future experiments with artificially prepared dots in topological insulators. It also provides a qualitative explanation of the resistance fluctuations observed in short HgTe quantum wells. In addition to the single-dot theory, we develop a statistical description of the helical edge resistivity introduced by random charge puddles in a long heterostructure carrying helical edge states. The presence of charge puddles in long samples may explain the observed coexistence of a high sample resistance with the propagation of electrons along the sample edges.


title: Revisiting the bound on axion-photon coupling from Globular Clusters

abstract: We derive a strong bound on the axion-photon coupling $g_{a\gamma}$ from the analysis of a sample of 39 Galactic Globular Clusters. As recognized long ago, the R parameter, i.e. the number ratio of stars in horizontal over red giant branch of old stellar clusters, would be reduced by the axion production from photon conversions occurring in stellar cores. In this regard we have compared the measured R with state-of-the-art stellar models obtained under different assumptions for $g_{a\gamma}$. We show that the estimated value of $g_{a\gamma}$ substantially depends on the adopted He mass fraction Y, an effect often neglected in previous investigations. Taking as benchmark for our study the most recent determinations of the He abundance in H II regions with O/H in the same range of the Galactic Globular Clusters, we obtain an upper bound $g_{a\gamma}<0.66\times 10^{-10}$ GeV$^{-1}$ at 95$\%$ confidence level. This result significantly improves the constraints from previous analyses and is currently the strongest limit on the axion-photon coupling in a wide mass range.


title: Higgs-flavon mixing and LHC phenomenology in a simplified model of broken flavor symmetry

abstract: The LHC phenomenology of a low-scale gauged flavor symmetry model with inverted hierarchy is studied, through introduction of a simplified model of broken flavor symmetry. A new scalar (a flavon) and a new neutral top-philic massive gauge boson emerge with mass in the TeV range along with a new heavy fermion associated with the standard model top quark. After checking constraints from electroweak precision observables, we investigate the influence of the model on Higgs boson physics, notably on its production cross section and decay branching fractions. Limits on the flavon $\varphi$ from heavy Higgs boson searches at the LHC at 7 and 8 TeV are presented. The branching fractions of the flavon are computed as a function of the flavon mass and the Higgs-flavon mixing angle. We also explore possible discovery of the flavon at 14 TeV, particularly via the $\varphi \rightarrow Z^0Z^0$ decay channel in the $2\ell2\ell'$ final state, and through standard model Higgs boson pair production $\varphi \rightarrow hh$ in the $b\bar{b}\gamma\gamma$ final state. We conclude that the flavon mass range up to $500$ GeV could probed down to quite small values of the Higgs-flavon mixing angle with 100 fb$^{-1}$ of integrated luminosity at 14 TeV.


title: Clues on the Missing Sources of Reionization from Self-consistent Modeling of Milky Way and Dwarf Galaxy Globular Clusters

abstract: Globular clusters are unique tracers of ancient star formation. We determine the formation efficiencies of globular clusters across cosmic time by modeling the formation and dynamical evolution of the globular cluster population of a Milky Way type galaxy in hierarchical cosmology, using the merger tree from the Via Lactea II simulation. All of the models are constrained to reproduce the observed specific frequency and initial mass function of globular clusters in isolated dwarfs. Globular cluster orbits are then computed in a time varying gravitational potential after they are either accreted from a satellite halo or formed in situ, within the Milky Way halo. We find that the Galactocentric distances and metallicity distribution of globular clusters are very sensitive to the formation efficiencies of globular clusters as a function of redshift and halo mass. Our most accurate models reveal two distinct peaks in the globular cluster formation efficiency at z~2 and z~7-12 and prefer a formation efficiency that is mildly increasing with decreasing halo mass, the opposite of what expected for feedback-regulated star formation. This model accurately reproduces the positions, velocities, mass function, metallicity distribution, and age distribution of globular clusters in the Milky Way and predicts that ~ 40% formed in situ, within the Milky Way halo, while the other ~ 60% were accreted from about 20 satellite dwarf galaxies with Vc > 30 km/s, and about 29% or all globular clusters formed at redshifts z > 7. These results further strengthen the notion that globular cluster formation was an important mode of star formation in high-redshift galaxies and likely played a significant role in the reionization of the intergalactic medium


title: The distribution of satellites around massive galaxies at 1<z<3 in ZFOURGE/CANDELS: dependence on star formation activity

abstract: We study the statistical distribution of satellites around star-forming and quiescent central galaxies at 1<z<3 using imaging from the FourStar Galaxy Evolution Survey (ZFOURGE) and the Cosmic Assembly Near-IR Deep Extragalactic Legacy Survey (CANDELS). The deep near-IR data select satellites down to $\log(M/M_\odot)>9$ at z<3. The radial satellite distribution around centrals is consistent with a projected NFW profile. Massive quiescent centrals, $\log(M/M_\odot)>10.78$, have $\sim$2 times the number of satellites compared to star-forming centrals with a significance of 2.7$\sigma$ even after accounting for differences in the centrals' stellar-mass distributions. We find no statistical difference in the satellite distributions of intermediate-mass quiescent and star-forming centrals, $10.48<\log(M/M_\odot)<10.78$. Comparing to the Guo2011 semi-analytic model, the excess number of satellites indicates that quiescent centrals have halo masses 0.3 dex larger than star-forming centrals, even when the stellar-mass distributions are fixed. We use a simple toy model that relates halo mass and quenching, which roughly reproduces the observed quenched fractions and the differences in halo mass between star-forming and quenched galaxies only if galaxies have a quenching probability that increases with halo mass from $\sim$0 for $\log(M_h/M_\odot)\sim$11 to $\sim$1 for $\log(M_h/M_\odot)\sim$13.5. A single halo-mass quenching threshold is unable to reproduce the quiescent fraction and satellite distribution of centrals. Therefore, while halo quenching may be an important mechanism, it is unlikely to be the only factor driving quenching. It remains unclear why a high fraction of centrals remain star-forming even in relatively massive halos.


title: Prospects for joint gravitational-wave and electromagnetic observations of neutron-star--black-hole coalescing binaries

abstract: Coalescing neutron-star-black-hole (NS-BH) binaries are a promising source of gravitational-wave (GW) signals detectable with large-scale laser interferometers such as Advanced LIGO and Virgo. They are also one of the main short gamma-ray burst (SGRB) progenitor candidates. If the BH tidally disrupts its companion, an SGRB may be ignited when a sufficiently massive accretion disk forms around the remnant BH. Detecting an NS-BH coalescence both in the GW and electromagnetic (EM) spectrum offers a wealth of information about the nature of the source. How much can actually be inferred from a joint detection is unclear, however, as a mass/spin degeneracy may reduce the GW measurement accuracy. To shed light on this problem and on the potential of joint EM+GW observations, we here combine recent semi-analytical predictions for the remnant disk mass with estimates of the parameter-space portion that is selected by a GW detection. We identify cases in which an SGRB ignition is supported, others in which it can be excluded, and finally others in which the outcome depends on the chosen model for the currently unknown NS equation of state. We pinpoint a range of systems that would allow us to place lower bounds on the equation of state stiffness if both the GW emission and its EM counterpart are observed. The methods we develop can broaden the scope of existing GW detection and parameter-estimation algorithms and could allow us to disregard about half of the templates in an NS-BH search following an SGRB trigger, increasing its speed and sensitivity.


title: Galactic conformity and central / satellite quenching, from the satellite profiles of M$^{\ast}$ galaxies at $0.4<z<1.9$ in the UKIDSS UDS

abstract: We explore the redshift evolution of a curious correlation between the star-formation properties of central galaxies and their satellites (`galactic conformity') at intermediate to high redshift ($0.4<z<1.9$). Using an extremely deep near-infrared survey, we study the distribution and properties of satellite galaxies with stellar masses, ${\rm log} ({\rm M}_*/{\rm M}_{\odot})>9.7$, around central galaxies at the characteristic Schechter function mass, ${\rm M} \sim {\rm M}^{\ast}$. We fit the radial profiles of satellite number densities with simple power laws, finding slopes in the range -1.1 to -1.4 for mass-selected satellites, and -1.3 to -1.6 for passive satellites. We confirm the tendency for passive satellites to be preferentially located around passive central galaxies at $3\sigma$ significance and show that it exists to at least $z\sim2$. Meanwhile, the quenched fraction of satellites around star-forming galaxies is consistent with field galaxies of equal stellar masses. We find no convincing evidence for a redshift-dependent evolution of these trends. One simple interpretation of these results is that only passive central galaxies occupy an environment that is capable of independently shutting off star-formation in satellite galaxies. By examining the satellites of higher stellar mass star-forming galaxies (${\rm log} ({\rm M}_*/{\rm M}_{\odot}) > 11$), we conclude that the origin of galactic conformity is unlikely to be exclusively due to the host dark-matter halo mass. A halo-mass-independent correlation could be established by either formation bias or a more physical connection between central and satellite star-formation histories. For the latter, we argue that a star-formation (or AGN) related outburst event from the central galaxy could establish a hot halo environment which is then capable of quenching both central and satellite galaxies.


title: Correcting for Telluric Absorption: Methods, Case Studies, and Release of the TelFit Code

abstract: Ground-based astronomical spectra are contaminated by the Earth's atmosphere to varying degrees in all spectral regions. We present a Python code that can accurately fit a model to the telluric absorption spectrum present in astronomical data, with residuals of $\sim 3-5\%$ of the continuum for moderately strong lines. We demonstrate the quality of the correction by fitting the telluric spectrum in a nearly featureless A0V star, HIP 20264, as well as to a series of dwarf M star spectra near the 819 nm sodium doublet. We directly compare the results to an empirical telluric correction of HIP 20264 and find that our model-fitting procedure is at least as good and sometimes more accurate. The telluric correction code, which we make freely available to the astronomical community, can be used as a replacement for telluric standard star observations for many purposes.


title: Circular Polarization of Pulsar Wind Nebulae and the Cosmic-Ray Positron Excess

abstract: Recent observations by the PAMELA and AMS-02 telescopes have uncovered an anomalous rise in the positron fraction at energies above 10~GeV. One possible explanation for this excess is the production of primary electron/positron pairs through electromagnetic cascades in pulsar magnetospheres. This process results in a high multiplicity of electron/positron pairs within the wind-termination shock of pulsar wind nebula (PWN). A consequence of this scenario is that no circular polarization should be observed within PWN, since the contributions from electrons and positrons exactly cancel. Here we note that current radio instruments are capable of setting meaningful limits on the circular polarization of synchrotron radiation in PWNs, which observationally test the model for pulsar production of the local positron excess. The observation of a PWN with detectable circular polarization would cast strong doubt on pulsar interpretations of the positron excess, while observations setting strong limits on the circular polarization of PWN would lend credence to these models. Finally, we indicate which pulsar wind nebulae are likely to provide the best targets for observational tests of the AMS-02 excess.


title: Renormalization Group Analysis of a Fermionic Hot Spot Model

abstract: We present a renormalization group (RG) analysis of a fermionic "hot spot" model of interacting electrons on the square lattice. We truncate the Fermi surface excitations to linearly dispersing quasiparticles in the vicinity of eight hot spots on the Fermi surface, with each hot spot separated from another by the wavevector $(\pi, \pi)$. This motivated by the importance of these Fermi surface locations to the onset of antiferromagnetic order; however, we allow for all possible quartic interactions between the fermions, and also for all possible ordering instabilities. We compute the RG equations for our model, which depend on whether the hot spots are perfectly nested or not, and relate our results to earlier models. We also compute the RG flow of the relevant order parameters for both Hubbard and $J$, $V$ interactions, and present our results for the dominant instabilities in the nested and non-nested cases. In particular, we find that non-nested hot spots with $J$, $V$ interactions have competing singlet $d_{x^2-y^2}$ superconducting and $d$-form factor incommensurate density wave instabilities. We also investigate the enhancement of incommensurate density waves near experimentally observed wavevectors, and find dominant $d$-form factor enhancement for a range of couplings.


title: Co-orbiting planes of sub-halos are similarly unlikely around paired and isolated hosts

abstract: Sub-halos in dark-matter-based cosmological simulations tend to be distributed approximately isotropically around their host. The existence of highly flattened, co-orbiting planes of satellite galaxies has therefore been identified as a possible problem for these cosmological models, but so far studies have not considered the hosts' environments. That satellite planes are now known around both major galaxies in the Local Group raises the question whether they are more likely around paired hosts. In a first attempt to investigate this possibility we focus on the flattening and orbital coherence of the 11 brightest satellite galaxies of the vast polar structure (VPOS) around the Milky Way (MW). We search for VPOS analogs in the ELVIS suite of cosmological simulations, which consist of 24 paired and 24 isolated host halos. We do not find significant differences between the properties of sub-halo distributions around paired and isolated hosts. The observed flattening and the observed orbital alignment are each reproduced by only 0.2 to 2 per cent of paired and isolated systems incorporating the obscuration of satellites by randomly oriented galactic discs. Only one of all 4800 analyzed realizations (0.02 per cent) reproduces both parameters simultaneously, but the average orbital pole of this sub-halo system does not align as well with the normal to the plane fit as observed. That the MW is part of a galaxy pair thus does not help in explaining the existence of the VPOS if the satellite galaxies are identified with sub-halos found in dissipationless simulations.


title: Milky Way Mass and Potential Recovery Using Tidal Streams in a Realistic Halo

abstract: We present a new method for determining the Galactic gravitational potential based on forward modeling of tidal stellar streams. We use this method to test the performance of smooth and static analytic potentials in representing realistic dark matter halos, which have substructure and are continually evolving by accretion. Our FAST-FORWARD method uses a Markov Chain Monte Carlo algorithm to compare, in 6D phase space, an "observed" stream to models created in trial analytic potentials. We analyze a large sample of streams evolved in the Via Lactea II (VL2) simulation, which represents a realistic Galactic halo potential. The recovered potential parameters are in agreement with the best fit to the global, present-day VL2 potential. However, merely assuming an analytic potential limits the dark matter halo mass measurement to an accuracy of 5 to 20%, depending on the choice of analytic parametrization. Collectively, mass estimates using streams from our sample reach this fundamental limit, but individually they can be highly biased. Individual streams can both under- and overestimate the mass, and the bias is progressively worse for those with smaller perigalacticons, motivating the search for tidal streams at galactocentric distances larger than 70 kpc. We estimate that the assumption of a static and smooth dark matter potential in modeling of the GD-1 and Pal5-like streams introduces an error of up to 50% in the Milky Way mass estimates.


title: HV2112, a Thorne-Zytkow Object or a Super Asymptotic Giant Branch Star

abstract: The very bright red star HV2112 in the Small Magellanic Cloud could be a massive Thorne-Zytkow Object, a supergiant-like star with a degenerate neutron core. With its luminosity of over $10^5\,\rm L_\odot$, it could also be a super asymptotic giant branch star, a star with an oxygen/neon core supported by electron degeneracy and undergoing thermal pulses with third dredge up. Both TZOs and SAGB stars are expected to be rare. Abundances of heavy elements in HV2112's atmosphere, as observed to date, do not allow us to distinguish between the two possibilities based on the latest models. Molybdenum and rubidium can be enhanced by both the irp-process in a TZO or by the s-process in SAGB stars. Lithium can be generated by hot bottom burning at the base of the convective envelope in either. HV2112's enhanced calcium could thus be the key determinant. A SAGB star is not able to synthesise its own calcium but it may be possible to produce this in the final stages of the process that forms a TZO, when the degenerate electron core of a giant star is tidally disrupted by a neutron star. Hence our calculations indicate that HV2112 is most likely a genuine TZO.


title: The AIMSS Project I: Bridging the Star Cluster - Galaxy Divide

abstract: We describe the structural and kinematic properties of the first compact stellar systems discovered by the AIMSS project. These spectroscopically confirmed objects have sizes ($\sim$6 $<$ R$_{\rm e}$ [pc] $<$ 500) and masses ($\sim$2$\times$10$^{6}$ $<$ M$_*$/M$_\odot$ $<$ 6$\times$10$^{9}$) spanning the range of massive globular clusters (GCs), ultra compact dwarfs (UCDs) and compact elliptical galaxies (cEs), completely filling the gap between star clusters and galaxies. Several objects are close analogues to the prototypical cE, M32. These objects, which are more massive than previously discovered UCDs of the same size, further call into question the existence of a tight mass-size trend for compact stellar systems, while simultaneously strengthening the case for a universal "zone of avoidance" for dynamically hot stellar systems in the mass-size plane. Overall, we argue that there are two classes of compact stellar systems: 1) massive star clusters and 2) a population closely related to galaxies. Our data provide indications for a further division of the galaxy-type UCD/cE population into two groups, one population that we associate with objects formed by the stripping of nucleated dwarf galaxies, and a second population that formed through the stripping of bulged galaxies or are lower-mass analogues of classical ellipticals. We find compact stellar systems around galaxies in low to high density environments, demonstrating that the physical processes responsible for forming them do not only operate in the densest clusters.


title: Dust and Gas in the Magellanic Clouds from the HERITAGE Herschel Key Project. I. Dust Properties and Insights into the Origin of the Submm Excess Emission

abstract: The dust properties in the Large and Small Magellanic Clouds are studied using the HERITAGE Herschel Key Project photometric data in five bands from 100 to 500 micron. Three simple models of dust emission were fit to the observations: a single temperature blackbody modified by a power- law emissivity (SMBB), a single temperature blackbody modified by a broken power-law emissivity (BEMBB), and two blackbodies with different temperatures, both modified by the same power-law emissivity (TTMBB). Using these models we investigate the origin of the submm excess; defined as the submillimeter (submm) emission above that expected from SMBB models fit to observations < 200 micron. We find that the BEMBB model produces the lowest fit residuals with pixel-averaged 500 micron submm excesses of 27% and 43% for the LMC and SMC, respectively. Adopting gas masses from previous works, the gas-to-dust ratios calculated from our the fitting results shows that the TTMBB fits require significantly more dust than are available even if all the metals present in the interstellar medium (ISM) were condensed into dust. This indicates that the submm excess is more likely to be due to emissivity variations than a second population of colder dust. We derive integrated dust masses of (7.3 +/- 1.7) x 10^5 and (8.3 +/- 2.1) times 10^4 M(sun) for the LMC and SMC, respectively. We find significant correlations between the submm excess and other dust properties; further work is needed to determine the relative contributions of fitting noise and ISM physics to the correlations.


title: Order complexes of coset posets of finite groups are not contractible

abstract: We show that the order complex of the poset of all cosets of all proper subgroups of a finite group $G$ is never $\mathbb{F}_{2}$-acyclic and therefore never contractible. This settles a question of K. S. Brown.


title: Soft X-ray characterisation of the long term properties of Supergiant Fast X-ray Transients

abstract: We perform the first high-sensitivity soft X-ray long-term monitoring with Swift/XRT of three relatively unexplored Supergiant Fast X-ray Transients (SFXTs), IGR J08408-4503, IGR J16328-4726, and IGR J16465-4507, whose hard X-ray duty cycles are the lowest measured among the SFXT sample, and compare their properties with those of the prototypical SFXTs. The behaviour of J08408 and J16328 resembles that of other SFXTs, and it is characterized by a relatively high inactivity duty cycle (IDC) and pronounced dynamic range (DR) in the X-ray luminosity. Like the SFXT prototypes, J08408 shows two distinct populations of flares, the first one associated with the brightest outbursts ($L_{\rm X}\gtrsim 10^{35-36}$ erg s$^{-1}$), the second one comprising less bright events with $L_{\rm X}\lesssim$10$^{35}$ erg s$^{-1}$. This double-peaked distribution seems to be a ubiquitous feature of the extreme SFXTs. The lower DR of J16328 suggests it is an intermediate SFXT. We find J16465 is characterized by IDC$\sim$5% and DR$\sim$40, reminiscent of classical supergiant HMXBs. The duty cycles measured with XRT are found to be comparable with those reported previously by BAT and INTEGRAL, when the higher limiting sensitivities of these instruments are taken into account and sufficiently long observational campaigns are available. We prove that no clear correlation exists between the duty cycles of the SFXTs and their orbital periods, which makes it difficult to interpret the SFXT peculiar variability by only using arguments related to the properties of supergiant star winds. Our findings favour the idea that a correct interpretation of the SFXT phenomenology requires a mechanism to strongly reduce the mass accretion rate onto the compact object during most of its orbit around the companion, as proposed in a number of theoretical works. [Abridged]


title: The mass-metallicity and fundamental metallicity relations at z>2 using VLT and Subaru near-infrared spectroscopy of zCOSMOS galaxies

abstract: In the local universe, there is good evidence that, at a given stellar mass M, the gas-phase metallicity Z is anti-correlated with the star formation rate (SFR) of the galaxies. It has also been claimed that the resulting Z(M,SFR) relation is invariant with redshift - the so-called Fundamental Metallicity Relation (FMR). Given a number of difficulties in determining metallicities, especially at higher redshifts, the form of the Z(M,SFR) relation and whether it is really independent of redshift is still very controversial. To explore this issue at z>2, we used VLT-SINFONI and Subaru-MOIRCS near-infrared spectroscopy of 20 zCOSMOS-deep galaxies at 2.1<z<2.5 to measure the strengths of up to five emission lines: [OII], Hbeta, [OIII], Halpha, and [NII]. This near-infrared spectroscopy enables us to derive O/H metallicities, and also SFRs from extinction corrected Halpha measurements. We find that the mass-metallicity relation (MZR) of these star-forming galaxies at z~2.3 is lower than the local SDSS MZR by a factor of three to five, a larger change than found by Erb et al. (2006) using [NII]/Halpha-based metallicities from stacked spectra. We discuss how the different selections of the samples and metallicity calibrations used may be responsible for this discrepancy. The galaxies show direct evidence that the SFR is still a second parameter in the mass-metallicity relation at these redshifts. However, determining whether the Z(M,SFR) relation is invariant with epoch depends on the choice of extrapolation used from local samples, because z>2 galaxies of a given mass have much higher SFRs than the local SDSS galaxies. We find that the zCOSMOS galaxies are consistent with a non-evolving FMR if we use the physically-motivated formulation of the Z(M,SFR) relation from Lilly et al. (2003), but not if we use the empirical formulation of Mannucci et al. (2010).


title: Precision Electroweak Analysis after the Higgs Boson Discovery

abstract: Until recently precision electroweak computations were fundamentally uncertain due to lack of knowledge about the existence of the Standard Model Higgs boson and its mass. For this reason substantial calculational machinery had to be carried along for each calculation that changed the Higgs boson mass and other parameters of the Standard Model. Now that the Higgs boson is discovered and its mass is known to within a percent, we are able to compute reliable semi-analytic expansions of electroweak observables. We present results of those computations in the form of expansion formulae. In addition to the convenience of having these expressions, we show how the approach makes investigating new physics contributions to precision electroweak observables much easier.


title: Towards the Standard Model in F-theory

abstract: This article explores possible embeddings of the Standard Model gauge group and its matter representations into F-theory. To this end we construct elliptic fibrations with gauge group SU(3)xSU(2)xU(1)xU(1) as suitable restrictions of a ${\rm Bl}_2{\mathbb P}^2$-fibration with rank-two Mordell-Weil group. We analyse the five inequivalent toric enhancements to gauge group SU(3)xSU(2) along two independent divisors W_3 and W_2 in the base. For each of the resulting smooth fibrations, the representation spectrum generically consists of a bifundamental (3,2), three types of (1,2) representations and five types of (3,1) representations (plus conjugates), in addition to charged singlet states. The precise spectrum of zero-modes in these representations depends on the 3-form background. We analyse the geometrically realised Yukawa couplings among all these states and find complete agreement with field theoretic expectations based on their U(1) charges. We classify possible identifications of the found representations with the Standard Model field content extended by right-handed neutrinos and extra singlets. The linear combination of the two abelian gauge group factors orthogonal to hypercharge acts as a selection rule which, depending on the specific model, can forbid dangerous dimension-four and -five proton decay operators.


title: On fast Generation of Cosmological Random Fields

abstract: The statistical translation invariance of cosmological random fields is broken by a finite survey boundary, correlating the observable Fourier modes. Standard methods for generating Gaussian fields either neglect these correlations, or are costly, or both. Here we report on a fast and exact simulation method applicable to a wide class of two-point statistics that requires the simulation of a periodic grid of only twice the survey side with fast Fourier transforms. Super-survey modes, dominating the covariance of power spectra beyond linear scales in galaxy surveys and causing the correlation of large and small scales, "beat coupling", or "super-sample" covariance, are precisely accounted for in non-linear transformations of the Gaussian field. As an application, we simulate the CFHTLS $\sim 7^\circ \times 7^\circ$ W1 galaxy density field, modeled as a Poisson sampling of a lognormal density field. We show that our simulations produce power spectra, $A^*$-power spectra, counts-in-cells probability distributions as well as covariances perfectly consistent with the data. In addition, our technique reproduces the information plateau beyond linear scales as observed previously in Sloan Digital Sky Survey galaxy catalogs and in $N$-body simulations. Our method is thus an efficient yet powerful simulation and prediction tool for galaxy survey data and covariances.


title: Testing New Ideas Regarding the Nature of Interstellar Extinction

abstract: The nature of Galactic interstellar extinction is tested using reddening line parameters for several fields in conjunction with equivalent widths $W(\lambda4430)$ for the diffuse interstellar band at $4430$ \AA. The Cardelli et al.$\;$relations [29] at infrared, optical, and ultraviolet wavelengths are inconsistent with the newly-derived quadratic variation of $R_V({\rm observed})$ on reddening slope $X$. A minimum of $R_V=2.82\pm0.06$ exists for $X=0.83\pm0.10$, and is argued to represent true Galactic extinction described by $A(\lambda)\propto \lambda^{-1.375}$. It matches expectations for a new description of extinction in the infrared, optical, and ultraviolet by Zagury [32]. Additional consequences, reddened stars with no 2175 \AA$\;$feature and a correlation of normalized $\lambda4430$ absorption with $X$, are not predicted by the Cardelli et al.$\;$relation [29]. Known variations in $X$ from 0.62 to 0.83, and corresponding variations in $R_V({\rm observed})$ from 4.0 to 2.8, presumably result from forward-scattered starlight in the ultraviolet contaminating optical light of stars affected by dust extinction. A new understanding of the true nature of interstellar extinction is important for establishing an accurate picture of the extragalactic distance scale, which in turn is related to our understanding of the nature of the Universe.


title: Formation of Warped Disks by Galactic Fly-by Encounters. I. Stellar Disks

abstract: Warped disks are almost ubiquitous among spiral galaxies. Here we revisit and test the `fly-by scenario' of warp formation, in which impulsive encounters between galaxies are responsible for warped disks. Based on N-body simulations, we investigate the morphological and kinematical evolution of the stellar component of disks when galaxies undergo fly-by interactions with adjacent dark matter halos. We find that the so-called `S'-shaped warps can be excited by fly-bys and sustained for even up to a few billion years, and that this scenario provides a cohesive explanation for several key observations. We show that disk warp properties are governed primarily by the following three parameters; (1) the impact parameter, i.e., the minimum distance between two halos, (2) the mass ratio between two halos, and (3) the incident angle of the fly-by perturber. The warp angle is tied up with all three parameters, yet the warp lifetime is particularly sensitive to the incident angle of the perturber. Interestingly, the modeled S-shaped warps are often non-symmetric depending on the incident angle. We speculate that the puzzling U- and L-shaped warps are geometrically superimposed S-types produced by successive fly-bys with different incident angles, including multiple interactions with a satellite on a highly elongated orbit.


title: Synchrotron pair halo and echo emission from blazars in the cosmic web: application to extreme TeV blazars

abstract: High frequency peaked high redshift blazars, are extreme in the sense that their spectrum is particularly hard and peaks at TeV energies. Standard leptonic scenarios require peculiar source parameters and/or a special setup in order to account for these observations. Electromagnetic cascades seeded by ultra-high energy cosmic rays (UHECR) in the intergalactic medium have also been invoked, assuming a very low intergalactic magnetic field (IGMF). Here we study the synchrotron emission of UHECR secondaries produced in blazars located in magnetised environments, and show that it can provide an alternative explanation to these challenged channels, for sources embedded in structured regions with magnetic field strengths of the order of $10^{-7}$ G. To demonstrate this, we focus on three extreme blazars: 1ES 0229+200, RGB J0710+591, and 1ES 1218+304. We model the expected gamma-ray signal from these sources through a combination of numerical Monte Carlo simulations and solving the kinetic equations of the particles in our simulations, and explore the UHECR source and intergalactic medium parameter space to test the robustness of the emission. We show that the generated synchrotron pair halo/echo flux at the peak energy is not sensitive to variations in the overall IGMF strength. This signal is unavoidable in contrast to the inverse Compton pair halo/echo intensity, which is appealing in view of the large uncertainties on the IGMF in voids of large scale structure. It is also shown that the variability of blazar gamma-ray emission can be accommodated by the synchrotron emission of secondary products of UHE neutral beams if these are emitted by UHECR accelerators inside magnetised regions.


title: Transport Coefficients of Black MQGP M3-Branes

abstract: The SYZ mirror, in the `delocalized limit' of [1], of (M)N (fractional)D3-branes, and wrapped N_f flavor D7-branes in the presence of a black-hole resulting in a non-Kaehler resolved warped deformed conifold (NKRWDC) in [2], was carried out in [3] and resulted in black M3-branes. The uplift, if valid globally(like [4] for fractional D3 branes in conifolds), asymptotes to M5-branes wrapping a two-cycle (homologously an (large) integer sum of two-spheres) in AdS_5xM_6. Interestingly, in the MQGP limit, assuming the deformation > resolution, by estimating the five SU(3) structure torsion (\tau) classes W_{1,2,3,4,5} we show that \tau\in W_4+W_5: 2/3 Re(W^3bar_5)=W^3bar_4 in the UV, implying the NKRWDC locally preserves SUSY. Further, the local T^3 of [3] in the large-r limit and the `MQGP' limit of [3], satisfies the same conditions as the maximal T^2-invariant special Lagrangian three-cycle of T^*(S^3) of [6], partly justifying use of local SYZ mirror symmetry in [3]. Using the Ouyang embedding in the DBI action of a D7-brane or by dimensionally reducing the 11-dimensional EH action to five (R^{1,3},r) dimensions, we then calculate a variety of gauge and metric-perturbation-modes' two-point functions using the prescription of [5], and show: (i) diffusion constant D~1/T, (ii) the electrical conductivity \sigma~T, (iii) the charge susceptibility \chi~ T^2, (iv) [using (i) - (iii)] the Einstein's relation \sigma/\chi=D, is indeed satisfied, (v) the R-charge diffusion constant D_R~1/T, and (vi) the possibility of generating \eta/s=1/4pi from solutions to the vector and tensor mode metric perturbations' EOMs, separately. All results are also valid in the limit of [2].


title: Kane-Mele-Hubbard model on the $\pi$-flux honeycomb lattice

abstract: We consider the Kane-Mele-Hubbard model with a magnetic $\pi$ flux threading each honeycomb plaquette. The resulting model has remarkably rich physical properties. In each spin sector, the noninteracting band structure is characterized by a total Chern number $C=\pm 2$. Fine-tuning of the intrinsic spin-orbit coupling $\lambda$ leads to a quadratic band crossing point associated with a topological phase transition. At this point, quantum Monte Carlo simulations reveal a magnetically ordered phase which extends to weak coupling. Although the spinful model has two Kramers doublets at each edge and is explicitly shown to be a $Z_{2}$ trivial insulator, the helical edge states are protected at the single-particle level by translation symmetry. Drawing on the bosonized low-energy Hamiltonian, we predict a correlation-induced gap as a result of umklapp scattering for half-filled bands. For strong interactions, this prediction is confirmed by quantum Monte Carlo simulations.


title: Flavor Symmetry Based MSSM (sMSSM): Theoretical Models and Phenomenological Analysis

abstract: We present a class of supersymmetric models in which symmetry considerations alone dictate the form of the soft SUSY breaking Lagrangian. We develop a class of minimal models, denoted as sMSSM -- for flavor symmetry-based minimal supersymmetric standard model, which respect a grand unified symmetry such as SO(10) and a non-Abelian flavor symmetry H which suppresses SUSY-induced flavor violation. Explicit examples are constructed with the flavor symmetry being gauged SU(2)_H and SO(3)_H with the three families transforming as 2 + 1 and 3 representations respectively. A simple solution is found in the case of SU(2)_H for suppressing the flavor violating D--terms based on an exchange symmetry. Explicit models based on SO(3)_H without the D--term problem are developed. In addition, models based on discrete non-Abelian flavor groups are presented which are automatically free from D--term issues. The permutation group S_3 with a 2 + 1 family assignment, as well as the tetrahedral group A_4 with a 3 assignment are studied. In all cases, a simple solution to the SUSY CP problem is found, based on spontaneous CP violation leading to a complex quark mixing matrix. We develop the phenomenology of the resulting sMSSM, which is controlled by seven soft SUSY breaking parameters for both the 2 + 1 assignment and the 3 assignment of fermion families. These models are special cases of the phenomenological MSSM (pMSSM), but with symmetry restrictions. We discuss the parameter space of sMSSM compatible with LHC searches, B physics constraints and dark matter relic abundance. Fine-tuning in these models is relatively mild, since all SUSY particles can have masses below about 3 TeV.


title: Dynamical models for the Sculptor dwarf spheroidal in a Lambda CDM universe

abstract: The Sculptor dwarf spheroidal galaxy contains two distinct stellar populations: one metal-rich and the other metal-poor. Several authors have argued that in order for these two populations to reside in the same gravitational potential, the dark matter halo must have a core similar to that observed in the stellar count profile. This would rule out the cuspy Navarro-Frenk-White (NFW) density profiles predicted for halos and subhalos by dark matter only simulations of the Lambda Cold Dark Matter (Lambda-CDM) cosmological model. We present a new theoretical framework to analyse stellar count and velocity observations in a self-consistent manner based on separable models, f(E,J)=g(J)h(E), for the distribution function of an equilibrium spherical system. We use this machinery to analyse available photometric and kinematic data for the two stellar populations in Sculptor. We find, contrary to some previous claims, that the data are consistent with populations in equilibrium within an NFW dark matter potential with structural parameters in the range expected in Lambda-CDM. Our solution gives a maximum circular velocity for Sculptor between 20 and 35 km/s. We discuss why some previous authors came to a different conclusion.


title: A reappraisal of the wrong-sign $hb\overline{b}$ coupling and the study of $h \rightarrow Z \gamma$

abstract: It has been pointed out recently that current experiments still allow for a two Higgs doublet model where the $h b \bar{b}$ coupling ($k_D m_b/v$) is negative; a sign opposite to that of the Standard Model. Due to the importance of delayed decoupling in the $h H^+ H^-$ coupling, $h \rightarrow \gamma \gamma$ improved measurements will have a strong impact on this issue. For the same reason, measurements or even bounds on $h \rightarrow Z \gamma$ are potentially interesting. In this article, we revisit this problem, highlighting the crucial importance of $h \rightarrow VV$, which can be understood with simple arguments. We show that the impacts on $k_D<0$ models of both $h \rightarrow b \bar{b}$ and $h \rightarrow \tau^+ \tau^-$ are very sensitive to input values for the gluon fusion production mechanism; in contrast, $h \rightarrow \gamma \gamma$ and $h \rightarrow Z \gamma$ are not. We also inquire if the search for $h \rightarrow Z \gamma$ and its interplay with $h \rightarrow \gamma \gamma$ will impact the sign of the $h b \bar{b}$ coupling. Finally, we study these issues in the context of the Flipped two Higgs doublet model.


title: A. Fluka Study of Underground Cosmogenic Neutron Production

abstract: Neutrons produced by cosmic muon interactions are important contributors to backgrounds in underground detectors when searching for rare events. Typically such neutrons can dominate the background, as they are particularly difficult to shield and detect. Since actual data is sparse and not well documented, simulation studies must be used to design shields and predict background rates. Thus validation of any simulation code is necessary to assure reliable results. This work compares in detail the predictions of the FLUKA simulation code to existing data, and uses this code to report a simulation of cosmogenic backgrounds for typical detectors embedded in a water tank with liquid scintillator shielding.


title: Spiral Galaxies - classical description of spiral arms and rotational velocity pattern - toy model

abstract: We propose an explanation of features of spiral galaxies: spiral arms and observed flat rotation curves, without the presence of an exotic form of matter. The formalism is based on Boltzmanns transport equation for the collisional matter and the very-low-velocity post-Newtonian approximation of the general relativity equations expressed in the Maxwell-like form. The Maxwell-like formulation provides the base for the explanation of the above phenomena in the language of dynamically created gravitoelectromagnetic fields by the movement of mass streams in the plane of the galaxy disc. In the model we use radical simplifications expressed as neglect of the gravitational interaction between neighbors and approximation of the incompressible mass flow. In this frame we show that if the galaxy disc is fuelled constantly by collisional mass carriers, then the amplification of the gravitomagnetic field can be large enough to create the rotational velocity pattern and spiral arms. In this framework the collisional part of the mass gas in the galaxy disc plane, i.e. molecules and atoms, is responsible for the creation of the gravitomagnetic field. According to the model the spiral pattern of arms is static and determined by a direction of the mass flow. The model reproduces qualitatively the observed spiral arms and reproduces well the shape of the rotational velocity pattern. As an example of the usability of the proposed mechanism, we reproduce qualitatively the above features for the IC 342 and NGC 4321 galaxies.


title: On the auto Igusa-zeta function of an Algebraic Curve

abstract: We study endomorphisms of complete Noetherian local rings in the context of motivic integration. Using the notion of an auto-arc space, we introduce the (reduced) auto-Igusa zeta series at a point, which appears to measure the degree to which a variety is not smooth that point. We conjecture a closed formula in the case of curves with one singular point, and we provide explicit formulas for this series in the case of the cusp and the node. Using the work of Denef and Loeser, one can show that this series will often be rational. These ideas were obtained through extensive calculations in Sage. Thus, we include a Sage script which was used in these calculations. It computes the affine arc spaces $\nabla_{\mathfrak{n}}X$ provided that $X$ is affine, $\mathfrak{n}$ is a fat point, and the ground field is of characteristic zero. Finally, we show that the auto Poincar\'e series will often be rational as well and connect this to questions concerning new types of motivic integrals.


title: From Black-Scholes to Online Learning: Dynamic Hedging under Adversarial Environments

abstract: We consider a non-stochastic online learning approach to price financial options by modeling the market dynamic as a repeated game between the nature (adversary) and the investor. We demonstrate that such framework yields analogous structure as the Black-Scholes model, the widely popular option pricing model in stochastic finance, for both European and American options with convex payoffs. In the case of non-convex options, we construct approximate pricing algorithms, and demonstrate that their efficiency can be analyzed through the introduction of an artificial probability measure, in parallel to the so-called risk-neutral measure in the finance literature, even though our framework is completely adversarial. Continuous-time convergence results and extensions to incorporate price jumps are also presented.


title: Spectrum Estimation: A Unified Framework for Covariance Matrix Estimation and PCA in Large Dimensions

abstract: Covariance matrix estimation and principal component analysis (PCA) are two cornerstones of multivariate analysis. Classic textbook solutions perform poorly when the dimension of the data is of a magnitude similar to the sample size, or even larger. In such settings, there is a common remedy for both statistical problems: nonlinear shrinkage of the eigenvalues of the sample covariance matrix. The optimal nonlinear shrinkage formula depends on unknown population quantities and is thus not available. It is, however, possible to consistently estimate an oracle nonlinear shrinkage, which is motivated on asymptotic grounds. A key tool to this end is consistent estimation of the set of eigenvalues of the population covariance matrix (also known as the spectrum), an interesting and challenging problem in its own right. Extensive Monte Carlo simulations demonstrate that our methods have desirable finite-sample properties and outperform previous proposals.


title: Role of electron-phonon interaction in a magnetically driven mechanism for superconductivity

abstract: We use the renormalization group method to examine the effect of phonon mediated interaction on d-wave superconductivity, as driven by spin fluctuations in a quasi-one-dimensional electron system. The influence of a tight-binding electron-phonon interaction on the spin-density-wave and d-wave superconducting instability lines is calculated for arbitrary temperature, phonon frequency and antinesting of the Fermi surface.The domain of electron-phonon coupling strength where spin-density-wave order becomes unstable against the formation of a bond-order-wave or Peierls state is determined at weak antinesting. We show the existence of a positive isotope effect for spin-density-wave and d-wave superconducting critical temperatures which scales with the antinesting distance from quantum critical point where the two instabilities merge. We single out a low phonon frequency zone where the bond-oder-wave ordering gives rise to triplet f-wave superconductivity under nesting alteration, with both orderings displaying a negative isotope effect. We also study the electron-phonon strengthening of spin fluctuations at the origin of extended quantum criticality in the metallic phase above superconductivity. The impact of our results on quasi-one-dimensional organic conductors like the Bechgaard salts where a Peierls distortion is absent and superconductivity emerges near a spin-density-wave state under pressure is emphasized.


title: Initialization Algorithms For Convolutional Network Coding

abstract: We present algorithms for initializing a convolutional network coding scheme in networks that may contain cycles. An initialization process is needed if the network is unknown or if local encoding kernels are chosen randomly. During the initialization process every source node transmits basis vectors and every sink node measures the impulse response of the network. The impulse response is then used to find a relationship between the transmitted and the received symbols, which is needed for a decoding algorithm and to find the set of all achievable rates. Unlike acyclic networks, for which it is enough to transmit basis vectors one after another, the initialization of cyclic networks is more involved, as pilot symbols interfere with each other and the impulse response is of infinite duration.


title: From Phase to Micro-Phase Separation in Flocking Models: The Essential Role of Non-Equilibrium Fluctuations

abstract: We show that the flocking transition in the Vicsek model is best understood as a liquid-gas transition, rather than an order-disorder one. The full phase separation observed in flocking models with Z2 rotational symmetry is, however, replaced by a microphase separation leading to a smectic arrangement of traveling ordered bands. Remarkably, continuous deterministic descriptions do not account for this difference, which is only recovered at the fluctuating hydrodynamics level. Scalar and vectorial order parameters indeed produce different types of number fluctuations, which we show to be essential in selecting the inhomogeneous patterns. This highlights an unexpected role of fluctuations in the selection of flock shapes.


title: Controlling Tokamak Geometry with 3D Magnetic Perturbations

abstract: It is shown that small externally applied magnetic perturbations can significantly alter important geometric properties of magnetic flux surfaces in tokamaks. Through 3D shaping, experimentally relevant perturbation levels are large enough to influence turbulent transport and MHD stability in the pedestal region. It is shown that the dominant pitch-resonant flux surface deformations are primarily induced by non-resonant 3D fields, particularly in the presence of significant axisymmetric shaping. The spectral content of the applied 3D field can be used to control these effects.


title: Semiclassical approximation in stochastic optimal control I. Portfolio construction problem

abstract: This is the first in a series of papers in which we study an efficient approximation scheme for solving the Hamilton-Jacobi-Bellman equation for multi-dimensional problems in stochastic control theory. The method is a combination of a WKB style asymptotic expansion of the value function, which reduces the second order HJB partial differential equation to a hierarchy of first order PDEs, followed by a numerical algorithm to solve the first few of the resulting first order PDEs. This method is applicable to stochastic systems with a relatively large number of degrees of freedom, and does not seem to suffer from the curse of dimensionality. Computer code implementation of the method using modest computational resources runs essentially in real time. We apply the method to solve a general portfolio construction problem.


title: 50 M31 black hole candidates identified by Chandra and XMM-Newton

abstract: Over the last ~5 years we have identified ~35 black hole candidates (BHCs) in M31 from their X-ray spectra. Our BHCs exhibited 0.3--10 keV spectra consistent with the X-ray binary (XB) hard state, at luminosities that are above the upper limit for neutron star (NS) XBs. When our BHC spectra were modeled with a disk blackbody + blackbody model, for comparison with bright NS XBs, we found that the BHCs inhabited a different parameter space to the NS XBs. However, BH XBs may also exhibit a thermally dominated (TD) state that has never been seen in NS XBs; this TD state is most often observed in X-ray transients. We examined the ~50 X-ray transients in our Chandra survey of M31, and found 13 with spectra suitable for analysis. We also examined 2 BHCs outside the field of view of our survey, in the globular clusters B045 and B375. We have 42 strong BHCs, and 8 plausible BHCs that may benefit from further observation. Of our 15 BHCs in globular clusters, 12 differ from NS spectra by >5 sigma. Due to improvements in our analysis, we have upgraded 10 previously identified plausible BHCs to strong BHCs. The mean maximum duty cycle of the 33 X-ray transients within 6' of M31* was 0.13; we estimate that >40% of the XBs in this region contain BH accretors. Remarkably, we estimate that BHCs contribute >90% of those XBs >1 E+38 erg/s.


title: The genetic code invariance: when Euler and Fibonacci meet

abstract: The number of atoms in the four ribonucleotides uridine monophosphate, cytidine monophosphate, adenine monophosphate and guanine monophosphate is taken as a key parameter. A mathematical relation describing the condensation of the three basic subunits a nucleobase, a ribose and a phosphate group, to form a ribonucleotide, is first obtained from this parameter. Next, the use of the latter and Euler totient function is shown to lead to the atom number content of the 64 codons and also to Rakocevic pattern. Finally, selected finite sums of Fibonacci numbers are shown to lead to the nucleon number content of the amino acids in various degeneracy patterns, and also to the multiplet structure of the 20 amino acids as well as to the degeneracy.


title: Interferometry of chemically peculiar stars: theoretical predictions vs. modern observing facilities

abstract: By means of numerical experiments we explore the application of interferometry to the detection and characterization of abundance spots in chemically peculiar (CP) stars using the brightest star eps~Uma as a case study. We find that the best spectral regions to search for spots and stellar rotation signatures are in the visual domain. The spots can clearly be detected already at a first visibility lobe and their signatures can be uniquely disentangled from that of rotation. The spots and rotation signatures can also be detected in NIR at low spectral resolution but baselines longer than 180~m are needed for all potential CP candidates. According to our simulations, an instrument like VEGA (or its successor e.g., FRIEND) should be able to detect, in the visual, the effect of spots and spots+rotation, provided that the instrument is able to measure $V^2\approx10^{-3}$, and/or closure phase. In infrared, an instrument like AMBER but with longer baselines than the ones available so far would be able to measure rotation and spots. Our study provides necessary details about strategies of spot detection and the requirements for modern and planned interferometric facilities essential for CP star research.


title: The M_BH - M_star relation for X-ray obscured, red QSOs at 1.2< z <2.6

abstract: We present near-infrared spectra, obtained with SINFONI and XShooter observations at ESO VLT, of nine dusty, red QSOs at 1.2<z<2.6. The sources are hard X-ray detected, characterized by cold absorption (N_H>10^{21} - 10^{22} cm^{-2}) and show a broad Ha component in the NIR spectra. We complement this sample with twelve additional sources taken from the literature with similar properties resulting in a total sample of 21 X-ray obscured, intermediate type (1.8-1.9), dusty reddened QSOs. From the broad Ha line we have computed the BH masses through the virial formula and derived Eddington ratios. Moreover, from optical/IR multi-component SED fitting we have derived the stellar mass of their host galaxies and their SFRs. We find that most of the sources in our sample are hosted in starburst and main sequence star-forming galaxies with Eddington ratios lambda>0.1. We find a strong trend with the BH mass i.e. less massive objects are scattered below and above the local relation while the most massive ones are mainly located above it. We also studied the evolution of these sources on the M_BH-M_star plane compared to a sample of optically blue type--1 QSOs and we find that obscured red QSOs show a ratio of M_BH to M_star that increases with redshift which is consistent with or slightly lower than what has been found for blue QSOs. These sources may represent the blow-out phase at the end of the rapid BH growth and immediately preceding the classical blue QSOs typically sampled in optical surveys. They in fact show evidence of outflows in the ionized gas component, but their BH has already fully formed.