title: Yukawa Unification in an SO(10) SUSY GUT: SUSY on the Edge

abstract: In this paper we analyze Yukawa unification in a three family SO(10) SUSY GUT. We perform a global $\chi^2$ analysis and show that SUSY effects do not decouple even though the universal scalar mass parameter at the GUT scale, $m_{16}$, is found to lie between 15 and 30 TeV with the best fit given for $m_{16} \approx 25$ TeV. Note, SUSY effects don't decouple since stops and bottoms have mass of order 5 TeV, due to RG running from $M_{GUT}$. The model has many testable predictions. Gauginos are the lightest sparticles and the light Higgs boson is very much Standard Model-like. The model is consistent with flavor and CP observables with the $BR(\mu \to e\gamma)$ close to the experimental upper bound. With such a large value of $m_{16}$ we clearly cannot be considered "natural" SUSY nor are we "Split" SUSY. We are thus in the region in between or "SUSY on the Edge."


title: Band-9 ALMA observations of the [NII] 122 micron line and FIR continuum in two high-z galaxies

abstract: We present Atacama Large Millimeter Array (ALMA) observations of two high-redshift systems (SMMJ02399-0136 and the Cloverleaf QSO) in their rest-frame 122 micron continuum (~650 GHz or ~450 micron on-sky) and [NII] 122 micron line emission. The continuum observations with a synthesized beam of ~0.25" resolve both sources and recover the expected flux. The Cloverleaf is resolved into a partial Einstein ring, while the SMMJ02399-0136 is unambiguously separated into two components; an AGN associated point source and an extend region at the location of a previously identified dusty starburst. We detect the [NII] line in both systems, though significantly weaker than our previous detections made with the 1st generation z(Redshift) and Early Universe Spectrometer. We show that this discrepancy is mostly explained if the line flux is resolved out due to significantly more extended emission and longer ALMA baselines than expected. Based on the ALMA observations we determine that greater than 75% of the total [NII] line flux in each source is produced via star formation. We use the [NII] line flux that is recovered by ALMA to constrain the N/H abundance, ionized gas mass, hydrogen ionizing photon rate, and star formation rate. In SMMJ02399-0136 we discover it contains a significant amount (~1000 solar masses per year) of unobscured star formation in addition to its dusty starburst and argue that SMMJ02399-0136 may be similar to the Antennae Galaxies (Arp 244) locally. In total these observations provide a new look at two well-studied systems while demonstrating the power and challenges of Band-9 ALMA observations of high-z systems.


title: Spectral Confusion for Cosmological Surveys of Redshifted CII Observations

abstract: Far infrared cooling lines are ubiquitous features in the spectra of star forming galaxies. Surveys of redshifted fine-structure lines provide a promising new tool to study structure formation and galactic evolution at redshifts including the epoch of reionization as well as the peak of star formation. Unlike neutral hydrogen surveys, where the 21 cm line is the only bright line, surveys of red-shifted fine-structure lines suffer from confusion generated by line broadening, spectral overlap of different lines, and the crowding of sources with redshift. We use simulations to investigate the resulting spectral confusion and derive observing parameters to minimize these effects in pencil-beam surveys of red-shifted far-IR line emission. We generate simulated spectra of the 17 brightest far-IR lines in galaxies, covering the 150 to 1300 micron wavelength region corresponding to redshifts 0 < z < 7, and develop a simple iterative algorithm that successfully identifies the 158 micron [CII] line and other lines. Although the [CII] line is a principal coolant for the interstellar medium, the assumption that the brightest observed lines in a given line of sight are always [CII] lines is a poor approximation to the simulated spectra once other lines are included. Blind line identification requires detection of fainter companion lines from the same host galaxies, driving survey sensitivity requirements. The observations require moderate spectral resolution 700 < R < 4000 with angular resolution between 20 arcsec and 10 armin, sufficiently narrow to minimize confusion yet sufficiently large to include a statistically meaningful number of sources.


title: Robust Neighbor Discovery in Multi-Hop Multi-Channel Heterogeneous Wireless Networks

abstract: An important first step when deploying a wireless ad hoc network is neighbor discovery in which every node attempts to determine the set of nodes it can communicate with in one wireless hop. In the recent years, cognitive radio (CR) technology has gained attention as an attractive approach to alleviate spectrum congestion. A cognitive radio transceiver can operate over a wide range of frequencies, possibly scanning multiple frequency bands. A cognitive radio node can opportunistically utilize unused wireless spectrum without interference from other wireless devices in its vicinity. Due to spatial variations in frequency usage and hardware variations in radio transceivers, different nodes in the network may perceive different subsets of frequencies available to them for communication. This heterogeneity in the available channel sets across the network increases the complexity of solving the neighbor discovery problem in a cognitive radio network. In this work, we design and analyze several randomized algorithms for neighbor discovery in such a (heterogeneous) network under a variety of assumptions (e.g. maximum node degree known or unknown) for both synchronous and asynchronous systems under minimal knowledge. We also show that our randomized algorithms are naturally suited to tolerate unreliable channels and adversarial attacks.


title: Circumnuclear Media of Quiescent Supermassive Black Holes

abstract: We calculate steady-state, one-dimensional hydrodynamic profiles of hot gas in slowly accreting ("quiescent") galactic nuclei for a range of central black hole masses $M_{\bullet}$, parametrized gas heating rates, and observationally-motivated stellar density profiles. Mass is supplied to the circumnuclear medium by stellar winds, while energy is injected primarily by stellar winds, supernovae, and black hole feedback. Analytic estimates are derived for the stagnation radius (where the radial velocity of the gas passes through zero) and the large scale gas inflow rate, $\dot{M}$, as a function of $M_{\bullet}$ and the gas heating efficiency, the latter being related to the star-formation history. We assess the conditions under which radiative instabilities develop in the hydrostatic region near the stagnation radius, both in the case of a single burst of star formation and for the average star formation history predicted by cosmological simulations. By combining a sample of measured nuclear X-ray luminosities, $L_x$, of nearby quiescent galactic nuclei with our results for $\dot{M}(M_{\bullet})$ we address whether the nuclei are consistent with accreting in a steady-state, thermally-stable manner for radiative efficiencies predicted for radiatively inefficiency accretion flows. We find thermally-stable accretion cannot explain the short average growth times of low mass black holes in the local Universe, which must instead result from gas being fed in from large radii, due either to gas inflows or thermal instabilities acting on larger, galactic scales. Our results have implications for attempts to constrain the occupation fraction of SMBHs in low mass galaxies using the mean $L_x-M_{\bullet}$ correlation, as well as the predicted diversity of the circumnuclear densities encountered by relativistic outflows from tidal disruption events.


title: Variability in the super-Earth 55 Cnc e

abstract: Considerable progress has been made in recent years in observations of atmospheric signatures of giant exoplanets, but processes in rocky exoplanets remain largely unknown due to major challenges in observing small planets. Numerous efforts to observe spectra of super-Earths, exoplanets with masses of 1-10 Earth masses, have thus far revealed only featureless spectra. In this paper we report a 4-$\sigma$ detection of variability in the dayside thermal emission from the transiting super-Earth 55 Cancri e. Dedicated space-based monitoring of the planet in the mid-infrared over eight eclipses revealed the thermal emission from its dayside atmosphere varying by a factor 3.7 between 2012 and 2013. The amplitude and trend of the variability are not explained by potential influence of star spots or by local thermal or compositional changes in the atmosphere over the short span of the observations. The possibility of large scale surface activity due to strong tidal interactions possibly similar to Io, or the presence of circumstellar/circumplanetary material appear plausible and motivate future long-term monitoring of the planet.


title: Ultra-stripped supernovae: progenitors and fate

abstract: The explosion of ultra-stripped stars in close binaries can lead to ejecta masses < 0.1 M_sun and may explain some of the recent discoveries of weak and fast optical transients. In Tauris et al. (2013), it was demonstrated that helium star companions to neutron stars (NSs) may experience mass transfer and evolve into naked ~1.5 M_sun metal cores, barely above the Chandrasekhar mass limit. Here we present a systematic investigation of the progenitor evolution leading to ultra-stripped supernovae (SNe). In particular, we examine the binary parameter space leading to electron-capture (EC SNe) and iron core-collapse SNe (Fe CCSNe), respectively, and determine the amount of helium ejected with applications to their observational classification as Type Ib or Type Ic. We mainly evolve systems where the SN progenitors are helium star donors of initial mass M_He = 2.5 - 3.5 M_sun in tight binaries with orbital periods of P_orb = 0.06 - 2.0 days, and hosting an accreting NS, but we also discuss the evolution of wider systems and of both more massive and lighter - as well as single - helium stars. In some cases we are able to follow the evolution until the onset of silicon burning, just a few days prior to the SN explosion. We find that ultra-stripped SNe are possible for both EC SNe and Fe CCSNe, and that the amount of helium ejected is correlated with P_orb - the tightest systems even having donors being stripped down to envelopes of less than 0.01 M_sun. We estimate the rise time of ultra-stripped SNe to be in the range 12 hr - 8 days, and light curve decay times between 1 and 50 days. Ultra-stripped SNe may produce NSs in the mass range 1.10 - 1.80 M_sun and are highly relevant for LIGO/VIRGO since most (possibly all) merging double NS systems have evolved through this phase. Finally, we discuss the low-velocity kicks which might be imparted on these resulting NSs at birth. [Abridged]


title: On the role of chalcogen vapor annealing in inducing bulk superconductivity in Fe$_{1+y}$Te$_{1-x}$Se$_{x}$

abstract: Recent investigations have shown that Fe$_{1+y}$Te$_{1-x}$Se$_{x}$ can be made superconducting by annealing it in Se and O vapors. The current lore is that these chalcogen vapors induce superconductivity by removing the magnetic excess Fe atoms. To investigate this phenomenon we performed a combination of magnetic susceptibility, specific heat and transport measurements together with scanning tunneling microscopy and spectroscopy and density functional theory calculations on Fe$_{1+y}$Te$_{1-x}$Se$_{x}$ treated with Te vapor. We conclude that the main role of the Te vapor is to quench the magnetic moments of the excess Fe atoms by forming FeTe$_{m}$ (m $\geq$ 1) complexes. We show that the remaining FeTe$_{m}$ complexes are still damaging to the superconductivity and therefore that their removal potentially could further improve superconductive properties in these compounds.


title: Pauli-Villars Regularization and Light Front Hamiltonian in (2+1)-dimensional Yang-Mills Theory

abstract: The renormalization problem of (2+1)-dimensional Yang-Mills theory quantized on the light front is considered. Extra fields analogous to those used in Pauli-Villars regularization are introduced to restore perturbative equivalence between such quantized theory and conventional formulation in Lorentz coordinates. These fields also provide necessary ultraviolet regularization to the theory. Obtained results allow to construct renormalized Hamiltonian of the theory on the light front.


title: A Survey of Stochastic Simulation and Optimization Methods in Signal Processing

abstract: Modern signal processing (SP) methods rely very heavily on probability and statistics to solve challenging SP problems. SP methods are now expected to deal with ever more complex models, requiring ever more sophisticated computational inference techniques. This has driven the development of statistical SP methods based on stochastic simulation and optimization. Stochastic simulation and optimization algorithms are computationally intensive tools for performing statistical inference in models that are analytically intractable and beyond the scope of deterministic inference methods. They have been recently successfully applied to many difficult problems involving complex statistical models and sophisticated (often Bayesian) statistical inference techniques. This survey paper offers an introduction to stochastic simulation and optimization methods in signal and image processing. The paper addresses a variety of high-dimensional Markov chain Monte Carlo (MCMC) methods as well as deterministic surrogate methods, such as variational Bayes, the Bethe approach, belief and expectation propagation and approximate message passing algorithms. It also discusses a range of optimization methods that have been adopted to solve stochastic problems, as well as stochastic methods for deterministic optimization. Subsequently, areas of overlap between simulation and optimization, in particular optimization-within-MCMC and MCMC-driven optimization are discussed.


title: Stick-Breaking Policy Learning in Dec-POMDPs

abstract: Expectation maximization (EM) has recently been shown to be an efficient algorithm for learning finite-state controllers (FSCs) in large decentralized POMDPs (Dec-POMDPs). However, current methods use fixed-size FSCs and often converge to maxima that are far from optimal. This paper considers a variable-size FSC to represent the local policy of each agent. These variable-size FSCs are constructed using a stick-breaking prior, leading to a new framework called \emph{decentralized stick-breaking policy representation} (Dec-SBPR). This approach learns the controller parameters with a variational Bayesian algorithm without having to assume that the Dec-POMDP model is available. The performance of Dec-SBPR is demonstrated on several benchmark problems, showing that the algorithm scales to large problems while outperforming other state-of-the-art methods.


title: Local Orthogonal Polynomial Expansion for Density Estimation

abstract: A Local Orthogonal Polynomial Expansion (LOrPE) of the empirical density function is proposed as a novel method to estimate the underlying density. The estimate is constructed by matching localized expectation values of orthogonal polynomials to the values observed in the sample. LOrPE is related to several existing methods, and generalizes straightforwardly to multivariate settings. By manner of construction, it is similar to Local Likelihood Density Estimation (LLDE). In the limit of small bandwidths, LOrPE functions as Kernel Density Estimation (KDE) with high-order (effective) kernels inherently free of boundary bias, a natural consequence of kernel reshaping to accommodate endpoints. Faster asymptotic convergence rates follow. In the limit of large bandwidths, LOrPE is equivalent to Orthogonal Series Density Estimation (OSDE) with Legendre polynomials. We compare the performance of LOrPE to KDE, LLDE, and OSDE, in a number of simulation studies. In terms of mean integrated squared error, the results suggest that with a proper balance of the two tuning parameters, bandwidth and degree, LOrPE generally outperforms these competitors when estimating densities with sharply truncated supports.


title: Joint Object and Part Segmentation using Deep Learned Potentials

abstract: Segmenting semantic objects from images and parsing them into their respective semantic parts are fundamental steps towards detailed object understanding in computer vision. In this paper, we propose a joint solution that tackles semantic object and part segmentation simultaneously, in which higher object-level context is provided to guide part segmentation, and more detailed part-level localization is utilized to refine object segmentation. Specifically, we first introduce the concept of semantic compositional parts (SCP) in which similar semantic parts are grouped and shared among different objects. A two-channel fully convolutional network (FCN) is then trained to provide the SCP and object potentials at each pixel. At the same time, a compact set of segments can also be obtained from the SCP predictions of the network. Given the potentials and the generated segments, in order to explore long-range context, we finally construct an efficient fully connected conditional random field (FCRF) to jointly predict the final object and part labels. Extensive evaluation on three different datasets shows that our approach can mutually enhance the performance of object and part segmentation, and outperforms the current state-of-the-art on both tasks.


title: Grounded Discovery of Coordinate Term Relationships between Software Entities

abstract: We present an approach for the detection of coordinate-term relationships between entities from the software domain, that refer to Java classes. Usually, relations are found by examining corpus statistics associated with text entities. In some technical domains, however, we have access to additional information about the real-world objects named by the entities, suggesting that coupling information about the "grounded" entities with corpus statistics might lead to improved methods for relation discovery. To this end, we develop a similarity measure for Java classes using distributional information about how they are used in software, which we combine with corpus statistics on the distribution of contexts in which the classes appear in text. Using our approach, cross-validation accuracy on this dataset can be improved dramatically, from around 60% to 88%. Human labeling results show that our classifier has an F1 score of 86% over the top 1000 predicted pairs.


title: Automatic Observer Script for StarCraft: Brood War Bot Games (technical report)

abstract: This short report describes an automated BWAPI-based script developed for live streams of a StarCraft Brood War bot tournament, SSCAIT. The script controls the in-game camera in order to follow the relevant events and improve the viewer experience. We enumerate its novel features and provide a few implementation notes.


title: Thermal evolution of the Schwinger model with Matrix Product Operators

abstract: We demonstrate the suitability of tensor network techniques for describing the thermal evolution of lattice gauge theories. As a benchmark case, we have studied the temperature dependence of the chiral condensate in the Schwinger model, using matrix product operators to approximate the thermal equilibrium states for finite system sizes with non-zero lattice spacings. We show how these techniques allow for reliable extrapolations in bond dimension, step width, system size and lattice spacing, and for a systematic estimation and control of all error sources involved in the calculation. The reached values of the lattice spacing are small enough to capture the most challenging region of high temperatures and the final results are consistent with the analytical prediction by Sachs and Wipf over a broad temperature range.


title: Thorium Abundances in Solar Twins and Analogues: Implications for the Habitability of Extrasolar Planetary Systems

abstract: We present the first investigation of Th abundances in Solar twins and analogues to understand the possible range of this radioactive element and its effect on rocky planet interior dynamics and potential habitability. The abundances of the radioactive elements Th and U are key components of a planet's energy budget, making up 30% to 50% of the Earth's (Korenaga 2008; All\`egre et al. 2001; Schubert et al. 1980; Lyubetskaya & Korenaga 2007; The KamLAND Collaboration 2011; Huang et al. 2013). Radiogenic heat drives interior mantle convection and surface plate tectonics, which sustains a deep carbon and water cycle and thereby aides in creating Earth's habitable surface. Unlike other heat sources that are dependent on the planet's specific formation history, the radiogenic heat budget is directly related to the mantle concentration of these nuclides. As a refractory element, the stellar abundance of Th is faithfully reflected in the terrestrial planet's concentration. We find that log eps Th varies from 59% to 251% that of Solar, suggesting extrasolar planetary systems may possess a greater energy budget with which to support surface to interior dynamics and thus increase their likelihood to be habitable compared to our Solar System.


title: The representation of integers by positive ternary quadratic polynomials

abstract: An integral quadratic polynomial is called regular if it represents every integer that is represented by the polynomial itself over the reals and over the $p$-adic integers for every prime $p$. It is called complete if it is of the form $Q({\mathbf x} + {\mathbf v})$, where $Q$ is an integral quadratic form in the variables ${\mathbf x} = (x_1, \ldots, x_n)$ and ${\mathbf v}$ is a vector in ${\mathbb Q}^n$. Its conductor is defined to be the smallest positive integer $c$ such that $c{\mathbf v} \in {\mathbb Z}^n$. We prove that for a fixed positive integer $c$, there are only finitely many equivalence classes of positive primitive ternary regular complete quadratic polynomials with conductor $c$. This generalizes the analogous finiteness results for positive definite regular ternary quadratic forms by Watson and for ternary triangular forms by Chan and Oh.


title: Which is conducive to high responsivity in a hybrid graphene-quantum dot transistor: Bulk- or layer- heterojunction

abstract: Heterojunction, a photoelectric conversion center, plays a critical role in photo detection. A compare study is performed on field effect phototransistors (FEpT) with two typical heterojunctions based on graphene and PbSe quantum dots (QDs) hybrid, including layer-heterojunction (LH) and bulk-heterojunction (BH). LH-FEpT exhibits a higher mobility ({\mu}n) 677 cm2V-1s-1, whereas {\mu}n of BH-FEpT is as low as 314 cm2V-1s-1, however the higher responsivity (106 A/W) is achieved by the former. Because R is proportional to the product of mobility and transfer rate, and high transfer rate of LH-FEpT overcompensates the shortage of low mobility. Although large area of heterojunction and high mobility benefit high density of photo-induced carriers, lack of transport mechanism becomes to the main constrain factor in BH-FEpT. Therefore, LH-FEpT is a better candidate of near infrared (NIR) photo detectors.


title: A new approach to codeword stabilized quantum codes using the algebraic structure of modules

abstract: In this work, we study the Codeword Stabilized Quantum Codes (CWS codes) a generalization of the stabilizers quantum codes using a new approach, the algebraic structure of modules, a generalization of linear spaces. We show then a new result that relates CWS codes with stabilizer codes generalizing results in the literature.


title: Bayesian Policy Reuse

abstract: A long-lived autonomous agent should be able to respond online to novel instances of tasks from a familiar domain. Acting online requires 'fast' responses, in terms of rapid convergence, especially when the task instance has a short duration, such as in applications involving interactions with humans. These requirements can be problematic for many established methods for learning to act. In domains where the agent knows that the task instance is drawn from a family of related tasks, albeit without access to the label of any given instance, it can choose to act through a process of policy reuse from a library, rather than policy learning from scratch. In policy reuse, the agent has prior knowledge of the class of tasks in the form of a library of policies that were learnt from sample task instances during an offline training phase. We formalise the problem of policy reuse, and present an algorithm for efficiently responding to a novel task instance by reusing a policy from the library of existing policies, where the choice is based on observed 'signals' which correlate to policy performance. We achieve this by posing the problem as a Bayesian choice problem with a corresponding notion of an optimal response, but the computation of that response is in many cases intractable. Therefore, to reduce the computation cost of the posterior, we follow a Bayesian optimisation approach and define a set of policy selection functions, which balance exploration in the policy library against exploitation of previously tried policies, together with a model of expected performance of the policy library on their corresponding task instances. We validate our method in several simulated domains of interactive, short-duration episodic tasks, showing rapid convergence in unknown task variations.


title: The temporal evolution of the energy flux across scales in homogeneous turbulence

abstract: A temporal study of energy transfer across length scales is performed in 3D numerical simulations of homogeneous shear flow and isotropic turbulence. The average time taken by perturbations in the energy flux to travel between scales is measured and shown to be additive. Our data suggests that the propagation of disturbances in the energy flux is independent of the forcing and that it defines a `velocity' that determines the energy flux itself. These results support that the cascade is, on average, a scale-local process where energy is continuously transmitted from one scale to the next in order of decreasing size.


title: An analytic regularisation scheme on curved spacetimes with applications to cosmological spacetimes

abstract: We develop a renormalisation scheme for time--ordered products in interacting field theories on curved spacetimes which consists of an analytic regularisation of Feynman amplitudes and a minimal subtraction of the resulting pole parts. This scheme is directly applicable to spacetimes with Lorentzian signature, manifestly generally covariant, invariant under any spacetime isometries present and constructed to all orders in perturbation theory. Moreover, the scheme captures correctly the non--geometric state--dependent contribution of Feynman amplitudes and it is well--suited for practical computations. To illustrate this last point, we compute explicit examples on a generic curved spacetime, and demonstrate how momentum space computations in cosmological spacetimes can be performed in our scheme. In this work, we discuss only scalar fields in four spacetime dimensions, but we argue that the renormalisation scheme can be directly generalised to other spacetime dimensions and field theories with higher spin, as well as to theories with local gauge invariance.


title: Matrix product formula for Macdonald polynomials

abstract: We derive a matrix product formula for symmetric Macdonald polynomials. Our results are obtained by constructing polynomial solutions of deformed Knizhnik--Zamolodchikov equations, which arise by considering representations of the Zamolodchikov--Faddeev and Yang--Baxter algebras in terms of $t$-deformed bosonic operators. These solutions form a basis of the ring of polynomials in $n$ variables, whose elements are indexed by compositions. For weakly increasing compositions (anti-dominant weights), these basis elements coincide with non-symmetric Macdonald polynomials. Our formulas imply a natural combinatorial interpretation in terms of solvable lattice models. They also imply that normalisations of stationary states of multi-species exclusion processes are obtained as Macdonald polynomials at $q=1$.


title: New copulas based on general partitions-of-unity and their applications to risk management

abstract: We construct new multivariate copulas on the basis of a generalized infinite partition-of-unity approach. This approach allows - in contrast to finite partition-of-unity copulas - for tail-dependence as well as for asymmetry. A possibility of fitting such copulas to real data from quantitative risk management is also pointed out.


title: Deep Remix: Remixing Musical Mixtures Using a Convolutional Deep Neural Network

abstract: Audio source separation is a difficult machine learning problem and performance is measured by comparing extracted signals with the component source signals. However, if separation is motivated by the ultimate goal of re-mixing then complete separation is not necessary and hence separation difficulty and separation quality are dependent on the nature of the re-mix. Here, we use a convolutional deep neural network (DNN), trained to estimate 'ideal' binary masks for separating voice from music, to perform re-mixing of the vocal balance by operating directly on the individual magnitude components of the musical mixture spectrogram. Our results demonstrate that small changes in vocal gain may be applied with very little distortion to the ultimate re-mix. Our method may be useful for re-mixing existing mixes.


title: Algorithms for Lipschitz Learning on Graphs

abstract: We develop fast algorithms for solving regression problems on graphs where one is given the value of a function at some vertices, and must find its smoothest possible extension to all vertices. The extension we compute is the absolutely minimal Lipschitz extension, and is the limit for large $p$ of $p$-Laplacian regularization. We present an algorithm that computes a minimal Lipschitz extension in expected linear time, and an algorithm that computes an absolutely minimal Lipschitz extension in expected time $\widetilde{O} (m n)$. The latter algorithm has variants that seem to run much faster in practice. These extensions are particularly amenable to regularization: we can perform $l_{0}$-regularization on the given values in polynomial time and $l_{1}$-regularization on the initial function values and on graph edge weights in time $\widetilde{O} (m^{3/2})$.


title: Dip or nothingness of a Higgs resonance from the interference with a complex phase

abstract: We show that new resonance shapes -- a pure dip, nothingness and an enhanced pure peak -- can be produced from the interference between resonance and continuum with a relative phase. Production conditions of those new shapes are derived based on a general parameterization of the interference. The narrow width approximation is modified to work with the non-zero imaginary part of interference, and the correction factor can characterize the resonance shape. We demonstrate that the new resonance shapes of heavy Higgs bosons, $H^0$ and $A^0$ in the Type II aligned two Higgs doublet model, generally show up in $gg \to H^0/A^0 \to t\bar{t}$ as well as $b\bar{b}$ and $\gamma \gamma$ channels. The pure $A^0$ resonance dip in the $t\bar{t}$ channel is a particularly interesting signal as it can be probed well by the current search techniques that do not even take into account interferences; the high-luminosity LHC 14 TeV can perhaps probe a large part of its parameter space.


title: Free-space quantum key distribution to a moving receiver

abstract: Technological realities limit terrestrial quantum key distribution (QKD) to single-link distances of a few hundred kilometers. One promising avenue for global-scale quantum communication networks is to use low-Earth-orbit satellites. Here we report the first demonstration of QKD from a stationary transmitter to a receiver platform traveling at an angular speed equivalent to a 600 km altitude satellite, located on a moving truck. We overcome the challenges of actively correcting beam pointing, photon polarization and time-of-flight. Our system generates an asymptotic secure key at 40 bits/s.


title: Oscillatory superfluid Ekman pumping in Helium II and neutron stars

abstract: The linear response of a superfluid, rotating uniformly in a cylindrical container and threaded with a large number of vortex lines, to an impulsive increase in the angular velocity of the container is investigated. At zero temperature and with perfect pinning of vortices to the top and bottom of the container, we demonstrate that the system oscillates persistently with a frequency proportional to the vortex line tension parameter to the quarter power. This low-frequency mode is generated by a secondary flow analogous to classical Ekman pumping that is periodically reversed by the vortex tension in the boundary layers. We compare analytic solutions to the two-fluid equations of Chandler & Baym (1986) with the spin-up experiments of Tsakadze & Tsakadze (1980) in helium II and find the frequency agrees within a factor of four, although the experiment is not perfectly suited to the application of the linear theory. We argue that this oscillatory Ekman pumping mode, and not Tkachenko modes provide a natural explanation for the observed oscillation. In neutron stars, the oscillation period depends on the pinning interaction between neutron vortices and flux tubes in the outer core. Using a simplified pinning model, we demonstrate that strong pinning can accommodate modes with periods of days to years, which are only weakly damped by mutual friction over longer timescales.


title: Monotonous (Semi-)Nonnegative Matrix Factorization

abstract: Nonnegative matrix factorization (NMF) factorizes a non-negative matrix into product of two non-negative matrices, namely a signal matrix and a mixing matrix. NMF suffers from the scale and ordering ambiguities. Often, the source signals can be monotonous in nature. For example, in source separation problem, the source signals can be monotonously increasing or decreasing while the mixing matrix can have nonnegative entries. NMF methods may not be effective for such cases as it suffers from the ordering ambiguity. This paper proposes an approach to incorporate notion of monotonicity in NMF, labeled as monotonous NMF. An algorithm based on alternating least-squares is proposed for recovering monotonous signals from a data matrix. Further, the assumption on mixing matrix is relaxed to extend monotonous NMF for data matrix with real numbers as entries. The approach is illustrated using synthetic noisy data. The results obtained by monotonous NMF are compared with standard NMF algorithms in the literature, and it is shown that monotonous NMF estimates source signals well in comparison to standard NMF algorithms when the underlying sources signals are monotonous.


title: Dense Optical Flow Prediction from a Static Image

abstract: Given a scene, what is going to move, and in what direction will it move? Such a question could be considered a non-semantic form of action prediction. In this work, we present a convolutional neural network (CNN) based approach for motion prediction. Given a static image, this CNN predicts the future motion of each and every pixel in the image in terms of optical flow. Our CNN model leverages the data in tens of thousands of realistic videos to train our model. Our method relies on absolutely no human labeling and is able to predict motion based on the context of the scene. Because our CNN model makes no assumptions about the underlying scene, it can predict future optical flow on a diverse set of scenarios. We outperform all previous approaches by large margins.


title: Object-Scene Convolutional Neural Networks for Event Recognition in Images

abstract: Event recognition from still images is of great importance for image understanding. However, compared with event recognition in videos, there are much fewer research works on event recognition in images. This paper addresses the issue of event recognition from images and proposes an effective method with deep neural networks. Specifically, we design a new architecture, called Object-Scene Convolutional Neural Network (OS-CNN). This architecture is decomposed into object net and scene net, which extract useful information for event understanding from the perspective of objects and scene context, respectively. Meanwhile, we investigate different network architectures for OS-CNN design, and adapt the deep (AlexNet) and very-deep (GoogLeNet) networks to the task of event recognition. Furthermore, we find that the deep and very-deep networks are complementary to each other. Finally, based on the proposed OS-CNN and comparative study of different network architectures, we come up with a solution of five-stream CNN for the track of cultural event recognition at the ChaLearn Looking at People (LAP) challenge 2015. Our method obtains the performance of 85.5% and ranks the $1^{st}$ place in this challenge.


title: Two-Dimensional Pursuit-Evasion in a Compact Domain with Piecewise Analytic Boundary

abstract: In a pursuit-evasion game, a team of pursuers attempt to capture an evader. The players alternate turns, move with equal speed, and have full information about the state of the game. We consider the most restictive capture condition: a pursuer must become colocated with the evader to win the game. We prove two general results about pursuit-evasion games in topological spaces. First, we show that one pursuer has a winning strategy in any CAT(0) space under this restrictive capture criterion. This complements a result of Alexander, Bishop and Ghrist, who provide a winning strategy for a game with positive capture radius. Second, we consider the game played in a compact domain in Euclidean two-space with piecewise analytic boundary and arbitrary Euler characteristic. We show that three pursuers always have a winning strategy by extending recent work of Bhadauria, Klein, Isler and Suri from polygonal environments to our more general setting.


title: Modular invariance and anomaly cancellation formulas in odd dimension

abstract: By studying modular invariance properties of some characteristic forms, we get some new anomaly cancellation formulas on $(4r-1)$ dimensional manifolds. As an application, we derive some results on divisibilities of the index of Toeplitz operators on $(4r-1)$ dimensional spin manifolds and some congruent formulas on characteristic number for $(4r-1)$ dimensional spin$^c$ manifolds.


title: Compressed Sensing Based Multi-User Millimeter Wave Systems: How Many Measurements Are Needed?

abstract: Millimeter wave (mmWave) systems will likely employ directional beamforming with large antenna arrays at both the transmitters and receivers. Acquiring channel knowledge to design these beamformers, however, is challenging due to the large antenna arrays and small signal-to-noise ratio before beamforming. In this paper, we propose and evaluate a downlink system operation for multi-user mmWave systems based on compressed sensing channel estimation and conjugate analog beamforming. Adopting the achievable sum-rate as a performance metric, we show how many compressed sensing measurements are needed to approach the perfect channel knowledge performance. The results illustrate that the proposed algorithm requires an order of magnitude less training overhead compared with traditional lower-frequency solutions, while employing mmWave-suitable hardware. They also show that the number of measurements need to be optimized to handle the trade-off between the channel estimate quality and the training overhead.


title: On the route to shear jamming, are fragile states real?

abstract: Starting from an unjammed initial state, applying shear to a granular material of a fixed packing fraction below $\phi_J$, i.e. the isotropic jamming density of frictionless spheres can produce shear jamming states, as have been discovered recently. In addition, it has also been discovered that the system will first experience a bulk fragile state before evolving into a shear jammed state. Due to the existence of friction between the system and the third dimension in the previous studies, it is unclear whether such fragile state would still exist on the route to shear jamming if the friction with the third dimension were completely eliminated and the background noise level were greatly reduced. Using a novel apparatus, we have completely eliminated the friction between particles and the third dimension by floating the particles on the surface of a shallow water layer thus revealing more details of the route of shear-jamming. We are able to measure weak boundary pressure of three orders of magnitude below the resolution of the photo-elastic method through the combination force-gauges of high sensitivity and a cantilever-like simple-beam apparatus. In a system with weak cohesion, we have indeed observed the bulk fragile states before the shear jamming; in addition, we have also discovered the boundary fragile states and regimes of negligible shear modus compared to bulk modulus. Most importantly, we now have a much better understanding of the bulk fragile state: its existence requires an auxiliary in some form, e.g. the weak cohesion in this experiment. By constructing a way to freely tune down the cohesion to more than six orders of magnitude smaller, to our surprise, we have discovered the complete vanishing of the bulk fragile states, which thus allows us to determine the plausible origin of the bulk fragile state.


title: Trace-distance correlations for X states and emergence of the pointer basis in Markovian and non-Markovian regimes

abstract: We provide analytical expressions for classical and total trace-norm (Schatten 1-norm) geometric correlations in the case of two-qubit X states. As an application, we consider the open-system dynamical behavior of such correlations under phase and generalized amplitude damping evolutions. Then, we show that geometric classical correlations can characterize the emergence of the pointer basis of an apparatus subject to decoherence in either Markovian or non-Markovian regimes. In particular, as a non-Markovian effect, we obtain a time delay for the information to be retrieved from the apparatus by a classical observer. Moreover, we show that the set of initial X states exhibiting sudden transitions in the geometric classical correlation has nonzero measure.


title: Maximal amenability and disjointness for the radial masa

abstract: We prove that the radial masa C in the free group factor is disjoint from other maximal amenable subalgebras in the following sense: any distinct maximal amenable subalgebra cannot have diffuse intersection with C.


title: Achievable Rates of Multi-User Millimeter Wave Systems with Hybrid Precoding

abstract: Millimeter wave (mmWave) systems will likely employ large antenna arrays at both the transmitters and receivers. A natural application of antenna arrays is simultaneous transmission to multiple users, which requires multi-user precoding at the transmitter. Hardware constraints, however, make it difficult to apply conventional lower frequency MIMO precoding techniques at mmWave. This paper proposes and analyzes a low complexity hybrid analog/digital beamforming algorithm for downlink multi-user mmWave systems. Hybrid precoding involves a combination of analog and digital processing that is motivated by the requirement to reduce the power consumption of the complete radio frequency and mixed signal hardware. The proposed algorithm configures hybrid precoders at the transmitter and analog combiners at multiple receivers with a small training and feedback overhead. For this algorithm, we derive a lower bound on the achievable rate for the case of single-path channels, show its asymptotic optimality at large numbers of antennas, and make useful insights for more general cases. Simulation results show that the proposed algorithm offers higher sum rates compared with analog-only beamforming, and approaches the performance of the unconstrained digital precoding solutions.


title: Proof of a conjecture of Granath on optimal bounds of the Landau constants

abstract: We study the asymptotic expansion for the Landau constants $G_n$, \begin{equation*} \pi G_{n}\sim \ln(16N)+\gamma+\sum^{\infty}_{k=1}\frac{\alpha_k}{N^k} ~~\mbox{as} ~ n\rightarrow\infty, \end{equation*} where $N=n+1$, and $\gamma$ is Euler's constant. We show that the signs of the coefficients $\alpha_{k}$ demonstrate a periodic behavior such that $(-1)^{\frac {l(l+1)} 2} \alpha_{l+1}< 0$ for all $l$. We further prove a conjecture of Granath which states that $(-1)^{\frac {l(l+1)} 2} \varepsilon_l(N)<0$ for $l=0,1,2,\cdots$ and $n=0,1,2,\cdots$, $\varepsilon_l(N)$ being the error due to truncation at the $l$-th order term. Consequently, we also obtain the sharp bounds up to arbitrary orders of the form \begin{equation*} \ln(16N)+\gamma+\sum_{k=1}^{p}\frac{\alpha_{k}}{N^{k}}<\pi G_{n}<\ln(16N)+\gamma+\sum_{k=1}^{q}\frac{\alpha_{k}}{N^{k}} \end{equation*} for all $n=0,1,2\cdots$, all $p=4s+1,\; 4s+2$ and $q=4m,\; 4m+3$, with $s=0,1,2,\cdots$ and $m=0, 1, 2,\cdots$.


title: Spectral Intensities of Antiprotons and the lifetime of Cosmic Rays in the Galaxy

abstract: In this paper we note that the spectral intensities of antiprotons observed in Galactic cosmic rays in the energy range ~ 1-100 GeV by BESS, PAMELA and AMS instruments display nearly the same spectral shape as that generated by primary cosmic rays through their interaction with matter in the interstellar medium, without any significant modifications. More importantly, a constant residence time of ~ 2.5 +/-0.7 million years in the Galactic volume, independent of the energy of cosmic rays, matches the observed intensities. A small additional component of secondary antiprotons in the energy below 10 GeV, generated in cocoon-like regions surrounding the cosmic-ray sources, seems to be present. We discuss this result in the context of observations of other secondary components like positrons and Boron, and conclude with general remarks about the origins and propagation of cosmic rays.


title: Diluting the inflationary axion fluctuation by a stronger QCD in the early Universe

abstract: We propose a new mechanism to suppress the axion isocurvature perturbation, while producing the right amount of axion dark matter, within the framework of supersymmetric axion models with the axion scale induced by supersymmetry breaking. The mechanism involves an intermediate phase transition to generate the Higgs \mu-parameter, before which the weak scale is comparable to the axion scale and the resulting stronger QCD yields an axion mass heavier than the Hubble scale over a certain period. Combined with that the Hubble-induced axion scale during the primordial inflation is well above the intermediate axion scale at present, the stronger QCD in the early Universe suppresses the axion fluctuation to be small enough even when the inflationary Hubble scale saturates the current upper bound, while generating an axion misalignment angle of order unity.


title: Optical spectroscopy study of three dimensional Dirac semimetal ZrTe$_5$

abstract: Three dimensional (3D) topological Dirac materials are under intensive study recently. The layered compound ZrTe$_5$ has been suggested to be one of them by transport and ARPES experiments. Here, we perform infrared reflectivity measurement to investigate the underlying physics of this material. The derived optical conductivity exhibits linear increasing with frequency below normal interband transitions, which provides the first optical spectroscopic proof of a 3D Dirac semimetal. Apart from that, the plasma edge shifts dramatically to lower energy upon temperature cooling, which might be associated with the consequence of lattice parameter shrinking. In addition, an extremely sharp peak shows up in the frequency dependent optical conductivity, indicating the presence of a Van Hove singularity in the joint density of state.


title: Multi-Object Classification and Unsupervised Scene Understanding Using Deep Learning Features and Latent Tree Probabilistic Models

abstract: Deep learning has shown state-of-art classification performance on datasets such as ImageNet, which contain a single object in each image. However, multi-object classification is far more challenging. We present a unified framework which leverages the strengths of multiple machine learning methods, viz deep learning, probabilistic models and kernel methods to obtain state-of-art performance on Microsoft COCO, consisting of non-iconic images. We incorporate contextual information in natural images through a conditional latent tree probabilistic model (CLTM), where the object co-occurrences are conditioned on the extracted fc7 features from pre-trained Imagenet CNN as input. We learn the CLTM tree structure using conditional pairwise probabilities for object co-occurrences, estimated through kernel methods, and we learn its node and edge potentials by training a new 3-layer neural network, which takes fc7 features as input. Object classification is carried out via inference on the learnt conditional tree model, and we obtain significant gain in precision-recall and F-measures on MS-COCO, especially for difficult object categories. Moreover, the latent variables in the CLTM capture scene information: the images with top activations for a latent node have common themes such as being a grasslands or a food scene, and on on. In addition, we show that a simple k-means clustering of the inferred latent nodes alone significantly improves scene classification performance on the MIT-Indoor dataset, without the need for any retraining, and without using scene labels during training. Thus, we present a unified framework for multi-object classification and unsupervised scene understanding.


title: Numerical Study of Crystal Size Distribution in Polynuclear Growth

abstract: The crystal size distribution in polynuclear growth is numerically studied using a coupled map lattice model. The width of the size distribution depends on c/D, where c is the growth rate at interface sites and $D$ is the diffusion constant. When c/D is sufficiently small, the width W increases linearly with c/D and saturates at large c/D. Monodisperse square and cubic crystals are obtained respectively on square and cubic lattices when c/D is sufficiently small for a small kinetic parameter b. The linear dependence of W on c/D in a parameter range of small c/D is explained by the eigenfunction for the first eigenvalue in a two-dimensional model and a mean-field model. For the mean-field model, the slope of the linear dependence is evaluated theoretically.


title: Progress in vacuum susceptibilities and their applications to the chiral phase transition of QCD

abstract: The QCD vacuum condensates and various vacuum susceptibilities are all important parameters which characterize the nonperturbative properties of the QCD vacuum. In the QCD sum rules external field formula, various QCD vacuum susceptibilities play important roles in determining the properties of hadrons. In this paper, we review the recent progress in studies of vacuum susceptibilities together with their applications to the chiral phase transition of QCD. The results of the tensor, the vector, the axial-vector, the scalar, and the pseudo-scalar vacuum susceptibilities are shown in detail in the framework of Dyson-Schwinger equations.


title: Radio Afterglows and Host Galaxies of Gamma-Ray Bursts

abstract: Considering the contribution of the emission from the host galaxies of gamma-ray bursts (GRBs) to the radio afterglows, we investigate the effect of host galaxies on observations statistically. For the three types of events, e.g. low-luminosity, standard and high-luminosity GRBs, it is found that a tight correlation exists between the ratio of the radio flux (RRF) of host galaxy to the total radio peak emission and the observational frequency. Especially, toward lower frequencies, the contribution from the host increases significantly. The correlation can be used to get a useful estimate for the radio brightness of those host galaxies which only have very limited radio afterglow data. Using this prediction, we re-considered the theoretical radio afterglow light curves for four kinds of events, i.e. high-luminosity, low-luminosity, standard and failed GRBs, taking into account the contribution from the host galaxies and aiming at exploring the detectability of these events by the Five-hundred-meter Aperture Spherical radio Telescope (FAST). Lying at a typical redshift of $z=1$, most of the events can be detected by FAST easily. For the less fierce low-luminosity GRBs, their radio afterglows are not strong enough to exceed the sensitivity limit of FAST at such distances. However, since a large number of low luminosity bursts actually happen very near to us, it is expected that FAST will still be able to detect many of them.


title: The Descent of Math

abstract: A perplexing problem in understanding physical reality is why the universe seems comprehensible, and correspondingly why there should exist physical systems capable of comprehending it. In this essay I explore the possibility that rather than being an odd coincidence arising due to our strange position as passive (and even more strangely, conscious) observers in the cosmos, these two problems might be related and could be explainable in terms of fundamental physics. The perspective presented suggests a potential unified framework where, when taken together, comprehenders and comprehensibility are part of causal structure of physical reality, which is considered as a causal graph (network) connecting states that are physically possible. I argue that in some local regions, the most probable states are those that include physical systems which contain information encodings - such as mathematics, language and art - because these are the most highly connected to other possible states in this causal graph. Such physical systems include life and - of particular interest for the discussion of the place of math in physical reality - comprehenders capable of making mathematical sense of the world. Within this framework, the descent of math is an undirected outcome of the evolution of the universe, which will tend toward states that are increasingly connected to other possible states of the universe, a process greatly facilitated if some physical systems know the rules of the game. I therefore conclude that our ability to use mathematics to describe, and more importantly manipulate, the natural world may not be an anomaly or trick, but instead could provide clues to the underlying causal structure of physical reality.


title: Compton scattering of twisted light: angular distribution and polarization of scattered photons

abstract: Compton scattering of twisted photons is investigated within a non-relativistic framework using first-order perturbation theory. We formulate the problem in the density matrix theory, which enables one to gain new insights into scattering processes of twisted particles by exploiting the symmetries of the system. In particular, we analyze how the angular distribution and polarization of the scattered photons are affected by the parameters of the initial beam such as the opening angle and the projection of orbital angular momentum. We present analytical and numerical results for the angular distribution and the polarization of Compton scattered photons for initially twisted light and compare them with the standard case of plane-wave light.

categories: physics.atom-ph hep-ph physics.optics quant-ph


title: Deconstructing Principal Component Analysis Using a Data Reconciliation Perspective

abstract: Data reconciliation (DR) and Principal Component Analysis (PCA) are two popular data analysis techniques in process industries. Data reconciliation is used to obtain accurate and consistent estimates of variables and parameters from erroneous measurements. PCA is primarily used as a method for reducing the dimensionality of high dimensional data and as a preprocessing technique for denoising measurements. These techniques have been developed and deployed independently of each other. The primary purpose of this article is to elucidate the close relationship between these two seemingly disparate techniques. This leads to a unified framework for applying PCA and DR. Further, we show how the two techniques can be deployed together in a collaborative and consistent manner to process data. The framework has been extended to deal with partially measured systems and to incorporate partial knowledge available about the process model.


title: Learning Temporal Embeddings for Complex Video Analysis

abstract: In this paper, we propose to learn temporal embeddings of video frames for complex video analysis. Large quantities of unlabeled video data can be easily obtained from the Internet. These videos possess the implicit weak label that they are sequences of temporally and semantically coherent images. We leverage this information to learn temporal embeddings for video frames by associating frames with the temporal context that they appear in. To do this, we propose a scheme for incorporating temporal context based on past and future frames in videos, and compare this to other contextual representations. In addition, we show how data augmentation using multi-resolution samples and hard negatives helps to significantly improve the quality of the learned embeddings. We evaluate various design decisions for learning temporal embeddings, and show that our embeddings can improve performance for multiple video tasks such as retrieval, classification, and temporal order recovery in unconstrained Internet video.


title: The chiral phase transition with a chiral chemical potential in the framework of Dyson-Schwinger equations

abstract: Within the framework of Dyson-Schwinger equations (DSEs), we discuss the chiral phase transition of QCD with a chiral chemical potential $\mu_5$ as an additional scale. We focus especially on the issues related to the widely accepted as well as interested critical end point (CEP). With the help of a scalar susceptibility, we find that there might be no CEP$_5$ in the $T-\mu_5$ plane, and the phase transition in the $T-\mu_5$ plane might be totally crossover when $\mu<50$ MeV, which has apparent consistency with the Lattice QCD calculation. Our study may also provide some useful hints to some other studies related to $\mu_5$.


title: A problem of Petersson about weight 0 meromorphic modular forms

abstract: In this paper, we provide an explicit construction of weight $0$ meromorphic modular forms. Following work of Petersson, we build these via Poincar\'e series. There are two main aspects of our investigation which differ from his approach. Firstly, the naive definition of the Poincar\'e series diverges and one must analytically continue via Hecke's trick. Hecke's trick is further complicated in our situation by the fact that the Fourier expansion does not converge everywhere due to singularities in the upper half-plane so it cannot solely be used to analytically continue the functions. To explain the second difference, we recall that Petersson constructed linear combinations from a family of meromorphic functions which are modular if a certain principal parts condition is satisfied. In contrast to this, we construct linear combinations from a family of non-meromorphic modular forms, known as polar harmonic Maass forms, which are meromorphic whenever the principal parts condition is satisfied.


title: Centralizer algebras of the primitive unitary reflection group of order $96$

abstract: Among the unitary reflection groups, the one on the title is singled out by its importance in, for example, coding theory and number theory. In this paper we start with describing the irreducible representations of this group and then examine the semi-simple structure of the centralizer algebra in the tensor representation.


title: On the Bloch Theorem and Orthogonality Relations

abstract: Bloch theorem for a periodic operator is being revisited here, and we notice extra orthogonality relationships. It is shown that solutions are bi-periodic, in the sense that eigenfunctions are periodic with respect to one argument, and pseudo-periodic with respect to the other. An additional kind of symmetry between r-space and k-space exists between the envelope and eignfunctions not apparently noticed before, which allows to define new invertible modified Wannier functions. As opposed to the Wannier functions in r-space, these modified Wannier functions are defined in the k-spaces, but satisfy similar basic properties. These could result in new algorithms and other novel applications in the computational tools of periodic structures.


title: Comment on 'Semiclassical Klein-Kramers and Smoluchowski equations for the Brownian motion of a particle in an external potential'

abstract: The range of validity of the semiclassical Smoluchowski equation derived recently by Coffey et al is discussed. The analysis is based on the quantum Smoluchowski equation derived by the present author before. A quantum generalization of the Einstein law of Brownian motion is also obtained.


title: On the Riemann-Hurwitz formula for graph coverings

abstract: The aim of this paper is to present a few versions of the Riemann-Hurwitz formula for a regular branched covering of graphs. By a graph, we mean a finite connected multigraph. The genus of a graph is defined as the rank of the first homology group. We consider a finite group acting on a graph, possibly with fixed and invertible edges, and the respective factor graph. Then, the obtained Riemann-Hurwitz formula relates genus of the graph with genus of the factor graph and orders of the vertex and edge stabilisers.


title: Using PCA to Efficiently Represent State Spaces

abstract: Reinforcement learning algorithms need to deal with the exponential growth of states and actions when exploring optimal control in high-dimensional spaces. This is known as the curse of dimensionality. By projecting the agent's state onto a low-dimensional manifold, we can represent the state space in a smaller and more efficient representation. By using this representation during learning, the agent can converge to a good policy much faster. We test this approach in the Mario Benchmarking Domain. When using dimensionality reduction in Mario, learning converges much faster to a good policy. But, there is a critical convergence-performance trade-off. By projecting onto a low-dimensional manifold, we are ignoring important data. In this paper, we explore this trade-off of convergence and performance. We find that learning in as few as 4 dimensions (instead of 9), we can improve performance past learning in the full dimensional space at a faster convergence rate.


title: Thermal-induced proteinquake in oxyhemoglobin

abstract: Oxygen is released to living tissues via conformational changes of hemoglobin from R-state (oxyhemoglobin) to T-state (desoxyhemoglobin). The detailed mechanism of this process is not yet fully understood. We have carried out micromechanical experiments on oxyhemoglobin crystals to determine the behavior of the Young's modulus and the internal friction for temperatures between 20 C and 70 C. We have found that around 49 C oxyhemoglobin crystal samples undergo a sudden and strong increase of their Young's modulus, accompanied by a sudden decrease of the internal friction. This sudden mechanical change (proteinquake) takes place in a partially unfolded state and precedes the full denaturation transition at higher temperatures. The hemoglobin crystals after the proteinquake has the same mechanical properies as the initial state at room temperatures. We conjecture that it can be relevant for explaining the oxygen-releasing function of native oxyhemoglobin when the temperature is increased, e.g. due to active sport. The effect is specific for the quaternary structure of hemoglobin, and is absent for myoglobin with only one peptide sequence.


title: Strong optomechanical interactions in a sliced photonic crystal nanobeam

abstract: Cavity optomechanical systems can be used for sensitive detection of mechanical motion and to control mechanical resonators, down to the quantum level. The strength with which optical and mechanical degrees of freedom interact is defined by the photon-phonon coupling rate $g_0$, which is especially large in nanoscale systems. Here, we demonstrate an optomechanical system based on a sliced photonic crystal nanobeam, that combines subwavelength optical confinement with a low-mass mechanical mode. Analyzing the transduction of motion and effects of radiation pressure we find a coupling rate $g_0$/2{\pi} = 11.5 MHz, exceeding previously reported values by an order of magnitude. Using this interaction we detect the resonator's motion with a noise imprecision below that at the standard quantum limit, even though the system has optical and mechanical quality factors smaller than $10^3$. The broad bandwidth is useful for application in miniature sensors, and for measurement-based control of the resonator's motional state.


title: Two-dimensional segmentation of small convective patterns in radiation hydrodynamics simulations

abstract: Recent results from high-resolution solar granulation observations indicate the existence of a population of small granular cells that are smaller than 600 km in diameter. These small convective cells strongly contribute to the total area of granules and are located in the intergranular lanes, where they form clusters and chains. We study high-resolution radiation hydrodynamics simulations of the upper convection zone and photosphere to detect small granular cells, define their spatial alignment, and analyze their physical properties. We developed an automated image-segmentation algorithm specifically adapted to high-resolution simulations to identify granules. The resulting segmentation masks were applied to physical quantities, such as intensity and vertical velocity profiles, provided by the simulation. A new clustering algorithm was developed to study the alignment of small granular cells. This study shows that small granules make a distinct contribution to the total area of granules and form clusters of chain-like alignments. The simulation profiles demonstrate a different nature for small granular cells because they exhibit on average lower intensities, lower horizontal velocities, and are located deeper inside of convective layers than regular granules. Their intensity distribution deviates from a normal distribution as known for larger granules, and follows a Weibull distribution.


title: Combining Existential Rules and Description Logics (Extended Version)

abstract: Query answering under existential rules -- implications with existential quantifiers in the head -- is known to be decidable when imposing restrictions on the rule bodies such as frontier-guardedness [BLM10, BLMS11]. Query answering is also decidable for description logics [Baa03], which further allow disjunction and functionality constraints (assert that certain relations are functions), however, they are focused on ER-type schemas, where relations have arity two. This work investigates how to get the best of both worlds: having decidable existential rules on arbitrary arity relations, while allowing rich description logics, including functionality constraints, on arity-two relations. We first show negative results on combining such decidable languages. Second, we introduce an expressive set of existential rules (frontier-one rules with a certain restriction) which can be combined with powerful constraints on arity-two relations (e.g. GC 2, ALCQIb) while retaining decidable query answering. Further, we provide conditions to add functionality constraints on the higher-arity relations.


title: pH/$T$ duality - wall properties and time evolution of plant cells

abstract: We examined the pH/$T$ (or $\mu$/$T$) duality of acidic pH and temperature ($T$) for the growth of grass shoots in order to determine the equation of state (EoS) for living plants. By considering non-meristematic growth as a dynamic series of 'state transitions' (STs) in the extending primary wall, we identified the critical (read: optimum) exponents for this phenomenon, which exhibit a singular behaviour at a critical temperature, critical pH and critical chemical potential ($\mu$) in the form of four power laws: $F(\tau)_\pi\propto|\tau|^{\beta-1}$, $F(\pi)_\tau\propto|\pi|^{1-\alpha}$, $G(\tau)_\mu\propto|\tau|^{-2-\alpha+2\beta}$ and $G(\mu)_\tau\propto|\mu|^{2-\alpha}$. The power-law exponents $\alpha$ and $\beta$ are numbers, which are independent of pH (or $\mu$) and $T$ that are known as critical exponents, while $\pi$ and $\tau$ represent a reduced pH and reduced temperature, respectively. Various 'scaling' predictions were obtained - a convexity relation $\alpha + \beta \ge 2$ for practical pH-based analysis and a $\beta \equiv 2$ identity in microscopic representation. In the presented scenario, the magnitude that is decisive is the chemical potential of H$^+$ ions (protons), enforcing subsequent STs and growth. The EoS span areas of the biological, physical, chemical and Earth sciences cross the borders with the language (adapted formalism) of phase transitions.


title: Profitability of contrarian strategies in the Chinese stock market

abstract: This paper reexamines the profitability of loser, winner and contrarian portfolios in the Chinese stock market using monthly data of all stocks traded on the Shanghai Stock Exchange and Shenzhen Stock Exchange covering the period from January 1997 to December 2012. We find evidence of short-term and long-term contrarian profitability in the whole sample period when the estimation and holding horizons are 1 month or longer than 12 months and the annualized returns of contrarian portfolios increases with the estimation and holding horizons. We perform subperiod analysis and find that the long-term contrarian effect is significant in both bullish and bearish states while the short-term contrarian effect disappears in bullish states. We compare the performance of contrarian portfolios based on different grouping manners in the estimation period and unveil that decile grouping outperforms quintile grouping and tertile grouping, which is more evident and robust in the long run. Generally, loser portfolios and winner portfolios have positive returns and loser portfolios perform much better than winner portfolios. Both loser and winner portfolios in bullish states perform better than those in the whole sample period. In contrast, loser and winner portfolios have smaller returns in bearish states in which loser portfolio returns are significant only in the long term and winner portfolio returns become insignificant. These results are robust to the one-month skipping between the estimation and holding periods and for the two stock exchanges. Our findings show that the Chinese stock market is not efficient in the weak form. These findings also have obvious practical implications for financial practitioners.


title: Detection of small single-cycle signals by stochastic resonance using a bistable superconducting quantum interference device

abstract: We propose and experimentally demonstrate detecting small single-cycle and few-cycle signals by using the symmetric double-well potential of a radio frequency superconducting quantum interference device (rf-SQUID). We show that the response of this bistable system to single- and few-cycle signals has a non-monotonic dependence on the noise strength. The response, measured by the probability of transition from initial potential well to the opposite one, becomes maximum when the noise-induced transition rate between the two stable states of the rf-SQUID is comparable to the signal frequency. Comparison to numerical simulations shows that the phenomenon is a manifestation of stochastic resonance.


title: Linear Precoding of Data and Artificial Noise in Secure Massive MIMO Systems

abstract: In this paper, we consider secure downlink transmission in a multi-cell massive multiple-input multiple-output (MIMO) system where the numbers of base station (BS) antennas, mobile terminals, and eavesdropper antennas are asymptotically large. The channel state information of the eavesdropper is assumed to be unavailable at the BS and hence, linear precoding of data and artificial noise (AN) are employed for secrecy enhancement. Four different data precoders (i.e., selfish zero-forcing (ZF)/regularized channel inversion (RCI) and collaborative ZF/RCI precoders) and three different AN precoders (i.e., random, selfish/collaborative null-space based precoders) are investigated and the corresponding achievable ergodic secrecy rates are analyzed. Our analysis includes the effects of uplink channel estimation, pilot contamination, multi-cell interference, and path-loss. Furthermore, to strike a balance between complexity and performance, linear precoders that are based on matrix polynomials are proposed for both data and AN precoding. The polynomial coefficients of the data and AN precoders are optimized respectively for minimization of the sum mean squared error of and the AN leakage to the mobile terminals in the cell of interest using tools from free probability and random matrix theory. Our analytical and simulation results provide interesting insights for the design of secure multi-cell massive MIMO systems and reveal that the proposed polynomial data and AN precoders closely approach the performance of selfish RCI data and null-space based AN precoders, respectively.


title: Evolution of structural ({\alpha}) relaxation-time anomalies in GexSe1-x Chalcogenide glasses

abstract: We examine enthalpy relaxation across the chalcogenide glass series GexSe1-x, prepared over close-by compositions using conventional melt-quenching technique. We estimate the timescale {\tau}(Tg) characterizing enthalpic relaxation near the kinetic glass transition temperature, using the non-reversing heat-flow data obtained from MDSC (modulated differential scanning calorimetry) measurements over a wide range of compositions (2.1 {\lessthanequalto} r {\lessthanequalto} 2.8, r = 2x+2). Anomalies in the enthalpy-relaxation characteristic-times {\tau}g(r) are identified as marking rigidity-transitions encountered in successive Ge-doping of polymeric selenium chains.


title: Black hole radiation in the presence of a universal horizon

abstract: In Ho\v{r}ava and Einstein-{\AE}ther theories of modified gravity, in spite of the violation of Lorentz invariance, spherically-symmetric stationary black hole solutions possess an inner universal horizon which separates field configurations into two disconnected classes. We compute the late time radiation emitted by a dispersive field propagating in such backgrounds. We fix the initial conditions on stationary modes by considering a regular collapsing geometry, and imposing that the state inside the infalling shell is vacuum. We find that the mode pasting across the shell is adiabatic at late time (large inside frequencies). This implies that large black holes emit a thermal flux with a temperature fixed by the surface gravity of the Killing horizon. In turn, this suggests that the universal horizon should play no role in the thermodynamical properties of these black holes.


title: 1D3V PIC simulation of propagation of relativistic electron beam in an inhomogeneous plasma

abstract: A recent experimental observation has shown efficient transport of Mega Ampere of electron currents through aligned carbon nanotube arrays [Phys. Rev Letts. 108, 235005 (2012)]. The result was subsequently interpreted on the basis of suppression of the filamentation instability in an inhomogeneous plasma [Phys. Plasmas 21, 012108 (2014)]. This inhomogeneity forms as a result of the ionization of the carbon nanotubes. In the present work a full 1D3V Particle-in-Cell (PIC) simulations have been carried out for the propagation of relativistic electron beams (REB) through an inhomogeneous background plasma. The suppression of the filamentation instability, responsible for beam divergence, is shown. The simulation also confirms that in the nonlinear regime too the REB propagation is better when it propagates through a plasma whose density is inhomogeneous transverse to the beam. The role of inhomogeneity scale length, its amplitude and the transverse beam temperature etc., in the suppression of the instability is studied in detail.


title: Dissipative Abelian Sandpile Models

abstract: We introduce a family of abelian sandpile models with two parameters $n, m \in {\bf N}$ defined on finite lattices on $d$-dimensional torus. Sites with $2dn+m$ or more grains of sand are unstable and topple, and in each toppling $m$ grains dissipate from the system. Because of dissipation in bulk, the models are well-defined on the shift-invariant lattices and the infinite-volume limit of systems can be taken. From the determinantal expressions, we obtain the asymptotic forms of the avalanche propagators and the height-$(0,0)$ correlations of sandpiles for large distances in the infinite-volume limit in any dimensions $d \geq 2$. We show that both of them decay exponentially with the correlation length $$ \xi(d, a)=(\sqrt{d} \sinh^{-1} \sqrt{a(a+2)} \ )^{-1}, $$ if the dissipation rate $a =m/(2dn)$ is positive. By considering a series of models with increasing $n$, we discuss the limit $a \downarrow 0$ and the critical exponent defined by $\nu_{a}=- \lim_{a \downarrow 0} \log \xi(d, a)/ \log a$ is determined as $$ \nu_{a}=1/2 $$ for all $d \geq 2$. Comparison with the $q \downarrow 0$ limit of $q$-state Potts model in external magnetic field is discussed.

categories: math-ph cond-mat.stat-mech math.MP math.PR nlin.SI


title: Cracking a hierarchical chaotic image encryption algorithm based on permutation

abstract: In year 2000, an efficient hierarchical chaotic image encryption (HCIE) algorithm was proposed, which divides a plain-image of size $M\times N$ with $T$ possible value levels into $K$ blocks of the same size and then operates position permutation on two levels: intra-block and inter-block. As a typical position permutation-only encryption algorithm, it has received intensive attention. The present paper analyzes specific security performance of HCIE against ciphertext-only attack and known/chosen-plaintext attack. It is found that only $O(\lceil\log_T(M\cdot N/K) \rceil)$ known/chosen plain-images are sufficient to achieve a good performance, and the computational complexity is $O(M\cdot N\cdot \lceil\log_T(M\cdot N/K) \rceil)$, which effectively demonstrates that hierarchical permutation-only image encryption algorithms are less secure than normal (i.e., non-hierarchical) ones. Detailed experiment results are given to verify the feasibility of the known-plaintext attack. In addition, it is pointed out that the security of HCIE against ciphertext-only attack was much overestimated.


title: Is independent choice possible?

abstract: This paper questions the generally accepted assumption that one can make a random choice that is independent of the rest of the universe. We give a general description of any setup that could be conceived to generate random numbers. Based on the fact that the initial state of such setup together with its environment cannot be known, we show that the independence of its generated output cannot be guaranteed. Some consequences of this theoretical limitation are discussed.


title: On the radial abundance gradients in disks of irregular galaxies

abstract: We determine the radial abundance distributions across the disks of fourteen irregular galaxies of the types Sm and Im (morphological T types T = 9 and T =10) as traced by their HII regions. The oxygen and nitrogen abundances in HII regions are estimated through the Te method or/and with the counterpart method (C method). Moreover, we examine the correspondence between the radial abundance gradient and the surface brightness profile. We find that irregular galaxies with a flat inner profile (flat or outwardly increasing surface brightness in the central region) show shallow (if any) radial abundance gradients. On the other hand, irregular galaxies with a steep inner profile (with or without a bulge or central star cluster) usually show rather steep radial abundance gradients. This is in contrast to the widely held belief that irregular galaxies do not usually show a radial abundance gradient.


title: Cosmology and stability in scalar tensor bigravity with non-minimal kinetic coupling gravity

abstract: We generalize the scalar tensor bigravity models to the non-minimal kinetic coupling scalar tensor bigravity models with two scalar fields whose kinetic terms are non-minimally coupled to two Einstein tensors constructed by two metrics. We show that a broad class of expanding universes can be explained by some solutions of this model. Then, we study the stability issue of the solutions by means of imposing homogeneous perturbation on the equations of motion and extract the stable solutions.


title: Electronic Structure of Oxide Interfaces: A Comparative Analysis of GdTiO$_3$/SrTiO$_3$ and LaAlO$_3$/SrTiO$_3$ Interfaces

abstract: Emergent phases in the two-dimensional electron gas (2DEG) formed at the interface between two insulating oxides have attracted great attention in the past decade. We present ab-initio electronic structure calculations for the interface between a Mott insulator GdTiO$_3$ (GTO) and a band insulator SrTiO$_3$ (STO) and compare our results with those for the widely studied LaAlO$_3$/SrTiO$_3$ (LAO/STO) interface between two band insulators. Our GTO/STO results are in excellent agreement with experiments, but qualitatively different from LAO/STO. We find an interface carrier density of 0.5$e^{-}$/Ti, independent of GTO thickness in both superlattice and thin film geometries, in contrast to LAO/STO. The superlattice geometry in LAO/STO offers qualitatively the same result as in GTO/STO. On the other hand, for a thin film geometry, the interface carrier density builds up only beyond a threshold thickness of LAO. The positive charge at the vacuum surface that compensates the 2DEG at the interface also exhibits distinct behaviors in the two systems. The top GTO layer is found to be insulating due to correlation-driven charge disproportionation, while the top LAO layer is metallic within band theory and may become insulating due to surface disorder or surface reconstruction.


title: Exploring positron characteristics utilizing two new positron-electron correlation schemes based on three electronic-structure calculation methods

abstract: We make a gradient correction to a new local density approximation form of positron-electron correlation. Then the positron lifetimes and affinities are probed by using these two approximation forms based on three electronic-structure calculation methods including the full-potential linearized augmented plane wave (FLAPW) plus local orbitals approach, the atomic superposition (ATSUP) approach and the projector augmented wave (PAW) approach. The differences between calculated lifetimes using the FLAPW and ATSUP methods are clearly interpreted in the view of positron and electron transfers. We further find that a well implemented PAW method can give near-perfect agreement on both the positron lifetimes and affinities with the FLAPW method, and the competitiveness of the ATSUP method against the FLAPW/PAW method is reduced within the best calculations. By comparing with experimental data, the new introduced gradient corrected correlation form is proved competitive for positron lifetime and affinity calculations.


title: Structure in Dichotomous Preferences

abstract: Many hard computational social choice problems are known to become tractable when voters' preferences belong to a restricted domain, such as those of single-peaked or single-crossing preferences. However, to date, all algorithmic results of this type have been obtained for the setting where each voter's preference list is a total order of candidates. The goal of this paper is to extend this line of research to the setting where voters' preferences are dichotomous, i.e., each voter approves a subset of candidates and disapproves the remaining candidates. We propose several analogues of the notions of single-peaked and single-crossing preferences for dichotomous profiles and investigate the relationships among them. We then demonstrate that for some of these notions the respective restricted domains admit efficient algorithms for computationally hard approval-based multi-winner rules.


title: On the Stability of Einstein Static Universe in Doubly General Relativity Scenario

abstract: By presenting a relation between average energy of the ensemble of probe photons and energy density of the Universe, in the context of {\it gravity's rainbow} or {\it doubly general relativity} scenario, we introduce a rainbow FRW Universe model. By analyzing the fixed points in flat FRW model modified by two well known rainbow functions, we find that the finite time singularity avoidance (i.e. Big-Bang) may still remain as a problem. Then, we follow the "Emergent Universe" scenario in which there is no beginning of time and consequently there is no Big-Bang singularity. Moreover, we study the impact of a high energy quantum gravity modifications related to the gravity's rainbow on the stability conditions of an "Einstein static Universe" (ESU). We find that independent of a particular rainbow function, the positive energy condition dictates a positive spatial curvature for the Universe. In fact, without raising a nonphysical energy condition in the quantum gravity regimes, we can address an agreement between gravity's rainbow scenario and basic assumption of modern version of "Emergent Universe". We show that in the absence and presence of an energy-dependent cosmological constant $\Lambda(\epsilon)$, a stable Einstein static solution is available versus the homogeneous and linear scalar perturbations under the variety of obtained conditions. Also, we explore the stability of ESU against the vector and tensor perturbations.


title: A first-order logic for string diagrams

abstract: Equational reasoning with string diagrams provides an intuitive means of proving equations between morphisms in a symmetric monoidal category. This can be extended to proofs of infinite families of equations using a simple graphical syntax called !-box notation. While this does greatly increase the proving power of string diagrams, previous attempts to go beyond equational reasoning have been largely ad hoc, owing to the lack of a suitable logical framework for diagrammatic proofs involving !-boxes. In this paper, we extend equational reasoning with !-boxes to a fully-fledged first order logic called with conjunction, implication, and universal quantification over !-boxes. This logic, called !L, is then rich enough to properly formalise an induction principle for !-boxes. We then build a standard model for !L and give an example proof of a theorem for non-commutative bialgebras using !L, which is unobtainable by equational reasoning alone.


title: Fireflies: New software for interactively exploring dynamical systems using GPU computing

abstract: In non-linear systems, where explicit analytic solutions usually can't be found, visualisation is a powerful approach which can give insights into the dynamical behaviour of models; it is also crucial for teaching this area of mathematics. In this paper we present new software, Fireflies, which exploits the power of graphical processing unit (GPU) computing to produce spectacular interactive visualisations of arbitrary systems of ordinary differential equations. In contrast to typical phase portraits, Fireflies draws the current position of trajectories (projected onto 2D or 3D space) as single points of light, which move as the system is simulated. Due to the massively parallel nature of GPU hardware, Fireflies is able to simulate millions of trajectories in parallel (even on standard desktop computer hardware), producing "swarms" of particles that move around the screen in real-time according to the equations of the system. Particles that move forwards in time reveal stable attractors (e.g. fixed points and limit cycles), while the option of integrating another group of trajectories backwards in time can reveal unstable objects (repellers). Fireflies allows the user to change the parameters of the system as it is running, in order to see the effect that they have on the dynamics and to observe bifurcations. We demonstrate the capabilities of the software with three examples: a two-dimensional "mean field" model of neuronal activity, the classical Lorenz system, and a 15-dimensional model of three interacting biologically realistic neurons.


title: Possible atomic structures for the sub-bandgap absorption of chalcogen hyperdoped silicon

abstract: Single-crystal silicon wafers were hyperdoped respectively by sulfur, selenium, and tellurium element using ion implantation and nanosecond laser melting. The hyperdoping of such chalcogen elements endowed the treated silicon with a strong and wide sub-bandgap light absorptance. When these hyperdoped silicons were thermally annealed even at low temperatures (such as 200~400 oC), however, this extra sub-bandgap absorptance began to attenuate. In order to explain this attenuation of absorptance, alternatively, we consider it corresponding to a chemical decomposition reaction from optically absorbing structure to non-absorbing structure, and obtain a very good fitting to the attenuated absorptances by using Arrhenius equation. Further, we extract the reaction activation energies from the fittings and they are 0.343(+/- 0.031) eV for S-, 0.426(+/-0.042) eV for Se-, and 0.317(+/-0.033) eV for Te-hyperdoped silicon, respectively. We discuss these activation energies in term of the bond energies of chalcogen-Si metastable bonds, and finally suggest that several high-energy interstitial sites instead of the substitutional site, are very possibly the atomic structures that are responsible for the sub-bandgap absorptance of chalcogen hyperdoped silicon.


title: Power Allocation and Measurement Matrix Design for Block CS-Based Distributed MIMO Radars

abstract: Multiple-input multiple-output (MIMO) radars offer higher resolution, better target detection, and more accurate target parameter estimation. Due to the sparsity of the targets in space-velocity domain, we can exploit Compressive Sensing (CS) to improve the performance of MIMO radars when the sampling rate is much less than the Nyquist rate. In distributed MIMO radars, block CS methods can be used instead of classical CS ones for more performance improvement, because the received signal in this group of MIMO radars is a block sparse signal in a basis. In this paper, two new methods are proposed to improve the performance of the block CS-based distributed MIMO radars. The first one is a new method for optimal energy allocation to the transmitters, and the other one is a new method for optimal design of the measurement matrix. These methods are based on the minimization of an upper bound of the sensing matrix block-coherence. Simulation results show an increase in the accuracy of multiple targets parameters estimation for both proposed methods.


title: Saturated Sodium Chloride Solution under an External Static Electric Field: a Molecular Dynamics Study

abstract: The behavior of saturated aqueous sodium chloride solutions under a constant external electric field (E) was studied by molecular dynamics (MD) simulation. Our dynamic MD simulations have indicated that the irreversible nucleation process towards crystallization is accelerated by a moderate E, but retarded or even prohibited under a stronger E, which can be understood by the competition between self-diffusion and drift motion. The former increases with E resulting in the acceleration of the nucleation process, and the latter tears oppositely charged ions more apart under a stronger E leading to the deceleration of nucleation. Moreover, our steady-state MD simulations have indicated that a first-order phase transition happens in saturated solutions only when the applied E is below a certain threshold Ec, and the ratio of crystallized ions does not change with the electric field. The magnitude of Ec increases with concentration, because larger clusters are easy to form in a more concentrated solution and require a stronger E to dissociate them.


title: Discrete Heisenberg group and its automorphism group

abstract: In this note we give more easy and short proof of a statement previously proved by P. Kahn that the automorphism group of the discrete Heisenberg group ${\rm Heis}(3, \mathbb{Z}) $ is isomorphic to the group $ (\mathbb{Z} \oplus \mathbb{Z}) \rtimes GL(2,\mathbb{Z})$. The method which we suggest to construct this isomorphism gives far more transparent picture of the structure of the automorphism group of the group ${\rm Heis}(3, \mathbb{Z}) $.


title: AC Stark effect in ThO $H^3\Delta_1$ for the electron EDM search

abstract: A method and code for calculations of diatomic molecules in the external variable electromagnetic field have been developed. Code applied for calculation of systematics in the electron's electric dipole moment search experiment on ThO $H^3\Delta_1$ state related to geometric phases, including dependence on $\Omega$-doublet, rotational level, and external static electric field. It is found that systematics decrease cubically with respect to the frequency of the rotating transverse component of the electric field. Calculation confirms that experiment on ThO $H^3\Delta_1$ state is very robust against systematic errors related to geometric phases.


title: Machine Learning for Quantum Mechanical Properties of Atoms in Molecules

abstract: We introduce machine learning models of quantum mechanical observables of atoms in molecules. Instant out-of-sample predictions for proton and carbon nuclear chemical shifts, atomic core level excitations, and forces on atoms reach accuracies on par with density functional theory reference. Locality is exploited within non-linear regression via local atom-centered coordinate systems. The approach is validated on a diverse set of 9k small organic molecules. Linear scaling of computational cost in system size is demonstrated for saturated polymers with up to sub-mesoscale lengths.


title: Slow Limb Movements Require Precise Control of Muscle Stiffness

abstract: Slow and accurate finger and limb movements are essential to daily activities, but their neural control and governing mechanics are relatively unexplored. We consider neuromechanical systems where slow movements are produced by neural commands that modulate muscle stiffness. This formulation based on strain-energy equilibria is in agreement with prior work on neural control of muscle and limb impedance. Slow limb movements are driftless in the sense that movement stops when neural commands stop. We demonstrate, in the context of two planar tendon-driven systems representing a finger and a leg, that the control of muscle stiffness suffices to produce stable and accurate limb postures and quasi-static (slow) transitions among them. We prove, however, that stable postures are achievable only when muscles are pre-tensioned, as is the case for natural muscle tone. Our results further indicate, in accordance with experimental findings, that slow movements are non-smooth. The non-smoothness arises because the precision with which individual muscle stiffnesses need to be controlled changes substantially throughout the limb's motion. These results underscore the fundamental roles of muscle tone and accurate neural control of muscle stiffness in producing stable limb postures and slow movements.


title: Volume and lattice points counting for the cyclopermutohedron

abstract: The face lattice of the permutohedron realizes the combinatorics of linearly ordered partitions of the set $[n]=\{1,...,n\}$. Similarly, the cyclopermutohedron is a virtual polytope that realizes the combinatorics of cyclically ordered partitions of $[n]$. It is known that the volume of the standard permutohedron equals the number of trees with $n$ labeled vertices multiplied by $\sqrt{n}$. The number of integer points of the standard permutohedron equals the number of forests on $n$ labeled vertices. In the paper we prove that the volume of the cyclopermutohedron also equals some weighted number of forests, which eventually reduces to zero. We also derive a combinatorial formula for the number of integer points in the cyclopermutohedron. Another object of the paper is the configuration space of a polygonal linkage $L$. It has a cell decomposition $\mathcal{K}(L)$ related to the face lattice of cyclopermutohedron. Using this relationship, we introduce and compute the volume $Vol(\mathcal{K}(L))$.


title: Joint Multi-Leaf Segmentation, Alignment and Tracking from Fluorescence Plant Videos

abstract: This paper proposes a novel framework for fluorescence plant video processing. The plant research community is interested in the leaf-level photosynthetic analysis within a plant. A prerequisite for such analysis is to segment all leaves, estimate their structures, and track them over time. We identify this as a joint multi-leaf segmentation, alignment, and tracking problem. First, leaf segmentation and alignment are applied on the last frame of a plant video to find a number of well-aligned leaf candidates. Second, leaf tracking is applied on the remaining frames with leaf candidate transformation from the previous frame. We form two optimization problems with shared terms in their objective functions for leaf alignment and tracking respectively. A quantitative evaluation framework is formulated to evaluate the performance of our algorithm with four metrics. Two models are learned to predict the alignment accuracy and detect tracking failure respectively in order to provide guidance for subsequent plant biology analysis. The limitation of our algorithm is also studied. Experimental results show the effectiveness, efficiency, and robustness of the proposed method.