title: Mutual emergence of noncausal optical response and nonclassicality in an optomechanical system

abstract: We show that single-mode nonclassicality of the output of an optomechanical cavity and the noncausal linear optical response of this cavity emerge at the same critical cavity-mechanical coupling. In other words, single-mode nonclassicality emerges when the barrier (in electromagnetism) avoiding faster-than-light communication is lifted off. The nature of the emergence of noncausal behavior does not depend on the length (boundary conditions) and the type of the cavity. Origin of the noncausal behavior is the temporal/frequency relations between the incident and reflected waves at the outer surface of the cavity. We further discuss the relations with the recent studies; (i) equivalence of the entanglement among identical particles to the nonclassicality of their quasiparticle excitations, (ii) necessity of superfluid behavior of vacuum, and (iii) entanglement-wormhole equivalence.


title: The set chromatic number of random graphs

abstract: In this paper we study the set chromatic number of a random graph $G(n,p)$ for a wide range of $p=p(n)$. We show that the set chromatic number, as a function of $p$, forms an intriguing zigzag shape.


title: (Almost) Dark HI Sources in the ALFALFA Survey: The Intriguing Case of HI1232+20

abstract: We report the discovery and follow-up observations of a system of three objects identified by the ALFALFA extragalactic HI survey, cataloged as (almost) dark extragalactic sources, i.e., extragalactic HI detections with no discernible counterpart in publicly available, wide-field, imaging surveys. We have obtained deep optical imaging with WIYN pODI and HI synthesis maps with WSRT of the HI1232+20 system. The source with the highest HI flux has a newly discovered ultra-low surface brightness (LSB) optical counterpart associated with it, while the other two sources have no detected optical counterparts in our images. Our optical observations show that the detected LSB optical counterpart has a peak surface brightness of ~26.4 mag/arcsec^2 in g', which is exceptionally faint. This source (AGC 229385) has the largest accurately measured HI mass-to-light ratio of an isolated object: MHI/Lg'=46 Msun/Lsun, and has an HI mass of 7.2*10^8 Msun. The other two HI sources (with HI masses 2.0*10^8 and 1.2*10^8 Msun) without optical counterparts have upper limit surface brightnesses of 27.9 and 27.8 mag/arcsec^2 in g', and lower limits on their gas mass-to-light ratio of MHI/Lg'>57 and >31 Msun/Lsun. This system lies relatively close in projection to the Virgo Cluster, but velocity flow models indicate that it is located at ~25 Mpc, substantially beyond Virgo. The system appears to be quite isolated, with no known object closer than 500 kpc. These HI sources may represent both sides of the threshold between "dark" star-less galaxies and galaxies with stellar populations. We discuss a variety of possible formation scenarios for the HI1232+20 system.


title: The equitable presentation of $\mathfrak{osp}_q(1|2)$ and a $q$-analog of the Bannai-Ito algebra

abstract: The equitable presentation of the quantum superalgebra $\mathfrak{osp}_q(1|2)$, in which all generators appear on an equal footing, is exhibited. It is observed that in their equitable presentations, the quantum algebras $\mathfrak{osp}_q(1|2)$ and $\mathfrak{sl}_q(2)$ are related to one another by the formal transformation $q\rightarrow -q$. A $q$-analog of the Bannai-Ito algebra is shown to arise as the covariance algebra of $\mathfrak{osp}_q(1|2)$.


title: Assessing the direction of climate interactions by means of complex networks and information theoretic tools

abstract: An estimate of the net direction of climate interactions in different geographical regions is made by constructing a directed climate network from a regular latitude-longitude grid of nodes, using a directionality index (DI) based on conditional mutual information. Two datasets of surface air temperature anomalies - one monthly-averaged and another daily-averaged - are analyzed and compared. The network links are interpreted in terms of known atmospheric tropical and extratropical variability patterns. Specific and relevant geographical regions are selected, the net direction of propagation of the atmospheric patterns is analyzed and the direction of the inferred links is validated by recovering some well-known climate variability structures. These patterns are found to be acting at various time-scales, such as atmospheric waves in the extra-tropics or longer range events in the tropics. This analysis demonstrates the capability of the DI measure to infer the net direction of climate interactions and may contribute to improve the present understanding of climate phenomena and climate predictability. The work presented here also stands out as an application of advanced tools to the analysis of empirical, real-world data.


title: Compound Chondrules fused Cold

abstract: About 4-5% of chondrules are compound: two separate chondrules stuck together. This is commonly believed to be the result of the two component chondrules having collided shortly after forming, while still molten. This allows high velocity impacts to result in sticking. However, at T ~ 1100K, the temperature below which chondrules collide as solids (and hence usually bounce), coalescence times for droplets of appropriate composition are measured in tens of seconds. Even at 1025K, at which temperature theory predicts that the chondrules must have collided extremely slowly to have stuck together, the coalescence time scale is still less than an hour. These coalescence time scales are too short for the collision of molten chondrules to explain the observed frequency of compound chondrules. We suggest instead a scenario where chondrules stuck together in slow collisions while fully solid; and the resulting chondrule pair was subsequently briefly heated to a temperature in the range of 900-1025K. In that temperature window the coalescence time is finite but long, covering a span of hours to a decade. This is particularly interesting because those temperatures are precisely the critical window for thermally ionized MRI activity, so compound chondrules provide a possible probe into that vital regime.


title: Towards a physics-based modelling of the electro-mechanical coupling in EAPs

abstract: Due to the increasing number of industrial applications of electro-active polymers (EAPs), there is a growing need for electromechanical models which accurately capture their behavior. To this end, we compare the predicted behavior of EAPs undergoing homogenous deformations according to three electromechanical models. The first model is a continuum based model composed of the mechanical Gent model and a linear relationship between the electric field and the polarization. The electrical and the mechanical responses according to the second model are based on the polymer microstructure, whereas the third model incorporates a neo-Hookean mechanical response and a microstructural based long-chains model for the electrical behavior. In the microstructural motivated models the integration from the microscopic to the macroscopic levels is accomplished by the micro-sphere technique. Four types of homogeneous boundary conditions are considered and the behaviors determined according to the three models are compared. The differences between the predictions of the models are discussed, highlighting the need for an in-depth investigation of the relations between the structure and the behaviors of the EAPs at microscopic level and their overall macroscopic response.


title: Sliding friction in the hydrodynamic lubrication regime for a power-law fluid

abstract: A scaling analysis is undertaken for the load balance in sliding friction in the hydrodynamic lubrication regime, with a particular emphasis on power-law shear-thinning typical of a structured liquid. It is argued that the shear-thinning regime is mechanically unstable if the power-law index n < 1/2, where n is the exponent that relates the shear stress to the shear rate. Consequently the Stribeck (friction) curve should be discontinuous, with possible hysteresis. Further analysis suggests that normal stress and flow transience (stress overshoot) do not destroy this basic picture, although they may provide stabilising mechanisms at higher shear rates. Extensional viscosity is also expected to be insignificant unless the Trouton ratio is large. A possible application to shear thickening in non-Brownian particulate suspensions is indicated.


title: Numerical Estimation of Frictional Torques with Rate and State Friction

abstract: In this paper, numerical estimation of frictional torques is carried out of a rotary elastic disc on a hard and rough surface under different rotating conditions. A one dimensional spring- mass rotary system is numerically solved under the quasistatic condition with the rate and state dependent friction model. It is established that torque of frictional strength as well as torque of steady dynamic stress increases with radius and found to be maximum at the periphery of the disc. Torque corresponding to frictional strength estimated using the analytical solution matches closely with the simulation only in the case of high stiffness of the connecting spring. In steady relaxation simulation, a steadily rotating disc is suddenly stopped and relaxational angular velocity and corresponding frictional torque decreases with both steady angular velocity and stiffness of the connecting spring in the velocity strengthening regime. In velocity weakening regime, in contrast, torque of relaxation stress deceases but relaxation velocity increases. The reason for the contradiction is explained.


title: A category mistake in observational claims regarding ultrashort-lived unstable particles

abstract: The physics literature contains many claims that ultrashort-lived unstable particles, such as a Higgs boson, have been observed. These claims are a matter of applying the $5\sigma$-convention in particle physics. This paper, however, shows that by applying this $5\sigma$-convention a category mistake is made, by which a pure reasoning is passed off as an observation. Not only are these two fundamentally different primitive notions at the very basis of science, but the pure reasoning in question is also weaker than an observation: what we have in each case is that the existence of the ultrashort-lived unstable particle is inferred to the best explanation, but that does absolutely not merit the stronger claim that the particle in question has been "observed". Consequently, the observational claims in question will thus have to be dismissed as overstatements. On a general note, this demonstrates that the empirical support for the Standard Model of particle physics is significantly less than hitherto thought.


title: Constraining the Physics of AM Canum Venaticorum Systems with the Accretion Disk Instability Model

abstract: Recent work by Levitan et al has expanded the long-term photometric database for AM CVn stars. In particular, their outburst properties are well-correlated with orbital period, and allow constraints to be placed on the secular mass transfer rate between secondary and primary if one adopts the disk instability model for the outbursts. We use the observed range of outbursting behavior for AM CVn systems as a function of orbital period to place a constraint on mass transfer rate versus orbital period P. We infer a rate ~5 x 10^{-9} Msun/yr (P/1000 s)^{-5.2}. We show the functional form so obtained is consistent with the recurrence time-orbital period relation found by Levitan et al using a simple theory for the recurrence time. Also, we predict their steep dependence of outburst duration on orbital period will flatten considerably once the longer orbital period systems have more complete observations.


title: Phase transitions of regular Schwarzschild-Anti-deSitter black holes

abstract: We study a solution of the Einstein's equations generated by a self-gravitating, anisotropic, static, non-singular matter fluid. The resulting Schwarzschild like solution is regular and accounts for smearing effects of noncommutative fluctuations of the geometry. We call this solution regular Schwarzschild spacetime. In the presence of an Anti-deSitter cosmological term, the regularized metric offers an extension of the Hawking-Page transition into a van der Waals-like phase diagram. Specifically the regular Schwarzschild-Anti-deSitter geometry undergoes a first order small/large black hole transition similar to the liquid/gas transition of a real fluid. In the present analysis we have considered the cosmological constant as a dynamical quantity and its variation is included in the first law of black hole thermodynamics.


title: Percolation on the stationary distributions of the voter model

abstract: The voter model on $\mathbb{Z}^d$ is a particle system that serves as a rough model for changes of opinions among social agents or, alternatively, competition between biological species occupying space. When $d \geq 3$, the set of (extremal) stationary distributions is a family of measures $\mu_\alpha$, for $\alpha$ between 0 and 1. A configuration sampled from $\mu_\alpha$ is a strongly correlated field of 0's and 1's on $\mathbb{Z}^d$ in which the density of 1's is $\alpha$. We consider such a configuration as a site percolation model on $\mathbb{Z}^d$. We prove that if $d \geq 5$, the probability of existence of an infinite percolation cluster of 1's exhibits a phase transition in $\alpha$. If the voter model is allowed to have sufficiently spread-out interactions, we prove the same result for $d \geq 3$.


title: A Comparison of Radial Intensity Profiles of Termination Shock Particles and Anomalous Cosmic Rays in the Outer North-South Heliosheaths Using CRS data from V1 and V2

abstract: Using data from the Voyager 1 and 2 CRS telescopes available on the web through the end of 2014 we have studied the intensity variations of termination shock H nuclei and anomalous cosmic ray H and He nuclei as a function of radial distance. In contrast to the inner part of the heliosheath where the intensity vs. radius profiles in the North and South heliosheaths are much different, these intensity vs. radius profiles, as well as the intensities themselves, are more similar in the outer North and South heliosheaths as measured by V1 and V2 respectively. In the N heliosheath, taken to be 27.6 AU thick beyond the HTS crossing distance of 94 AU, the intensities of termination shock particles and anomalous cosmic rays reach a maximum at between 110-112 AU or at a location ~halfway between the termination shock and the heliopause. They then decrease more or less continuously to an intensity ~2/3 of the maximum for each component just before the final dropout at ~121 AU, just inside the heliopause. This intensity-time behavior favors theories in which the acceleration and loss effects lead to a maximum near the middle of the heliosheath.


title: Hybridization and the origin of the Fano resonances in symmetric nanoparticle trimers

abstract: We study the light scattering by plasmonic and dielectric symmetric trimers to investigate the existence of polarization-independent Fano resonances. Plasmonic hybridization theory is revealed to hide simple physics, and we instead provide a simplified model for hybridization to derive a plasmonic trimer's eigenmodes analytically. This approach is demonstrated to accurately recreate full wave simulations of plasmonic trimers and their Fano resonances. We are subsequently able to deduce the grounds for modal interference in plasmonic trimers and the related formation of Fano resonances. However, by taking advantage of the generality of our simplified hybridization approach, we are also able to investigate the eigenmodes of all-dielectric trimers. We show that bianisotropic coupling channels between high-index dielectric nanoparticles are able to increase the capacity for Fano resonances, even at normal incidence. We finally provide the first experimental measurements of sharp, polarization-independent Fano resonances from a symmetric all-dielectric trimer, with very good agreement to the predicted response from our simplified hybridization theory.


title: Entropy determination based on the ordinal structure of a dynamical system

abstract: The ordinal approach to evaluate time series due to innovative works of Bandt and Pompe has increasingly established itself among other techniques of nonlinear time series analysis. In this paper, we summarize and generalize the theory of determining the Kolmogorov-Sinai entropy of a measure-preserving dynamical system via increasing sequences of order generated partitions of the state space. Our main focus are measuring processes without information loss. Particularly, we consider the question of the minimal necessary number of measurements related to the properties of a given dynamical system.


title: Correlations and synchronization in a Bose-Fermi mixture

abstract: We study a Bose-Fermi mixture within the framework of the mean-field theory, including three possible regimes for the fermionic species: fully polarized, BCS, and unitarity. Starting from the 3D description and using the variational approximation (VA), we derive 1D and 2D systems of equations, under the corresponding confining potentials. This method produces a pair of nonlinear Schr\"{o}dinger (NLS) equations coupled to algebraic equations for the transverse widths of the confined state. The equations incorporate interactions between atoms of the same species and between the species, assuming that the latter can be manipulated by means of the Feshbach resonance (FR). As an application, we explore spatial density correlations in the ground state (GS) between the species, concluding that they strongly depend on the sign and strength of the inter-species interaction. Also studied are the dynamics of the mixture in a vicinity of the GS and the corresponding spatiotemporal inter-species correlation. The correlations are strongly affected by the fermionic component, featuring the greatest variation in the unitary regime. Results produced by the VA are verified by comparison with full numerical solutions.


title: The effect of reconnection on the structure of the Sun's open-closed-flux boundary

abstract: Global magnetic field extrapolations are now revealing the huge complexity of the Sun's corona, and in particular the structure of the boundary between open and closed magnetic flux. Moreover, recent developments indicate that magnetic reconnection in the corona likely occurs in highly fragmented current layers, and that this typically leads to a dramatic increase in the topological complexity beyond that of the equilibrium field. In this paper we investigate the consequences of reconnection at the open-closed flux boundary ("interchange reconnection") in a fragmented current layer. We demonstrate that it leads to a situation in which magnetic flux (and therefore plasma) from open and closed field regions is efficiently mixed together. This corresponds to an increase in the length and complexity of the open-closed boundary. Thus, whenever reconnection occurs at a null point or separator of the open-closed boundary, the associated separatrix arc of the so-called "S-web" in the high corona becomes not a single line but a band of finite thickness within which the open-closed flux boundary is highly structured. This has significant implications for the acceleration of the slow solar wind, for which the interaction of open and closed field is thought to be important, and may also explain the coronal origins of certain solar energetic particles. The topological structures examined contain magnetic null points, separatrices and separators, and include a model for a pseudo-streamer. The potential for understanding both the large scale morphology and fine structure observed in flare ribbons associated with coronal nulls is also discussed.


title: Vivace: a collaborative live coding language and platform

abstract: Live coding is a performance and creative technique based on improvised and interactive coding. Many recent endeavors have focused in live coding both because of aesthetics and as a way to alleviate performance drawbacks when the musical instrument is a computer. This paper describes the principles and the design of Vivace, a live coding language and environment built with Web technologies to be executed on web browsers. The approach is compelling by 1) allowing many performers to code simultaneously, 2) the synthesis of audio and video, 3) a very simple syntax, 4) being a multiplatform software. We also strive to contextualize Vivace by means of historical and usage summaries including a live coding sub-genre.


title: Wedge-local fields in integrable models with bound states

abstract: Recently, large families of two-dimensional quantum field theories with factorizing S-matrices have been constructed by the operator-algebraic methods, by first showing the existence of observables localized in wedge-shaped regions. However, these constructions have been limited to the class of S-matrices whose components are analytic in rapidity in the physical strip. In this work, we construct candidates for observables in wedges for scalar factorizing S-matrices with poles in the physical strip and show that they weakly commute on a certain domain. We discuss some technical issues concerning further developments, especially the self-adjointness of the candidate operators here and strong commutativity between them.


title: Twisted index theory on orbifold symmetric products and the fractional quantum Hall effect

abstract: We extend the noncommutative geometry model of the fractional quantum Hall effect, previously developed by Mathai and the first author, to orbifold symmetric products. It retains the same properties of quantization of the Hall conductance at integer multiples of the fractional Satake orbifold Euler characteristics. We show that it also allows for interesting composite fermions and anyon representations, and possibly for Laughlin type wave functions.


title: Photonic Nanojets in Optical Tweezers

abstract: Photonic nanojets have been brought into attention ten years ago for potential application in ultramicroscopy, because of its sub-wavelength resolution that can enhance detection and interaction with matter. For these novel applications under development, the optical trapping of a sphere acts as an ideal framework to employ photonic nanojets. In the present study, we generated nanojets by using a highly focused incident beam, in contrast to traditional plane waves. The method inherits the advantage of optical trapping, especially for intracellular applications, with the microsphere in equilibrium on the beam propagation axis and positioned arbitrarily in space. Moreover, owing to optical scattering forces, when the sphere is in equilibrium, its center shifts with respect to the focal point of the incident beam. However, when the system is in stable equilibrium with a configuration involving optical tweezers, photonic nanojets cannot be formed. To overcome this issue, we employed double optical tweezers in an unorthodox configuration involving two collinear and co-propagating beams, the precise positioning of which would turn on/off the photonic nanojets, thereby improving the applicability of photonic nanojets.


title: Gradient systems on coupled cell networks

abstract: For networks of coupled dynamical systems we characterize admissible functions, that is, functions whose gradient is an admissible vector field. The schematic representation of a gradient network dynamical system is of an undirected cell graph, and we use tools from graph theory to deduce the general form of such functions, relating it to the topological structure of the graph defining the network. The coupling of pairs of dynamical systems cells is represented by edges of the graph, and from spectral graph theory we detect the existence and nature of equilibria of the gradient system from the critical points of the coupling function. In particular, we study fully synchronous and 2-state patterns of equilibria on regular graphs.These are two special types of equilibrium configurations for gradient networks. We also investigate equilibrium configurations of S1-invariant admissible functions on a ring of cells.


title: Euler characteristics of centralizer subcategories

abstract: Let p be a prime number, G a finite group, and A a finite group acting on G. The Brown poset of nonidentity p-subgroups of G is then an A-poset. We investigate the equivariant subposet and the equivariant Euler characteristics and establish a global relation between locally defined Euler characteristics and the number of p-elements of G centralized by A. It is a consequence of this relation that the equivariant version of Brown's theorem holds: The reduced Euler characteristic of the A-equivariant Brown poset is divisible by the p-part of the order of the centralizer of A. The second equivariant Euler characteristic for the conjugation of G on the Brown poset for G is especially intriguing because of its relation to the Knorr-Robinson conjecture and we carry out a concrete numerical verification of the conjecture in case of the smallest simple Mathieu group.


title: New RR Lyrae variables in binary systems

abstract: Despite their importance, very few RR Lyrae (RRL) stars have been known to reside in binary systems. We report on a search for binary RRL in the OGLE-III Galactic bulge data. Our approach consists in the search for evidence of the light-travel time effect in so-called observed minus calculated ($O-C$) diagrams. Analysis of 1952 well-observed fundamental-mode RRL in the OGLE-III data revealed an initial sample of 29 candidates. We used the recently released OGLE-IV data to extend the baselines up to 17 years, leading to a final sample of 12 firm binary candidates. We provide $O-C$ diagrams and binary parameters for this final sample, and also discuss the properties of 8 additional candidate binaries whose parameters cannot be firmly determined at present. We also estimate that $\gtrsim 4$ per cent of the RRL reside in binary systems.


title: Computation of Local Time of Reflecting Brownian Motion and Probabilistic Representation of the Neumann Problem

abstract: In this paper, we propose numerical methods for computing the boundary local time of reflecting Brownian motion (RBM) in R3 and its use in the probabilistic representation of the solution of the Laplace equation with the Neumann boundary condition. Approximations of the RBM based on a walk-on-spheres (WOS) and random walk on lattices are discussed and tested for sampling the RBM paths and their applicability in finding accurate approximation of the local time and discretization of the probabilistic formula. Numerical tests for several types of domains (cube, sphere, and ellipsoid) have shown the convergence of the numerical methods as the length of the RBM path and number of paths sampled increase.


title: TASI 2014 Lectures: The Hunt for Dark Matter

abstract: These lectures, given at the 2014 Theoretical Advanced Study Institute (TASI), are an introduction to what we know at present about dark matter and the major current experimental and observational efforts to identify what it consists of. They attempt to present the complexities of the subject, making clear common simplifying assumptions, to better understand the reach of dark matter searches.


title: Numerical Solution of Fuzzy Stochastic Differential Equation

abstract: In this paper an alternative approach to solve uncertain Stochastic Differential Equation (SDE) is proposed. This uncertainty occurs due to the involved parameters in system and these are considered as Triangular Fuzzy Numbers (TFN). Here the proposed fuzzy arithmetic in [2] is used as a tool to handle Fuzzy Stochastic Differential Equation (FSDE). In particular, a system of Ito stochastic differential equations is analysed with fuzzy parameters. Further exact and Euler Maruyama approximation methods with fuzzy values are demonstrated and solved some standard SDE.


title: Sensor Control for Multi-Object Tracking Using Labeled Multi-Bernoulli Filter

abstract: The recently developed labeled multi-Bernoulli (LMB) filter uses better approximations in its update step, compared to the unlabeled multi-Bernoulli filters, and more importantly, it provides us with not only the estimates for the number of targets and their states, but also with labels for existing tracks. This paper presents a novel sensor-control method to be used for optimal multi-target tracking within the LMB filter. The proposed method uses a task-driven cost function in which both the state estimation errors and cardinality estimation errors are taken into consideration. Simulation results demonstrate that the proposed method can successfully guide a mobile sensor in a challenging multi-target tracking scenario.


title: The Pierre Auger Cosmic Ray Observatory

abstract: The Pierre Auger Observatory, located on a vast, high plain in western Argentina, is the world's largest cosmic ray observatory. The objectives of the Observatory are to probe the origin and characteristics of cosmic rays above $10^{17}$ eV and to study the interactions of these, the most energetic particles observed in nature. The Auger design features an array of 1660 water-Cherenkov particle detector stations spread over 3000 km$^2$ overlooked by 24 air fluorescence telescopes. In addition, three high elevation fluorescence telescopes overlook a 23.5 km$^2$, 61-detector infilled array with 750 m spacing. The Observatory has been in successful operation since completion in 2008 and has recorded data from an exposure exceeding 40,000 km$^2$ sr yr. This paper describes the design and performance of the detectors, related subsystems and infrastructure that make up the Auger Observatory.


title: Subsystem real-time Time Dependent Density Functional Theory

abstract: We present the extension of Frozen Density Embedding (FDE) theory to real-time Time Dependent Density Functional Theory (rt-TDDFT). FDE a is DFT-in-DFT embedding method that allows to partition a larger Kohn-Sham system into a set of smaller, coupled Kohn-Sham systems. Additional to the computational advantage, FDE provides physical insight into the properties of embedded systems and the coupling interactions between them. The extension to rt-TDDFT is done straightforwardly by evolving the Kohn-Sham subsystems in time simultaneously, while updating the embedding potential between the systems at every time step. Two main applications are presented: the explicit excitation energy transfer in real time between subsystems is demonstrated for the case of the Na$_4$ cluster and the effect of the embedding on optical spectra of coupled chromophores. In particular, the importance of including the full dynamic response in the embedding potential is demonstrated.


title: Independent Manipulation of Electric and Thermal Fields with Bilayer Structure

abstract: Recently, increasing attention has been focused on the employment of transformation and metamaterial for manipulation of various physical fields, which requires complicated configuration and usually limits in single field. Here, for the first time, we propose and experimentally demonstrated bilayer structure to achieve simultaneously independent manipulation of multi-physics field (dc electric fields and thermal) by directly solving the dc electric/ thermal field equations. This structure is composed of two layers: the outer layer is made of isotropic and homogeneous material, while the inner layer is fan-shape layer. Since it is not based on TO, it can be readily experimentally fabricated with naturally occurring materials. Experimentally, we has designed, fabricated and characterized two structures simultaneously behaving as dc electric cloak/ thermal concentrator and dc electric concentrator/ thermal cloak, respectively. The simulation results agree well with the experiment ones, thus confirming the feasibility of our methodology. This job provides a novel avenue to manipulation of multiphysics fields, thus might find potential applications in various areas.


title: Predicting the extinction of Ebola spreading in Liberia due to mitigation strategies

abstract: The Ebola virus is spreading throughout West Africa and is causing thousands of deaths. In order to quantify the effectiveness of different strategies for controlling the spread, we develop a mathematical model in which the propagation of the Ebola virus through Liberia is caused by travel between counties. For the initial months in which the Ebola virus spreads, we find that the arrival times of the disease into the counties predicted by our model are compatible with World Health Organization data, but we also find that reducing mobility is insufficient to contain the epidemic because it delays the arrival of Ebola virus in each county by only a few weeks. We study the effect of a strategy in which safe burials are increased and effective hospitalisation instituted under two scenarios: (i) one implemented in mid-July 2014 and (ii) one in mid-August---which was the actual time that strong interventions began in Liberia. We find that if scenario (i) had been pursued the lifetime of the epidemic would have been three months shorter and the total number of infected individuals 80\% less than in scenario (ii). Our projection under scenario (ii) is that the spreading will stop by mid-spring 2015.


title: Unknotting number for Lorenz knots

abstract: The unknotting number of a positive braid with n strands and k intersections is known to be equal to (k-n+1)/2. We consider Lorenz knots (which are positive braids) and, using a different method, find their unknotting numbers in terms of their positions on the Lorenz attractor.


title: A layman's note on a class of frequentist hypothesis testing problems

abstract: It is observed that for testing between simple hypotheses where the cost of Type I and Type II errors can be quantified, it is better to let the optimization choose the test size.


title: A proposal for regularly updated review/survey articles: "Perpetual Reviews"

abstract: We advocate the publication of review/survey articles that will be updated regularly, both in traditional journals and novel venues. We call these "perpetual reviews." This idea naturally builds on the dissemination and archival capabilities present in the modern internet, and indeed perpetual reviews exist already in some forms. Perpetual review articles allow authors to maintain over time the relevance of non-research scholarship that requires a significant investment of effort. Further, such reviews published in a purely electronic format without space constraints can also permit more pedagogical scholarship and clearer treatment of technical issues that remain obscure in a brief treatment.


title: Propagation and dispersion of sausage wave trains in magnetic flux tubes

abstract: A localized perturbation of a magnetic flux tube produces a pair of wave trains that propagate in opposite directions along the tube. These wave packets disperse as they propagate, where the extent of dispersion depends on the physical properties of the magnetic structure, on the length of the initial excitation, and on its nature (e.g., transverse or axisymmetric). In Oliver et al. (2014) we considered a transverse initial perturbation, whereas the temporal evolution of an axisymmetric one is examined here. In both papers we use a method based on Fourier integrals to solve the initial value problem. Previous studies on wave propagation in magnetic wave guides have emphasized that the wave train dispersion is influenced by the particular dependence of the group velocity on the longitudinal wavenumber. Here we also find that long initial perturbations result in low amplitude wave packets and that large values of the magnetic tube to environment density ratio yield longer wave trains. To test the detectability of propagating transverse or axisymmetric wave packets in magnetic tubes of the solar atmosphere (e.g., coronal loops, spicules, or prominence threads) a forward modelling of the perturbations must be carried out. This is left for a future work.


title: Azimuthal anisotropies in p+Pb collisions from classical Yang-Mills dynamics

abstract: We compute single and double inclusive gluon distributions in classical Yang-Mills simulations of proton-lead collisions and extract the associated transverse momentum dependent Fourier harmonics $v_2(p_T)$ and $v_3(p_T)$. Gluons have a large $v_2$ in the initial state, while odd harmonics such as $v_3$ vanish identically at the initial time $\tau=0^{+}$. By the time $\tau \lesssim 0.4\,{\rm fm/c}$ final state effects in the classical Yang-Mills evolution generate a non-zero $v_{3}$ and only mildly modify the gluon $v_{2}$. Unlike hydrodynamic flow, these momentum space anisotropies are uncorrelated with the global spatial anisotropy of the collision. A principal ingredient for the generation of $v_2$ and $v_3$ in this framework is the event-by-event breaking of rotational invariance in domains the size of the inverse of the saturation scale $Q_s$. In contrast to our findings in p+Pb collisions Yang-Mills simulations of lead-lead collisions generate much smaller values of $v_{2,3} (p_T)$ and additional collective flow effects are needed to explain experimental data. This is because the locally generated anisotropy due to the breaking of rotational invariance is depleted with the increase in the number of uncorrelated domains.


title: Effects of liquid pore water on acoustic wave propagation in snow as a Biot-type porous material

abstract: A method to estimate phase velocity and attenuation of acoustic waves in the presence of liquid water in a snowpack is presented. The method is based on Biot's theory of wave propagation in porous materials. Empirical relations and a priori information is used to characterize snow as a porous material as a function of porosity. Plane wave theory and an equivalent pore fluid are used to solve Biot's differential equations and to asses the impact of the air and water in the pore space. The liquid water in the pore space of a snow pack reduces the velocity of the first compressional wave by roughly 300 m/s for every 0.1 increase in liquid water saturation. Also the attenuation of the compressional waves is increased with increasing liquid water content. Two end member models for compaction are evaluated to asses the importance of an independent density measurement for an estimate of liquid pore water saturation in snow with acoustic waves. The two end members correspond to no compaction at all and to a melting sphere packing where the grains remain in contact. The change of velocity for the first compressional wave was found to strongly depend on the compaction model, while the velocity of the second compressional wave mainly depends on liquid water saturation. Also the attenuation for both compressional waves depends mainly on the liquid water saturation and only little on the compaction model. Finally, a simple field experiment illustrates the potential use of the method to estimate the liquid water content of snow with acoustic sensors.


title: On Piterbarg's max-discretisation theorem for homogeneous Gaussian random fields

abstract: Motivated by the papers of Piterbarg (2004) and H\"{u}sler (2004), in this paper the asymptotic relation between the maximum of a continuous dependent homogeneous Gaussian random field and the maximum of this field sampled at discrete time points is studied. It is shown that, for the weakly dependent case, these two maxima are asymptotically independent, dependent and coincide when the grid of the discrete time points is a sparse grid, Pickands grid and dense grid, respectively, while for the strongly dependent case, these two maxima are asymptotically totally dependent if the grid of the discrete time points is sufficiently dense, and asymptotically dependent if the the grid points are sparse or Pickands grids.


title: Dynamically Induced Planck Scale and Inflation

abstract: Theories where the Planck scale is dynamically generated from dimensionless interactions provide predictive inflationary potentials and super-Planckian field variations. We first study the minimal single-field realisation in the low-energy effective field theory limit, finding the predictions $n_s \approx 0.96$ for the spectral index and $r \approx 0.13$ for the tensor-to-scalar ratio, which can be reduced down to $\approx 0.04$ in presence of large couplings. Next we consider agravity as a dimensionless quantum gravity theory finding a multi-field inflation that converges towards an attractor trajectory that predicts $n_s\approx 0.96$ and $0.003<r<0.13$, interpolating between the quadratic and Starobinsky inflation. These theories relate the smallness of the weak scale to the smallness of inflationary perturbations: both arise naturally because of small couplings, implying a reheating temperature of $10^{7-9}$ GeV. A measurement of $r$ by Keck/Bicep3 would give us information on quantum gravity in the dimensionless scenario.


title: Approximately Counting H-Colourings is #BIS-Hard

abstract: We consider the problem of counting H-colourings from an input graph G to a target graph H. We show that if H is any fixed graph without trivial components, then the problem is as hard as the well-known problem #BIS, which is the problem of (approximately) counting independent sets in a bipartite graph. #BIS is a complete problem in an important complexity class for approximate counting, and is believed not to have an FPRAS. If this is so, then our result shows that for every graph H without trivial components, the H-colouring counting problem has no FPRAS. This problem was studied a decade ago by Goldberg, Kelk and Paterson. They were able to show that approximately sampling H-colourings is #BIS-hard, but it was not known how to get the result for approximate counting. Our solution builds on non-constructive ideas using the work of Lovasz.


title: Renormalized vacuum polarization of rotating black holes

abstract: Quantum field theory on rotating black hole spacetimes is plagued with technical difficulties. Here, we describe a general method to renormalize and compute the vacuum polarization of a quantum field in the Hartle-Hawking state on rotating black holes. We exemplify the technique with a massive scalar field on the warped AdS3 black hole solution to topologically massive gravity, a deformation of (2+1)-dimensional Einstein gravity. We use a "quasi-Euclidean" technique, which generalizes the Euclidean techniques used for static spacetimes, and we subtract the divergences by matching to a sum over mode solutions on Minkowski spacetime. This allows us, for the first time, to have a general method to compute the renormalized vacuum polarization (and, more importantly, the renormalized stress-energy tensor), for a given quantum state, on a rotating black hole, such as the physically relevant case of the Kerr black hole in four dimensions.


title: Counterintuitive effect of gravity on the heat capacity of a metal sphere: re-examination of a well-known problem

abstract: A well-known high-school problem asking the final temperature of two spheres that are given the same amount of heat, one lying on a table and the other hanging from a thread, is re-examined. The conventional solution states that the sphere on the table ends up colder, since thermal expansion raises its center of mass. It is found that this solution violates the second law of thermodynamics and is therefore incorrect. Two different new solutions are proposed. The first uses statistical mechanics, while the second is based on purely classical thermodynamical arguments. It is found that gravity produces a counterintuitive effect on the heat capacity, and the new answer to the problem goes in the opposite direction of what traditionally thought.


title: Revealing a hard X-ray spectral component that reverberates within one light hour of the central supermassive black hole in Ark 564

abstract: Ark 564 (z=0.0247) is an X-ray-bright NLS1. By using advanced X-ray timing techniques, an excess of "delayed" emission in the hard X-ray band (4-7.5 keV) following about 1000 seconds after "flaring" light in the soft X-ray band (0.4-1 keV) was recently detected. We report on the X-ray spectral analysis of eight XMM-Newton and one Suzaku observation of Ark 564. High-resolution spectroscopy was performed with the RGS in the soft X-ray band, while broad-band spectroscopy was performed with the EPIC-pn and XIS/PIN instruments. We analysed time-averaged, flux-selected, and time-resolved spectra. Despite the strong variability in flux during our observational campaign, the broad-band spectral shape of Ark 564 does not vary dramatically and can be reproduced either by a superposition of a power law and a blackbody emission or by a Comptonized power-law emission model. High-resolution spectroscopy revealed ionised gas along the line of sight at the systemic redshift of the source, with a low column density and a range of ionisation states. Broad-band spectroscopy revealed a very steep intrinsic continuum and a rather weak emission feature in the iron K band; modelling this feature with a reflection component requires highly ionised gas. A reflection-dominated or an absorption-dominated model are similarly able to well reproduce the time-averaged data from a statistical point of view, in both cases requiring contrived geometries and/or unlikely physical parameters. Finally, through time-resolved analysis we spectroscopically identified the "delayed" emission as a spectral hardening above ~4 keV; the most likely interpretation for this component is a reprocessing of the "flaring" light by gas located at 10-100 r_g from the central supermassive black hole that is so hot that it can Compton-upscatter the flaring intrinsic continuum emission.


title: The merger rate of galaxies in the Illustris Simulation: a comparison with observations and semi-empirical models

abstract: We have constructed merger trees for galaxies in the Illustris Simulation by directly tracking the baryonic content of subhalos. These merger trees are used to calculate the galaxy-galaxy merger rate as a function of descendant stellar mass, progenitor stellar mass ratio, and redshift. We demonstrate that the most appropriate definition for the mass ratio of a galaxy-galaxy merger consists in taking both progenitor masses at the time when the secondary progenitor reaches its maximum stellar mass. Additionally, we avoid effects from `orphaned' galaxies by allowing some objects to `skip' a snapshot when finding a descendant, and by only considering mergers which show a well-defined `infall' moment. Adopting these definitions, we obtain well-converged predictions for the galaxy-galaxy merger rate with the following main features, which are qualitatively similar to the halo-halo merger rate except for the last one: a strong correlation with redshift that evolves as $\sim (1+z)^{2.4-2.8}$, a power law with respect to mass ratio, and an increasing dependence on descendant stellar mass, which steepens significantly for descendant stellar masses greater than $\sim 2 \times 10^{11} \, {\rm M_{\odot}}$. These trends are consistent with observational constraints for medium-sized galaxies ($M_{\ast} \gtrsim 10^{10} \, {\rm M_{\odot}}$), but in tension with some recent observations of the close pair fraction for massive galaxies ($M_{\ast} \gtrsim 10^{11} \, {\rm M_{\odot}}$), which report a nearly constant or decreasing evolution with redshift. Finally, we provide a fitting function for the galaxy-galaxy merger rate which is accurate over a wide range of stellar masses, progenitor mass ratios, and redshifts.


title: Competing Abelian and non-Abelian topological orders in $\nu = 1/3+1/3$ quantum Hall bilayers

abstract: Bilayer quantum Hall systems, realized either in two separated wells or in the lowest two sub-bands of a wide quantum well, provide an experimentally realizable way to tune between competing quantum orders at the same filling fraction. Using newly developed density matrix renormalization group techniques combined with exact diagonalization, we return to the problem of quantum Hall bilayers at filling $\nu = 1/3 + 1/3$. We first consider the Coulomb interaction at bilayer separation $d$, bilayer tunneling energy $\Delta_\textrm{SAS}$, and individual layer width $w$, where we find a phase diagram which includes three competing Abelian phases: a bilayer-Laughlin phase (two nearly decoupled $\nu = 1/3$ layers); a bilayer-spin singlet phase; and a bilayer-symmetric phase. We also study the order of the transitions between these phases. A variety of non-Abelian phases have also been proposed for these systems. While absent in the simplest phase diagram, by slightly modifying the interlayer repulsion we find a robust non-Abelian phase which we identify as the "interlayer-Pfaffian" phase. In addition to non-Abelian statistics similar to the Moore-Read state, it exhibits a novel form of bilayer-spin charge separation. Our results suggest that $\nu = 1/3 + 1/3$ systems merit further experimental study.


title: On the serendipitous discovery of a Li-rich giant in the globular cluster NGC 362

abstract: We have serendipitously identified the first lithium-rich giant star located close to the red giant branch bump in a globular cluster. Through intermediate-resolution FLAMES spectra we derived a lithium abundance of A(Li)=2.55 (assuming local thermodynamical equilibrium), which is extremely high considering the star's evolutionary stage. Kinematic and photometric analysis confirm the object as a member of the globular cluster NGC 362. This is the fourth Li-rich giant discovered in a globular cluster but the only one known to exist at a luminosity close to the bump magnitude. The three previous detections are clearly more evolved, located close to, or beyond the tip of their red giant branch. Our observations are able to discard the accretion of planets/brown dwarfs, as well as an enhanced mass-loss mechanism as a formation channel for this rare object. Whilst the star sits just above the cluster bump luminosity, its temperature places it towards the blue side of the giant branch in the colour-magnitude diagram. We require further dedicated observations to unambiguously identify the star as a red giant: we are currently unable to confirm whether Li production has occurred at the bump of the luminosity function or if the star is on the pre zero-age horizontal branch. The latter scenario provides the opportunity for the star to have synthesised Li rapidly during the core helium flash or gradually during its red giant branch ascent via some extra mixing process.


title: Using Galaxy Pairs to Probe Star Formation During Major Halo Mergers

abstract: Currently-proposed galaxy quenching mechanisms predict very different behaviours during major halo mergers, ranging from significant quenching enhancement (e.g., clump-induced gravitational heating models) to significant star formation enhancement (e.g., gas starvation models). To test real galaxies' behaviour, we present an observational galaxy pair method for selecting galaxies whose host haloes are preferentially undergoing major mergers. Applying the method to central L* (10^10 Msun < M_* < 10^10.5 Msun) galaxies in the Sloan Digital Sky Survey (SDSS) at z<0.06, we find that major halo mergers can at most modestly reduce the star-forming fraction, from 59% to 47%. Consistent with past research, however, mergers accompany enhanced specific star formation rates for star-forming L* centrals: ~10% when a paired galaxy is within 200 kpc (approximately the host halo's virial radius), climbing to ~70% when a paired galaxy is within 30 kpc. No evidence is seen for even extremely close pairs (<30 kpc separation) rejuvenating star formation in quenched galaxies. For galaxy formation models, our results suggest: (1) quenching in L* galaxies likely begins due to decoupling of the galaxy from existing hot and cold gas reservoirs, rather than a lack of available gas or gravitational heating from infalling clumps, (2) state-of-the-art semi-analytic models currently over-predict the effect of major halo mergers on quenching, and (3) major halo mergers can trigger enhanced star formation in non-quenched central galaxies.


title: Effective-one-body Hamiltonian with next-to-leading order spin-spin coupling

abstract: We propose a way of including the next-to-leading (NLO) order spin-spin coupling into an effective-one-body (EOB) Hamiltonian. This work extends [S. Balmelli and P. Jetzer, Phys. Rev. D 87, 124036 (2013)], which is restricted to the case of equatorial orbits and aligned spins, to general orbits with arbitrary spin orientations. This is done applying appropriate canonical phase-space transformations to the NLO spin-spin Hamiltonian in Arnowitt-Deser-Misner (ADM) coordinates, and systematically adding "effectiv" quantities at NLO to all spin-squared terms appearing in the EOB Hamiltonian. As required by consistency, the introduced quantities reduce to zero in the test- mass limit. We expose the result both in a general gauge and in a gauge-fixed form. The last is chosen such as to minimize the number of new coefficients that have to be inserted into the effective spin squared. As a result, the 25 parameters that describe the ADM NLO spin-spin dynamics get condensed into only 12 EOB terms.


title: Periodograms for Multiband Astronomical Time Series

abstract: This paper introduces the multiband periodogram, a general extension of the well-known Lomb-Scargle approach for detecting periodic signals in time-domain data. In addition to advantages of the Lomb-Scargle method such as treatment of non-uniform sampling and heteroscedastic errors, the multiband periodogram significantly improves period finding for randomly sampled multiband light curves (e.g., Pan-STARRS, DES and LSST). The light curves in each band are modeled as arbitrary truncated Fourier series, with the period and phase shared across all bands. The key aspect is the use of Tikhonov regularization which drives most of the variability into the so-called base model common to all bands, while fits for individual bands describe residuals relative to the base model and typically require lower-order Fourier series. This decrease in the effective model complexity is the main reason for improved performance. We use simulated light curves and randomly subsampled SDSS Stripe 82 data to demonstrate the superiority of this method compared to other methods from the literature, and find that this method will be able to efficiently determine the correct period in the majority of LSST's bright RR Lyrae stars with as little as six months of LSST data. A Python implementation of this method, along with code to fully reproduce the results reported here, is available on GitHub.


title: A Multiwavelength Study of the Relativistic Tidal Disruption Candidate Sw J2058+05 at Late Times

abstract: We report a multiwavelength (X-ray, ultraviolet/optical/infrared, radio) analysis of the relativistic tidal disruption event candidate Sw J2058+05 from 3 months to 3 yr post-discovery in order to study its properties and compare its behavior with that of Sw J1644+57. Our main results are as follows. (1) The long-term X-ray light curve of Sw J2058+05 shows a remarkably similar trend to that of Sw J1644+57. After a prolonged power-law decay, the X-ray flux drops off rapidly by a factor of $\gtrsim 160$ within a span of $\Delta$$t$/$t$ $\le$ 0.95. Associating this sudden decline with the transition from super-Eddington to sub-Eddington accretion, we estimate the black hole mass to be in the range of $10^{4-6}$ M$_{\odot}$. (2) We detect rapid ($\lesssim 500$ s) X-ray variability before the dropoff, suggesting that, even at late times, the X-rays originate from close to the black hole (ruling out a forward-shock origin). (3) We confirm using HST and VLBA astrometry that the location of the source coincides with the galaxy's center to within $\lesssim 400$ pc (in projection). (4) We modeled Sw J2058+05's ultraviolet/optical/infrared spectral energy distribution with a single-temperature blackbody and find that while the radius remains more or less constant at a value of $63.4 \pm 4.5$ AU ($\sim 10^{15}$ cm) at all times during the outburst, the blackbody temperature drops significantly from $\sim$ 30,000 K at early times to a value of $\sim$ 15,000 K at late times (before the X-ray dropoff). Our results strengthen Sw J2058+05's interpretation as a tidal disruption event similar to Sw J1644+57.


title: The Brown-dwarf Atmosphere Monitoring (BAM) Project II: Multi-epoch monitoring of extremely cool brown dwarfs

abstract: With the discovery of Y dwarfs by the WISE mission, the population of field brown dwarfs now extends to objects with temperatures comparable to those of Solar System planets. To investigate the atmospheres of these newly identified brown dwarfs, we have conducted a pilot study monitoring an initial sample of three late T-dwarfs (T6.5, T8 and T8.5) and one Y-dwarf (Y0) for infrared photometric variability at multiple epochs. With J-band imaging, each target was observed for a period of 1.0h to 4.5h per epoch, which covers a significant fraction of the expected rotational period. These measurements represent the first photometric monitoring for these targets. For three of the four targets (2M1047, Ross 458C and WISE0458), multi-epoch monitoring was performed, with the time span between epochs ranging from a few hours to ~2 years. During the first epoch, the T8.5 target WISE0458 exhibited variations with a remarkable min-to-max amplitude of 13%, while the second epoch light curve taken ~2 years later did not note any variability to a 3% upper limit. With an effective temperature of ~600 K, WISE0458 is the coldest variable brown dwarf published to-date, and combined with its high and variable amplitude makes it a fascinating target for detailed follow-up. The three remaining targets showed no significant variations, with a photometric precision between 0.8% and 20.0%, depending on the target brightness. Combining the new results with previous multi-epoch observations of brown dwarfs with spectral types of T5 or later, the currently identified variables have locations on the colour-colour diagram better matched by theoretical models incorporating cloud opacities rather than cloud-free atmospheres. This preliminary result requires further study to determine if there is a definitive link between variability among late-T dwarfs and their location on the colour-colour diagram.


title: Measuring star formation with resolved observations: the test case of M33

abstract: Context. Measuring star formation at a local scale is important to constrain star formation laws. Yet, it is not clear whether and how the measure of star formation is affected by the spatial scale at which a galaxy is observed. Aims. We want to understand the impact of the resolution on the determination of the spatially resolved star formation rate (SFR) and other directly associated physical parameters such as the attenuation. Methods. We have carried out a multi-scale, pixel-by-pixel study of the nearby galaxy M33. Assembling FUV, Halpha, 8, 24, 70, and 100 micron maps, we have systematically compared the emission in individual bands with various SFR estimators from a resolution of 33 pc to 2084 pc. Results. We have found that there are strong, scale-dependent, discrepancies up to a factor 3 between monochromatic SFR estimators and Halpha+24 micron. The scaling factors between individual IR bands and the SFR show a strong dependence on the spatial scale and on the intensity of star formation. Finally, strong variations of the differential reddening between the nebular emission and the stellar continuum are seen, depending on the specific SFR (sSFR) and on the resolution. At the finest spatial scales, there is little differential reddening at high sSFR. The differential reddening increases with decreasing sSFR. At the coarsest spatial scales the differential reddening is compatible with the canonical value found for starburst galaxies. Conclusions. Our results confirm that monochromatic estimators of the SFR are unreliable at scales smaller than 1 kpc. Furthermore, the extension of local calibrations to high redshift galaxies presents non-trivial challenges as the properties of these systems may be poorly known.


title: PRIMUS: The Effect of Physical Scale on the Luminosity-Dependence of Galaxy Clustering via Cross-Correlations

abstract: We report small-scale clustering measurements from the PRIMUS spectroscopic redshift survey as a function of color and luminosity. We measure the real-space cross-correlations between 62,106 primary galaxies with PRIMUS redshifts and a tracer population of 545,000 photometric galaxies over redshifts from z=0.2 to z=1. We separately fit a power-law model in redshift and luminosity to each of three independent color-selected samples of galaxies. We report clustering amplitudes at fiducial values of z=0.5 and L=1.5 L*. The clustering of the red galaxies is ~3 times as strong as that of the blue galaxies and ~1.5 as strong as that of the green galaxies. We also find that the luminosity dependence of the clustering is strongly dependent on physical scale, with greater luminosity dependence being found between r=0.0625 Mpc/h and r=0.25 Mpc/h, compared to the r=0.5 Mpc/h to r=2 Mpc/h range. Moreover, over a range of two orders of magnitude in luminosity, a single power-law fit to the luminosity dependence is not sufficient to explain the increase in clustering at both the bright and faint ends at the smaller scales. We argue that luminosity-dependent clustering at small scales is a necessary component of galaxy-halo occupation models for blue, star-forming galaxies as well as for red, quenched galaxies.


title: On the Diffuse Lyman-alpha Halo Around Lyman-alpha Emitting Galaxies

abstract: Ly$\alpha$ photons scattered by neutral hydrogen atoms in the circumgalactic media or produced in the halos of star-forming galaxies are expected to lead to extended Ly$\alpha$ emission around galaxies. Such low surface brightness Ly$\alpha$ halos (LAHs) have been detected by stacking Ly$\alpha$ images of high-redshift star-forming galaxies. We study the origin of LAHs by performing radiative transfer modeling of nine $z=3.1$ Lyman-Alpha Emitters (LAEs) in a high resolution hydrodynamic cosmological galaxy formation simulation. We develop a method of computing the mean Ly$\alpha$ surface brightness profile of each LAE by effectively integrating over many different observing directions. Without adjusting any parameters, our model yields an average Ly$\alpha$ surface brightness profile in remarkable agreement with observations. We find that observed LAHs cannot be accounted for solely by photons originating from the central LAE and scattered to large radii by hydrogen atoms in the circumgalactic gas. Instead, Ly$\alpha$ emission from regions in the outer halo is primarily responsible for producing the extended LAHs seen in observations, which potentially includes both star-forming and cooling radiation. With the limit on the star formation contribution set by the ultra-violet (UV) halo measurement, we find that cooling radiation can play an important role in forming the extended LAHs. We discuss the implications and caveats of such a picture.


title: Extragalactic Transients in the Era of Wide-Field Radio Surveys. I. Detection Rates and Light Curve Characteristics

abstract: The impending era of wide-field radio surveys has the potential to revolutionize our understanding of astrophysical transients. Here we evaluate the prospects of a wide range of planned and hypothetical radio surveys using the properties and volumetric rates of known and hypothetical classes of extragalactic synchrotron radio transients (e.g., on- and off-axis gamma-ray bursts [GRB], supernovae, tidal disruption events [TDE], compact object mergers). Utilizing these sources and physically motivated considerations we assess the allowed phase-space of radio luminosity and peak timescale for extragalactic transients. We also include for the first time effects such as redshift evolution of the rates, K-corrections, and non-Euclidean luminosity distance, which affect the detection rates of the most sensitive surveys. The number of detected events is calculated by means of a Monte Carlo method, using the various survey properties (depth, cadence, area) and realistic detection criteria that include a cut on the minimum variability of the transients during the survey and an assessment of host galaxy contamination. Near-term GHz frequency surveys (ASKAP/VAST, Very Large Array Sky Survey) will detect few events: <~30-50 on- and off-axis long GRBs and off-axis tidal disruption events, and ~10-20 neutron star binary mergers if ~1% of the mergers result in a stable millisecond magnetar. Low-frequency surveys (e.g., LOFAR) are unlikely to detect any transients, while a hypothetical large-scale mm survey may detect ~40 on-axis long GRBs. On the other hand, SKA surveys at ~0.1-1 GHz have the potential to uncover thousands of transients, mainly on- and off-axis long GRBs, on-axis short GRBs, off-axis TDEs, and neutron star binary mergers with magnetar remnants.


title: Rational design of self-assembly pathways for complex multicomponent structures

abstract: The field of complex self-assembly is moving toward the design of multi-particle structures consisting of thousands of distinct building blocks. To exploit the potential benefits of structures with such `addressable complexity,' we need to understand the factors that optimize the yield and the kinetics of self-assembly. Here we use a simple theoretical method to explain the key features responsible for the unexpected success of DNA-brick experiments, which are currently the only demonstration of reliable self-assembly with such a large number of components. Simulations confirm that our theory accurately predicts the narrow temperature window in which error-free assembly can occur. Even more strikingly, our theory predicts that correct assembly of the complete structure may require a time-dependent experimental protocol. Furthermore, we predict that low coordination numbers result in non-classical nucleation behavior, which we find to be essential for achieving optimal nucleation kinetics under mild growth conditions. We also show that, rather surprisingly, the use of heterogeneous bond energies improves the nucleation kinetics and in fact appears to be necessary for assembling certain intricate three-dimensional structures. This observation makes it possible to sculpt nucleation pathways by tuning the distribution of interaction strengths. These insights not only suggest how to improve the design of structures based on DNA bricks, but also point the way toward the creation of a much wider class of chemical or colloidal structures with addressable complexity.


title: Electroweak naturalness in three-flavour Type I see-saw and implications for leptogenesis

abstract: In the Type I see-saw model, the naturalness requirement that corrections to the electroweak $\mu$ parameter not exceed 1 TeV results in a rough bound on the lightest right-handed neutrino mass, $M_{N_1}\lesssim 3\times 10^7$ GeV. In this letter we derive generic bounds applicable in any three-flavour Type I see-saw model. We find $M_{N_1}\lesssim 4\times 10^7$ GeV and $M_{N_2}\lesssim 7\times 10^7$ GeV. In the limit of one massless neutrino, there is no naturalness bound on $M_{N_3}$ in the Poincare protected decoupling limit. Our results confirm that no Type I see-saw model can explain the observed neutrino masses and baryogenesis via hierarchical ($N_1$-, $N_2$-, or $N_3$-dominated) thermal leptogenesis while remaining completely natural.


title: Tailoring Strong Lensing Cosmographic Observations

abstract: Strong lensing time delay cosmography has excellent complementarity with other dark energy probes, and will soon have abundant systems detected. We investigate two issues in the imaging and spectroscopic followup required to obtain the time delay distance. The first is optimization of spectroscopic resources. We develop a code to optimize the cosmological leverage under the constraint of constant spectroscopic time, and find that sculpting the lens system redshift distribution can deliver a 40% improvement in dark energy figure of merit. The second is the role of systematics, correlated between different quantities of a given system or model errors common to all systems. We show how the levels of different systematics affect the cosmological parameter estimation, and derive guidance for the fraction of double image vs quad image systems to follow as a function of differing systematics between them.


title: Cold parsec-scale gas in a zabs~0.1 sub-DLA with disparate H2 and 21-cm absorption

abstract: We present a detailed analysis of a H2-bearing metal-rich sub-damped Lyman-alpha system at zabs = 0.10115 towards the radio-loud quasar J0441-4313, at a projected separation of ~7.6 kpc from a star-forming galaxy. The H2, C I and Na I absorption are much stronger in the redder of the two components seen in the Hubble Space Telescope / Cosmic Origins Spectrograph spectrum. The best single component fit to the strong H2 component gives log N(H2) = 16.61 +/- 0.05. However, possible hidden saturation in the medium resolution spectrum can allow for log N(H2) to be as high as 18.9. The rotational excitation temperature of H2 in this component is 133 +33/-22 K. Photoionization models suggest 30-80% of the total N(H I) is associated with the strong H2 component, that has a density <= 100 cm^-3 and is subject to a radiation field that is <= 0.5 times the Galactic mean field. The Very Large Baseline Array 1.4 GHz continuum image of the radio source contains only 27% of the arcsecond scale emission. Using a previously published spectrum, no 21-cm absorption is found to be associated with the strong H2 component. This suggests that either the N(H I) associated with this component is <= 50% of the total N(H I) or the gas covering factor is <= 0.27. This is consistent with the results of the photoionization model that uses UV radiation due to stars in the associated galaxy. The 21-cm absorption previously reported from the weaker H2 component suggests a spin temperature of <= 90 K, at odds with the weakness of H2, C I and Na I absorption in this component. From the inferred physical and chemical conditions, we suggest that the gas may be tracing a recent metal-rich outflow from the host-galaxy.


title: Symmetry-protected Topological Phases at Finite Temperature

abstract: We have applied the recently developed theory of topological Uhlmann numbers to a representative model of a topological insulator in two dimensions, the Qi-Wu-Zhang model. We have found a stable symmetry-protected topological (SPT) phase under external thermal fluctuations in two-dimensions. A complete phase diagram for this model is computed as a function of temperature and coupling constants in the original Hamiltonian. It shows the appearance of large stable phases of matter with topological properties compatible with thermal fluctuations or external noise and the existence of critical lines separating abruptly trivial phases from topological phases. These novel critical temperatures represent thermal topological phase transitions. The initial part of the paper comprises a self-contained explanation of the Uhlmann geometric phase needed to understand the topological properties that it may acquire when applied to topological insulators and superconductors.


title: Center symmetry and Hagedorn spectrum

abstract: This paper explores the conjecture that large Nc gauge theories have a Hagedorn spectrum, if, and only if, they are confining and posses an explicit or emergent center symmetry. Evidence in support of this conjecture is presented. Many classes of large Nc gauge theories are considered. In all cases, we find that theories for which there exists a strong plausibility argument for a Hagedorn spectrum at large Nc are also believed to be confining and possess either an explicit center symmetric or have a strong plausibility argument for the existence of an emergent center symmetry at large Nc . Conversely, all theories we considered which are believed not to have a Hagedorn spectrum at large Nc, either were known not to be confining or else were believed to lack an emergent center symmetry. This is consistent with expectations based on the conjecture.


title: Asymmetric MHD Outflows/Jets from Accreting T Tauri Stars

abstract: Observations of jets from young stellar objects reveal the asymmetric outflows from some sources. A large set of $2.5$D MHD simulations has been carried out for axisymmetric viscous/diffusive disc accretion to rotating magnetized stars for the purpose of assessing the conditions where the outflows or jets are asymmetric relative to the equatorial plane. We consider initial magnetic fields that are symmetric about the equatorial plane and consist of a radially distributed field threading the disc (disc-field) and a stellar dipole field.({\bf 1}). For pure disc-fields the symmetry or asymmetry of the outflows is affected by the midplane plasma $\beta$ of the disc (where $\beta$ is the ratio of the plasma pressure to the magnetic pressure). For the low density discs with small plasma $\beta$ values, outflows are observed to be symmetric to within $10\%$ over timescales of hundreds of inner disc orbits. For the denser higher $\beta$ discs, the coupling of the upper and lower coronal plasmas is broken, and quasi-periodic field motion in the two hemispheres becomes different. This asymmetry leads to asymmetric episodic outflows. ({\bf 2.}) Accreting stars with a stellar dipole field and no disc-field exhibit episodic, two component outflows - a magnetospheric wind and an inner disc wind from somewhat larger radial distances. Both are characterized by similar velocity profiles but the magnetospheric wind has densities $\gtrsim 10$ times that of the disc wind. ({\bf 3}.)Adding a disc-field parallel to the stellar dipole field acts to enhance the magnetospheric winds but suppress the disc wind. ({\bf 4}.) In contrast, adding a disc-field which is anti-parallel to the stellar dipole field in the disc acts to suppress the magnetospheric and disc winds. Our simulations reproduce some key features of observations of asymmetric outflows of T Tauri stars.


title: Clustered Cepheid Variables 90 kiloparsec from the Galactic Center

abstract: Distant regions close to the plane of our Galaxy are largely unexplored by optical surveys as they are hidden by dust. We have used near-infrared data (that minimizes dust obscuration) from the ESO Public survey VISTA Variables of the Via Lactea (VVV) (Minniti et al. 2011; Saito et al. 2012; henceforth S12) to search for distant stars at low latitudes. We have discovered four Cepheid variables within an angular extent of one degree centered at Galactic longitude of $l = -27.4^\circ$ and Galactic latitude of $b = -1.08 ^\circ$. We use the tightly constrained period-luminosity relationship that these pulsating stars obey (Persson et al. 2004; Matsunaga et al. 2011) to derive distances. We infer an average distance to these Cepheid variables of 90 kpc. The Cepheid variables are highly clustered in angle (within one degree) and in distance (the standard deviation of the distances is 12 kpc). They are at an average distance of $\sim 2~\rm kpc$ from the plane and their maximum projected separation is $\sim 1~ \rm kpc$. These young ($\sim$ 100 Myr old), pulsating stars (Bono et al. 2005) are unexpected at such large distances from the Galactic disk, which terminates at $\sim$ 15 kpc (Minniti et al. 2011). The highly clustered nature in distance and angle of the Cepheid variables suggests that the stars may be associated with a dwarf galaxy, one that was earlier predicted by a dynamical analysis (Chakrabarti \& Blitz 2009).


title: The Three-Terminal Interactive Lossy Source Coding Problem

abstract: The three-node multiterminal lossy source coding problem is investigated. We derive an inner bound to the general rate-distortion region of this problem which is a natural extension of the seminal work by Kaspi'85 on the interactive two-terminal source coding problem. It is shown that this (rather involved) inner bound contains several rate-distortion regions of some relevant source coding settings. In this way, besides the non-trivial extension of the interactive two terminal problem, our results can be seen as a generalization and hence unification of several previous works in the field. Specializing to particular cases we obtain novel rate-distortion regions for several lossy source coding problems. We finish by describing some of the open problems and challenges. However, the general three-node multiterminal lossy source coding problem seems to offer a formidable mathematical complexity.


title: Properties of the Virial Expansion and Equation of State of Ideal Quantum Gases in Arbitrary Dimensions

abstract: The virial expansion of ideal quantum gases reveals some interesting and amusing properties when considered as a function of dimensionality $d$. In particular, the convergence radius $\rho_c(d)$ of the expansion is particulary large at {\em exactly\/} $d=3$ dimensions, $\rho_c(3) = 7.1068\ldots \times \lim_{d\to3} \rho_c(d)$. The same phenomenon occurs in a few other special (non-integer) dimensions. We explain the origin of these facts, and discuss more generally the structure of singularities governing the asymptotic behavior of the ideal gas virial expansion.


title: A second Higgs from the Higgs portal

abstract: In the Higgs portal framework, the Higgs field generally mixes with the Standard Model singlet leading to the existence of two states, one of which is identified with the 125 GeV scalar observed at the LHC. In this work, we analyse direct and indirect constraints on the second mass eigenstate and the corresponding mixing angle. The existence of the additional scalar can be beneficial as it can stabilise the otherwise-metastable electroweak vacuum. We find parameter regions where all of the bounds, including the stability constraints, are satisfied. We also study prospects for observing the decay of the heavier state into a pair of the 125 GeV Higgs-like scalars.


title: Visualizing gauge unification with high-scale thresholds

abstract: Tests of gauge coupling unification require knowledge of thresholds between the weak scale and the high scale of unification. If these scales are far separated, as is the case in most unification scenarios considered in the literature, the task can be factorized into IR and UV analyses. We advocate "$\Delta\lambda$ plots" as an efficient IR analysis projected to the high scale. The data from these plots gives an immediate qualitative guide to the size of threshold corrections needed at the high scale (e.g., the indices of high-scale representations) and provides precise quantitative data needed to test the viability of hypothesized high-scale unification theories. Such an approach shows more clearly the reasonable prospects of non-supersymmetric grand unification in large rank groups, and also shows the low summed values of high-scale threshold corrections required for supersymmetric unification. The latter may imply tuned cancellations of high-scale thresholds in theories based on weak-scale supersymmetry. For that reason we view non-supersymmetric unification to be just as viable as supersymmetric unification when confining ourselves only to the question of reasonable high-scale threshold corrections needed for exact unification. We illustrate these features for a non-supersymmetric $SO(10)$ graund unified theory and a supersymmetric $SU(5)$ theory.


title: Radiation Skyshine Calculation with MARS15 for the mu2e Experiment at Fermilab

abstract: The Fermilab Antiproton source is to be repurposed to provide an 8 kW proton beam to the Mu2e experiment by 1/3 integer, slow resonant extraction. Shielding provided by the existing facility must be supplemented with in-tunnel shielding to limit the radiation effective dose rate above the shield in the AP30 service building. In addition to the nominal radiation shield calculations, radiation skyshine calculations were required to ensure compliance with Fermilab Radiological Control Manual. A complete model of the slow resonant extraction system including magnets, electrostatic septa, magnetic fields, tunnel enclosure with shield, and a nearby exit stairway are included in the model. The skyshine model extends above the beam enclosure surface to 10 km vertically and 5 km radially.


title: On the Atiyah problem on hyperbolic configurations of four points

abstract: Given a configuration $\mathbf{x}$ of $n$ distinct points in hyperbolic $3$-space $H^3$, Michael Atiyah associated $n$ polynomials $p_1,\ldots,p_n$ of a variable $t \in \mathbb{C}P^1$, of degree $n-1$, and conjectured that they are linearly independent over $\mathbb{C}$, no matter which configuration $\mathbf{x}$ one starts with. We prove this conjecture for $n=4$ in two cases: in case the $4$ points are non-coplanar, and in case one of the points lies in the hyperbolic convex hull of the other three.


title: Enhancing non-melonic triangulations: A tensor model mixing melonic and planar maps

abstract: Ordinary tensor models of rank $D\geq 3$ are dominated at large $N$ by tree-like graphs, known as melonic triangulations. We here show that non-melonic contributions can be enhanced consistently, leading to different types of large $N$ limits. We first study the most generic quartic model at $D=4$, with maximally enhanced non-melonic interactions. The existence of the $1/N$ expansion is proved and we further characterize the dominant triangulations. This combinatorial analysis is then used to define a non-quartic, non-melonic class of models for which the large $N$ free energy and the relevant expectations can be calculated explicitly. They are matched with random matrix models which contain multi-trace invariants in their potentials: they possess a branched polymer phase and a 2D quantum gravity phase, and a transition between them whose entropy exponent is positive. Finally, a non-perturbative analysis of the generic quartic model is performed, which proves analyticity in the coupling constants in cardioid domains.


title: Gaussian Approximation Potentials: a brief tutorial introduction

abstract: We present a swift walk-through of our recent work that uses machine learning to fit interatomic potentials based on quantum mechanical data. We describe our Gaussian Approximation Potentials (GAP) framework, discussing a variety of descriptors, how to train the model on total energies and derivatives and the simultaneous use of multiple models. We also show a small example using QUIP, the software sandbox implementation of GAP that is available for non-commercial use.


title: Visualizing Marden's theorem with Scilab

abstract: A theorem which is named after the American Mathematician Moris Marden states a very surprising and interesting fact concerning the relationship between the points of a triangle in the complex plane and the zeros of two complex polynomials related to this triangle: "Suppose the zeroes z1, z2, and z3 of a third-degree polynomial p(z) are non-collinear. There is a unique ellipse inscribed in the triangle with vertices z1, z2, z3 and tangent to the sides at their midpoints: the Steiner in-ellipse. The foci of that ellipse are the zeroes of the derivative p'(z)." (Wikipedia contributors, "Marden's theorem", http://en.wikipedia.org/wiki/Marden%27s_theorem). This document describes how Scilab, a popular and powerful open source alternative to MATLAB, can be used to visualize the above stated theorem for arbitrary complex numbers z1, z2, and z3 which are not collinear. It is further demonstrated how the equations of the Steiner ellipses of a triangle in the complex plane can be calculated and plotted by applying this theorem.


title: Sparse Representation Classification Beyond L1 Minimization and the Subspace Assumption

abstract: The sparse representation classifier (SRC) has been utilized in various classification problems, which makes use of L1 minimization and works well for image recognition satisfying a subspace assumption. In this paper we propose a new implementation of SRC via screening, establish its equivalence to the original SRC under regularity conditions, and prove its classification consistency under a latent subspace model and contamination. The results are demonstrated via simulations and real data experiments, where the new algorithm achieves comparable numerical performance and significantly faster.


title: On the partial transpose of fermionic Gaussian states

abstract: We consider Gaussian states of fermionic systems and study the action of the partial transposition on the density matrix. It is shown that, with a suitable choice of basis, these states are transformed into a linear combination of two Gaussian operators that are uniquely defined in terms of the covariance matrix of the original state. In case of a reflection symmetric geometry, this result can be used to efficiently calculate a lower bound for a well-known entanglement measure, the logarithmic negativity. Furthermore, exact expressions can be derived for traces involving integer powers of the partial transpose. The method can also be applied to the quantum Ising chain and the results show perfect agreement with the predictions of conformal field theory.


title: Sufficient and Necessary Conditions for Limit Theorems for Quadratic Variations of Gaussian Sequences

abstract: Quadratic variations of Gaussian processes play important role in both stochastic analysis and in applications such as estimation of model parameters, and for this reason the topic has been extensively studied in the literature. In this article we study the problem for general Gaussian processes and we provide sufficient and necessary conditions for different types of convergence which include convergence in probability, almost sure convergence, $L^p$-convergence as well as convergence in law. Furthermore, we study general Gaussian vectors from which different interesting cases including first or second order quadratic variations can be studied by appropriate choice of the underlying vector. Finally, we provide a practical and simple approach to attack the problem which simplifies the existing methodology considerably.


title: Inflationary scenarios in Starobinsky model with higher order corrections

abstract: We consider the Starobinsky inflation with a set of higher order corrections parametrised by two real coefficients $\lambda_1, \lambda_2$. In the Einstein frame we have found a potential with the Starobinsky plateau, steep slope and possibly with an additional minimum, local maximum or a saddle point. We have identified three types of inflationary behaviour that may be generated in this model: i) inflation on the plateau, ii) at the local maximum (topological inflation), iii) at the saddle point. We have found limits on parameters $\lambda_i$ and initial conditions at the Planck scale which enable successful inflation and disable eternal inflation at the plateau. We have checked that the local minimum away from the GR vacuum is stable and that the field cannot leave it neither via quantum tunnelling nor via thermal corrections.


title: TASI Lectures on Flavor Physics

abstract: These notes overlap with lectures given at the TASI summer schools in 2014 and 2011, as well as at the European School of High Energy Physics in 2013. This is primarily an attempt at transcribing my hand-written notes, with emphasis on topics and ideas discussed in the lectures. It is not a comprehensive introduction or review of the field, nor does it include a complete list of references. I hope, however, that some may find it useful to better understand the reasons for excitement about recent progress and future opportunities in flavor physics.


title: Fractal Butterflies of Chiral Fermions in Bilayer Graphene: Phase Transitions and Emergent Properties

abstract: We report on our studies of fractal butterflies in biased bilayer graphene in the fractional quantum Hall effect (FQHE) regime. We have considered the case when the external periodic potential is present in one layer and have illustrated the effect of varying both the periodic potential strength and the bias voltage on the FQHE and the butterfly energy gaps. Interestingly, the butterfly spectra exhibits remarkable phase transitions between the FQHE gap and the butterfly gap for chiral electrons in bilayer graphene, by varying either the periodic potential strength or the bias voltage. We also find that, in addition to those phase transitions, by varying the bias voltage one can essentially control the periodic potential strength experienced by the electrons.


title: Assortative and disassortative mixing investigated using the spectra of graphs

abstract: We investigate the impact of degree-degree correlations on the spectra of networks. Even though density distributions exhibit drastic changes depending on the (dis)assortative mixing and the network architecture, the short range correlations in eigenvalues exhibit universal RMT predictions. The long range correlations turn out to be a measure of randomness in (dis)assortative networks. The analysis further provides insight in to the origin of high degeneracy at the zero eigenvalue displayed by majority of the biological networks.


title: Moving in a crowd: human perception as a multiscale process

abstract: The strategic behaviour of pedestrians is largely determined by how they perceive and react to neighbouring people. This issue is addressed in this paper by a model which combines, in a time and space-dependent way, discrete and continuous effects of pedestrian interactions. Numerical simulations and qualitative analysis suggest that human perception, and its impact on crowd dynamics, can be effectively modelled as a multiscale process based on a dual microscopic/macroscopic representation of groups of agents.


title: Mitigation of numerical Cerenkov radiation and instability using a hybrid finite difference-FFT Maxwell solver and a local charge conserving current deposit

abstract: A hybrid Maxwell solver for fully relativistic and electromagnetic (EM) particle-in-cell (PIC) codes is described. In this solver, the EM fields are solved in $k$ space by performing an FFT in one direction, while using finite difference operators in the other direction(s). This solver eliminates the numerical Cerenkov radiation for particles moving in the preferred direction. Moreover, the numerical Cerenkov instability (NCI) induced by the relativistically drifting plasma and beam can be eliminated using this hybrid solver by applying strategies that are similar to those recently developed for pure FFT solvers. A current correction is applied for the charge conserving current deposit to correctly account for the EM calculation in hybrid Yee-FFT solver. A theoretical analysis of the dispersion properties in vacuum and in a drifting plasma for the hybrid solver is presented, and compared with PIC simulations with good agreement obtained. This hybrid solver is applied to both 2D and 3D Cartesian and quasi-3D (in which the fields and current are decomposed into azimuthal harmonics) geometries. Illustrative results for laser wakefield accelerator simulation in a Lorentz boosted frame using the hybrid solver in the 2D Cartesian geometry are presented, and compared against results from 2D UPIC-EMMA simulation which uses a pure spectral Maxwell solver, and from OSIRIS 2D lab frame simulation using the standard Yee solver. Very good agreement is obtained which demonstrates the feasibility of using the hybrid solver for high fidelity simulation of relativistically drifting plasma with no evidence of the numerical Cerenkov instability.


title: The Arithmetic Cosine Transform: Exact and Approximate Algorithms

abstract: In this paper, we introduce a new class of transform method --- the arithmetic cosine transform (ACT). We provide the central mathematical properties of the ACT, necessary in designing efficient and accurate implementations of the new transform method. The key mathematical tools used in the paper come from analytic number theory, in particular the properties of the Riemann zeta function. Additionally, we demonstrate that an exact signal interpolation is achievable for any block-length. Approximate calculations were also considered. The numerical examples provided show the potential of the ACT for various digital signal processing applications.


title: Inflation in maximal gauged supergravities

abstract: We discuss the dynamics of multiple scalar fields and the possibility of realistic inflation in the maximal gauged supergravity. In this paper, we address this problem in the framework of recently discovered 1-parameter deformation of ${\rm SO}(4,4)$ and ${\rm SO}(5,3)$ dyonic gaugings, for which the base point of the scalar manifold corresponds to an unstable de Sitter critical point. In the gauge-field frame where the embedding tensor takes the value in the sum of the {\bf 36} and {\bf 36'} representations of ${\rm SL}(8)$, we present a scheme that allows us to derive an analytic expression for the scalar potential. With the help of this formalism, we derive the full potential and gauge coupling functions in analytic forms for the ${\rm SO}(3)\times {\rm SO}(3)$-invariant subsectors of ${\rm SO}(4,4)$ and ${\rm SO}(5,3)$ gaugings, and argue that there exist no new critical points in addition to those discovered so far. For the ${\rm SO}(4,4)$ gauging, we also study the behavior of 6-dimensional scalar fields in this sector near the Dall'Agata-Inverso de Sitter critical point at which the negative eigenvalue of the scalar mass square with the largest modulus goes to zero as the deformation parameter approaches a critical value. We find that when the deformation parameter is taken sufficiently close to the critical value, inflation lasts more than 60 e-folds even if the initial point of the inflaton allows an $O(0.1)$ deviation in Planck units from the Dall'Agata-Inverso critical point. It turns out that the spectral index $n_s$ of the curvature perturbation at the time of the 60 e-folding number is always about 0.96 and within the $1\sigma$ range $n_s=0.9639\pm0.0047$ obtained by Planck, irrespective of the value of the $\eta$ parameter at the critical saddle point. The tensor-scalar ratio predicted by this model is around $10^{-3}$ and is close to the value in the Starobinsky model.


title: Butterfly Factorization

abstract: The paper introduces the butterfly factorization as a data-sparse approximation for the matrices that satisfy a complementary low-rank property. The factorization can be constructed efficiently if either fast algorithms for applying the matrix and its adjoint are available or the entries of the matrix can be sampled individually. For an $N \times N$ matrix, the resulting factorization is a product of $O(\log N)$ sparse matrices, each with $O(N)$ non-zero entries. Hence, it can be applied rapidly in $O(N\log N)$ operations. Numerical results are provided to demonstrate the effectiveness of the butterfly factorization and its construction algorithms.


title: Artificial neural networks in calibration of nonlinear mechanical models

abstract: Rapid development in numerical modelling of materials and the complexity of new models increases quickly together with their computational demands. Despite the growing performance of modern computers and clusters, calibration of such models from noisy experimental data remains a nontrivial and often computationally exhaustive task. The layered neural networks thus represent a robust and efficient technique to overcome the time-consuming simulations of a calibrated model. The potential of neural networks consists in simple implementation and high versatility in approximating nonlinear relationships. Therefore, there were several approaches proposed to accelerate the calibration of nonlinear models by neural networks. This contribution reviews and compares three possible strategies based on approximating (i) model response, (ii) inverse relationship between the model response and its parameters and (iii) error function quantifying how well the model fits the data. The advantages and drawbacks of particular strategies are demonstrated on the calibration of four parameters of the affinity hydration model from simulated data as well as from experimental measurements. This model is highly nonlinear, but computationally cheap thus allowing its calibration without any approximation and better quantification of results obtained by the examined calibration strategies. The paper can be thus viewed as a guide intended for the engineers to help them select an appropriate strategy in their particular calibration problems.


title: Behavior of canonical divisors under purely inseparable base changes

abstract: Let $k$ be an imperfect field. Let $X$ be a regular variety over $k$ and set $Y$ to be the normalization of $(X \times_k k^{1/p^{\infty}})_{{\rm red}}$. In this paper, we show that $K_Y+C=f^*K_X$ for some effective divisor $C$ on $Y$. We obtain the following three applications. First, we show that a $K_X$-trivial fiber space with non-normal fibers is uniruled. Second, we prove that general fibers of Mori fiber spaces are rationally chain connected. Third, we obtain a weakening of the cone theorem for surfaces and threefolds defined over an imperfect field.


title: New knotted solutions of Maxwell's equations

abstract: In this note we have further developed the study of topologically non-trivial solutions of vacuum electrodynamics. We have discovered a novel method of generating such solutions by applying conformal transformations with complex parameters on known solutions expressed in terms of Bateman's variables. This has enabled us to get a wide class of solutions from the basic configuration like constant electromagnetic fields and plane-waves. We have introduced a covariant formulation of the Bateman's construction and discussed the conserved charges associated with the conformal group as well as a set of four types of conserved helicities. We have also given a formulation in terms of quaternions. This led to a simple map between the electromagnetic knotted and linked solutions into flat connections of $SU(2)$ gauge theory. We have computed the corresponding CS charge in a class of solutions and it takes integer values.

categories: hep-th math-ph math.MP physics.class-ph physics.optics


title: Abundance theorem for surfaces over imperfect fields

abstract: In this paper, we show the abundance theorem for log canonical surfaces over fields of positive characteristic.