Statistical mechanics of general discrete nonlinear Schr{\"o}dinger models: Localization transition and its relevance for Klein-Gordon lattices

We extend earlier work [Phys.Rev.Lett. 84, 3740 (2000)] on the statistical mechanics of the cubic one-dimensional discrete nonlinear Schrodinger (DNLS) equation to a more general class of models, including higher dimensionalities and nonlinearities of arbitrary degree. These extensions are physically motivated by the desire to describe situations with an excitation threshold for creation of localized excitations, as well as by recent work suggesting non-cubic DNLS models to describe Bose-Einstein condensates in deep optical lattices, taking into account the effective condensate dimensionality. Considering ensembles of initial conditions with given values of the two conserved quantities, norm and Hamiltonian, we calculate analytically the boundary of the 'normal' Gibbsian regime corresponding to infinite temperature, and perform numerical simulations to illuminate the nature of the localization dynamics outside this regime for various cases. Furthermore, we show quantitatively how this DNLS localization transition manifests itself for small-amplitude oscillations in generic Klein-Gordon lattices of weakly coupled anharmonic oscillators (in which energy is the only conserved quantity), and determine conditions for existence of persistent energy localization over large time scales.Comment: to be published in Physical Review

Classification of Big Data with Application to Imaging Genetics

Big data applications, such as medical imaging and genetics, typically generate datasets that consist of few observations n on many more variables p, a scenario that we denote as p>>n. Traditional data processing methods are often insufficient for extracting information out of big data. This calls for the development of new algorithms that can deal with the size, complexity, and the special structure of such datasets. In this paper, we consider the problem of classifying p>>n data and propose a classification method based on linear discriminant analysis (LDA). Traditional LDA depends on the covariance estimate of the data, but when p>>n the sample covariance estimate is singular. The proposed method estimates the covariance by using a sparse version of noisy principal component analysis (nPCA). The use of sparsity in this setting aims at automatically selecting variables that are relevant for classification. In experiments, the new method is compared to state-of-the art methods for big data problems using both simulated datasets and imaging genetics datasets

Energetically stable singular vortex cores in an atomic spin-1 Bose-Einstein condensate

We analyze the structure and stability of singular singly quantized vortices in a rotating spin-1 Bose-Einstein condensate. We show that the singular vortex can be energetically stable in both the ferromagnetic and polar phases despite the existence of a lower-energy nonsingular coreless vortex in the ferromagnetic phase. The spin-1 system exhibits energetic hierarchy of length scales resulting from different interaction strengths and we find that the vortex cores deform to a larger size determined by the characteristic length scale of the spin-dependent interaction. We show that in the ferromagnetic phase the resulting stable core structure, despite apparent complexity, can be identified as a single polar core with everywhere nonvanishing axially symmetric density profile. In the polar phase, the energetically favored core deformation leads to a splitting of a singly quantized vortex into a pair of half-quantum vortices that preserves the topology of the vortex outside the extended core region, but breaks the axial symmetry of the core. The resulting half-quantum vortices exhibit nonvanishing ferromagnetic cores.<br/

A Verified Certificate Checker for Finite-Precision Error Bounds in Coq and HOL4

Being able to soundly estimate roundoff errors of finite-precision computations is important for many applications in embedded systems and scientific computing. Due to the discrepancy between continuous reals and discrete finite-precision values, automated static analysis tools are highly valuable to estimate roundoff errors. The results, however, are only as correct as the implementations of the static analysis tools. This paper presents a formally verified and modular tool which fully automatically checks the correctness of finite-precision roundoff error bounds encoded in a certificate. We present implementations of certificate generation and checking for both Coq and HOL4 and evaluate it on a number of examples from the literature. The experiments use both in-logic evaluation of Coq and HOL4, and execution of extracted code outside of the logics: we benchmark Coq extracted unverified OCaml code and a CakeML-generated verified binary

Existence and Stability of Non-Trivial Scalar Field Configurations in Orbifolded Extra Dimensions

We consider the existence and stability of static configurations of a scalar field in a five dimensional spacetime in which the extra spatial dimension is compactified on an $S^1/Z_2$ orbifold. For a wide class of potentials with multiple minima there exist a finite number of such configurations, with total number depending on the size of the orbifold interval. However, a Sturm-Liouville stability analysis demonstrates that all such configurations with nodes in the interval are unstable. Nodeless static solutions, of which there may be more than one for a given potential, are far more interesting, and we present and prove a powerful general criterion that allows a simple determination of which of these nodeless solutions are stable. We demonstrate our general results by specializing to a number of specific examples, one of which may be analyzed entirely analytically.Comment: 23 pages, 7 figures, references added, factor of two corrected in kink energy definition, submitted to PR

Topological interface engineering and defect crossing in ultracold atomic gases

We propose an experimentally feasible scheme for topological interface engineering and show how it can be used for studies of dynamics of topologically nontrivial interfaces and perforation of defects and textures across such interfaces. The method makes use of the internal spin structure of the atoms together with locally applied control of interaction strengths to create many-particle states with highly complex topological properties. In particular, we consider a constructed coherent interface between topologically distinct phases of spinor Bose-Einstein condensates.Comment: 9 pages, 7 figure

Functional Big-step Semantics

When doing an interactive proof about a piece of software, it is important that the underlying programming language’s semantics does not make the proof unnecessarily difficult or unwieldy. Both smallstep and big-step semantics are commonly used, and the latter is typically given by an inductively defined relation. In this paper, we consider an alternative: using a recursive function akin to an interpreter for the language. The advantages include a better induction theorem, less duplication, accessibility to ordinary functional programmers, and the ease of doing symbolic simulation in proofs via rewriting. We believe that this style of semantics is well suited for compiler verification, including proofs of divergence preservation. We do not claim the invention of this style of semantics: our contribution here is to clarify its value, and to explain how it supports several language features that might appear to require a relational or small-step approach. We illustrate the technique on a simple imperative language with C-like for-loops and a break statement, and compare it to a variety of other approaches. We also provide ML and lambda-calculus based examples to illustrate its generality