title: Long-Lived Non-Equilibrium Interstitial-Solid-Solutions in Binary Mixtures

abstract: We perform particle resolved experimental studies on the heterogeneous crystallisation process of two compo- nent mixtures of hard spheres. The components have a size ratio of 0.39. We compared these with molecular dynamics simulations of homogenous nucleation. We find for both experiments and simulations that the final assemblies are interstitial solid solutions, where the large particles form crystalline close-packed lattices, whereas the small particles occupy random interstitial sites. This interstitial solution resembles that found at equilibrium when the size ratios are 0.3 [Filion et al., Phys. Rev. Lett. 107, 168302 (2011)] and 0.4 [Filion, PhD Thesis, Utrecht University (2011)]. However, unlike these previous studies, for our system sim- ulations showed that the small particles are trapped in the octahedral holes of the ordered structure formed by the large particles, leading to long-lived non-equilibrium structures in the time scales studied and not the equilibrium interstitial solutions found earlier. Interestingly, the percentage of small particles in the crystal formed by the large ones rapidly reaches a maximum of around 14% for most of the packing fractions tested, unlike previous predictions where the occupancy of the interstitial sites increases with the system concentration. Finally, no further hopping of the small particles was observed.


title: Generalized Stieltjes constants and integrals involving the log-log function: Kummer's Theorem in action

abstract: In this note, we recall Kummer's Fourier series expansion of the 1-periodic function that coincides with the logarithm of the Gamma function on the unit interval $(0,1)$, and we use it to find closed forms for some numerical series related to the generalized Stieltjes constants, and some integrals involving the function $x\mapsto \ln \ln(1/x)$.


title: Simple modules in the Auslander-Reiten quiver of principal blocks with abelian defect groups

abstract: Given an odd prime $p$, we investigate the position of simple modules in the stable Auslander-Reiten quiver of the principal block of a finite group with non-cyclic abelian Sylow $p$-subgroups. In particular, we prove a reduction to finite simple groups. In the case that the characteristic is $3$, we prove that simple modules in the principal block all lie at the end of their components


title: Toward Abstraction from Multi-modal Data: Empirical Studies on Multiple Time-scale Recurrent Models

abstract: The abstraction tasks are challenging for multi- modal sequences as they require a deeper semantic understanding and a novel text generation for the data. Although the recurrent neural networks (RNN) can be used to model the context of the time-sequences, in most cases the long-term dependencies of multi-modal data make the back-propagation through time training of RNN tend to vanish in the time domain. Recently, inspired from Multiple Time-scale Recurrent Neural Network (MTRNN), an extension of Gated Recurrent Unit (GRU), called Multiple Time-scale Gated Recurrent Unit (MTGRU), has been proposed to learn the long-term dependencies in natural language processing. Particularly it is also able to accomplish the abstraction task for paragraphs given that the time constants are well defined. In this paper, we compare the MTRNN and MTGRU in terms of its learning performances as well as their abstraction representation on higher level (with a slower neural activation). This was done by conducting two studies based on a smaller data- set (two-dimension time sequences from non-linear functions) and a relatively large data-set (43-dimension time sequences from iCub manipulation tasks with multi-modal data). We conclude that gated recurrent mechanisms may be necessary for learning long-term dependencies in large dimension multi-modal data-sets (e.g. learning of robot manipulation), even when natural language commands was not involved. But for smaller learning tasks with simple time-sequences, generic version of recurrent models, such as MTRNN, were sufficient to accomplish the abstraction task.


title: An infinitely differentiable function with compact support: Definition and properties

abstract: This is the English translation of my old paper 'Definici\'on y estudio de una funci\'on indefinidamente diferenciable de soporte compacto', Rev. Real Acad. Ciencias 76 (1982) 21-38. In it a function (essentially Fabius function) is defined and given its main properties, including: unicity, interpretation as a probability, partition of unity with its translates, formulas for its $n$-th derivates, rationality of its values at dyadic points, formulas for the effective computation of these values, and some arithmetical properties of these values. Since I need it now for a reference, I have translated it.


title: How close are the eigenvectors and eigenvalues of the sample and actual covariance matrices?

abstract: How many samples are sufficient to guarantee that the eigenvectors and eigenvalues of the sample covariance matrix are close to those of the actual covariance matrix? For a wide family of distributions, including distributions with finite second moment and distributions supported in a centered Euclidean ball, we prove that the inner product between eigenvectors of the sample and actual covariance matrices decreases proportionally to the respective eigenvalue distance. Our findings imply non-asymptotic concentration bounds for eigenvectors, eigenspaces, and eigenvalues. They also provide conditions for distinguishing principal components based on a constant number of samples.


title: A catalog of optical to X-ray spectral energy distributions of z~2 quasars observed with Swift. I: First results

abstract: We present the Swift optical to X-ray spectral energy distributions (SEDs) of 44 quasars at redshifts z~2 observed by Swift, part of a larger program to establish and characterize the optical through X-ray SEDs of moderate-redshift quasars. Here we outline our analysis approach and present preliminary analysis and results for the first third of the full quasar sample. Not all quasars in the sample are detected in X-rays; all of the X-ray detected objects so far are radio loud. As expected for radio loud objects, they are X-ray bright relative to radio-quiet quasars of comparable optical luminosities, with an average alpha_ox = 1.39 +/- 0.03 (where alpha_ox is the power-law slope connecting the monochromatic flux at 2500 Ang and at 2 keV), and display hard X-ray spectra. We find integrated 3000 Ang - 25 keV accretion luminosities of between 0.7*10^(46) erg s^(-1) and 5.2*10^(47) erg s^(-1). Based on single-epoch spectroscopic virial black hole mass estimates, we find that these quasars are accreting at substantial Eddington fractions, 0.1 \le L/LEdd \le 1.


title: Ergodicity in umbrella billiards

abstract: We investigate a three-parameter family of billiard tables with circular arc boundaries. These umbrella-shaped billiards may be viewed as a generalization of two-parameter moon and asymmetric lemon billiards, in which the latter classes comprise instances where the new parameter is $0$. Like those two previously studied classes, for certain parameters umbrella billiards exhibit evidence of chaotic behavior despite failing to meet certain criteria for defocusing or dispersing, the two most well understood mechanisms for generating ergodicity and hyperbolicity. For some parameters corresponding to non-ergodic lemon and moon billiards, small increases in the new parameter transform elliptic $2$-periodic points into a cascade of higher order elliptic points. These may either stabilize or dissipate as the new parameter is increased. We characterize the periodic points and present evidence of new ergodic examples.


title: Network Flow Based Post Processing for Sales Diversity

abstract: Collaborative filtering is a broad and powerful framework for building recommendation systems that has seen widespread adoption. Over the past decade, the propensity of such systems for favoring popular products and thus creating echo chambers have been observed. This has given rise to an active area of research that seeks to diversify recommendations generated by such algorithms. We address the problem of increasing diversity in recommendation systems that are based on collaborative filtering that use past ratings to predicting a rating quality for potential recommendations. Following our earlier work, we formulate recommendation system design as a subgraph selection problem from a candidate super-graph of potential recommendations where both diversity and rating quality are explicitly optimized: (1) On the modeling side, we define a new flexible notion of diversity that allows a system designer to prescribe the number of recommendations each item should receive, and smoothly penalizes deviations from this distribution. (2) On the algorithmic side, we show that minimum-cost network flow methods yield fast algorithms in theory and practice for designing recommendation subgraphs that optimize this notion of diversity. (3) On the empirical side, we show the effectiveness of our new model and method to increase diversity while maintaining high rating quality in standard rating data sets from Netflix and MovieLens.


title: Counting edge-injective homomorphisms and matchings on restricted graph classes

abstract: We consider the $\#\mathsf{W}[1]$-hard problem of counting all matchings with exactly $k$ edges in a given input graph $G$; we prove that it remains $\#\mathsf{W}[1]$-hard on graphs $G$ that are line graphs or bipartite graphs with degree $2$ on one side. In our proofs, we use that $k$-matchings in line graphs can be equivalently viewed as edge-injective homomorphisms from the disjoint union of $k$ length-$2$ paths into (arbitrary) host graphs. Here, a homomorphism from $H$ to $G$ is edge-injective if it maps any two distinct edges of $H$ to distinct edges in $G$. We show that edge-injective homomorphisms from a pattern graph $H$ can be counted in polynomial time if $H$ has bounded vertex-cover number after removing isolated edges. For hereditary classes $\mathcal{H}$ of pattern graphs, we complement this result: If the graphs in $\mathcal{H}$ have unbounded vertex-cover number even after deleting isolated edges, then counting edge-injective homomorphisms with patterns from $\mathcal{H}$ is $\#\mathsf{W}[1]$-hard. Our proofs rely on an edge-colored variant of Holant problems and a delicate interpolation argument; both may be of independent interest.


title: Learning to Detect Human-Object Interactions

abstract: We study the problem of detecting human-object interactions (HOI) in static images, defined as predicting a human and an object bounding box with an interaction class label that connects them. HOI detection is a fundamental problem in computer vision as it provides semantic information about the interactions among the detected objects. We introduce HICO-DET, a new large benchmark for HOI detection, by augmenting the current HICO classification benchmark with instance annotations. To solve the task, we propose Human-Object Region-based Convolutional Neural Networks (HO-RCNN). At the core of our HO-RCNN is the Interaction Pattern, a novel DNN input that characterizes the spatial relations between two bounding boxes. Experiments on HICO-DET demonstrate that our HO-RCNN, by exploiting human-object spatial relations through Interaction Patterns, significantly improves the performance of HOI detection over baseline approaches.


title: IGR J17062-6143 is an Accreting Millisecond X-ray Pulsar

abstract: We present the discovery of 163.65 Hz X-ray pulsations from IGR J17062-6143 in the only observation obtained from the source with the Rossi X-ray Timing Explorer. This detection makes IGR J17062-6143 the lowest-frequency accreting millisecond X-ray pulsar presently known. The pulsations are detected in the 2 - 12 keV band with an overall significance of 4.3 sigma, and an observed pulsed amplitude of 5.54 +- 0.67 % (in this band). Both dynamic power spectral and coherent phase timing analysis indicate that the pulsation frequency is decreasing during the 1.2 ks observation in a manner consistent with orbital motion of the neutron star. Because the observation interval is short, we cannot precisely measure the orbital period; however, periods shorter than 17 minutes are excluded at 90 % confidence. For the range of acceptable circular orbits the inferred binary mass function substantially overlaps the observed range for the AMXP population as a whole.


title: Work drives time evolution

abstract: We propose the idea that time evolution of quantum systems is driven by work. The formalism presented here falls within the scope of a recently proposed theory of gravitating quantum matter where extractible work, and not energy, is responsible for gravitation. Our main assumption is that extractible work, and not the Hamiltonian, dictates dynamics. We find that expectation values of meaningful quantities, such as the occupation number, deviate from those predicted by standard quantum mechanics. The scope, applications and validity of this proposal are also discussed.


title: Stress Response of Granular Systems

abstract: We develop a framework for stress response in two dimensional granular media, with and without friction, that respects vector force balance at the microscopic level. We introduce local gauge degrees of freedom that determine the response of contact forces between constituent grains on a given, disordered, contact network, to external perturbations. By mapping this response to the spectral properties of the graph Laplacian corresponding to the underlying contact network, we show that this naturally leads to spatial localization of forces. We present numerical evidence for localization using exact diagonalization studies of network Laplacians of soft disk packings. Finally, we discuss the role of other constraints, such as torque balance, in determining the stability of a granular packing to external perturbations.


title: Random walks on Ramanujan complexes and digraphs

abstract: The cutoff phenomenon was recently confirmed for random walks on Ramanujan graphs by the first author and Peres. In this work, we obtain analogs in higher dimensions, for random walk operators on any Ramanujan complex associated with a simple group $G$ over a local field $F$. We show that if $T$ is any $k$-regular $G$-equivariant operator on the Bruhat-Tits building with a simple combinatorial property (collision-free), the associated random walk on the $n$-vertex Ramanujan complex has cutoff at time $\log_k n$. The high dimensional case, unlike that of graphs, requires tools from non-commutative harmonic analysis and the infinite-dimensional representation theory of $G$. Via these, we show that operators $T$ as above on Ramanujan complexes give rise to Ramanujan digraphs with a special property ($r$-normal), implying cutoff. Applications include geodesic flow operators, geometric implications, and a confirmation of the Riemann Hypothesis for the associated zeta functions over every group $G$, previously known for groups of type $\widetilde A_n$ and $\widetilde C_2$.


title: On the Four Dimensional Holoraumy of the 4D, $\cal N$ = 1 Complex Linear Supermultiplet

abstract: We present arguments to support the existence of weight spaces for supersymmetric field theories and identify the calculations of information about supermultiplets to define such spaces via the concept of "holoraumy." For the first time this is extended to the complex linear supermultiplet by a calculation of the commutator of supercovariant derivatives on all of its component fields.


title: Cache-Aided Content Delivery over Erasure Broadcast Channels

abstract: A cache-aided broadcast network is studied, in which a server delivers contents to a group of receivers over a packet erasure broadcast channel (BC). The receivers are divided into two sets with regards to their channel qualities: the weak and strong receivers, where all the weak receivers have statistically worse channel qualities than all the strong receivers. The weak receivers, in order to compensate for the high erasure probability they encounter over the channel, are equipped with cache memories of equal size, while the receivers in the strong set have no caches. Data can be pre-delivered to weak receivers' caches over the off-peak traffic period before the receivers reveal their demands. Allowing arbitrary erasure probabilities for the weak and strong receivers, a joint caching and channel coding scheme, which divides each file into several subfiles, and applies a different caching and delivery scheme for each subfile, is proposed. It is shown that all the receivers, even those without any cache memories, benefit from the presence of caches across the network. An information theoretic trade-off between the cache size and the achievable rate is formulated. It is shown that the proposed scheme improves upon the state-of-the-art in terms of the achievable trade-off.


title: Improving the upper bound on the length of the shortest reset words

abstract: We improve the best known upper bound on the length of the shortest reset words of synchronizing automata. The new bound is slightly better than $114 n^3 / 685 + O(n^2)$. The \v{C}ern\'y conjecture states that $(n-1)^2$ is an upper bound. So far, the best general upper bound was $(n^3-n)/6-1$ obtained by J.-E.~Pin and P.~Frankl in 1982. Despite a number of efforts, it remained unchanged for about 35 years. To obtain the new upper bound we utilize avoiding words. A word is avoiding for a state $q$ if after reading the word the automaton cannot be in $q$. We obtain upper bounds on the length of the shortest avoiding words, and using the approach of Trahtman from 2011 combined with the well known Frankl theorem from 1982, we improve the general upper bound on the length of the shortest reset words. For all the bounds, there exist polynomial algorithms finding a word of length not exceeding the bound.


title: LCL problems on grids

abstract: LCLs or locally checkable labelling problems (e.g. maximal independent set, maximal matching, and vertex colouring) in the LOCAL model of computation are very well-understood in cycles (toroidal 1-dimensional grids): every problem has a complexity of $O(1)$, $\Theta(\log^* n)$, or $\Theta(n)$, and the design of optimal algorithms can be fully automated. This work develops the complexity theory of LCL problems for toroidal 2-dimensional grids. The complexity classes are the same as in the 1-dimensional case: $O(1)$, $\Theta(\log^* n)$, and $\Theta(n)$. However, given an LCL problem it is undecidable whether its complexity is $\Theta(\log^* n)$ or $\Theta(n)$ in 2-dimensional grids. Nevertheless, if we correctly guess that the complexity of a problem is $\Theta(\log^* n)$, we can completely automate the design of optimal algorithms. For any problem we can find an algorithm that is of a normal form $A' \circ S_k$, where $A'$ is a finite function, $S_k$ is an algorithm for finding a maximal independent set in $k$th power of the grid, and $k$ is a constant. Finally, partially with the help of automated design tools, we classify the complexity of several concrete LCL problems related to colourings and orientations.


title: Symmetric bi-skew maps and symmetrized motion planning in projective spaces

abstract: This work is motivated by the question of whether there are spaces $X$ for which the Farber-Grant symmetric topological complexity $TC^S(X)$ differs from the Basabe-Gonz\'alez-Rudyak-Tamaki symmetric topological complexity $TC^\Sigma(X)$. It is known that, for a projective space $RP^m$, $TC^S(RP^m)$ captures, with a few potentially exceptional cases, the Euclidean embedding dimension of $RP^m$. We now show that, for all $m\geq1$, $TC^\Sigma(RP^m)$ is characterized as the smallest positive integer $n$ for which there is a symmetric $\mathbb{Z}_2$-biequivariant map $S^m\times S^m\to S^n$ with a "monoidal" behavior on the diagonal. This result thus lies at the core of the efforts in the 1970's to characterize the embedding dimension of real projective spaces in terms of the existence of symmetric axial maps. Together with Nakaoka's description of the cohomology ring of symmetric squares, this allows us to compute both $TC$ numbers in the case of $RP^{2^e}$ for $e\geq1$. In particular, this leaves the torus $S^1\times S^1$ as the only closed surface whose symmetric (symmetrized) $TC^S$ ($TC^\Sigma$) -invariant is currently unknown.


title: Luminous, pc-scale CO 6-5 emission in the obscured nucleus of NGC1377

abstract: High resolution submm observations are important in probing the morphology, column density and dynamics of obscured active galactic nuclei (AGNs). With high resolution (0.06 x 0.05) ALMA 690 GHz observations we have found bright (TB >80 K) and compact (FWHM 10x7 pc) CO 6-5 line emission in the nucleus of the extremely radio-quiet galaxy NGC1377. The integrated CO 6-5 intensity is aligned with the previously discovered jet/outflow of NGC1377 and is tracing the dense (n>1e4 cm-3), hot gas at the base of the outflow. The velocity structure is complex and shifts across the jet/outflow are discussed in terms of jet-rotation or separate, overlapping kinematical components. High velocity gas (deltaV +-145 km/s) is detected inside r<2-3 pc and we suggest that it is emerging from an inclined rotating disk or torus of position angle PA=140+-20 deg with a dynamical mass of approx 3e6 Msun. This mass is consistent with that of a supermassive black hole (SMBH), as inferred from the M-sigma relation. The gas mass of the proposed disk/torus constitutes <3% of the nuclear dynamical mass. In contrast to the intense CO 6-5 line emission, we do not detect dust continuum with an upper limit of S(690GHz)<2mJy. The corresponding, 5 pc, H2 column density is estimated to N(H2)<3e23 cm-2, which is inconsistent with a Compton Thick (CT) source. We discuss the possibility that CT obscuration may be occuring on small (subparsec) or larger scales. From SED fitting we suggest that half of the IR emission of NGC1377 is nuclear and the rest (mostly the far-infrared (FIR)) is more extended. The extreme radio quietness, and the lack of emission from other star formation tracers, raise questions on the origin of the FIR emission. We discuss the possibility that it is arising from the dissipation of shocks in the molecular jet/outflow or from irradiation by the nuclear source along the poles.


title: Communication Reducing Algorithms for Distributed Hierarchical N-Body Problems with Boundary Distributions

abstract: Reduction of communication and efficient partitioning are key issues for achieving scalability in hierarchical $N$-Body algorithms like FMM. In the present work, we propose four independent strategies to improve partitioning and reduce communication. First of all, we show that the conventional wisdom of using space-filling curve partitioning may not work well for boundary integral problems, which constitute about 50% of FMM's application user base. We propose an alternative method which modifies orthogonal recursive bisection to solve the cell-partition misalignment that has kept it from scaling previously. Secondly, we optimize the granularity of communication to find the optimal balance between a bulk-synchronous collective communication of the local essential tree and an RDMA per task per cell. Finally, we take the dynamic sparse data exchange proposed by Hoefler et al. and extend it to a hierarchical sparse data exchange, which is demonstrated at scale to be faster than the MPI library's MPI_Alltoallv that is commonly used.


title: Shifts and Splittings of the Hole Bands in the Nematic Phase of FeSe

abstract: We report a high-resolution laser-based angle-resolved photoemission spectroscopy (laser-ARPES) study of single crystals of FeSe, focusing on the temperature-dependence of the hole-like bands around the ${\rm \Gamma}$ point. As the system cools through the tetragonal-orthorhombic "nematic" structural transition at 90~K, the splitting of the $d_{xz}$/$d_{yz}$ bands is observed to increase by a magnitude of 13 meV. Moreover, the onset of a $\sim$10 meV downward shift of the $d_{xy}$ band is also at 90~K. These measurements provide clarity on the nature, magnitude and temperature-dependence of the band shifts at the ${\rm \Gamma}$ point in the nematic phase of FeSe.


title: Quantum interference effects in resonant Raman spectroscopy of single- and triple-layer MoTe$_2$ from first principles

abstract: We present a combined experimental and theoretical study of resonant Raman spectroscopy in single- and triple-layer MoTe$_2$. Raman intensities are computed entirely from first principles by calculating finite differences of the dielectric susceptibility. In our analysis, we investigate the role of quantum interference effects and the electron-phonon coupling. With this method, we explain the experimentally observed intensity inversion of the $A^\prime_1$ vibrational modes in triple-layer MoTe2 with increasing laser photon energy. Finally, we show that a quantitative comparison with experimental data requires the proper inclusion of excitonic effects.


title: Objective Bayesian Analysis for Change Point Problems

abstract: In this paper we present a loss-based approach to change point analysis. In particular, we look at the problem from two perspectives. The first focuses on the definition of a prior when the number of change points is known a priori. The second contribution aims to estimate the number of change points by using a loss-based approach recently introduced in the literature. The latter considers change point estimation as a model selection exercise. We show the performance of the proposed approach on simulated data and real data sets.

categories: stat.ME math.ST stat.AP stat.CO stat.ML stat.TH


title: A Performance Analysis of Simple Runtime System for Actor Programming in C++

abstract: In this paper, we propose the Templet -- a runtime system for actor programming of high performance computing in C++. We provide a compact source code of the runtime system, which uses standard library of C++11 only. We demonstrate how it differs from the classic implementations of the actor model. The practical significance of the Templet was examined by comparative study on the performance of three applications: the reference code in C++, managed by the OpenMP; the actor code in C++, managed by the Templet; the actor code in Java, managed by the Akka. As a test problem we used a numerical algorithm for solving the heat equation.


abstract: Adversarial learning methods are a promising approach to training robust deep networks, and can generate complex samples across diverse domains. They also can improve recognition despite the presence of domain shift or dataset bias: several adversarial approaches to unsupervised domain adaptation have recently been introduced, which reduce the difference between the training and test domain distributions and thus improve generalization performance. Prior generative approaches show compelling visualizations, but are not optimal on discriminative tasks and can be limited to smaller shifts. Prior discriminative approaches could handle larger domain shifts, but imposed tied weights on the model and did not exploit a GAN-based loss. We first outline a novel generalized framework for adversarial adaptation, which subsumes recent state-of-the-art approaches as special cases, and we use this generalized view to better relate the prior approaches. We propose a previously unexplored instance of our general framework which combines discriminative modeling, untied weight sharing, and a GAN loss, which we call Adversarial Discriminative Domain Adaptation (ADDA). We show that ADDA is more effective yet considerably simpler than competing domain-adversarial methods, and demonstrate the promise of our approach by exceeding state-of-the-art unsupervised adaptation results on standard cross-domain digit classification tasks and a new more difficult cross-modality object classification task.


title: A Large Systematic Search for Close Supermassive Binary and Rapidly Recoiling Black Holes - III. Radial Velocity Variations

abstract: We have been spectroscopically monitoring 88 quasars selected to have broad H$\beta$ emission lines offset from their systemic redshift by thousands of km s$^{-1}$. By analogy with single-lined spectroscopic binary stars, we consider these quasars to be candidates for hosting supermassive black hole binaries (SBHBs). In this work we present new radial velocity measurements, typically 3-4 per object over a time period of up to 12 years in the observer's frame. In 29/88 of the SBHB candidates no variability of the shape of the broad H$\beta$ profile is observed, which allows us to make reliable measurements of radial velocity changes. Among these, we identify three objects that have displayed systematic and monotonic velocity changes by several hundred km s$^{-1}$ and are prime targets for further monitoring. Because the periods of the hypothetical binaries are expected to be long, we cannot hope to observe many orbital cycles during our lifetimes. Instead, we seek to evaluate the credentials of the SBHB candidates by attempting to rule out the SBHB hypothesis. In this spirit, we present a method for placing a lower limit on the period, and thus the mass, of the SBHBs under the assumption that the velocity changes we observe are due to orbital motion. Given the duration of our monitoring campaign and the uncertainties in the radial velocities, we were able to place a lower limit on the total mass in the range $4.7\times10^4-3.8\times10^8$ $M_{\scriptscriptstyle \odot}$, which does not yet allow us to rule out the SBHB hypothesis for any candidates.


title: On affine Tverberg-type results without continuous generalization

abstract: Recent progress building on the groundbreaking work of Mabillard and Wagner has shown that there are important differences between the affine and continuous theory for Tverberg-type results. These results aim to describe the intersection pattern of convex hulls of point sets in Euclidean space and continuous relaxations thereof. Here we give additional examples of an affine-continuous divide, but our deductions are almost elementary and do not build on the technical work of Mabillard and Wagner. Moreover, these examples show a difference between the affine and continuous theory even asymptotically for arbitrarily large complexes. Along the way we settle the Tverberg admissible-prescribable problem (or AP conjecture) in the negative, give a new, short and elementary proof of the balanced case of the AP conjecture which was recently proven by Joji\'c, Vre\'cica, and \v{Z}ivaljevi\'c in a series of two papers, provide examples of Tverberg-type results that hold affinely but not continuously without divisibility conditions on the intersection multiplicity, extend a result of Sober\'on, and show that this extension has a topological generalization if and only if the intersection multiplicity is a prime.


title: Smoothing closed gridded surfaces embedded in ${\mathbb R}^4$

abstract: We say that a topological $n$-manifold $N$ is a cubical $n$-manifold if it is contained in the $n$-skeleton of the canonical cubulation $\mathcal{C}$ of ${\mathbb{R}}^{n+k}$ ($k\geq1$). In this paper, we prove that any closed, oriented cubical $2$-manifold has a transverse field of 2-planes in the sense of Whitehead and therefore it is smoothable by a small ambient isotopy.


title: Rigorous bounds on the stationary distributions of the chemical master equation via mathematical programming

abstract: The stochastic dynamics of biochemical networks are usually modelled with the chemical master equation (CME). The stationary distributions of CMEs are seldom solvable analytically, and numerical methods typically produce estimates with uncontrolled errors. Here, we introduce mathematical programming approaches that yield approximations of these distributions with computable error bounds which enable the verification of their accuracy. First, we use semidefinite programming to compute increasingly tighter upper and lower bounds on the moments of the stationary distributions for networks with rational propensities. Second, we use these moment bounds to formulate linear programs that yield convergent upper and lower bounds on the stationary distributions themselves, their marginals and stationary averages. The bounds obtained also provide a computational test for the uniqueness of the distribution. In the unique case, the bounds form an approximation of the stationary distribution with a computable bound on its error. In the non-unique case, our approach yields converging approximations of the ergodic distributions. We illustrate our methodology through several biochemical examples taken from the literature: Schl\"ogl's model for a chemical bifurcation, a two-dimensional toggle switch, a model for bursty gene expression, and a dimerisation model with multiple stationary distributions.

categories: math.PR math.OC q-bio.MN q-bio.PE q-bio.QM


title: On biconservative hypersurfaces in 4-dimensional Riemannian space forms

abstract: In this paper, we study biconservative hypersurfaces in $\mathbb S^{n}$ and $\mathbb H^{n}$. Further, we obtain complete explicit classification of biconservative hypersurfaces in $4$-dimensional Riemannian space form with exactly three distinct principal curvatures.


title: Wilsonian RG-flow approach to holographic transport with momentum dissipation

abstract: We systematically present a new approach for studying the coupled linear transport of holographic systems. In this approach, the set of equations for the linear perturbations can be reduced to a first-order nonlinear ordinary differential equation expressed as the radial (renormalization group) flow equation of the transport matrices. As an important application, we use this approach to compute the DC and AC conductivities of a holographic model with momentum dissipation, which can be easily read off from the nonlinear flow equations. This method also works for transport in the presence of an external magnetic field.


abstract: In the era of big data, reducing data dimensionality is critical in many areas of science. Widely used Principal Component Analysis (PCA) addresses this problem by computing a low dimensional data embedding that maximally explain variance of the data. However, PCA has two major weaknesses. Firstly, it only considers linear correlations among variables (features), and secondly it is not suitable for categorical data. We resolve these issues by proposing Maximally Correlated Principal Component Analysis (MCPCA). MCPCA computes transformations of variables whose covariance matrix has the largest Ky Fan norm. Variable transformations are unknown, can be nonlinear and are computed in an optimization. MCPCA can also be viewed as a multivariate extension of Maximal Correlation. For jointly Gaussian variables we show that the covariance matrix corresponding to the identity (or the negative of the identity) transformations majorizes covariance matrices of non-identity functions. Using this result we characterize global MCPCA optimizers for nonlinear functions of jointly Gaussian variables for every rank constraint. For categorical variables we characterize global MCPCA optimizers for the rank one constraint based on the leading eigenvector of a matrix computed using pairwise joint distributions. For a general rank constraint we propose a block coordinate descend algorithm and show its convergence to stationary points of the MCPCA optimization. We compare MCPCA with PCA and other state-of-the-art dimensionality reduction methods including Isomap, LLE, multilayer autoencoders (neural networks), kernel PCA, probabilistic PCA and diffusion maps on several synthetic and real datasets. We show that MCPCA consistently provides improved performance compared to other methods.


title: Threshold Constraints with Guarantees for Parity Objectives in Markov Decision Processes

abstract: The beyond worst-case synthesis problem was introduced recently by Bruy\`ere et al. [BFRR14]: it aims at building system controllers that provide strict worst-case performance guarantees against an antagonistic environment while ensuring higher expected performance against a stochastic model of the environment. Our work extends the framework of [BFRR14] and follow-up papers, which focused on quantitative objectives, by addressing the case of $\omega$-regular conditions encoded as parity objectives, a natural way to represent functional requirements of systems. We build strategies that satisfy a main parity objective on all plays, while ensuring a secondary one with sufficient probability. This setting raises new challenges in comparison to quantitative objectives, as one cannot easily mix different strategies without endangering the functional properties of the system. We establish that, for all variants of this problem, deciding the existence of a strategy lies in ${\sf NP} \cap {\sf coNP}$, the same complexity class as classical parity games. Hence, our framework provides additional modeling power while staying in the same complexity class. [BFRR14] V\'eronique Bruy\`ere, Emmanuel Filiot, Mickael Randour, and Jean-Fran\c{c}ois Raskin. Meet your expectations with guarantees: Beyond worst-case synthesis in quantitative games. In Ernst W. Mayr and Natacha Portier, editors, 31st International Symposium on Theoretical Aspects of Computer Science, STACS 2014, March 5-8, 2014, Lyon, France, volume 25 of LIPIcs, pages 199-213. Schloss Dagstuhl - Leibniz - Zentrum fuer Informatik, 2014.


title: Costas cubes

abstract: A Costas array is a permutation array for which the vectors joining pairs of $1$s are all distinct. We propose a new three-dimensional combinatorial object related to Costas arrays: an order $n$ Costas cube is an array $(d_{i,j,k})$ of size $n \times n \times n$ over $\mathbb{Z}_2$ for which each of the three projections of the array onto two dimensions, namely $(\sum_i d_{i,j,k})$ and $(\sum_j d_{i,j,k})$ and $(\sum_k d_{i,j,k})$, is an order $n$ Costas array. We determine all Costas cubes of order at most $29$, showing that Costas cubes exist for all these orders except $18$ and $19$ and that a significant proportion of the Costas arrays of certain orders occur as projections of Costas cubes. We then present constructions for four infinite families of Costas cubes.


title: A numerical study of the F-model with domain-wall boundaries

abstract: We perform a numerical study of the F-model with domain-wall boundary conditions. Various exact results are known for this particular case of the six-vertex model, including closed expressions for the partition function for any system size as well as its asymptotics and leading finite-size corrections. To complement this picture we use a full lattice multi-cluster algorithm to study equilibrium properties of this model for systems of moderate size, up to L=512. We compare the energy to its exactly known large-L asymptotics. We investigate the model's infinite-order phase transition by means of finite-size scaling for an observable derived from the staggered polarization in order to test the method put forward in our recent joint work with Duine and Barkema. In addition we analyse local properties of the model. Our data are perfectly consistent with analytical expressions for the arctic curves. We investigate the structure inside the temperate region of the lattice, confirming the oscillations in vertex densities that were first observed by Sylju{\aa}sen and Zvonarev, and recently studied by Lyberg et al. We point out '(anti)ferroelectric' oscillations close to the corresponding frozen regions as well as 'higher-order' oscillations forming an intricate pattern with saddle-point-like features.

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


title: Large Fluctuations of Radiation in Stochastically Activated Two-Level Systems

abstract: We study the large fluctuations of emitted radiations in the system of $N$ non-interacting two-level atoms. Two methods are used to calculate the probability of the large fluctuations and the time dependence of the excitation and emission. The first method is based on the large deviation principle for Markov processes. The second one uses an analogue of the quantum formalism for classical probability problems. Particularly we prove that in a large fluctuation limit approximately a half of the atoms are excited. This fact is independent on the fraction of the excited atoms in the stationary state.


title: Narrow band emission from Au-film on Cu- substrate

abstract: The experiments reveal the strong narrow band in the spectra of the emission from Au-film on the Cu substrate. The emission was excited by red He-Ne laser in standard geometry for Raman back-scattering experiments. Some additional experiments were made. The set of obtained data did not give enough arguments to make the final conclusion about the nature of this emission. The details of experiments, results, and possible explanations are discussed.


title: A new concept multi-stage Zeeman decelerator

abstract: We present a new concept for a multi-stage Zeeman decelerator that is optimized particularly for applications in molecular beam scattering experiments. The decelerator consists of a series of alternating hexapoles and solenoids, that effectively decouple the transverse focusing and longitudinal deceleration properties of the decelerator. It can be operated in a deceleration and acceleration mode, as well as in a hybrid mode that makes it possible to guide a particle beam through the decelerator at constant speed. The deceleration features phase stability, with a relatively large six-dimensional phase-space acceptance. The separated focusing and deceleration elements result in an unequal partitioning of this acceptance between the longitudinal and transverse directions. This is ideal in scattering experiments, which typically benefit from a large longitudinal acceptance combined with narrow transverse distributions. We demonstrate the successful experimental implementation of this concept using a Zeeman decelerator consisting of an array of 25 hexapoles and 24 solenoids. The performance of the decelerator in acceleration, deceleration and guiding modes is characterized using beams of metastable Helium ($^3S$) atoms. Up to 60 % of the kinetic energy was removed for He atoms that have an initial velocity of 520 m/s. The hexapoles consist of permanent magnets, whereas the solenoids are produced from a single hollow copper capillary through which cooling liquid is passed. The solenoid design allows for excellent thermal properties, and enables the use of readily available and cheap electronics components to pulse high currents through the solenoids. The Zeeman decelerator demonstrated here is mechanically easy to build, can be operated with cost-effective electronics, and can run at repetition rates up to 10 Hz


title: Experimental realization of Floquet PT-symmetric systems

abstract: We provide an experimental framework where periodically driven PT-symmetric systems can be investigated. The set-up, consisting of two UHF oscillators coupled by a time-dependent capacitance, demonstrates a cascade of PT-symmetric broken domains bounded by exceptional point degeneracies. These domains are analyzed and understood using an equivalent Floquet frequency lattice with local PT-symmetry. Management of these PT-phase transition domains is achieved through the amplitude and frequency of the drive.


title: Violation of space-time Bell-CHSH inequality beyond Tsirelson bound via post selection and quantum cryptography

abstract: Here we show that, if we insert context dependent unitary evolutions (which can be achieved via post selection) into spatial (i.e., normal) Bell-CHSH test, then it is possible to violate space-time Bell-CHSH inequality maximally (i.e., up to $4$). However this does not contradict Tsirelson quantum bound ($2\sqrt{2}$), as the latter has been derived without taking into consideration context dependent unitary evolutions and/or post selection. As an important application, this leads to a more efficient quantum key distribution protocol.


title: Computing toric degenerations of flag varieties

abstract: We compute toric degenerations arising from the tropicalization of the full flag varieties $\mathcal{F}\ell_4$ and $\mathcal{F}\ell_5$ embedded in a product of Grassmannians. For $\mathcal{F}\ell_4$ and $\mathcal{F}\ell_5$ we compare toric degenerations arising from string polytopes and the FFLV polytope with those obtained from the tropicalization of the flag varieties. We also present a general procedure to find toric degenerations in the cases where the initial ideal arising from a cone of the tropicalization of a variety is not prime.


title: Understanding the importance of transient resonances in extreme mass ratio inspirals

abstract: Extreme mass ratio inspirals (EMRIs) occur when a compact object orbits a much larger one, like a solar-mass black hole around a supermassive black hole. The orbit has 3 frequencies which evolve through the inspiral. If the orbital radial frequency and polar frequency become commensurate, the system passes through a transient resonance. Evolving through resonance causes a jump in the evolution of the orbital parameters. We study these jumps and their impact on EMRI gravitational-wave detection. Jumps are smaller for lower eccentricity orbits; since most EMRIs have small eccentricities when passing through resonances, we expect that the impact on detection will be small. Neglecting the effects of transient resonances leads to a loss of ~4% of detectable signals for an astrophysically motivated population of EMRIs.


title: Dark Catalysis

abstract: Recently it was shown that dark matter with mass of order the weak scale can be charged under a new long-range force, decoupled from the Standard Model, with only weak constraints from early Universe cosmology. Here we consider the implications of an additional charged particle $C$ that is light enough to lead to significant dissipative dynamics on galactic times scales. We highlight several novel features of this model, which can be relevant even when the $C$ particle constitutes only a small fraction of the number density (and energy density). We assume a small asymmetric abundance of the $C$ particle whose charge is compensated by a heavy $X$ particle so that the relic abundance of dark matter consists mostly of symmetric $X$ and $\bar{X}$, with a small asymmetric component made up of $X$ and $C$. As the universe cools, it undergoes asymmetric recombination binding the free $C$s into $(XC)$ dark atoms efficiently. Even with a tiny asymmetric component, the presence of $C$ particles catalyzes tight coupling between the heavy dark matter $X$ and the dark photon plasma that can lead to a significant suppression of the matter power spectrum on small scales and lead to some of the strongest bounds on such dark matter theories. We find a viable parameter space where structure formation constraints are satisfied and significant dissipative dynamics can occur in galactic haloes but show a large region is excluded. Our model shows that subdominant components in the dark sector can dramatically affect structure formation.


title: Once in a blue moon: detection of 'bluing' during debris transits in the white dwarf WD1145+017

abstract: The first transiting planetesimal orbiting a white dwarf was recently detected in K2 data of WD1145+017 and has been followed up intensively. The multiple, long, and variable transits suggest the transiting objects are dust clouds, probably produced by a disintegrating asteroid. In addition, the system contains circumstellar gas, evident by broad absorption lines, mostly in the u'-band, and a dust disc, indicated by an infrared excess. Here we present the first detection of a change in colour of WD1145+017 during transits, using simultaneous multi-band fast-photometry ULTRACAM measurements over the u'g'r'i'-bands. The observations reveal what appears to be 'bluing' during transits; transits are deeper in the redder bands, with a u'-r' colour difference of up to ~-0.05 mag. We explore various possible explanations for the bluing. 'Spectral' photometry obtained by integrating over bandpasses in the spectroscopic data in- and out-of-transit, compared to the photometric data, shows that the observed colour difference is most likely the result of reduced circumstellar absorption in the spectrum during transits. This indicates that the transiting objects and the gas share the same line-of-sight, and that the gas covers the white dwarf only partially, as would be expected if the gas, the transiting debris, and the dust emitting the infrared excess, are part of the same general disc structure (although possibly at different radii). In addition, we present the results of a week-long monitoring campaign of the system.


title: Deciphering the MSSM Higgs Mass at Future Hadron Colliders

abstract: Future hadron colliders will have a remarkable capacity to discover massive new particles, but their capabilities for precision measurements of couplings that can reveal underlying mechanisms have received less study. In this work we study the capability of future hadron colliders to shed light on a precise, focused question: is the higgs mass of 125 GeV explained by the MSSM? If supersymmetry is realized near the TeV scale, a future hadron collider could produce huge numbers of gluinos and electroweakinos. We explore whether precision measurements of their properties could allow inference of the scalar masses and $\tan \beta$ with sufficient accuracy to test whether physics beyond the MSSM is needed to explain the higgs mass. We also discuss dark matter direct detection and precision higgs physics as complementary probes of $\tan \beta$. For concreteness, we focus on the mini-split regime of MSSM parameter space at a 100 TeV $pp$ collider, with scalar masses ranging from 10s to about 1000 TeV.


title: KINETyS: Constraining spatial variations of the stellar initial mass function in early-type galaxies

abstract: The heavyweight stellar initial mass function (IMF) observed in the cores of massive early-type galaxies (ETGs) has been linked to formation of their cores in an initial swiftly-quenched rapid starburst. However, the outskirts of ETGs are thought to be assembled via the slow accumulation of smaller systems in which the star formation is less extreme; this suggests the form of the IMF should exhibit a radial trend in ETGs. Here we report radial stellar population gradients out to the half-light radii of a sample of eight nearby ETGs. Spatially resolved spectroscopy at 0.8-1.35{\mu}m from the VLT's KMOS instrument was used to measure radial trends in the strengths of a variety of IMF-sensitive absorption features (including some which are previously unexplored). We find weak or no radial variation in some of these which, given a radial IMF trend, ought to vary measurably, e.g. for the Wing-Ford band we measure a gradient of +0.06$\pm$0.04 per decade in radius. Using stellar population models to fit stacked and individual spectra, we infer that the measured radial changes in absorption feature strengths are primarily accounted for by abundance gradients which are fairly consistent across our sample (e.g. we derive an average [Na/H] gradient of -0.53$\pm$0.07). The inferred contribution of dwarf stars to the total light typically corresponds to a bottom heavy IMF, but we find no evidence for radial IMF variations in the majority of our sample galaxies.


title: Non-toric Cones and Chern-Simons Quivers

abstract: We obtain an integral formula for the volume of non-toric tri-Sasaki Einstein manifolds arising from nonabelian hyperkahler quotients. The derivation is based on equivariant localization and generalizes existing formulas for Abelian quotients, which lead to toric manifolds. The formula is particularly valuable in the context of AdS$_{4}\times Y_{7}$ vacua of M-theory and their field theory duals. As an application, we consider 3d $\mathcal N=3$ Chern-Simons theories with affine ADE quivers. While the $\widehat A$ series corresponds to toric $Y_{7}$, the $\widehat D$ and $\widehat E$ series are non-toric. We compute the volumes of the corresponding seven-manifolds and compare to the prediction from supersymmetric localization in field theory, finding perfect agreement. This is the first test of an infinite number of non-toric AdS$_4$/CFT$_3$ dualities.


title: Scalable explicit implementation of anisotropic diffusion with Runge-Kutta-Legendre super-time-stepping

abstract: An important ingredient in numerical modelling of high temperature magnetised astrophysical plasmas is the anisotropic transport of heat along magnetic field lines from higher to lower temperatures.Magnetohydrodynamics (MHD) typically involves solving the hyperbolic set of conservation equations along with the induction equation. Incorporating anisotropic thermal conduction requires to also treat parabolic terms arising from the diffusion operator. An explicit treatment of parabolic terms will considerably reduce the simulation time step due to its dependence on the square of the grid resolution ($\Delta x$) for stability. Although an implicit scheme relaxes the constraint on stability, it is difficult to distribute efficiently on a parallel architecture. Treating parabolic terms with accelerated super-time stepping (STS) methods has been discussed in literature but these methods suffer from poor accuracy (first order in time) and also have difficult-to-choose tuneable stability parameters. In this work we highlight a second order (in time) Runge Kutta Legendre (RKL) scheme (first described by Meyer et. al. 2012) that is robust, fast and accurate in treating parabolic terms alongside the hyperbolic conversation laws. We demonstrate its superiority over the first order super time stepping schemes with standard tests and astrophysical applications. We also show that explicit conduction is particularly robust in handling saturated thermal conduction. Parallel scaling of explicit conduction using RKL scheme is demonstrated up to more than $10^4$ processors.


title: An update to the EVEREST K2 pipeline: Short cadence, saturated stars, and Kepler-like photometry down to Kp = 15

abstract: We present an update to the EVEREST K2 pipeline that addresses various limitations in the previous version and improves the photometric precision of the de-trended light curves. We develop a fast regularization scheme for third order pixel level decorrelation (PLD) and adapt the algorithm to include the PLD vectors of neighboring stars to enhance the predictive power of the model and minimize overfitting, particularly for faint stars. We also modify PLD to work for saturated stars and improve its performance on extremely variable stars. On average, EVEREST 2.0 light curves have 10-20% higher photometric precision than those in the previous version, yielding the highest precision light curves at all Kp magnitudes of any publicly available K2 catalog. For most K2 campaigns, we recover the original Kepler precision to at least Kp = 14, and to at least Kp = 15 for campaigns 1, 5, and 6. We also de-trend all short cadence targets observed by K2, obtaining even higher photometric precision for these stars. All light curves for campaigns 0-8 are available online in the EVEREST catalog, which will be continuously updated with future campaigns. EVEREST 2.0 is open source and is coded in a general framework that can be applied to other photometric surveys, including Kepler and the upcoming TESS mission.


title: The Odd story of $\alpha'$-corrections

abstract: The $\alpha'$-deformed frame-like Double Field Theory (DFT) is a T-duality and gauge invariant extension of DFT in which generalized Green-Schwarz transformations provide a gauge principle that fixes the higher-derivative corrections. It includes all the first order $\alpha'$-corrections of the bosonic and heterotic string low energy effective actions and of the Hohm-Siegel-Zwiebach $\alpha'$-geometry. Here we gauge this theory and parameterize it in terms of a frame, a two-form, a dilaton, gauge vectors and scalar fields. This leads to a unified framework that extends the previous construction by including all duality constrained interactions in generic (gauged/super)gravity effective field theories in arbitrary number of dimensions, to first order in $\alpha'$.


title: Violating weak cosmic censorship in AdS$_4$

abstract: We consider time-dependent solutions of the Einstein-Maxwell equations using anti-de Sitter (AdS) boundary conditions, and provide the first counterexample to the weak cosmic censorship conjecture in four spacetime dimensions. Our counterexample is entirely formulated in the Poincar\'e patch of AdS. We claim that our results have important consequences for quantum gravity, most notably to the weak gravity conjecture.


title: Identification of partially resolved binaries in Pan-STARRS1 data

abstract: Using shape measurement techniques developed for weak lensing surveys we have identified three new ultracool binaries in the Pan-STARRS1 survey. Binary companions which are not completely resolved can still alter the shapes of stellar images. These shape distortions can be measured if PSF anisotropy caused by the telescope is properly accounted for. We show using both a sample of known binary stars and simulated binaries that we can reliably recover binaries wider than around 0.3" and with flux ratios greater than around 0.1. We then applied our method to a sample of ultracool dwarfs within 30pc with 293 objects having sufficient Pan-STARRS1 data for our method. In total we recovered all but one of the 11 binaries wider than 0.3" in this sample. Our one failure was a true binary detected with a significant but erroneously high ellipticity which led it to be rejected in our analysis. We identify three new binaries, one a simultaneous discovery, with primary spectral types M6.5, L1 and T0.5. These latter two were confirmed with Keck/NIRC2 follow-up imaging. This technique will be useful for identifying large numbers of stellar and substellar binaries in the upcoming LSST and DES sky surveys.


title: Quantum Simulation of the Abelian-Higgs Lattice Gauge Theory with Ultracold Atoms

abstract: We present a quantum simulation scheme for the Abelian-Higgs lattice gauge theory using ultracold bosonic atoms in optical lattices. The model contains both gauge and Higgs scalar fields, and exhibits interesting phases related to confinement and the Higgs mechanism. The model can be simulated by an atomic Hamiltonian, by first mapping the local gauge symmetry to an internal symmetry of the atomic system, the conservation of hyperfine angular momentum in atomic collisions. By including auxiliary bosons in the simulation, we show how the Abelian-Higgs Hamiltonian emerges effectively. We analyze the accuracy of our method in terms of different experimental parameters, as well as the effect of the finite number of bosons on the quantum simulator. Finally, we propose possible experiments for studying the ground state of the system in different regimes of the theory, and measuring interesting high energy physics phenomena in real time.


title: Bimodal radio variability in OVRO-40m-monitored blazars

abstract: Blazars are known to show periods of quiescence followed by outbursts visible throughout the electromagnetic spectrum. We present a novel maximum likelihood approach to capture this bimodal behavior by examining blazar radio variability in the flux-density domain. We separate quiescent and flaring components of a source's light curve by modeling its flux-density distribution as a series of "off" and "on" states. Our modeling allows us to extract information regarding the flaring ratio, duty cycle, and the modulation index in the "off"-state, in the "on"-state, as well as throughout the monitoring period of each blazar. We apply our method to a flux-density-limited subsample from the Owens Valley Radio observatory's 15 GHz blazar monitoring program, and explore differences in the variability characteristics between BL Lacs and FSRQs as well as between $\gamma$-ray detected and non-detected sources. We find that: (1) BL Lacs are more variable and have relatively larger outbursts than the FSRQs, (2) unclassified blazar candidates in our sample show similar variability characteristics as the FSRQs, and (3) $\gamma$-ray detected differ from the $\gamma$-ray non-detected sources in all their variability properties, suggesting a link between the production of $\gamma$-rays and the mechanism responsible for the radio variability. Finally, we fit distributions for blazar flaring ratios, duty cycles, and on- and off- modulation indices that can be used in population studies of variability-dependent blazar properties.


title: Time-resolved polarimetry of the superluminous SN 2015bn with the Nordic Optical Telescope

abstract: We present imaging polarimetry of the superluminous supernova SN 2015bn, obtained over nine epochs between $-$20 and $+$46 days with the Nordic Optical Telescope. This was a nearby, slowly-evolving Type I superluminous supernova that has been studied extensively and for which two epochs of spectropolarimetry are also available. Based on field stars, we determine the interstellar polarisation in the Galaxy to be negligible. The polarisation of SN 2015bn shows a statistically significant increase during the last epochs, confirming previous findings. Our well-sampled imaging polarimetry series allows us to determine that this increase (from $\sim 0.54\%$ to $\gtrsim 1.10\%$) coincides in time with rapid changes that took place in the optical spectrum. We conclude that the supernova underwent a `phase transition' at around $+$20 days, when the photospheric emission shifted from an outer layer, dominated by natal C and O, to a more aspherical inner core, dominated by freshly nucleosynthesized material. This two-layered model might account for the characteristic appearance and properties of Type I superluminous supernovae.


title: Darboux theory of integrability for real polynomial vector fields on $\sss^n$

abstract: This is a survey on the Darboux theory of integrability for polynomial vector fields, first in $\R^n$ and second in the $n$-dimensional sphere $\sss^n$. We also provide new results about the maximum number of parallels and meridians that a polynomial vector field $\X$ on $\sss^n$ can have in function of its degree. These results in some sense extend the known result on the maximum number of hyperplanes that a polynomial vector field $\Y$ in $\R^n$ can have in function of the degree of $\Y$.


title: Behavioral Analytics for Myopic Agents

abstract: Many multi-agent systems have the structure of a single coordinator providing behavioral or financial incentives to a large number of agents. Two challenges faced by the coordinator are a finite budget from which to allocate incentives, and an initial lack of knowledge about the utility function of the agents. Here, we present a behavioral analytics approach to solve the coordinator's problem when agents make decisions by maximizing utility functions that depend on prior system states, inputs, and other parameters that are initially unknown and subject to temporal dynamics. Our behavioral analytics framework involves three steps: first, we develop a behavioral model that describes the decision-making process of an agent; second, we use data to estimate behavioral model parameters for each agent and then use these estimates to predict future decisions of each agent; and third, we use the estimated behavioral model parameters to optimize a set of costly incentives to provide to each agent. In this paper, we describe a specific set of tools, models, and approaches that fit into this framework, and that adapt models and incentives as new information is collected by repeating the second and third steps of this framework. Furthermore, we prove that incentives computed by this adaptive approach are asymptotically optimal with respect to a given loss function that describes the coordinator's objective. We optimize incentives using a decomposition scheme, where each sub-problem solves the coordinator's problem for a single agent, and the master problem is a pure integer program. We conclude with a simulation study to evaluate the effectiveness of our behavioral analytics approach in designing personalized treatment plans for a weight loss program. The results show our approach maintains efficacy of the program while reducing costs by up to 60%, while adaptive heuristics provide less savings.


title: Multiple D3-instantons and mock modular forms II

abstract: We analyze the modular properties of D3-brane instanton corrections to the hypermultiplet moduli space in type IIB string theory compactified on a Calabi-Yau threefold. In Part I, we found a necessary condition for the existence of an isometric action of S-duality on this moduli space: the generating function of DT invariants in the large volume attractor chamber must be a vector-valued mock modular form with specified modular properties. In this work, we prove that this condition is also sufficient at two-instanton order. This is achieved by producing a holomorphic action of SL(2,Z) on the twistor space which preserves the holomorphic contact structure. The key step is to cancel the anomalous modular variation of the Darboux coordinates by a local holomorphic contact transformation, which is generated by a suitable indefinite theta series. For this purpose we introduce a new family of theta series of signature (2,n-2), find their modular completion, and conjecture sufficient conditions for their convergence, which may be of independent mathematical interest.


title: $B_s\to\mu^+\mu^-$ as current and future probe of new physics

abstract: The rare flavour-changing neutral current decay $B_s\to\mu^+\mu^-$ is among the most important indirect probes of new physics at the LHC, as it is strongly suppressed in the Standard Model, very sensitive to new physics effects, and theoretically exceptionally clean. We present a thorough state-of-the-art analysis of the constraints on new physics from present and future measurements of this decay, focusing on scalar operators. We show model-independently and in concrete new physics models, namely the MSSM and two leptoquark scenarios, that a future precise measurement of the mass-eigenstate rate asymmetry in $B_s\to\mu^+\mu^-$ would allow to disentangle new physics scenarios that would be indistinguishable based on measurements of the branching ratio alone. We also highlight the complementarity between $B_s\to\mu^+\mu^-$ and direct searches in both model classes. Our numerics is based on the open source code flavio and is made publicly available.


title: When is a Network a Network? Multi-Order Graphical Model Selection in Pathways and Temporal Networks

abstract: We introduce a framework for the modeling of sequential data capturing pathways of varying lengths observed in a network. Such data are important, e.g., when studying click streams in information networks, travel patterns in transportation systems, information cascades in social networks, biological pathways or time-stamped social interactions. While it is common to apply graph analytics and network analysis to such data, recent works have shown that temporal correlations can invalidate the results of such methods. This raises a fundamental question: when is a network abstraction of sequential data justified? Addressing this open question, we propose a framework which combines Markov chains of multiple, higher orders into a multi-layer graphical model that captures temporal correlations in pathways at multiple length scales simultaneously. We develop a model selection technique to infer the optimal number of layers of such a model and show that it outperforms previously used Markov order detection techniques. An application to eight real-world data sets on pathways and temporal networks shows that it allows to infer graphical models which capture both topological and temporal characteristics of such data. Our work highlights fallacies of network abstractions and provides a principled answer to the open question when they are justified. Generalizing network representations to multi-order graphical models, it opens perspectives for new data mining and knowledge discovery algorithms.

categories: cs.SI cond-mat.dis-nn physics.data-an physics.soc-ph


title: Natural and Artificial Spectral Edges in Exoplanets

abstract: Technological civilizations may rely upon large-scale photovoltaic arrays to harness energy from their host star. Photovoltaic materials, such as silicon, possess distinctive spectral features, including an "artificial edge" that is characteristically shifted in wavelength shortwards of the "red edge" of vegetation. Future observations of reflected light from exoplanets would be able to detect both natural and artificial edges photometrically, if a significant fraction of the planet's surface is covered by vegetation or photovoltaic arrays respectively. The stellar energy thus tapped can be utilized for terraforming activities by transferring heat and light from the day side to the night side on tidally locked exoplanets, thereby producing detectable artifacts.


title: Limits on the heralding efficiencies and spectral purities of spectrally-filtered single photons from photon pair sources

abstract: Photon pairs produced by parametric down-conversion or four-wave mixing can interfere with each other in multiport interferometers, or carry entanglement between distant nodes for use in entanglement swapping. This requires the photons be spectrally pure to ensure good interference, and have high heralding efficiency to know accurately the number of photons involved and to maintain high rates as the number of photons grows. Spectral filtering is often used to remove noise and define spectral properties. For heralded single photons high purity and heralding efficiency is possible by filtering the heralding arm, but when both photons in typical pair sources are filtered, we show that the heralding efficiency of one or both of the photons is strongly reduced even by ideal spectral filters with 100% transmission in the passband: any improvement in reduced-state spectral purity from filtering comes at the cost of lowered heralding efficiency. We consider the fidelity to a pure, lossless single photon, symmetrize it to include both photons of the pair, and show this quantity is intrinsically limited for sources with spectral correlation. We then provide a framework for this effect for benchmarking common photon pair sources, and present an experiment where we vary the photon filter bandwidths and measure the increase in purity and corresponding reduction in heralding efficiency.


title: Inclusive Production Through AdS/CFT

abstract: It has been shown that AdS/CFT calculations can reproduce certain exclusive 2->2 cross sections in QCD at high energy, both for near-forward and for fixed-angle scattering. In this paper, we extend prior treatments by using AdS/CFT to calculate the inclusive single-particle production cross section in QCD at high center-of-mass energy. We find that conformal invariance in the UV restricts the cross section to have a characteristic power-law falloff in the transverse momentum of the produced particle, with the exponent given by twice the conformal dimension of the produced particle, independent of incoming particle types. We conclude by comparing our findings to recent LHC experimental data from ATLAS and ALICE, and find good agreement.


title: Elliptic theory for sets with higher co-dimensional boundaries

abstract: Many geometric and analytic properties of sets hinge on the properties of harmonic measure, notoriously missing for sets of higher co-dimension. The aim of this manuscript is to develop a version of elliptic theory, associated to a linear PDE, which ultimately yields a notion analogous to that of the harmonic measure, for sets of codimension higher than 1. To this end, we turn to degenerate elliptic equations. Let $\Gamma \subset \mathbb R^n$ be an Ahlfors regular set of dimension $d<n-1$ (not necessarily integer) and $\Omega = \mathbb R^n \setminus \Gamma$. Let $L = - {\rm div} A\nabla$ be a degenerate elliptic operator with measurable coefficients such that the ellipticity constants of the matrix $A$ are bounded from above and below by a multiple of ${\rm dist}(\cdot, \Gamma)^{d+1-n}$. We define weak solutions; prove trace and extension theorems in suitable weighted Sobolev spaces; establish the maximum principle, De Giorgi-Nash-Moser estimates, the Harnack inequality, the H\"older continuity of solutions (inside and at the boundary). We define the Green function and provide the basic set of pointwise and/or $L^p$ estimates for the Green function and for its gradient. With this at hand, we define harmonic measure associated to $L$, establish its doubling property, non-degeneracy, change-of-the-pole formulas, and, finally, the comparison principle for local solutions. In another article to appear, we will prove that when $\Gamma$ is the graph of a Lipschitz function with small Lipschitz constant, we can find an elliptic operator $L$ for which the harmonic measure given here is absolutely continuous with respect to the $d$-Hausdorff measure on $\Gamma$ and vice versa. It thus extends Dahlberg's theorem to some sets of codimension higher than 1.


title: QCD Axion Star Collapse with the Chiral Potential

abstract: In a previous work, we analyzed collapsing axion stars using the low-energy instanton potential, showing that the total energy is always bounded and that collapsing axion stars do not form black holes. In this paper, we provide a proof that the conclusions are unchanged when using instead the more general chiral potential for QCD axions.


title: Equisingular and Equinormalizable Deformations of Isolated Non-Normal Singularities

abstract: We present new results on equisingularity and equinormalizability of families with isolated non-normal singularities (INNS) of arbitrary dimension. We define a $\delta$-invariant and a $\mu$-invariant for an INNS and prove necessary and sufficient numerical conditions for equinormalizability and weak equinormalizability using $\delta$ and $\mu$. Moreover, we determine the number of connected components of the Milnor fibre of an arbitrary INNS. For families of generically reduced curves, we investigate the topological behavior of the Milnor fibre and characterize topological triviality of such families. Finally we state some open problems and conjectures. In addition we give a survey of classical results about equisingularity and equinormalizability so that the article may be useful as a reference source.


title: An Unsupervised Approach for Overlapping Cervical Cell Cytoplasm Segmentation

abstract: The poor contrast and the overlapping of cervical cell cytoplasm are the major issues in the accurate segmentation of cervical cell cytoplasm. This paper presents an automated unsupervised cytoplasm segmentation approach which can effectively find the cytoplasm boundaries in overlapping cells. The proposed approach first segments the cell clumps from the cervical smear image and detects the nuclei in each cell clump. A modified Otsu method with prior class probability is proposed for accurate segmentation of nuclei from the cell clumps. Using distance regularized level set evolution, the contour around each nucleus is evolved until it reaches the cytoplasm boundaries. Promising results were obtained by experimenting on ISBI 2015 challenge dataset.


title: Magnetars in Ultra-Long Gamma-Ray Bursts

abstract: Supernova 2011kl, associated with the ultra-long gamma-ray burst (ULGRB) 111209A, exhibited a higher-than-normal peak luminosity, placing it in the parameter space between regular supernovae and super-luminous supernovae. Its light curve can only be matched by an abnormally high fraction of $^{56}$Ni that appears inconsistent with the observed spectrum, and as a result it has been suggested that the supernova, and by extension the gamma-ray burst, are powered by the spin-down of a highly magnetised millisecond pulsar, known as a magnetar. We investigate the broadband observations of ULGRB 111209A, and find two independent measures that suggest a high density circumburst environment. However, the light curve of the GRB afterglow shows no evidence of a jet break (the steep decline that would be expected as the jet slows due to the resistance of the external medium) out to three weeks after trigger, implying a wide jet. Combined with the high isotropic energy of the burst, this implies that only a magnetar with a spin period of ~ 1 ms or faster can provide enough energy to power both ULGRB 111209A and Supernova 2011kl.


title: Ultraluminous X-ray sources

abstract: The origin of Ultraluminous X-ray Sources (ULXs) in external galaxies whose X-ray luminosities exceed those of the brightest black holes in our Galaxy by hundreds and thousands of times is mysterious. The most popular models for the ULXs involve either intermediate mass black holes (IMBHs) or stellar-mass black holes accreting at super-Eddington rates. Here we review the ULX properties. Their X-ray spectra indicate a presence of hot winds in their accretion disks supposing the supercritical accretion. In recent years, new surprising results were discovered in X-ray data, ULX-pulsars and high-velocity outflows up to 0.2 c. They are also in accordance with the super-Eddington accretion. However, the strongest evidences come from optical spectroscopy. The spectra of the ULX counterparts are very similar to that of SS433, the only known supercritical accretor in our Galaxy. The spectra are apparently of WNL type (late nitrogen Wolf-Rayet stars) or LBV (luminous blue variables) in their hot state, which are very scarce stellar objects. We find that the spectra do not originate from WNL/LBV type donors but from very hot winds from the accretion disks, which have similar physical conditions as the stellar winds from these stars. The results suggest that bona-fide ULXs must constitute a homogeneous class of objects, which most likely have supercritical accretion disks.


title: New lower bounds for hypergraph Ramsey numbers

abstract: The Ramsey number $r_k(s,n)$ is the minimum $N$ such that for every red-blue coloring of the $k$-tuples of $\{1,\ldots, N\}$, there are $s$ integers such that every $k$-tuple among them is red, or $n$ integers such that every $k$-tuple among them is blue. We prove the following new lower bounds for 4-uniform hypergraph Ramsey numbers: $$r_4(5,n) > 2^{n^{c\log n}} \qquad \hbox{ and } \qquad r_4(6,n) > 2^{2^{cn^{1/5}}},$$ where $c$ is an absolute positive constant. This substantially improves the previous best bounds of $2^{n^{c\log\log n}}$ and $2^{n^{c\log n}}$, respectively. Using previously known upper bounds, our result implies that the growth rate of $r_4(6,n)$ is double exponential in a power of $n$. As a consequence, we obtain similar bounds for the $k$-uniform Ramsey numbers $r_k(k+1, n)$ and $r_k(k+2, n)$ where the exponent is replaced by an appropriate tower function. This almost solves the question of determining the tower growth rate for {\emph {all}} classical off-diagonal hypergraph Ramsey numbers, a question first posed by Erd\H os and Hajnal in 1972. The only problem that remains is to prove that $r_4(5,n)$ is double exponential in a power of $n$.


title: Java Code Analysis and Transformation into AWS Lambda Functions

abstract: Software developers are faced with the issue of either adapting their programming model to the execution model (e.g. cloud platforms) or finding appropriate tools to adapt the model and code automatically. A recent execution model which would benefit from automated enablement is Function-as-a-Service. Automating this process requires a pipeline which includes steps for code analysis, transformation and deployment. In this paper, we outline the design and runtime characteristics of Podilizer, a tool which implements the pipeline specifically for Java source code as input and AWS Lambda as output. We contribute technical and economic metrics about this concrete 'FaaSification' process by observing the behaviour of Podilizer with two representative Java software projects.


title: A Concurrent Perspective on Smart Contracts

abstract: In this paper, we explore remarkable similarities between multi-transactional behaviors of smart contracts in cryptocurrencies such as Ethereum and classical problems of shared-memory concurrency. We examine two real-world examples from the Ethereum blockchain and analyzing how they are vulnerable to bugs that are closely reminiscent to those that often occur in traditional concurrent programs. We then elaborate on the relation between observable contract behaviors and well-studied concurrency topics, such as atomicity, interference, synchronization, and resource ownership. The described contracts-as-concurrent-objects analogy provides deeper understanding of potential threats for smart contracts, indicate better engineering practices, and enable applications of existing state-of-the-art formal verification techniques.


title: soc2seq: Social Embedding meets Conversation Model

abstract: While liking or upvoting a post on a mobile app is easy to do, replying with a written note is much more difficult, due to both the cognitive load of coming up with a meaningful response as well as the mechanics of entering the text. Here we present a novel textual reply generation model that goes beyond the current auto-reply and predictive text entry models by taking into account the content preferences of the user, the idiosyncrasies of their conversational style, and even the structure of their social graph. Specifically, we have developed two types of models for personalized user interactions: a content-based conversation model, which makes use of location together with user information, and a social-graph-based conversation model, which combines content-based conversation models with social graphs.


title: Rethinking High Performance Computing Platforms: Challenges, Opportunities and Recommendations

abstract: A new class of Second generation high-performance computing applications with heterogeneous, dynamic and data-intensive properties have an extended set of requirements, which cover application deployment, resource allocation, -control, and I/O scheduling. These requirements are not met by the current production HPC platform models and policies. This results in a loss of opportunity, productivity and innovation for new computational methods and tools. It also decreases effective system utilization for platform providers due to unsupervised workarounds and rogue resource management strategies implemented in application space. In this paper we critically discuss the dominant HPC platform model and describe the challenges it creates for second generation applications because of its asymmetric resource view, interfaces and software deployment policies. We present an extended, more symmetric and application-centric platform model that adds decentralized deployment, introspection, bidirectional control and information flow and more comprehensive resource scheduling. We describe cHPC: an early prototype of a non-disruptive implementation based on Linux Containers (LXC). It can operate alongside existing batch queuing systems and exposes a symmetric platform API without interfering with existing applications and usage modes. We see our approach as a viable, incremental next step in HPC platform evolution that benefits applications and platform providers alike. To demonstrate this further, we layout out a roadmap for future research and experimental evaluation.


title: Quantum dimer models emerging from large-spin ultracold atoms

abstract: We propose an experimental protocol for using cold atoms to create and probe quantum dimer models, thereby exploring the Pauling-Anderson vision of a macroscopic collection of resonating bonds. This process can allow the study of exotic crystalline phases, fractionalization, topological spin liquids, and the relationship between resonating dimers and superconductivity subjects which have been challenging to address in solid-state experiments. Our key technical development is considering the action of an off-resonant photoassociation laser on large spin atoms localized at the sites of a deep optical lattice. The resulting superexchange interaction favors nearest-neighbor singlets. We derive an effective Hamiltonian in terms of these dimer degrees of freedom, finding that it is similar to well-known quantum dimer models, which boast a rich variety of valence bond crystal and spin liquid phases. We numerically study the ground state, explain how to tune the parameters, and develop a protocol to directly measure the dimers and their resonating patterns.


title: Overview: Generalizations of Multi-Agent Path Finding to Real-World Scenarios

abstract: Multi-agent path finding (MAPF) is well-studied in artificial intelligence, robotics, theoretical computer science and operations research. We discuss issues that arise when generalizing MAPF methods to real-world scenarios and four research directions that address them. We emphasize the importance of addressing these issues as opposed to developing faster methods for the standard formulation of the MAPF problem.


title: Quantum Identities for the Action

abstract: In this paper we derive various identities involving the {\it action} functional which enters the path-integral formulation of quantum mechanics. They provide some kind of generalisations of the Ehrenfest theorem giving correlations between powers of the action and its functional derivatives.


title: The evolution of gravitons in accelerating cosmologies: the case of extended gravity

abstract: We discuss the production and evolution of cosmological gravitons showing how the cosmological background affects their dynamics. Besides, the detection of cosmological gravitons could constitute an extremely important signature to discriminate among different cosmological models. Here we consider the cases of scalar-tensor gravity and $f(R)$ gravity where it is demonstrated the amplification of graviton amplitude changes if compared with General Relativity. Possible observational constraints are discussed.


title: Sampling strategies for fast updating of Gaussian Markov random fields

abstract: Gaussian Markov random fields (GMRFs) are popular for modeling dependence in large areal datasets due to their ease of interpretation and computational convenience afforded by the sparse precision matrices needed for random variable generation. Typically in Bayesian computation, GMRFs are updated jointly in a block Gibbs sampler or componentwise in a single-site sampler via the full conditional distributions. The former approach can speed convergence by updating correlated variables all at once, while the latter avoids solving large matrices. We consider a sampling approach in which the underlying graph can be cut so that conditionally independent sites are updated simultaneously. This algorithm allows a practitioner to parallelize updates of subsets of locations or to take advantage of `vectorized' calculations in a high-level language such as R. Through both simulated and real data, we demonstrate computational savings that can be achieved versus both single-site and block updating, regardless of whether the data are on a regular or an irregular lattice. The approach provides a good compromise between statistical and computational efficiency and is accessible to statisticians without expertise in numerical analysis or advanced computing.


title: The peculiar light-curve of J1415+1320: A case study in extreme scattering events

abstract: The radio light-curve of J1415+1320 (PKS 1413+135) shows time-symmetric and recurring U-shaped features across the cm-wave and mm-wave bands. The symmetry of these features points to lensing by an intervening object as the cause. U-shaped events in radio light curves in the cm-wave band have previously been attributed to Extreme scattering events (ESE). ESEs are thought to be the result of lensing by compact plasma structures in the Galactic interstellar medium, but the precise nature of these plasma structures remains unknown. Since the strength of a plasma lens evolves with wavelength $\lambda$ as $\lambda^2$, the presence of correlated variations at over a wide wavelength range casts doubt on the canonical ESE interpretation for J1415+1320. In this paper, we critically examine the evidence for plasma lensing in J1415+1320. We compute limits on the lensing strength, and the associated free-free opacity of the putative plasma lenses. We compare the observed and model ESE light curves, and also derive a lower limit on the lens distance based on the effects of parallax due to the Earth's orbit around the Sun. We conclude that plasma lensing is not a viable interpretation for J1415+1320's light curves and that symmetric U-shaped features in the radio light curves of extragalactic sources do not present {\em prima facie} evidence for ESEs. The methodology presented here is generic enough to be applicable to any plasma lensing candidate.


title: The role of disorder in NaO$_2$ and its implications for Na-O$_2$ batteries

abstract: We present a DFT study utilizing the Hubbard U correction to probe structural and magnetic disorder in $\mathrm{NaO_{2}}$, primary discharge product of Na-O$_2$ batteries. We show that $\mathrm{NaO_{2}}$ exhibits a large degree of rotational and magnetic disorder; a 3-body Ising Model is necessary to capture the subtle interplay of this disorder. MC simulations demonstrate that energetically favorable, FM phases near room temperature consist of alternating bands of orthogonally-oriented $\mathrm{O_{2}}$ dimers. We find that bulk structures are insulating, with a subset of FM structures showing a moderate gap ($<2$ eV) in one spin channel.


title: Hollow-Core Terahertz Optical Waveguides With Hyperuniform Disordered Reflectors

abstract: Novel hollow-core THz waveguides featuring hyperuniform disordered reflectors are proposed, fabricated, and characterized. The reflector comprise aperiodically positioned dielectric cylinders connected with dielectric bridges. The proposed waveguides are fabricated using a 3D stereolithography printer. Optical properties of the fabricated waveguide are investigated numerically using finite element method, as well as experimentally using THz-TDS spectroscopy. The results confirm that proposed waveguides exhibit sizable photonic band gaps (20%) even when relatively low refractive index contrast used (resin/air). Position of the band gaps can be easily tuned by varying reflector geometrical parameters.


title: Reflexive Polytopes and Lattice-Polarized K3 Surfaces

abstract: In this expository note, we review the standard formulation of mirror symmetry for Calabi-Yau hypersurfaces in toric varieties, and compare this construction to a description of mirror symmetry for K3 surfaces which relies on a sublattice of the Picard lattice. We then show how to combine information about the Picard group of a toric ambient space with data about automorphisms of the toric variety to identify families of K3 surfaces with high Picard rank.


title: The new modified Vlasov equation for the systems with dissipative processes

abstract: A new modified Vlasov equation has been obtained in this paper for the systems with dissipative phenomena such as, for example, plasma with irradiation.


title: Inequalities for the modified k-Bessel function

abstract: The article considers the generalized k-Bessel functions and represents it as Wright functions. Then we study the monotonicity properties of the ratio of two different orders k- Bessel functions, and the ratio of the k-Bessel and the m-Bessel functions. The log-convexity with respect to the order of the k-Bessel also given. An investigation regarding the monotonicity of the ratio of the k-Bessel and k-confluent hypergeometric functions are discussed.


title: NRQCD Confronts LHCb Data on Quarkonium Production within Jets

abstract: We analyze the recent LHCb measurement of the distribution of the fraction of the transverse momentum, $z(J/\psi)$, carried by the $J/\psi$ within a jet. LHCb data is compared to analytic calculations using the fragmenting jet function (FJF) formalism for studying $J/\psi$ in jets. Logarithms in the FJFs are resummed using DGLAP evolution. We also convolve hard QCD partonic cross sections, showered with PYTHIA, with leading order Non Relativistic Quantum Chromodynamics (NRQCD) fragmentation functions and obtain consistent results. Both approaches use Madgraph to calculate the hard process that creates the jet initiating parton. These calculations give reasonable agreement with the $z(J/\psi)$ distribution that was shown to be poorly described by default PYTHIA simulations in the LHCb paper. We compare our predictions for the $J/\psi$ distribution using various extractions of nonperturbative NRQCD long-distance matrix elements (LDMEs) in the literature. NRQCD calculations agree with LHCb data better than default PYTHIA regardless of which fit to the LDMEs is used. LDMEs from fits that focus exclusively on high transverse momentum data from colliders are in good agreement with the LHCb measurement.


title: Chen invariants for Riemannian submersions and their applications in meteorology

abstract: In this paper, an optimal inequality involving the delta curvature is exposed. An application of Riemannian submersions dealing meteorology is presented. Some characterizations about the vertical motion and the horizontal divergence are obtained.


title: Local Redundancy in SAT: Generalizations of Blocked Clauses

abstract: Clause-elimination procedures that simplify formulas in conjunctive normal form play an important role in modern SAT solving. Before or during the actual solving process, such procedures identify and remove clauses that are irrelevant to the solving result. These simplifications usually rely on so-called redundancy properties that characterize cases in which the removal of a clause does not affect the satisfiability status of a formula. One particularly successful redundancy property is that of blocked clauses, because it generalizes several other redundancy properties. To find out whether a clause is blocked---and therefore redundant---one only needs to consider its resolution environment, i.e., the clauses with which it can be resolved. For this reason, we say that the redundancy property of blocked clauses is local. In this paper, we show that there exist local redundancy properties that are even more general than blocked clauses. We present a semantic notion of blocking and prove that it constitutes the most general local redundancy property. We furthermore introduce the syntax-based notions of set-blocking and super-blocking, and show that the latter coincides with our semantic blocking notion. In addition, we show how semantic blocking can be alternatively characterized via Davis and Putnam's rule for eliminating atomic formulas. Finally, we perform a detailed complexity analysis and relate our novel redundancy properties to prominent redundancy properties from the literature.


title: Degrees of Freedom in Cached MIMO Relay Networks With Multiple Base Stations

abstract: The ability of physical layer relay caching to increase the degrees of freedom (DoF) of a single cell was recently illustrated. In this paper, we extend this result to the case of multiple cells in which a caching relay is shared among multiple non-cooperative base stations (BSs). In particular, we show that a large DoF gain can be achieved by exploiting the benefits of having a shared relay that cooperates with the BSs. We first propose a cache-assisted relaying protocol that improves the cooperation opportunity between the BSs and the relay. Next, we consider the cache content placement problem that aims to design the cache content at the relay such that the DoF gain is maximized. We propose an optimal algorithm and a near-optimal low-complexity algorithm for the cache content placement problem. Simulation results show significant improvement in the DoF gain using the proposed relay-caching protocol.


title: Optimal Non-uniform Deployments in Ultra-Dense Finite-Area Cellular Networks

abstract: Network densification and heterogenisation through the deployment of small cellular access points (picocells and femtocells) are seen as key mechanisms in handling the exponential increase in cellular data traffic. Modelling such networks by leveraging tools from Stochastic Geometry has proven particularly useful in understanding the fundamental limits imposed on network coverage and capacity by co-channel interference. Most of these works however assume infinite sized and uniformly distributed networks on the Euclidean plane. In contrast, we study finite sized non-uniformly distributed networks, and find the optimal non-uniform distribution of access points which maximises network coverage for a given non-uniform distribution of mobile users, and vice versa.


title: Surface creasing of soft elastic continua as a Kosterlitz-Thouless transition

abstract: Harnessing a model from composite materials science, we show how point-like cusped surface features arise as quasi-particle excitations, termed "ghost fibers", on the surface of a homogeneous soft elastic material. These deformations appear above a critical compressive strain at which ghost fiber dipoles unbind, analogous to vortices in the Kosterlitz-Thouless transition. Finite-length creases can be described in the same framework. Our predictions for crease surface profiles and onset strain agree with previous experiments and simulations, and further experimental tests are proposed.


title: Analysis and Optimization of fastText Linear Text Classifier

abstract: The paper [1] shows that simple linear classifier can compete with complex deep learning algorithms in text classification applications. Combining bag of words (BoW) and linear classification techniques, fastText [1] attains same or only slightly lower accuracy than deep learning algorithms [2-9] that are orders of magnitude slower. We proved formally that fastText can be transformed into a simpler equivalent classifier, which unlike fastText does not have any hidden layer. We also proved that the necessary and sufficient dimensionality of the word vector embedding space is exactly the number of document classes. These results help constructing more optimal linear text classifiers with guaranteed maximum classification capabilities. The results are proven exactly by pure formal algebraic methods without attracting any empirical data.