title: How gas-dynamic flare models powered by Petschek reconnection differ from those with ad hoc energy sources

abstract: Aspects of solar flare dynamics, such as chromospheric evaporation and flare light-curves, have long been studied using one-dimensional models of plasma dynamics inside a static flare loop, subjected to some energy input. While extremely successful at explaining the observed characteristics of flares, all such models so far have specified energy input ad hoc, rather than deriving it self-consistently. There is broad consensus that flares are powered by magnetic energy released through reconnection. Recent work has generalized Petschek's basic reconnection scenario, topological change followed by field line retraction and shock heating, to permit its inclusion into a one-dimensional flare loop model. Here we compare the gas dynamics driven by retraction and shocking to those from more conventional static loop models energized by ad hoc source terms. We find significant differences during the first minute, when retraction leads to larger kinetic energies and produces higher densities at the loop top, while ad hoc heating tends to rarify the loop top. The loop-top density concentration is related to the slow magnetosonic shock, characteristic of Petschek's model, but persists beyond the retraction phase occurring in the outflow jet. This offers an explanation for observed loop-top sources of X-ray and EUV emission, with advantages over that provided by ad hoc heating scenarios. The cooling phases of the two models are, however, notably similar to one another, suggesting observations at that stage will yield little information on the nature of energy input.


title: Springer correspondence, hyperelliptic curves, and cohomology of Fano varieties

abstract: In \cite{CVX3}, we have established a Springer theory for the symmetric pair $(\operatorname{SL}(N),\operatorname{SO}(N))$. In this setting we obtain representations of (the Tits extension) of the braid group rather than just Weyl group representations. These representations arise from cohomology of families of certain (Hessenberg) varieties. In this paper we determine the Springer correspondence explicitly for IC sheaves supported on order 2 nilpotent orbits. In this process we encounter universal families of hyperelliptic curves. As an application we calculate the cohomolgy of Fano varieties of $k$-planes in the smooth intersection of two quadrics in an even dimensional projective space.


title: Additive triples of bijections, or the toroidal semiqueens problem

abstract: We prove an asymptotic for the number of additive triples of bijections $\{1,\dots,n\}\to\mathbb{Z}/n\mathbb{Z}$, that is, the number of pairs of bijections $\pi_1,\pi_2\colon \{1,\dots,n\}\to\mathbb{Z}/n\mathbb{Z}$ such that the pointwise sum $\pi_1+\pi_2$ is also a bijection. This problem is equivalent to counting the number of orthomorphisms or complete mappings of $\mathbb{Z}/n\mathbb{Z}$, to counting the number of arrangements of $n$ mutually nonattacking semiqueens on an $n\times n$ toroidal chessboard, and to counting the number of transversals in a cyclic Latin square. The method of proof is a version of the Hardy--Littlewood circle method from analytic number theory, adapted to the group $(\mathbb{Z}/n\mathbb{Z})^n$.


title: Actinide chemistry using singlet-paired coupled cluster and its combinations with density functionals

abstract: Singlet-paired coupled cluster doubles (CCD0) is a simplification of CCD that relinquishes a fraction of dynamic correlation in order to be able to describe static correlation. Combinations of CCD0 with density functionals that recover specifically the dynamic correlation missing in the former have also been developed recently. Here, we assess the accuracy of CCD0 and CCD0+DFT (and variants of these using Brueckner orbitals) as compared to well-established quantum chemical methods for describing ground-state properties of singlet actinide molecules. The $f^0$ actinyl series (UO$_2^{2+}$, NpO$_2^{2+}$, PuO$_2^{2+}$), the isoelectronic NUN, and Thorium (ThO, ThO$^{2+}$) and Nobelium (NoO, NoO$_2$) oxides are studied.


title: Hot subdwarfs in the Galactic halo - Tracers of prominent events in late stellar evolution

abstract: Hot subdwarf stars (sdO/Bs) are the stripped cores of red giants located at the bluest extension of the horizontal branch. They constitute the dominant population of UV-bright stars in old stellar environments and are most likely formed by binary interactions. We perform the first systematic, spectroscopic analysis of a sample of those stars in the Galactic halo based on data from SDSS. In the course of this project we discovered 177 close binary candidates. A significant fraction of the sdB binaries turned out to have close substellar companions, which shows that brown dwarfs and planets can significantly influence late stellar evolution. Close hot subdwarf binaries with massive white dwarf companions on the other hand are good candidates for the progenitors of type Ia supernovae. We discovered a hypervelocity star, which not only turned out to be the fastest unbound star known in our Galaxy, but also the surviving companion of such a supernova explosion.


title: Holomorph of generalized Bol loops II

abstract: The notion of the holomorph of a generalized Bol loop (GBL) is characterized afresh. The holomorph of a right inverse property loop (RIPL) is shown to be a GBL if and only if the loop is a GBL and some bijections of the loop are right (middle) regular. The holomorph of a RIPL is shown to be a GBL if and only if the loop is a GBL and some elements of the loop are right (middle) nuclear. Necessary and sufficient condition for the holomorph of a RIPL to be a Bol loop are deduced. Some algebraic properties and commutative diagrams are established for a RIPL whose holomorph is a GBL.


title: One-point concentration of the clique and chromatic numbers of the random Cayley graph on F_2^n

abstract: Green showed that there exist constants $C_1,C_2>0$ such that the clique number $\omega$ of the random Cayley graph on $\mathbb{F}_2^n$ satisfies $\lim_{n\to\infty}\mathbb{P}(C_1n\log n < \omega < C_2n\log n)=1$. In this paper we find the best possible $C_1$ and $C_2$. Moreover, we prove that for $n$ in a set of density $1$, clique number is actually concentrated on a single value. As a simple consequence of these results, we also prove the one-point concentration result for the chromatic number, thus proving the $\mathbb{F}_2^n$ analogue of the famous conjecture by Bollob\'{a}s and giving almost the complete answer to the question by Green.


title: Seeing lens imaging as a superposition of multiple views

abstract: In the conventional approach to lens imaging, rays are used to map object points to image points. However, many students have a need to think of the image as a whole. To answer this need, lens imaging is reinterpreted as a superposition of sharp images from different viewpoints. These so-called elemental images are uncovered by covering the lens with a pinhole array. Rays are introduced to connect elemental images. Lens ray diagrams are constructed based on bundles of elemental images. The conventional construction method is included as a special case. The proposed approach proceeds from concrete images to abstract rays.


title: R\'enyi Entropy for the $\sun1$ WZW model on the torus

abstract: The $\sun1$ WZW model is constructed on a n-sheeted branched torus, which allows the investigation of the R\'enyi entropy for a single interval at finite temperature. The small and large interval limits, as well as the low temperature expansion are presented for this theory.


title: Adapting a Fourier pseudospectral method to Dirichlet boundary conditions for Rayleigh--B\'enard convection

abstract: We present the adaptation to non--free boundary conditions of a pseudospectral method based on the (complex) Fourier transform. The method is applied to the numerical integration of the Oberbeck--Boussinesq equations in a Rayleigh--B\'enard cell with no-slip boundary conditions for velocity and Dirichlet boundary conditions for temperature. We show the first results of a 2D numerical simulation of dry air convection at high Rayleigh number ($R\sim10^9$). These results are the basis for the later study, by the same method, of wet convection in a solar still.


title: Frog model wakeup time on the complete graph

abstract: The frog model is a system of random walks where active particles set sleeping particles in motion. On the complete graph with n vertices it is equivalent to a well-understood rumor spreading model. We given an alternate and elementary proof that the wake-up time, i.e. the expected time for every particle to be activated, is Theta(log n). Additionally, we give an explicit distributional equation for the wakeup time as a weighted sum of geometric random variables. This project was part of the University of Washington Research Experience for Undergraduates program.


title: Concerning the variability of beta-decay measurements

abstract: Many experiments have been carried out to study the beta-decay rates of a variety of nuclides, and many - but not all - of these experiments yield evidence of variability of these rates. While there is as yet no accepted theory to explain patterns in the results, a number of conjectures have been proposed. We discuss three prominent conjectures (which are not mutually exclusive) - that variability of beta-decay rates may be due to (a) environmental influences, (b) solar neutrinos, and (c) cosmic neutrinos. We find evidence in support of each of these conjectures.


title: Synthesis of Electromagnetic Metasurfaces: Principles and Illustrations

abstract: The paper presents partial overview of the mathematical synthesis and the physical realization of metasurfaces, and related illustrative examples. The synthesis consists in determining the exact tensorial surface susceptibility functions of the metasurface, based on generalized sheet transition conditions, while the realization deals with both metallic and dielectric scattering particle structures. The examples demonstrate the capabilities of the synthesis and realization techniques, thereby showing the plethora of possible metasurface field transmission and subsequent applications. The first example is the design of two diffraction engineering birefringent metasurfaces performing polarization beam splitting and orbital angular momentum multiplexing, respectively. Next, we discuss the concept of the "transistor" metasurface, which is an electromagnetic linear switch based on destructive interferences. Then, we introduce a non-reciprocal non-gyrotropic metasurface using a pick-up circuit radiator (PCR) architecture. Finally, the implementation of all-dielectric metasurfaces for spatial dispersion engineering is discussed.


title: Extractors in Paley graphs: a random model

abstract: A well-known conjecture in analytic number theory states that for every pair of sets $X,Y\subset\mathbb{Z}/p\mathbb{Z}$, each of size at least $\log ^C p$ (for some constant $C$) we have that the number of pairs $(x,y)\in X\times Y$ such that $x+y$ is a quadratic residue modulo $p$ differs from $\frac12|X||Y|$ by $o\left(|X||Y|\right)$. We address the probabilistic analogue of this question, that is for every fixed $\delta>0$, given a finite group $G$ and $A\subset G$ a random subset of density $\frac12$, we prove that with high probability for all subsets $|X|,|Y|\geq \log ^{2+\delta} |G|$, the number of pairs $(x,y)\in X\times Y$ such that $xy\in A$ differs from $\frac12|X||Y|$ by $o\left(|X||Y|\right)$.


title: Dynamical ambiguities in models with spontaneous Lorentz violation

abstract: Spontaneous Lorentz violation is a viable mechanism to look for Planck scale physics. In this work, we study spontaneous Lorentz violation models, in flat spacetime, where a vector field produces such a violation and matter is modeled by a complex scalar field. We show that it is possible to construct a Hamilton density for which the evolution respects the dynamical constraints. However, we also find that the initial data, as required by standard field theory, does not determine the fields evolution in a unique way. In addition, we present some examples where the physical effects of such ambiguities can be recognized. As a consequence, the proposals in which the electromagnetic and gravitational interactions emerge from spontaneous Lorentz violation are challenged.


title: Spin-polarized neutron matter: Critical unpairing and BCS-BEC precursor

abstract: We obtain the critical magnetic field required for complete destruction of $S$-wave pairing in neutron matter, thereby setting limits on the pairing and superfluidity of neutrons in the crust and outer core of magnetars. We find that for fields $B \ge 10^{17}$ G the neutron fluid is non-superfluid -- if weaker spin-1 superfluidity does not intervene -- a result with profound consequences for the thermal, rotational, and oscillatory behavior of magnetars. Because the dineutron is not bound in vacuum, cold dilute neutron matter cannot exhibit a proper BCS-BEC crossover. Nevertheless, owing to the strongly resonant behavior of the $nn$ interaction at low densities, neutron matter shows a precursor of the BEC state, as manifested in Cooper-pair correlation lengths {being} comparable to the interparticle distance. We make a systematic quantitative study of this type of BCS-BEC crossover in the presence of neutron fluid spin-polarization induced by an ultra-strong magnetic field. We evaluate the Cooper pair wave-function, quasiparticle occupation numbers, and quasiparticle spectra for densities and temperatures spanning the BCS-BEC crossover region. The phase diagram of spin-polarized neutron matter is constructed and explored at different polarizations.


title: A Weak Galerkin Finite Element Method for A Type of Fourth Order Problem Arising From Fluorescence Tomography

abstract: In this paper, a new and efficient numerical algorithm by using weak Galerkin (WG) finite element methods is proposed for a type of fourth order problem arising from fluorescence tomography(FT). Fluorescence tomography is an emerging, in vivo non-invasive 3-D imaging technique which reconstructs images that characterize the distribution of molecules that are tagged by fluorophores. Weak second order elliptic operator and its discrete version are introduced for a class of discontinuous functions defined on a finite element partition of the domain consisting of general polygons or polyhedra. An error estimate of optimal order is derived in an $H^2$-equivalent norm for the WG finite element solutions. Error estimates in the usual $L^2$ norm are established, yielding optimal order of convergence for all the WG finite element algorithms except the one corresponding to the lowest order (i.e., piecewise quadratic elements). Some numerical experiments are presented to illustrate the efficiency and accuracy of the numerical scheme.


title: Fast and Scalable Structural SVM with Slack Rescaling

abstract: We present an efficient method for training slack-rescaled structural SVM. Although finding the most violating label in a margin-rescaled formulation is often easy since the target function decomposes with respect to the structure, this is not the case for a slack-rescaled formulation, and finding the most violated label might be very difficult. Our core contribution is an efficient method for finding the most-violating-label in a slack-rescaled formulation, given an oracle that returns the most-violating-label in a (slightly modified) margin-rescaled formulation. We show that our method enables accurate and scalable training for slack-rescaled SVMs, reducing runtime by an order of magnitude compared to previous approaches to slack-rescaled SVMs.


title: Root-counting measures of Jacobi polynomials and topological types and critical geodesics of related quadratic differentials

abstract: Two main topics of this paper are asymptotic distributions of zeros of Jacobi polynomials and topology of critical trajectories of related quadratic differentials. First, we will discuss recent developments and some new results concerning the limit of the root-counting measures of these polynomials. In particular, we will show that the support of the limit measure sits on the critical trajectories of a quadratic differential of the form Q(z)dz^2=(az^2+bz+c)dz^2/(z^2-1)^2. Then we will give a complete classification, in terms of complex parameters a, b, and c, of possible topological types of critical geodesics for the quadratic differential of this type.


title: Anticommutativity of Symmetric Elements under Generalized Oriented Involutions

abstract: Let $R$ be a ring with $char(R)\neq2$ whose unit group are denoted by $\mathcal{U}(R)$, $G$ a group with involution $*$, and $\sigma:G\rightarrow\mathcal{U}(R)$ a nontrivial group homomorphism, with $ker\ \sigma=N$, satisfying $xx^*\in N$ for all $x\in G$. Let $RG$ be the group ring of $G$ over $R$ and define the involution $\sigma*$ in $RG$ by $\left( \sum_{x\in G}\alpha_xx\right)^{\sigma*}=\sum_{x\in G}\sigma(x)\alpha_xx^*$. In this paper, we will classify the group rings $RG$ such that $\mathcal{S}$ is anticommutative, where $\mathcal{S}$ is the largest subset of $(RG)^+=\left\{ \alpha\in RG: \alpha^{\sigma *}=\alpha\right\}$ that can satisfy anticommutativity under $char(R)\neq2$.


title: Almost Primes in Almost All Short Intervals

abstract: Let $E_k$ be the set of positive integers having exactly $k$ prime factors. We show that almost all intervals $[x,x+\log^{1+\varepsilon} x]$ contain $E_3$ numbers, and almost all intervals $[x,x+\log^{3.51} x]$ contain $E_2$ numbers. By this we mean that there are only $o(X)$ integers $1\leq x\leq X$ for which the mentioned intervals do not contain such numbers. The result for $E_3$ numbers is optimal up to the $\varepsilon$ in the exponent. The theorem on $E_2$ numbers improves a result of Harman, which had the exponent $7+\varepsilon$ in place of $3.51$. We will also consider general $E_k$ numbers, and find them on intervals whose lengths approach $\log x$ as $k\to \infty$.


title: Stellar orbital properties as diagnostics of the origin of the stellar halo

abstract: We examine metallicities, ages and orbital properties of halo stars in a Milky-Way like disk galaxy formed in the cosmological hydrodynamical MaGICC simulations. Halo stars were either accreted from satellites or they formed in situ in the disk or bulge of the galaxy and were then kicked up into the halo ("in situ/ kicked-up" stars). Regardless of where they formed both types show surprisingly similar orbital properties: the majority of stars of both types are on short-axis tubes with the same sense of rotation as the disk -- implying that a large fraction of satellites are accreted onto the halo with the same sense of angular momentum as the disk.


title: Photoconversion in the HIT solar cells: Theory vs experiment

abstract: We obtain theoretical expressions for the photocurrent in the Heterojunction solar cells with Intrinsic Thin layer (HIT cells). Our calculations take into account tunneling of electrons and holes through wide-bandgap layers of $\alpha$-Si:H or $\alpha$-SiC:H. We introduce the criteria, under which tunneling does not lead to the deterioration of solar cell characteristics, in particular, to the reduction of the short-circuit current and open-circuit voltage. We propose an algorithm to compute the photoconversion efficiency of HIT elements, taking into account the peculiarities of the open-circuit voltage generation, in particular, its rather high values. We test our theoretical predictions against the experimental results. For this, we fabricate HIT elements with the efficiency of about $20\,\%$. We measured the temperature dependence of the short-circuit current, open-circuit voltage, photoconversion power, and fill factor of the current-voltage curve of these elements in a wide temperature range from 80 to 420\,K. In the low-temperature range, the open-circuit voltage and the photoconversion power decrease on cooling. At $T \ge 200$\,K, the theoretical expressions and the experimental curves agree rather well. The behavior of the fill factor and output power at low temperatures is explained by the increase of the series resistance on cooling. We discuss the reasons behind the reduction of the power temperature coefficient in HIT elements. We show that they are related to the low value of the combined surface and volume recombination rate. Finally, we derive a theoretical expression for the HIT element's operation temperature under natural working conditions.


title: A large scale prediction of bacteriocin gene blocks suggests a wide functional spectrum for bacteriocins

abstract: Bacteriocins are peptide-derived molecules produced by bacteria, whose recently-discovered functions include virulence factors and signalling molecules as well as their better known roles as antibiotics. To date, close to five hundred bacteriocins have been identified and classified. Recent discoveries have shown that bacteriocins are highly diverse and widely distributed among bacterial species. Given the heterogeneity of bacteriocin compounds, many tools struggle with identifying novel bacteriocins due to their vast sequence and structural diversity. Many bacteriocins undergo post-translational processing or modifications necessary for the biosynthesis of the final mature form. Enzymatic modification of bacteriocins as well as their export is achieved by proteins whose genes are often located in a discrete gene cluster proximal to the bacteriocin precursor gene, referred to as \textit{context genes} in this study. Although bacteriocins themselves are structurally diverse, context genes have been shown to be largely conserved across unrelated species. Using this knowledge, we set out to identify new candidates for context genes which may clarify how bacteriocins are synthesized, and identify new candidates for bacteriocins that bear no sequence similarity to known toxins. To achieve these goals, we have developed a software tool, Bacteriocin Operon and gene block Associator (BOA) that can identify homologous bacteriocin associated gene clusters and predict novel ones. We discover that several phyla have a strong preference for bactericon genes, suggesting distinct functions for this group of molecules. Availability: https://github.com/idoerg/BOA


title: ManeParse: Mathematica Toolbox for PDF Uncertainties and Application to New Physics Searches

abstract: As the LHC begins Run 2 at an even higher energy, one of the top priorities will be to search for new particles (possibly from SUSY) at the highest energy scales. In addition to direct production of new particles, they can mix with Standard Model (SM) particles to yield discrepancies from the usual predictions. To distinguish this new physics from old uncertainties, we need tools to precisely quantify the uncertainties of the SM predictions. Here we present a versatile set of utility functions built with Mathematica that is capable of working with a variety of PDF formats including the recent LHAPDF6 format. This software can perform PDF calculations within the Mathematica framework and compare results from different PDF collaborations. The package includes both the central PDF value as well as the full error sets needed for PDF uncertainty analysis; a variety of sample error definitions are implemented. We demonstrate this package for the case of a new heavy scalar particle production at the LHC.


title: Suppression of type I migration by disk winds

abstract: Planets less massive than Saturn tend to rapidly migrate inward in protoplanetary disks. This is the so-called type I migration. Simulations attempting to reproduce the observed properties of exoplanets show that type I migration needs to be significantly reduced over a wide region of the disk for a long time. However, the mechanism capable of suppressing type I migration over a wide region has remained elusive. The recently found turbulence-driven disk winds offer new possibilities. We investigate the effects of disk winds on the disk profile and type I migration for a range of parameters that describe the strength of disk winds. We also examine the in situ formation of close-in super-Earths in disks that evolve through disk winds. The disk profile, which is regulated by viscous diffusion and disk winds, was derived by solving the diffusion equation. We carried out a number of simulations and plot here migration maps that indicate the type I migration rate. We also performed N-body simulations of the formation of close-in super-Earths from a population of planetesimals and planetary embryos. We define a key parameter, Kw, which determines the ratio of strengths between the viscous diffusion and disk winds. For a wide range of Kw, the type I migration rate is presented in migration maps. These maps show that type I migration is suppressed over the whole close-in region when the effects of disk winds are relatively strong (Kw < 100). From the results of N-body simulations, we see that type I migration is significantly slowed down assuming Kw = 40. We also show that the results of N-body simulations match statistical orbital distributions of close-in super-Earths.


title: The black hole spin influence on accretion disk neutrino detection

abstract: Neutrinos are copiously emitted from black hole accretion disks playing a fundamental role in their evolution, as well as in the production of gamma ray bursts and r-process nucleosynthesis. The black hole generates a strong gravitational field able to change the properties of the emerging neutrinos. We study the influence of the black hole spin on the structure of the neutrino surfaces, neutrino luminosities, average neutrino energies, and event counts at SuperK. We consider several disk models and provide estimates that cover different black hole efficiency scenarios. We discuss the influence of the detector's inclination with respect to the axis of the torus on neutrino properties. We find that tori around spinning black holes have larger luminosities, energies and rates compared to tori around static black holes, and that the inclination of the observer causes a reduction in the luminosities and detection rates but an increase in the average energies.


title: Quasinormal modes of Four Dimensional Topological Nonlinear Charged Lifshitz Black Holes

abstract: We study scalar perturbations of four dimensional topological nonlinear charged Lifshitz black holes with spherical and plane transverse sections, and we find numerically the quasinormal modes for scalar fields. Then, we study the stability of these black holes under massive and massless scalar field perturbations. We focus our study on the dependence of the dynamical exponent, the nonlinear exponent, the angular momentum and the mass of the scalar field in the modes. It is found that the modes are overdamped depending strongly on the dynamical exponent and the angular momentum of the scalar field for a spherical transverse section. In constrast, for plane transverse sections the modes are always overdamped.


title: Size biased couplings and the spectral gap for random regular graphs

abstract: Let $\lambda$ be the second largest eigenvalue in absolute value of a uniform random $d$-regular graph on $n$ vertices. It was famously conjectured by Alon and proved by Friedman that if $d$ is fixed independent of $n$, then $\lambda=2\sqrt{d-1} +o(1)$ with high probability. In the present work we show that $\lambda=O(\sqrt{d})$ continues to hold with high probability as long as $d=O(n^{2/3})$, making progress towards a conjecture of Vu that the bound holds for all $1\le d\le n/2$. Prior to this work the best result was obtained by Broder, Frieze, Suen and Upfal (1999) using the configuration model, which hits a barrier at $d=o(n^{1/2})$. We are able to go beyond this barrier by proving concentration of measure results directly for the uniform distribution on $d$-regular graphs. These come as consequences of advances we make in the theory of concentration by size biased couplings. Specifically, we obtain Bennett-type tail estimates for random variables admitting certain unbounded size biased couplings.


title: On reduced models for gravity waves generated by moving bodies

abstract: In 1982, Marshall P. Tulin published a report proposing a framework for reducing the equations for gravity waves generated by moving bodies into a single nonlinear differential equation solvable in closed form [Proc. 14th Symp. on Naval Hydrodynamics, 1982, pp.19-51]. Several new and puzzling issues were highlighted by Tulin, notably the existence of weak and strong wave-making regimes, and the paradoxical fact that the theory seemed to be applicable to flows at low speeds, "but not too low speeds". These important issues were left unanswered, and despite the novelty of the ideas, Tulin's report fell into relative obscurity. Now thirty years later, we will revive Tulin's observations, and explain how an asymptotically consistent framework allows us to address these concerns. Most notably, we will explain, using the asymptotic method of steepest descents, how the production of free-surface waves can be related to the arrangement of integration contours connected to the shape of the moving body. This approach provides an intuitive and visual procedure for studying nonlinear wave-body interactions.


title: On the illusory connection between cross-conjugation and quantum interference

abstract: Quantum interference, be it destructive or constructive, has a substantial influence on the magnitude of molecular conductance and consequently there is significant interest in predicting these effects. It is commonly thought that cross-conjugated paths result in suppressed conductance due to destructive quantum interference. Using H\"uckel theory and DFT calculations we investigate systems that break this cross-conjugation rule of thumb. We predict and rationalize how a class of conjugated molecules containing closed loops can exhibit destructive interference despite being linearly conjugated and exhibit constructive interference despite being cross-conjugated. The arguments build on the graphical rules derived by Markussen et al. [1] and the hitherto neglected effects of closed loops in the molecular structure. Finally, we identify the 1,3 connected Azulene molecule as belonging to the closed-loop class and argue that this explains recent measurements of its electrical conductance.


title: Luminescent Ge-related centre in high-pressure synthesized diamond

abstract: We report on the high-pressure synthesis of novel nano- and microcrystalline high-quality diamonds with luminescent Ge-related centers. Observation of the four-line fine structure in luminescence at 2 eV (602 nm) at temperatures below 80 K manifests a high quality of diamonds. We demonstrate germanium and carbon isotope shifts in the fine structure of luminescence at 602 nm and its vibrational sideband which allows us to unambiguously associate the center with the germanium impurity entering into the diamond lattice. We show that there are two ground-state energy levels with the separation of 0.7 meV and two excited-state levels separated by 4.6 meV in the electronic structure of the center and suggest a split-vacancy structure of this center.


title: Two-Locus Likelihoods under Variable Population Size and Fine-Scale Recombination Rate Estimation

abstract: Two-locus sampling probabilities have played a central role in devising an efficient composite likelihood method for estimating fine-scale recombination rates. Due to mathematical and computational challenges, these sampling probabilities are typically computed under the unrealistic assumption of a constant population size, and simulation studies have shown that resulting recombination rate estimates can be severely biased in certain cases of historical population size changes. To alleviate this problem, we develop here new methods to compute the sampling probability for variable population size functions that are piecewise constant. Our main theoretical result, implemented in a new software package called LDpop, is a novel formula for the sampling probability that can be evaluated by numerically exponentiating a large but sparse matrix. This formula can handle moderate sample sizes ($n \leq 50$) and demographic size histories with a large number of epochs ($\mathcal{D} \geq 64$). In addition, LDpop implements an approximate formula for the sampling probability that is reasonably accurate and scales to hundreds in sample size ($n \geq 256$). Finally, LDpop includes an importance sampler for the posterior distribution of two-locus genealogies, based on a new result for the optimal proposal distribution in the variable-size setting. Using our methods, we study how a sharp population bottleneck followed by rapid growth affects the correlation between partially linked sites. Then, through an extensive simulation study, we show that accounting for population size changes under such a demographic model leads to substantial improvements in fine-scale recombination rate estimation. LDpop is freely available for download at https://github.com/popgenmethods/ldpop


title: Quasars Probing Quasars VIII. The Physical Properties of the Cool Circumgalactic Medium Surrounding z ~ 2-3 Massive Galaxies Hosting Quasars

abstract: We characterize the physical properties of the cool T ~10^4 K circumgalactic medium surrounding z ~2-3 quasar host galaxies, which are predicted to evolve into present day massive ellipticals. Using a statistical sample of 14 quasar pairs with projected separation < 300 kpc and high dispersion, high S/N spectra, we find extreme kinematics with low metal ion lines typically spanning ~ 500 km/s, exceeding any previously studied galactic population. The CGM is significantly enriched, even beyond the virial radius, with a median metallicity [M/H] ~ -0.6. The alpha/Fe abundance ratio is enhanced, suggesting that halo gas is primarily enriched by core-collapse supernovae. The projected cool gas mass within the virial radius is estimated to be 1.9*10^11 M_sun (R_\perp/160 kpc)^2, accounting for ~ 1/3 of the galaxy halo baryonic budget. The ionization state of CGM gas increases with projected distance from the foreground quasars, contrary to expectation if the quasar dominates the ionizing radiation flux. However, we also found peculiarities not exhibited in the CGM of other galaxy populations. In one absorption system, we may be detecting unresolved fluorescent Ly-alpha emission, and another system shows strong NV lines. Taken together these anomalies suggest that transverse sightlines are at least in some cases possibly illuminated. We also discovered a peculiar case where detection of the CII* fine structure line implies an electron density > 100 cm^-3 and subparsec scale gas clumps.


title: Comment on: Doubly hybrid density functional xDH-PBE0 from a parameter-free global hybrid model PBE0 (J. Chem. Phys. 136, 174103 (2012))

abstract: We have compared the performance of Grimme style DH/DSD and Zhang-Xu-Goddard type xDH/xDSD forms for double hybrids. In the DH and DSD forms, KS orbitals with elevated HF exchange and damped DFT correlation are used, while in the xDH and xDSD forms, the KS orbitals are obtained from a conventional hybrid functional with undamped DFT correlation. Generally, the difference in performance between DSD and xDSD functionals is small, slightly favoring xDSD. Augmentation of the xDH form with either same-spin MP2 correlation or a dispersion correction markedly improves performance. Best xDSD results appear to be obtained for orbitals obtained with `exact exchange' fractions in the 50-70% range. The orbitals for xDSD appear to be fairly transferable between different correlation functionals.


title: Petri Nets and its Polynomials

abstract: For every finite Petri net, we construct a commutative polynomial in two variables and with coefficients from the semiring of natural numbers. We also present an inverse construction and show that multiplication of polynomials correspondence to the product of the corresponding Petri nets in the category of Petri nets with Winskel's morphisms. We endow the set of all Petri nets with Zariski topology.


title: Lightning Strikes and Attribution of Climatic Change

abstract: Using lightning strikes as an example, two possible schemes are discussed for the attribution of changes in event frequency to climate change, and estimating the cost associated with them. The schemes determine the fraction of events that should be attributed to climatic change, and the fraction that should be attributed to natural chance. They both allow for the expected increase in claims and the fluctuations about this expected value. Importantly, the attribution fraction proposed in the second of these schemes is necessarily different to that found in epidemiological studies. This ensures that the statistically expected fraction of attributed claims is correctly equal to the expected increase in claims. The analysis of lightning data highlights two particular difficulties with data-driven, as opposed to modeled, attribution studies. The first is the possibility of unknown "confounding" variables that can influence the strike frequency. This is partly accounted for here by considering the influence of temperature changes within a given month, so as to standardise the data to allow for cyclical climatic influences. The second is the possibility suggested by the data presented here, that climate change may lead to qualitatively different climate patterns, with a different relationship between e.g. strike frequency and temperature.


title: Sequences from zero entropy noncommutative toral automorphisms and Sarnak Conjecture

abstract: In this paper we study $C^*$-algebra version of Sarnak Conjecture for noncommutative toral automorphisms. Let $A_\Theta$ be a noncommutative torus and $\alpha_\Theta$ be the noncommutative toral automorphism arising from a matrix $S\in GL(d,\mathbb{Z})$. We show that if the Voiculescu-Brown entropy of $\alpha_{\Theta}$ is zero, then the sequence $\{\rho(\alpha_{\Theta}^nu)\}_{n\in \mathbb{Z}}$ is a sum of a nilsequence and a zero-density-sequence, where $u\in A_\Theta$ and $\rho$ is any state on $A_\Theta$. Then by a result of Green and Tao, this sequence is linearly disjoint from the M\"obius function.


title: Comments on complete actions for open superstring field theory

abstract: We clarify a Wess-Zumino-Wtten-like structure including Ramond fields and propose one systematic way to construct gauge invariant actions: Wess-Zumino-Witten-like complete action $S_{\rm WZW}$. We show that Kunitomo-Okawa's action proposed in arXiv:1508.00366 can obtain a topological parameter dependence of Ramond fields and belongs to our WZW-like framework. In this framework, once a WZW-like functional $\mathcal{A}_{\eta } = \mathcal{A}_{\eta } [\Psi ]$ of a dynamical string field $\Psi $ is constructed, we obtain one realization of $S_{\rm WZW}$ parametrized by $\Psi $. On the basis of this way, we construct an action $\widetilde{S}$ whose on-shell condition is equivalent to the Ramond equations of motion proposed in arXiv:1506.05774. Using these results, we provide the equivalence of two these theories: arXiv:1508.00366 and arXiv:1506.05774.


title: Robust Semi-Supervised Classification for Multi-Relational Graphs

abstract: Graph-regularized semi-supervised learning has been used effectively for classification when (i) instances are connected through a graph, and (ii) labeled data is scarce. If available, using multiple relations (or graphs) between the instances can improve the prediction performance. On the other hand, when these relations have varying levels of veracity and exhibit varying relevance for the task, very noisy and/or irrelevant relations may deteriorate the performance. As a result, an effective weighing scheme needs to be put in place. In this work, we propose a robust and scalable approach for multi-relational graph-regularized semi-supervised classification. Under a convex optimization scheme, we simultaneously infer weights for the multiple graphs as well as a solution. We provide a careful analysis of the inferred weights, based on which we devise an algorithm that filters out irrelevant and noisy graphs and produces weights proportional to the informativeness of the remaining graphs. Moreover, the proposed method is linearly scalable w.r.t. the number of edges in the union of the multiple graphs. Through extensive experiments we show that our method yields superior results under different noise models, and under increasing number of noisy graphs and intensity of noise, as compared to a list of baselines and state-of-the-art approaches.


title: Probing Final Stages of Stellar Evolution with X-Ray Observations of SN 2013ej

abstract: Massive stars shape their surroundings with mass loss from winds during their lifetimes. Fast ejecta from supernovae, from these massive stars, shocks this circumstellar medium. Emission generated by this interaction provides a window into the final stages of stellar evolution, by probing the history of mass loss from the progenitor. Here we use Chandra and Swift x-ray observations of the type II-P/L SN 2013ej to probe the history of mass loss from its progenitor. We model the observed x-rays as emission from both heated circumstellar matter and supernova ejecta. The circumstellar density profile probed by the supernova shock reveals a history of steady mass loss during the final 400 years. The inferred mass loss rate of $3 \times 10^{-6} {\rm \; M_\odot \; yr^{-1}}$ points back to a 14 $M_\odot$ progenitor. Soon after the explosion we find significant absorption of reverse shock emission by a cooling shell. The column depth of this shell observed in absorption provides an independent and consistent measurement of the circumstellar density seen in emission. We also determine the efficiency of cosmic ray acceleration from x-rays produced by Inverse Compton scattering of optical photons by relativistic electrons. Only about 1 percent of the thermal energy is used to accelerate electrons. Our x-ray observations and modeling provides stringent tests for models of massive stellar evolution and micro-physics of shocks.


title: Cosmological constraints to dark matter with two- and many-body decays

abstract: We present a study of cosmological implications of generic dark matter decays. We consider two-body and many-body decaying scenarios. In the two-body case the massive particle has a possibly relativistic kick velocity and thus possesses a dynamical equation of state. This has implications to the expansion history of the universe. We use recent observational data from the cosmic microwave background, baryon acoustic oscillations and supernovae Type Ia to obtain constraints on the lifetime of the dark matter particle. We find that for an energy splitting where more than 40% of the dark matter particle energy is transferred to massless, relativistic particles in the two-body case, or more than 50% in the many-body case, lifetimes less than the age of the universe are excluded at more than 95% confidence. When the energy splitting falls to 10% the lifetime is constrained to be more than roughly half the age.


title: An Observational Perspective of the IMF: Progress and Challenges

abstract: The stellar initial mass function (IMF) is a fundamental astrophysical quantity that impacts a wide range of astrophysical problems from heavy element distribution to galactic evolution to planetary system formation. However, the origin and universality of the IMF are hotly debated both observationally and theoretically. I review recent observations of the IMF across a variety of environments. These suggest the IMF is surprisingly invariant between star-forming regions, star clusters, and spiral galaxies but that it may also vary under extreme conditions, including within the Galactic center and early type galaxies.


title: Simulated {\Lambda}CDM analogues of the thin Plane of Satellites around the Andromeda galaxy are not kinematically coherent structures

abstract: A large fraction of the dwarf satellites orbiting the Andromeda galaxy are surprisingly aligned in a thin, extended and apparently kinematically coherent planar structure. Such a structure is not easily found in simulations based on the Cold Dark Matter model ({\Lambda}CDM). Using 21 high resolution cosmological simulations we analyse the kinematics of planes of satellites similar to the one around Andromeda. We find good agreement when co-rotation is characterized by the line-of-sight velocity. At the same time, when co-rotation is inferred by the angular momenta of the satellites, the planes are in agreement with the plane around our Galaxy. We find such planes to be common in our high concentration haloes. The number of co-rotating satellites obtained from the sign of the line-of-sight velocity shows large variations depending on the viewing angle and is consistent with that obtained from a sample with random velocities. We find that the clustering of angular momentum vectors of the satellites in the plane is a better measure of the kinematic coherence. Thus we conclude that the line-of- sight velocity is not well suited as a proxy for the kinematical coherence of the plane. Analysis of the kinematics of our planes shows a fraction of $\sim$30% chance aligned satellites. Tracking the satellites in the plane back in time reveals that these planes are a transient feature and not kinematically coherent as would appear at first sight. Thus we expect some of the satellites in the plane around Andromeda to have high velocities perpendicular to the plane.


title: Tidal capture formation of Low Mass X-Ray Binaries from wide binaries in the field

abstract: We present a potentially efficient dynamical formation scenario for Low Mass X-ray Binaries (LMXBs) in the field, focusing on black-hole (BH) LMXBs. In this formation channel LMXBs are formed from wide binaries $(>1000$ AU) with a BH component and a stellar companion. The wide binary is perturbed by fly-by's of field stars and its orbit random-walks and changes over time. This diffusion process can drive the binary into a sufficiently eccentric orbit such that the binary components tidally interact at peri-center and the binary evolves to become a short period binary, which eventually evolves into an LMXB. The formation rate of LMXBs through this channel mostly depends on the number of such BH wide binaries progenitors, which in turn depends on the velocity kicks imparted to BHs (or NSs) at birth. We consider several models for the formation and survival of such wide binaries, and calculate the LMXB formation rates for each model. We find that models where BHs form through direct collapse with no/little natal kicks can give rise to high formation rates comparable with those inferred from observations. This formation scenario had several observational signatures: (1) the number density of LMXBs generally follows the background stellar density and (2) the mass function of the BH stellar companion should be comparable to the mass function of the background stellar population, likely peaking at $0.4-0.6$ M$_{\odot}$. The latter aspect, in particular, is unique to this model compared with previously suggested LMXB formation models following common envelope binary stellar evolution. We note that NS LMXBs can similarly form from wide binaries, but their formation rate through this channel is likely significantly smaller due to their much higher natal kicks.


title: Initial conditions of formation of starburst clusters: constraints from stellar dynamics

abstract: How starburst clusters form out of molecular clouds is still an open question. In this article, I highlight some of the key constraints in this regard, that one can get from the dynamical evolutionary properties of dense stellar systems. I particularly focus on secular expansion of massive star clusters and hierarchical merging of sub-clusters, and discuss their implications vis-a-vis the observed properties of young massive clusters. The analysis suggests that residual gas expulsion is necessary for shaping these clusters as we see them today, irrespective of their monolithic or hierarchical mode of formation.


title: Optical spectroscopy of the high-mass gamma-ray binary 1FGL J1018.6-5856: A probable neutron star primary

abstract: We present medium-resolution optical spectroscopy with the SOAR telescope of the O star secondary of the high-mass gamma-ray binary 1FGL J1018.6-5856 to help determine whether the primary is a neutron star or black hole. We find that the secondary has a low radial velocity semi-amplitude of 11-12 km/s, with consistent values obtained for H and He absorption lines. This low value strongly favors a neutron star primary: while a black hole cannot be excluded if the system is close to face on, such inclinations are disallowed by the observed rotation of the secondary. We also find the high-energy (X-ray and gamma-ray) flux maxima occur when the star is behind the compact object along our line of sight, inconsistent with a simple model of anisotropic inverse Compton scattering for the gamma-ray photons.


title: Self-dual Quantum Electrodynamics as Boundary State of the three dimensional Bosonic Topological Insulator

abstract: Inspired by the recent developments of constructing novel Dirac liquid boundary states of the $3d$ topological insulator, we propose one possible $2d$ boundary state of the $3d$ bosonic symmetry protected topological state with $U(1)_e \rtimes Z_2^T \times U(1)_s $ symmetry. This boundary theory is described by a $(2+1)d$ quantum electrodynamics (QED$_3$) with two flavors of Dirac fermions ($N_f = 2$) coupled with a noncompact U(1) gauge field: $ \mathcal{L} = \sum_{j = 1}^2 \bar{\psi}_j \gamma_\mu (\partial_\mu - i a_\mu) \psi_j - i A^{s}_\mu \bar{\psi_i} \gamma_\mu \tau^z_{ij} \psi_j + \frac{i}{2\pi} \epsilon_{\mu\nu\rho} a_\mu \partial_\nu A^{e}_\rho $, where $a_\mu$ is the internal noncompact U(1) gauge field, $A^s_\mu$ and $A^e_\mu$ are two external gauge fields that couple to $U(1)_s$ and $U(1)_e$ global symmetries respectively. We demonstrate that this theory has a "self-dual" structure, which is a fermionic analogue of the self-duality of the noncompact CP$^1$ theory with easy plane anisotropy. Under the self-duality, the boundary action takes exactly the same form except for an exchange between $A^s_\mu$ and $A^e_\mu$. The self-duality may still hold after we break one of the U(1) symmetries (which makes the system a bosonic topological insulator), with some subtleties that will be discussed.


title: Intrinsic Diophantine approximation in Carnot groups and in the Siegel model of the Heisenberg group

abstract: We initiate the study of an intrinsic notion of Diophantine approximation on a rational Carnot group $G$. If $G$ has Hausdorff dimension $Q$, we show that its Diophantine exponent is equal to $(Q+1)/Q$, generalizing the case $G=\mathbb R^n$. We furthermore obtain a precise asymptotic on the count of rational approximations. We then focus on the case of the Heisenberg group $\mathbb H^n$, distinguishing between two notions of Diophantine approximation by rational points in $\mathbb H^n$: Carnot Diophantine approximation and Siegel Diophantine approximation. After computing the Siegel Diophantine exponent (surprisingly, equal to 1 for all $\mathbb H^n$), we consider Siegel-badly approximable points to show that Siegel approximation is linked to both Heisenberg continued fractions and to geodesics in complex hyperbolic space. We conclude by showing that Carnot and Siegel approximation are qualitatively different: Siegel-badly approximable points are Schmidt winning in any complete Ahlfors regular subset of $\mathbb H^n$, while the set of Carnot-badly approximable points does not have this property.


title: A new method to quantify the effects of baryons on the matter power spectrum

abstract: Future large-scale galaxy surveys have the potential to become leading probes for cosmology provided the influence of baryons on the total mass distribution is understood well enough. As hydrodynamical simulations strongly depend on details in the feedback implementations, no unique and robust predictions for baryonic effects currently exist. In this paper we propose a baryonic correction model that modifies the density field of dark-matter-only $N$-body simulations to mimic the effects of baryons from any underlying adopted feedback recipe. The model assumes haloes to consist of 4 components: 1- hot gas in hydrostatical equilibrium, 2- ejected gas from feedback processes, 3- central galaxy stars, and 4- adiabatically relaxed dark matter, which all modify the initial dark-matter-only density profiles. These altered profiles allow to define a displacement field for particles in $N$-body simulations and to modify the total density field accordingly. The main advantage of the baryonic correction model is to connect the total matter density field to the observable distribution of gas and stars in haloes, making it possible to parametrise baryonic effects on the matter power spectrum. We show that the most crucial quantities are the mass fraction of ejected gas and its corresponding ejection radius. The former controls how strongly baryons suppress the power spectrum, while the latter provides a measure of the scale where baryonic effects become important. A comparison with X-ray and Sunyaev-Zel'dovich cluster observations suggests that baryons suppress wave modes above $k\sim0.5$ h/Mpc with a maximum suppression of 10-25 percent around $k\sim 2$ h/Mpc. More detailed observations of the gas in the outskirts of groups and clusters are required to decrease the large uncertainties of these numbers.


title: Aperiodic Weak Topological Superconductors

abstract: Weak topological phases are usually described in terms of protection by the lattice translation symmetry. Their characterization explicitly relies on periodicity since weak invariants are expressed in terms of the momentum-space torus. We prove the compatibility of weak topological superconductors with aperiodic systems, such as quasicrystals. We go beyond usual descriptions of weak topological phases and introduce a novel, real-space formulation of the weak invariant, based on the Clifford pseudospectrum. A non-trivial value of this index implies a non-trivial bulk phase, which is robust against disorder and hosts localized zero-energy modes at the edge. Our recipe for determining the weak invariant is directly applicable to any finite-sized system, including disordered lattice models. This direct method enables a quantitative analysis of the level of disorder the topological protection can withstand.


title: Semicontinuity of 4d N=2 spectrum under renormalization group flow

abstract: We study renormalization group flow of four dimensional N=2 SCFTs defined by isolated hypersurface three-fold singularities. We define the spectrum of N=2 theory as the set of scaling dimensions of the parameters on the Coulomb branch, which include Coulomb branch moduli, mass parameters and coupling constants. We prove that the spectrum of those theories is semicontinous under the RG flow on the Coulomb branch using the mathematical result about the singularity spectra under deformation. The semicontinuity behavior of N=2 spectrum implies a theorem under relevant and Coulomb branch moduli deformation, the absence of dangerous irrelevant deformations and can be taken as the necessary condition for the ending point of a RG flow. This behavior is also true for (c,c) ring deformation of two dimensional Landau-Ginzburg model with (2,2) supersymmetry.


title: The sign of active galactic nucleus quenching in a merger remnant with radio jets

abstract: We investigate optical, infrared, and radio active galactic nucleus (AGN) signs in the merger remnant Arp 187, which hosts luminous jets launched in the order of $10^5$ yr ago but whose present-day AGN activity is still unknown. We find AGN signs from the optical BPT diagram and infrared [OIV]25.89 $\mu$m line, originating from the narrow line regions of AGN. On the other hand, Spitzer/IRS show the host galaxy dominated spectra, suggesting that the thermal emission from the AGN torus is considerably small or already diminished. Combining the black hole mass, the upper limit of radio luminosity of the core, and the fundamental plane of the black hole enable us to estimate X-ray luminosity, which gives $<10^{40}$ erg s$^{-1}$. Those results suggest that the AGN activity of Arp 187 has already been quenched, but the narrow line region is still alive owing to the time delay of emission from the past AGN activity.


title: The SAMI Galaxy Survey: Can we trust aperture corrections to predict star formation?

abstract: In the low redshift Universe (z<0.3), our view of galaxy evolution is primarily based on fibre optic spectroscopy surveys. Elaborate methods have been developed to address aperture effects when fixed aperture sizes only probe the inner regions for galaxies of ever decreasing redshift or increasing physical size. These aperture corrections rely on assumptions about the physical properties of galaxies. The adequacy of these aperture corrections can be tested with integral-field spectroscopic data. We use integral-field spectra drawn from 1212 galaxies observed as part of the SAMI Galaxy Survey to investigate the validity of two aperture correction methods that attempt to estimate a galaxy's total instantaneous star formation rate. We show that biases arise when assuming that instantaneous star formation is traced by broadband imaging, and when the aperture correction is built only from spectra of the nuclear region of galaxies. These biases may be significant depending on the selection criteria of a survey sample. Understanding the sensitivities of these aperture corrections is essential for correct handling of systematic errors in galaxy evolution studies.


title: Roy-Steiner-equation analysis of pion-nucleon scattering

abstract: We review the structure of Roy-Steiner equations for pion-nucleon scattering, the solution for the partial waves of the t-channel process $\pi\pi\to \bar N N$, as well as the high-accuracy extraction of the pion-nucleon S-wave scattering lengths from data on pionic hydrogen and deuterium. We then proceed to construct solutions for the lowest partial waves of the s-channel process $\pi N\to \pi N$ and demonstrate that accurate solutions can be found if the scattering lengths are imposed as constraints. Detailed error estimates of all input quantities in the solution procedure are performed and explicit parameterizations for the resulting low-energy phase shifts as well as results for subthreshold parameters and higher threshold parameters are presented. Furthermore, we discuss the extraction of the pion-nucleon $\sigma$-term via the Cheng-Dashen low-energy theorem, including the role of isospin-breaking corrections, to obtain a precision determination consistent with all constraints from analyticity, unitarity, crossing symmetry, and pionic-atom data. We perform the matching to chiral perturbation theory in the subthreshold region and detail the consequences for the chiral convergence of the threshold parameters and the nucleon mass.


title: A Bayesian Calibration-Prediction Method for Reducing Model-Form Uncertainties with Application in RANS Simulations

abstract: Model-form uncertainties in complex mechanics systems are a major obstacle for predictive simulations. Reducing these uncertainties is critical for stake-holders to make risk-informed decisions based on numerical simulations. For example, Reynolds-Averaged Navier-Stokes (RANS) simulations are increasingly used in mission-critical systems involving turbulent flows. However, for many practical flows the RANS predictions have large model-form uncertainties originating from the uncertainty in the modeled Reynolds stresses. Recently, a physics-informed Bayesian framework has been proposed to quantify and reduce model-form uncertainties in RANS simulations by utilizing sparse observation data. However, in the design stage of engineering systems, measurement data are usually not available. In the present work we extend the original framework to scenarios where there are no available data on the flow to be predicted. In the proposed method, we first calibrate the model discrepancy on a related flow with available data, leading to a statistical model for the uncertainty distribution of the Reynolds stress discrepancy. The obtained distribution is then sampled to correct the RANS-modeled Reynolds stresses for the flow to be predicted. The extended framework is a Bayesian calibration-prediction method. The merits of the proposed method are demonstrated on two flows that are challenging to standard RANS models. By not requiring observation data on the flow to be predicted, the present calibration-prediction method will gain wider acceptance in practical engineering application. While RANS modeling is chosen to demonstrate the merits of the proposed framework, the methodology is generally applicable to other complex mechanics models involving solids, fluids flows, or the coupling between the two, where model-form uncertainties are present in the constitutive relations.


title: Towards higher dimensional black rings

abstract: Despite all attempts, exact solutions for black rings in more than five dimensions remain elusive. In this paper we clarify some of the reasons for that, in particular we show that a peculiar symmetry of the five--dimensional black ring - separability of the base - cannot occur in dimensions higher than five. We also construct supersymmetric solutions that have symmetries of 5D supersymmetric black ring and show that they do not have regular horizons.


title: The Quest for B Modes from Inflationary Gravitational Waves

abstract: The search for the curl component (B mode) in the cosmic microwave background (CMB) polarization induced by inflationary gravitational waves is described. The canonical single-field slow-roll model of inflation is presented, and we explain the quantum production of primordial density perturbations and gravitational waves. It is shown how these gravitational waves then give rise to polarization in the CMB. We then describe the geometric decomposition of the CMB polarization pattern into a curl-free component (E mode) and curl component (B mode) and show explicitly that gravitational waves induce B modes. We discuss the B modes induced by gravitational lensing and by Galactic foregrounds and show how both are distinguished from those induced by inflationary gravitational waves. Issues involved in the experimental pursuit of these B modes are described, and we summarize some of the strategies being pursued. We close with a brief discussion of some other avenues toward detecting/characterizing the inflationary gravitational-wave background.


title: Entropy-continuity for interval maps with holes

abstract: We study the dependence of the topological entropy of piecewise monotonic maps with holes under perturbations, for example sliding a hole of fixed size at uniform speed or expanding a hole with uniform expansion. We show that under suitable conditions the topological entropy varies locally Hoelder continuously with the local Hoelder exponent depending itself on the value of the topological entropy.


title: Evidence for Enrichment by Supernovae in the Globular Cluster NGC 6273

abstract: In our recent investigation (Lim et al. 2015), we have shown that narrow-band photometry can be combined with low-resolution spectroscopy to effectively search for globular clusters (GCs) with supernovae (SNe) enrichments. Here we apply this technique to the metal-poor bulge GC NGC 6273, and find that the red giant branch stars in this GC are clearly divided into two distinct subpopulations having different calcium abun- dances. The Ca rich subpopulation in this GC is also enhanced in CN and CH, showing a positive correlation between them. This trend is identical to the result we found in M22, suggesting that this might be a ubiquitous nature of GCs more strongly affected by SNe in their chemical evolution. Our results suggest that NGC 6273 was massive enough to retain SNe ejecta which would place this cluster in the growing group of GCs with Galactic building block characteristics, such as {\omega} Centauri and Terzan 5.


title: HE 0017+0055 : A probable pulsating CEMP-rs star and long-period binary

abstract: A radial-velocity monitoring of the Carbon-Enhanced Metal-Poor (CEMP) star HE 0017+0055 over 8 years with the Nordic Optical Telescope and Mercator telescopes reveals variability with a period of 384 d and amplitude of 540$\pm27$ m s$^{-1}$, superimposed on a nearly linear long-term decline of $\sim$1 m s$^{-1}$ day$^{-1}$. High-resolution HERMES/Mercator and Keck/HIRES spectra have been used to derive elemental abundances using 1-D LTE MARCS models. A metallicity of [Fe/H] $\sim -2.4$ is found, along with s-process overabundances on the order of 2 dex (with the exception of [Y/Fe] $\sim+0.5$), and most notably overabundances of r-process elements like Sm, Eu, Dy, and Er in the range 0.9 - 2.0 dex. With [Ba/Fe] $ > 1.9$ dex and [Eu/Fe] = 2.3 dex, HE 0017+0055 is a CEMP-rs star. It appears to be a giant star below the tip of the red giant branch (RGB). The s-process pollution must therefore originate from mass transfer from a companion formerly on the AGB, now a carbon-oxygen white dwarf (WD). If the 384 d velocity variations are attributed to the WD companion, its orbit must be seen almost face-on, with $i \sim 2.3^\circ$, because the mass function is very small: $f(M_1,M_2) = (6.1\pm1.1)\times10^{-6}$ Msun. Alternatively, the WD orbital motion could be responsible for the long-term velocity variations, with a period of several decades. The 384 d variations should then be attributed either to a low-mass inner companion (perhaps a brown dwarf, depending on the orbital inclination), or to stellar pulsations. The latter possibility is made likely by the fact that similar low-amplitude velocity variations, with periods close to 1 yr, have been reported for other CEMP stars in a companion paper (Jorissen et al., 2015). A definite conclusion about the origin of the 384 d velocity variations should however await the detection of synchronous low-amplitude photometric variations.


title: Nonlinear stochastic heat equation driven by spatially colored noise: moments and intermittency

abstract: In this paper, we study the stochastic heat equation in the spatial domain $\mathbb{R}^d$ subject to a Gaussian noise which is white in time and colored in space. The spatial correlation can be any symmetric, nonnegative and nonnegative-definite function that satisfies {\it Dalang's condition}. We establish the existence and uniqueness of a random field solution starting from measure-valued initial data. We find the upper and lower bounds for the second moment. As a first application of these moments bounds, we find the necessary and sufficient conditions for the solution to have phase transition for the second moment Lyapunov exponents. As another application, we prove a localization result for the intermittency fronts.


title: Polarization microlensing in the quadruply imaged broad absorption line quasar H1413+117

abstract: We have obtained spectropolarimetric observations of the four images of the gravitationally lensed broad absorption line quasar H1413+117. The polarization of the microlensed image D is significantly different both in the continuum and in the broad lines from the polarization of image A which is essentially unaffected by microlensing. The observations suggest that the continuum is scattered off two regions spatially separated and producing roughly perpendicular polarizations. These results are compatible with a model in which the microlensed polarized continuum comes from a compact region located in the equatorial plane close to the accretion disk and the non-microlensed continuum from an extended region located along the polar axis.


title: Cubic Single Crystal Representations in Classical and Size-dependent Couple Stress Elasticity

abstract: Beginning with Cosserat theory in the early 20th century, there have been several different formulations for size-dependent elastic response. In this paper, we concentrate on the application of classical Cauchy theory and the recent parsimonious consistent couple stress theory to model a homogeneous linear elastic solid, exemplified by a pure single crystal with cubic structure. The focus is on an examination of elastodynamic response based upon wave velocities from ultrasonic excitation and phonon dispersion curves, along with adiabatic bulk moduli measurements. In particular, we consider in detail elastic parameter estimation within classical elasticity and consistent couple stress theory for four different cubic single crystals (NaCl, KCl, Cu, CuZn). The classical theory requires the estimation of three independent material parameters, while only one additional parameter relating skew-symmetric mean curvature to skew-symmetric couple-stress is needed for the size-dependent consistent couple stress theory. This additional parameter can be defined for cubic crystals in terms of a material length scale, which is found to be on the order of tens of microns for the four materials studied here. Furthermore, a detailed statistical investigation provides strong to very strong evidence that couple stress theory is superior to classical Cauchy elasticity for representing the wave velocities and adiabatic bulk moduli for all four single crystals.


title: What can we learn from the dynamics of entanglement and quantum discord in the Tavis-Cummings model?

abstract: We revisit the problem of the dynamics of quantum correlations in the exact Tavis-Cummings model. We show that many of the dynamical features of quantum discord attributed to dissipation are already present in the exact framework and are due to the well known non-linearities in the model and to the choice of initial conditions. Through a comprehensive analysis, supported by explicit analytical calculations, we find that the dynamics of entanglement and quantum discord are far from being trivial or intuitive. In this context, we find states that are indistinguishable from the point of view of entanglement and distinguishable from the point of view of quantum discord, states where the two quantifiers give opposite information and states where they give roughly the same information about correlations at a certain time. Depending on the initial conditions, this model exhibits a fascinating range of phenomena that can be used for experimental purposes such as: Robust states against change of manifold or dissipation, tunable entanglement states and states with a counterintuitive sudden birth as the number of photons increase. We furthermore propose an experiment called quantum discord gates where discord is zero or non-zero depending on the number of photons.


title: Volume and Determinant Densities of Hyperbolic Rational Links

abstract: The volume density of a hyperbolic link is defined as the ratio of hyperbolic volume to crossing number. We study its properties and a closely-related invariant called the determinant density. It is known that the sets of volume densities and determinant densities of links are dense in the interval [0,v_{oct}]. We construct sequences of alternating knots whose volume and determinant densities both converge to any x in [0,v_{oct}]. We also investigate the distributions of volume and determinant densities for hyperbolic rational links, and establish upper bounds and density results for these invariants.


title: The Complexity of Pattern Matching for $321$-Avoiding and Skew-Merged Permutations

abstract: The Permutation Pattern Matching problem, asking whether a pattern permutation $\pi$ is contained in a permutation $\tau$, is known to be NP-complete. In this paper we present two polynomial time algorithms for special cases. The first algorithm is applicable if both $\pi$ and $\tau$ are $321$-avoiding; the second is applicable if $\pi$ and $\tau$ are skew-merged. Both algorithms have a runtime of $O(kn)$, where $k$ is the length of $\pi$ and $n$ the length of $\tau$.


title: A method to calculate Franck-Condon factors in terms of the tomographic probability representation

abstract: We introduce a new method to calculate Franck-Condon factors in polyatomic molecules that is based on the tomographic probability approach to quantum mechanics. This approach is implemented to calculate transition probabilities in various systems under an instantaneous change of frequency and equilibrium position of nuclei in a molecule by an external force. The problem is considered for different types of the Dushinsky matrix and for any quantity of atoms in a molecule.


title: Scalable posterior approximations for large-scale Bayesian inverse problems via likelihood-informed parameter and state reduction

abstract: Two major bottlenecks to the solution of large-scale Bayesian inverse problems are the scaling of posterior sampling algorithms to high-dimensional parameter spaces and the computational cost of forward model evaluations. Yet incomplete or noisy data, the state variation and parameter dependence of the forward model, and correlations in the prior collectively provide useful structure that can be exploited for dimension reduction in this setting--both in the parameter space of the inverse problem and in the state space of the forward model. To this end, we show how to jointly construct low-dimensional subspaces of the parameter space and the state space in order to accelerate the Bayesian solution of the inverse problem. As a byproduct of state dimension reduction, we also show how to identify low-dimensional subspaces of the data in problems with high-dimensional observations. These subspaces enable approximation of the posterior as a product of two factors: (i) a projection of the posterior onto a low-dimensional parameter subspace, wherein the original likelihood is replaced by an approximation involving a reduced model; and (ii) the marginal prior distribution on the high-dimensional complement of the parameter subspace. We present and compare several strategies for constructing these subspaces using only a limited number of forward and adjoint model simulations. The resulting posterior approximations can rapidly be characterized using standard sampling techniques, e.g., Markov chain Monte Carlo. Two numerical examples demonstrate the accuracy and efficiency of our approach: inversion of an integral equation in atmospheric remote sensing, where the data dimension is very high; and the inference of a heterogeneous transmissivity field in a groundwater system, which involves a partial differential equation forward model with high dimensional state and parameters.


title: When is a network epidemic hard to eliminate?

abstract: We consider the propagation of a contagion process (epidemic) on a network and study the problem of dynamically allocating a fixed curing budget to the nodes of the graph, at each time instant. For bounded degree graphs, we provide a lower bound on the expected time to extinction under any such dynamic allocation policy, in terms of a combinatorial quantity that we call the resistance of the set of initially infected nodes, the available budget, and the number of nodes n. Specifically, we consider the case of bounded degree graphs, with the resistance growing linearly in n. We show that if the curing budget is less than a certain multiple of the resistance, then the expected time to extinction grows exponentially with n. As a corollary, if all nodes are initially infected and the CutWidth of the graph grows linearly, while the curing budget is less than a certain multiple of the CutWidth, then the expected time to extinction grows exponentially in n. The combination of the latter with our prior work establishes a fairly sharp phase transition on the expected time to extinction (sub-linear versus exponential) based on the relation between the CutWidth and the curing budget.


title: A lower bound on the performance of dynamic curing policies for epidemics on graphs

abstract: We consider an SIS-type epidemic process that evolves on a known graph. We assume that a fixed curing budget can be allocated at each instant to the nodes of the graph, towards the objective of minimizing the expected extinction time of the epidemic. We provide a lower bound on the optimal expected extinction time as a function of the available budget, the epidemic parameters, the maximum degree, and the CutWidth of the graph. For graphs with large CutWidth (close to the largest possible), and under a budget which is sublinear in the number of nodes, our lower bound scales exponentially with the size of the graph.


title: The Slice Spectral Sequence for certain $RO(C_{p^n})$-graded Suspensions of $H\underline{\mathbb Z}$

abstract: We study the slice filtration and associated spectral sequence for a family of $RO(C_{p^{n}})$-graded suspensions of the Eilenberg-MacLane spectrum for the constant Mackey functor $\underline{\mathbb Z}$. Since $H\underline{\mathbb Z}$ is the zero slice of the sphere spectrum, this begins an analysis of how one can describe the slices of a suspension in terms of the original slices.


title: Refined Tur\'an numbers and Ramsey numbers for the loose 3-uniform path of length three

abstract: Let $P$ denote a 3-uniform hypergraph consisting of 7 vertices $a,b,c,d,e,f,g$ and 3 edges $\{a,b,c\}, \{c,d,e\},$ and $\{e,f,g\}$. It is known that the $r$-color Ramsey number for $P$ is $R(P;r)=r+6$ for $r\le 7$. The proof of this result relies on a careful analysis of the Tur\'an numbers for $P$. In this paper, we refine this analysis further and compute, for all $n$, the third and fourth order Tur\'an numbers for $P$. With the help of the former, we confirm the formula $R(P;r)=r+6$ for $r\in\{8,9\}$.


title: Polarization dependent tuning of the Hanle effect in the ground state of Cs

abstract: We demonstrate that the Hanle effect can be tuned between magnetically induced absorption (MIA) and magnetically induced transmission (MIT) simply by changing the polarization of the input laser beam. The experiments are done using closed hyperfine transitions of the $ \rm D_2 $ line of ${\rm ^{133}Cs}$ ---$ F_g = 3 \rightarrow F_e = 2 $ and $ F_g =4 \rightarrow F_e = 5 $. The former shows a transformation from MIT to MIA, while the latter shows the opposite behavior. A qualitative explanation based on optical pumping and coherences among the magnetic sublevels of the ground state is borne out by a detailed density-matrix calculation. To increase the coherence time, the experiments are done in a Cs vapor cell with paraffin coating on the walls. The observed linewidth is extremely narrow ($\sim 0.1$ mG) compared to previous work in this area, making this a promising technique for all kinds of precision measurements.


title: IndIGO and LIGO-India: Scope and Plans for Gravitational Wave Research and Precision Metrology in India

abstract: Initiatives by the IndIGO (Indian Initiative in Gravitational Wave Observations) Consortium during the past three years have materialized into concrete plans and project opportunities for instrumentation and research based on advanced interferometer detectors . With the LIGO-India opportunity, this initiative has a taken a promising path towards significant participation in gravitational wave (GW) astronomy and research, and in developing and nurturing precision fabrication and measurement technologies in India. The proposed LIGO-India detector will foster integrated development of frontier GW research in India and will provide opportunity for substantial contributions to global GW research and astronomy. Widespread interest and enthusiasm about these developments in premier research and educational institutions in India lead to the expectation that there will be a grand surge of activity in precision metrology, instrumentation, data handling and computation etc. in the context of LIGO-India. I discuss the scope of such research in the backdrop of the IndIGO action plan and the LIGO-India project.


title: Staircase diagrams and enumeration of smooth Schubert varieties

abstract: We enumerate smooth and rationally smooth Schubert varieties in the classical finite types A, B, C, and D, extending Haiman's enumeration for type A. To do this enumeration, we introduce a notion of staircase diagrams on a graph. These combinatorial structures are collections of steps of irregular size, forming interconnected staircases over the given graph. Over a Dynkin-Coxeter graph, the set of "nearly-maximally labelled" staircase diagrams is in bijection with the set of Schubert varieties with a complete Billey-Postnikov (BP) decomposition. We can then use an earlier result of the authors showing that all finite-type rationally smooth Schubert varieties have a complete BP decomposition to finish the enumeration.


title: Rigidity and Curvature Estimates for Graphical Self-shrinkers

abstract: Self-shrinkers are hypersurfaces that shrink homothetically under mean curvature flow; these solitons model the singularities of the flow. It it presently known that an entire self-shrinking graph must be a hyperplane. In this paper we show that the hyperplane is rigid in an even stronger sense, namely: For $2 \leq n \leq 6$, any smooth, complete self-shrinker $\Sigma^n\subset\mathbf{R}^{n+1}$ that is graphical inside a large, but compact, set must be a hyperplane. In fact, this rigidity holds within a larger class of almost stable self-shrinkers. A key component of this paper is the procurement of linear curvature estimates for almost stable shrinkers, and it is this step that is responsible for the restriction on $n$. Our methods also yield uniform curvature bounds for translating solitons of the mean curvature flow.


title: Motion Blur Filtering: A Statistical Approach for Extracting Confinement Forces and Diffusivity from a Single Blurred Trajectory

abstract: Single Particle Tracking (SPT) can aid in understanding complex spatio-temporal processes. However, quantifying diffusivity and forces from individual live cell trajectories is complicated by inter- & intra-trajectory kinetic heterogeneity, thermal fluctuations, and statistical temporal dependence inherent to the underlying molecule's time correlated confined dynamics experienced in the cell. Experimental artifacts such as localization uncertainty and motion blur also obscure the data. We introduce a new maximum likelihood estimation (MLE) technique that decouples the above noise sources and systematically treats temporal correlation via a likelihood function (permitting more reliable extraction of effective forces from position vs. time data). Our estimator is demonstrated to be consistent over a wide range of exposure times, diffusion coefficients, and confinement "radii". The algorithm and corresponding software can reliably extract motion parameters independent of exposure time in trajectories exhibiting confined and/or non-stationary dynamics and will aid in directly comparing trajectories obtained from different imaging modalities.


title: Compact Picture in Extended Superconformal Field Theories

abstract: There is a complex conformal transformation, which maps the $D$ - dimensional real Minkowski space on a bounded set in the $D$ - dimensional complex vector space. It generalizes the Cayley map from $D=1$ dimensions to higher space-time dimensions. This transformation provides a very convenient coordinate picture for Conformal Field Theories called compact picture. In this paper we extend the compact picture coordinates for superconformal field theories in four space-time dimensions.


title: Forbidden integer ratios of consecutive power sums

abstract: Let $S_k(m):=1^k+2^k+\cdots+(m-1)^k$ denote a power sum. In 2011 Bernd Kellner formulated the conjecture that for $m\ge 4$ the ratio $S_k(m+1)/S_k(m)$ of two consecutive power sums is never an integer. We will develop some techniques that allow one to exclude many integers $\rho$ as a ratio and combine them to exclude the integers $3\le \rho\le 1501$ and, assuming a conjecture on irregular primes to be true, a set of density $1$ of ratios $\rho$. To exclude a ratio $\rho$ one has to show that the Erd\H{o}s-Moser type equation $(\rho-1)S_k(m)=m^k$ has no non-trivial solutions.


title: Compact Low-Voltage, High-Power, Multi-beam Klystron for ILC: Initial Test Results

abstract: Initial test results of an L-band multi-beam klystron with parameters relevant for ILC are presented. The chief distinction of this tube from MBKs already developed for ILC is its low operating voltage of 60 kV, a virtue that implies considerable technological simplifications in the accelerator complex. To demonstrate the concept underlying the tubes design, a six-beamlet quadrant (a 54 inch high one-quarter portion of the full 1.3 GHz tube) was built and recently underwent initial tests, with main goals of demonstrating rated gun perveance, rated gain, and at least one-quarter of the full 10-MW rated power. Our initial three-day conditioning campaign without RF drive (140 microsec pulses @ 60 Hz) was stopped at 53% of full rated duty because of time-limits at the test-site; no signs appeared that would seem to prevent achieving full duty operation (i.e., 1.6 msec pulses @ 10 Hz). The subsequent tests with 10-15 microsec RF pulses confirmed the rated gain, produced output powers of up to 2.86 MW at 60 kV with high efficiency and 56 dB gain, and showed acceptable beam interception. These results suggest that a full version of the tube should be able to produce up to 11.5 MW. Follow-on tests are planned for later in 2015.


title: The far-infrared emission of the radio-loud quasar 3C318

abstract: 3C318, a radio-loud quasar at z=1.574, is a subgalactic-sized radio source, and a good test-bed for the interplay between black hole and galaxy growth in the high-z Universe. Based on its IRAS, ISO, and SCUBA detections, it has long been considered as one of the most intrinsically luminous (L$_{\mathrm{IR}}$ > 10$^{13}$ L$_{\odot}$) infrared sources in the Universe. Recent far-infrared data from the Herschel Space Observatory reveal that most of the flux associated with 3C318 measured with earlier instruments in fact comes from a bright nearby source. Optical imaging and spectroscopy show that this infrared-bright source is a strongly star-forming pair of interacting galaxies at z=0.35. Adding existing Spitzer and SDSS photometry, we perform a spectral energy distribution analysis of the pair, and find that it has a combined infrared luminosity of L$_{\mathrm{IR}}$ = 1.5 $\times$ 10$^{12}$ L$_{\odot}$, comparable to other intermediate-redshift ultra-luminous infrared galaxies studied with Herschel. Isolating the emission from 3C318's host, we robustly constrain the level of star formation to a value a factor of three lower than that published earlier, which is more in line with the star formation activity found in other Herschel-detected 3CR objects at similar redshift.


title: Stability of the semi-tamed and tamed Euler schemes for stochastic differential equations with jumps under non-global Lipschitz condition

abstract: Under non-global Lipschitz condition, Euler Explicit method fails to converge strongly to the exact solution, while Euler implicit method converges but requires much computational efforts. Tamed scheme was first introduced in [2] to overcome this failure of the standard explicit method. This technique is extended to SDEs driven by Poisson jump in [3] where several schemes were analyzed. In this work, we investigate their nonlinear stability under non-global Lipschitz and their linear stability. Numerical simulations to sustain the theoretical results are also provided.


title: An apparatus for studying electrical breakdown in liquid helium at 0.4 K and testing electrode materials for the SNS nEDM experiment

abstract: We have constructed an apparatus to study DC electrical breakdown in liquid helium at temperatures as low as 0.4 K and at pressures between the saturated vapor pressure and $\sim$600 torr. The apparatus can house a set of electrodes that are 12 cm in diameter with a gap of $1-2$ cm between them, and a potential up to $\pm 50$ kV can be applied to each electrode. Initial results demonstrated that it is possible to apply fields exceeding 100 kV/cm in a 1 cm gap between two electropolished stainless steel electrodes 12 cm in diameter for a wide range of pressures at 0.4 K. We also measured the current between two electrodes. Our initial results, $I<1$ pA at 45 kV, correspond to a lower bound on the effective volume resistivity of LHe of $\rho_V > 5\times10^{18}$ $\Omega\cdot$cm. This lower bound is 5 times larger than the bound previously measured. We report the design, construction, and operational experience of the apparatus, as well as initial results.


title: The spontaneous $\mathbb{Z}_2$ breaking Twin Higgs

abstract: The Twin Higgs model seeks to address the little hierarchy problem by making the Higgs a pseudo-Goldstone of a global $SU(4)$ symmetry that is spontaneously broken to $SU(3)$. Gauge and Yukawa couplings, which explicitly break $SU(4)$, enjoy a discrete $\mathbb{Z}_2$ symmetry that accidentally maintains $SU(4)$ at the quadratic level and therefore keeps the Higgs light. Contrary to most beyond the Standard Model theories, the quadratically divergent corrections to the Higgs mass are cancelled by a mirror sector, which is uncharged under the Standard Model groups. However, the Twin Higgs with an exact $\mathbb{Z}_2$ symmetry leads to equal vevs in the Standard Model and mirror sectors, which is phenomenologically unviable. An explicit $\mathbb{Z}_2$ breaking potential must then be introduced and tuned against the $SU(4)$ breaking terms to produce a hierarchy of vevs between the two sectors. This leads to a moderate but non-negligible tuning. We propose a model to alleviate this tuning, without the need for an explicit $\mathbb{Z}_2$ breaking sector. The model consists of two $SU(4)$ fundamental Higgses, one whose vacuum preserves $\mathbb{Z}_2$ and one whose vacuum breaks it. As the interactions between the two Higgses are turned on, the $\mathbb{Z}_2$ breaking is transmitted from the broken to the unbroken sector and a small hierarchy of vevs is naturally produced. The presence of an effective tadpole and feedback between the two Higgses lead to a sizable improvement of the tuning. The resulting Higgs boson is naturally very Standard Model like.


title: Jacobi no-core shell model for p-shell nuclei

abstract: We introduce an algorithm to obtain coefficients of fractional parentage for light $p$-shell nuclei. The coefficients enable to use Jacobi coordinates in no-core shell model calculations separating off the center-of-mass motion. Fully antisymmetrized basis states are given together with recoupling coefficients that allow one to apply two- and three-nucleon operators. As an example, we study the dependence on the harmonic oscillator frequency of $^3$H, $^4$He, $^6$He, $^6$Li and $^7$Li and extract their binding and excitation energies. The coefficients will be made openly accessible as HDF5 data files.


title: Augmentations and Rulings of Legendrian Links in $\#^k(S^1\times S^2)$

abstract: Given a Legendrian link in $\#^k(S^1\times S^2)$, we extend the definition of a normal ruling from $J^1(S^1)$ given by Lavrov and Rutherford and show that the existence of an augmentation to any field of the Chekanov-Eliashberg differential graded algebra over $\mathbb{Z}[t,t^{-1}]$ is equivalent to the existence of a normal ruling of the front diagram. For Legendrian knots, we also show that any even graded augmentation must send $t$ to $-1$. We use the correspondence to give nonvanishing results for the symplectic homology of certain Weinstein $4$-manifolds. We show a similar correspondence for the related case of Legendrian links in $J^1(S^1)$, the solid torus.


title: The SFR-M* Relation and Empirical Star-Formation Histories from ZFOURGE at 0.5 < z < 4

abstract: We explore star-formation histories (SFHs) of galaxies based on the evolution of the star-formation rate stellar mass relation (SFR-M*). Using data from the FourStar Galaxy Evolution Survey (ZFOURGE) in combination with far-IR imaging from the Spitzer and Herschel observatories we measure the SFR-M* relation at 0.5 < z < 4. Similar to recent works we find that the average infrared SEDs of galaxies are roughly consistent with a single infrared template across a broad range of redshifts and stellar masses, with evidence for only weak deviations. We find that the SFR-M* relation is not consistent with a single power-law of the form SFR ~ M*^a at any redshift; it has a power-law slope of a~1 at low masses, and becomes shallower above a turnover mass (M_0) that ranges from 10^9.5 - 10^10.8 Msol, with evidence that M_0 increases with redshift. We compare our measurements to results from state-of-the-art cosmological simulations, and find general agreement in the slope of the SFR-M* relation albeit with systematic offsets. We use the evolving SFR-M* sequence to generate SFHs, finding that typical SFRs of individual galaxies rise at early times and decline after reaching a peak. This peak occurs earlier for more massive galaxies. We integrate these SFHs to generate mass-growth histories and compare to the implied mass-growth from the evolution of the stellar mass function. We find that these two estimates are in broad qualitative agreement, but that there is room for improvement at a more detailed level. At early times the SFHs suggest mass-growth rates that are as much as 10x higher than inferred from the stellar mass function. However, at later times the SFHs under-predict the inferred evolution, as is expected in the case of additional growth due to mergers.


title: Input Sparsity and Hardness for Robust Subspace Approximation

abstract: In the subspace approximation problem, we seek a k-dimensional subspace F of R^d that minimizes the sum of p-th powers of Euclidean distances to a given set of n points a_1, ..., a_n in R^d, for p >= 1. More generally than minimizing sum_i dist(a_i,F)^p,we may wish to minimize sum_i M(dist(a_i,F)) for some loss function M(), for example, M-Estimators, which include the Huber and Tukey loss functions. Such subspaces provide alternatives to the singular value decomposition (SVD), which is the p=2 case, finding such an F that minimizes the sum of squares of distances. For p in [1,2), and for typical M-Estimators, the minimizing $F$ gives a solution that is more robust to outliers than that provided by the SVD. We give several algorithmic and hardness results for these robust subspace approximation problems. We think of the n points as forming an n x d matrix A, and letting nnz(A) denote the number of non-zero entries of A. Our results hold for p in [1,2). We use poly(n) to denote n^{O(1)} as n -> infty. We obtain: (1) For minimizing sum_i dist(a_i,F)^p, we give an algorithm running in O(nnz(A) + (n+d)poly(k/eps) + exp(poly(k/eps))), (2) we show that the problem of minimizing sum_i dist(a_i, F)^p is NP-hard, even to output a (1+1/poly(d))-approximation, answering a question of Kannan and Vempala, and complementing prior results which held for p >2, (3) For loss functions for a wide class of M-Estimators, we give a problem-size reduction: for a parameter K=(log n)^{O(log k)}, our reduction takes O(nnz(A) log n + (n+d) poly(K/eps)) time to reduce the problem to a constrained version involving matrices whose dimensions are poly(K eps^{-1} log n). We also give bicriteria solutions, (4) Our techniques lead to the first O(nnz(A) + poly(d/eps)) time algorithms for (1+eps)-approximate regression for a wide class of convex M-Estimators.


title: Bounded Control for Double Integrator in Quadrotor Dynamics

abstract: We construct a trajectory tracking controller for a quadrotor system by finding a coordinate change which transforms the quadrotor's vector field into that of a thrust propelled system. In a thrust propelled system, the goal is to stabilize its position around the origin, while the system is actuated by a one dimensional acceleration/thrust along a direction vector, by a time-varying gravity, and by the angular acceleration of the direction vector. For this system, a solution has been proposed in a companion article, submitted to ECC 2016, based on the implicit knowledge of a bounded controller for a double integrator system, and on the implicit knowledge of a Lyapunov function that guarantees the origin is asymptotically stable for the double integrator controlled by the bounded controller. We present two alternative bounded controllers for a double integrator system, and corresponding Lyapunov functions.


title: The cup subalgebra has the absorbing amenability property

abstract: Consider an inclusion of diffuse von Neumann algebras A c M . We say that A c M has the absorbing amenability property if for any diffuse subalgebra B c A and any amenable intermediate algebra B c D c M we have that D is contained in A. We prove that the cup subalgebra associated to any subfactor planar algebra has the absorbing amenability property.


title: Chebyshev approximation for multivariate functions

abstract: In this paper, we derive optimality conditions (Chebyshev approximation) for multivariate functions. The theory of Chebyshev (uniform) approximation for univariate functions is very elegant. The optimality conditions are based on the notion of alternance (maximal deviation points with alternating deviation signs). It is not very straightforward, however, how to extend the notion of alternance to the case of multivariate functions. There have been several attempts to extend the theory of Chebyshev approximation to the case of multivariate functions. We propose an alternative approach, which is based on the notion of convexity and nonsmooth analysis.


title: Coagulation-fragmentation model for animal group-size statistics

abstract: We study coagulation-fragmentation equations inspired by a simple model proposed in fisheries science to explain data for the size distribution of schools of pelagic fish. Although the equations lack detailed balance and admit no $H$-theorem, we are able to develop a rather complete description of equilibrium profiles and large-time behavior, based on recent developments in complex function theory for Bernstein and Pick functions. In the large-population continuum limit, a scaling-invariant regime is reached in which all equilibria are determined by a single scaling profile. This universal profile exhibits power-law behavior crossing over from exponent $-\frac23$ for small size to $-\frac32$ for large size, with an exponential cut-off.