title: Penning-trap mass spectrometry of highly charged, neutron-rich Rb and Sr isotopes in the vicinity of $A\approx100$

abstract: The neutron-rich mass region around $A\approx100$ presents challenges for modeling the astrophysical $r$-process because of rapid shape transitions. We report on mass measurements using the TITAN Penning trap at TRIUMF-ISAC to attain more reliable theoretical predictions of $r$-process nucleosynthesis paths in this region. A new approach using highly charged ($q=15+$) ions has been applied which considerably saves measurement time and preserves accuracy. New mass measurements of neutron-rich $^{94,97,98}$Rb and $^{94,97-99}$Sr have uncertainties of less than 4 keV and show deviations of up to 11$\sigma$ to previous measurements. An analysis using a parameterized $r$-process model is performed and shows that mass uncertainties for the A=90 abundance region are eliminated.


title: Beyond Random Walk and Metropolis-Hastings Samplers: Why You Should Not Backtrack for Unbiased Graph Sampling

abstract: Graph sampling via crawling has been actively considered as a generic and important tool for collecting uniform node samples so as to consistently estimate and uncover various characteristics of complex networks. The so-called simple random walk with re-weighting (SRW-rw) and Metropolis-Hastings (MH) algorithm have been popular in the literature for such unbiased graph sampling. However, an unavoidable downside of their core random walks -- slow diffusion over the space, can cause poor estimation accuracy. In this paper, we propose non-backtracking random walk with re-weighting (NBRW-rw) and MH algorithm with delayed acceptance (MHDA) which are theoretically guaranteed to achieve, at almost no additional cost, not only unbiased graph sampling but also higher efficiency (smaller asymptotic variance of the resulting unbiased estimators) than the SRW-rw and the MH algorithm, respectively. In particular, a remarkable feature of the MHDA is its applicability for any non-uniform node sampling like the MH algorithm, but ensuring better sampling efficiency than the MH algorithm. We also provide simulation results to confirm our theoretical findings.

categories: stat.ME cs.DS cs.NI cs.SI physics.data-an physics.soc-ph


title: Analysis of a Natural Gradient Algorithm on Monotonic Convex-Quadratic-Composite Functions

abstract: In this paper we investigate the convergence properties of a variant of the Covariance Matrix Adaptation Evolution Strategy (CMA-ES). Our study is based on the recent theoretical foundation that the pure rank-mu update CMA-ES performs the natural gradient descent on the parameter space of Gaussian distributions. We derive a novel variant of the natural gradient method where the parameters of the Gaussian distribution are updated along the natural gradient to improve a newly defined function on the parameter space. We study this algorithm on composites of a monotone function with a convex quadratic function. We prove that our algorithm adapts the covariance matrix so that it becomes proportional to the inverse of the Hessian of the original objective function. We also show the speed of covariance matrix adaptation and the speed of convergence of the parameters. We introduce a stochastic algorithm that approximates the natural gradient with finite samples and present some simulated results to evaluate how precisely the stochastic algorithm approximates the deterministic, ideal one under finite samples and to see how similarly our algorithm and the CMA-ES perform.


title: Comments on "Translation-invariant bipolarons and the problem of high-temperature superconductivity"

abstract: We comment on the recent results of Refs. [1, 2] on the bipolaron problem derived using an approximation of Gross - Tulub. It is proved that, contrary to the claim made in Refs. [1, 2], the bipolaron ground state energy calculated there in the strong-coupling approximation has not been shown to constitute a variational upper bound.


title: Qualitative properties of certain piecewise deterministic Markov processes

abstract: We study a class of Piecewise Deterministic Markov Processes with state space Rd x E where E is a finite set. The continuous component evolves according to a smooth vector field that is switched at the jump times of the discrete coordinate. The jump rates may depend on the whole position of the process. Working under the general assumption that the process stays in a compact set, we detail a possible construction of the process and characterize its support, in terms of the solutions set of a differential inclusion. We establish results on the long time behaviour of the process, in relation to a certain set of accessible points, which is shown to be strongly linked to the support of invariant measures. Under H\"ormander-type bracket conditions, we prove that there exists a unique invariant measure and that the processes converges to equilibrium in total variation. Finally we give examples where the bracket condition does not hold, and where there may be one or many invariant measures, depending on the jump rates between the flows.


title: A new variational formulation based on discontinuous Galerkin technique for a reaction-diffusion problem

abstract: In this paper, a new variational formulation based on discontinuous Galerkin technique for a reaction-diffusion problem is introduced, and the discontinuous Galerkin technique of this work is different from the general discontinuous Galerkin methods. The well posedness of the new formulation is given. Finally, it is pointed that the new variational formulation will be helpful to design better hybrid numerical methods which will not only strongly stable in spatial variable and absolutely stable in temporal variable but also be optimally convergent.


title: Learning From An Optimization Viewpoint

abstract: In this dissertation we study statistical and online learning problems from an optimization viewpoint.The dissertation is divided into two parts : I. We first consider the question of learnability for statistical learning problems in the general learning setting. The question of learnability is well studied and fully characterized for binary classification and for real valued supervised learning problems using the theory of uniform convergence. However we show that for the general learning setting uniform convergence theory fails to characterize learnability. To fill this void we use stability of learning algorithms to fully characterize statistical learnability in the general setting. Next we consider the problem of online learning. Unlike the statistical learning framework there is a dearth of generic tools that can be used to establish learnability and rates for online learning problems in general. We provide online analogs to classical tools from statistical learning theory like Rademacher complexity, covering numbers, etc. We further use these tools to fully characterize learnability for online supervised learning problems. II. In the second part, for general classes of convex learning problems, we provide appropriate mirror descent (MD) updates for online and statistical learning of these problems. Further, we show that the the MD is near optimal for online convex learning and for most cases, is also near optimal for statistical convex learning. We next consider the problem of convex optimization and show that oracle complexity can be lower bounded by the so called fat-shattering dimension of the associated linear class. Thus we establish a strong connection between offline convex optimization problems and statistical learning problems. We also show that for a large class of high dimensional optimization problems, MD is in fact near optimal even for convex optimization.


title: Restoration of Rotational Symmetry in the Continuum Limit of Lattice Field Theories

abstract: We explore how rotational invariance is systematically recovered from calculations on hyper-cubic lattices through the use of smeared lattice operators that smoothly evolve into continuum operators with definite angular momentum as the lattice-spacing is reduced. Perturbative calculations of the angular momentum violation associated with such operators at tree-level and at one-loop are presented in phi^4 theory and QCD. Contributions from these operators that violate rotational invariance occur at tree-level, with coefficients that are suppressed by O(a^2) in the continuum limit. Quantum loops do not modify this behavior in phi^4, nor in QCD if the gauge-fields are smeared over a comparable spatial region. Consequently, the use of this type of operator should, in principle, allow for Lattice QCD calculations of the higher moments of the hadron structure functions.


title: Controlling trapping potentials and stray electric fields in a microfabricated ion trap through design and compensation

abstract: Recent advances in quantum information processing with trapped ions have demonstrated the need for new ion trap architectures capable of holding and manipulating chains of many (>10) ions. Here we present the design and detailed characterization of a new linear trap, microfabricated with scalable complementary metal-oxide-semiconductor (CMOS) techniques, that is well-suited to this challenge. Forty-four individually controlled DC electrodes provide the many degrees of freedom required to construct anharmonic potential wells, shuttle ions, merge and split ion chains, precisely tune secular mode frequencies, and adjust the orientation of trap axes. Microfabricated capacitors on DC electrodes suppress radio-frequency pickup and excess micromotion, while a top-level ground layer simplifies modeling of electric fields and protects trap structures underneath. A localized aperture in the substrate provides access to the trapping region from an oven below, permitting deterministic loading of particular isotopic/elemental sequences via species-selective photoionization. The shapes of the aperture and radio-frequency electrodes are optimized to minimize perturbation of the trapping pseudopotential. Laboratory experiments verify simulated potentials and characterize trapping lifetimes, stray electric fields, and ion heating rates, while measurement and cancellation of spatially-varying stray electric fields permits the formation of nearly-equally spaced ion chains.


title: Algorithm for multivariate data standardization up to third moment

abstract: An algorithm for transforming multivariate data to a form with normalized first, second and third moments is presented.


title: Quantum phase transition in a multicomponent anyonic Lieb-Liniger model

abstract: We study a one-dimensional multicomponent anyon model that reduces to a multicomponent Lieb-Liniger gas of impenetrable bosons (Tonks-Girardeau gas) for vanishing statistics parameter. At fixed component densities, the coordinate Bethe ansatz gives a family of quantum phase transitions at special values of the statistics parameter. We show that the ground state energy changes extensively between different phases. Special regimes are studied and a general classification for the transition points is given. An interpretation in terms of statistics of composite particles is proposed.


title: The parsec-scale jet of PKS 1749+096

abstract: PKS 1749+096 is a BL Lac object showing weak extended jet emission to the northeast of the compact VLBI core on parsec scales. We aim at better understanding the jet kinematics and variability of this source and finding clues that may applicable to other BL Lac objects. The jet was studied with multi-epoch multi-frequency high-resolution VLBI observations. The jet is characterized by a one-sided curved morphology at all epochs and all frequencies. The VLBI core, located at the southern end of the jet, was identified based on its spectral properties. The equipartition magnetic field of the core was investigated, through which we derived a Doppler factor of 5, largely consistent with that derived from kinematics (component C5). The study of the detailed jet kinematics at 22 and 15 GHz, spanning a period of more than 10 years, indicates the possible existence of a bimodal distribution of the jet apparent speed. Ballistic and non-ballistic components are found to coexist in the jet. Superluminal motions in the range of 5-21 c were measured in 11 distinct components. We estimated the physical jet parameters with the minimum Lorentz factor of 10.2 and Doppler factors in the range of 10.2-20.4 (component C5). The coincidence in time of the component's ejection and flares supports the idea that, at least in PKS 1749+096, ejection of new jet components is connected with major outbursts in flux density. For the best-traced component (C5) we found that the flux density decays rapidly as it travels downstream the jet, accompanied by a steepening of its spectra, which argues in favor of a contribution of inverse Compton cooling. These properties make PKS 1749+096 a suitable target for an intensive monitoring to decipher the variability phenomenon of BL Lac objects.


title: Ultra Low Complexity Soft Output Detector for Non-Binary LDPC Coded Large MIMO Systems

abstract: The theoretic results of MIMO capacity tell us that the higher the number of antennas are employed, the higher the transmission rate is. This makes MIMO systems with hundreds of antennas very attractive but one of the major problems that obstructs such large dimensional MIMO systems from the practical realization is a high complexity of the MIMO detector. We present in this paper the new soft output MIMO detector based on matched filtering that can be applied to the large MIMO systems which are coded by the powerful non-binary LDPC codes. The per-bit complexity of the proposed detector is just 0.28% to that of low complexity soft output MMSE detector and scales only linearly with a number of antennas. Furthermore, the coded performances with small information length 800 bits are within 4.2 dB from the associated MIMO capacity.


title: Numerical evaluation of multi-loop integrals for arbitrary kinematics with SecDec 2.0

abstract: We present the program SecDec 2.0 which contains various new features: First, it allows the numerical evaluation of multi-loop integrals with no restriction on the kinematics. Dimensionally regulated ultraviolet and infrared singularities are isolated via sector decomposition, while threshold singularities are handled by a deformation of the integration contour in the complex plane. As an application we present numerical results for various massive two-loop four-point diagrams. SecDec 2.0 also contains new useful features for the calculation of more general parameter integrals, related e.g. to phase space integrals.


title: Multilevel Sparse Kernel-Based Interpolation

abstract: A multilevel kernel-based interpolation method, suitable for moderately high-dimensional function interpolation problems, is proposed. The method, termed multilevel sparse kernel-based interpolation (MLSKI, for short), uses both level-wise and direction-wise multilevel decomposition of structured (or mildly unstructured) interpolation data sites in conjunction with the application of kernel-based interpolants with different scaling in each direction. The multilevel interpolation algorithm is based on a hierarchical decomposition of the data sites, whereby at each level the detail is added to the interpolant by interpolating the resulting residual of the previous level. On each level, anisotropic radial basis functions are used for solving a number of small interpolation problems, which are subsequently linearly combined to produce the interpolant. MLSKI can be viewed as an extension of $d$-boolean interpolation (which is closely related to ideas in sparse grid and hyperbolic crosses literature) to kernel-based functions, within the hierarchical multilevel framework to achieve accelerated convergence. Numerical experiments suggest that the new algorithm is numerically stable and efficient for the reconstruction of large data in $\mathbb{R}^{d}\times \mathbb{R}$, for $d = 2, 3, 4$, with tens or even hundreds of thousands data points. Also, MLSKI appears to be generally superior over classical radial basis function methods in terms of complexity, run time and convergence at least for large data sets.


title: The Artificial Regression Market

abstract: The Artificial Prediction Market is a recent machine learning technique for multi-class classification, inspired from the financial markets. It involves a number of trained market participants that bet on the possible outcomes and are rewarded if they predict correctly. This paper generalizes the scope of the Artificial Prediction Markets to regression, where there are uncountably many possible outcomes and the error is usually the MSE. For that, we introduce the reward kernel that rewards each participant based on its prediction error and we derive the price equations. Using two reward kernels we obtain two different learning rules, one of which is approximated using Hermite-Gauss quadrature. The market setting makes it easy to aggregate specialized regressors that only predict when an observation falls into their specialization domain. Experiments show that regression markets based on the two learning rules outperform Random Forest Regression on many UCI datasets and are rarely outperformed.


title: Inflation from the Higgs field false vacuum with hybrid potential

abstract: We have recently suggested [1,2] that Inflation could have started in a local minimum of the Higgs potential at field values of about $10^{15}-10^{17}$ GeV, which exists for a narrow band of values of the top quark and Higgs masses and thus gives rise to a prediction on the Higgs mass to be in the range 123-129 GeV, together with a prediction on the the top mass and the cosmological tensor-to-scalar ratio $r$. Inflation can be achieved provided there is an additional degree of freedom which allows the transition to a radiation era. In [1] we had proposed such field to be a Brans-Dicke scalar. Here we present an alternative possibility with an additional subdominant scalar very weakly coupled to the Higgs, realizing an (inverted) hybrid Inflation scenario. Interestingly, we show that such model has an additional constraint $m_H<125.3 \pm 3_{th}$, where $3_{th}$ is the present theoretical uncertainty on the Standard Model RGEs. The tensor-to-scalar ratio has to be within the narrow range $10^{-4}\lesssim r<0.007$, and values of the scalar spectral index compatible with the observed range can be obtained. Moreover, if we impose the model to have subplanckian field excursion, this selects a narrower range $10^{-4} \lesssim r<0.001$ and an upper bound on the Higgs mass of about $m_H <124 \pm 3_{th}$.


title: Collective excitations and marginal stability of quantum Ising spin glasses

abstract: We solve the Sherrington-Kirkpatrick (SK) model in a transverse field \Gamma deep in its quantum glass phase at zero temperature. We show that the glass phase is critical everywhere, exhibiting collective excitations with a gapless Ohmic spectral function. Using an effective potential approach, we interpret the latter as being due to disordered spin waves which behave as weakly coupled, underdamped harmonic oscillators. In the limit of small transverse field the low frequency tail of the spectrum tends to a universal, \Gamma-independent limit.


title: Supersymmetric BCS

abstract: We implement relativistic BCS superconductivity in N=1 supersymmetric field theories with a U(1)_R symmetry. The simplest model contains two chiral superfields with a Kahler potential modified by quartic terms. We study the phase diagram of the gap as a function of the temperature and the specific heat. The superconducting phase transition turns out to be first order, due to the scalar contribution to the one-loop potential. By virtue of supersymmetry, the critical curves depend logarithmically with the UV cutoff, rather than quadratically as in standard BCS theory. We comment on the difficulties in having fermion condensates when the chemical potential is instead coupled to a baryonic U(1)_B current. We also discuss supersymmetric models of BCS with canonical Kahler potential constructed by "integrating-in" chiral superfields.


title: Small generators of function fields

abstract: Let $K/k$ be a finite extension of a global field. Such an extension can be generated over $k$ by a single element. The aim of this article is to prove the existence of a "small" generator in the function field case. This answers the function field version of a question of Ruppert on small generators of number fields.


title: Connectivity Oracles for Planar Graphs

abstract: We consider dynamic subgraph connectivity problems for planar graphs. In this model there is a fixed underlying planar graph, where each edge and vertex is either "off" (failed) or "on" (recovered). We wish to answer connectivity queries with respect to the "on" subgraph. The model has two natural variants, one in which there are $d$ edge/vertex failures that precede all connectivity queries, and one in which failures/recoveries and queries are intermixed. We present a $d$-failure connectivity oracle for planar graphs that processes any $d$ edge/vertex failures in $sort(d,n)$ time so that connectivity queries can be answered in $pred(d,n)$ time. (Here $sort$ and $pred$ are the time for integer sorting and integer predecessor search over a subset of $[n]$ of size $d$.) Our algorithm has two discrete parts. The first is an algorithm tailored to triconnected planar graphs. It makes use of Barnette's theorem, which states that every triconnected planar graph contains a degree-3 spanning tree. The second part is a generic reduction from general (planar) graphs to triconnected (planar) graphs. Our algorithm is, moreover, provably optimal. An implication of Patrascu and Thorup's lower bound on predecessor search is that no $d$-failure connectivity oracle (even on trees) can beat $pred(d,n)$ query time. We extend our algorithms to the subgraph connectivity model where edge/vertex failures (but no recoveries) are intermixed with connectivity queries. In triconnected planar graphs each failure and query is handled in $O(\log n)$ time (amortized), whereas in general planar graphs both bounds become $O(\log^2 n)$.


title: On Shape Dependence and RG Flow of Entanglement Entropy

abstract: We use a mix of field theoretic and holographic techniques to elucidate various properties of quantum entanglement entropy. In (3+1)-dimensional conformal field theory we study the divergent terms in the entropy when the entangling surface has a conical or a wedge singularity. In (2+1)-dimensional field theory with a mass gap we calculate, for an arbitrary smooth entanglement contour, the expansion of the entropy in inverse odd powers of the mass. We show that the shape-dependent coefficients that arise are even powers of the extrinsic curvature and its derivatives. A useful dual construction of a (2+1)-dimensional theory, which allows us to exhibit these properties, is provided by the CGLP background. This smooth warped throat solution of 11-dimensional supergravity describes renormalization group flow from a conformal field theory in the UV to a gapped one in the IR. For this flow we calculate the recently introduced renormalized entanglement entropy and confirm that it is a monotonic function.


title: Azimuthal distributions of charged hadrons, pions, and kaons produced in deep-inelastic scattering off unpolarized protons and deuterons

abstract: The azimuthal cos{\phi} and cos2{\phi} modulations of the distribution of hadrons produced in unpolarized semi-inclusive deep-inelastic scattering of electrons and positrons off hydrogen and deuterium targets have been measured in the HERMES experiment. For the first time these modulations were determined in a four-dimensional kinematic space for positively and negatively charged pions and kaons separately, as well as for unidentified hadrons. These azimuthal dependences are sensitive to the transverse motion and polarization of the quarks within the nucleon via, e.g., the Cahn, Boer-Mulders and Collins effects.


title: Evaluations for zeta(2),zeta(4),...,zeta(2k)based on the WZ method

abstract: Based on the framework of the WZ theory, a new evaluation for $\varsigma (2) = \frac{\pi ^2}{6}$ and $\varsigma (4) = \frac{\pi ^4}{90}$ was given respectively, finally, a new recurrence formula for $\varsigma (2k)$ was given.


title: A Self-Consistent Model for Positronium Formation from Helium Atoms

abstract: The differential and total cross sections for electron capture by positrons from helium atoms are calculated using a first-order distorted wave theory satisfying the Coulomb boundary conditions. In this formalism a parametric potential is used to describe the electron screening in a consistent and realistic manner. The present procedure is self consistent because (i) it satisfies the correct boundary conditions and post-prior symmetry, and (ii) the potential and the electron binding energies appearing in the transition amplitude are consistent with the wave functions describing the collision system. The results are compared with the other theories and with the available experimental measurements. At the considered range of collision energies, the results agree reasonably well with recent experiments and theories. [Note: This paper will be published on volume 42 of the Brazilian Journal of Physics]


title: Confronting Seiberg's Duality with r Duality in N=1 Supersymmetric QCD

abstract: Systematizing our results on r duality obtained previously we focus on comparing r duality with the generalized Seiberg duality in the r vacua of N=2 and N=1 super-Yang-Mills theories with the U(N) gauge group and N_f matter flavors (N_f>N). The number of condensed (s)quarks r is assumed to be in the interval 2N_f/3 < r \le N. To pass to N=1 we introduce an N=2-breaking deformation, a mass term \mu for the adjoint matter, eventually decoupling the adjoint matter in the limit of large \mu. If one starts from a large value of the parameter \xi\sim\mu m, where the original theory is at weak coupling, and decreases \xi one hits a a crossover transition from weak to strong coupling (here m is a typical value of the quark masses). Below this transition the original theory is described in terms of a weakly coupled infrared-free r dual theory with the U(N_f-r) gauge group and N_f light quark-like dyon flavors. Dyon condensation leads to confinement of monopoles, defying a naive expectation of quark confinement. The quarks and gauge bosons of the original theory are in an "instead-of-confinement" phase. The r and Seiberg dualities are demonstrated to coincide in the r=N vacua. In the (2/3) N_f <r<N vacua two dualities do not match. In this window Seiberg's dual is at strong coupling while our r-dual model is at weak coupling. Thus, we can speak of triality. Seiberg's dual solution at weak coupling reappears again at r<N_f-N< (1/3) N_f$.


title: r Duality and "Instead-of-Confinement" Mechanism in N=1 Supersymmetric QCD

abstract: We consider N=2 SQCD with the U(N) gauge group and N_f flavors (N_f>N) perturbed by an N=2 breaking deformation - a small mass term \mu for the adjoint matter. We study r-vacua, with the constraint $2N_f/3 < r \le N. At large values of the parameter \xi\sim\mu m (m is a typical value of the quark masses) r quark flavors condense, by construction. The effective low-energy theory with the gauge group U(r)\times U(1)^{N-r} is at weak coupling. Upon reducing \xi the original theory undergoes a crossover transition from weak to strong coupling. As the original theory becomes strongly coupled, at low energies it is described by a weakly coupled infrared-free dual theory with the gauge group U(N_f-r)\times U(1)^{N-N_f+r} and N_f light dyon flavors. These dyons condense triggering formation of non-Abelian strings which still confine monopoles, rather than quarks, contrary to naive duality arguments. "Instead-of-confinement" mechanism for quarks and gauge bosons of the original theory takes place: screened quarks and gauge bosons of the original theory decay, on curves of the marginal stability (CMS), into confined monopole-antimonopole pairs that form stringy mesons. Next, we increase the deformation parameter \mu thus decoupling the adjoint fields. Then our theory flows to N=1 SQCD. The gauge group of the dual theory becomes U(N_f-r). We show that the dual theory is weakly coupled if we are sufficiently close to the Argyres-Douglas point. The "instead-of-confinement" mechanism for quarks and gauge bosons survives in the limit of large \mu. It determines low-energy non-Abelian dynamics in the r-vacua of N=1 SQCD.


title: Message passing with relaxed moment matching

abstract: Bayesian learning is often hampered by large computational expense. As a powerful generalization of popular belief propagation, expectation propagation (EP) efficiently approximates the exact Bayesian computation. Nevertheless, EP can be sensitive to outliers and suffer from divergence for difficult cases. To address this issue, we propose a new approximate inference approach, relaxed expectation propagation (REP). It relaxes the moment matching requirement of expectation propagation by adding a relaxation factor into the KL minimization. We penalize this relaxation with a $l_1$ penalty. As a result, when two distributions in the relaxed KL divergence are similar, the relaxation factor will be penalized to zero and, therefore, we obtain the original moment matching; In the presence of outliers, these two distributions are significantly different and the relaxation factor will be used to reduce the contribution of the outlier. Based on this penalized KL minimization, REP is robust to outliers and can greatly improve the posterior approximation quality over EP. To examine the effectiveness of REP, we apply it to Gaussian process classification, a task known to be suitable to EP. Our classification results on synthetic and UCI benchmark datasets demonstrate significant improvement of REP over EP and Power EP--in terms of algorithmic stability, estimation accuracy and predictive performance.


title: The stability to instability transition in the structure of large scale networks

abstract: We examine phase transitions between the easy, hard, and the unsolvable phases when attempting to identify structure in large complex networks (community detection) in the presence of disorder induced by network noise (spurious links that obscure structure), heat bath temperature $T$, and system size $N$. When present, transitions at low temperature or low noise correspond to entropy driven (or "order by disorder") annealing effects wherein stability may initially increase as temperature or noise is increased before becoming unsolvable at sufficiently high temperature or noise. Additional transitions between contending viable solutions (such as those at different natural scales) are also possible. When analyzing community structure via a dynamical approach, "chaotic-type" transitions were previously identified [Phil. Mag. {\bf 92} 406 (2012)]. The correspondence between the spin-glass-type complexity transitions and transitions into chaos in dynamical analogs might extend to other hard computational problems. In this work, we examine large networks that have a large number of communities. We infer that large systems at a constant ratio of $q$ to the number of nodes $N$ asymptotically tend toward insolvability in the limit of large $N$ for any positive $T$. The asymptotic behavior of temperatures below which structure identification might be possible, $T_\times =O[1/\log q]$, decreases slowly, so for practical system sizes, there remains an accessible, and generally easy, global solvable phase at low temperature. We further employ multivariate Tutte polynomials to show that increasing $q$ emulates increasing $T$ for a general Potts model, leading to a similar stability region at low $T$. Given the relation between Tutte and Jones polynomials, our results further suggest a link between the above complexity transitions and transitions associated with random knots.


title: D-modules with finite support are semi-simple

abstract: Let $(R, \mf, k_R)$ be regular local $k$-algebra satisfying the weak Jacobian criterion, such that $k_R/k$ is an algebraic field extension. Let $D_R$ be the ring of $k$-linear differential operators of $R$. We give an explicit decomposition of the $D_R$-module $D_R/D_R \mf_R^{n+1}$ as a direct sum of simple modules, all isomorphic to $D_R/D_R \mf$, where certain "Pochhammer" differential operators are used to describe generators of the simple components.


title: Apparent Fracture in Polymeric Fluids under Step Shear

abstract: Recent step strain experiments in well-entangled polymeric liquids demonstrated a bulk fracture-like phenomenon. We have studied this instability using a modern version of the Doi-Edwards theory for entangled polymers, and we find close quantitative agreement with the experiments. The phenomenon occurs because the viscoelastic liquid is sheared into a rubbery state that possesses an elastic constitutive instability (Marrucci and Grizzuti, 1983). The fracture is a transient manifestation of this instability, which relies on the amplification of spatially inhomogeneous fluctuations. This mechanism differs from fracture in glassy materials and dense suspensions.


title: Quantum Field Lagrangian model for charge density waves in one-dimensional systems at finite temperature : A Thermofield Dynamics Approach

abstract: We consider the Thermofield Dynamics bosonization to perform a field theory analysis of the effective Lagrangian model for incommensurate charge density waves (ICDW) in one-dimensional systems at finite temperature. The phonon degree of freedom is carryied by a dynamical phase field, contributing to the quantum dynamics and symmetry related features of the ICDW phenomenon. The electron chiral density and the phase of the phonon field condensate as a thermal soliton, carrying the symmetry under the linked electron-phonon $U_e^5(1) \otimes U_{ph}(1)$ global transformations. Using the Gell'Mann-Low formula for finite temperature, the perturbative series of the phonon thermal correlation function is obtained. Due to the electron-phonon charge selection rule we obtain for the thermal vacuum expectation value for the order parameter $<0 (\beta) |\Phi| 0 (\beta)> = 0$, in accordance with the cluster decomposition property of the corresponding correlation function. This reflects the fact that the quantum description of the ICDW corresponds to a local charge transport through the lattice which is accomplished by an electron-lattice energy redistribution, which accounts for a thermal dynamical mass gap generation. The electron-phonon coupling can be rewritten in terms of a mass operator for the "physical" fermion operator $\Psi$ such that $<0 (\beta)| \bar\Psi (x)\Psi (x) | 0 (\beta)> \neq 0$, without the breakdown of the linked electron-phonon symmetry.


title: Time-Reversal-Invariant Hofstadter-Hubbard Model with Ultracold Fermions

abstract: We consider the time-reversal-invariant Hofstadter-Hubbard model which can be realized in cold atom experiments. In these experiments, an additional staggered potential and an artificial Rashba--type spin-orbit coupling are available. Without interactions, the system exhibits various phases such as topological and normal insulator, metal as well as semi--metal phases with two or even more Dirac cones. Using a combination of real-space dynamical mean-field theory and analytical techniques, we discuss the effect of on-site interactions and determine the corresponding phase diagram. In particular, we investigate the semi--metal to antiferromagnetic insulator transition and the stability of different topological insulator phases in the presence of strong interactions. We compute spectral functions which allow us to study the edge states of the strongly correlated topological phases.


title: Extremal vs. Non-Extremal Correlators with Giant Gravitons

abstract: We consider extremal and non-extremal three-point functions of two giant gravitons and one point-like graviton using Schur polynomials in N=4 super Yang-Mills theory and holographically, using a semiclassical Born-Infeld analysis as well as bubbling geometries. For non-extremal three-point functions our computations using all three approaches are in perfect agreement. For extremal correlators we find that our results from the bubbling geometry analysis agree with existing results from the gauge theory. The semiclassical Born-Infeld computation for the extremal case is known to give a different answer, which we interpret as a manifestation of the known subtlety of holography for extremal correlators.


title: Rudin-Kisler ordering on the P-hierarchy

abstract: M. E. Rudin in "Partial orders on the types of $\beta {N}$", (Trans. Amer. Math. Soc., 155, 1971, 353-362) proved (under CH) that for each P-point $u$ there is a P-point $v$ such that $v>_{RK}u$. A. Blass in "Rudin - Kisler ordering on P-points" (Trans. Amer. Math. Soc. 179, 1973, 145-166) improved that theorem assuming MA in the place of CH, in that paper he also proved that under MA each RK-increasing sequence of P-points is upper bounded by a P-point. We improve Blass results simultaneously in 3 directions - we prove it for each class of index $\geq 2$ of P-hierarchy (P-points coincidence with a class $P_2$ of P-hierarchy), assuming $\mathfrak{b}=\mathfrak{c}$ in the place of MA and we show that there are at least $\mathfrak{b}$ many Rudin-Kisler incomparable such upper bounds.


title: The inexorable resistance of inertia determines the initial regime of drop coalescence

abstract: Drop coalescence is central to diverse processes involving dispersions of drops in industrial, engineering and scientific realms. During coalescence, two drops first touch and then merge as the liquid neck connecting them grows from initially microscopic scales to a size comparable to the drop diameters. The curvature of the interface is infinite at the point where the drops first make contact, and the flows that ensue as the two drops coalesce are intimately coupled to this singularity in the dynamics. Conventionally, this process has been thought to have just two dynamical regimes: a viscous and an inertial regime with a crossover region between them. We use experiments and simulations to reveal that a third regime, one that describes the initial dynamics of coalescence for all drop viscosities, has been missed. An argument based on force balance allows the construction of a new coalescence phase diagram.


title: Ab-initio non-equilibrium quantum transport and forces with the real space projector augmented wave method

abstract: We present an efficient implemention of a non-equilibrium Green function (NEGF) method for self-consistent calculations of electron transport and forces in nanostructured materials. The electronic structure is described at the level of density functional theory (DFT) using the projector augmented wave method (PAW) to describe the ionic cores and an atomic orbital basis set for the valence electrons. External bias and gate voltages are treated in a self-consistent manner and the Poisson equation with appropriate boundary conditions is solved in real space. Contour integration of the Green function and parallelization over k-points and real space makes the code highly efficient and applicable to systems containing several hundreds of atoms. The method is applied to a number of different systems demonstrating the effects of bias and gate voltages, multiterminal setups, non-equilibrium forces, and spin transport.


title: Deterministic Function Computation with Chemical Reaction Networks

abstract: Chemical reaction networks (CRNs) formally model chemistry in a well-mixed solution. CRNs are widely used to describe information processing occurring in natural cellular regulatory networks, and with upcoming advances in synthetic biology, CRNs are a promising language for the design of artificial molecular control circuitry. Nonetheless, despite the widespread use of CRNs in the natural sciences, the range of computational behaviors exhibited by CRNs is not well understood. CRNs have been shown to be efficiently Turing-universal when allowing for a small probability of error. CRNs that are guaranteed to converge on a correct answer, on the other hand, have been shown to decide only the semilinear predicates. We introduce the notion of function, rather than predicate, computation by representing the output of a function f:N^k --> N^l by a count of some molecular species, i.e., if the CRN starts with n_1,...,n_k molecules of some "input" species X1,...,Xk, the CRN is guaranteed to converge to having f(n_1,...,n_k) molecules of the "output" species Y1,...,Yl. We show that a function f:N^k --> N^l is deterministically computed by a CRN if and only if its graph {(x,y) | f(x) = y} is a semilinear set. Furthermore, each semilinear function f can be computed on input x in expected time O(polylog(|x|)).


title: Transconductance fluctuations as a probe for interaction induced quantum Hall states in graphene

abstract: Transport measurements normally provide a macroscopic, averaged view of the sample, so that disorder prevents the observation of fragile interaction induced states. Here, we demonstrate that transconductance fluctuations in a graphene field effect transistor reflect charge localization phenomena on the nanometer scale due to the formation of a dot network which forms near incompressible quantum states. These fluctuations give access to fragile broken-symmetry and fractional quantum Hall states even though these states remain hidden in conventional magnetotransport quantities.


title: Accessing surface Brillouin zone and band structure of picene single crystals

abstract: We have experimentally revealed the band structure and the surface Brillouin zone of insulating picene single crystals (SCs), the mother organic system for a recently discovered aromatic superconductor, with ultraviolet photoelectron spectroscopy (UPS) and low-energy electron diffraction with laser for photoconduction. A hole effective mass of 2.24 m_0 and the hole mobility mu_h >= 9.0 cm^2/Vs (298 K) were deduced in Gamma-Y direction. We have further shown that some picene SCs did not show charging during UPS even without the laser, which indicates that pristine UPS works for high-quality organic SCs.


title: Decoherence in the chemical compass: The role of decoherence for avian magnetoreception

abstract: Contrary to the usual picture that decoherence destroys quantum properties and causes the quantum-to-classical transition, we argue that decoherence can also play a constructive role in driving quantum dynamics and amplifying its results to macroscopic scales. We support this perspective by presenting an example system from spin chemistry, which is also of importance for biological systems, e.g. in avian magnetoreception.

categories: physics.bio-ph physics.chem-ph physics.hist-ph quant-ph


title: Construction of a Calibrated Probabilistic Classification Catalog: Application to 50k Variable Sources in the All-Sky Automated Survey

abstract: With growing data volumes from synoptic surveys, astronomers must become more abstracted from the discovery and introspection processes. Given the scarcity of follow-up resources, there is a particularly sharp onus on the frameworks that replace these human roles to provide accurate and well-calibrated probabilistic classification catalogs. Such catalogs inform the subsequent follow-up, allowing consumers to optimize the selection of specific sources for further study and permitting rigorous treatment of purities and efficiencies for population studies. Here, we describe a process to produce a probabilistic classification catalog of variability with machine learning from a multi-epoch photometric survey. In addition to producing accurate classifications, we show how to estimate calibrated class probabilities, and motivate the importance of probability calibration. We also introduce a methodology for feature-based anomaly detection, which allows discovery of objects in the survey that do not fit within the predefined class taxonomy. Finally, we apply these methods to sources observed by the All Sky Automated Survey (ASAS), and unveil the Machine-learned ASAS Classification Catalog (MACC), which is a 28-class probabilistic classification catalog of 50,124 ASAS sources. We estimate that MACC achieves a sub-20% classification error rate, and demonstrate that the class posterior probabilities are reasonably calibrated. MACC classifications compare favorably to the classifications of several previous domain-specific ASAS papers and to the ASAS Catalog of Variable Stars, which had classified only 24% of those sources into one of 12 science classes. The MACC is publicly available at http://www.bigmacc.info.


title: Discovery of Bright Galactic R Coronae Borealis and DY Persei Variables: Rare Gems Mined from ACVS

abstract: We present the results of a machine-learning (ML) based search for new R Coronae Borealis (RCB) stars and DY Persei-like stars (DYPers) in the Galaxy using cataloged light curves from the All-Sky Automated Survey (ASAS) Catalog of Variable Stars (ACVS). RCB stars - a rare class of hydrogen-deficient carbon-rich supergiants - are of great interest owing to the insights they can provide on the late stages of stellar evolution. DYPers are possibly the low-temperature, low-luminosity analogs to the RCB phenomenon, though additional examples are needed to fully establish this connection. While RCB stars and DYPers are traditionally identified by epochs of extreme dimming that occur without regularity, the ML search framework more fully captures the richness and diversity of their photometric behavior. We demonstrate that our ML method can use newly discovered RCB stars to identify additional candidates within the same data set. Our search yields 15 candidates that we consider likely RCB stars/DYPers: new spectroscopic observations confirm that four of these candidates are RCB stars and four are DYPers. Our discovery of four new DYPers increases the number of known Galactic DYPers from two to six; noteworthy is that one of the new DYPers has a measured parallax and is m ~ 7 mag, making it the brightest known DYPer to date. Future observations of these new DYPers should prove instrumental in establishing the RCB connection. We consider these results, derived from a machine-learned probabilistic classification catalog, as an important proof-of-concept for the efficient discovery of rare sources with time-domain surveys.


title: SSOS: A Moving Object Image Search Tool for Asteroid Precovery

abstract: It is very difficult to find archival images of solar system objects. While regular archive searches can find images at a fixed location, they cannot find images of moving targets. Archival images have become increasingly useful to galactic and stellar astronomers the last few years but, until now, solar system researchers have been at a disadvantage in this respect. The Solar System Object Search (SSOS) at the Canadian Astronomy Data Centre allows users to search for images of moving objects. SSOS accepts as input either a list of observations, an object designation, a set of orbital elements, or a user-generated ephemeris for an object. It then searches for images containing that object over a range of dates. The user is then presented with a list of images containing that object from a variety of archives. Initially created to search the CFHT MegaCam archive, SSOS has been extended to other telescope archives including Gemini, Subaru/SuprimeCam, HST, several ESO instruments and the SDSS for a total of 6.5 million images. The SSOS tool is located on the web at: http://www.cadc.hia.nrc.gc.ca/ssos


title: Distinguishing Dynamical Dark Matter at the LHC

abstract: Dynamical dark matter (DDM) is a new framework for dark-matter physics in which the dark sector comprises an ensemble of individual component fields which collectively conspire to act in ways that transcend those normally associated with dark matter. Because of its non-trivial structure, this DDM ensemble --- unlike most traditional dark-matter candidates --- cannot be characterized in terms of a single mass, decay width, or set of scattering cross-sections, but must instead be described by parameters which describe the collective behavior of its constituents. Likewise, the components of such an ensemble need not be stable so long as lifetimes are balanced against cosmological abundances across the ensemble as a whole. In this paper, we investigate the prospects for identifying a DDM ensemble at the LHC and for distinguishing such a dark-matter candidate from the candidates characteristic of traditional dark-matter models. In particular, we focus on DDM scenarios in which the component fields of the ensemble are produced at colliders alongside some number of Standard-Model particles via the decays of additional heavy fields. The invariant-mass distributions of these Standard-Model particles turn out to possess several unique features that cannot be replicated in most traditional dark-matter models. We demonstrate that in many situations it is possible to differentiate between a DDM ensemble and a traditional dark-matter candidate on the basis of such distributions. Moreover, many of our results also apply more generally to a variety of other extensions of the Standard Model which involve multiple stable or metastable neutral particles.


title: A fundamental problem in our understanding of low mass galaxy evolution

abstract: Recent studies have found a dramatic difference between the observed number density evolution of low mass galaxies and that predicted by semi-analytic models. While models accurately reproduce the z=0 number density, they require that the evolution occurs rapidly at early times, which is incompatible with the strong late evolution found observationally. We report here the same discrepancy in two state-of-the-art cosmological hydrodynamical simulations, which is evidence that the problem is fundamental. We search for the underlying cause of this problem using two complementary methods. Firstly, we look for evidence of a different history of today's low mass galaxies in models and observations and we find that the models yield too few young, strongly star-forming galaxies. Secondly, we construct a toy model to link the observed evolution of specific star formation rates (sSFR) with the evolution of the galaxy stellar mass function. We infer from this model that a key problem in both semi-analytic and hydrodynamical models is the presence of a positive instead of a negative correlation between sSFR and stellar mass. A similar positive correlation is found between the specific dark matter halo accretion rate and the halo mass, indicating that model galaxies are growing in a way that follows the growth of their host haloes too closely. It therefore appears necessary to find a mechanism that decouples the growth of low mass galaxies, which occurs primarily at late times, from the growth of their host haloes, which occurs primarily at early times. We argue that the current form of star-formation driven feedback implemented in most galaxy formation models is unlikely to achieve this goal, owing to its fundamental dependence on host halo mass and time. [Abridged]


title: Stable counteralignment of a circumbinary disc

abstract: In general, when gas accretes on to a supermassive black hole binary it is likely to have no prior knowledge of the binary angular momentum. Therefore a circumbinary disc forms with a random inclination angle, theta, to the binary. It is known that for theta < 90 degrees the disc will coalign wrt the binary. If theta > 90 degrees the disc wholly counteraligns if it satisfies cos(theta) < -J_d/2J_b, where J_d and J_b are the magnitudes of the disc and binary angular momentum vectors respectively. If however theta > 90 degrees and this criterion is not satisfied the same disc may counteralign its inner regions and, on longer timescales, coalign its outer regions. I show that for typical disc parameters, describing an accretion event on to a supermassive black hole binary, a misaligned circumbinary disc is likely to wholly co-- or counter--align with the binary plane. This is because the binary angular momentum dominates the disc angular momentum. However with extreme parameters (binary mass ratio M_2/M_1 << 1 or binary eccentricity e ~ 1) the same disc may simultaneously co- and counter-align. It is known that coplanar prograde circumbinary discs are stable. I show that coplanar retrograde circumbinary discs are also stable. A chaotic accretion event on to an SMBH binary will therefore result in a coplanar circumbinary disc that is either prograde or retrograde with respect to the binary plane.


title: Matter and Antimatter in the Universe

abstract: We review observational evidence for a matter-antimatter asymmetry in the early universe, which leads to the remnant matter density we observe today. We also discuss observational bounds on the presence of antimatter in the present day universe, including the possibility of a large lepton asymmetry in the cosmic neutrino background. We briefly review the theoretical framework within which baryogenesis, the dynamical generation of a matter-antimatter asymmetry, can occur. As an example, we discuss a testable minimal model that simultaneously explains the baryon asymmetry of the universe, neutrino oscillations and dark matter.


title: Stability of FFLO states in optical lattices with bilayer structure

abstract: We investigate the stability of the superfluid state in a bilayer fermionic optical lattice system with a confining potential, using the Bogoliubov de-Gennes equations. It is clarified that in the imbalanced case, the introduction of the interlayer hopping stabilizes the radial Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) state, while makes the angular FFLO state unstable. We also discuss the system size dependence of the superfluid ground state. It is clarified that in a certain ring region the A-FFLO state is indeed realized in a large system.


title: Metal-Insulator transition in optical lattice system with site-dependent interactions

abstract: We investigate the half-filled Hubbard model with spatially alternating interactions by means of the two-site dynamical mean-field theory. It is found that a single Mott transition occurs when two kinds of interactions are increased. This implies that the different interactions are essentially irrelevant at the critical point. The nature of the Mott states is also addressed.


title: Holographic Description of Finite Size Effects in Strongly Coupled Superconductors

abstract: Despite its fundamental and practical interest, the understanding of mesoscopic effects in strongly coupled superconductors is still limited. Here we address this problem by studying holographic superconductivity in a disk and a strip of typical size $\ell$. For $\ell < \ell_c$, where $\ell_c$ depends on the chemical potential and temperature, we have found that the order parameter vanishes. The superconductor-metal transition at $\ell = \ell_c$ is controlled by mean-field critical exponents which suggests that quantum and thermal fluctuations induced by finite size effects are suppressed in holographic superconductors. Intriguingly, the effective interactions that bind the order parameter increases as $\ell$ decreases. Most of these results are consistent with experimental observations in Pb nanograins at low temperature and qualitatively different from the ones expected in a weakly coupled superconductor.


title: Milky Way Star Forming Complexes and the Turbulent Motion of the Galaxy's Molecular Gas

abstract: We analyze Spitzer GLIMPSE, MSX, and WMAP images of the Milky Way to identify 8 micron and free-free sources in the Galaxy. Seventy-two of the eighty-eight WMAP sources have coverage in the GLIMPSE and MSX surveys suitable for identifying massive star forming complexes (SFC). We measure the ionizing luminosity functions of the SFCs and study their role in the turbulent motion of the Galaxy's molecular gas. We find a total Galactic free-free flux f_{\nu} = 46177.6 Jy; the 72 WMAP sources with full 8 micron coverage account for 34263.5 Jy (~75%), with both measurements made at \nu=94GHz (W band). We find a total of 280 SFCs, of which 168 have unique kinematic distances and free-free luminosities. We use a simple model for the radial distribution of star formation to estimate the free-free and ionizing luminosity for the sources lacking distance determinations. The total dust-corrected ionizing luminosity is Q = 2.9 \pm 0.5 x 10^53 photons s^-1, which implies a galactic star formation rate of 1.2 \pm 0.2 M_{\sun} yr^-1. We present the (ionizing) luminosity function of the SFCs, and show that 24 sources emit half the ionizing luminosity of the Galaxy. The SFCs appear as bubbles in GLIMPSE or MSX images; the radial velocities associated with the bubble walls allow us to infer the expansion velocity of the bubbles. We calculate the kinetic luminosity of the bubble expansion and compare it to the turbulent luminosity of the inner molecular disk. SFCs emitting 80% of the total galactic free-free luminosity produce a kinetic luminosity equal to 65% of the turbulent luminosity in the inner molecular disk. This suggests that the expansion of the bubbles is a major driver of the turbulent motion of the inner Milky Way molecular gas.


title: Universality of modulation length (and time) exponents

abstract: We study systems with a crossover parameter lambda, such as the temperature T, which has a threshold value lambda* across which the correlation function changes from exhibiting fixed wavelength (or time period) modulations to continuously varying modulation lengths (or times). We report on a new exponent, nuL, characterizing the universal nature of this crossover. These exponents, similar to standard correlation length exponents, are obtained from motion of the poles of the momentum (or frequency) space correlation functions in the complex k-plane (or omega-plane) as the parameter lambda is varied. Near the crossover, the characteristic modulation wave-vector KR on the variable modulation length "phase" is related to that on the fixed modulation length side, q via |KR-q|\propto|T-T*|^{nuL}. We find, in general, that nuL=1/2. In some special instances, nuL may attain other rational values. We extend this result to general problems in which the eigenvalue of an operator or a pole characterizing general response functions may attain a constant real (or imaginary) part beyond a particular threshold value, lambda*. We discuss extensions of this result to multiple other arenas. These include the ANNNI model. By extending our considerations, we comment on relations pertaining not only to the modulation lengths (or times) but also to the standard correlation lengths (or times). We introduce the notion of a Josephson timescale. We comment on the presence of "chaotic" modulations in "soft-spin" and other systems. These relate to glass type features. We discuss applications to Fermi systems - with particular application to metal to band insulator transitions, change of Fermi surface topology, divergent effective masses, Dirac systems, and topological insulators. Both regular periodic and glassy (and spatially chaotic behavior) may be found in strongly correlated electronic systems.


title: The Omega Deformation From String and M-Theory

abstract: We present a string theory construction of Omega-deformed four-dimensional gauge theories with generic values of \epsilon_1 and \epsilon_2. Our solution gives an explicit description of the geometry in the core of Nekrasov and Witten's realization of the instanton partition function, far from the asymptotic region of their background. This construction lifts naturally to M-theory and corresponds to an M5-brane wrapped on a Riemann surface with a selfdual flux. Via a 9-11 flip, we finally reinterpret the Omega deformation in terms of non-commutative geometry. Our solution generates all modified couplings of the \Omega-deformed gauge theory, and also yields a geometric origin for the quantum spectral curve of the associated quantum integrable system.


title: Quasi-isometric co-Hopficity of non-uniform lattices in rank-one semi-simple Lie groups

abstract: We prove that if $G$ is a non-uniform lattice in a rank-one semi-simple Lie group $\ne Isom(\H^2_\R)$ then $G$ is quasi-isometrically co-Hopf. This means that every quasi-isometric embedding $G\to G$ is coarsely onto and thus is a quasi-isometry.


title: CANDELS: Correlations of SEDs and Morphologies with Star-formation Status for Massive Galaxies at z ~ 2

abstract: We present a study on Spectral Energy Distributions, Morphologies, and star formation for an IRAC-selected extremely red object sample in the GOODS Chandra Deep Field-South. This work was enabled by new HST/WFC3 near-IR imaging from the CANDELS survey as well as the deepest available X-ray data from Chandra 4 Ms observations. This sample consists of 133 objects with the 3.6um limiting magnitude of [3.6] = 21.5, and is approximately complete for galaxies with M >10^{11}M_sun at 1.5 < z < 2.5. We classify this sample into two types, quiescent and star-forming galaxies, in the observed infrared color-color ([3.6]-[24] vs K-[3.6]) diagram. The further morphological study of this sample show a consistent result with the observed color classification. The classified quiescent galaxies are bulge dominated and star-forming galaxies in the sample have disk or irregular morphologies. Our observed infrared color classification is also consistent with the rest-frame color (U-V vs V-J) classification. We also found that quiescent and star-forming galaxies are well separated in the nonparametric morphology parameter (Gini vs M_{20}) diagram measuring their concentration and clumpiness: quiescent galaxies have Gini coefficient higher than 0.58 and star forming galaxies have Gini coefficient lower that 0.58. We argue that the star formation quenching process must lead to or be accompanied by the increasing galaxy concentration. One prominent morphological feature of this sample is that disks are commonly seen in this massive galaxy sample at 1.5 < z < 2.5: 30% of quiescent galaxies and 70% of star forming galaxies with M >10^{11}M_sun have disks in their rest-frame optical morphologies. The prevalence of these extended, relatively undisturbed disks challenges the merging scenario as the main mode of massive galaxy formation.


title: Six-dimensional Supergravity and Projective Superfields

abstract: We propose a superspace formulation of N=(1,0) conformal supergravity in six dimensions. The corresponding superspace constraints are invariant under super-Weyl transformations generated by a real scalar parameter. The known variant Weyl super-multiplet is recovered by coupling the geometry to a super-3-form tensor multiplet. Isotwistor variables are introduced and used to define projective superfields. We formulate a locally supersymmetric and super-Weyl invariant action principle in projective superspace. Some families of dynamical supergravity-matter systems are presented.


title: The Nature of the Compton-thick X-ray Reprocessor in NGC 4945

abstract: We present an exhaustive methodology for fitting Compton-thick X-ray reprocessor models to obscured AGNs and for interpreting the results. We focus on the MYTORUS model but also utilize other models. We apply the techniques to Suzaku, BeppoSAX, and Swift BAT spectra of the Sy 2 galaxy NGC 4945, but the methods are applicable to other AGNs including Compton-thin sources. The models overcome a major restriction of disk-reflection models, namely the assumption of an infinite column density. Finite column-density models produce a richer variety of spectral shapes and characteristics, even for Compton-thin AGNs. Although NGC 4945 is one of the brightest AGNs above 10 keV, the models span nearly a factor of 3 in column density (~2 to 6 x 10^{24} cm^{-2}) and 2 orders of magnitude in the intrinsic 2-195 keV luminosity. Models in which the continuum above 10 keV is dominated by the direct (unscattered) continuum or Compton-scattered continuum give the highest and lowest intrinsic luminosities respectively. Variability properties favor solutions in which the unscattered continuum dominates above 10 keV. The data require that the Compton-scattered continuum and Fe Kalpha line emission come predominantly from the illuminated surfaces of the X-ray reprocessor, implying a clumpy medium with a global covering factor that is small enough that the Compton-scattered continuum does not dominate the spectrum above 10 keV. This can be identified with the ~30 pc region spatially resolved by Chandra. The implied intrinsic bolometric luminosity is close to, or greater than, the Eddington luminosity. However, a strongly beamed AGN embedded in a shell of Compton-thick (but clumpy) matter requires less fine-tuning of the covering factor. Beaming is consistent with recent radio and Fermi results. Such beamed Compton-thick AGNs would be preferentially selected in surveys over unbeamed Compton-thick AGNs.


title: Magneto-seismology of solar atmospheric loops by means of longitudinal oscillations

abstract: There is increasingly strong observational evidence that slow magnetoacoustic modes arise in the solar atmosphere. Solar magneto-seismology is a novel tool to derive otherwise directly un-measurable properties of the solar atmosphere when magnetohydrodynamic (MHD) wave theory is compared to wave observations. Here, MHD wave theory is further developed illustrating how information about the magnetic and density structure along coronal loops can be determined by measuring the frequencies of the slow MHD oscillations. The application to observations of slow magnetoacoustic waves in coronal loops is discused.


title: A remark on representations of infinite symmetric groups

abstract: We simplify construction of Thoma representations of an infinite symmetric group


title: Constrained Supersymmetry after two years of LHC data: a global view with Fittino

abstract: We perform global fits to the parameters of the Constrained Minimal Supersymmetric Standard Model (CMSSM) and to a variant with non-universal Higgs masses (NUHM1). In addition to constraints from low-energy precision observables and the cosmological dark matter density, we take into account the LHC exclusions from searches in jets plus missing transverse energy signatures with about 5\,fb$^{-1}$ of integrated luminosity. We also include the most recent upper bound on the branching ratio $B_s\to\mu\mu$ from LHCb. Furthermore, constraints from and implications for direct and indirect dark matter searches are discussed. The best fit of the CMSSM prefers a light Higgs boson just above the experimentally excluded mass. We find that the description of the low-energy observables, $(g-2)_{\mu}$ in particular, and the non-observation of SUSY at the LHC become more and more incompatible within the CMSSM. A potential SM-like Higgs boson with mass around 126 GeV can barely be accommodated. Values for ${\cal B}(B_s\to\mu\mu)$ just around the Standard Model prediction are naturally expected in the best fit region. The most-preferred region is not yet affected by limits on direct WIMP searches, but the next generation of experiments will probe this region. Finally, we discuss implications from fine-tuning for the best fit regions.


title: Discrete Dynamical Genetic Programming in XCS

abstract: A number of representation schemes have been presented for use within Learning Classifier Systems, ranging from binary encodings to neural networks. This paper presents results from an investigation into using a discrete dynamical system representation within the XCS Learning Classifier System. In particular, asynchronous random Boolean networks are used to represent the traditional condition-action production system rules. It is shown possible to use self-adaptive, open-ended evolution to design an ensemble of such discrete dynamical systems within XCS to solve a number of well-known test problems.


title: Longitudinal oscillations in density stratified and expanding solar waveguides

abstract: Waves and oscillations can provide vital information about the internal structure of waveguides they propagate in. Here, we analytically investigate the effects of density and magnetic stratification on linear longitudinal magnetohydrodynamic (MHD) waves. The focus of this paper is to study the eigenmodes of these oscillations. It is our specific aim is to understand what happens to these MHD waves generated in flux tubes with non-constant (e.g., expanding or magnetic bottle) cross-sectional area and density variations. The governing equation of the longitudinal mode is derived and solved analytically and numerically. In particular, the limit of the thin flux tube approximation is examined. The general solution describing the slow longitudinal MHD waves in an expanding magnetic flux tube with constant density is found. Longitudinal MHD waves in density stratified loops with constant magnetic field are also analyzed. From analytical solutions, the frequency ratio of the first overtone and fundamental mode is investigated in stratified waveguides. For small expansion, a linear dependence between the frequency ratio and the expansion factor is found. From numerical calculations it was found that the frequency ratio strongly depends on the density profile chosen and, in general, the numerical results are in agreement with the analytical results. The relevance of these results for solar magneto-seismology is discussed.


title: Fuzzy Dynamical Genetic Programming in XCSF

abstract: A number of representation schemes have been presented for use within Learning Classifier Systems, ranging from binary encodings to Neural Networks, and more recently Dynamical Genetic Programming (DGP). This paper presents results from an investigation into using a fuzzy DGP representation within the XCSF Learning Classifier System. In particular, asynchronous Fuzzy Logic Networks are used to represent the traditional condition-action production system rules. It is shown possible to use self-adaptive, open-ended evolution to design an ensemble of such fuzzy dynamical systems within XCSF to solve several well-known continuous-valued test problems.


title: Integrated high frequency aluminum nitride optomechanical resonators

abstract: Aluminum nitride (AlN) has been widely used in microeletromechanical resonators for its excellent electromechanical properties. Here we demonstrate the use of AlN as an optomechanical material that simultaneously offer low optical and mechanical loss. Integrated AlN microring resonators in the shape of suspended rings exhibit high optical quality factor (Q) with loaded Q up to 125,000. Optomechanical transduction of the Brownian motion of a GHz contour mode yields a displacement sensitivity of 6.2\times10^(-18)m/Hz^(1/2) in ambient air.


title: Tilings with $n$-Dimensional Chairs and their Applications to Asymmetric Codes

abstract: An $n$-dimensional chair consists of an $n$-dimensional box from which a smaller $n$-dimensional box is removed. A tiling of an $n$-dimensional chair has two nice applications in coding for write-once memories. The first one is in the design of codes which correct asymmetric errors with limited-magnitude. The second one is in the design of $n$ cells $q$-ary write-once memory codes. We show an equivalence between the design of a tiling with an integer lattice and the design of a tiling from a generalization of splitting (or of Sidon sequences). A tiling of an $n$-dimensional chair can define a perfect code for correcting asymmetric errors with limited-magnitude. We present constructions for such tilings and prove cases where perfect codes for these type of errors do not exist.


title: Observation of eta_b(2S) in Upsilon(2S) -> gamma eta_b(2S), eta_b(2S) -> hadrons, and Confirmation of eta_b(1S)

abstract: The data for 9.3 million Upsilon(2S) and 20.9 million Upsilon(1S) taken with the CLEO III detector has been used to study the radiative population of states identified by their decay into twenty six different exclusive hadronic final states. In the Upsilon(2S) decays an enhancement is observed at a ~5 sigma level at a mass of 9974.6+-2.3(stat)+-2.1(syst) MeV. It is attributed to eta_b(2S), and corresponds to the Upsilon(2S) hyperfine splitting of 48.7+-2.3(stat)+-2.1(syst) MeV. In the Upsilon(1S) decays, the identification of eta_b(1S) is confirmed at a ~3 sigma level with M(eta_b(1S)) in agreement with its known value.


title: Revealing companions to nearby stars with astrometric acceleration

abstract: A subset of 51 Hipparcos astrometric binaries among FG dwarfs within 67pc has been surveyed with the NICI adaptive optics system at Gemini-S, directly resolving for the first time 17 sub-arcsecond companions and 7 wider ones. Using these data together with published speckle interferometry of 57 stars, we compare the statistics of resolved astrometric companions with those of a simulated binary population. The fraction of resolved companions is slightly lower than expected from binary statistics. About 10% of astrometric companions could be "dark" (white dwarfs and close pairs of late M-dwarfs). To our surprise, several binaries are found with companions too wide to explain the acceleration. Re-analysis of selected intermediate astrometric data shows that some acceleration solutions in the original Hipparcos catalog are spurious.


title: On Soft Limits of Inflationary Correlation Functions

abstract: Soft limits of inflationary correlation functions are both observationally relevant and theoretically robust. Various theorems can be proven about them that are insensitive to detailed model-building assumptions. In this paper, we re-derive several of these theorems in a universal way. Our method makes manifest why soft limits are such an interesting probe of the spectrum of additional light fields during inflation. We illustrate these abstract results with a detailed case study of the soft limits of quasi-single-field inflation.


title: Maximal R-symmetry violating amplitudes in type IIB superstring theory

abstract: On-shell superspace techniques are used to quantify R-symmetry violation in type IIB superstring theory amplitudes in a flat background in ten dimensions. This shows the existence of a particularly simple class of non-vanishing amplitudes in this theory which violate R-symmetry maximally. General properties of the class and some of its extensions are established which at string tree level are shown to determine the first three non-trivial effective field theory contributions to all multiplicity. This leads to a natural conjecture for the exact analytic part of the first two of these.


title: Folded Codes from Function Field Towers and Improved Optimal Rate List Decoding

abstract: We give a new construction of algebraic codes which are efficiently list decodable from a fraction $1-R-\eps$ of adversarial errors where $R$ is the rate of the code, for any desired positive constant $\eps$. The worst-case list size output by the algorithm is $O(1/\eps)$, matching the existential bound for random codes up to constant factors. Further, the alphabet size of the codes is a constant depending only on $\eps$ - it can be made $\exp(\tilde{O}(1/\eps^2))$ which is not much worse than the lower bound of $\exp(\Omega(1/\eps))$. The parameters we achieve are thus quite close to the existential bounds in all three aspects - error-correction radius, alphabet size, and list-size - simultaneously. Our code construction is Monte Carlo and has the claimed list decoding property with high probability. Once the code is (efficiently) sampled, the encoding/decoding algorithms are deterministic with a running time $O_\eps(N^c)$ for an absolute constant $c$, where $N$ is the code's block length. Our construction is based on a linear-algebraic approach to list decoding folded codes from towers of function fields, and combining it with a special form of subspace-evasive sets. Instantiating this with the explicit "asymptotically good" Garcia-Stichtenoth tower of function fields yields the above parameters. To illustrate the method in a simpler setting, we also present a construction based on Hermitian function fields, which offers similar guarantees with a list and alphabet size polylogarithmic in the block length $N$. Along the way, we shed light on how to use automorphisms of certain function fields to enable list decoding of the folded version of the associated algebraic-geometric codes.


title: On the Link Between Central Black Holes, Bar Dynamics, and Dark Matter Halos in Spiral Galaxies

abstract: The discovery of a relationship between supermassive black hole (SMBH) mass and spiral arm pitch angle (P) is evidence that SMBHs are tied to the overall secular evolution of a galaxy. The discovery of SMBHs in late-type galaxies with little or no bulge suggests that an underlying correlation between the dark matter halo concentration and SMBH mass (MBH) exists, rather than between the bulge mass and MBH. In this paper we measure P using a two-dimensional fast fourier transform and estimate the bar pattern speeds of 40 barred spiral galaxies from the Carnegie-Irvine Galaxy Survey. The pattern speeds were derived by estimating the gravitational potentials of our galaxies from Ks-band images and using them to produce dynamical simulation models. The pattern speeds allow us to identify those galaxies with low central dark halo densities, or fast rotating bars, while P provides an estimate of MBH. We find that a wide range of MBH exists in galaxies with low central dark matter halo densities, which appears to support other theoretical results. We also find that galaxies with low central dark halo densities appear to follow more predictable trends in P versus de Vaucouleurs morphological type (T) and bar strength versus T than barred galaxies in general. The empirical relationship between MBH and total gravitational mass of a galaxy (Mtot) allows us to predict the minimum Mtot that will be observationally measured of our fast bar galaxies. These predictions will be investigated in a subsequent paper.


title: Wave Functions, Quantum Diffusion, and Scaling Exponents in Golden-Mean Quasiperiodic Tilings

abstract: We study the properties of wave functions and the wave-packet dynamics in quasiperiodic tight-binding models in one, two, and three dimensions. The atoms in the one-dimensional quasiperiodic chains are coupled by weak and strong bonds aligned according to the Fibonacci sequence. The associated d-dimensional quasiperiodic tilings are constructed from the direct product of d such chains, which yields either the hypercubic tiling or the labyrinth tiling. This approach allows us to consider rather large systems numerically. We show that the wave functions of the system are multifractal and that their properties can be related to the structure of the system in the regime of strong quasiperiodic modulation by a renormalization group (RG) approach. We also study the dynamics of wave packets to get information about the electronic transport properties. In particular, we investigate the scaling behaviour of the return probability of the wave packet with time. Applying again the RG approach we show that in the regime of strong quasiperiodic modulation the return probability is governed by the underlying quasiperiodic structure. Further, we also discuss lower bounds for the scaling exponent of the width of the wave packet and propose a modified lower bound for the absolute continuous regime.


title: Observation of the fractional ac Josephson effect: the signature of Majorana particles

abstract: Topological superconductors which support Majorana fermions are thought to be realized in one-dimensional semiconducting wires coupled to a superconductor. Such excitations are expected to exhibit non-Abelian statistics and can be used to realize quantum gates that are topologically protected from local sources of decoherence. Here we report the observation of the fractional a.c. Josephson effect in a hybrid semiconductor/superconductor InSb/Nb nanowire junction, a hallmark of topological matter. When the junction is irradiated with a radio-frequency f in the absence of an external magnetic field, quantized voltage steps (Shapiro steps) with a height hf/2e are observed, as is expected for conventional superconductor junctions, where the supercurrent is carried by charge-2e Cooper pairs. At high magnetic fields the height of the first Shapiro step is doubled to hf/e, suggesting that the supercurrent is carried by charge-e quasiparticles. This is a unique signature of Majorana fermions, elusive particles predicted ca. 80 years ago.


title: First Evidence of Globular Cluster Formation from the Ejecta of Prompt Type Ia Supernovae

abstract: Recent spectroscopic observations of globular clusters (GCs) in the Large Magellanic Cloud (LMC) have discovered that one of the intermediate-age GC, NGC 1718 with [Fe/H]=-0.7 has an extremely low [Mg/Fe] ratio of ~-0.9. We propose that NGC 1718 was formed from the ejecta of type Ia supernovae (SNe Ia) mixed with very metal-poor ([Fe/H] <-1.3) gas about ~ 2 Gyr ago. The proposed scenario is shown to be consistent with the observed abundances of Fe-group elements such as Cr, Mn, and Ni. In addition, compelling evidence for asymptotic giant branch stars playing a role in chemical enrichment during this GC formation is found. We suggest that the origin of the metal-poor gas is closely associated with the efficient gas-transfer from the outer gas disk of the Small Magellanic Cloud to the LMC disk. We anticipate that the outer part of the LMC disk contains field stars exhibiting significantly low [Mg/Fe] ratios, formed through the same process as NGC 1718.


title: Outbreak size distributions in epidemics with multiple stages

abstract: Multiple-type branching processes that model the spread of infectious diseases are investigated. In these stochastic processes, the disease goes through multiple stages before it eventually disappears. We mostly focus on the critical multistage Susceptible-Infected-Recovered (SIR) infection process. In the infinite population limit, we compute the outbreak size distributions and show that asymptotic results apply to more general multiple-type critical branching processes. Finally using heuristic arguments and simulations we establish scaling laws for a multistage SIR model in a finite population.


title: A low phase noise microwave source for atomic spin squeezing experiments in 87Rb

abstract: We describe and characterize a simple, low cost, low phase noise microwave source that operates near 6.800 GHz for agile, coherent manipulation of ensembles of 87Rb. Low phase noise is achieved by directly multiplying a low phase noise 100 MHz crystal to 6.8 GHz using a non-linear transmission line and filtering the output with custom band-pass filters. The fixed frequency signal is single sideband modulated with a direct digital synthesis frequency source to provide the desired phase, amplitude, and frequency control. Before modulation, the source has a single sideband phase noise near -140 dBc/Hz in the range of 10 kHz to 1 MHz offset from the carrier frequency and -130 dBc/Hz after modulation. The resulting source is estimated to contribute added spin-noise variance 16 dB below the quantum projection noise level during quantum nondemolition measurements of the clock transition in an ensemble 7x10^5 87Rb atoms.


title: Fast detection of single-charge tunneling to a graphene quantum dot in a multi-level regime

abstract: In situ tunable radio-frequency charge detection is used for the determination of the tunneling rates into and out of a graphene single quantum dot connected to only one lead. An analytical model for calculating these rates in the multi-level tunneling regime is presented and found to correspond very well to our experimental observations.


title: Probing the Neutrino Mass Hierarchy with Super-Kamiokande

abstract: We show that for recently discovered large values of theta(13), a superbeam with an average neutrino energy of ~ 5 GeV, such as those being proposed at CERN, if pointing to Super-Kamiokande (L = 8770 km), could reveal the neutrino mass hierarchy at 5 sigma in less than two years irrespective of the true hierarchy and CP phase. The measurement relies on the near resonant matter effect in the numu to nue oscillation channel, and can be done counting the total number of appearance events with just a neutrino beam.


title: Optimal strategies for continuous gravitational wave detection in pulsar timing arrays

abstract: Supermassive black hole binaries (SMBHBs) are expected to emit continuous gravitational waves in the pulsar timing array (PTA) frequency band ($10^{-9}$--$10^{-7}$ Hz). The development of data analysis techniques aimed at efficient detection and characterization of these signals is critical to the gravitational wave detection effort. In this paper we leverage methods developed for LIGO continuous wave gravitational searches, and explore the use of the $\mathcal{F}$-statistic for such searches in pulsar timing data. Babak & Sesana 2012 have already used this approach in the context of PTAs to show that one can resolve multiple SMBHB sources in the sky. Our work improves on several aspects of prior continuous wave search methods developed for PTA data analysis. The algorithm is implemented fully in the time domain, which naturally deals with the irregular sampling typical of PTA data and avoids spectral leakage problems associated with frequency domain methods. We take into account the fitting of the timing model, and have generalized our approach to deal with both correlated and uncorrelated colored noise sources. We also develop an incoherent detection statistic that maximizes over all pulsar dependent contributions to the likelihood. To test the effectiveness and sensitivity of our detection statistics, we perform a number of monte-carlo simulations. We produce sensitivity curves for PTAs of various configurations, and outline an implementation of a fully functional data analysis pipeline. Finally, we present a derivation of the likelihood maximized over the gravitational wave phases at the pulsar locations, which results in a vast reduction of the search parameter space.


title: An Absence of Neutrinos Associated with Cosmic Ray Acceleration in Gamma-Ray Bursts

abstract: Gamma-Ray Bursts (GRBs) have been proposed as a leading candidate for acceleration of ultra high-energy cosmic rays, which would be accompanied by emission of TeV neutrinos produced in proton-photon interactions during acceleration in the GRB fireball. Two analyses using data from two years of the IceCube detector produced no evidence for this neutrino emission, placing strong constraints on models of neutrino and cosmic-ray production in these sources.


title: Fully pseudospectral time evolution and its application to 1+1 dimensional physical problems

abstract: It was recently demonstrated that time-dependent PDE problems can numerically be solved with a fully pseudospectral scheme, i.e. using spectral expansions with respect to both spatial and time directions (Hennig and Ansorg, 2009 [15]). This was done with the example of simple scalar wave equations in Minkowski spacetime. Here we show that the method can be used to study interesting physical problems that are described by systems of nonlinear PDEs. To this end we consider two 1+1 dimensional problems: radial oscillations of spherically symmetric Newtonian stars and time evolution of Gowdy spacetimes as particular cosmological models in general relativity.


title: Magic-state distillation with the four-qubit code

abstract: The distillation of magic states is an often-cited technique for enabling universal quantum computing once the error probability for a special subset of gates has been made negligible by other means. We present a routine for magic-state distillation that reduces the required overhead for a range of parameters of practical interest. Each iteration of the routine uses a four-qubit error-detecting code to distill the +1 eigenstate of the Hadamard gate at a cost of ten input states per two improved output states. Use of this routine in combination with the 15-to-1 distillation routine described by Bravyi and Kitaev allows for further improvements in overhead.


title: The Holographic Dual of 2+1 Dimensional QFTs with N=1 SUSY and Massive Fundamental Flavours

abstract: The Maldacena Nastase solution is generalised to include massive fundamental matter through the addition of a flavour profile. This gives a holographic dual to N=1 SYM-CS with massive fundamental matter with a singularity free IR. We study this solution in some detail confirming confinement and asymptotic freedom. A recently proposed solution generating technique is then applied which results in a new type-IIA supergravity solution. In a certain limit the geometry of this solution is asymptotically AdS_4X Y, where Y is the metric at the base of the Bryant-Salamon G_2 cone, which has topology S^3XS^3.


title: Improved Quantum LDPC Decoding Strategies For The Misidentified Quantum Depolarizing Channel

abstract: Quantum cryptography via key distribution mechanisms that utilize quantum entanglement between sender-receiver pairs will form the basis of future large-scale quantum networks. A key engineering challenge in such networks will be the ability to correct for decoherence effects in the distributed entanglement resources. It is widely believed that sophisticated quantum error correction codes, such as quantum low-density parity-check (LDPC) codes, will be pivotal in such a role. However, recently the importance of the channel mismatch effect in degrading the performance of deployed quantum LDPC codes has been pointed out. In this work we help remedy this situation by proposing new quantum LDPC decoding strategies that can significantly reduce performance degradation by as much as $50\%$. Our new strategies for the quantum LDPC decoder are based on previous insights from classical LDPC decoders in mismatched channels, where an asymmetry in performance is known as a function of the estimated channel noise. We show how similar asymmetries carry over to the quantum depolarizing channel, and how an estimate of the depolarization flip parameter weighted to larger values leads to significant performance improvement.


title: Software Mutational Robustness

abstract: Neutral landscapes and mutational robustness are believed to be important enablers of evolvability in biology. We apply these concepts to software, defining mutational robustness to be the fraction of random mutations that leave a program's behavior unchanged. Test cases are used to measure program behavior and mutation operators are taken from genetic programming. Although software is often viewed as brittle, with small changes leading to catastrophic changes in behavior, our results show surprising robustness in the face of random software mutations. The paper describes empirical studies of the mutational robustness of 22 programs, including 14 production software projects, the Siemens benchmarks, and 4 specially constructed programs. We find that over 30% of random mutations are neutral with respect to their test suite. The results hold across all classes of programs, for mutations at both the source code and assembly instruction levels, across various programming languages, and are only weakly related to test suite coverage. We conclude that mutational robustness is an inherent property of software, and that neutral variants (i.e., those that pass the test suite) often fulfill the program's original purpose or specification. Based on these results, we conjecture that neutral mutations can be leveraged as a mechanism for generating software diversity. We demonstrate this idea by generating a population of neutral program variants and showing that the variants automatically repair unknown bugs with high probability. Neutral landscapes also provide a partial explanation for recent results that use evolutionary computation to automatically repair software bugs.


title: Gauge and spacetime connections in the Plebanski formulation of complex general relativity

abstract: The Plebanski formulation of complex general relativity is given in terms of variables valued in the complexification of the $so(3)$ Lie algebra. Therefore, it is genuinely a gauge theory that is also diffeomorphism-invariant. For this reason, the way that the Levi-Civita connection emerges from this formulation is not direct because both the internal (gauge) and the spacetime connections are geometrical structures a priori not related, there is not a natural link between them. Any possible relationship must be put in by hand or must come from extra hypotheses. In this paper, we analyze the correct relationship between these connections and show how they are related.


title: Majorana Fermions and Odd-frequency Cooper Pairs in a Nano Wire

abstract: We discuss a strong relationship between Majorana fermions and odd-frequency Cooper pairs which appear at a disordered normal (N) nano wire attached to a topologically nontrivial superconducting (S) one. The transport properties in superconducting nano wire junctions show universal behaviors irrespective of the degree of disorder: the quantized zero-bias differential conductance at 2e^{2}/h in NS junction and the fractional current-phase (J-\varphi) relationship of the Josephson effect in SNS junction J\propto \sin(\varphi/2). Such behaviors are exactly the same as those in the anomalous proximity effect of odd-parity spin-triplet superconductors. We show that odd-frequency Cooper pairs support the universal transport properties. The odd-frequency pairs exist wherever the Majorana fermions stay.


title: Estimating Unknown Sparsity in Compressed Sensing

abstract: In the theory of compressed sensing (CS), the sparsity ||x||_0 of the unknown signal x\in\R^p is commonly assumed to be a known parameter. However, it is typically unknown in practice. Due to the fact that many aspects of CS depend on knowing ||x||_0, it is important to estimate this parameter in a data-driven way. A second practical concern is that ||x||_0 is a highly unstable function of x. In particular, for real signals with entries not exactly equal to 0, the value ||x||_0=p is not a useful description of the effective number of coordinates. In this paper, we propose to estimate a stable measure of sparsity s(x):=||x||_1^2/||x||_2^2, which is a sharp lower bound on ||x||_0. Our estimation procedure uses only a small number of linear measurements, does not rely on any sparsity assumptions, and requires very little computation. A confidence interval for s(x) is provided, and its width is shown to have no dependence on the signal dimension p. Moreover, this result extends naturally to the matrix recovery setting, where a soft version of matrix rank can be estimated with analogous guarantees. Finally, we show that the use of randomized measurements is essential to estimating s(x). This is accomplished by proving that the minimax risk for estimating s(x) with deterministic measurements is large when n<<p.

categories: cs.IT math.IT math.ST stat.ME stat.ML stat.TH


title: Fixed-smoothing asymptotics for time series

abstract: In this paper, we derive higher order Edgeworth expansions for the finite sample distributions of the subsampling-based t-statistic and the Wald statistic in the Gaussian location model under the so-called fixed-smoothing paradigm. In particular, we show that the error of asymptotic approximation is at the order of the reciprocal of the sample size and obtain explicit forms for the leading error terms in the expansions. The results are used to justify the second-order correctness of a new bootstrap method, the Gaussian dependent bootstrap, in the context of Gaussian location model.


title: Compressibility of an Ultracold Fermi Gas with Repulsive Interactions

abstract: Fermi gases with repulsive interactions are characterized by measuring their compressibility as a function of interaction strength. The compressibility is obtained from in-trap density distributions monitored by phase contrast imaging. For interaction parameters k_F a > 0.25 fast decay of the gas prevents the observation of equilibrium profiles. For smaller interaction parameters, the results are adequately described by first-order perturbation theory. A novel phase contrast imaging method compensates for dispersive distortions of the images.


title: Planar F-Deletion: Approximation, Kernelization and Optimal FPT Algorithms

abstract: Let F be a finite set of graphs. In the F-Deletion problem, we are given an n-vertex graph G and an integer k as input, and asked whether at most k vertices can be deleted from G such that the resulting graph does not contain a graph from F as a minor. F-Deletion is a generic problem and by selecting different sets of forbidden minors F, one can obtain various fundamental problems such as Vertex Cover, Feedback Vertex Set or Treewidth t-Deletion. In this paper we obtain a number of generic algorithmic results about Planar F-Deletion, when F contains at least one planar graph. The highlights of our work are - A constant factor approximation algorithm for the optimization version of Planar F-Deletion. - A randomized linear time and single exponential parameterized algorithm, that is, an algorithm running in time O(2^{O(k)} n), for the parameterized version of Planar F-deletion where all graphs in F are connected. The algorithm can be made deterministic at the cost of making the polynomial factor in the running time n*log^2 n rather than linear. - A polynomial kernel for parameterized Planar F-deletion These algorithms unify, generalize, and improve a multitude of results in the literature. Our main results have several direct applications, but also the methods we develop on the way have applicability beyond the scope of this paper. Our results -- constant factor approximation, polynomial kernelization and FPT algorithms -- are stringed together by a common theme of polynomial time preprocessing.


title: Neutrino Capture on $^{13}$C

abstract: We present neutrino cross sections on $^{13}$C. The charged-current cross sections leading to various states in the daughter $^{13}N$ and the neutral-current cross sections leading to various states in the daughter $^{13}$C are given. We also provide simple polynomial fits to those cross sections for quick estimates of the reaction rates. We briefly discuss possible implications for the current and future scintillator-based experiments.


title: Schauder Bases and Operator Theory III: Schauder Spectrums

abstract: In this paper, we study spectrums of Schauder operators. We show that we always can choose a Schauder operator in a given orbit such that the Schauder spectrum of it is empty.