title: pth moment noise-to-state stability of stochastic differential equations with persistent noise

abstract: This paper studies the stability properties of stochastic differential equations subject to persistent noise (including the case of additive noise), which is noise that is present even at the equilibria of the underlying differential equation and does not decay with time. The class of systems we consider exhibit disturbance attenuation outside a closed, not necessarily bounded, set. We identify conditions, based on the existence of Lyapunov functions, to establish the noise-to-state stability in probability and in pth moment of the system with respect to a closed set. As part of our analysis, we study the concept of two functions being proper with respect to each other formalized via a pair of inequalities with comparison functions. We show that such inequalities define several equivalence relations for increasingly strong refinements on the comparison functions. We also provide a complete characterization of the properties that a pair of functions must satisfy to belong to the same equivalence class. This characterization allows us to provide checkable conditions to determine whether a function satisfies the requirements to be a strong NSS-Lyapunov function in probability or a pth moment NSS-Lyapunov function. Several examples illustrate our results.


title: Comment on two recent arXiv postings [arXiv:1412.2200 and arXiv:1412.8244]

abstract: This is a comment on two recent arXiv postings.


title: Occurrence of high-speed solar wind streams over the Grand Modern Maximum

abstract: In the declining phase of the solar cycle, when the new-polarity fields of the solar poles are strengthened by the transport of same-signed magnetic flux from lower latitudes, the polar coronal holes expand and form non-axisymmetric extensions toward the solar equator. These extensions enhance the occurrence of high-speed solar wind streams (HSS) and related co-rotating interaction regions in the low-latitude heliosphere, and cause moderate, recurrent geomagnetic activity in the near-Earth space. Here, using a novel definition of geomagnetic activity at high (polar cap) latitudes and the longest record of magnetic observations at a polar cap station, we calculate the annually averaged solar wind speeds as proxies for the effective annual occurrence of HSS over the whole Grand Modern Maximum (GMM) from 1920s onwards. We find that a period of high annual speeds (frequent occurrence of HSS) occurs in the declining phase of each solar cycle 16-23. For most cycles the HSS activity clearly maximizes during one year, suggesting that typically only one strong activation leading to a coronal hole extension is responsible for the HSS maximum. We find that the most persistent HSS activity occurred in the declining phase of solar cycle 18. This suggests that cycle 19, which marks the sunspot maximum period of the GMM, was preceded by exceptionally strong polar fields during the previous sunspot minimum. This gives interesting support for the validity of solar dynamo theory during this dramatic period of solar magnetism.


title: Two-dimensional volume-frozen percolation: exceptional scales

abstract: We study a percolation model on the square lattice, where clusters "freeze" (stop growing) as soon as their volume (i.e. the number of sites they contain) gets larger than N, the parameter of the model. A model where clusters freeze when they reach diameter at least N was studied in earlier papers. Using volume as a way to measure the size of a cluster - instead of diameter - leads, for large N, to a quite different behavior (contrary to what happens on the binary tree, where the volume model and the diameter model are "asymptotically the same"). In particular, we show the existence of a sequence of "exceptional" length scales.


title: A 3D Search for Companions to 12 Nearby M-Dwarfs

abstract: We present a carefully vetted equatorial ($\pm$ 30$^\circ$ Decl.) sample of all known single (within 4'') mid M-dwarfs (M2.5V-M8.0V) extending out to 10 pc; their proximity and low masses make them ideal targets for planet searches. For this sample of 58 stars, we provide V$_J$, R$_{KC}$, I$_{KC}$ photometry, new low dispersion optical ($6000 - 9000$\AA) spectra from which uniform spectral types are determined, multi-epoch H$\alpha$ equivalent widths, and gravity sensitive $Na\,I$ indices. For 12 of these 58 stars, strict limits are placed on the presence of stellar and sub-stellar companions, based on a pioneering program described here that utilizes precise infrared radial velocities and optical astrometric measurements in an effort to search for Jupiter-mass, brown dwarf and stellar-mass companions. Our infrared radial velocity precision using CSHELL at NASA's IRTF is $\sim$90 m s$^{-1}$ over timescales from 13 days to 5 years. With our spectroscopic results the mean companion masses that we rule out of existence are 1.5 M$_{JUP}$ or greater in 10 day orbital periods and 7 M$_{JUP}$ or greater in 100 day orbital periods. We use these spectra to determine rotational velocities and absolute radial velocities of these twelve stars. Our mean astrometric precision using RECONS data from 0.9-m telescope at Cerro Tololo Inter-American Observatory is $\sim$3 milli-arcseconds over baselines ranging from 9 to 13 years. With our astrometric results the mean companion masses that we rule out of existence are greater than 11.5 M$_{JUP}$ with an orbital period of 4 years and greater than 7.5 M$_{JUP}$ with an orbital period of 8 years. Although we do not detect companions around our sub-sample of 12 stars, we demonstrate that our two techniques probe a regime that is commonly missed in other companion searches of late type stars.


title: Current status of the Standard Model CKM fit and constraints on $\Delta F=2$ New Physics

abstract: This letter summarises the status of the global fit of the CKM parameters within the Standard Model performed by the CKMfitter group. Special attention is paid to the inputs for the CKM angles $\alpha$ and $\gamma$ and the status of $B_s\to\mu\mu$ and $B_d\to \mu\mu$ decays. We illustrate the current situation for other unitarity triangles. We also discuss the constraints on generic $\Delta F=2$ New Physics. All results have been obtained with the CKMfitter analysis package, featuring the frequentist statistical approach and using Rfit to handle theoretical uncertainties.


title: Experimental Simulation of Closed Timelike Curves

abstract: Closed timelike curves are among the most controversial features of modern physics. As legitimate solutions to Einstein's field equations, they allow for time travel, which instinctively seems paradoxical. However, in the quantum regime these paradoxes can be resolved leaving closed timelike curves consistent with relativity. The study of these systems therefore provides valuable insight into non-linearities and the emergence of causal structures in quantum mechanics-essential for any formulation of a quantum theory of gravity. Here we experimentally simulate the non-linear behaviour of a qubit interacting unitarily with an older version of itself, addressing some of the fascinating effects that arise in systems traversing a closed timelike curve. These include perfect discrimination of non-orthogonal states and, most intriguingly, the ability to distinguish nominally equivalent ways of preparing pure quantum states. Finally, we examine the dependence of these effects on the initial qubit state, the form of the unitary interaction, and the influence of decoherence.


title: Damping of long wavelength collective modes in spinor Bose-Fermi mixtures

abstract: Using an effective field theory we describe the low energy bosonic excitations in a three dimensional ultra-cold mixture of spin-1 bosons and spin-1/2 fermions. We establish an interesting fermionic excitation induced generic damping of the usual undamped long wavelength bosonic collective Goldstone modes. Two states with bosons forming either a ferromagnetic or polar superfluid are studied. The linear dispersion of the bosonic Bogoliubov excitations is preserved with a renormalized sound velocity. For the polar superfluid we find both gapless modes (density and spin) are damped, whereas in the ferromagnetic superfluid we find the density (spin) mode is (not) damped. We argue quite generally that this holds for any mixture of bosons and fermions that are coupled through at least a density-density interaction. We discuss the implications of our many-body interaction results for experiments on Bose-Fermi mixtures.


title: A Categorical Semantics for Linear Logical Frameworks

abstract: A type theory is presented that combines (intuitionistic) linear types with type dependency, thus properly generalising both intuitionistic dependent type theory and full linear logic. A syntax and complete categorical semantics are developed, the latter in terms of (strict) indexed symmetric monoidal categories with comprehension. Various optional type formers are treated in a modular way. In particular, we will see that the historically much-debated multiplicative quantifiers and identity types arise naturally from categorical considerations. These new multiplicative connectives are further characterised by several identities relating them to the usual connectives from dependent type theory and linear logic. Finally, one important class of models, given by families with values in some symmetric monoidal category, is investigated in detail.


title: Antenna array characterization via radio interferometry observation of astronomical sources

abstract: We present an in-situ antenna characterization method and results for a "low-frequency" radio astronomy engineering prototype array, characterized over the 75-300 MHz frequency range. The presence of multiple cosmic radio sources, particularly the dominant Galactic noise, makes in-situ characterization at these frequencies challenging; however, it will be shown that high quality measurement is possible via radio interferometry techniques. This method is well-known in the radio astronomy community but seems less so in antenna measurement and wireless communications communities, although the measurement challenges involving multiple undesired sources in the antenna field-of-view bear some similarities. We discuss this approach and our results with the expectation that this principle may find greater application in related fields.


title: Supersymmetric non conservative systems

abstract: We give the generalization of a recent variational formulation for nonconservative classical mechanics, for fermionic and sypersymmetric systems. Both cases require slightly modified boundary conditions. The supersymmetric version is given in the superfield formalism. The corresponding Noether theorem is formulated. As expected, like the energy, the supersymmetric charges are not conserved.


title: Robust Hypothesis Testing with $\alpha$-Divergence

abstract: A robust minimax test for two composite hypotheses, which are determined by the neighborhoods of two nominal distributions with respect to a set of distances - called $\alpha-$divergence distances, is proposed. Sion's minimax theorem is adopted to characterize the saddle value condition. Least favorable distributions, the robust decision rule and the robust likelihood ratio test are derived. If the nominal probability distributions satisfy a symmetry condition, the design procedure is shown to be simplified considerably. The parameters controlling the degree of robustness are bounded from above and the bounds are shown to be resulting from a solution of a set of equations. The simulations performed evaluate and exemplify the theoretical derivations.


title: Layouts of Expander Graphs

abstract: Bourgain and Yehudayoff recently constructed $O(1)$-monotone bipartite expanders. By combining this result with a generalisation of the unraveling method of Kannan, we construct 3-monotone bipartite expanders, which is best possible. We then show that the same graphs admit 3-page book embeddings, 2-queue layouts, 4-track layouts, and have simple thickness 2. All these results are best possible.


title: Stochastic Block Model and Community Detection in the Sparse Graphs: A spectral algorithm with optimal rate of recovery

abstract: In this paper, we present and analyze a simple and robust spectral algorithm for the stochastic block model with $k$ blocks, for any $k$ fixed. Our algorithm works with graphs having constant edge density, under an optimal condition on the gap between the density inside a block and the density between the blocks. As a co-product, we settle an open question posed by Abbe et. al. concerning censor block models.


title: The uTPC Method: Improving the Position Resolution of Neutron Detectors Based on MPGDs

abstract: Due to the Helium-3 crisis, alternatives to the standard neutron detection techniques are becoming urgent. In addition, the instruments of the European Spallation Source (ESS) require advances in the state of the art of neutron detection. The instruments need detectors with excellent neutron detection efficiency, high-rate capabilities and unprecedented spatial resolution. The Macromolecular Crystallography instrument (NMX) requires a position resolution in the order of 200 um over a wide angular range of incoming neutrons. Solid converters in combination with Micro Pattern Gaseous Detectors (MPGDs) are proposed to meet the new requirements. Charged particles rising from the neutron capture have usually ranges larger than several millimetres in gas. This is apparently in contrast with the requirements for the position resolution. In this paper, we present an analysis technique, new in the field of neutron detection, based on the Time Projection Chamber (TPC) concept. Using a standard Single-GEM with the cathode coated with 10B4C, we extract the neutron interaction point with a resolution of better than sigma = 200 um.


title: Extremal dichotomy for uniformly hyperbolic systems

abstract: We consider the extreme value theory of a hyperbolic toral automorphism $T: \mathbb{T}^2 \to \mathbb{T}^2$ showing that if a H\"older observation $\phi$ which is a function of a Euclidean-type distance to a non-periodic point $\zeta$ is strictly maximized at $\zeta$ then the corresponding time series $\{\phi\circ T^i\}$ exhibits extreme value statistics corresponding to an iid sequence of random variables with the same distribution function as $\phi$ and with extremal index one. If however $\phi$ is strictly maximized at a periodic point $q$ then the corresponding time-series exhibits extreme value statistics corresponding to an iid sequence of random variables with the same distribution function as $\phi$ but with extremal index not equal to one. We give a formula for the extremal index (which depends upon the metric used and the period of $q$). These results imply that return times are Poisson to small balls centered at non-periodic points and compound Poisson for small balls centered at periodic points.


title: Hydrogen Fluoride Capture by Imidazolium Acetate Ionic Liquid

abstract: Extraction of hydrofluoric acid (HF) from oils is a drastically important problem in petroleum industry, since HF causes quick corrosion of pipe lines and brings severe health problems to humanity. Some ionic liquids (ILs) constitute promising scavenger agents thanks to strong binding to polar compounds and tunability. PM7-MD simulations and hybrid density functional theory are employed here to consider HF capture ability of ILs. Discussing the effects and impacts of the cation and the anion separately and together, I will evaluate performance of imidazolium acetate and outline systematic search guidelines for efficient adsorption and extraction of HF.


title: The long-lasting optical afterglow plateau of short burst GRB 130912A

abstract: The short burst GRB 130912A was detected by Swift, Fermi satellites and several ground-based optical telescopes. Its X-ray light curve decayed with time normally. The optical emission, however, displayed a long term plateau, which is the longest one in current short GRB observations. In this work we examine the physical origin of the X-ray and optical emission of this peculiar event. We find that the canonical forward shock afterglow emission model can account for the X-ray and optical data self-consistently and the energy injection model that has been widely adopted to interpret the shallowly-decaying afterglow emission is not needed. We also find that the burst was born in a very-low density interstellar medium, consistent with the compact object merger model. Significant fractions of the energy of the forward shock have been given to accelerate the non-thermal electrons and amplify the magnetic fields (i.e., $\epsilon_{\rm e}\sim 0.37$ and $\epsilon_{\rm B}\sim 0.16$, respectively), which are much larger than those inferred in most short burst afterglow modeling and can explain why the long-lasting optical afterglow plateau is rare in short GRBs.


title: Two-photon coherent control of atomic collisions by light with entangled polarization

abstract: We describe a new method of coherent optical control of internal dynamics of atomic collisions by means of two correlated light beams having entangled polarizations. We show that if excitation of a colliding pair of atoms is by two photons having entangled polarizations, it is possible to redirect the output fragments of the collision into certain channels with a selected type of internal transition symmetry. The transition symmetry is defined in the body-fixed coordinate frame which has random and originally unknown orientation in space.


title: On the stable hom relation and stable degenerations of Cohen-Macaulay modules

abstract: We study the stable hom relation for Cohen-Macaulay modules over Gorenstein local algebras. We give the sufficient condition to make the stable hom relation a partial order when the base algebra is of finite representation type. As an application, we give the description of stable degenerations of Cohen-Macaulay modules over simple singularities of several types by using the stable hom relation.


title: Quantum Spectrum Testing

abstract: In this work, we study the problem of testing properties of the spectrum of a mixed quantum state. Here one is given $n$ copies of a mixed state $\rho\in\mathbb{C}^{d\times d}$ and the goal is to distinguish whether $\rho$'s spectrum satisfies some property $\mathcal{P}$ or is at least $\epsilon$-far in $\ell_1$-distance from satisfying $\mathcal{P}$. This problem was promoted in the survey of Montanaro and de Wolf under the name of testing unitarily invariant properties of mixed states. It is the natural quantum analogue of the classical problem of testing symmetric properties of probability distributions. Here, the hope is for algorithms with subquadratic copy complexity in the dimension $d$. This is because the "empirical Young diagram (EYD) algorithm" can estimate the spectrum of a mixed state up to $\epsilon$-accuracy using only $\widetilde{O}(d^2/\epsilon^2)$ copies. In this work, we show that given a mixed state $\rho\in\mathbb{C}^{d\times d}$: (i) $\Theta(d/\epsilon^2)$ copies are necessary and sufficient to test whether $\rho$ is the maximally mixed state, i.e., has spectrum $(\frac1d, ..., \frac1d)$; (ii) $\Theta(r^2/\epsilon)$ copies are necessary and sufficient to test with one-sided error whether $\rho$ has rank $r$, i.e., has at most $r$ nonzero eigenvalues; (iii) $\widetilde{\Theta}(r^2/\Delta)$ copies are necessary and sufficient to distinguish whether $\rho$ is maximally mixed on an $r$-dimensional or an $(r+\Delta)$-dimensional subspace; and (iv) The EYD algorithm requires $\Omega(d^2/\epsilon^2)$ copies to estimate the spectrum of $\rho$ up to $\epsilon$-accuracy, nearly matching the known upper bound. In addition, we simplify part of the proof of the upper bound. Our techniques involve the asymptotic representation theory of the symmetric group; in particular Kerov's algebra of polynomial functions on Young diagrams.


title: Integrability Conditions for Killing-Yano Tensors and Maximally Symmetric Spaces in the Presence of Torsion

abstract: The integrability conditions for the existence of Killing-Yano tensors or, equivalently, covariantly closed conformal Killing-Yano tensors, in the presence of torsion are worked out. As an application, all metrics and torsions compatible with the existence of a Killing-Yano tensor of order n-1 are obtained. Finally, the issue of defining a maximally symmetric space with respect to connections with torsion is addressed.


title: Transverse Wobbling in $^{135}$Pr

abstract: A pair of transverse wobbling bands has been observed in the nucleus $^{135}$Pr. The wobbling is characterized by $\Delta I$ =1, E2 transitions between the bands, and a decrease in the wobbling energy confirms its transverse nature. Additionally, a transition from transverse wobbling to a three-quasiparticle band comprised of strong magnetic dipole transitions is observed. These observations conform well to results from calculations with the Tilted Axis Cranking (TAC) model and the Quasiparticle Triaxial Rotor (QTR) Model.


title: Maximin Safety: When Failing to Lose is Preferable to Trying to Win

abstract: We present a new decision rule, \emph{maximin safety}, that seeks to maintain a large margin from the worst outcome, in much the same way minimax regret seeks to minimize distance from the best. We argue that maximin safety is valuable both descriptively and normatively. Descriptively, maximin safety explains the well-known \emph{decoy effect}, in which the introduction of a dominated option changes preferences among the other options. Normatively, we provide an axiomatization that characterizes preferences induced by maximin safety, and show that maximin safety shares much of the same behavioral basis with minimax regret.


title: Edge reconstruction in armchair phosphorene nanoribbons revealed by discontinuous Galerkin density functional theory

abstract: With the help of our recently developed massively parallel DGDFT (Discontinuous Galerkin Density Functional Theory) methodology, we perform large-scale Kohn-Sham density functional theory calculations on phosphorene nanoribbons with armchair edges (ACPNRs) containing a few thousands to ten thousand atoms. The use of DGDFT allows us to systematically achieve conventional plane wave basis set type of accuracy, but with a much smaller number (about 15) of adaptive local basis (ALB) functions per atom for this system. The relatively small number degrees of freedom required to represent the Kohn-Sham Hamiltonian, together with the use of the pole expansion the selected inversion (PEXSI) technique that circumvents the need to diagonalize the Hamiltonian, result in a highly efficient and scalable computational scheme for analyzing the electronic structures of ACPNRs as well as its dynamics. The total wall clock time for calculating the electronic structures of large-scale ACPNRs containing 1080-10800 atoms is only 10-25 s per self-consistent field (SCF) iteration, with accuracy fully comparable to that obtained from conventional planewave DFT calculations. For the ACPNR system, we observe that the DGDFT methodology can scale to 5,000-50,000 processors. We use DGDFT based ab-initio molecular dynamics (AIMD) calculations to study the thermodynamic stability of ACPNRs. Our calculations reveal that a 2 * 1 edge reconstruction appears in ACPNRs at room temperature.


title: New simultaneous methods for finding all zeros of a polynomial

abstract: The purpose of this paper is to present three new methods for finding all simple zeros of polynomials simultaneously. First, we give a new method for finding simultaneously all simple zeros of polynomials constructed by applying the Weierstrass method to the zero in the trapezoidal Newton's method, and prove the convergence of the method. We also present two modified Newton's methods combined with the derivative-free method, which are constructed by applying the derivative-free method to the zero in the trapezoidal Newton's method and the midpoint Newton's method, respectively. Finally, we give a numerical comparison between various simultaneous methods for finding zeros of a polynomial.


title: Simplifying the Reinsch algorithm for the Baker-Campbell-Hausdorff series

abstract: The Baker-Campbell-Hausdorff series computes the quantity \begin{equation*} Z(X,Y)=\ln\left( e^X e^Y \right) = \sum_{n=1}^\infty z_n(X,Y), \end{equation*} where $X$ and $Y$ are not necessarily commuting, in terms of homogeneous multinomials $z_n(X,Y)$ of degree $n$. (This is essentially equivalent to computing the so-called Goldberg coefficients.) The Baker-Campbell-Hausdorff series is a general purpose tool of wide applicability in mathematical physics, quantum physics, and many other fields. The Reinsch algorithm for the truncated series permits one to calculate up to some fixed order $N$ by using $(N+1)\times(N+1)$ matrices. We show how to further simplify the Reinsch algorithm, making implementation (in principle) utterly straightforward. This helps provide a deeper understanding of the Goldberg coefficients and their properties. For instance we establish strict bounds (and some equalities) on the number of non-zero Goldberg coefficients. Unfortunately, we shall see that the number of terms in the multinomial $z_n(X,Y)$ often grows very rapidly (in fact exponentially) with the degree $n$. We also present some closely related results for the symmetric product \begin{equation*} S(X,Y)=\ln\left( e^{X/2} e^Y e^{X/2} \right) = \sum_{n=1}^\infty s_n(X,Y). \end{equation*} Variations on these themes are straightforward. For instance, one can just as easily consider the series \begin{equation*} L(X,Y)=\ln\left( e^{X} e^Y e^{-X} e^{-Y}\right) = \sum_{n=1}^\infty \ell_n(X,Y). \end{equation*} This type of series is of interest, for instance, when considering parallel transport around a closed curve. Several other related series are investigated.


title: Musings on the theory that variation in cancer risk among tissues can be explained by the number of divisions of normal stem cells

abstract: This manuscript has been written to address questions related to our recent publication (Science 347:78-81, 2015). We appreciate the many reactions to this paper that have been communicated to us, either privately or publicly. The following addresses several of the most important statistical and technical issues related to our analysis and conclusions. Our responses to non-technical questions are available at http://www.hopkinsmedicine.org/news/media/releases/bad_luck_of_random_mutations_plays_predominant_role_in_cancer_study_shows


title: Hyperbolic Lagrangian coherent structures align with contours of path-averaged scalars

abstract: While inequality (9) is mathematically correct, it does not imply alignment between path-averaged scalars and the hyperbolic LCSs.


title: Extensions of isometric embeddings of Pseudo-Euclidean metric polyhedra

abstract: We extend the results of B. Minemyer by showing that any indefinite metric polyhedron (either compact or not) with the vertex degree bounded from above admits an isometric simplicial embedding into a Minkowski space of the lowest possible dimension. We provide a simple algorithm of constructing such embeddings. We also show that every partial simplicial isometric embedding of such space in general position extends to a simplicial isometric embedding of the whole space.


title: Light curve solutions of the ultrashort-period $Kepler$ binaries

abstract: We carried out light curve solutions of the ultrashort-period binaries with MS components observed by $Kepler$. All six targets turned out almost in thermal contact with contact or slightly overcontact configurations. Two of them, KID 4921906 and KID 6309193, are not eclipsing but reveal ellipsoidal and spot variability. One of the components of KID 8108785 exhibits inherent, quasi-sinusoidal, small-amplitude variability. KID 12055255 turned out a very rare case of ultrashort-period overcontact binary consisting of two M dwarfs. Our modeling indicated that the variability of KID 9532219 is due to eclipses but not to $\delta$ Sct pulsations as it was previously supposed.


title: Defining Data Science

abstract: Data science is gaining more and more and widespread attention, but no consensus viewpoint on what data science is has emerged. As a new science, its objects of study and scientific issues should not be covered by established sciences. Data in cyberspace have formed what we call datanature. In the present paper, data science is defined as the science of exploring datanature.


title: Modular Dynamics of Financial Market Networks

abstract: The financial market is a complex dynamical system composed of a large variety of intricate relationships between several entities, such as banks, corporations and institutions. At the heart of the system lies the stock exchange mechanism, which establishes a time-evolving network of trades among companies and individuals. Such network can be inferred through correlations between time series of companies stock prices, allowing the overall system to be characterized by techniques borrowed from network science. Here we study the presence of communities in the inferred stock market network, and show that the knowledge about the communities alone can provide a nearly complete representation of the system topology. This is done by defining a simple null model, a randomized version of the studied network sharing only the sizes and interconnectivity between communities observed. We show that many topological characteristics of the inferred networks are carried over the networks generated by the null model. In particular, we find that in periods of instability, such as during a financial crisis, the network strays away from a state of well-defined community structure to a much more uniform topological organization. We show that the framework presented here provides a good null model representation of topological variations taking place in the market during crises. Also, the general approach used in this work can be extended to other systems.


title: Pilot-Abstraction: A Valid Abstraction for Data-Intensive Applications on HPC, Hadoop and Cloud Infrastructures?

abstract: HPC environments have traditionally been designed to meet the compute demand of scientific applications and data has only been a second order concern. With science moving toward data-driven discoveries relying more on correlations in data to form scientific hypotheses, the limitations of HPC approaches become apparent: Architectural paradigms such as the separation of storage and compute are not optimal for I/O intensive workloads (e.g. for data preparation, transformation and SQL). While there are many powerful computational and analytical libraries available on HPC (e.g. for scalable linear algebra), they generally lack the usability and variety of analytical libraries found in other environments (e.g. the Apache Hadoop ecosystem). Further, there is a lack of abstractions that unify access to increasingly heterogeneous infrastructure (HPC, Hadoop, clouds) and allow reasoning about performance trade-offs in this complex environment. At the same time, the Hadoop ecosystem is evolving rapidly and has established itself as de-facto standard for data-intensive workloads in industry and is increasingly used to tackle scientific problems. In this paper, we explore paths to interoperability between Hadoop and HPC, examine the differences and challenges, such as the different architectural paradigms and abstractions, and investigate ways to address them. We propose the extension of the Pilot-Abstraction to Hadoop to serve as interoperability layer for allocating and managing resources across different infrastructures. Further, in-memory capabilities have been deployed to enhance the performance of large-scale data analytics (e.g. iterative algorithms) for which the ability to re-use data across iterations is critical. As memory naturally fits in with the Pilot concept of retaining resources for a set of tasks, we propose the extension of the Pilot-Abstraction to in-memory resources.


title: Spatially resolved resonant tunneling on single atoms in silicon

abstract: The ability to control single dopants in solid-state devices has opened the way towards reliable quantum computation schemes. In this perspective it is essential to understand the impact of interfaces and electric fields, inherent to address coherent electronic manipulation, on the dopants atomic scale properties. This requires both fine energetic and spatial resolution of the energy spectrum and wave-function, respectively. Here we present an experiment fulfilling both conditions: we perform transport on single donors in silicon close to a vacuum interface using a scanning tunneling microscope (STM) in the single electron tunneling regime. The spatial degrees of freedom of the STM tip provide a versatility allowing a unique understanding of electrostatics. We obtain the absolute energy scale from the thermal broadening of the resonant peaks, allowing to deduce the charging energies of the donors. Finally we use a rate equations model to derive the current in presence of an excited state, highlighting the benefits of the highly tunable vacuum tunnel rates which should be exploited in further experiments. This work provides a general framework to investigate dopant-based systems at the atomic scale.


title: Ultrashort-period MS eclipsing systems. New observations and light curve solutions of six NSVS binaries

abstract: We carried out photometric and low-resolution spectral observations of six eclipsing ultrashort-period binaries with MS components. The light curve solutions of the Rozhen observations show that all targets are overcontact systems. We found well-defined empirical relation "period -- semi-major axis" for the short-period binaries and used it for estimation of the global parameters of the targets. Our results revealed that NSVS 925605 is quite interesting target: (a) it is one of a few contact binaries with M components; (b) it exhibits high activity (emission in H$\alpha$ line, X-ray emission, large cool spots, non-Planck energy distribution); (c) its components differ in temperature by 700 K. All appearances of high magnetic activity and huge fillout factor (0.7) of NSVS 925605 might be assumed as a precursor of the predicted merging of close magnetic binaries. Another unusual binary is NSVS 2700153 which reveals considerable long-term variability.


title: Quantum Noises, Physical Realizability and Coherent Quantum Feedback Control

abstract: Physical Realizability addresses the question of whether it is possible to implement a given linear time invariant (LTI) system as a quantum system. A given synthesized quantum controller described by a set of stochastic differential equations does not necessarily correspond to a physically meaningful quantum system. However, if additional quantum noises are permitted in the implementation, it is always possible to implement an arbitrary LTI system as a quantum system. In this paper, we give an expression for the number of introduced noise channels required to implement a given LTI system as a quantum system. We then consider the special case where only the transfer function to be implemented is of interest. We give results showing when it is possible to implement a transfer function as a quantum system by introducing the same number of quantum noises as there are system outputs. Finally, we demonstrate the utility of these results by providing an algorithm for obtaining a suboptimal solution to a coherent quantum LQG control problem.


title: A note on mass-minimising extensions

abstract: A conjecture related to the Bartnik quasilocal mass, is that the infimum of the ADM energy, over an appropriate space of extensions to a compact 3-manifold with boundary, is realised by a static metric. It was shown by Corvino [Comm. Math. Phys. 214(1), (2000)] that if the infimum is indeed achieved, then it is achieved by a static metric; however, the more difficult question of whether or not the infimum is achieved, is still an open problem. Bartnik [Comm. Anal. Geom. 13(5), (2005)] then proved that critical points of the ADM mass, over the space of solutions to the Einstein constraints on an asymptotically flat manifold without boundary, correspond to stationary solutions. In that article, he stated that it should be possible to use a similar construction to provide a more natural proof of Corvino's result. In the first part of this note, we discuss the required modifications to Bartnik's argument to adapt it to include a boundary. Assuming that certain results concerning a Hilbert manifold structure for the space of solutions carry over to the case considered here, we then demonstrate how Bartnik's proof can be modified to consider the simpler case of scalar-flat extensions and obtain Corvino's result. In the second part of this note, we consider a space of extensions in a fixed conformal class. Sufficient conditions are given to ensure that the infimum is realised within this class.


title: Vertex-primitive digraphs having vertices with almost equal neighbourhoods

abstract: We consider vertex-primitive digraphs having two vertices with almost equal neighbourhoods (that is, the set of vertices that are neighbours of one but not the other is small). We prove a structural result about such digraphs and then apply it to answer a question of Ara\'ujo and Cameron about synchronising groups.


title: Antiferromagnetism of Ni2NbBO6 with S = 1 dimer quasi-one-dimensional armchair chains

abstract: Long range antiferromagnetic (AFM) ordering of Ni spins in Ni2NbBO6 has been studied with single crystal from spin susceptibility measurement and comparedwith the ab initio calculation results consistently. Below TN = 23.5 K, the S = 1 spins align along the a direction for edge-shared NiO6 octahedra which form crystallographic armchair chains along the b direction. The isothermal magnetization M(H) below TN shows spin-flop transition for magnetic field above 36 kOe along the a axis,which indicates the spin anisotropy is along the a direction. The electronic and magnetic structures of Ni2NbBO6 have also been explored theoretically using density functional theory with generalized gradient approximation plus on-site Coulomb interaction (U). These calculations support the experimentally observed antiferromagnetism of Ni2NbBO6. In particular, the long range AFM ordering below TN can be dissected into armchair chains which consists of S = 1 dimers of J2 = 2.43 meV with ferromagnetic (FM) intrachain and interchain couplings of size half |J2|.


title: "Dust, Ice, and Gas In Time" (DIGIT) Herschel Observations of GSS30-IRS1 in Ophiuchus

abstract: As a part of the "Dust, Ice, and Gas In Time" (DIGIT) key program on Herschel, we observed GSS30-IRS1, a Class I protostar located in Ophiuchus (d = 120 pc), with Herschel/Photodetector Array Camera and Spectrometer (PACS). More than 70 lines were detected within a wavelength range from 50 micron to 200 micron, including CO, H2O, OH, and two atomic [O I] lines at 63 and 145 micron. The [C II] line, known as a tracer of externally heated gas by the interstellar radiation field, is also detected at 158 micron. All lines, except [O I] and [C II], are detected only at the central spaxel of 9.4" X 9.4". The [O I] emissions are extended along a NE-SW orientation, and the [C II] line is detected over all spaxels, indicative of external PDR. The total [C II] intensity around GSS30 reveals that the far-ultraviolet radiation field is in the range of 3 to 20 G0, where G0 is in units of the Habing Field, 1.6 X 10^{-3} erg cm^{-2} s^{-1}. This enhanced external radiation field heats the envelope of GSS30-IRS1, causing the continuum emission to be extended, unlike the molecular emission. The best-fit continuum model of GSS30-IRS1 with the physical structure including flared disk, envelope, and outflow shows that the internal luminosity is 10 Lsun, and the region is externally heated by a radiation field enhanced by a factor of 130 compared to the standard local interstellar radiation field.


title: Bowling ball representation of virtual string links

abstract: In this paper we investigate the virtual string links via a probabilistic interpretation. This representation can be used to distinguish some virtual string links from classical string links. In order to study the algebraic structure behind this probabilistic interpretation we introduce the notion of virtual flat biquandle. The cocycle invariants associated with virtual flat biquandle is discussed.


title: Electrical injection and detection of spin-polarized currents in topological insulator Bi2Te2Se

abstract: Topological insulators (TIs) are an unusual phase of quantum matter with nontrivial spin-momentum locked topological surface states (TSS). The electrical detection of spin-momentum locking of the TSS in 3D TIs has been lacking till very recently. Many of the results are measured on samples with significant bulk conduction, such as metallic Bi2Se3, where it can be challenging to separate the surface and bulk contribution to the measured spin signal. Here, we report spin potentiometric measurements in thin flakes exfoliated from bulk insulating 3D TI Bi2Te2Se (BTS221) crystals, using two outside nonmagnetic (Au) contacts for driving a DC spin helical current and a middle ferromagnetic (FM)-Al2O3 tunneling contact for detecting spin polarization. The voltage measured by the FM electrode exhibits a hysteretic step-like change when sweeping an in-plane magnetic field between opposite directions along the easy axis of the FM contact to switch its magnetization. Importantly, the direction of this step-like voltage change can be reversed by reversing the direction of the DC current, and the amplitude of the change as measured by the difference in the detector voltage between opposite FM magnetization increases linearly with increasing bias current, consistent with the current-induced spin polarization of spin-momentum-locked TSS. Our work directly demonstrates the electrical injection and detection of spin polarization in TI and may enable utilization of spin-helical TSS for future applications in nanoelectronics and spintronics.


title: Resolving the Merging Planck Cluster PLCK G147.3-16.6 with GISMO

abstract: The $Planck$ satellite has recently completed an all-sky galaxy cluster survey exploiting the thermal Sunyaev-Zel'dovich (SZ) effect to locate some of the most massive systems observable. With a median redshift of $\left<z\right>=0.22$, the clusters found by $Planck$ at $z>0.3$ are proving to be exceptionally massive and/or disturbed systems. One notable $Planck$ discovery at $z=0.645$, PLCK G147.3-16.6, has a dual core and hosts a radio halo, indicating it is likely in the process of merging. We present a 16."5 resolution SZ observation of this high-$z$ merger using the Goddard-IRAM Superconducting 2 Millimeter Observer (GISMO), and compare it to X-ray follow-up observations with XMM-$Newton$. We find the SZ pressure substructure is offset from the core components seen in X-ray. We interpret this offset as possible line of sight temperature or density substructure due to the on-going merger.


title: Theory of angular dispersive imaging hard x-ray spectrographs

abstract: A spectrograph is an optical instrument that disperses photons of different energies into distinct directions and space locations, and images photon spectra on a position-sensitive detector. Spectrographs consist of collimating, angular dispersive, and focusing optical elements. Bragg reflecting crystals arranged in an asymmetric scattering geometry are used as the dispersing elements. A ray-transfer matrix technique is applied to propagate x-rays through the optical elements. Several optical designs of hard x-ray spectrographs are proposed and their performance is analyzed. Spectrographs with an energy resolution of 0.1 meV and a spectral window of imaging up to a few tens of meVs are shown to be feasible for inelastic x-ray scattering (IXS) spectroscopy applications. In another example, a spectrograph with a 1-meV spectral resolution and 85-meV spectral window of imaging is considered for Cu K-edge resonant IXS (RIXS).


title: Generalized quasiisometry on smooth riemannian manifolds

abstract: We study the boundary behavior of finitely Bi-Lipschitz mappings on smooth Riemannian manifolds


title: Summary of Long-Baseline Systematics Session at CETUP*2014

abstract: A session studying systematics in long-baseline neutrino oscillation physics was held July 14-18, 2014 as part of CETUP* 2014. Systematic effects from flux normalization and modeling, modeling of cross sections and nuclear interactions, and far detector effects were addressed. Experts presented the capabilities of existing and planned tools. A program of study to determine estimates of and requirements for the size of these effects was designed. This document summarizes the results of the CETUP* systematics workshop and the current status of systematic uncertainty studies in long-baseline neutrino oscillation measurements.


title: Some Estimates Regarding Integrated density of States for Random Schr\"{o}dinger Operator with decaying Random Potentials

abstract: We investigate some bounds for the density of states in the pure point regime for the random Schr\"{o}dinger operators $H^{\omega}=-\Delta+\displaystyle\sum_{n\in\mathbb{Z}^d}a_nq_n(\omega)$, acting on $\ell^2(\mathbb{Z}^d)$, where $\{q_n\}$ are iid random variables and $a_n\simeq|n|^{-\alpha}~~\alpha>0$.


title: Electron Thermionic Emission from Graphene and a Thermionic Energy Converter

abstract: In this paper, we propose a model to investigate the electron thermionic emission from a single-layer graphene (ignoring the effects of substrate) and to explore its application as the emitter of thermionic energy convertor (TIC). An analytical formula has been derived, which is a function of temperature, work function and Fermi energy level. The formula is significantly different from the traditional Richardson-Dushman (RD) law for which it is independent of mass to account for the supply function of the electrons in the graphene behaving like massless Fermion quasiparticles. By comparing with a recent experiment [Kaili Jiang et al., Nano Research 7, 553 (2014)] measuring electron thermionic emission from a suspended single layer graphene, our model predicts that the intrinsic work function of a single-layer graphene is about 4.514 eV with a Fermi energy level of 0.083 eV. For a given work function, a new scaling of $T^{3}$ is predicted, which is different from the traditional RD scaling of $T^2$. If the work function of the graphene is lowered to 2.5 to 3 eV, and the Fermi energy level is increased to 0.8 to 0.9 eV, it is possible to design a graphene cathode based TIC operating at around 900 K or lower, as compared with the metal-based cathode TIC (operating at about 1500 K). With a graphene based cathode (work function = 4.514 eV) at 900 K, and a metallic based anode (work function = 2.5 eV) like LaB$_6$ at 425 K, the efficiency of our proposed-TIC is about 45$\%$


title: Optimal Energy Allocation Policy for Wireless Networks in the Sky

abstract: Google's Project Loon was launched in 2013 with the aim of providing Internet access to rural and remote areas. In the Loon network, balloons travel around the Earth and bring access points to the users who cannot connect directly to the global wired Internet. The signals from the users will be transmitted through the balloon network to the base stations connected to the Internet service provider (ISP) on Earth. The process of transmitting and receiving data consume a certain amount of energy from the balloon, while the energy on balloons cannot be supplied by stable power source or by replacing batteries frequently. Instead, the balloons can harvest energy from natural energy sources, e.g., solar energy, or from radio frequency energy by equipping with appropriate circuits. However, such kinds of energy sources are often dynamic and thus how to use this energy efficiently is the main goal of this paper. In this paper, we study the optimal energy allocation problem for the balloons such that network performance is optimized and the revenue for service providers is maximized. We first formulate the stochastic optimization problem as a Markov decision process and then apply a learning algorithm based on simulation-based method to obtain optimal policies for the balloons. Numerical results obtained by extensive simulations clearly show the efficiency and convergence of the proposed learning algorithm.


title: Charged Lepton Flavor Violation in the Semi-Constrained NMSSM with Right-Handed Neutrinos

abstract: We study the \mu \to e \gamma decay in the Z_3-invariant next-to-minimal supersymmetric (SUSY) Standard Model (NMSSM) with superheavy right-handed neutrinos. We assume that the soft SUSY breaking parameters are generated at the GUT scale, not universally as in the minimal supergravity scenario but in such a way that those soft parameters which are specific to the NMSSM can differ from the soft parameters which involve only the MSSM fields while keeping the universality at the GUT scale within the soft parameters for the MSSM and right-handed neutrino fields. We call this type of boundary conditions "semi-constrained". In this model, the lepton-flavor-violating off-diagonal elements of the slepton mass matrix are induced by radiative corrections from the neutrino Yukawa couplings, just like as in the MSSM extended with the right-handed neutrinos, and these off-diagonal elements induce sizable rates of \mu \to e \gamma depending on the parameter space. Since this model has more free parameters than the MSSM, the parameter region favored from the Higgs boson mass can slightly differ from that in the MSSM. We show that there is a parameter region in which the \mu \to e \gamma decay can be observable in the near future even if the SUSY mass scale is about 4 TeV.


title: Energetics, Charge Transfer and Magnetism of Small Molecules Physisorbed on Phosphorene

abstract: First-principles calculations are performed to investigate the interaction of physisorbed small molecules, including CO, H2, H2O, NH3, NO, NO2, and O2, with phosphorene, and their energetics, charge transfer, and magnetic moment are evaluated on the basis of dispersion corrected density functional theory. Our calculations reveal that CO, H2, H2O and NH3 molecules act as a weak donor, whereas O2 and NO2 act as a strong acceptor. While NO molecule donates electrons to graphene, it receives electrons from phosphorene. Among all the investigated molecules, NO2 has the strongest interaction through hybridizing its frontier orbitals with the 3p orbital of phosphorene. The nontrivial and distinct charge transfer occurring between phosphorene and these physisorbed molecules not only renders phosphorene promising for application as a gas sensor, but also provides an effective route to modulating the polarity and density of carriers in phosphorene. In addition, the binding energy of H2 on phosphorene is found to be 0.13 eV/H2, indicating that phosphorene is suitable for both stable room-temperature hydrogen storage and its subsequent facile release.


title: Error Correcting Index Codes and Matroids

abstract: The connection between index coding and matroid theory have been well studied in the recent past. El Rouayheb et al. established a connection between multilinear representation of matroids and wireless index coding. Muralidharan and Rajan showed that a vector linear solution to an index coding problem exists if and only if there exists a representable discrete polymatroid satisfying certain conditions. Recently index coding with erroneous transmission was considered by Dau et al.. Error correcting index codes in which all receivers are able to correct a fixed number of errors was studied. In this paper we consider a more general scenario in which each receiver is able to correct a desired number of errors, calling such index codes differential error correcting index codes. A link between differential error correcting index codes and certain matroids is established. We define matroidal differential error correcting index codes and we show that a scalar linear differential error correcting index code exists if and only if it is matroidal differential error correcting index code associated with a representable matroid.


title: Symmetry and Explicit Marking of the Critical Phase in Two-Bath Spin-Boson Model

abstract: The spin-boson model is a paradigm for studying decoherence, relaxation, entanglement and other effects that arise in a quantum system coupled to environmental degrees of freedom. At zero temperature, a localization-delocalization phase transition is known to exist in the sub-Ohmic regime, where the standard density matrix renormalization group algorithm is inadequate due to the divergence in the number of low-frequency modes. This limitation is circumvented in this work by symmetrically optimizing the phonon basis and introducing an order parameter accounting for the U(1) symmetry for a two-bath spin-boson model, by which we are able to determine the classification and criticality of the phase transition explicitly. Compared with variational results, the critical phase is characterized by spontaneous vanishing of boson displacements in both the baths, resulting in an accurate phase diagram with three model parameters.


title: The focusing cubic NLS on exterior domains in three dimensions

abstract: We consider the focusing cubic NLS in the exterior $\Omega$ of a smooth, compact, strictly convex obstacle in three dimensions. We prove that the threshold for global existence and scattering is the same as for the problem posed on Euclidean space. Specifically, we prove that if $E(u_0)M(u_0)<E(Q)M(Q)$ and $\|\nabla u_0\|_2\|u_0\|_2<\|\nabla Q\|_2\|Q\|_2$, the corresponding solution to the initial-value problem with Dirichlet boundary conditions exists globally and scatters to linear evolutions asymptotically in the future and in the past. Here, $Q(x)$ denotes the ground state for the focusing cubic NLS in $\mathbb{R}^3$.


title: Quantum walk and Anderson localization of rotational excitations in disordered ensembles of polar molecules

abstract: We consider the dynamics of rotational excitations placed on a single molecule in spatially disordered 1D, 2D and 3D ensembles of ultracold molecules trapped in optical lattices. The disorder arises from incomplete populations of optical lattices with molecules. This leads to a model corresponding to a quantum particle with long-range tunnelling amplitudes moving on a lattice with the same on-site energy but with forbidden access to random sites (vacancies). We examine the time and length scales of Anderson localization for this type of disorder with realistic experimental parameters in the Hamiltonian. We show that for an experimentally realized system of KRb molecules on an optical lattice this type of disorder leads to disorder-induced localization in 1D and 2D systems on a time scale $t \sim 1$ sec. For 3D lattices with $55$ sites in each dimension and vacancy concentration $ 90~\%$, the rotational excitations diffuse to the edges of the lattice and show no signature of Anderson localization. We examine the role of the long-range tunnelling amplitudes allowing for transfer of rotational excitations between distant lattice sites. Our results show that the long-range tunnelling has little impact on the dynamics in the diffusive regime but affects significantly the localization dynamics in lattices with large concentrations of vacancies, enhancing the width of the localized distributions in 2D lattices by more than a factor of 2. Our results raise a general question whether quantum particles with long-range tunnelling can undergo quantum localization in 3D lattices with substitutional disorder.


title: Determination of the stretch tensor for structural transformations

abstract: The transformation stretch tensor plays an essential role in the evaluation of conditions of compatibility between phases and the use of the Cauchy-Born rule. This tensor is difficult to measure directly from experiment. We give an algorithm for the determination of the transformation stretch tensor from x-ray measurements of structure and lattice parameters. When evaluated on some traditional and emerging phase transformations the algorithm gives unexpected results.


title: Operator-Valued Tensors on Manifolds: A Framework for Field Quantization

abstract: In this paper we try to prepare a framework for field quantization. To this end, we aim to replace the field of scalars R by self-adjoint elements of a commutative C-algebra, and reach an appropriate generalization of geometrical concepts on manifolds. First, we put forward the concept of operator-valued tensors and extend semi-Riemannian metrics to operator valued metrics. Then, in this new geometry, some essential concepts of Riemannian geometry such as curvature ten- sor, Levi-Civita connection, Hodge star operator, exterior derivative, divergence,... will be considered.


title: A Property of Random Walks on a Cycle Graph

abstract: We analyze the Hunter vs Rabbit game on graph, which is a kind of model of communication in an adhoc mobile network. Let $G$ be a cycle graph with $N$ nodes. The hunter can move from a vertex to another vertex on the graph along an edge. The rabbit can move to any vertex on graph at once. We formalized the game using the random walk framework. The strategy of the rabbit is formalized using a one dimensional random walk over $\mathbb{Z}$. We classify strategies using the order $O(k^{-\beta-1})$ of their Fourier transformation. We investigate lower bounds and upper bounds of a probability that the hunter catches the rabbit. We found a constant lower bound if $\beta \in (0,1)$. That is there is not depend on the size $N$ of the graph. We show the order is equivalent to $O(1/\log N)$ if $\beta=1$ and a lower bound is $1/N^{(\beta-1)/\beta}$ if $\beta \in (1,2]$. Those results assist to choose the parameter $\beta$ of a rabbit strategy according to the size $N$ of the given graph. We introduce a formalization of strategies using a random walk, theoretical estimation of bounds of a probability that the hunter catches the rabbit, and we also show computing simulation results.


title: Visualizing the elongated vortices in $\gamma$-Ga nanostrips

abstract: We study the magnetic response of superconducting $\gamma$-Ga via low temperature scanning tunneling microscopy and spectroscopy. The magnetic vortex cores rely substantially on the Ga geometry, and exhibit an unexpectedly-large axial elongation with aspect ratio up to 40 in rectangular Ga nano-strips (width $l$ $<$ 100 nm). This is in stark contrast with the isotropic circular vortex core in a larger round-shaped Ga island. We suggest that the unusual elongated vortices in Ga nanostrips originate from geometric confinement effect probably via the strong repulsive interaction between the vortices and Meissner screening currents at the sample edge. Our finding provides novel conceptual insights into the geometrical confinement effect on magnetic vortices and forms the basis for the technological applications of superconductors.


title: Difficulties applying recent blind source separation techniques to EEG and MEG

abstract: High temporal resolution measurements of human brain activity can be performed by recording the electric potentials on the scalp surface (electroencephalography, EEG), or by recording the magnetic fields near the surface of the head (magnetoencephalography, MEG). The analysis of the data is problematic due to the fact that multiple neural generators may be simultaneously active and the potentials and magnetic fields from these sources are superimposed on the detectors. It is highly desirable to un-mix the data into signals representing the behaviors of the original individual generators. This general problem is called blind source separation and several recent techniques utilizing maximum entropy, minimum mutual information, and maximum likelihood estimation have been applied. These techniques have had much success in separating signals such as natural sounds or speech, but appear to be ineffective when applied to EEG or MEG signals. Many of these techniques implicitly assume that the source distributions have a large kurtosis, whereas an analysis of EEG/MEG signals reveals that the distributions are multimodal. This suggests that more effective separation techniques could be designed for EEG and MEG signals.


title: Convergent Bayesian formulations of blind source separation and electromagnetic source estimation

abstract: We consider two areas of research that have been developing in parallel over the last decade: blind source separation (BSS) and electromagnetic source estimation (ESE). BSS deals with the recovery of source signals when only mixtures of signals can be obtained from an array of detectors and the only prior knowledge consists of some information about the nature of the source signals. On the other hand, ESE utilizes knowledge of the electromagnetic forward problem to assign source signals to their respective generators, while information about the signals themselves is typically ignored. We demonstrate that these two techniques can be derived from the same starting point using the Bayesian formalism. This suggests a means by which new algorithms can be developed that utilize as much relevant information as possible. We also briefly mention some preliminary work that supports the value of integrating information used by these two techniques and review the kinds of information that may be useful in addressing the ESE problem.


title: New trigonometric and hyperbolic inequalities

abstract: The aim of this paper is to prove new trigonometric and hyperbolic inequalities, which constitute among others refinements or analogs of famous Cusa-Huygens, Wu-Srivastava, and related inequalities. In most cases, the obtained results are sharp.


title: Non-probabilistic odds and forecasting with imperfect models

abstract: Probability forecasts are intended to account for the uncertainties inherent in forecasting. It is suggested that from an end-user's point of view probability is not necessarily sufficient to reflect uncertainties that are not simply the result of complexity or randomness, for example, probability forecasts may not adequately account for uncertainties due to model error. It is suggested that an alternative forecast product is to issue non-probabilistic odds forecasts, which may be as useful to end-users, and give a less distorted account of the uncertainties of a forecast. Our analysis of odds forecasts derives from game theory using the principle that if forecasters truly believe their forecasts, then they should take bets at the odds they offer and not expect to be bankrupted. Despite this game theoretic approach, it is not a market or economic evaluation; it is intended to be a scientific evaluation. Illustrative examples are given of the calculation of odds forecasts and their application to investment, loss mitigation and ensemble weather forecasting.


title: Estimation of Reliability in the Two-Parameter Geometric Distribution

abstract: In this article, the reliabilities $R(t)=P(X\geq t)$, when $X$ follows two-parameter geometric distribution and $R=P(X\leq Y)$, arises under stress-strength setup, when X and Y assumed to follow two-parameter geometric independently have been found out. Maximum Likelihood Estimator (MLE) and an Unbiased Estimator (UE) of these have been derived. MLE and UE of the reliability of k-out-of-m system have also been derived. The estimators have been compared through simulation study.


title: Simple model for the Darwinian transition in early evolution

abstract: It has been hypothesized that in the era just before the last universal common ancestor emerged, life on earth was fundamentally collective. Ancient life forms shared their genetic material freely through massive horizontal gene transfer (HGT). At a certain point, however, life made a transition to the modern era of individuality and vertical descent. Here we present a minimal model for this hypothesized "Darwinian transition." The model suggests that HGT-dominated dynamics may have been intermittently interrupted by selection-driven processes during which genotypes became fitter and decreased their inclination toward HGT. Stochastic switching in the population dynamics with three-point (hypernetwork) interactions may have destabilized the HGT-dominated collective state and led to the emergence of vertical descent and the first well-defined species in early evolution. A nonlinear analysis of a stochastic model dynamics covering key features of evolutionary processes (such as selection, mutation, drift and HGT) supports this view. Our findings thus suggest a viable route from early collective evolution to the start of individuality and vertical Darwinian evolution, enabling the emergence of the first species.

categories: q-bio.PE cond-mat.stat-mech nlin.AO physics.bio-ph


title: Duality and Modularity in Elliptic Integrable Systems and Vacua of N=1* Gauge Theories

abstract: We study complexified elliptic Calogero-Moser integrable systems. We determine the value of the potential at isolated extrema, as a function of the modular parameter of the torus on which the integrable system lives. We calculate the extrema for low rank B,C,D root systems using a mix of analytical and numerical tools. For so(5) we find convincing evidence that the extrema constitute a vector valued modular form for a congruence subgroup of the modular group. For so(7) and so(8), the extrema split into two sets. One set contains extrema that make up vector valued modular forms for congruence subgroups, and a second set contains extrema that exhibit monodromies around points in the interior of the fundamental domain. The former set can be described analytically, while for the latter, we provide an analytic value for the point of monodromy for so(8), as well as extensive numerical predictions for the Fourier coefficients of the extrema. Our results on the extrema provide a rationale for integrality properties observed in integrable models, and embed these into the theory of vector valued modular forms. Moreover, using the data we gather on the modularity of complexified integrable system extrema, we analyse the massive vacua of mass deformed N=4 supersymmetric Yang-Mills theories with low rank gauge group of type B,C and D. We map out their transformation properties under the infrared electric-magnetic duality group as well as under triality for N=1* with gauge algebra so(8). We find several intriguing properties of the quantum gauge theories.


title: Calculations relating to some special Harmonic numbers

abstract: We report on the results of a computer search for primes $p$ which divide an Harmonic number $H_{\lfloor p/N \rfloor}$ with small $N > 1$.


title: Statistical properties of the localization measure in a finite-dimensional model of the quantum kicked rotator

abstract: We study the quantum kicked rotator in the classically fully chaotic regime $K=10$ and for various values of the quantum parameter $k$ using Izrailev's $N$-dimensional model for various $N \le 3000$, which in the limit $N \rightarrow \infty$ tends to the exact quantized kicked rotator. By numerically calculating the eigenfunctions in the basis of the angular momentum we find that the localization length ${\cal L}$ for fixed parameter values has a certain distribution, in fact its inverse is Gaussian distributed, in analogy and in connection with the distribution of finite time Lyapunov exponents of Hamilton systems. However, unlike the case of the finite time Lyapunov exponents, this distribution is found to be independent of $N$, and thus survives the limit $N=\infty$. This is different from the tight-binding model of Anderson localization. The reason is that the finite bandwidth approximation of the underlying Hamilton dynamical system in the Shepelyansky picture (D.L. Shepelyansky, {\em Phys. Rev. Lett.} {\bf 56}, 677 (1986)) does not apply rigorously. This observation explains the strong fluctuations in the scaling laws of the kicked rotator, such as e.g. the entropy localization measure as a function of the scaling parameter $\Lambda={\cal L}/N$, where $\cal L$ is the theoretical value of the localization length in the semiclassical approximation. These results call for a more refined theory of the localization length in the quantum kicked rotator and in similar Floquet systems, where we must predict not only the mean value of the inverse of the localization length $\cal L$ but also its (Gaussian) distribution, in particular the variance. In order to complete our studies we numerically analyze the related behavior of finite time Lyapunov exponents in the standard map and of the 2$\times$2 transfer matrix formalism. This paper is extending our recent work.


title: Planar Markovian Holonomy Fields

abstract: We study planar random holonomy fields which are processes indexed by paths on the plane which behave well under the concatenation and orientation-reversing operations on paths. We define the Planar Markovian Holonomy Fields as planar random holonomy fields which satisfy some independence and invariance by area-preserving homeomorphisms properties. We use the theory of braids in the framework of classical probabilities: for finite and infinite random sequences the notion of invariance by braids is defined and we prove a new version of the de-Finetti's Theorem. This allows us to construct a family of Planar Markovian Holonomy Fields, the Yang-Mills fields, and we prove that any regular Planar Markovian Holonomy Field is a planar Yang-Mills field. This family of planar Yang-Mills fields can be partitioned into three categories according to the degree of symmetry: we study some equivalent conditions in order to classify them. Finally, we recall the notion of Markovian Holonomy Fields and construct a bridge between the planar and non-planar theories. Using the results previously proved in the article, we compute, for any Markovian Holonomy Field, the "law" of any family of contractible loops drawn on a surface.


title: Renormalization procedure for random tensor networks and the canonical tensor model

abstract: We discuss a renormalization procedure for random tensor networks, and show that the corresponding renormalization-group flow is given by the Hamiltonian vector flow of the canonical tensor model, which is a discretized model of quantum gravity. The result is the generalization of the previous one concerning the relation between the Ising model on random networks and the canonical tensor model with N=2. We also prove a general theorem which relates discontinuity of the renormalization-group flow and the phase transitions of random tensor networks.


title: European VLBI Network: Present and Future

abstract: The European VLBI Network is a collaboration of the major radio astronomical institutes in Europe, Asia, South Africa and Puerto Rico. Established four decades ago, since then it has constantly improved its performance in terms made using resolution, data bit-rate and image fidelity with improvements in performance, and the addition of new stations and observing capabilities. The EVN provides open skies access and has over time become a common-user facility. In this contribution we discuss the present status and perspectives for the array in a continuously changing environment, especially in the era of ALMA and with the Square Kilometre Array ante portas.


title: Enabling High-Level Application Development for the Internet of Things

abstract: Application development in the Internet of Things (IoT) is challenging because it involves dealing with a wide range of related issues such as lack of separation of concerns, and lack of high-level of abstractions to address both the large scale and heterogeneity. Moreover, stakeholders involved in the application development have to address issues that can be attributed to different life-cycles phases. when developing applications. First, the application logic has to be analyzed and then separated into a set of distributed tasks for an underlying network. Then, the tasks have to be implemented for the specific hardware. Apart from handling these issues, they have to deal with other aspects of life-cycle such as changes in application requirements and deployed devices. Several approaches have been proposed in the closely related fields of wireless sensor network, ubiquitous and pervasive computing, and software engineering in general to address the above challenges. However, existing approaches only cover limited subsets of the above mentioned challenges when applied to the IoT. This paper proposes an integrated approach for addressing the above mentioned challenges. The main contributions of this paper are: (1) a development methodology that separates IoT application development into different concerns and provides a conceptual framework to develop an application, (2) a development framework that implements the development methodology to support actions of stakeholders. The development framework provides a set of modeling languages to specify each development concern and abstracts the scale and heterogeneity related complexity. It integrates code generation, task-mapping, and linking techniques to provide automation. Code generation supports the application development phase by producing a programming framework that allows stakeholders to focus on the application logic, while our mapping and linking techniques together support the deployment phase by producing device-specific code to result in a distributed system collaboratively hosted by individual devices. Our evaluation based on two realistic scenarios shows that the use of our approach improves the productivity of stakeholders involved in the application development.


title: Improving the Description of Nonmagnetic and Magnetic Molecular Crystals via the van der Waals Density Functional

abstract: We have derived and implemented a stress tensor formulation for the van derWaals density functional (vdW-DF) with spin-polarization-dependent gradient correction (GC) recently proposed by the authors [J. Phys. Soc. Jpn. 82, 093701 (2013)] and applied it to nonmagnetic and magnetic molecular crystals under ambient condition. We found that the cell parameters of the molecular crystals obtained with vdW-DF show an overall improvement compared with those obtained using local density and generalized gradient approximations. In particular, the original vdW-DF with GC gives the equilibrium structural parameters of solid oxygen in the {\alpha}-phase, which are in good agreement with the experiment.


title: Entropy and drift in word hyperbolic groups

abstract: The fundamental inequality of Guivarc'h relates the entropy and the drift of random walks on groups. It is strict if and only if the random walk does not behave like the uniform measure on balls. We prove that, in any nonelementary hyperbolic group which is not virtually free, endowed with a word distance, the fundamental inequality is strict for symmetric measures with finite support, uniformly for measures with a given support. This answers a conjecture of S. Lalley. For admissible measures, this is proved using previous results of Ancona and Blach{\`e}re-Ha{\"i}ssinsky-Mathieu. For non-admissible measures, this follows from a counting result, interesting in its own right: we show that, in any infinite index subgroup, the number of non-distorted points is exponentially small. The uniformity is obtained by studying the behavior of measures that degenerate towards a measure supported on an elementary subgroup.


title: Certifying isolated singular points and their multiplicity structure

abstract: This paper presents two new constructions related to singular solutions of polynomial systems. The first is a new deflation method for an isolated singular root. This construc-tion uses a single linear differential form defined from the Jacobian matrix of the input, and defines the deflated system by applying this differential form to the original system. The advantages of this new deflation is that it does not introduce new variables and the increase in the number of equations is linear instead of the quadratic increase of previous methods. The second construction gives the coefficients of the so-called inverse system or dual basis, which defines the multiplicity structure at the singular root. We present a system of equations in the original variables plus a relatively small number of new vari-ables. We show that the roots of this new system include the original singular root but now with multiplicity one, and the new variables uniquely determine the multiplicity structure. Both constructions are "exact", meaning that they permit one to treat all conjugate roots simultaneously and can be used in certification procedures for singular roots and their multiplicity structure with respect to an exact rational polynomial system.


title: Strong $\Sigma_bNB$ and $\Sigma_cND$ coupling constants in QCD

abstract: We study the strong interactions among the heavy bottom spin--1/2 $\Sigma_b$ baryon, nucleon and $B$ meson as well as the heavy charmed spin--1/2 $\Sigma_c$ baryon, nucleon and $D$ meson in the context of QCD sum rules. We calculate the corresponding strong coupling form factors defining these vertices by using a three point correlation function. We obtain the numerical values of the corresponding strong coupling constants via the most prominent structure entering the calculations.


title: A Note on a Question of Erd\H{o}s & Graham

abstract: Erd\H{o}s & Graham ask whether the equation $x^2+y^2=z^2$ is partition regular, i.e. whether it has a finite Rado number. This note provides a lower bound and also states results in the affirmative for two similar quadratic equations.


title: Multi-resource Energy-efficient Routing in Cloud Data Centers with Networks-as-a-Service

abstract: With the rapid development of software defined networking and network function virtualization, researchers have proposed a new cloud networking model called Network-as-a-Service (NaaS) which enables both in-network packet processing and application-specific network control. In this paper, we revisit the problem of achieving network energy efficiency in data centers and identify some new optimization challenges under the NaaS model. Particularly, we extend the energy-efficient routing optimization from single-resource to multi-resource settings. We characterize the problem through a detailed model and provide a formal problem definition. Due to the high complexity of direct solutions, we propose a greedy routing scheme to approximate the optimum, where flows are selected progressively to exhaust residual capacities of active nodes, and routing paths are assigned based on the distributions of both node residual capacities and flow demands. By leveraging the structural regularity of data center networks, we also provide a fast topology-aware heuristic method based on hierarchically solving a series of vector bin packing instances. Our simulations show that the proposed routing scheme can achieve significant gain on energy savings and the topology-aware heuristic can produce comparably good results while reducing the computation time to a large extent.


title: The 2012 flare of PG 1553+113 seen with H.E.S.S. and Fermi-LAT

abstract: Very high energy (VHE, $E>$100 GeV) $\gamma$-ray flaring activity of the high-frequency peaked BL Lac object \pg\ has been detected by the \hess\ telescopes. The flux of the source increased by a factor of 3 during the nights of 2012 April 26 and 27 with respect to the archival measurements with hint of intra-night variability. No counterpart of this event has been detected in the \fla\ data. This pattern is consistent with VHE $\gamma$ ray flaring being caused by the injection of ultrarelativistic particles, emitting $\gamma$ rays at the highest energies. The dataset offers a unique opportunity to constrain the redshift of this source at \bestz\ using a novel method based on Bayesian statistics. The indication of intra-night variability is used to introduce a novel method to probe for a possible Lorentz Invariance Violation (LIV), and to set limits on the energy scale at which Quantum Gravity (QG) effects causing LIV may arise. For the subluminal case, the derived limits are $\textrm{E}_{\rm QG,1}>4.10\times 10^{17}$ GeV and $\textrm{E}_{\rm QG,2}>2.10\times 10^{10}$ GeV for linear and quadratic LIV effects, respectively.


title: GRB Spectra in the complex of synchrotron and Compton processes

abstract: Under the steady state condition, the spectrum of electrons is investigated by solving the continuity equation under the complex radiation of both the synchrotron and Compton processes. The resulted GRB spectrum is a broken power law in both the fast and slow cooling phases. On the basis of this electron spectrum, the spectral indices of the Band function in four different phases are presented. In the complex radiation frame, the detail investigation on physical parameters reveals that both the reverse shock photosphere model and the forward shock with strong coupling model can answer the $\alpha \sim -1$ problem. A possible marginal to fast cooling phase transition in GRB 080916C is discussed. The time resolved spectra in different pulses of GRB 100724B, GRB 100826A and GRB 130606B are investigated. We found that the flux is proportional to the peak energy in almost all pulses. The phases for different pulses are determined according to the spectral index revolution. We found the strong correlations between spectral indices and the peak energy in GRB 100826A, which can be explained by the Compton effect in the fast cooling phase. However, the complex scenario predicts a steeper index for the injected electrons, which challenges the acceleration mechanism in GRBs.


title: Walk-powers and homomorphism bound of planar graphs

abstract: As an extension of the Four-Color Theorem it is conjectured that every planar graph of odd-girth at least $2k+1$ admits a homomorphism to $PC_{2k}=(\mathbb{Z}_2^{2k}, \{e_1, e_2, ...,e_{2k}, J\})$ where $e_i$'s are standard basis and $J$ is all 1 vector. Noting that $PC_{2k}$ itself is of odd-girth $2k+1$, in this work we show that if the conjecture is true, then $PC_{2k}$ is an optimal such a graph both with respect to number of vertices and number of edges. The result is obtained using the notion of walk-power of graphs and their clique numbers. An analogous result is proved for bipartite signed planar graphs of unbalanced-girth $2k$. The work is presented on a uniform frame work of planar consistent signed graphs.


title: Many-Body electronic structure calculations of Eu doped ZnO

abstract: The formation energies and electronic structure of europium doped zinc oxide has been determined using DFT and many-body $GW$ methods. In the absence of intrisic defects we find that the europium-$f$ states are located in the ZnO band gap with europium possessing a formal charge of 2+. On the other hand, the presence of intrinsic defects in ZnO allows intraband $f-f$ transitions otherwise forbidden in atomic europium. This result coorroborates with recently observed photoluminescence in the visible red region [1].


title: Mapping the $q$-voter model: From a single chain to complex networks

abstract: We propose and compare six different ways of mapping the modified $q$-voter model to complex networks. Considering square lattices, Barab\'asi-Albert, Watts-Strogatz and real Twitter networks, we ask the question if always a particular choice of the group of influence of a fixed size $q$ leads to different behavior at the macroscopic level. Using Monte Carlo simulations we show that the answer depends on the relative average path length of the network and for real-life topologies the differences between the considered mappings may be negligible.


title: Uniqueness of p(f) and P[f]

abstract: Let f be a non constant meromorphic function and a(not identically zero or infinity) be a meromorphic function satisfying T(r,a) = o(T(r,f)) as r tends to infinity, and p(z) be a polynomial of degree n greater than or equal to 1 with p(0) = 0. Let P[f] be a non constant differential polynomial of f. Under certain essential conditions, we prove the uniqueness of p(f) and P[f] when p(f) and P[f] share a with weight l greater than or equal to zero. Our result generalizes the results due to Zang and Lu, Banerjee and Majumder, Bhoosnurmath and Kabbur and answers a question of Zang and Lu.


title: Analysis of the heavy mesons in the nuclear matter with the QCD sum rules

abstract: In this article, I calculate the contributions of the nuclear matter induced condensates up to dimension 5, take into account the next-to-leading order contributions of the nuclear matter induced quark condensate, study the properties of the scalar, pseudoscalar, vector and axial-vector heavy mesons in the nuclear matter with the QCD sum rules in a systematic way, and obtain the shifts of the masses and decay constants. Furthermore, I study the heavy-meson-nucleon scattering lengths as a byproduct, and obtain the conclusion qualitatively about the possible existence of heavy-meson-nucleon bound states.


title: Orbifold construction of holomorphic vertex operator algebras associated to inner automorphisms

abstract: In this article, we construct three new holomorphic vertex operator algebras of central charge $24$ using the $\mathbb{Z}_2$-orbifold construction associated to inner automorphisms. Their weight one subspaces has the Lie algebra structures $D_{7,3}A_{3,1}G_{2,1}$, $E_{7,3}A_{5,1}$, and $A_{8,3}A_{2,1}^2$. In addition, we discuss the constructions of holomorphic vertex operator algebras with Lie algebras $A_{5,6}C_{2,3}A_{1,2}$ and $D_{6,5}A_{1,1}^2$ from holomorphic vertex operator algebras with Lie algebras $C_{5,3}G_{2,2}A_{1,1}$ and $A_{4,5}^2$, respectively.


title: Period integrals and mutation

abstract: Let $f$ be a Laurent polynomial in two variables, whose Newton polygon strictly contains the origin and whose vertices are primitive lattice points, and let $L_f$ be the minimal-order differential operator that annihilates the period integral of $f$. We prove several results about $f$ and $L_f$ in terms of the Newton polygon of $f$ and the combinatorial operation of *mutation*, in particular we give an in principle complete description of the monodromy of $L_f$ around the origin. Special attention is given to the class of *maximally mutable* Laurent polynomials, which has applications to the conjectured classification of Fano manifolds via mirror symmetry.


title: Experimental realisation of generalised qubit measurements based on quantum walks

abstract: We report an experimental implementation of a single-qubit generalised measurement scenario(POVM) based on a quantum walk model. The qubit is encoded in a single-photon polarisation. The photon performs a quantum walk on an array of optical elements, where the polarisation-dependent translation is performed via birefringent beam displacers and a change of the polarisation is implemented with the help of wave-plates. We implement: (i) Trine-POVM, i.e., the POVM elements uniformly distributed on an equatorial plane of the Bloch sphere; (ii) Symmetric-Informationally- Complete (SIC) POVM; and (iii) Unambiguous Discrimination of two non-orthogonal qubit states.


title: Discrete-Time Models for Implicit Port-Hamiltonian Systems

abstract: Implicit representations of finite-dimensional port-Hamiltonian systems are studied from the perspective of their use in numerical simulation and control design. Implicit representations arise when a system is modeled in Cartesian coordinates and when the system constraints are applied in the form of additional algebraic equations (the system model is in a DAE form). Such representations lend themselves better to sample-data approximations. An implicit representation of a port-Hamiltonian system is given and it is shown how to construct a sampled-data model that preserves the port-Hamiltonian structure under sample and hold.


title: Better Answers to Real Questions

abstract: We consider existential problems over the reals. Extended quantifier elimination generalizes the concept of regular quantifier elimination by providing in addition answers, which are descriptions of possible assignments for the quantified variables. Implementations of extended quantifier elimination via virtual substitution have been successfully applied to various problems in science and engineering. So far, the answers produced by these implementations included infinitesimal and infinite numbers, which are hard to interpret in practice. We introduce here a post-processing procedure to convert, for fixed parameters, all answers into standard real numbers. The relevance of our procedure is demonstrated by application of our implementation to various examples from the literature, where it significantly improves the quality of the results.


title: Superradiant cascade emissions in an atomic ensemble via four-wave mixing

abstract: We investigate superradiant cascade emissions from an atomic ensemble driven by two-color classical fields. The correlated pair of photons (signal and idler) is generated by adiabatically driving the system with large-detuned light fields via four-wave mixing. The signal photon from the upper transition of the diamond-type atomic levels is followed by the idler one which can be superradiant due to light-induced dipole-dipole interactions. We then calculate the cooperative Lamb shift (CLS) of the idler photon, which is a cumulative effect of interaction energy.We study its dependence on a cylindrical geometry, a conventional setup in cold atom experiments, and estimate the maximum CLS which can be significant and observable. Manipulating the CLS of cascade emissions enables frequency qubits that provide alternative robust elements in quantum network.


title: Efficient Linear Scaling Approach for Computing the Kubo Hall Conductivity

abstract: We report an order-N approach to compute the Kubo Hall conductivity for disorderd two-dimensional systems reaching tens of millions of orbitals, and realistic values of the applied external magnetic fields (as low as a few Tesla). A time-evolution scheme is employed to evaluate the Hall conductivity $\sigma_{xy}$ using a wavepacket propagation method and a continued fraction expansion for the computation of diagonal and off-diagonal matrix elements of the Green functions. The validity of the method is demonstrated by comparison of results with brute-force diagonalization of the Kubo formula, using (disordered) graphene as system of study. This approach to mesoscopic system sizes is opening an unprecedented perspective for so-called reverse engineering in which the available experimental transport data are used to get a deeper understanding of the microscopic structure of the samples. Besides, this will not only allow addressing subtle issues in terms of resistance standardization of large scale materials (such as wafer scale polycrystalline graphene), but will also enable the discovery of new quantum transport phenomena in complex two-dimensional materials, out of reach with classical methods.


title: Direct laser acceleration of electrons in free-space

abstract: Compact laser-driven accelerators are versatile and powerful tools of unarguable relevance on societal grounds for the diverse purposes of science, health, security, and technology because they bring enormous practicality to state-of-the-art achievements of conventional radio-frequency accelerators. Current benchmarking laser-based technologies rely on a medium to assist the light-matter interaction, which impose material limitations or strongly inhomogeneous fields. The advent of few cycle ultra-intense radially polarized lasers has materialized an extensively studied novel accelerator that adopts the simplest form of laser acceleration and is unique in requiring no medium to achieve strong longitudinal energy transfer directly from laser to particle. Here we present the first observation of direct longitudinal laser acceleration of non-relativistic electrons that undergo highly-directional multi-GeV/m accelerating gradients. This demonstration opens a new frontier for direct laser-driven particle acceleration capable of creating well collimated and relativistic attosecond electron bunches and x-ray pulses.


title: Averaged controllability for random evolution partial differential equations

abstract: We analyze the averaged controllability properties of random evolution Partial Differential Equations. We mainly consider heat and Schr\"odinger equations with random parameters, although the problem is also formulated in an abstract frame. We show that the averages of parabolic equations lead to parabolic-like dynamics that enjoy the null-controllability properties of solutions of heat equations in an arbitrarily short time and from arbitrary measurable sets of positive measure. In the case of Schr\"odinger equations we show that, depending on the probability density governing the random parameter, the average may behave either as a conservative or a parabolic-like evolution, leading to controllability properties, in average, of very different kind.


title: On the large time behavior of the solutions of a nonlocal ordinary differential equation with mass conservation

abstract: We consider an initial value problem for a nonlocal differential equation with a bistable nonlinearity in several space dimensions. The equation is an ordinary differential equation with respect to the time variable t, while the nonlocal term is expressed in terms of spatial integration. We discuss the large time behavior of solutions and prove, among other things, the convergence to steady-states. The proof that the solution orbits are relatively compact is based upon the rearrangement theory.


title: Memoization for Unary Logic Programming: Characterizing PTIME

abstract: We give a characterization of deterministic polynomial time computation based on an algebraic structure called the resolution semiring, whose elements can be understood as logic programs or sets of rewriting rules over first-order terms. More precisely, we study the restriction of this framework to terms (and logic programs, rewriting rules) using only unary symbols. We prove it is complete for polynomial time computation, using an encoding of pushdown automata. We then introduce an algebraic counterpart of the memoization technique in order to show its PTIME soundness. We finally relate our approach and complexity results to complexity of logic programming. As an application of our techniques, we show a PTIME-completeness result for a class of logic programming queries which use only unary function symbols.


title: An Application of the Separator of Subsets of Semigroups in the Number Theory

abstract: In this paper we give a construction for a special type of congruences on commutative semigroups. We apply our result for the multiplicative semigroup of all positive integers.


title: Topology of strongly polar weighted homogeneous links

abstract: We consider a canonical $S^1$ action on $S^3$ which is defined by $(\rho,(z_1,z_2))\mapsto (z_1\rho^p,z_2\rho^q)$ for $\rho\in S^1$ and $(z_1,z_2)\in S^3\subset {\mathbb C}^2$. We consider a link consisting of finite orbits of this action, which some of the orbits are reversely oriented. Such a link appears as a link of a certain type of mixed polynomials. We study the space of such links and show smooth degeneration relations.


title: Simple Polyominoes are Prime

abstract: In this paper we show that polyomino ideal of a simple polyomino coincides with the toric ideal of a weakly chordal bipartite graph and hence it has a quadratic Gr\"obner basis with respect to a suitable monomial order.


title: BDgraph: An R Package for Bayesian Structure Learning in Graphical Models

abstract: Graphical models provide powerful tools to uncover complicated patterns in multivariate data and are commonly used in Bayesian statistics and machine learning. In this paper, we introduce the R package BDgraph which performs Bayesian structure learning for general undirected graphical models (decomposable and non-decomposable) with continuous, discrete, and mixed variables. The package efficiently implements recent improvements in the Bayesian literature, including that of Mohammadi and Wit (2015) and Dobra and Mohammadi (2018). To speed up computations, the computationally intensive tasks have been implemented in C++ and interfaced with R, and the package has parallel computing capabilities. In addition, the package contains several functions for simulation and visualization, as well as several multivariate datasets taken from the literature and used to describe the package capabilities. The paper includes a brief overview of the statistical methods which have been implemented in the package. The main part of the paper explains how to use the package. Furthermore, we illustrate the package's functionality in both real and artificial examples.