title: Global existence for a coupled wave system related to the Strauss conjecture

abstract: A coupled system of semilinear wave equations is considered, and a small data global existence result related to the Strauss conjecture is proved. Previous results have shown that one of the powers may be reduced below the critical power for the Strauss conjecture provided the other power sufficiently exceeds such. The stability of such results under asymptotically flat perturbations of the space-time where an integrated local energy decay estimate is available is established.


title: A note on primitive $1-$normal elements over finite fields

abstract: Let $q$ be a prime power of a prime $p$, $n$ a positive integer and $\mathbb F_{q^n}$ the finite field with $q^n$ elements. The $k-$normal elements over finite fields were introduced and characterized by Huczynska et al (2013). Under the condition that $n$ is not divisible by $p$, they obtained an existence result on primitive $1-$normal elements of $\mathbb F_{q^n}$ over $\mathbb F_q$ for $q>2$. In this note, we extend their result to the excluded case $q=2$.


title: Rare Disease Physician Targeting: A Factor Graph Approach

abstract: In rare disease physician targeting, a major challenge is how to identify physicians who are treating diagnosed or underdiagnosed rare diseases patients. Rare diseases have extremely low incidence rate. For a specified rare disease, only a small number of patients are affected and a fractional of physicians are involved. The existing targeting methodologies, such as segmentation and profiling, are developed under mass market assumption. They are not suitable for rare disease market where the target classes are extremely imbalanced. The authors propose a graphical model approach to predict targets by jointly modeling physician and patient features from different data spaces and utilizing the extra relational information. Through an empirical example with medical claim and prescription data, the proposed approach demonstrates better accuracy in finding target physicians. The graph representation also provides visual interpretability of relationship among physicians and patients. The model can be extended to incorporate more complex dependency structures. This article contributes to the literature of exploring the benefit of utilizing relational dependencies among entities in healthcare industry.


title: Exact Diagonalization library for quantum electron models

abstract: We present an exact diagonalization C++ template library (EDLib) for solving quantum electron models, including single-band finite Hubbard cluster and multi-orbital impurity Anderson model. The observables that can be computed using EDLib are single particle Green's functions and spin-spin correlation functions. This code provides three different types of Hamiltonian matrix storage that can be chosen based on the model.


title: Effect of nuclear saturation parameters on possible maximum mass of neutron stars

abstract: In order to systematically examine the possible maximum mass of neutron stars, which is one of the important properties characterizing the physics in high-density region, I construct neutron star models by adopting phenomenological equations of state with various values of nuclear saturation parameters for low-density region, which are connected to the equation of state for high-density region characterized by the possible maximum sound velocity in medium. I derive an empirical formula for the possible maximum mass of neutron star. If massive neutron stars are observed, it could be possible to get a constraint on the possible maximum sound velocity for high-density region.


title: Simultaneous Confidence Band for Partially Linear Panel Data Models with Fixed Effects

abstract: In this paper, we construct the simultaneous confidence band (SCB) for the nonparametric component in partially linear panel data models with fixed effects. We remove the fixed effects, and further obtain the estimators of parametric and nonparametric components, which do not depend on the fixed effects. We establish the asymptotic distribution of their maximum absolute deviation between the estimated nonparametric component and the true nonparametric component under some suitable conditions, and hence the result can be used to construct the simultaneous confidence band of the nonparametric component. Based on the asymptotic distribution, it becomes difficult for the construction of the simultaneous confidence band. The reason is that the asymptotic distribution involves the estimators of the asymptotic bias and conditional variance, and the choice of the bandwidth for estimating the second derivative of nonparametric function. Clearly, these will cause computational burden and accumulative errors. To overcome these problems, we propose a Bootstrap method to construct simultaneous confidence band. Simulation studies indicate that the proposed Bootstrap method exhibits better performance under the limited samples.


title: NLP2Code: Code Snippet Content Assist via Natural Language Tasks

abstract: Developers increasingly take to the Internet for code snippets to integrate into their programs. To save developers the time required to switch from their development environments to a web browser in the quest for a suitable code snippet, we introduce NLP2Code, a content assist for code snippets. Unlike related tools, NLP2Code integrates directly into the source code editor and provides developers with a content assist feature to close the vocabulary gap between developers' needs and code snippet meta data. Our preliminary evaluation of NLP2Code shows that the majority of invocations lead to code snippets rated as helpful by users and that the tool is able to support a wide range of tasks.


title: On ordinal ranks of Baire class functions

abstract: The theory of ordinal ranks on Baire class 1 functions developed by Kechris and Loveau was recently extended by Elekes, Kiss and Vidny\'{a}nszky to Baire class $\xi$ functions for any countable ordinal $\xi\geq1$. In this paper, we answer two of the questions raised by them in their paper (Ranks on the Baire class $\xi$ functions, Trans. Amer. Math. Soc. 368(2016), 8111-8143). Specifically, we show that for any countable ordinal $\xi\geq1,$ the ranks $\beta_{\xi}^{\ast}$ and $\gamma_{\xi}^{\ast}$ are essentially equivalent, and that neither of them is essentially multiplicative. Since the rank $\beta$ is not essentially multiplicative, we investigate further the behavior of this rank with respect to products. We characterize the functions $f$ so that $\beta(fg)\leq \omega^{\xi}$ whenever $\beta(g)\leq\omega^{\xi}$ for any countable ordinal $\xi.$


title: Demand-Driven Pointer Analysis with Strong Updates via Value-Flow Refinement

abstract: We present a new demand-driven flow- and context-sensitive pointer analysis with strong updates for C programs, called SUPA, that enables computing points-to information via value-flow refinement, in environments with small time and memory budgets such as IDEs. We formulate SUPA by solving a graph reachability problem on an inter-procedural value-flow graph representing a program's def-use chains, which are pre-computed efficiently but over-approximately. To answer a client query (a request for a variable's points-to set), SUPA reasons about the flow of values along the pre-computed def-use chains sparsely (rather than across all program points), by performing only the work necessary for the query (rather than analyzing the whole program). In particular, strong updates are performed to filter out spurious def-use chains through value-flow refinement as long as the total budget is not exhausted. SUPA facilitates efficiency and precision tradeoffs by applying different pointer analyses in a hybrid multi-stage analysis framework. We have implemented SUPA in LLVM (3.5.0) and evaluate it by choosing uninitialized pointer detection as a major client on 18 open-source C programs. As the analysis budget increases, SUPA achieves improved precision, with its single-stage flow-sensitive analysis reaching 97.4% of that achieved by whole-program flow-sensitive analysis by consuming about 0.18 seconds and 65KB of memory per query, on average (with a budget of at most 10000 value-flow edges per query). With context-sensitivity also considered, SUPA's two- stage analysis becomes more precise for some programs but also incurs more analysis times. SUPA is also amenable to parallelization. A parallel implementation of its single-stage flow-sensitive analysis achieves a speedup of up to 6.9x with an average of 3.05x a 8-core machine with respect its sequential version.


title: Nucleon Charges from 2+1+1-flavor HISQ and 2+1-flavor clover lattices

abstract: Precise estimates of the nucleon charges $g_A$, $g_S$ and $g_T$ are needed in many phenomenological analyses of SM and BSM physics. In this talk, we present results from two sets of calculations using clover fermions on 9 ensembles of 2+1+1-flavor HISQ and 4 ensembles of 2+1-flavor clover lattices. We show that high statistics can be obtained cost-effectively using the truncated solver method with bias correction and the coherent source sequential propagator technique. By performing simulations at 4--5 values of the source-sink separation $t_{\rm sep}$, we demonstrate control over excited-state contamination using 2- and 3-state fits. Using the high-precision 2+1+1-flavor data, we perform a simultaneous fit in $a$, $M_\pi$ and $M_\pi L$ to obtain our final results for the charges.


title: Dual Recovery Network with Online Compensation for Image Super-Resolution

abstract: Image super-resolution (SR) methods essentially lead to a loss of some high-frequency (HF) information when predicting high-resolution (HR) images from low-resolution (LR) images without using external references. To address this issue, we additionally utilize online retrieved data to facilitate image SR in a unified deep framework. A novel dual high-frequency recovery network (DHN) is proposed to predict an HR image with three parts: an LR image, an internal inferred HF (IHF) map (HF missing part inferred solely from the LR image) and an external extracted HF (EHF) map. In particular, we infer the HF information based on both the LR image and similar HR references which are retrieved online. For the EHF map, we align the references with affine transformation and then in the aligned references, part of HF signals are extracted by the proposed DHN to compensate for the HF loss. Extensive experimental results demonstrate that our DHN achieves notably better performance than state-of-the-art SR methods.


title: Rigorous Dynamics of Expectation-Propagation-Based Signal Recovery from Unitarily Invariant Measurements

abstract: This paper investigates sparse signal recovery based on expectation propagation (EP) from unitarily invariant measurements. A rigorous analysis is presented for the state evolution (SE) of an EP-based message-passing algorithm in the large system limit, where both input and output dimensions tend to infinity at an identical speed. The main result is the justification of an SE formula conjectured by Ma and Ping.


title: Bayesian Network Learning via Topological Order

abstract: We propose a mixed integer programming (MIP) model and iterative algorithms based on topological orders to solve optimization problems with acyclic constraints on a directed graph. The proposed MIP model has a significantly lower number of constraints compared to popular MIP models based on cycle elimination constraints and triangular inequalities. The proposed iterative algorithms use gradient descent and iterative reordering approaches, respectively, for searching topological orders. A computational experiment is presented for the Gaussian Bayesian network learning problem, an optimization problem minimizing the sum of squared errors of regression models with L1 penalty over a feature network with application of gene network inference in bioinformatics.


title: The Theory of Optical Black Hole Lasers

abstract: The event horizon of black holes and white holes can be achieved in the context of analogue gravity. It was proven for a sonic case that if these two horizons are close to each other their dynamics resemble a laser, a black hole laser, where the analogue of Hawking radiation is trapped and amplified. Optical analogues are also very successful and a similar system can be achieved there. In this work we develop the theory of optical black hole lasers and prove that the amplification is also possible. Then, we study the optical system by determining the forward propagation of modes, obtaining an approximation for the phase difference which governs the amplification, and performing numerical simulations of the pulse propagation of our system.


title: A Two-Step Geometric Framework For Density Modeling

abstract: We introduce a novel two-step approach for estimating a probability density function (pdf) given its samples, with the second and important step coming from a geometric formulation. The procedure involves obtaining an initial estimate of the pdf and then transforming it via a warping function to reach the final estimate. The initial estimate is intended to be computationally fast, albeit suboptimal, but its warping creates a larger, flexible class of density functions, resulting in substantially improved estimation. The search for optimal warping is accomplished by mapping diffeomorphic functions to the tangent space of a Hilbert sphere, a vector space whose elements can be expressed using an orthogonal basis. Using a truncated basis expansion, we estimate the optimal warping under a (penalized) likelihood criterion and, thus, the optimal density estimate. This framework is introduced for univariate, unconditional pdf estimation and then extended to conditional pdf estimation. The approach avoids many of the computational pitfalls associated with classical conditional-density estimation methods, without losing on estimation performance. We derive asymptotic convergence rates of the density estimator and demonstrate this approach using both synthetic datasets and real data, the latter relating to the association of a toxic metabolite on preterm birth.


title: Miura transformations for discrete Painlev\'e equations coming from the affine E$_8$ Weyl group

abstract: We derive integrable equations starting from autonomous mappings with a general form inspired by the additive systems associated to the affine Weyl group E$_8^{(1)}$. By deautonomisation we obtain two hitherto unknown systems, one of which turns out to be a linearisable one, and we show that both these systems arise from the deautonomisation of a non-QRT mapping. In order to unambiguously prove the integrability of these nonautonomous systems, we introduce a series of Miura transformations which allows us to prove that one of these systems is indeed a discrete Painlev\'e equation, related to the affine Weyl group E$_7^{(1)}$, and to cast it in canonical form. A similar sequence of Miura transformations allows us to effectively linearise the second system we obtain. An interesting off-shoot of our calculations is that the series of Miura transformations, when applied at the autonomous limit, allows one to transform a non-QRT invariant into a QRT one.


title: Minimal surfaces in the 3-sphere by desingularizing intersecting Clifford tori

abstract: For each integer $k \geq 2$, we apply gluing methods to construct sequences of minimal surfaces embedded in the round $3$-sphere. We produce two types of sequences, all desingularizing collections of intersecting Clifford tori. Sequences of the first type converge to a collection of $k$ Clifford tori intersecting with maximal symmetry along these two circles. Near each of the circles, after rescaling, the sequences converge smoothly on compact subsets to a Karcher-Scherk tower of order $k$. Sequences of the second type desingularize a collection of the same $k$ Clifford tori supplemented by an additional Clifford torus equidistant from the original two circles of intersection, so that the latter torus orthogonally intersects each of the former $k$ tori along a pair of disjoint orthogonal circles, near which the corresponding rescaled sequences converge to a singly periodic Scherk surface. The simpler examples of the first type resemble surfaces constructed by Choe and Soret \cite{CS} by different methods where the number of handles desingularizing each circle is the same. There is a plethora of new examples which are more complicated and on which the number of handles for the two circles differs. Examples of the second type are new as well.


title: Ordered qausi-two-dimensional structure of nanoparticles in semiflexiblering polymer brushes under compression

abstract: Molecular Dynamics (MD) simulations are presented for a coarse-grained bead-spring model of ring polymer brushes under compression. Flexible polymer brushes are always disordered during compression, whereas semiflexible brushes tend to be ordered under sufficiently strong compression. Besides, the polymer monomer density of semiflexible polymer brush is very high near the polymer brush surface, inducing a peak value of free energy near the polymer brush surface. Therefore, by compressing nanoparticles (NPs) in semiflexible ring brush system, NPs tend to exhibit a closely packed single layer structure between the brush surface and the impenetrable wall, which provide a new access of designing responsive applications.


title: Fragmentation to a jet in the large $z$ limit

abstract: We consider the fragmentation of a parton into a jet with small radius $R$ in the large $z$ limit, where $z$ is the ratio of the jet energy to the mother parton energy. In this region of phase space, large logarithms of both $R$ and $1-z$ can appear, requiring resummation in order to have a well defined perturbative expansion. Using soft-collinear effective theory, we study the fragmentation function to a jet (FFJ) in this endpoint region. We derive a factorization theorem for this object, separating collinear and collinear-soft modes. This allows for the resummation using renormalization group evolution of the logarithms $\ln R$ and $\ln(1-z)$ simultaneously. We show results valid to next-to-leading logarithmic order for the global Sudakov logarithms. We also discuss the possibility of non-global logarithms that should appear at two-loops and give an estimate of their size.


title: Quasi-periodic two-scale homogenisation and effective spatial dispersion in high-contrast media

abstract: The convergence of spectra via two-scale convergence for double-porosity models is well known. A crucial assumption in these works is that the stiff component of the body forms a connected set. We show that under a relaxation of this assumption the (periodic) two-scale limit of the operator is insufficient to capture the full asymptotic spectral properties of high-contrast periodic media. Asymptotically, waves of all periods (or quasi-momenta) are shown to persist and an appropriate extension of the notion of two-scale convergence is introduced. As a result, homogenised limit equations with none trivial quasimomentum dependence are found as resolvent limits of the original operator family, resulting in limiting spectral behaviour with a rich dependence on quasimomenta.


title: Detection and tracking of chemical trails by local sensory systems

abstract: Many aquatic organisms exhibit remarkable abilities to detect and track chemical signals when foraging, mating and escaping. For example, the male copepod { \em T. longicornis} identifies the female in the open ocean by following its chemically-flavored trail. Here, we develop a mathematical framework in which a local sensory system is able to detect the local concentration field and adjust its orientation accordingly. We show that this system is able to detect and track chemical trails without knowing the trail's global or relative position.


title: Nonstandard finite difference schemes for a general predator-prey system

abstract: In this paper we transform a continuous-time predator-prey system with general functional response and recruitment for both species into a discrete-time model by nonstandard finite difference scheme (NSFD). The NSFD model shows complete dynamic consistency with its continuous counterpart for any step size. Especially, the global stability of a non-hyperbolic equilibrium point in a particular case of parameters is proved by the Lyapunov stability theorem. The performed numerical simulations confirmed the validity of the obtained theoretical results.


title: Eccentricity in planetary systems and the role of binarity: sample definition, initial results, and the system of HD 211847

abstract: We explore the multiplicity of exoplanet host stars with high-resolution images obtained with VLT/SPHERE. Two different samples of systems were observed: one containing low-eccentricity outer planets, and the other containing high-eccentricity outer planets. We find that 10 out of 34 stars in the high-eccentricity systems are members of a binary, while the proportion is 3 out of 27 for circular systems. Eccentric-exoplanet hosts are, therefore, significantly more likely to have a stellar companion than circular-exoplanet hosts. The median magnitude contrast over the 68 data sets is 11.26 and 9.25, in H and K, respectively, at 0.30 arcsec. The derived detection limits reveal that binaries with separations of less than 50au are rarer for exoplanet hosts than for field stars. Our results also imply that the majority of high-eccentricity planets are not embedded in multiple stellar systems (24 out of 34), since our detection limits exclude the presence of a stellar companion. We detect the low-mass stellar companions of HD 7449 and HD 211847, both members of our high-eccentricity sample. HD 7449B was already detected by Rodigas et al (2016) and our independent observation is in agreement with this earlier work. HD 211847's substellar companion, previously detected by the radial velocity method, is actually a low-mass star seen face-on. The role of stellar multiplicity in shaping planetary systems is confirmed by this work, although it does not appear as the only source of dynamical excitation.


title: Exponential Stability of Nonlinear Differential Repetitive Processes with Applications to Iterative Learning Control

abstract: This paper studies exponential stability properties of a class of two-dimensional (2D) systems called differential repetitive processes (DRPs). Since a distinguishing feature of DRPs is that the problem domain is bounded in the "time" direction, the notion of stability to be evaluated does not require the nonlinear system defining a DRP to be stable in the typical sense. In particular, we study a notion of exponential stability along the discrete iteration dimension of the 2D dynamics, which requires the boundary data for the differential pass dynamics to converge to zero as the iterations evolve. Our main contribution is to show, under standard regularity assumptions, that exponential stability of a DRP is equivalent to that of its linearized dynamics. In turn, exponential stability of this linearization can be readily verified by a spectral radius condition. The application of this result to Picard iterations and iterative learning control (ILC) is discussed. Theoretical findings are supported by a numerical simulation of an ILC algorithm.


title: A New Family of Error Distributions for Bayesian Quantile Regression

abstract: We propose a new family of error distributions for model-based quantile regression, which is constructed through a structured mixture of normal distributions. The construction enables fixing specific percentiles of the distribution while, at the same time, allowing for varying mode, skewness and tail behavior. It thus overcomes the severe limitation of the asymmetric Laplace distribution -- the most commonly used error model for parametric quantile regression -- for which the skewness of the error density is fully specified when a particular percentile is fixed. We develop a Bayesian formulation for the proposed quantile regression model, including conditional lasso regularized quantile regression based on a hierarchical Laplace prior for the regression coefficients, and a Tobit quantile regression model. Posterior inference is implemented via Markov Chain Monte Carlo methods. The flexibility of the new model relative to the asymmetric Laplace distribution is studied through relevant model properties, and through a simulation experiment to compare the two error distributions in regularized quantile regression. Moreover, model performance in linear quantile regression, regularized quantile regression, and Tobit quantile regression is illustrated with data examples that have been previously considered in the literature.


title: Lyapunov direct method for investigating stability of nonstandard finite difference schemes for metapopulation models

abstract: In this paper nonstandard finite difference (NSFD) schemes of two metapopulation models are constructed. The stability properties of the discrete models are investigated by the use of a generalization of Lyapunov stability theorem. Due to this result we have proved that the NSFD schemes preserve all properties of the metapopulation models. Numerical examples confirm the obtained theoretical results of the properties of the constructed difference schemes. The method of Lyapunov functions proves to be much simpler than the standard method for studying stability of the discrete metapopulation model in our very recent paper.


title: How wrong intuitions about weak* topology and completions of a normed space pose serious problems

abstract: We learn mathematics subjectively and must apply it objectively. But sometimes, we apply it subjectively by using wrong intuitions which may be elusive to our eyes. The aim of this note is to disclose the secretes of two kinds of these false intuitions and the opportunities they may provide. We first discuss the wrong assumption which says that each topology is uniquely determined by studying a very bad phenomenon happening in dealing with weak* topology. Then, we consider the problem of completing a normed space under two comparable norms, one being smaller than the other. Here, we show that the common belief contending that smaller norms give rise to larger completions is wrong. We then pose some serious questions arising from these wrong intuitions. As we will see finding fallacies are as important as major mathematical activities like proving and disproving.


title: Parametric Solution of a Small-Large Black Hole Coexistence Curve

abstract: We consider the first-order phase transition of a charged anti-de Sitter black hole, and find that the equation of state with the conditions of the two coexisting phases, leads to the two coupled equations about the thermodynamic volumes of small black hole and large black hole. By solving the equations, it is found that each reduced volume is only a function of the parameter $\omega$ . All properties of the coexistence curve can be studied from the two volume functions. In particular, each thermodynamic quantity is described by a piecewise analytic function. The demarcation point is located at $\omega_{d}=12(2\sqrt{3}-3)$. The thermodynamic function but not its derivative, is continuous at the point. This property is completely different from that of the ven der Waals fluid. Moreover, the thermodynamic behaviors as $\omega\rightarrow0$ are discussed. From which one can easily obtain some critical exponents and amplitudes for small-large black hole phase transitions.


title: Oscillators that sync and swarm

abstract: Synchronization occurs in many natural and technological systems, from cardiac pacemaker cells to coupled lasers. In the synchronized state, the individual cells or lasers coordinate the timing of their oscillations, but they do not move through space. A complementary form of self-organization occurs among swarming insects, flocking birds, or schooling fish; now the individuals move through space, but without conspicuously altering their internal states. Here we explore systems in which both synchronization and swarming occur together. Specifically, we consider oscillators whose phase dynamics and spatial dynamics are coupled. We call them swarmalators, to highlight their dual character. A case study of a generalized Kuramoto model predicts five collective states as possible long-term modes of organization. These states may be observable in groups of sperm, Japanese tree frogs, colloidal suspensions of magnetic particles, and other biological and physical systems in which self-assembly and synchronization interact.


title: Diversity of chimera-like patterns from a model of 2D arrays of neurons with nonlocal coupling

abstract: Chimera states have been studied in 1D arrays, and a variety of different chimera states have been found using different models. Research has recently been extended to 2D arrays but only to phase models of them. Here, we extend it to a nonphase model of 2D arrays of neurons and focus on the influence of nonlocal coupling. Using extensive numerical simulations, we find, surprisingly, that this system can show most types of previously observed chimera states, in contrast to previous models, where only one or a few types of chimera states can be observed in each model. We also find that this model can show some special chimera-like patterns such as gridding and multicolumn patterns, which were previously observed only in phase models. Further, we present an effective approach, i.e., removing some of the coupling links, to generate heterogeneous coupling, which results in diverse chimera-like patterns and even induces transformations from one chimera-like pattern to another.


title: Stability Enhanced Large-Margin Classifier Selection

abstract: Stability is an important aspect of a classification procedure because unstable predictions can potentially reduce users' trust in a classification system and also harm the reproducibility of scientific conclusions. The major goal of our work is to introduce a novel concept of classification instability, i.e., decision boundary instability (DBI), and incorporate it with the generalization error (GE) as a standard for selecting the most accurate and stable classifier. Specifically, we implement a two-stage algorithm: (i) initially select a subset of classifiers whose estimated GEs are not significantly different from the minimal estimated GE among all the candidate classifiers; (ii) the optimal classifier is chosen as the one achieving the minimal DBI among the subset selected in stage (i). This general selection principle applies to both linear and nonlinear classifiers. Large-margin classifiers are used as a prototypical example to illustrate the above idea. Our selection method is shown to be consistent in the sense that the optimal classifier simultaneously achieves the minimal GE and the minimal DBI. Various simulations and real examples further demonstrate the advantage of our method over several alternative approaches.


title: A modified Newton method for multilinear PageRank

abstract: When studying the multilinear PageRank problem, a system of polynomial equations needs to be solved. In this paper, we develop convergence theory for a modified Newton method in a particular parameter regime. The sequence of vectors produced by Newton-like method is monotonically increasing and converges to the nonnegative solution. Numerical results illustrate the effectiveness of this procedure.


title: Improved Algorithms For Structured Sparse Recovery

abstract: It is known that certain structures of the signal in addition to the standard notion of sparsity (called structured sparsity) can improve the sample complexity in several compressive sensing applications. Recently, Hegde et al. proposed a framework, called approximation-tolerant model-based compressive sensing, for recovering signals with structured sparsity. Their framework requires two oracles, the head- and the tail-approximation projection oracles. The two oracles should return approximate solutions in the model which is closest to the query signal. In this paper, we consider two structured sparsity models and obtain improved projection algorithms. The first one is the tree sparsity model, which captures the support structure in the wavelet decomposition of piecewise-smooth signals. We propose a linear time $(1-\epsilon)$-approximation algorithm for head-approximation projection and a linear time $(1+\epsilon)$-approximation algorithm for tail-approximation projection. The best previous result is an $\tilde{O}(n\log n)$ time bicriterion approximation algorithm (meaning that their algorithm may return a solution of sparsity larger than $k$) by Hegde et al. Our result provides an affirmative answer to the open problem mentioned in the survey of Hegde and Indyk. As a corollary, we can recover a constant approximate $k$-sparse signal. The other is the Constrained Earth Mover Distance (CEMD) model, which is useful to model the situation where the positions of the nonzero coefficients of a signal do not change significantly as a function of spatial (or temporal) locations. We obtain the first single criterion constant factor approximation algorithm for the head-approximation projection. The previous best known algorithm is a bicriterion approximation. Using this result, we can get a faster constant approximation algorithm with fewer measurements for the recovery problem in CEMD model.


title: Analysis of Scattering from Archival Pulsar Data using a CLEAN-based Method

abstract: In this work we adopted a CLEAN-based method to determine the scatter time, $\tau$, from archived pulsar profiles under both the thin screen and uniform medium scattering models and to calculate the scatter time frequency scale index $\alpha$ where $\tau\propto\nu^{\alpha}$. The value of $\alpha$ is $-4.4$, if a Kolmogorov spectrum of the interstellar medium turbulence is assumed. We deconvolved 1342 profiles from 347 pulsars over a broad range of frequencies and dispersion measures. In our survey, in the majority of cases the scattering effect was not significant compared to pulse profile widths. For a subset of 21 pulsars scattering at the lowest frequencies was large enough to be measured. Because reliable scatter time measurements were determined only for the lowest frequency, we were limited to using upper limits on scatter times at higher frequencies for the purpose of our scatter time frequency slope estimation. We scaled the deconvolved scatter time to 1~GHz assuming $\alpha=-4.4$ and considered our results in the context of other observations which yielded a broad relation between scatter time and dispersion measure.


title: Efficient Feature Matching by Progressive Candidate Search

abstract: We present a novel feature matching algorithm that systematically utilizes the geometric properties of features such as position, scale, and orientation, in addition to the conventional descriptor vectors. In challenging scenes with the presence of repetitive patterns or with a large viewpoint change, it is hard to find the correct correspondences using feature descriptors only, since the descriptor distances of the correct matches may not be the least among the candidates due to appearance changes. Assuming that the layout of the nearby features does not changed much, we propose the bidirectional transfer measure to gauge the geometric consistency of a pair of feature correspondences. The feature matching problem is formulated as a Markov random field (MRF) which uses descriptor distances and relative geometric similarities together. The unmatched features are explicitly modeled in the MRF to minimize its negative impact. For speed and stability, instead of solving the MRF on the entire features at once, we start with a small set of confident feature matches, and then progressively search the candidates in nearby features and expand the MRF with them. Experimental comparisons show that the proposed algorithm finds better feature correspondences, i.e. more matches with higher inlier ratio, in many challenging scenes with much lower computational cost than the state-of-the-art algorithms.


title: RVP-FLMS : A Robust Variable Power Fractional LMS Algorithm

abstract: In this paper, we propose an adaptive framework for the variable power of the fractional least mean square (FLMS) algorithm. The proposed algorithm named as robust variable power FLMS (RVP-FLMS) dynamically adapts the fractional power of the FLMS to achieve high convergence rate with low steady state error. For the evaluation purpose, the problems of system identification and channel equalization are considered. The experiments clearly show that the proposed approach achieves better convergence rate and lower steady-state error compared to the FLMS. The MATLAB code for the related simulation is available online at https://goo.gl/dGTGmP.


title: Study of radially excited $D_s(2^1S_0)$ and $D^+_{sJ}(2632)$

abstract: The $J^P=0^-$ radial excitation $D_s(2^1S_0)$ is anticipated to have mass $2650$ MeV (denoted with $D_s(2650)$) though it has not been observed. $D_s(2650)$ is anticipated to be observed in inclusive $e^+e^-$ and $pp$ collisions in $D^*K$ channel. $D_s(2650)$ can be produced from hadronic decays of higher excited resonances. In a $^3P_0$ model, hadronic production of $D_s(2650)$ has been studied, relevant decay widths have been estimated. The hadronic decays of $D_s(3P)$ have been explored in detail, and the dominant decay channels have been pointed out. Two hadronic decay channels of $D^+_{sJ}(2632)$ were observed by SELEX, but have not been observed by any other experiment. Hadronic decays of $D^+_{sJ}(2632)$ in different assignments have been explored in the $^3P_0$ model. In view of the hadronic decay manner of $D^+_{sJ}(2632)$ , we conclude that $D^+_{sJ}(2632)$ cannot be a conventional charmed strange $c\bar s$ meson.


title: A Computational Study of Yttria-Stabilized Zirconia: II. Cation Diffusion

abstract: Cubic yttria-stabilized zirconia is widely used in industrial electrochemical devices. While its fast oxygen ion diffusion is well understood, why cation diffusion is much slower-its activation energy (~5 eV) is 10 times that of anion diffusion-remains a mystery. Indeed, all previous computational studies predicted more than 5 eV is needed for forming a cation defect, and another 5 eV for moving one. In contrast, our ab initio calculations have correctly predicted the experimentally observed cation diffusivity. We found Schottky pairs are the dominant defects that provide cation vacancies, and their local environments and migrating path are dictated by packing preferences. As a cation exchanges position with a neighboring vacancy, it passes by an empty interstitial site and severely displaces two oxygen neighbors with shortened Zr-O distances. This causes a short-range repulsion against the migrating cation and a long-range disturbance of the surrounding, which explains why cation diffusion is relatively difficult. In comparison, cubic zirconia's migrating oxygen only minimally disturbs neighboring Zr, which explains why it is a fast oxygen conductor.


title: Strong Convergence Rate of Splitting Schemes for Stochastic Nonlinear Schr\"odinger Equations

abstract: We prove the optimal strong convergence rate of a fully discrete scheme, based on a splitting approach, for a stochastic nonlinear Schr\"odinger (NLS) equation. The main novelty of our method lies on the uniform a priori estimate and exponential integrability of a sequence of splitting processes which are used to approximate the solution of the stochastic NLS equation. We show that the splitting processes converge to the solution with strong order $1/2$. Then we use the Crank--Nicolson scheme to temporally discretize the splitting process and get the temporal splitting scheme which also possesses strong order $1/2$. To obtain a full discretization, we apply this splitting Crank--Nicolson scheme to the spatially discrete equation which is achieved through the spectral Galerkin approximation. Furthermore, we establish the convergence of this fully discrete scheme with optimal strong convergence rate $\mathcal{O}(N^{-2}+\tau^\frac12)$, where $N$ denotes the dimension of the approximate space and $\tau$ denotes the time step size. To the best of our knowledge, this is the first result about strong convergence rates of temporally numerical approximations and fully discrete schemes for stochastic NLS equations, or even for stochastic partial differential equations (SPDEs) with non-monotone coefficients. Numerical experiments verify our theoretical result.


title: Git Blame Who?: Stylistic Authorship Attribution of Small, Incomplete Source Code Fragments

abstract: Program authorship attribution has implications for the privacy of programmers who wish to contribute code anonymously. While previous work has shown that complete files that are individually authored can be attributed, we show here for the first time that accounts belonging to open source contributors containing short, incomplete, and typically uncompilable fragments can also be effectively attributed. We propose a technique for authorship attribution of contributor accounts containing small source code samples, such as those that can be obtained from version control systems or other direct comparison of sequential versions. We show that while application of previous methods to individual small source code samples yields an accuracy of about 73% for 106 programmers as a baseline, by ensembling and averaging the classification probabilities of a sufficiently large set of samples belonging to the same author we achieve 99% accuracy for assigning the set of samples to the correct author. Through these results, we demonstrate that attribution is an important threat to privacy for programmers even in real-world collaborative environments such as GitHub. Additionally, we propose the use of calibration curves to identify samples by unknown and previously unencountered authors in the open world setting. We show that we can also use these calibration curves in the case that we do not have linking information and thus are forced to classify individual samples directly. This is because the calibration curves allow us to identify which samples are more likely to have been correctly attributed. Using such a curve can help an analyst choose a cut-off point which will prevent most misclassifications, at the cost of causing the rejection of some of the more dubious correct attributions.


title: Laser opacity in underdense preplasma of solid targets due to quantum electrodynamics effects

abstract: We investigate how next-generation laser pulses at 10 PW $-$ 200 PW interact with a solid target in the presence of a relativistically underdense preplasma produced by amplified spontaneous emission (ASE). Laser hole boring and relativistic transparency are strongly restrained due to the generation of electron-positron pairs and $\gamma$-ray photons via quantum electrodynamics (QED) processes. A pair plasma with a density above the initial preplasma density is formed, counteracting the electron-free channel produced by the hole boring. This pair-dominated plasma can block the laser transport and trigger an avalanche-like QED cascade, efficiently transfering the laser energy to photons. This renders a 1-$\rm\mu m$-scalelength, underdense preplasma completely opaque to laser pulses at this power level. The QED-induced opacity therefore sets much higher contrast requirements for such pulse in solid-target experiments than expected by classical plasma physics. Our simulations show for example, that proton acceleration from the rear of a solid with a preplasma would be strongly impaired.


title: Black hole mass measurement using molecular gas kinematics: what ALMA can do

abstract: We study the limits of the spatial and velocity resolution of radio interferometry to infer the mass of supermassive black holes (SMBHs) in galactic centres using the kinematics of circum-nuclear molecular gas, by considering the shapes of the galaxy surface brightness profile, signal-to-noise ratios (S/Ns) of the position-velocity diagram (PVD) and systematic errors due to the spatial and velocity structure of the molecular gas. We argue that for fixed galaxy stellar mass and SMBH mass, the spatial and velocity scale that need to be resolved increase and decrease, respectively, with decreasing \sersic\ index of the galaxy surface brightness profile. We validate our arguments using simulated PVDs for varying beam size and velocity channel width. Furthermore, we consider the systematic effects on the inference of the SMBH mass by simulating PVDs including the spatial and velocity structure of the molecular gas, which demonstrates that their impacts are not significant for a PVD with good S/N unless the spatial and velocity scale associated with the systematic effects are comparable to or larger than the angular resolution and velocity channel width of the PVD from pure circular motion. Also, we caution that a bias in a galaxy surface brightness profile owing to the poor resolution of a galaxy photometric image can largely bias the SMBH mass by an order of magnitude. This study shows the promise and the limit of ALMA observations for measuring SMBH mass using molecular gas kinematics and provides a useful technical justification for an ALMA proposal with the science goal of measuring SMBH mass.


title: Two-magnon excitations in resonant inelastic x-ray scattering studied by spin-density-wave formalism

abstract: We study two-magnon excitations in resonant inelastic x-ray scattering (RIXS) at the transition-metal $K$-edge. Instead of working with effective Heisenberg spin models, we work with a Hubbard-type model ($d$-$p$ model) for a typical insulating cuprate La$_2$CuO$_4$. For the antiferromagnetic ground state within the spin-density-wave (SDW) mean-field formalism, we calculate the dynamical correlation function within the random-phase approximation (RPA), and then obtain two-magnon excitation spectra by calculating the convolution of it. Coupling between the $K$-shell hole and the magnons in the intermediate state is calculated by means of diagrammatic perturbation expansion in the Coulomb interaction. Calculated momentum dependence of RIXS spectra agrees well with that of experiments. A notable difference from previous calculations based on the Heisenberg spin models is that RIXS spectra have a large two-magnon weight near the zone center, which may be confirmed by further careful high-resolution experiments.


title: Quantitative modeling and analysis of bifurcation-induced bursting

abstract: Modeling and parameter estimation for neuronal dynamics are often challenging because many parameters can range over orders of magnitude and are difficult to measure experimentally. Moreover, selecting a suitable model complexity requires a sufficient understanding of the model's potential use, such as highlighting essential mechanisms underlying qualitative behavior or precisely quantifying realistic dynamics. We present a novel approach that can guide model development and tuning to achieve desired qualitative and quantitative solution properties. Our approach relies on the presence of disparate time scales and employs techniques of separating the dynamics of fast and slow variables, which are well known in the analysis of qualitative solution features. We build on these methods to show how it is also possible to obtain quantitative solution features by imposing designed dynamics for the slow variables in the form of specified two-dimensional paths in a bifurcation-parameter landscape.


title: High Rate LDPC Codes from Difference Covering Arrays

abstract: This paper presents a combinatorial construction of low-density parity-check (LDPC) codes from difference covering arrays. While the original construction by Gallagher was by randomly allocating bits in a sparse parity-check matrix, over the past 20 years researchers have used a variety of more structured approaches to construct these codes, with the more recent constructions of well-structured LDPC coming from balanced incomplete block designs (BIBDs) and from Latin squares over finite fields. However these constructions have suffered from the limited orders for which these designs exist. Here we present a construction of LDPC codes of length $4n^2 - 2n$ for all $n$ using the cyclic group of order $2n$. These codes achieve high information rate (greater than 0.8) for $n \geq 8$, have girth at least 6 and have minimum distance 6 for $n$ odd.


title: Scalar extensions of categorical resolutions of singularities

abstract: Let $X$ be a quasi-compact, separated scheme over a field k and we can consider the categorical resolution of singularities of $X$. In this paper let $k^{\prime}/k$ be a field extension and we study the scalar extension of a categorical resolution of singularities of $X$ and we show how it gives a categorical resolution of the base change scheme $X_{k^{\prime}}$. Our construction involves the scalar extension of derived categories of DG-modules over a DG algebra. As an application we use the technique of scalar extension developed in this paper to prove the non-existence of full exceptional collections of categorical resolutions for a projective curve of genus $\geq 1$ over a non-algebraically closed field.


title: Verified measurement-based quantum computing with hypergraph states

abstract: Hypergraph states are generalizations of graph states where controlled-$Z$ gates on edges are replaced with generalized controlled-$Z$ gates on hyperedges. Hypergraph states have several advantages over graph states. For example, certain hypergraph states, such as the Union Jack states, are universal resource states for measurement-based quantum computing with only Pauli measurements, while graph state measurement-based quantum computing needs non-Clifford basis measurements. Furthermore, it is impossible to classically efficiently sample measurement results on hypergraph states with a constant $L$1-norm error unless the polynomial hierarchy collapses to the third level. Although several protocols have been proposed to verify graph states with only sequential single-qubit Pauli measurements, there was no verification method for hypergraph states. In this paper, we propose a method for verifying hypergraph states with only sequential single-qubit Pauli measurements. As applications, we consider verified blind quantum computing with hypergraph states, and quantum supremacy demonstrations with hypergraph states.


title: Photometric classification and redshift estimation of LSST Supernovae

abstract: Supernova (SN) classification and redshift estimation using photometric data only have become very important for the Large Synoptic Survey Telescope (LSST), given the large number of SNe that LSST will observe and the impossibility of spectroscopically following up all the SNe. We investigate the performance of a SN classifier that uses SN colors to classify LSST SNe with the Random Forest classification algorithm. Our classifier results in an AUC of 0.98 which represents excellent classification. We are able to obtain a photometric SN sample containing 99$\%$ SNe Ia by choosing a probability threshold. We estimate the photometric redshifts (photo-z) of SNe in our sample by fitting the SN light curves using the SALT2 model with nested sampling. We obtain a mean bias ($\left<z_\mathrm{phot}-z_\mathrm{spec}\right>$) of 0.012 with $\sigma\left( \frac{z_\mathrm{phot}-z_\mathrm{spec}}{1+z_\mathrm{spec}}\right) = 0.0294$ without using a host-galaxy photo-z prior, and a mean bias ($\left<z_\mathrm{phot}-z_\mathrm{spec}\right>$) of 0.0017 with $\sigma\left( \frac{z_\mathrm{phot}-z_\mathrm{spec}}{1+z_\mathrm{spec}}\right) = 0.0116$ using a host-galaxy photo-z prior. Assuming a flat $\Lambda CDM$ model with $\Omega_m=0.3$, we obtain $\Omega_m$ of $0.305\pm0.008$ (statistical errors only), using the simulated LSST sample of photometric SNe Ia (with intrinsic scatter $\sigma_\mathrm{int}=0.11$) derived using our methodology without using host-galaxy photo-z prior. Our method will help boost the power of SNe from the LSST as cosmological probes.


title: Dark matter statistics for large galaxy catalogs: power spectra and covariance matrices

abstract: Upcoming and existing large-scale surveys of galaxies require accurate theoretical predictions of the dark matter clustering statistics for thousands of mock galaxy catalogs. We demonstrate that this goal can be achieve with our new Parallel Particle-Mesh (PM) Nbody code (PPM-GLAM) at a very low computational cost. We run about 15,000 simulations with ~2 billion particles that provide ~1% accuracy of the dark matter power spectra P(k) for wave-numbers up to k~ 1h/Mpc. Using this large data-set we study the power spectrum covariance matrix, the stepping stone for producing mock catalogs. In contrast to many previous analytical and numerical results, we find that the covariance matrix normalised to the power spectrum C(k,k')/P(k)P(k') has a complex structure of non-diagonal components. It has an upturn at small k, followed by a minimum at k=0.1-0.2h/Mpc. It also has a maximum at k=0.5-0.6h/Mpc. The normalised covariance matrix strongly evolves with redshift: C(k,k')~delta(t)^alpha P(k)P(k'), where delta is the linear growth factor and alpha ~ 1-1.25, which indicates that the covariance matrix depends on cosmological parameters. We also show that waves longer than 1Gpc have very little impact on the power spectrum and covariance matrix. This significantly reduces the computational costs and complexity of theoretical predictions: relatively small volume ~ (1Gpc)^3 simulations capture the necessary properties of dark matter clustering statistics. All the power spectra obtained from many thousands of our simulations are publicly available.


title: Real-time Traffic Accident Risk Prediction based on Frequent Pattern Tree

abstract: Traffic accident data are usually noisy, contain missing values, and heterogeneous. How to select the most important variables to improve real-time traffic accident risk prediction has become a concern of many recent studies. This paper proposes a novel variable selection method based on the Frequent Pattern tree (FP tree) algorithm. First, all the frequent patterns in the traffic accident dataset are discovered. Then for each frequent pattern, a new criterion, called the Relative Object Purity Ratio (ROPR) which we proposed, is calculated. This ROPR is added to the importance score of the variables that differentiate one frequent pattern from the others. To test the proposed method, a dataset was compiled from the traffic accidents records detected by only one detector on interstate highway I-64 in Virginia in 2005. This dataset was then linked to other variables such as real-time traffic information and weather conditions. Both the proposed method based on the FP tree algorithm, as well as the widely utilized, random forest method, were then used to identify the important variables or the Virginia dataset. The results indicate that there are some differences between the variables deemed important by the FP tree and those selected by the random forest method. Following this, two baseline models (i.e. a nearest neighbor (k-NN) method and a Bayesian network) were developed to predict accident risk based on the variables identified by both the FP tree method and the random forest method. The results show that the models based on the variable selection using the FP tree performed better than those based on the random forest method for several versions of the k-NN and Bayesian network models.The best results were derived from a Bayesian network model using variables from FP tree. That model could predict 61.11% of accidents accurately while having a false alarm rate of 38.16%.


title: Finiteness theorems for holomorphic mappings from products of hyperbolic Riemann surfaces

abstract: We prove that the space of dominant/non-constant holomorphic mappings from a product of hyperbolic Riemann surfaces of finite type into certain hyperbolic manifolds with universal cover a bounded domain is a finite set.


title: Gravitational Lensing of Rays through the Levitating Atmospheres of Compact Objects

abstract: Electromagnetic rays travel on curved paths under the influence of gravity. When a dispersive optical medium is included, these trajectories are frequency-dependent. In this work we consider the behaviour of rays when a spherically symmetric, luminous compact object described by the Schwarzschild metric is surrounded by an optically thin shell of plasma supported by radiation pressure. Such levitating atmospheres occupy a position of stable radial equilibrium, where radiative flux and gravitational effects are balanced. Using general relativity and an inhomogeneous plasma we find the existence of a stable circular orbit within the atmospheric shell for low-frequency rays. We explore families of bound orbits that exist between the shell and the compact object, and identify sets of novel periodic orbits. Finally, we examine conditions necessary for the trapping and escape of low-frequency radiation.


title: Efficient and accurate numerical schemes for a hydrodynamically coupled phase field diblock copolymer model

abstract: In this paper, we consider numerical approximations of a hydrodynamically coupled phase field diblock copolymer model, in which the free energy contains a kinetic potential, a gradient entropy, a Ginzburg-Landau double well potential, and a long range nonlocal type potential. We develop a set of second order time marching schemes for this system using the "Invariant Energy Quadratization" approach for the double well potential, the projection method for the Navier-Stokes equation, and a subtle implicit-explicit treatment for the stress and convective term. The resulting schemes are linear and lead to symmetric positive definite systems at each time step, thus they can be efficiently solved. We further prove that these schemes are unconditionally energy stable. Various numerical experiments are performed to validate the accuracy and energy stability of the proposed schemes.


title: The Value of Timing Risk

abstract: The aim of this paper is to provide a mathematical contribution on the semi-static hedge of timing risk associated to positions in American-style options under a multi-dimensional market model. Barrier options are considered in the paper and semi-static hedges are studied and discussed for a fairly large class of underlying price dynamics. Timing risk is identified with the uncertainty associated to the time at which the payoff payment of the barrier option is due. Starting from the work by Carr and Picron (1999), where the authors show that the timing risk can be hedged via static positions in plain vanilla options, the present paper extends the static hedge formula proposed in Carr and Picron (1999) by giving sufficient conditions to decompose a generalized timing risk into an integral of knock-in options in a multi-dimensional market model. A dedicated study of the semi-static hedge is then conducted by defining the corresponding strategy based on positions in barrier options. The proposed methodology allows to construct not only first order hedges but also higher order semi-static hedges, that can be interpreted as asymptotic expansions of the hedging error. The convergence of these higher order semi-static hedges to an exact hedge is shown. An illustration of the main theoretical results is provided for i) a symmetric case, ii) a one dimensional case, where the first order and second order hedging errors are derived in analytic closed form. The materiality of the hedging benefit gain of going from order one to order two by re-iterating the timing risk hedging strategy is discussed through numerical evidences by showing that order two can bring to more than 90% reduction of the hedging 'cost' w.r.t. order one (depending on the specific barrier option characteristics).


title: Two-step pulse observation to improve resonance contrast for coherent population trapping atomic clock

abstract: We study resonance contrast by a two-step pulse observation method to enhance the frequency stability of coherent population trapping (CPT) atomic clocks. The proposed method is a two-step Raman--Ramsey scheme with low intensity during resonance observation and high intensity after the observation. This method reduces the frequency variation in the light intensity and maintains a high signal-to-noise ratio. The resonance characteristics were calculated by density matrix analysis of a $\Lambda$-type three-level system that was modeled on the $^{133}$Cs D1 line, and the characteristics were also measured using a vertical-cavity surface-emitting laser and a Cs vapor cell.


title: High order approximation for the Boltzmann equation without angular cutoff

abstract: In order to solve the Boltzmann equation numerically, in the present work, we propose a new model equation to approximate the Boltzmann equation without angular cutoff. Here the approximate equation incorporates Boltzmann collision operator with angular cut-off and the Landau collision operator. As a first step, we prove the well-posedness theory for our approximate equation. Then in the next step we show the error estimate between the solutions to the approximate equation and the original equation. Compared to the standard angular cut-off approximation method, our method results in higher order of accuracy.


title: The AT-LESS CO(1-0) survey of submillimetre galaxies in the Extended Chandra Deep Field South: First results on cold molecular gas in galaxies at z ~ 2

abstract: We present the first results from our on-going Australia Telescope Compact Array survey of CO(1-0) in ALMA-identified submillimetre galaxies in the Extended Chandra Deep Field South. Strong detections of CO(1-0) emission from two submillimetre galaxies, ALESS 122.1 (z = 2.0232) and ALESS 67.1 (z = 2.1230), were obtained. We estimate gas masses of M_gas ~ 1.3 \times 10^{11} M_odot and M_gas ~ 1.0 \times 10^{11} M_\odot for ALESS 122.1 and ALESS 67.1, respectively, adopting alpha_CO = 1.0. Dynamical mass estimates from the kinematics of the CO(1-0) line yields M_dyn (sin i)^2 = 2.1 +- 1.1 \times 10^{11} M_odot and (3.2 +- 0.9) \times 10^{11} M_\odot for ALESS 122.1 and ALESS 67.1, respectively. This is consistent with the total baryonic mass estimates of these two systems. We examine star formation efficiency using the L_FIR versus L'_CO(1-0) relation for samples of local ULIRGs and LIRGs, and more distant star-forming galaxies, with CO(1-0) detections. We find some evidence of a shallower slope for ULIRGs and SMGs compared to less luminous systems, but a larger sample is required for definite conclusions. We determine gas-to-dust ratios of 170 +- 30 and 140 +- 30 for ALESS 122.1 and ALESS 67.1, respectively, showing ALESS 122.1 has an unusually large gas reservoir. By combining the 38.1 GHz continuum detection of ALESS 122.1 with 1.4 and 5.5 GHz data, we estimate that the free-free contribution to radio emission at 38.1 GHz is 34 +- 17 microJy, yielding a star formation rate (1400 +- 700 M_\odot yr^{-1}) consistent with that from the infrared luminosity.


title: Efficiently Computing Provenance Graphs for Queries with Negation

abstract: Explaining why an answer is in the result of a query or why it is missing from the result is important for many applications including auditing, debugging data and queries, and answering hypothetical questions about data. Both types of questions, i.e., why and why-not provenance, have been studied extensively. In this work, we present the first practical approach for answering such questions for queries with negation (first-order queries). Our approach is based on a rewriting of Datalog rules (called firing rules) that captures successful rule derivations within the context of a Datalog query. We extend this rewriting to support negation and to capture failed derivations that explain missing answers. Given a (why or why-not) provenance question, we compute an explanation, i.e., the part of the provenance that is relevant to answer the question. We introduce optimizations that prune parts of a provenance graph early on if we can determine that they will not be part of the explanation for a given question. We present an implementation that runs on top of a relational database using SQL to compute explanations. Our experiments demonstrate that our approach scales to large instances and significantly outperforms an earlier approach which instantiates the full provenance to compute explanations.


title: Antenna Deployment Method for MIMO Radar under the Situation of Multiple Interference Regions

abstract: In this paper, considering multiple interference regions simultaneously, an optimal antenna deployment problem for distributed Multi-Input Multi-Output (MIMO) radar is investigated. The optimal antenna deployment problem is solved by proposing an antenna deployment method based on Multi-Objective Particle Swarm Optimization (MOPSO). Firstly, we construct a multi-objective optimization problem for MIMO radar antenna deployment by choosing the interference power densities of different regions as objective functions. Then, to obtain the optimal deployment result without wasting time and computational resources, an iteration convergence criterion based on interval distance is proposed. The iteration convergence criterion can be used to stop the MOPSO optimization process efficiently when the optimal antenna deployment algorithm reaches the desired convergence level. Finally, numerical results are provided to verify the validity of the proposed algorithm.


title: Analytic evaluation of some three- and four- electron atomic integrals involving s STO's and exponential correlation with unlinked $r_{ij}$'s

abstract: The method of evaluation outlined in a previous work has been utilized here to evaluate certain other three- electron and four- electron atomic integrals involving s Slater-type orbitals and exponential correlation with unlinked $r_{ij}$'s. Limiting expressions for various such integrals have been derived, which has not been done earlier. Closed-form expressions for $<r_{12} r_{13} / r_{14}>$, $<r_{12}r_{34}/r_{23}>$, $<r_{12}r_{23}/r_{34}>$, $<r_{12}r_{13}/r_{34}>$ and $<r_{12}r_{34}/r_{13}>$ have been obtained.


title: A review of the past and present MAGIC dark matter search program and a glimpse at the future

abstract: The MAGIC TeV gamma-ray telescopes have devoted several hundreds hour of observation time in about a decade, to hunt for particle dark matter indirect signatures in gamma rays, from various candidate targets of interest in the sky: the galactic center, satellite galaxies, galaxy clusters and unidentified objects in other bands. Despite the effort, no hints are present in MAGIC data. These observation are nevertheless not unusable. MAGIC indeed derived the most robust upper limits in the TeV range than any other instrument. These results, for the time being, only mildly constrain some classic dark matter models, but are of use in the construction of dark matter models for the next searches, that consider also the negative results from accelerator and direct-detection experiments. In the contribution, we discuss and review MAGIC results, putting them into context, and in perspective with the next generation of ground-based Cherenkov telescopes. We will briefly inform about future MAGIC projects regarding dark matter searches.


title: Automatic Generation of Typographic Font from a Small Font Subset

abstract: This paper addresses the automatic generation of a typographic font from a subset of characters. Specifically, we use a subset of a typographic font to extrapolate additional characters. Consequently, we obtain a complete font containing a number of characters sufficient for daily use. The automated generation of Japanese fonts is in high demand because a Japanese font requires over 1,000 characters. Unfortunately, professional typographers create most fonts, resulting in significant financial and time investments for font generation. The proposed method can be a great aid for font creation because designers do not need to create the majority of the characters for a new font. The proposed method uses strokes from given samples for font generation. The strokes, from which we construct characters, are extracted by exploiting a character skeleton dataset. This study makes three main contributions: a novel method of extracting strokes from characters, which is applicable to both standard fonts and their variations; a fully automated approach for constructing characters; and a selection method for sample characters. We demonstrate our proposed method by generating 2,965 characters in 47 fonts. Objective and subjective evaluations verify that the generated characters are similar to handmade characters.


title: Some functional inequalities on non-reversible Finsler manifolds

abstract: We continue our study of geometric analysis on (possibly non-reversible) Finsler manifolds, based on the Bochner inequality established by the author and Sturm. Following the approach of the $\Gamma$-calculus a la Bakry et al, we show the dimensional versions of the Poincare--Lichnerowicz inequality, the logarithmic Sobolev inequality, and the Sobolev inequality. In the reversible case, these inequalities were obtained by Cavalletti--Mondino in the framework of curvature-dimension condition by means of the localization method. We show that the same (sharp) estimates hold also for non-reversible metrics.


title: Construction and nonexistence of strong external difference families

abstract: Strong external difference families (SEDFs) were introduced by Paterson and Stinson as a more restrictive version of external difference families. SEDFs can be used to produce optimal strong algebraic manipulation detection codes. We characterize the parameters $(v, m, k, \lambda)$ of a nontrivial SEDF that is near-complete (satisfying $v=km+1$). We construct the first known nontrivial example of a $(v, m, k, \lambda)$ SEDF having $m > 2$. The parameters of this example are $(243,11,22,20)$, giving a near-complete SEDF, and its group is $\mathbb{Z}_3^5$. We provide a comprehensive framework for the study of SEDFs using character theory and algebraic number theory, showing that the cases $m=2$ and $m>2$ are fundamentally different. We prove a range of nonexistence results, greatly narrowing the scope of possible parameters of SEDFs.


title: Discrete Spectrum Reconstruction using Integral Approximation Algorithm

abstract: An inverse problem in spectroscopy is considered. The objective is to restore the discrete spectrum from observed spectrum data, taking into account the spectrometer's line spread function. The problem is reduced to solution of a system of linear-nonlinear equations (SLNE) with respect to intensities and frequencies of the discrete spectral lines. The SLNE is linear with respect to lines' intensities and nonlinear with respect to the lines' frequencies. The integral approximation algorithm is proposed for the solution of this SLNE. The algorithm combines solution of linear integral equations with solution of a system of linear algebraic equations and avoids nonlinear equations. Numerical examples of the application of the technique, both to synthetic and experimental spectra, demonstrate the efficacy of the proposed approach in enabling an effective enhancement of the spectrometer's resolution.


title: Dynamics of Dirac solitons in networks

abstract: We study dynamics of Dirac solitons in prototypical networks modeling them by the nonlinear Dirac equation on metric graphs. Soliton solutions of the nonlinear Dirac equation on simple metric graphs are obtained. It is shown that these solutions provide reflectionless vertex transmission of the Dirac solitons under suitable conditions. The constraints for bond nonlinearity coefficients, allowing reflectionless transmission over a Y-junction are derived. The analytical results are confirmed by direct numerical simulations.

categories: nlin.PS cond-mat.mes-hall hep-th quant-ph


title: The regularization theory of the Krylov iterative solvers LSQR and CGLS for linear discrete ill-posed problems, part I: the simple singular value case

abstract: For the large-scale linear discrete ill-posed problem $\min\|Ax-b\|$ or $Ax=b$ with $b$ contaminated by a white noise, the Lanczos bidiagonalization based LSQR method and its mathematically equivalent Conjugate Gradient (CG) method for $A^TAx=A^Tb$ are most commonly used. They have intrinsic regularizing effects, where the number $k$ of iterations plays the role of regularization parameter. However, there has been no answer to the long-standing fundamental concern by Bj\"{o}rck and Eld\'{e}n in 1979: for which kinds of problems LSQR and CGLS can find best possible regularized solutions? Here a best possible regularized solution means that it is at least as accurate as the best regularized solution obtained by the truncated singular value decomposition (TSVD) method or standard-form Tikhonov regularization. In this paper, assuming that the singular values of $A$ are simple, we analyze the regularization of LSQR for severely, moderately and mildly ill-posed problems. We establish accurate estimates for the 2-norm distance between the underlying $k$-dimensional Krylov subspace and the $k$-dimensional dominant right singular subspace of $A$. For the first two kinds of problems, we then prove that LSQR finds a best possible regularized solution at semi-convergence occurring at iteration $k_0$ and that, for $k=1,2,\ldots,k_0$, (i) the $k$-step Lanczos bidiagonalization always generates a near best rank $k$ approximation to $A$; (ii) the $k$ Ritz values always approximate the first $k$ large singular values in natural order; (iii) the $k$-step LSQR always captures the $k$ dominant SVD components of $A$. For the third kind of problem, we prove that LSQR generally cannot find a best possible regularized solution. Numerical experiments confirm our theory.


title: A Unified Stochastic Formulation of Dissipative Quantum Dynamics. I. Generalized Hierarchical Equations

abstract: We extend a standard stochastic theory to study open quantum systems coupled to generic quantum environments including the three fundamental classes of noninteracting particles: bosons, fermions and spins. In this unified stochastic approach, the generalized stochastic Liouville equation (SLE) formally captures the exact quantum dissipations when noise variables with appropriate statistics for different bath models are applied. Anharmonic effects of a non-Gaussian bath are precisely encoded in the bath multi-time correlation functions that noise variables have to satisfy. Staring from the SLE, we devise a family of generalized hierarchical equations by averaging out the noise variables and expand bath multi-time correlation functions in a complete basis of orthonormal functions. The general hiearchical equations constitute systems of linear equations that provide numerically exact simulations of quantum dynamics. For bosonic bath models, our general hierarchical equation of motion reduces exactly to an extended version of hierarchical equation of motion which allows efficient simulation for arbitrary spectral densities and temperature regimes. Similar efficiency and exibility can be achieved for the fermionic bath models within our formalism. The spin bath models can be simulated with two complementary approaches in the presetn formalism. (I) They can be viewed as an example of non-Gaussian bath models and be directly handled with the general hierarchical equation approach given their multi-time correlation functions. (II) Alterantively, each bath spin can be first mapped onto a pair of fermions and be treated as fermionic environments within the present formalism.


title: High-Contrast Chirped-Pulse Amplification Enabled by In-Band Noise Filtering

abstract: Lasers that generate ultra-intense light pulses are under development for experiments in high-field and high-energy-density physics, as well as for applications such as particle acceleration. Extensions to even higher powers are being considered for future investigations that can only be imagined today, such as the quantum electrodynamics of plasmas and isolated attosecond-pulse generation with solid targets. For all of these areas, it is vital to produce high-contrast pulses, so that no pre-plasma is created in the target before the arrival of the main pulse. However, noise is unavoidable in high-gain amplification, and is manifested in the form of background light that accompanies pulses generated by chirped-pulse amplification (CPA). Here, we introduce a linear filtering technique based on spatio-spectral coupling, which allows in-band filtering of amplified pulses for the first time. Experiments demonstrate approximately 40 times contrast enhancement in optical parametric chirped-pulse amplification (OPCPA) and provide a foundation for scaling to much higher performance. The simplicity, efficiency, and direct compatibility with existing techniques for short-pulse generation will make spatio-spectral filtering attractive to a wide range of applications in ultrafast optics and time-resolved spectroscopy, and may open new directions in noise reduction.


title: Age-Optimal Information Updates in Multihop Networks

abstract: The problem of reducing the age-of-information has been extensively studied in the single-hop networks. In this paper, we minimize the age-of-information in general multihop networks. If the packet transmission times over the network links are exponentially distributed, we prove that a preemptive Last Generated First Served (LGFS) policy results in smaller age processes at all nodes of the network (in a stochastic ordering sense) than any other causal policy. In addition, for arbitrary general distributions of packet transmission times, the non-preemptive LGFS policy is shown to minimize the age processes at all nodes of the network among all non-preemptive work-conserving policies (again in a stochastic ordering sense). It is surprising that such simple policies can achieve optimality of the joint distribution of the age processes at all nodes even under arbitrary network topologies, as well as arbitrary packet generation and arrival times. These optimality results not only hold for the age processes, but also for any non-decreasing functional of the age processes.


title: Generation of ultra broadband coherent supercontinuum in tapered and dispersion managed silicon nanophotonic waveguides

abstract: Tapered and dispersion managed (DM) silicon nanophotonic waveguides are investigated for the generation of optimal ultra broadband supercontinuum (SC). DM waveguides are structures showing a longitudinally dependent group velocity dispersion that results from the variation of the waveguide width with the propagation distance. For the generation of optimal SC, a genetic algorithm has been used to find the best dispersion map. This allows for the generation of highly coherent supercontinuums that span over 1.14 octaves from 1300 nm to 2860 nm and 1.25 octaves from 1200 nm to 2870 nm at -20 dB level for the tapered and DM waveguides respectively, for a 2 $\mu$m, 200 fs and 6.4 pJ input pulse. The comparison of these two structures with the usually considered optimal fixed width waveguide shows that the SC is broader and flatter in the more elaborated DM waveguide, while the high coherence is ensured by the varying dispersion.


title: A Unified Stochastic Formulation of Dissipative Quantum Dynamics. II. Beyond Linear Response of Spin Baths

abstract: We use the "generalized hierarchical equation of motion" proposed in Paper I to study decoherence in a system coupled to a spin bath. The present methodology allows a systematic incorporation of higher order anharmonic effects of the bath in dynamical calculations. We investigate the leading order corrections to the linear response approximations for spin bath models. Two types of spin-based environments are considered: (1) a bath of spins discretized from a continuous spectral density and (2) a bath of physical spins such as nuclear or electron spins. The main difference resides with how the bath frequency and the system-bath coupling parameters are chosen to represent an environment. When discretized from a continuous spectral density, the system-bath coupling typically scales as $\sim 1/\sqrt{N_B}$ where $N_B$ is the number of bath spins. This scaling suppresses the non-Gaussian characteristics of the spin bath and justify the linear response approximations in the thermodynamic limit. For the physical spin bath models, system-bath couplings are directly deduced from spin-spin interactions with no reason to obey the $1/\sqrt{N_B}$ scaling. It is not always possible to justify the linear response approximations. Furthermore, if the spin-spin Hamiltonian and/or the bath parameters are highly symmetrical, these additional constraints generate non-Markovian and persistent dynamics that is beyond the linear response treatments.


title: A geometric Iwatsuka type effect in quantum layers

abstract: We study motion of a charged particle confined to Dirichlet layer of a fixed width placed into a homogeneous magnetic field. If the layer is planar and the field is perpendicular to it the spectrum consists of infinitely degenerate eigenvalues. We consider translationally invariant geometric perturbations and derive several sufficient conditions under which a magnetic transport is possible, that is, the spectrum, in its entirety or a part of it, becomes absolutely continuous.

categories: math-ph cond-mat.mes-hall math.MP math.SP quant-ph


title: Linear Bound for Majority Colourings of Digraphs

abstract: Given $\eta \in [0, 1]$, a colouring $C$ of $V(G)$ is an $\eta$-majority colouring if at most $\eta d^+(v)$ out-neighbours of $v$ have colour $C(v)$, for any $v \in V(G)$. We show that every digraph $G$ equipped with an assignment of lists $L$, each of size at least $k$, has a $2/k$-majority $L$-colouring. For even $k$ this is best possible, while for odd $k$ the constant $2/k$ cannot be replaced by any number less than $2/(k+1)$. This generalizes a result of Anholcer, Bosek and Grytczuk, who proved the cases $k=3$ and $k=4$ and gave a weaker result for general $k$.


title: New bi-Hamiltonian systems on the plane

abstract: We discuss several new bi-Hamiltonian integrable systems on the plane with integrals of motion of third, fourth and sixth order in momenta. The corresponding variables of separation, separated relations, compatible Poisson brackets and recursion operators are also presented in the framework of the Jacobi method.


title: Controllability of impulse controlled systems of heat equations coupled by constant matrices

abstract: This paper studies the approximate and null controllability for impulse controlled systems of heat equations coupled by a pair (A,B) of constant matrices. We present a necessary and sufficient condition for the approximate controllability, which is exactly Kalman's controllability rank condition of (A,B). We prove that when such a system is approximately controllable, the approximate controllability over an interval [0,T] can be realized by adding controls at arbitrary n different control instants 0<\tau_1<\tau_2<\cdots<\tau_n<T, provided that \tau_n-\tau_1<d_A, where d_A=\min\{\pi/|Im \lambda| : \lambda\in \sigma(A)\}. We also show that in general, such systems are not null controllable.


title: Semihierarchical quantum repeaters based on moderate lifetime quantum memories

abstract: The construction of large-scale quantum networks relies on the development of practical quantum repeaters. Many approaches have been proposed with the goal of outperforming the direct transmission of photons, but most of them are inefficient or difficult to implement with current technology. Here, we present a protocol that uses a semi-hierarchical structure to improve the entanglement distribution rate while reducing the requirement of memory time to a range of tens of milliseconds. This protocol can be implemented with a fixed distance of elementary links and fixed requirements on quantum memories, which are independent of the total distance. This configuration is especially suitable for scalable applications in large-scale quantum networks.


title: Binary systems with an RR Lyrae component - progress in 2016

abstract: In this contribution, we summarize the progress made in the investigation of binary candidates with an RR Lyrae component in 2016. We also discuss the actual status of the RRLyrBinCan database.


title: Lithium in stellar atmospheres: observations and theory

abstract: Of all the light elements, lithium is the most sensitive indicator of stellar evolution. This review discusses current data on the abundance of lithium in the atmospheres of A-, F-, G-, and K-stars of different types, as well as the consistency of these data with theoretical predictions. The variety of observed Li abundances is illustrated by the following objects in different stages of evolution: 1) Old stars in the galactic halo, which have a lithium abundance log e(Li)=2.2 (the "lithium plateau") that appears to be 0.5 dex lower than the primordial abundance predicted by cosmological models. 2) Young stars in the galactic disk, which have been used to estimate the contemporary initial lithium abundance log e(Li)=3.2 +/- 0.1 for stars in the Main sequence. Possible sources of lithium enrichment in the interstellar medium during evolution of the galaxy are discussed. 3) Evolving FGK dwarfs in the galactic disk, which have lower log e(Li) for lower effective temperature Teff and mass M. The "lithium dip" near Teff~6600 K in the distribution of log e(Li) with respect to Teff in old clusters is discussed. 4) FGK giants and supergiants, of which most have no lithium at all. This phenomenon is consistent with rotating star model calculations. 5) Lithium rich cold giants with log e(Li) > 2.0 , which form a small, enigmatic group. Theoretical models with rotation can explain the existence of these stars only in the case of low initial rotation velocities V0 < 50 km/s. In all other cases it is necessary to assume recent synthesis of lithium (capture of a giant planet is an alternative). 6) Magnetic Ap-stars, where lithium is concentrated in spots located at the magnetic poles. There the lithium abundance reaches log e(Li)=6. Discrepancies between observations and theory are noted for almost all the stars discussed in this review.


title: Irreducible components of the eigencurve of finite degree are finite over the weight space

abstract: Let p be a rational prime and N a positive integer which is prime to p. Let W be the p-adic weight space for GL_{2,Q}. Let C_N be the p-adic Coleman-Mazur eigencurve of tame level N. In this paper, we prove that any irreducible component of C_N which is of finite degree over W is in fact finite over W. Combined with an argument of Chenevier and a conjecture of Coleman-Mazur-Buzzard-Kilford (which has been proven in special cases, and for general quaternionic eigencurves) this shows that the only finite degree components of the eigencurve are the ordinary components.


title: The VMC Survey - XXIV. Signatures of tidally stripped stellar populations from the inner Small Magellanic Cloud

abstract: We study the luminosity function of intermediate-age red-clump stars using deep, near-infrared photometric data covering $\sim$ 20 deg$^2$ located throughout the central part of the Small Magellanic Cloud (SMC), comprising the main body and the galaxy's eastern wing, based on observations obtained with the VISTA Survey of the Magellanic Clouds (VMC). We identified regions which show a foreground population ($\sim$11.8 $\pm$ 2.0 kpc in front of the main body) in the form of a distance bimodality in the red-clump distribution. The most likely explanation for the origin of this feature is tidal stripping from the SMC rather than the extended stellar haloes of the Magellanic Clouds and/or tidally stripped stars from the Large Magellanic Cloud. The homogeneous and continuous VMC data trace this feature in the direction of the Magellanic Bridge and, particularly, identify (for the first time) the inner region ($\sim$ 2 -- 2.5 kpc from the centre) from where the signatures of interactions start becoming evident. This result provides observational evidence of the formation of the Magellanic Bridge from tidally stripped material from the SMC.


title: Hyperbolic Dispersion Dominant Regime Identified through Spontaneous Emission Variations near Metamaterial Interfaces

abstract: Surface plasmon polariton, hyberbolic dispersion of energy and momentum, and emission interference provide opportunities to control photoluminescence properties. However, the interplays between these regimes need to be understood to take advantage of them in optoelectronic applications. Here, we investigate broadband variations induced by hyperbolic metamaterial (HMM) multilayer nanostructures on the spontaneous emission of selected organic chromophores. Experimental and calculated spontaneous emission lifetimes are shown to vary non-monotonously near HMM interfaces. With the SPP and interference dominant regimes. With the HMM number of pairs used as the analysis parameter, the lifetime is shown to be independent of the number of pairs in the surface plasmon polaritons, and emission interference dominant regimes, while it decreases in the Hyperbolic Dispersion dominant regime. We also show that the spontaneous emission lifetime is similarly affected by transverse positive and transverse negative HMMs. This work has broad implications on the rational design of functional photonic surfaces to control the luminescence of semiconductor chromophores.


title: Logical Inferences with Contexts of RDF Triples

abstract: Logical inference, an integral feature of the Semantic Web, is the process of deriving new triples by applying entailment rules on knowledge bases. The entailment rules are determined by the model-theoretic semantics. Incorporating context of an RDF triple (e.g., provenance, time, and location) into the inferencing process requires the formal semantics to be capable of describing the context of RDF triples also in the form of triples, or in other words, RDF contextual triples about triples. The formal semantics should also provide the rules that could entail new contextual triples about triples. In this paper, we propose the first inferencing mechanism that allows context of RDF triples, represented in the form of RDF triples about triples, to be the first-class citizens in the model-theoretic semantics and in the logical rules. Our inference mechanism is well-formalized with all new concepts being captured in the model-theoretic semantics. This formal semantics also allows us to derive a new set of entailment rules that could entail new contextual triples about triples. To demonstrate the feasibility and the scalability of the proposed mechanism, we implement a new tool in which we transform the existing knowledge bases to our representation of RDF triples about triples and provide the option for this tool to compute the inferred triples for the proposed rules. We evaluate the computation of the proposed rules on a large scale using various real-world knowledge bases such as Bio2RDF NCBI Genes and DBpedia. The results show that the computation of the inferred triples can be highly scalable. On average, one billion inferred triples adds 5-6 minutes to the overall transformation process. NCBI Genes, with 20 billion triples in total, took only 232 minutes for the transformation of 12 billion triples and added 42 minutes for inferring 8 billion triples to the overall process.


title: An infinite dimensional KAM theorem with application to two dimensional completely resonant beam equation

abstract: In this paper we consider the completely resonant beam equation on \T^2 with cubic nonlinearity on a subspace of L^2 (\T^2) which will be explained later. We establish an abstract infinite dimensional KAM theorem and apply it to the completely resonant beam equation. We prove the existence of a class of Whitney smooth small amplitude quasi-periodic solutions corresponding to finite dimensional tori.


title: Sto\"ilow's theorem revisited

abstract: Sto\"ilow's theorem from 1928 states that a continuous, open, and light map between surfaces is a discrete map with a discrete branch set. This result implies that such maps between orientable surfaces are locally modeled by power maps $z\mapsto z^k$ and admit a holomorphic factorization. The purpose of this expository article is to give a proof of this classical theorem having readers in mind that are interested in continuous, open and discrete maps.


title: An abstract infinite dimensional KAM theorem with application to nonlinear higher dimensional Schrodinger equation systems

abstract: In this paper we consider nonlinear Schrodinger systems with periodic boundary condition in high dimension. We establish an abstract infinite dimensional KAM theorem and apply it to the nonlinear Schrodinger equation systems with real Fourier Multiplier. By establishing a block-diagonal normal form, We prove the existence of a class of Whitney smooth small amplitude quasi-periodic solutions corresponding to finite dimensional invariant tori of an associated infinite dimensional dynamical system.


title: In-medium parton branching beyond eikonal approximation

abstract: The description of the in-medium modifications of partonic showers has been at the forefront of current theoretical and experimental efforts in heavy-ion collisions. It provides a unique laboratory to extend our knowledge frontier of the theory of the strong interactions, and to assess the properties of the hot and dense medium (QGP) that is produced in ultra-relativistic heavy-ion collisions at RHIC and the LHC. The theory of jet quenching, a commonly used alias for the modifications of the parton branching resulting from the interactions with the QGP, has been significantly developed over the last years. Within a weak coupling approach, several elementary processes that build up the parton shower evolution, such as single gluon emissions, interference effects between successive emissions and corrections to radiative energy loss of massive quarks, have been addressed both at eikonal accuracy and beyond by taking into account the Brownian motion that high-energy particles experience when traversing a hot and dense medium. In this work, by using the setup of single gluon emission from a color correlated quark-antiquark pair in a singlet state ($q\bar{q}$ antenna), we calculate the in-medium gluon radiation spectrum beyond the eikonal approximation. The results show that we are able to factorize broadening effects from the modifications of the radiation process itself. This constitutes the final proof that a probabilistic picture of the parton shower evolution holds even in the presence of a QGP.


title: Cyclotomic p-adic multi-zeta values

abstract: The cyclotomic $p$-adic multi-zeta values are the $p$-adic periods of $\pi_{1}(\mathbb{G}_{m} \setminus \mu_{M},\cdot),$ the unipotent fundamental group of the multiplicative group minus the $M$-th roots of unity. In this paper, we compute the cyclotomic $p$-adic multi-zeta values at all depths. This paper generalizes the results in [6] and [7]. Since the main result gives explicit formulas we expect it to be useful in proving non-vanishing and transcendence results for these $p$-adic periods and also, through the use of $p$-adic Hodge theory, in proving non-triviality results for the corresponding $p$-adic Galois representations.


title: Using LLVM-based JIT Compilation in Genetic Programming

abstract: The paper describes an approach to implementing genetic programming, which uses the LLVM library to just-in-time compile/interpret the evolved abstract syntax trees. The solution is described in some detail, including a parser (based on FlexC++ and BisonC++) that can construct the trees from a simple toy language with C-like syntax. The approach is compared with a previous implementation (based on direct execution of trees using polymorphic functors) in terms of execution speed.


title: Non-minimally coupled gravity and vacuum stability

abstract: We investigate the properties of vacuum decay taking into account a non-minimal coupling to gravity. We extend the standard thin-wall solution to include the non-minimal coupling and verify its validity by comparison with a full numerical study. We also investigate the implications of a~large cosmological constant whose influence on the geometry boosts the tunneling rate. Our analysis shows that the influence of the non-minimal coupling differs significantly between the cases of Minkowski and deSitter backgrounds.


title: Fast vortex wall motion in wide Permalloy strips from double switching of the vortex core

abstract: We study vortex domain wall dynamics in wide Permalloy strips driven by applied magnetic fields and spin-polarized electric currents. As recently reported [V. Est\'evez and L. Laurson, Phys. Rev. B, 93, 064403 (2016)], for sufficiently wide strips and above a threshold field, periodic dynamics of the vortex core are localized in the vicinity of one of the strip edges, and the velocity drop typically observed for narrow strips is replaced by a high-velocity plateau. Here, we analyze this behavior in more detail by means of micromagnetic simulations. We show that the high-velocity plateau originates from a repeated double switching of the magnetic vortex core, underlying the periodic vortex core dynamics in the vicinity of the strip edge, i.e., the "attraction-repulsion" effect. We also discuss the corresponding dynamics driven by spin-polarized currents, as well as the effect of including quenched random structural disorder to the system.