title: Ab-Initio Theory of Superconductivity in a Magnetic Field II. : Numerical solution

abstract: We numerically investigate the Spin Density Functional theory for superconductors (SpinSCDFT) and the approximated exchange-correlation functional, derived and presented in the preceding paper I. As a test system we employ a free electron gas featuring an exchange-splitting, a phononic pairing field and a Coulomb repulsion. SpinSCDFT results are compared with Sarma, the Bardeen Cooper and Schrieffer theory and with an Eliashberg type of approach. We find that the spectrum of the superconducting Kohn-Sham SpinSCDFT system is not in agreement with the true quasi particle structure. Therefore, starting from the Dyson equation, we derive a scheme that allows to compute the many body excitations of the superconductor and represents the extension to superconductivity of the G0W0 method in band structure theory. This superconducting G0 W0 method vastly improves the predicted spectra.


title: Widening the Axion Window via Kinetic and St\"uckelberg Mixings

abstract: We point out that kinetic and St\"uckelberg mixings that are generically present in the low energy effective action of axions can significantly widen the window of axion decay constants. We show that an effective super-Planckian decay constant can be obtained even when the axion kinetic matrix has only sub-Planckian entries. Our minimal model involves only two axions, a St\"uckelberg U(1) and a modest rank instanton generating non-Abelian group. Below the mass of the St\"uckelberg U(1), there is only a single axion with a non-perturbatively generated potential. In contrast to previous approaches, the enhancement of the axion decay constant is not tied to the number of degrees of freedom introduced. We also discuss how kinetic mixings can lower the decay constant to the desired axion dark matter window. String theory embeddings of this scenario and their phenomenological features are briefly discussed.


title: Bifurcations in models of a society of reasonable contrarians and conformists

abstract: We study models of a society composed of a mixture of conformist and reasonable contrarian agents that at any instant hold one of two opinions. Conformists tend to agree with the average opinion of their neighbors and reasonable contrarians to disagree, but revert to a conformist behavior in the presence of an overwhelming majority, in line with psychological experiments. The model is studied in the mean field approximation and on small-world and scale-free networks. In the mean field approximation, a large fraction of conformists triggers a polarization of the opinions, a pitchfork bifurcation, while a majority of reasonable contrarians leads to coherent oscillations, with an alternation of period-doubling and pitchfork bifurcations up to chaos. Similar scenarios are obtained by changing the fraction of long-range rewiring and the parameter of scale-free networks related to the average connectivity.


title: Special Relativity in a Discrete Quantum Universe

abstract: The hypothesis of a discrete fabric of the universe--the "Planck scale"--is always on stage, since it solves mathematical and conceptual problems in the infinitely small. However, it clashes with special relativity, which is designed for the continuum. Here we show how the clash can be overcome within a discrete quantum theory where the evolution of fields is described by a quantum cellular automaton. The reconciliation is achieved by defining the change of observer as a change of representation of the dynamics, without any reference to space-time. We use the relativity principle, i.e. the invariance of dynamics under change of inertial observer, to identify a change of inertial frame with a symmetry of the dynamics. We consider the full group of such symmetries, and recover the usual Lorentz group in the relativistic regime of low energies, while at the Planck scale the covariance is nonlinearly distorted.


title: Poisson distribution of a prime counting function corresponding to elliptic curves

abstract: Let $E$ be an elliptic curve defined over rational field $\mathbb{Q}$ and $N$ be a positive integer. Now, $M_E(N)$ denotes the number of primes $p$, such that the group $E_p(\mathbb{F}_p)$ is of order $N$. We show that $M_E(N)$ follows Poisson distribution when an average is taken over a large class of curves.


title: An experimental apparatus for observing deterministic structure formation in plate-on-pedestal ice crystal growth

abstract: We describe an experimental apparatus for making detailed morphological observations of the growth of isolated plate-like ice crystals from water vapor. Each crystal develops a plate-on-pedestal (POP) geometry, in which a large, thin, plate-like crystal grows out from the top edge of an initially prismatic seed crystal resting on a substrate. With the POP geometry, the substrate is not in contact with the growing plate (except at its center), so substrate interactions do not adversely affect the crystal growth. By controlling the temperature and supersaturation around the crystal, we can manipulate the resulting ice growth behavior in predictable ways, producing morphologies spanning the full range from simple faceted hexagonal plates to complex dendritic structures. We believe that the experimental apparatus described here will allow unprecedented investigations of ice crystal growth behaviors under controlled conditions, identifying and exploring robust morphological features in detail. Such investigations will provide valuable observational inputs for developing numerical modeling techniques that can accurately reproduce the faceted and branched structures that frequently emerge during diffusion-limited crystal growth.


title: Tangent lines and Lipschitz differentiability spaces

abstract: We study the existence of tangent lines, i.e. subsets of the tangent space isometric to the real line, in tangent spaces of metric spaces. We first revisit the almost everywhere metric differentiability of Lipschitz continuous curves. We then show that any blow-up done at a point of metric differentiability and of density one for the domain of the curve gives a tangent line. Metric differentiability enjoys a Borel measurability property and this will permit us to use it in the framework of Lipschitz differentiability spaces. We show that any tangent space of a Lipschitz differentiability space contains at least $n$ distinct tangent lines, obtained as the blow-up of $n$ Lipschitz curves, where $n$ is the dimension of the local measurable chart. Under additional assumptions on the space, such as curvature lower bounds, these $n$ distinct tangent lines span an $n$-dimensional part of the tangent space.


title: A necessary condition for lower semicontinuity of line energies

abstract: We are interested in some energy functionals concentrated on the discontinuity lines of divergence-free 2D vector fields valued in the circle $\mathbb{S}^1$. This kind of energy has been introduced first by P. Aviles and Y. Giga. They show in particular that, with the cubic cost function $f(t)=t^3$, this energy is lower semicontinuous. In this paper, we construct a counter-example which excludes the lower semicontinuity of line energies for cost functions of the form $t^p$ with $0<p<1$. We also show that, in this case, the viscosity solution corresponding to a certain convex domain is not a minimizer.


title: On the double points of operator stable L\'evy processes

abstract: We determine the Hausdorff dimension of the set of double points for a symmetric operator stable L\'evy process in terms of the eigenvalues of its stability exponent.


title: Cylindrical confinement of semiflexible polymers

abstract: Equilibrium states of a closed semiflexible polymer binding to a cylinder are described. This may be either by confinement or by constriction. Closed completely bound states are labeled by two integers: the number of oscillations, $n$, and the number of times it winds the cylinder, $p$, the latter a topological invariant. We examine the behavior of these states as the length of the loop is increased by evaluating the energy, the conserved axial torque and the contact force. The ground state for a given $p$ is the state with $n=1$; a short loop with $p=1$ is an elliptic deformation of a parallel; as its length increases it elongates along the cylinder axis, with two hairpin ends. Excited states with $n\geq2$ and $p=1$ possess $n$-fold axial symmetry. Short (long) loops possess energies $\approx p E_0$ ($nE_0$), with $E_0$ the energy of a circular loop with same radius as the cylinder; in long loops the axial torque vanishes. Confined bound excited states are initially unstable; however, above a critical length each $n$-fold state becomes stable: the folded hairpin cannot be unfolded. The ground state for each $p$ is also initially unstable with respect to deformations rotating the loop off the surface into the interior. A closed planar elastic curve aligned along the cylinder axis making contact with the cylinder on its two sides is identified as the ground state of a confined loop. Exterior bound states behave very differently, if free to unbind, as signaled by the reversal in the sign of the contact force. If $p=1$, all such states are unstable. If $p \geq 2$, however, a topological obstruction to complete unbinding exists. If the loop is short, the bound state with $p=2$ and $n=1$ provides a stable constriction of the cylinder, partially unbinding as the length is increased. This motif could be relevant to an understanding of the process of membrane fission mediated by dynamin rings.


title: Constructing Galileons

abstract: In this plenary talk delivered at the DISCRETE 2014 conference in London, I briefly summarize the ideas behind and attractive properties of the Galileon field theories, and describe a broad class of scalar field theories that share these properties. After describing how Galileons arise, and commenting on their fascinating properties, in the latter half of the talk I focus on novel ways of constructing Galileon-like theories, using both the probe brane construction, and the coset construction.


title: Ideal clocks - a convenient fiction

abstract: We show that no device built according to the rules of quantum field theory can measure proper time along its path. Highly accelerated quantum clocks experience the Unruh effect, which inevitably influences their time rate. This contradicts the concept of an ideal clock, whose rate should only depend on the instantaneous velocity.


title: A new general normal mode approach to dynamic tides in rotating stars with realistic structure and its applications

abstract: We review our recent results on a unified normal mode approach to dynamic tides proposed in Ivanov, Papaloizou $\&$ Chernov (2013) and Chernov, Papaloizou $\&$ Ivanov (2013). Our formalism can be used whenever the tidal interactions are mainly determined by normal modes of a star with identifiable regular spectrum of low frequency modes. We provide in the text basic expressions for tidal energy and angular momentum transfer valid both for periodic and parabolic orbits, and different assumptions about efficiency of normal mode damping due to viscosity and/or non-linear effects and discuss applications to binary stars and close orbiting extrasolar planets.


title: Large Deviations for the Langevin equation with strong damping

abstract: We study large deviations in the Langevin dynamics, with damping of order $\e^{-1}$ and noise of order $1$, as $\e\downarrow 0$. The damping coefficient is assumed to be state dependent. We proceed first with a change of time and then, we use a weak convergence approach to large deviations and their equivalent formulation in terms of the Laplace principle, to determine the good action functional. Some applications of these results to the exit problem from a domain and to the wave front propagation for a suitable class of reaction diffusion equations are considered.


title: Effect of Longitudinally Varying Cloud Coverage on Visible Wavelength Reflected-Light Exoplanet Phase Curves

abstract: We use a planetary albedo model to investigate variations in visible wavelength phase curves of exoplanets. The presence of clouds on these exoplanets significantly alters their planetary albedo spectra. We confirm that non-uniform cloud coverage on the dayside of tidally locked exoplanets will manifest as changes to the magnitude and shift of the phase curve. In this work, we first investigate a test case of our model using a Jupiter-like planet, at temperatures consistent to 2.0 AU insolation from a solar type star, to consider the effect of H2O clouds. We then extend our application of the model to the exoplanet Kepler-7b and consider the effect of varying cloud species, sedimentation efficiency, particle size, and cloud altitude. We show that, depending on the observational filter, the largest possible shift of the phase curve maximum will be 2-10 deg for a Jupiter-like planet, and up to 30 deg (0.08 in fractional orbital phase) for hot-Jupiter exoplanets at visible wavelengths as a function of dayside cloud distribution with a uniformly averaged thermal profile. Finally, we tailor our model for comparison with, and confirmation of, the recent optical phase-curve observations of Kepler-7b with the Kepler space telescope. The average planetary albedo can vary between 0.1-0.6 for the 1300 cloud scenarios that were compared to the observations. We observe that smaller particle size and increasing cloud altitude have a strong effect on increasing albedo. In particular, we show that a set of models where Kepler-7b has roughly half of its dayside covered in small-particle clouds high in the atmosphere, made of bright minerals like MgSiO3 and Mg2SiO4, provide the best fits to the observed offset and magnitude of the phase-curve, whereas Fe clouds are found to have too dark to fit the observations.


title: Discrete Malliavin-Stein method: Berry-Esseen bounds for random graphs and percolation

abstract: A new Berry-Esseen bound for non-linear functionals of non-symmetric and non-homogeneous infinite Rademacher sequences is established. It is based on a discrete version of the Malliavin-Stein method and an analysis of the discrete Ornstein-Uhlenbeck semigroup. The result is applied to sub-graph counts and to the number of vertices having a prescribed degree in the Erd\H{o}s-Renyi random graph. A further application deals with a percolation problem on trees.


title: Magneto-optical Kramers-Kronig analysis

abstract: We describe a simple magneto-optical experiment and introduce a magneto-optical Kramers-Kronig analysis (MOKKA) that together allow extracting the complex dielectric function for left- and right-handed circular polarizations in a broad range of frequencies without actually generating circularly polarized light. The experiment consists of measuring reflectivity and Kerr rotation, or alternatively transmission and Faraday rotation, at normal incidence using only standard broadband polarizers without retarders or quarter-wave plates. In a common case, where the magneto-optical rotation is small (below $\sim$ 0.2 rad), a fast measurement protocol can be realized, where the polarizers are fixed at 45$^\circ$ with respect to each other. Apart from the time-effectiveness, the advantage of this protocol is that it can be implemented at ultra-high magnetic fields and in other situations, where an \emph{in-situ} polarizer rotation is difficult. Overall, the proposed technique can be regarded as a magneto-optical generalization of the conventional Kramers-Kronig analysis of reflectivity on bulk samples and the Kramers-Kronig constrained variational analysis of more complex types of spectral data. We demonstrate the application of this method to the textbook semimetals bismuth and graphite and also use it to obtain handedness-resolved magneto-absorption spectra of graphene on SiC.


title: Gears, Pregears and Related Domains

abstract: We study conformal mappings from the unit disk to one-toothed gear-shaped planar domains from the point of view of the Schwarzian derivative. Gear-shaped (or "gearlike") domains fit into a more general category of domains we call "pregears" (images of gears under Mobius transformations), which aid in the study of the conformal mappings for gears and which we also describe in detail. Such domains being bounded by arcs of circles, the Schwarzian derivative of the Riemann mapping is known to be a rational function of a specific form. One accessory parameter of these mappings is naturally related to the conformal modulus of the gear (or pregear) and we prove several qualitative results relating it to the principal remaining accessory parameter. The corresponding region of univalence (parameters for which the rational function is the Schwarzian derivative of a conformal mapping) is determined precisely.


title: The power conjugacy problem in Higman-Thompson groups

abstract: An introduction to the universal algebra approach to Higman-Thompson groups (including Thompson's group $V$) is given, following a series of lectures by Graham Higman in 1973. In these talks, Higman outlined an algorithm for the conjugacy problem; which although essentially correct fails in certain cases, as we show here. A revised and complete version of the algorithm is written out explicitly. From this, we construct an algorithm for the power conjugacy problem in these groups. Python implementations of these algorithms can be found at [26].


title: On the sharp regularity for arbitrary actions of nilpotent groups on the interval: the case of N_4

abstract: We determine the largest {\alpha} for which the nilpotent group of 4-by-4 triangular matrices with integer coefficients and 1 in the diagonal embeds into the group of C1+{\alpha} diffeomorphism of the closed interval.


title: Quantomatic: A Proof Assistant for Diagrammatic Reasoning

abstract: Monoidal algebraic structures consist of operations that can have multiple outputs as well as multiple inputs, which have applications in many areas including categorical algebra, programming language semantics, representation theory, algebraic quantum information, and quantum groups. String diagrams provide a convenient graphical syntax for reasoning formally about such structures, while avoiding many of the technical challenges of a term-based approach. Quantomatic is a tool that supports the (semi-)automatic construction of equational proofs using string diagrams. We briefly outline the theoretical basis of Quantomatic's rewriting engine, then give an overview of the core features and architecture and give a simple example project that computes normal forms for commutative bialgebras.


title: A Fit to the Galactic Cosmic Ray Hydrogen and Helium Spectra at Voyager 1 at Low Energies and Earth Based Measurements at Much Higher Energies with Identical Rigidity Independent Source Spectra for the Hydrogen and Helium Nuclei

abstract: Voyager 1 data from beyond the heliopause provide the first direct measurements of the interstellar cosmic ray spectra below 1 GeVnuc. In this paper we combine these Voyager measurements of H and He nuclei from 3-600 MeVnuc with higher energy measurements at 1 AU from the BESS and PAMELA experiments up to 100 GeVnuc. Using a Weighted Leaky Box Model for propagation in the galaxy, we obtain an excellent fit to these new Voyager observations and the much higher energy spectra up to 100 GeVnuc by using source spectra which are P-2.28, with the exponent independent of rigidity from low to high rigidities; along with a rigidity dependence of the diffusion path length which is P-0.5 at rigidities 1.00 GV, and possibly changing to P1.0 at lower rigidities.


title: Amorphic complexity

abstract: We introduce amorphic complexity as a new topological invariant that measures the complexity of dynamical systems in the regime of zero entropy. Its main purpose is to detect the very onset of disorder in the asymptotic behaviour. For instance, it gives positive value to Denjoy examples on the circle and Sturmian subshifts, while being zero for all isometries and Morse-Smale systems. After discussing basic properties and examples, we show that amorphic complexity and the underlying asymptotic separation numbers can be used to distinguish almost automorphic minimal systems from equicontinuous ones. For symbolic systems, amorphic complexity equals the box dimension of the associated Besicovitch space. In this context, we concentrate on regular Toeplitz flows and give a detailed description of the relation to the scaling behaviour of the densities of the p-skeletons. Finally, we take a look at strange non-chaotic attractors appearing in so-called pinched skew product systems. Continuous-time systems, more general group actions and the application to cut and project quasicrystals will be treated in subsequent work.


title: Investigations of the coordinates in Ptolemy's Geographike Hyphegesis Book 8

abstract: In Book 8 of his Geographike Hyphegesis Ptolemy gives coordinates for ca. 360 so-called noteworthy cities. These coordinates are the time difference to Alexandria, the length of the longest day, and partly the ecliptic distance from the summer solstice. The supposable original conversions between the coordinates in Book 8 and the geographical coordinates in the location catalogue of Books 2-7 including the underlying parameters and tabulations are here reconstructed. The results document the differences between the $\Omega$- and $\Xi$-recension. The known difference in the longitude of Alexandria underlying the conversion of the longitudes is examined more closely. For the ecliptic distances from the summer solstice of the $\Omega$-recension it is revealed that they were originally computed by means of a so far undiscovered approximate, linear conversion. Further it is shown that the lengths of the longest day could be based on a linear interpolation of the data in the Mathematike Syntaxis 2.6.


title: Symmetry: a bridge between nature and culture

abstract: Are symmetries discovered or rather invented by humans ? The stand you may take firmly here reveals a lot of your epistemological position. Conversely, the arguments you may forge for answering to this question, or to one of its numerous narrower or broader variations, shape your whole philosophical thoughts; not specifically about science, by the way. I will try to show how physics helps to (re)consider this issue.


title: Indications of r-mode Oscillations in SOHO/MDI Solar Radius Measurements

abstract: Analysis of solar radius measurements acquired by the Michelson Doppler Imager on the SOHO spacecraft supports previously reported evidence of solar internal r-mode oscillations in Mt Wilson radius data and in nuclear-decay data acquired at the Lomonosov Moscow State University. The frequencies of these oscillations are compatible with oscillations in a putative inner tachocline that separates a slowly rotating core from the radiative envelope.


title: Possible role of gravity in collapse of the wave-function: a brief survey of some ideas

abstract: This article is a brief survey of some approaches to implementing the suggestion that collapse of the wave function is mediated by gravity. These approaches include: a possible connection between the problem of time and problem of quantum measurement, decoherence models based on space-time uncertainty, the Schr\"{o}dinger-Newton equation, attempts to introduce gravity into collapse models such as CSL, ideas based on the black hole - elementary particle complementarity, and the possible role of a complex space-time metric.


title: Numerical Conformal Mapping to One-Tooth Gear-Shaped Domains and Applications

abstract: We study conformal mappings from the unit disk (or a rectangle) to one-tooth gear-shaped planar domains from the point of view of the Schwarzian derivative, with emphasis on numerical considerations. Applications are given to evaluation of a singular integral, mapping to the complement of an annular rectangle, and symmetric multitooth domains.


title: Magnetic Fields of Uranus and Neptune: Metallic Fluid Hydrogen

abstract: The magnetic fields of the Ice Giant Planets Uranus and Neptune (U/N) are unique in the solar system. Based on a substantial database measured on Earth for representative planetary fluids at representative dynamic pressures up to 200 GPa (2 Mbar) and a few 1000 K, the complex magnetic fields of U/N are (i) probably made primarily by degenerate metallic fluid H (MFH) at or near the crossover from the H-He envelopes to Ice cores at ~100 GPa (Mbar) pressures and normalized radii of ~90% of the radii of U/N; (ii) because those magnetic fields are made relatively close to the surfaces of U/N, non-dipolar fields can be expected; (iii) the Ice cores are most probably a heterogeneous fluid mixture of H, N, O, C, Fe/Ni and silicate-oxides and their mutual reaction products at high pressures and temperatures,as discussed elsewhere. Ironically, there is probably little nebular Ice in the Ice Giant Planets.


title: Varieties of elementary subalgebras of maximal dimension for modular Lie algebras

abstract: Motivated by questions in modular representation theory, Carlson, Friedlander, and the first author introduced the varieties E(r, g) of r-dimensional abelian p-nilpotent subalgebras of a p-restricted Lie algebra g. In this paper, we identify the varieties E(r, g) for a reductive restricted Lie algebra g and r the maximal dimension of an abelian p-nilpotent subalgebra of g.


title: Resolution of SU(2) monopole singularities by oxidation

abstract: We show how "colored" SU(2) BPS monopoles (that is: SU(2) monopoles satisfying the Bogomol'nyi equation whose Higgs field and magnetic charge vanish at infinity and which are singular at the origin) can be obtained from the BPST instanton by a singular dimensional reduction, explaining the origin of the singularity and implying that the singularity can be cured by the oxidation of the solution. We study the oxidation of other monopole solutions in this scheme.


title: Unsteady self-sustained detonation waves in flake aluminum dust/air mixtures

abstract: Self-sustained detonation waves in flake aluminum dust/air mixtures have been studied in a tube of diameter 199 mm and length 32.4 m. A pressure sensor array of 32 sensors mounted around certain circumferences of the tube was used to measure the shape of the detonation front in the circumferential direction and pressure histories of the detonation wave. A two-head spin detonation wave front was observed for the aluminum dust/air mixtures, and the cellular structure resulting from the spinning movement of the triple point was analyzed. The variations in velocity and overpressure of the detonation wave with propagation distance in a cell were studied. The interactions of waves in triple-point configurations were analyzed and the flow-field parameters were calculated. Three types of triple-point configuration exist in the wave front of the detonation wave of an aluminum dust/air mixture. Both strong and weak transverse waves exist in the unstable self-sustained detonation wave.


title: Single-cell protein dynamics reproduce universal fluctuations in cell populations

abstract: Protein variability in single cells has been studied extensively in populations, but little is known about temporal protein fluctuations in a single cell over extended times. We present here traces of protein copy number measured in individual bacteria over multiple generations and investigate their statistical properties, comparing them to previously measured population snapshots. We find that temporal fluctuations in individual traces exhibit the same universal features as those previously observed in populations. Scaled fluctuations around the mean of each trace exhibit the same universal distribution shape as found in populations measured under a wide range of conditions and in two distinct microorganisms. Additionally, the mean and variance of the traces over time obey the same quadratic relation. Analyzing the temporal features of the protein traces in individual cells, reveals that within a cell cycle protein content increases as an exponential function with a rate that varies from cycle to cycle. This leads to a compact description of the protein trace as a 3-variable stochastic process - the exponential rate, the cell-cycle duration and the value at the cycle start - sampled once each cell cycle. This compact description is sufficient to preserve the universal statistical properties of the protein fluctuations, namely, the protein distribution shape and the quadratic relationship between variance and mean. Our results show that the protein distribution shape is insensitive to sub-cycle intracellular microscopic details and reflects global cellular properties that fluctuate between generations.


title: Generalized ballistic deposition in 2 dimensions : scaling of surface width, porosity and conductivity

abstract: A deposition process with particles having realistic intermediate stickiness is studied in 2+1 dimensions. At each stage of the deposition process, for any given configuration, a newly depositing particle gives rise to allowed set of configurations that are vastly larger than those for deposition of a mixture of purely non-sticky (random like) and purely sticky (ballistic like) particles. We obtain scaling behavior and demonstrate collapse of scaled data for surface width and porosity. Scaling of conductivity, when a porous structure thus formed, is saturated with conductive fluid, e.g. brine, is studied. The results obtained are in good agreement with Archie's law for porous sedimentary rocks.


title: 2-Swappability and the Edge-Reconstruction Number of Regular Graphs

abstract: The edge-reconstruction number of graph $G$, denoted $ern(G)$,is the size of the smallest multiset of edge-deleted, unlabeled subgraphs of $G$, from which the structure of $G$ can be uniquely determined. That there was some connection between the areas of edge reconstruction and swappability has been known since the swapping number of a graph was first introduced by Froncek, Rosenberg, and Hlavacek in 2013. This paper illustrates the depth of that connection by proving several bridging results between those areas; in particular, when the graphs in question are both regular and 2-swappable. These connections led to the discovery of four infinite families of $r\geq 3$ regular graphs with $ern(G) \geq 3$, contradicting the formerly conjectured upper bound.


title: Charge dynamics and spin blockade in a hybrid double quantum dot in silicon

abstract: Electron spin qubits in silicon, whether in quantum dots or in donor atoms, have long been considered attractive qubits for the implementation of a quantum computer due to the semiconductor vacuum character of silicon and its compatibility with the microelectronics industry. While donor electron spins in silicon provide extremely long coherence times and access to the nuclear spin via the hyperfine interaction, quantum dots have the complementary advantages of fast electrical operations, tunability and scalability. Here we present an approach to a novel hybrid double quantum dot by coupling a donor to a lithographically patterned artificial atom. Using gate-based rf reflectometry, we probe the charge stability of this double quantum dot system and the variation of quantum capacitance at the interdot charge transition. Using microwave spectroscopy, we find a tunnel coupling of 2.7 GHz and characterise the charge dynamics, which reveals a charge T2* of 200 ps and a relaxation time T1 of 100 ns. Additionally, we demonstrate spin blockade at the inderdot transition, opening up the possibility to operate this coupled system as a singlet-triplet qubit or to transfer a coherent spin state between the quantum dot and the donor electron and nucleus.


title: Parametric oscillator in a Kerr medium: evolution of coherent states

abstract: We study the temporal evolution of a coherent state under the action of a parametric oscillator and a nonlinear Kerr-like medium. We make use of the interaction picture representation and use an exact time evolution operator for the time independent part of the Hamiltonian. We approximate the interaction picture Hamiltonian in such a way as to make it a member of a Lie algebra. The corresponding time evolution operator behaves like a squeezing operator due to the temporal dependence of the oscillator's frequency. We analyze the probability amplitude and the auto correlation function for different Hamiltonian parameters and we find a very good agreement between our approximate results and converged numerical calculations.


title: Combining Probabilistic, Causal, and Normative Reasoning in CP-logic

abstract: In recent years the search for a proper formal definition of actual causation -- i.e., the relation of cause-effect as it is instantiated in specific observations, rather than general causal relations -- has taken on impressive proportions. In part this is due to the insight that this concept plays a fundamental role in many different fields, such as legal theory, engineering, medicine, ethics, etc. Because of this diversity in applications, some researchers have shifted focus from a single idealized definition towards a more pragmatic, context-based account. For instance, recent work by Halpern and Hitchcock draws on empirical research regarding people's causal judgments, to suggest a graded and context-sensitive notion of causation. Although we sympathize with many of their observations, their restriction to a merely qualitative ordering runs into trouble for more complex examples. Therefore we aim to improve on their approach, by using the formal language of CP-logic (Causal Probabilistic logic), and the framework for defining actual causation that was developed by the current authors using it. First we rephrase their ideas into our quantitative, probabilistic setting, after which we modify it to accommodate a greater class of examples. Further, we introduce a formal distinction between statistical and normative considerations.


title: Estimating the Benefits of Electric Vehicle Smart Charging at Non-Residential Locations: A Data-Driven Approach

abstract: In this paper, we use data collected from over 2000 non-residential electric vehicle supply equipments (EVSEs) located in Northern California for the year of 2013 to estimate the potential benefits of smart electric vehicle (EV) charging. We develop a smart charging framework to identify the benefits of non-residential EV charging to the load aggregators and the distribution grid. Using this extensive dataset, we aim to improve upon past studies focusing on the benefits of smart EV charging by relaxing the assumptions made in these studies regarding: (i) driving patterns, driver behavior and driver types; (ii) the scalability of a limited number of simulated vehicles to represent different load aggregation points in the power system with different customer characteristics; and (iii) the charging profile of EVs. First, we study the benefits of EV aggregations behind-the-meter, where a time-of-use pricing schema is used to understand the benefits to the owner when EV aggregations shift load from high cost periods to lower cost periods. For the year of 2013, we show a reduction of up to 24.8% in the monthly bill is possible. Then, following a similar aggregation strategy, we show that EV aggregations decrease their contribution to the system peak load by approximately 40% when charging is controlled within arrival and departure times. Our results also show that it could be expected to shift approximately 0.25kWh (~2.8%) of energy per non-residential EV charging session from peak periods (12PM-6PM) to off-peak periods (after 6PM) in Northern California for the year of 2013.


title: Increment entropy as a measure of complexity for time series

abstract: Entropy has been a common index to quantify the complexity of time series in a variety of fields. Here, we introduce increment entropy to measure the complexity of time series in which each increment is mapped into a word of two letters, one letter corresponding to direction and the other corresponding to magnitude. The Shannon entropy of the words is termed as increment entropy (IncrEn). Simulations on synthetic data and tests on epileptic EEG signals have demonstrated its ability of detecting the abrupt change, regardless of energetic (e.g. spikes or bursts) or structural changes. The computation of IncrEn does not make any assumption on time series and it can be applicable to arbitrary real-world data.


title: High temperature limits for $(1+1)$-dimensional directed polymer with heavy-tailed disorder

abstract: The directed polymer model at intermediate disorder regime was introduced by Alberts-Khanin-Quastel~\cite{AKQ12}. It was proved that at inverse temperature $\beta n^{-\gamma}$ with $\gamma=1/4$ the partition function, centered appropriately, converges in distribution and the limit is given in terms of the solution of the stochastic heat equation. This result was obtained under the assumption that the disorder variables posses exponential moments, but its universality was also conjectured under the assumption of six moments. We show that this conjecture is valid and we further extend it by exhibiting classes of different universal limiting behaviors in the case of less than six moments. We also explain the behavior of the scaling exponent for the log-partition function under different moment assumptions and values of $\gamma$.


title: The Bernstein-Sato $b$-function of the Vandermonde determinant

abstract: The Bernstein-Sato polynomial, or the $b$-function, is an important invariant of singularities of hypersurfaces that is difficult to compute in general. We describe a few different results towards computing the $b$-function of the Vandermonde determinant $\xi$. We use a result of Opdam to produce a lower bound for the $b$-function of $\xi$. This bound proves a conjecture of Budur, Musta\c{t}\u{a}, and Teitler for the case of finite Coxeter hyperplane arrangements, proving the Strong Monodromy Conjecture in this case. In our second set of results, we show the duality of two $\mathcal{D}$-modules, and conclude that the roots of the $b$-function of $\xi$ are symmetric about $-1$. We then use some results about jumping coefficients to prove an upper bound for the $b$-function of $\xi$, and finally we conjecture a formula for the $b$-function of $\xi$.


title: Secrecy Transmit Beamforming for Heterogeneous Networks

abstract: In this paper, we pioneer the study of physical-layer security in heterogeneous networks (HetNets). We investigate secure communications in a two-tier downlink HetNet, which comprises one macrocell and several femtocells. Each cell has multiple users and an eavesdropper attempts to wiretap the intended macrocell user. Firstly, we consider an orthogonal spectrum allocation strategy to eliminate co-channel interference, and propose the secrecy transmit beamforming only operating in the macrocell (STB-OM) as a partial solution for secure communication in HetNet. Next, we consider a secrecy-oriented non-orthogonal spectrum allocation strategy and propose two cooperative STBs which rely on the collaboration amongst the macrocell base station (MBS) and the adjacent femtocell base stations (FBSs). Our first cooperative STB is the STB sequentially operating in the macrocell and femtocells (STB-SMF), where the cooperative FBSs individually design their STB matrices and then feed their performance metrics to the MBS for guiding the STB in the macrocell. Aiming to improve the performance of STB-SMF, we further propose the STB jointly designed in the macrocell and femtocells (STB-JMF), where all cooperative FBSs feed channel state information to the MBS for designing the joint STB. Unlike conventional STBs conceived for broadcasting or interference channels, the three proposed STB schemes all entail relatively sophisticated optimizations due to QoS constraints of the legitimate users. In order to efficiently use these STB schemes, the original optimization problems are reformulated and convex optimization techniques, such as second-order cone programming and semidefinite programming, are invoked to obtain the optimal solutions. Numerical results demonstrate that the proposed STB schemes are highly effective in improving the secrecy rate performance of HetNet.


title: Kernel Interpolation for Scalable Structured Gaussian Processes (KISS-GP)

abstract: We introduce a new structured kernel interpolation (SKI) framework, which generalises and unifies inducing point methods for scalable Gaussian processes (GPs). SKI methods produce kernel approximations for fast computations through kernel interpolation. The SKI framework clarifies how the quality of an inducing point approach depends on the number of inducing (aka interpolation) points, interpolation strategy, and GP covariance kernel. SKI also provides a mechanism to create new scalable kernel methods, through choosing different kernel interpolation strategies. Using SKI, with local cubic kernel interpolation, we introduce KISS-GP, which is 1) more scalable than inducing point alternatives, 2) naturally enables Kronecker and Toeplitz algebra for substantial additional gains in scalability, without requiring any grid data, and 3) can be used for fast and expressive kernel learning. KISS-GP costs O(n) time and storage for GP inference. We evaluate KISS-GP for kernel matrix approximation, kernel learning, and natural sound modelling.


title: An algorithm for multiplication of split-octonions

abstract: In this paper we introduce efficient algorithm for the multiplication of split-octonions. The direct multiplication of two split-octonions requires 64 real multiplications and 56 real additions. More effective solutions still do not exist. We show how to compute a product of the split-octonions with 28 real multiplications and 92 real additions. During synthesis of the discussed algorithm we use the fact that product of two split-octonions may be represented as vector-matrix product. The matrix that participates in the product calculating has unique structural properties that allow performing its advantageous decomposition. Namely this decomposition leads to significant reducing of the multiplicative complexity of split-octonions multiplication.


title: Neutrino Oscillation Studies with Reactors

abstract: Nuclear reactors are one of the most intense, pure, controllable, cost-effective, and well-understood sources of neutrinos. Reactors have played a major role in the study of neutrino oscillations, a phenomenon that indicates that neutrinos have mass and that neutrino flavors are quantum mechanical mixtures. Over the past several decades reactors were used in the discovery of neutrinos, were crucial in solving the solar neutrino puzzle, and allowed the determination of the smallest mixing angle $\theta_{13}$. In the near future, reactors will help to determine the neutrino mass hierarchy and to solve the puzzling issue of sterile neutrinos.


title: Constraints on the off-shell Higgs boson signal strength in the high-mass $ZZ$ and $WW$ final states with the ATLAS detector

abstract: Measurements of the $ZZ$ and $WW$ final states in the mass range above the $2m_Z$ and $2m_W$ thresholds provide a unique opportunity to measure the off-shell coupling strength of the Higgs boson. This paper presents constraints on the off-shell Higgs boson event yields normalised to the Standard Model prediction (signal strength) in the $ZZ \to 4\ell$, $ZZ\to2\ell2\nu$ and $WW\to e\nu\mu\nu$ final states. The result is based on $pp$ collision data collected by the ATLAS experiment at the LHC, corresponding to an integrated luminosity of $20.3$ fb$^\text{-1}$ at a collision energy of $\sqrt{s}$=8 TeV. Using the $CL_s$ method, the observed 95% confidence level (CL) upper limit on the off-shell signal strength is in the range 5.1-8.6, with an expected range of 6.7-11.0. In each case the range is determined by varying the unknown $gg\to ZZ$ and $gg\to WW$ background K-factor from higher-order QCD corrections between half and twice the value of the known signal K-factor. Assuming the relevant Higgs boson couplings are independent of the energy scale of the Higgs production, a combination with the on-shell measurements yields an observed (expected) 95\% CL upper limit on $\Gamma_H/\Gamma_H^{\mathrm{SM}}$ in the range 4.5- -7.5 (6.5- -11.2) using the same variations of the background K-factor. Assuming that the unknown $gg\rightarrow VV$ background K-factor is equal to the signal K-factor,this translates into an observed (expected) 95% CL upper limit on the Higgs boson total width of 22.7 (33.0) MeV.


title: Distributed Hierarchical Control versus an Economic Model for Cloud Resource Management

abstract: We investigate a hierarchically organized cloud infrastructure and compare distributed hierarchical control based on resource monitoring with market mechanisms for resource management. The latter do not require a model of the system, incur a low overhead, are robust, and satisfy several other desiderates of autonomic computing. We introduce several performance measures and report on simulation studies which show that a straightforward bidding scheme supports an effective admission control mechanism, while reducing the communication complexity by several orders of magnitude and also increasing the acceptance rate compared to hierarchical control and monitoring mechanisms. Resource management based on market-based mechanisms can be seen as an intermediate step towards cloud self-organization, an ideal alternative to current mechanisms for cloud resource management.


title: Broadband detuned Sagnac interferometer for future generation gravitational wave astronomy

abstract: Broadband suppression of quantum noise below the Standard Quantum Limit (SQL) becomes a top-priority problem for the future generation of large-scale terrestrial detectors of gravitational waves, as the interferometers of the Advanced LIGO project, predesigned to be quantum-noise-limited in the almost entire detection band, are phased in. To this end, among various proposed methods of quantum noise suppression or signal amplification, the most elaborate approach implies a so-called *xylophone* configuration of two Michelson interferometers, each optimised for its own frequency band, with a combined broadband sensitivity well below the SQL. Albeit ingenious, it is a rather costly solution. We demonstrate that changing the optical scheme to a Sagnac interferometer with weak detuned signal recycling and frequency dependent input squeezing can do almost as good a job, as the xylophone for significantly lower spend. We also show that the Sagnac interferometer is more robust to optical loss in filter cavity, used for frequency dependent squeezed vacuum injection, than an analogous Michelson interferometer, thereby reducing building cost even more.


title: On the Capacity of Wireless Networks with Random Transmission Delay

abstract: In this paper, we introduce novel coding schemes for wireless networks with random transmission delays. These coding schemes obviate the need for synchronicity, reduce the number of transmissions and achieve the optimal rate region in the corresponding wired model for both multiple unicast and multicast cases with up to three users under the equal rate constraint. The coding schemes are presented in two phases; first, coding schemes for line, star and line-star topologies with random transmission delays are provided. Second, any general topology with multiple bidirectional unicast and multicast sessions is shown to be decomposable into these canonical topologies to reduce the number of transmissions without rate redundancy. As a result, the coding schemes developed for the line, star and line-star topologies serve as building blocks for the construction of more general coding schemes for all networks. The proposed schemes are proved to be Real Time (RT) for wireless networks in the sense that they achieve the minimal decoding delay. With a negligible size header, these coding schemes are shown to be applicable to unsynchronized networks, i.e., networks with random transmission delays. Finally, we demonstrate the applicability of these schemes by extensive simulations. The implementation of such coding schemes on a wireless network with random transmission delay can improve performance and power efficiency.


title: Existence and uniqueness of the global solution to the Navier-Stokes boundary problem

abstract: A proof is given of the global existence and uniqueness of a weak solution to Navier-Stokes boundary problem. The proof is short and essentially self-contained.


title: Collaboration Tools and Patterns for Creative Thinking

abstract: Many creativity methods follow similar structures and principles. Design Patterns capture such invariants of proven good practices and discuss why, when and how creative thinking methods match various situations of collaboration. Moreover patterns connect different forms with each other. Once we understand the underlying structures of creative thinking processes we can facilitate digital tools to support them. While such tools can foster the effective application of established methods and even change their properties, tools can also enable new patterns of collaboration.


title: Collective achievement of making in cosplay culture

abstract: This paper analyzes peer-based learning and the concept of Scaffolding represented in ethnographic case studies of ten female informants aged 20-25 participating in the cosplay community. Cosplay is a female-dominated niche subculture of extreme fans and mavens, who are devoted to dressing up as characters from manga, games, and anime. Cosplayers are highly conscious of quality standards for costumes, makeup, and accessories. Cosplay events and dedicated SNSs for cosplayers are a valuable venue for exchanging information about costume making. First we frame this work as an effort to think about their learning environment using the concept of connected learning by Ito et al(2013). Then we share an overview of cosplay culture in Japan and our methodologies based on interviews and fieldwork. The interview transcripts were analyzed according to the Steps for coding and Theorization method, a qualitative data analysis technique by Otani (2008). In this study, We focus on their reciprocal learning and expanding the concept of Scaffolding.


title: Exploring Cultures through Pattern Mining - Practices from Generative Beauty Workshops

abstract: This paper presents a method for understanding personal ways of thinking and doing in daily lives among different countries by mining their ways as patterns in a sense of pattern language. Pattern language is a methodology of describing tacit practical knowledge, where each pattern consists of context, problem, and solution. In this paper, patterns mined from the workshops we held in the following three countries: Japan, Korea, and the United States, are analysed. The results demonstrate similarities and reflect characteristics of the patterns of each country. We anticipate that this workshop can be used as a method for better understanding of cultural similarities and features in the light of practical knowledge in daily lives.


title: An approach to computing downward closures

abstract: The downward closure of a word language is the set of all (not necessarily contiguous) subwords of its members. It is well-known that the downward closure of any language is regular. While the downward closure appears to be a powerful abstraction, algorithms for computing a finite automaton for the downward closure of a given language have been established only for few language classes. This work presents a simple general method for computing downward closures. For language classes that are closed under rational transductions, it is shown that the computation of downward closures can be reduced to checking a certain unboundedness property. This result is used to prove that downward closures are computable for (i) every language class with effectively semilinear Parikh images that are closed under rational transductions, (ii) matrix languages, and (iii) indexed languages (equivalently, languages accepted by higher-order pushdown automata of order 2).


title: The spectral index of signed Laplacians and their structural stability

abstract: Given a graph Laplacian with positively and negatively weighted edges we are interested in characterizing the set of weights that give a particular spectral index, i.e.~give a prescribed number of positive, zero, and negative eigenvalues. One of the main results of this paper is that the set of signed Laplacians that exhibit multiple zero eigenvalues is "small", and that eigenvalue crossings are nongeneric --- specifically, eigenvalues repel each other near zero in a sense that can be made precise. We exhibit an algebraic discriminant that measures the level of repulsion, and show that this discriminant admits a combinatorial interpretation. Conversely, we exhibit a constructive method for finding the sets of Laplacians that exhibit a large degree of degeneracy (many eigenvalues at or near zero) in terms of these discriminants.


title: Using Descriptive Video Services to Create a Large Data Source for Video Annotation Research

abstract: In this work, we introduce a dataset of video annotated with high quality natural language phrases describing the visual content in a given segment of time. Our dataset is based on the Descriptive Video Service (DVS) that is now encoded on many digital media products such as DVDs. DVS is an audio narration describing the visual elements and actions in a movie for the visually impaired. It is temporally aligned with the movie and mixed with the original movie soundtrack. We describe an automatic DVS segmentation and alignment method for movies, that enables us to scale up the collection of a DVS-derived dataset with minimal human intervention. Using this method, we have collected the largest DVS-derived dataset for video description of which we are aware. Our dataset currently includes over 84.6 hours of paired video/sentences from 92 DVDs and is growing.


title: Asymptotic optimality of Tailored Base-Surge policies in dual-sourcing inventory systems

abstract: Dual-sourcing inventory systems, in which one supplier is faster (i.e. express) and more costly, while the other is slower (i.e. regular) and cheaper, arise naturally in many real-world supply chains. These systems are notoriously difficult to optimize due to the complex structure of the optimal solution and the curse of dimensionality, having resisted solution for over 40 years. Recently, so-called Tailored Base-Surge (TBS) policies have been proposed as a heuristic for the dual-sourcing problem. Under such a policy, a constant order is placed at the regular source in each period, while the order placed at the express source follows a simple order-up-to rule. Numerical experiments by several authors have suggested that such policies perform well as the lead time difference between the two sources grows large, which is exactly the setting in which the curse of dimensionality leads to the problem becoming intractable. However, providing a theoretical foundation for this phenomenon has remained a major open problem. In this paper, we provide such a theoretical foundation by proving that a simple TBS policy is indeed asymptotically optimal as the lead time of the regular source grows large, with the lead time of the express source held fixed. Our main proof technique combines novel convexity and lower-bounding arguments, an explicit implementation of the vanishing discount factor approach to analyzing infinite-horizon Markov decision processes, and ideas from the theory of random walks and queues, significantly extending the methodology and applicability of a novel framework for analyzing inventory models with large lead times recently introduced by Goldberg and co-authors in the context of lost-sales models with positive lead times.


title: Frobenius-Schur indicators for some fusion categories associated to symmetric and alternating groups

abstract: We calculate Frobenius-Schur indicator values for some fusion categories obtained from inclusions of finite groups $H\subset G$, where more concretely $G$ is symmetric or alternating, and $H$ is a symmetric, alternating or cyclic group. Our work is strongly related to earlier results by Kashina-Mason-Montgomery, Jedwab-Montgomery, and Timmer for bismash product Hopf algebras obtained from exact factorizations of groups. We can generalize some of their results, settle some open questions and offer shorter proofs; this already pertains to the Hopf algebra case, while our results also cover fusion categories not associated to Hopf algebras.


title: T3PS: Tool for Parallel Processing in Parameter Scans

abstract: T3PS is a program that can be used to quickly design and perform parameter scans while easily taking advantage of the multi-core architecture of current processors. It takes an easy to read and write parameter scan definition file format as input. Based on the parameter ranges and other options contained therein, it distributes the calculation of the parameter space over multiple processes and possibly computers. The derived data is saved in a plain text file format readable by most plotting software. The supported scanning strategies include: grid scan, random scan, Markov Chain Monte Carlo, numerical optimization. Several example parameter scans are shown and compared with results in the literature.


title: Fundamental limitations to high-precision tests of the universality of free fall by dropping atoms

abstract: Tests of the universality of free fall and the weak equivalence principle probe the foundations of General Relativity. Evidence of a violation may lead to the discovery of a new force. The best torsion balance experiments have ruled it out to 10^-13. Cold-atom drop tests have reached 10^-7 and promise to do 7 to 10 orders of magnitude better, on the ground or in space. They are limited by the random shot noise, which depends on the number N of atoms in the clouds. As mass-dropping experiments in the non-uniform gravitational field of Earth, they are sensitive to the initial conditions. Random accelerations due to initial condition errors of the clouds are designed to be at the same level as shot noise, so that they can be reduced with the number of drops along with it. This sets the requirements for the initial position and velocity spreads of the clouds with given N. In the STE-QUEST space mission proposal aiming at 2x10^-15 they must be about a factor 8 above Heisenberg's principle limit, and the integration time required to reduce both errors is 3 years, with a mission duration of 5 years. Instead, offset errors at release between different atom clouds are systematic and give rise to a systematic effect which mimics a violation. Such offsets must be demonstrated to be as small as required in all drops, must be small by design and must be measured. For STE-QUEST to meet its goal they must be several orders of magnitude smaller than the size of each individual cloud, which in its turn must be at most 8 times larger than the uncertainty principle limit. Even if all technical problems are solved and the clouds are released with negligible systematic errors, still they must be measured. Then, Heisenberg's principle dictates that the measurement lasts as long as the experiment and the systematic nature of the effect requires many measurements for it to be ruled out as a source of violation.


title: Two Interesting Properties of the Exponential Distribution

abstract: Let $X_1, X_2,\ldots, X_n$ be $n$ independent and identically distributed random variables, here $n \geq 2.$ Let $X_{(1)}, X_{(2)}, \ldots, X_{(n)}$ be the order statistics of $X_1, X_2,..., X_n.$ In this note we proved that: (I) If $X_1, X_2,..., X_n$ are exponential random variables with parameter $c > 0,$ then the "correlation coefficient" between $X_{(k)}$ and $X_{(k+t)}$ is strictly increasing in $k$ from $1$ to $m,$ and then is strictly decreasing in $k$ from $m$ to $n - t,$ here $t$ is a fixed integer between $1$ and $n - 3,$ and $m = (n - t)/2$ if $n - t$ is even, $m = (n - t + 1)/2$ if $n - t$ is odd. We also proved that if $t = n - 2$, then the "correlation coefficient" between $X_{(1)}$ and $X_{(n-1)}$ is greater than the "correlation coefficient" between $X_{(2)}$ and$X_{(n)}.$ (II) The "correlation coefficient" between $X_{(k)}$ and $X_{(k+t)}$ for the exponential random variables is always less than the "correlation coefficient" between $X_{(k)}$ and $X_{(k+t)}$ for the uniform random variables for all $k$ and $t$ such that $k + t \leq n.$ A combinatorial identity is also given as a bi-product. \vs


title: Trapping of classical particles by an electromagnetic radiation amplifying with time

abstract: We analyze new possible applications of the trapping mechanism of sufficiently slow-speed particles by an electromagnetic potential well deepening with time (up to a certain limit) which was recently established by author from basic relations of classical mechanics. It is assumed that given particles are under conditions of the high vacuum and forces acting on these particles are not dissipative. Such wells may be created by means of an electromagnetic field (in particular radiation) with a fixed spatial distribution but with a nondecreasing strength in time. Trapping and localization of particles in such electromagnetic traps are analyzed on example of typical laser beams. Obtained results may be used in high resolution spectroscopy of various particles (including also atoms and molecules in definite cases).


title: Variational Approach to Projected Entangled Pair States at Finite Temperature

abstract: The projected entangled pair state (PEPS) ansatz can represent a thermal state in a strongly correlated system. We introduce a novel variational algorithm to optimize this tensor network. Since full tensor environment is taken into account, then with increasing bond dimension the optimized PEPS becomes the exact Gibbs state. Our presentation opens with a 1D version for a matrix product state (MPS) and then generalizes to PEPS in 2D. Benchmark results in the quantum Ising model are presented.


title: $\mu\to e \gamma$ Decay in an MSSM Extension

abstract: An analysis is given of the decay $\mu \to e+ \gamma$ in an MSSM extension with a vectorlike generation. Here mixing with the mirrors allows the possibility of this decay. The analysis is done at the one loop level with the exchange of charginos and neutralinos and of sleptons and mirror sleptons in the loops. A one loop analysis with W and Z boson exchange and mirror leptons and neutrinos is also considered. The effects of CP violating phases from the new sector on the decay $\mu\to e \gamma$ are analyzed in detail. The constraints arising from the current upper limit on the branching ratio ${\cal{B}}(\mu\to {e} \gamma)$ from the MEG experiment of $2.4\times 10^{-12}$ (at 90\% CL) on the parameter space of SUSY models and on vectorlike models are explored. Further, the MEG experiment is likely to improve the upper limit by an order of magnitude in the coming years. The improved limits will allow one to probe a much larger domain of the parameter space of the extended models.


title: Duality results in Banach and quasi-Banach spaces of homogeneous polynomials and applications

abstract: Spaces of homogeneous polynomials on a Banach space are frequently equipped with quasinorms instead of norms. In this paper we develop a technique to replace the original quasi-norm by a norm in a dual preserving way, in the sense that the dual of the space with the new norm coincides with the dual of the space with the original quasi-norm. Applications to problems on the existence and approximation of solutions of convolution equations and on hypercyclic convolution operators on spaces of entire functions are provided.


title: A Generalization of Wantzel's Theorem, m-sectable angles, and the density of certain Chebyshev-polynomial images

abstract: The eponymous theorem of P.L. Wantzel presents a necessary and sufficient criterion for angle trisectability in terms of the third Chebyshev polynomial $T_3$, thus making it easy to prove that there exist non-trisectable angles. We generalize this theorem to the case of all Chebyshev polynomials $T_m$ . We also study the set \textbf{m-Sect} consisting of all cosines of $m$-sectable angles (see \S 1), showing that, when $m$ is not a power of two, \textbf{m-Sect} contains only algebraic numbers . We then introduce a notion of density based on the diophantine-geometric concept of height of an algebraic number and obtain a result on the density of certain polynomial images. Using this in conjunction with the Generalized Wantzel Theorem, we obtain our main result: for every real algebraic number field $K$, the set \textbf{m-Sect}\ $\cap\ K$ has density zero in $[-1,1]\ \cap\ K$ when $m$ is not a power of two. (To appear in the Journal of Pure and Applied Algebra.)


title: Genome-wide modelling of transcription kinetics reveals patterns of RNA production delays

abstract: Genes with similar transcriptional activation kinetics can display very different temporal mRNA profiles due to differences in transcription time, degradation rate and RNA processing kinetics. Recent studies have shown that a splicing-associated RNA production delay can be significant. We introduce a joint model of transcriptional activation and mRNA accumulation which can be used for inference of transcription rate, RNA production delay and degradation rate given genome-wide data from high-throughput sequencing time course experiments. We combine a mechanistic differential equation model with a non-parametric statistical modelling approach allowing us to capture a broad range of activation kinetics, and use Bayesian parameter estimation to quantify the uncertainty in the estimates of the kinetic parameters. We apply the model to data from estrogen receptor (ER-{\alpha}) activation in the MCF-7 breast cancer cell line. We use RNA polymerase II (pol-II) ChIP-Seq time course data to characterise transcriptional activation and mRNA-Seq time course data to quantify mature transcripts. We find that 11% of genes with a good signal in the data display a delay of more than 20 minutes between completing transcription and mature mRNA production. The genes displaying these long delays are significantly more likely to be short. We also find a statistical association between high delay and late intron retention in pre-mRNA data, indicating significant splicing-associated production delays in many genes.


title: Assessing a human mediated current awareness service

abstract: In this paper, we present an approach for analyzing the behavior of editors in the large current awareness service "NEP: New Economics Papers". We processed data from more than 38,000 issues derived from 90 different NEP reports over the past ten years. The aim of our analysis was to gain an inside to the editor behaviour when creating an issue and to look for factors that influence the success of a report. In our study we looked at the following features: average editing time, the average number of papers in an issue and the editor effort measured on presorted issues as relative search length (RSL). We found an average issue size of 12.4 documents per issue. The average editing time is rather low with 14.5 minute. We get to the point that the success of a report is mainly driven by its topic and the number of subscribers, as well as proactive action by the editor to promote the report in her community.


title: A Performance Estimator for Quantum Annealers: Gauge selection and Parameter Setting

abstract: With the advent of large-scale quantum annealing devices, several challenges have emerged. For example, it has been shown that the performance of a device can be significantly affected by several degrees of freedom when programming the device; a common example being gauge selection. To date, no experimentally-tested strategy exists to select the best programming specifications. We developed a score function that can be calculated from a number of readouts much smaller than the number of readouts required to find the desired solution. We show how this performance estimator can be used to guide, for example, the selection of the optimal gauges out of a pool of random gauge candidates and how to select the values of parameters for which we have no a priori knowledge of the optimal value. For the latter, we illustrate the concept by applying the score function to set the strength of the parameter intended to enforce the embedding of the logical graph into the hardware architecture, a challenge frequently encountered in the implementation of real-world problem instances. Since the harder the problem instances, the more useful the strategies proposed in this work are, we expect the programming strategies proposed to significantly reduce the time of future benchmark studies and in help finding the solution of hard-to-solve real-world applications implemented in the next generation of quantum annealing devices.


title: Clues for flavor from rare lepton and quark decays

abstract: Flavor symmetries successfully explain lepton and quark masses and mixings yet it is usually hard to distinguish different models that predict the same mixing angles. Further experimental input could be available, if the agents of flavor breaking are sufficiently low in mass and detectable or if new physics with non-trivial flavor charges is sufficiently low in mass and detectable. The recent hint for lepton-nonuniversality in the ratio of branching fractions $B \to K \mu \mu$ over $B \to K e e$, $R_K$, suggests the latter, at least for indirect detection via rare decays. We demonstrate the discriminating power of the rare decay data on flavor model building taking into account viable leptonic mixings and show how correlations with other observables exist in leptoquark models. We give expectations for branching ratios $B \to K \ell \ell^\prime, B_{(s)} \to \ell \ell^\prime$ and $\ell \to \ell^\prime \gamma$, and Higgs decays $h \to \ell \ell^\prime$.


title: Environmental impact of Supernova Remnants

abstract: The explosion of a supernovae (SN) represents the sudden injection of about 10^51 ergs of thermal and mechanical energy in a small region of space, causing the formation of powerful shock waves that propagate through the interstellar medium at speeds of several thousands of km/s. These waves sweep, compress and heat the interstellar material that they encounter, forming the supernova remnants. Their evolution over thousands of years change forever, irreversibly, not only the physical but also the chemical properties of a vast region of space that can span hundreds of parsecs. This contribution briefly analyzes the impact of these explosions, discussing the relevance of some phenomena usually associated with SNe and their remnants in the light of recent theoretical and observational results.


title: Some properties of a sequence defined with the aid of prime numbers

abstract: For every integer $n\ge 1$ let $a_n$ be the smallest positive integer such that $n+a_n$ is prime. We investigate the behavior of the sequence $(a_n)_{n\ge 1}$, and prove asymptotic results for the sums $\sum_{n\le x} a_n$, $\sum_{n\le x} 1/a_n$ and $\sum_{n\le x} \log a_n$.


title: Fundamental properties of High-Mass X-ray Binaries

abstract: The aim of this PhD Thesis is to characterize a representative sample of Supergiant X-ray Binaries (SGXBs) formed by 4 sources: XTE J1855-026, a classical SGXB with long-term stable X-ray flux; AX J1841.0-0535 and AX J1845.0-0433, two supergiant fast X-ray transients (SFXTs) with the X-ray emission mostly dominated by flaring; and IGR J00370+6122, something in between these 2 sub-groups. The physical processes that produce these observable differences are still a matter of debate. In this PhD Thesis I performed a study of these 4 different systems to provide new data to constrain the models. This study consists of:(i) the determination of the orbital solution,(ii) a systematic study of the wind behavior along the orbit by the measure of Halpha variations,(iii) a model of stellar atmospheres of the donor star,(iv) establish whether there are X-ray flux variations modulated by the orbital period. The study of the wind shows that Halpha variations are dominated by intrinsic wind processes. The stellar atmospheres study shows that the supergiant stars that harbor these binaries have a higher projected rotational velocity, higher He abundance and higher N/C ratio than that of isolated supergiant stars. The results of this study show that the eccentricity of the binary does not have a simple correlation with the differences in the X-ray flux of the different sub-groups. Therefore, the idea of a more complex scenario is consolidated. The discovery of a new type of system, IGR J00370+6122, which properties do not fit in any of the established sub-groups, reinforces the idea of a continuum in the observed properties more than a strict classification of systems. Furthermore, I have developed a pipeline to reduce spectra of the FRODOSpec spectrograph at the Liverpool Telescope optimized for the reduction of the red spectra of the obscured supergiants that we find in these SGXBs.


title: On jet substructure methods for signal jets

abstract: We carry out simple analytical calculations and Monte Carlo studies to better understand the impact of QCD radiation on some well-known jet substructure methods for jets arising from the decay of boosted Higgs bosons. Understanding differences between taggers for these signal jets assumes particular significance in situations where they perform similarly on QCD background jets. As an explicit example of this we compare the Y-splitter method to the more recently proposed Y-pruning technique. We demonstrate how the insight we gain can be used to significantly improve the performance of Y-splitter by combining it with trimming and show that this combination outperforms the other taggers studied here, at high $p_T$. We also make analytical estimates for optimal parameter values, for a range of methods and compare to results from Monte Carlo studies.


title: Dark Matter Balls Help Supernovae to Explode

abstract: As a solution to the well-known problem that the shock wave potentially responsible for the explosion of a supernova actually tends to stall, we propose a new energy source arising from our model for dark matter. Our earlier model proposed that dark matter should consist of cm-large white dwarf-like objects kept together by a skin separating two different sorts of vacua. These dark matter balls or pearls will collect in the middle of any star throughout its lifetime. At some stage during the development of a supernova the balls will begin to take in neutrons and then other surrounding material. By passing into a ball nucleons fall through a potential of order 10 MeV, causing a severe production of heat - of order 10 foe for a solar mass of material eaten by the balls. The temperature in the iron core will thereby be raised, splitting up the iron into smaller nuclei. This provides a mechanism for reviving the shock wave when it arrives and making the supernova explosion really occur. The onset of the heating due to the dark matter balls would at first stop the collapse of the supernova progenitor. This opens up the possibility of there being {\em two} collapses giving two neutrino outbursts, as apparently seen in the supernova SN1987A - one in Mont Blanc, and one 4 hours 43 minutes later in both IMB and Kamiokande.


title: Super compact pairwise model for SIS epidemic on heterogeneous networks

abstract: In this paper we provide the derivation of a super compact pairwise model with only 4 equations in the context of describing susceptible-infected-susceptible (SIS) epidemic dynamics on heterogenous networks. The super compact model is based on a new closure relation that involves not only the average degree but also the second and third moments of the degree distribution. Its derivation uses an a priori approximation of the degree distribution of susceptible nodes in terms of the degree distribution of the network. The new closure gives excellent agreement with heterogeneous pairwise models that contain significantly more differential equations.


title: Lower bounds for interior nodal sets of Steklov eigenfunctions

abstract: We study the interior nodal sets, $Z_\lambda$ of Steklov eigenfunctions in an $n$-dimensional relatively compact manifolds $M$ with boundary and show that one has the lower bounds $|Z_\lambda|\ge c\lambda^{\frac{2-n}2}$ for the size of its $(n-1)$-dimensional Hausdorff measure. The proof is based on a Dong-type identity and estimates for the gradient of Steklov eigenfunctions, similar to those in \cite{SZ1} and \cite{SZ2}, respectively.


title: Never Judge a Black Hole by Its Area

abstract: Christodoulou and Rovelli have shown that black holes have large interiors that grow asymptotically linearly in advanced time, and speculated that this may be relevant to the information loss paradox. We show that there is no simple relation between the interior volume of an arbitrary black hole and its horizon area. That is, the volume enclosed is not necessarily a monotonically increasing function of the surface area.


title: A note on the longest common Abelian factor problem

abstract: Abelian string matching problems are becoming an object of considerable interest in last years. Very recently, Alatabbi et al. \cite{AILR2015} presented the first solution for the longest common Abelian factor problem for a pair of strings, reaching $O(\sigma n^2)$ time with $O(\sigma n \log n)$ bits of space, where $n$ is the length of the strings and $\sigma$ is the alphabet size. In this note we show how the time complexity can be preserved while the space is reduced by a factor of $\sigma$, and then how the time complexity can be improved, if the alphabet is not too small, when superlinear space is allowed.


title: Universal logarithmic corrections to entanglement entropies in two dimensions with spontaneously broken continuous symmetries

abstract: We explore the R\'enyi entanglement entropies of a one-dimensional (line) subsystem of length $L$ embedded in two-dimensional $L\times L$ square lattice for quantum spin models whose ground-state breaks a continuous symmetry in the thermodynamic limit. Using quantum Monte Carlo simulations, we first study the $J_1 - J_2$ Heisenberg model with antiferromagnetic nearest-neighbor $J_1>0$ and ferromagnetic second-neighbor couplings $J_2\le 0$. The signature of SU(2) symmetry breaking on finite size systems, ranging from $L=4$ up to $L=40$ clearly appears as a universal additive logarithmic correction to the R\'enyi entanglement entropies: $l_q \ln L$ with $l_q\simeq 1$, independent of the R\'enyi index and values of $J_2$. We confirm this result using a high precision spin-wave analysis (with restored spin rotational symmetry) on finite lattices up to $10^5\times 10^5$ sites, allowing to explore further non-universal finite size corrections and study in addition the case of U(1) symmetry breaking. Our results fully agree with the prediction $l_q=n_G/2$ where $n_G$ is the number of Goldstone modes, by Metlitski and Grover [arXiv:1112.5166].


title: libRoadRunner: A High Performance SBML Simulation and Analysis Library

abstract: This paper presents libRoadRunner, an extensible, high-performance, cross-platform, open-source software library for the simulation and analysis of models \ expressed using Systems Biology Markup Language (SBML). SBML is the most widely used standard for representing dynamic networks, especially biochemical networks. libRoadRunner supports solution of both large models and multiple replicas of a single model on desktop, mobile and cluster computers. libRoadRunner is a self-contained library, able to run both as a component inside other tools via its C++ and C bindings andnteractively through its Python interface. The Python Application Programming Interface (API) is similar to the APIs of Matlab and SciPy, making it fast and easy to learn, even for new users. libRoadRunner uses a custom Just-In-Time (JIT) compiler built on the widely-used LLVM JIT compiler framework to compile SBML-specified models directly into very fast native machine code for a variety of processors, making it appropriate for solving very large models or multiple replicas of smaller models. libRoadRunner is flexible, supporting the bulk of the SBML specification (except for delay and nonlinear algebraic equations) and several of its extensions. It offers multiple deterministic and stochastic integrators, as well as tools for steady-state, stability analyses and flux balance analysis. We regularly update libRoadRunner binary distributions for Mac OS X, Linux and Windows and license them under Apache License Version 2.0. http://www.libroadrunner.org provides online documentation, full build instructions, binaries and a git source repository.


title: Similarity Is Not Enough: Tipping Points of Ebola Zaire Mortalities

abstract: In early 2014 an outbreak of a slightly mutated Zaire Ebola subtype appeared in West Africa which is less virulent than 1976 and 1994 strains. The numbers of cases per year appear to be ~ 1000 times larger than the earlier strains, suggesting a greatly enhanced transmissibility. Although the fraction of the 2014 spike glycoprotein mutations is very small (~3%), the mortality is significantly reduced, while the transmission appears to have increased strongly. Bioinformatic scaling had previously shown similar inversely correlated trends in virulence and transmission in N1 (H1N1) and N2 (H3N2) influenza spike glycoprotein mutations. These trends appear to be related to various external factors (migration, availability of pure water, and vaccination programs). The molecular mechanisms for Ebola's mutational response involve mainly changes in the disordered mucin-like domain (MLD) of its spike glycoprotein amino acids. The MLD has been observed to form the tip of an oligomeric amphiphilic wedge that selectively pries apart cell-cell interfaces via an oxidative mechanism.


title: Negative differential resistance and characteristic nonlinear electromagnetic response of a Topological Insulator

abstract: Materials exhibiting negative differential resistance have important applications in technologies involving microwave generation, which range from motion sensing to radio astronomy. Despite their usefulness, there has been few physical mechanisms giving rise to materials with such properties, i.e. GaAs employed in the Gunn diode. In this work, we show that negative differential resistance also generically arise in Dirac ring systems, an example of which has been experimentally observed in the surface states of Topological Insulators. This novel realization of negative differential resistance is based on a completely different physical mechanism from that of the Gunn effect, relying on the characteristic non-monotonicity of the response curve that remains robust in the presence of nonzero temperature, chemical potential, mass gap and impurity scattering. As such, it opens up new possibilities for engineering applications, such as frequency upconversion devices which are highly sought for terahertz signal generation. Our results may be tested with thin films of Bi2Se3 Topological Insulators, and are expected to hold qualitatively even in the absence of a strictly linear Dirac dispersion, as will be the case in more generic samples of Bi2Se3 and other materials with topologically nontrivial Fermi sea regions.


title: Recognizing k-equistable graphs in FPT time

abstract: A graph $G = (V,E)$ is called equistable if there exist a positive integer $t$ and a weight function $w : V \to \mathbb{N}$ such that $S \subseteq V$ is a maximal stable set of $G$ if and only if $w(S) = t$. Such a function $w$ is called an equistable function of $G$. For a positive integer $k$, a graph $G = (V,E)$ is said to be $k$-equistable if it admits an equistable function which is bounded by $k$. We prove that the problem of recognizing $k$-equistable graphs is fixed parameter tractable when parameterized by $k$, affirmatively answering a question of Levit et al. In fact, the problem admits an $O(k^5)$-vertex kernel that can be computed in linear time.


title: Singularities of secant varieties

abstract: We study the singularities of the secant variety $\Sigma(X,L)$ associated to a smooth variety $X$ embedded by a sufficiently positive adjoint bundle $L$. We show that $\Sigma(X,L)$ is always Du Bois singular. Examples of secant varieties with worse singularities when $L$ has weak positivity are provided. We also give a necessary and sufficient condition for $\Sigma(X, L)$ to have rational singularities.


title: Multi-Agent Decentralized Network Interdiction Games

abstract: In this work, we introduce decentralized network interdiction games, which model the interactions among multiple interdictors with differing objectives operating on a common network. As a starting point, we focus on decentralized shortest path interdiction (DSPI) games, where multiple interdictors try to increase the shortest path lengths of their own adversaries, who all attempt to traverse a common network. We first establish results regarding the existence of equilibria for DSPI games under both discrete and continuous interdiction strategies. To compute such an equilibrium, we present a reformulation of the DSPI games, which leads to a generalized Nash equilibrium problem (GNEP) with non-shared constraints. While such a problem is computationally challenging in general, we show that under continuous interdiction actions, a DSPI game can be formulated as a linear complementarity problem and solved by Lemke's algorithm. In addition, we present decentralized heuristic algorithms based on best response dynamics for games under both continuous and discrete interdiction strategies. Finally, we establish theoretical bounds on the worst-case efficiency loss of equilibria in DSPI games, with such loss caused by the lack of coordination among noncooperative interdictors, and use the decentralized algorithms to empirically study the average-case efficiency loss.


title: Antideuterons from Decaying Gravitino Dark Matter

abstract: We study the possibility of improving the constraints on the lifetime of gravitino dark matter in scenarios with bilinear R-parity violation by estimating the amount of cosmic-ray antideuterons that can be produced in gravitino decays. Taking into account all different sources of theoretical uncertainties, we find that the margin of improvement beyond the limits already set by cosmic-ray antiproton data are quite narrow and unachievable for the next generation of experiments. However, we also identify more promising energy ranges for future experiments.


title: Cooperative Base Station Coloring for Pair-wise Multi-Cell Coordination

abstract: This paper proposes a method for designing BS clusters and cluster patterns for pair-wise BS coordination. The key idea is that each BS cluster is formed by using the 2nd-order Voronoi region, and the BS clusters are assigned to a specific cluster pattern by using edge-coloring for a graph drawn by Delaunay triangulation. The main advantage of the proposed method is that the BS selection conflict problem is prevented, while selected users are guaranteed to communicate with their two closest BSs in any irregular BS topology. With the proposed coordination method, analytical expressions for the rate distribution and the ergodic spectral efficiency are derived as a function of relevant system parameters in a fixed irregular network model. In a random network model with a homogeneous Poisson point process, a lower bound on the ergodic spectral efficiency is characterized. Through system level simulations, the performance of the proposed method is compared with that of conventional coordination methods: dynamic clustering and static clustering. Our major finding is that, when users are dense enough in a network, the proposed method provides the same level of coordination benefit with dynamic clustering to edge users.


title: On uniform estimates for Laplace equation in balls with small holes

abstract: In this paper, we consider the Dirichlet problem of the three-dimensional Laplace equation in the unit ball with a shrinking hole. The problem typically arises from homogenization problems in domains perforated with tiny holes. We give an almost complete description concerning the uniform $W^{1,p}$ estimates: for any $3/2<p<3$ there hold the uniform $W^{1,p}$ estimates; for any $1<p<3/2$ or $3<p<\infty $, there are counterexamples indicating that the uniform $W^{1,p}$ estimates do not hold. The results can be generalized to higher dimensions.


title: You Can Run, You Can Hide: The Epidemiology and Statistical Mechanics of Zombies

abstract: We use a popular fictional disease, zombies, in order to introduce techniques used in modern epidemiology modelling, and ideas and techniques used in the numerical study of critical phenomena. We consider variants of zombie models, from fully connected continuous time dynamics to a full scale exact stochastic dynamic simulation of a zombie outbreak on the continental United States. Along the way, we offer a closed form analytical expression for the fully connected differential equation, and demonstrate that the single person per site two dimensional square lattice version of zombies lies in the percolation universality class. We end with a quantitative study of the full scale US outbreak, including the average susceptibility of different geographical regions.


title: ROSA: Robust sparse adaptive channel estimation in the presence of impulsive noises

abstract: Based on the assumption of Gaussian noise model, conventional adaptive filtering algorithms for reconstruction sparse channels were proposed to take advantage of channel sparsity due to the fact that broadband wireless channels usually have the sparse nature. However, state-of-the-art algorithms are vulnerable to deteriorate under the assumption of non-Gaussian noise models (e.g., impulsive noise) which often exist in many advanced communications systems. In this paper, we study the problem of RObust Sparse Adaptive channel estimation (ROSA) in the environment of impulsive noises using variable step-size affine projection sign algorithm (VSS-APSA). Specifically, standard VSS-APSA algorithm is briefly reviewed and three sparse VSS-APSA algorithms are proposed to take advantage of channel sparsity with different sparse constraints. To fairly evaluate the performance of these proposed algorithms, alpha-stable noise is considered to approximately model the realistic impulsive noise environments. Simulation results show that the proposed algorithms can achieve better performance than standard VSS-APSA algorithm in different impulsive environments.


title: Test of a fine pitch SOI pixel detector with laser beam

abstract: A silicon pixel detector with fine pitch size of 19x19 um, developed base on SOI (silicon on insulator) technology, was tested under the illumination of infrared laser pulses. As an alternative way to particle beam tests, the laser pulses were tuned to very short duration and small transverse profile to simulate the tracks of MIPs (minimum ionization particles) in silicon. Hit cluster sizes were measured with focused laser pulses propagating through the SOI detector perpendicular to its surface and most of the induced charge was found to be collected inside the seed pixel. For the first time, the signal amplitude as a function of the applied bias voltage was measured for this SOI detector, deepening understanding of its depletion characteristics.


title: Hot Gaseous Coronae around Spiral Galaxies: Probing the Illustris Simulation

abstract: The presence of hot gaseous coronae around present-day massive spiral galaxies is a fundamental prediction of galaxy formation models. However, our observational knowledge remains scarce, since to date only four gaseous coronae were detected around spirals with massive stellar bodies ($\gtrsim2\times10^{11} \ \rm{M_{\odot}}$). To explore the hot coronae around lower mass spiral galaxies, we utilized Chandra X-ray observations of a sample of eight normal spiral galaxies with stellar masses of $(0.7-2.0)\times10^{11} \ \rm{M_{\odot}}$. Although statistically significant diffuse X-ray emission is not detected beyond the optical radii ($\sim20$ kpc) of the galaxies, we derive $3\sigma$ limits on the characteristics of the coronae. These limits, complemented with previous detections of NGC 1961 and NGC 6753, are used to probe the Illustris Simulation. The observed $3\sigma$ upper limits on the X-ray luminosities and gas masses exceed or are at the upper end of the model predictions. For NGC 1961 and NGC 6753 the observed gas temperatures, metal abundances, and electron density profiles broadly agree with those predicted by Illustris. These results hint that the physics modules of Illustris are broadly consistent with the observed properties of hot coronae around spiral galaxies. However, a shortcoming of Illustris is that massive black holes, mostly residing in giant ellipticals, give rise to powerful radio-mode AGN feedback, which results in under luminous coronae for ellipticals.


title: The 21-SPONGE HI Absorption Survey I: Techniques and Initial Results

abstract: We present methods and results from "21-cm Spectral Line Observations of Neutral Gas with the EVLA" (21-SPONGE), a large survey for Galactic neutral hydrogen (HI) absorption with the Karl G. Jansky Very Large Array (VLA). With the upgraded capabilities of the VLA, we reach median root-mean-square (RMS) noise in optical depth of $\sigma_{\tau}=9\times 10^{-4}$ per $0.42\rm\,km\,s^{-1}$ channel for the 31 sources presented here. Upon completion, 21-SPONGE will be the largest HI absorption survey with this high sensitivity. We discuss the observations and data reduction strategies, as well as line fitting techniques. We prove that the VLA bandpass is stable enough to detect broad, shallow lines associated with warm HI, and show that bandpass observations can be combined in time to reduce spectral noise. In combination with matching HI emission profiles from the Arecibo Observatory ($\sim3.5'$ angular resolution), we estimate excitation (or spin) temperatures ($\rm T_s$) and column densities for Gaussian components fitted to sightlines along which we detect HI absorption (30/31). We measure temperatures up to $\rm T_s\sim1500\rm\,K$ for individual lines, showing that we can probe the thermally unstable interstellar medium (ISM) directly. However, we detect fewer of these thermally unstable components than expected from previous observational studies. We probe a wide range in column density between $\sim10^{16}$ and $>10^{21}\rm\,cm^{-2}$ for individual HI clouds. In addition, we reproduce the trend between cold gas fraction and average $\rm T_s$ found by synthetic observations of a hydrodynamic ISM simulation by Kim et al. (2014). Finally, we investigate methods for estimating HI $\rm T_s$ and discuss their biases.


title: Exploring the liminality: properties of haloes and subhaloes in borderline $f(R)$ gravity

abstract: We investigate the properties of dark matter haloes and subhaloes in an $f(R)$ gravity model with $|f_{R0}|=10^{-6}$, using a very high-resolution N-body simulation. The model is a borderline between being cosmologically interesting and yet still consistent with current data. We find that the halo mass function in this model has a maximum 20% enhancement compared with the $\Lambda$CDM predictions between $z=1$ and $z=0$. Because of the chameleon mechanism which screens the deviation from standard gravity in dense environments, haloes more massive than $10^{13}h^{-1}M_\odot$ in this $f(R)$ model have very similar properties to haloes of similar mass in $\Lambda$CDM, while less massive haloes, such as that of the Milky Way, can have steeper inner density profiles and higher velocity dispersions due to their weaker screening. The halo concentration is remarkably enhanced for low-mass haloes in this model due to a deepening of the total gravitational potential. Contrary to the naive expectation, the halo formation time $z_f$ is later for low-mass haloes in this model, a consequence of these haloes growing faster than their counterparts in $\Lambda$CDM at late times and the definition of $z_f$. Subhaloes, especially those less massive than $10^{11}h^{-1}M_\odot$, are substantially more abundant in this $f(R)$ model for host haloes less massive than $10^{13}h^{-1}M_\odot$. We discuss the implications of these results for the Milky Way satellite abundance problem. Although the overall halo and subhalo properties in this borderline $f(R)$ model are close to their $\Lambda$CDM predictions, our results suggest that studies of the Local Group and astrophysical systems, aided by high-resolution simulations, can be valuable for further tests of it.


title: Probing The Ultraviolet Luminosity Function of the Earliest Galaxies with the Renaissance Simulations

abstract: In this paper, we present the first results from the Renaissance Simulations, a suite of extremely high-resolution and physics-rich AMR calculations of high redshift galaxy formation performed on the Blue Waters supercomputer. These simulations contain hundreds of well-resolved galaxies at $z \sim 25-8$, and make several novel, testable predictions. Most critically, we show that the ultraviolet luminosity function of our simulated galaxies is consistent with observations of high-z galaxy populations at the bright end of the luminosity function (M$_{1600} \leq -17$), but at lower luminosities is essentially flat rather than rising steeply, as has been inferred by Schechter function fits to high-z observations, and has a clearly-defined lower limit in UV luminosity. This behavior of the luminosity function is due to two factors: (i) the strong dependence of the star formation rate on halo virial mass in our simulated galaxy population, with lower-mass halos having systematically lower star formation rates and thus lower UV luminosities; and (ii) the fact that halos with virial masses below $\simeq 2 \times 10^8$ M$_\odot$ do not universally contain stars, with the fraction of halos containing stars dropping to zero at $\simeq 7 \times 10^6$ M$_\odot$. Finally, we show that the brightest of our simulated galaxies may be visible to current and future ultra-deep space-based surveys, particularly if lensed regions are chosen for observation.