title: Optimal shortcut approach based on an easily obtained intermediate Hamiltonian

abstract: We present a general approach to speed up the adiabatic process without adding the traditional counterdiabatic driving (CD) Hamiltonian. The strategy is to design an easy-to-get intermediate Hamiltonian to connect the original Hamiltonian and final transitionless Hamiltonian. With final transitionless Hamiltonian, the same target can be achieved as in the adiabatic process governed by the original Hamiltonian, but in a shorter time. We apply the present approach to a three-level system, and the result shows that the final transitionless Hamiltonian usually has the same structure as the original Hamiltonian but with different time-dependent coefficients, allowing speedup to be achieved in a much easier way compared to previous methods.


title: Lie theory of multiplicative tensors

abstract: We study tensors on Lie groupoids suitably compatible with the groupoid structure, called {\em multiplicative}. Our main result gives a complete description of these objects only in terms of infinitesimal data. Special cases include the infinitesimal counterparts of multiplicative forms, multivector fields and holomorphic structures, obtained through a unifying and conceptual method. We also give a full treatment of multiplicative vector-valued forms, particularly Nijenhuis operators and related structures.


title: Provable Estimation of the Number of Blocks in Block Models

abstract: Community detection is a fundamental unsupervised learning problem for unlabeled networks which has a broad range of applications. Many community detection algorithms assume that the number of clusters $r$ is known apriori. In this paper, we propose an approach based on semi-definite relaxations, which does not require prior knowledge of model parameters like many existing convex relaxation methods and recovers the number of clusters and the clustering matrix exactly under a broad parameter regime, with probability tending to one. On a variety of simulated and real data experiments, we show that the proposed method often outperforms state-of-the-art techniques for estimating the number of clusters.


title: The $\imath \epsilon$ prescription in the SYK model

abstract: We introduce an $\imath \epsilon$ prescription for the SYK model both at finite and at zero temperature. This prescription regularizes all the naive ultraviolet divergences of the model. As expected the prescription breaks the conformal invariance, but the latter is restored in the $\epsilon \to 0$ limit. We prove rigorously that the Schwinger Dyson equation of the resummed two point function at large $N$ and low momentum is recovered in this limit. Based on this $\imath \epsilon$ prescription we introduce an effective field theory Lagrangian for the infrared SYK model.


title: On the multiply robust estimation of the mean of the g-functional

abstract: We study multiply robust (MR) estimators of the longitudinal g-computation formula of Robins (1986). In the first part of this paper we review and extend the recently proposed parametric multiply robust estimators of Tchetgen-Tchetgen (2009) and Molina, Rotnitzky, Sued and Robins (2017). In the second part of the paper we derive multiply and doubly robust estimators that use non-parametric machine-learning (ML) estimators of nuisance functions in lieu of parametric models. We use sample splitting to avoid the need for Donsker conditions, thereby allowing an analyst to select the ML algorithms of their choosing. We contrast the asymptotic behavior of our non-parametric doubly robust and multiply robust estimators. In particular, we derive formulas for their asymptotic bias. Examining these formulas we conclude that although, under certain data generating laws, the rate at which the bias of the MR estimator converges to zero can exceed that of the DR estimator, nonetheless, under most laws, the bias of the DR and MR estimators converge to zero at the same rate.


title: Sequence Summarization Using Order-constrained Kernelized Feature Subspaces

abstract: Representations that can compactly and effectively capture temporal evolution of semantic content are important to machine learning algorithms that operate on multi-variate time-series data. We investigate such representations motivated by the task of human action recognition. Here each data instance is encoded by a multivariate feature (such as via a deep CNN) where action dynamics are characterized by their variations in time. As these features are often non-linear, we propose a novel pooling method, kernelized rank pooling, that represents a given sequence compactly as the pre-image of the parameters of a hyperplane in an RKHS, projections of data onto which captures their temporal order. We develop this idea further and show that such a pooling scheme can be cast as an order-constrained kernelized PCA objective; we then propose to use the parameters of a kernelized low-rank feature subspace as the representation of the sequences. We cast our formulation as an optimization problem on generalized Grassmann manifolds and then solve it efficiently using Riemannian optimization techniques. We present experiments on several action recognition datasets using diverse feature modalities and demonstrate state-of-the-art results.


title: MMD GAN: Towards Deeper Understanding of Moment Matching Network

abstract: Generative moment matching network (GMMN) is a deep generative model that differs from Generative Adversarial Network (GAN) by replacing the discriminator in GAN with a two-sample test based on kernel maximum mean discrepancy (MMD). Although some theoretical guarantees of MMD have been studied, the empirical performance of GMMN is still not as competitive as that of GAN on challenging and large benchmark datasets. The computational efficiency of GMMN is also less desirable in comparison with GAN, partially due to its requirement for a rather large batch size during the training. In this paper, we propose to improve both the model expressiveness of GMMN and its computational efficiency by introducing adversarial kernel learning techniques, as the replacement of a fixed Gaussian kernel in the original GMMN. The new approach combines the key ideas in both GMMN and GAN, hence we name it MMD GAN. The new distance measure in MMD GAN is a meaningful loss that enjoys the advantage of weak topology and can be optimized via gradient descent with relatively small batch sizes. In our evaluation on multiple benchmark datasets, including MNIST, CIFAR- 10, CelebA and LSUN, the performance of MMD-GAN significantly outperforms GMMN, and is competitive with other representative GAN works.


title: Elliptic function solutions in Jackiw-Teitelboim dilaton gravity

abstract: We present a new family of solutions for the Jackiw-Teitelboim model of two-dimensional gravity with a negative cosmological constant. Here, a metric of constant Ricci scalar curvature is constructed, and explicit linearly independent solutions of the corresponding dilaton field equations are determined. The metric is transformed to a black hole metric, and the dilaton solutions are expressed in terms of Jacobi elliptic functions. Using these solutions we compute, for example, Killing vectors for the metric.


title: Self-organized Segregation on the Grid

abstract: We consider an agent-based model in which two types of agents interact locally over a graph and have a common intolerance threshold $\tau$ for changing their types with exponentially distributed waiting times. The model is equivalent to an unperturbed Schelling model of self-organized segregation, an Asynchronous Cellular Automata (ACA) with extended Moore neighborhoods, or a zero-temperature Ising model with Glauber dynamics, and has applications in the analysis of social and biological networks, and spin glasses systems. Some rigorous results were recently obtained in the theoretical computer science literature, and this work provides several extensions. We enlarge the intolerance interval leading to the formation of large segregated regions of agents of a single type from the known size $\epsilon>0$ to size $\approx 0.134$. Namely, we show that for $0.433 < \tau < 1/2$ (and by symmetry $1/2<\tau<0.567$), the expected size of the largest segregated region containing an arbitrary agent is exponential in the size of the neighborhood. We further extend the interval leading to large segregated regions to size $\approx 0.312$ considering "almost segregated" regions, namely regions where the ratio of the number of agents of one type and the number of agents of the other type vanishes quickly as the size of the neighborhood grows. In this case, we show that for $0.344 < \tau \leq 0.433$ (and by symmetry for $0.567 \leq \tau<0.656$) the expected size of the largest almost segregated region containing an arbitrary agent is exponential in the size of the neighborhood. The exponential bounds that we provide also imply that complete segregation, where agents of a single type cover the whole grid, does not occur with high probability for $p=1/2$ and the range of tolerance considered.


title: Ages and structural and dynamical parameters of two globular clusters in the M81 group

abstract: GC-1 and GC-2 are two globular clusters (GCs) in the remote halo of M81 and M82 in the M81 group discovered by Jang et al. using the {\it Hubble Space Telescope} ({\it HST}) images. These two GCs were observed as part of the Beijing--Arizona--Taiwan--Connecticut (BATC) Multicolor Sky Survey, using 14 intermediate-band filters covering a wavelength range of 4000--10000 \AA. We accurately determine these two clusters' ages and masses by comparing their spectral energy distributions (from 2267 to 20000~{\AA}, comprising photometric data in the near-ultraviolet of the {\it Galaxy Evolution Explorer}, 14 BATC intermediate-band, and Two Micron All Sky Survey near-infrared $JHK_{\rm s}$ filters) with theoretical stellar population-synthesis models, resulting in ages of $15.50\pm3.20$ for GC-1 and $15.10\pm2.70$ Gyr for GC-2. The masses of GC-1 and GC-2 obtained here are $1.77-2.04\times 10^6$ and $5.20-7.11\times 10^6 \rm~M_\odot$, respectively. In addition, the deep observations with the Advanced Camera for Surveys and Wide Field Camera 3 on the {\it HST} are used to provide the surface brightness profiles of GC-1 and GC-2. The structural and dynamical parameters are derived from fitting the profiles to three different models; in particular, the internal velocity dispersions of GC-1 and GC-2 are derived, which can be compared with ones obtained based on spectral observations in the future. For the first time, in this paper, the $r_h$ versus $M_V$ diagram shows that GC-2 is an ultra-compact dwarf in the M81 group.


title: A universal scaling law of black hole activity including gamma-ray bursts

abstract: Previous works show that a correlation among radio luminosity, X-ray luminosity, and black hole (BH) mass from stellar-mass BHs in X-ray binaries to supermassive BHs in active galactic nuclei (AGNs), which leads to the so-called fundamental plane of BH activity. However, there are two competing explanations for this fundamental plane, including the jet-dominated model and the disk-jet model. Thus, the physical origin of this fundamental plane remains unknown. In this paper, we show that the X-ray luminosities, radio luminosities and BH masses of gamma-ray bursts (GRBs) and M82 X-1 also show a similar distribution. The universal scaling law among stellar-mass, intermediate and supermassive BH systems, together with the fact that radio and X-ray emission of GRBs originates from relativistic jets, reveals that the fundamental plane of BH activity is controlled by a jet, i.e., the radio and X-ray emission is mainly from the jet. Our work also suggests that the jets are scale-invariant with respect to the BH mass.


title: Metal surface structuring with spatiotemporally focused femtosecond laser pulses

abstract: Femtosecond laser micromachining provides high precision and less thermal diffusion in surface structuring as a result of the ultrashort temporal duration and ultrahigh peak intensity of the femtosecond laser pulses. To increase the throughput of surface patterning, the focal spot size can be expanded with loose focusing, which, however, could lead to nonlinear self-focusing of the pulses when the pulses propagate in air. We solve the problem by use of spatiotemporally focused femtosecond laser pulses for ablation of metal surfaces, which gives rise to improved surface quality as compared with that obtained with the conventional focusing scheme.


title: Generative Model with Coordinate Metric Learning for Object Recognition Based on 3D Models

abstract: Given large amount of real photos for training, Convolutional neural network shows excellent performance on object recognition tasks. However, the process of collecting data is so tedious and the background are also limited which makes it hard to establish a perfect database. In this paper, our generative model trained with synthetic images rendered from 3D models reduces the workload of data collection and limitation of conditions. Our structure is composed of two sub-networks: semantic foreground object reconstruction network based on Bayesian inference and classification network based on multi-triplet cost function for avoiding over-fitting problem on monotone surface and fully utilizing pose information by establishing sphere-like distribution of descriptors in each category which is helpful for recognition on regular photos according to poses, lighting condition, background and category information of rendered images. Firstly, our conjugate structure called generative model with metric learning utilizing additional foreground object channels generated from Bayesian rendering as the joint of two sub-networks. Multi-triplet cost function based on poses for object recognition are used for metric learning which makes it possible training a category classifier purely based on synthetic data. Secondly, we design a coordinate training strategy with the help of adaptive noises acting as corruption on input images to help both sub-networks benefit from each other and avoid inharmonious parameter tuning due to different convergence speed of two sub-networks. Our structure achieves the state of the art accuracy of over 50\% on ShapeNet database with data migration obstacle from synthetic images to real photos. This pipeline makes it applicable to do recognition on real images only based on 3D models.


title: Invariant-based pulse design for three-level systems without the rotating-wave approximation

abstract: In this paper, a scheme is put forward to design pulses which drive a three-level system based on the reverse engineering with Lewis-Riesenfeld invariant theory. The scheme can be applied to a three-level system even when the rotating-wave approximation (RWA) can not be used. The amplitudes of pulses and the maximal values of detunings in the system could be easily controlled by adjusting control parameters. We analyze the dynamics of the system by an invariant operator, so additional couplings are unnecessary. Moreover, the approaches to avoid singularity of pulses are studied and several useful results are obtained. We hope the scheme could contribute to fast quantum information processing without RWA.


title: Sufficient conditions for the existence of a path-factor which are related to odd components

abstract: In this paper, we are concerned with sufficient conditions for the existence of a $\{P_{2},P_{2k+1}\}$-factor. We prove that for $k\geq 3$, there exists $\varepsilon_{k}>0$ such that if a graph $G$ satisfies $\sum_{0\leq j\leq k-1}c_{2j+1}(G-X)\leq \varepsilon_{k}|X|$ for all $X\subseteq V(G)$, then $G$ has a $\{P_{2},P_{2k+1}\}$-factor, where $c_{i}(G-X)$ is the number of components $C$ of $G-X$ with $|V(C)|=i$. On the other hand, we construct infinitely many graphs $G$ having no $\{P_{2},P_{2k+1}\}$-factor such that $\sum_{0\leq j\leq k-1}c_{2j+1}(G-X)\leq \frac{32k+141}{72k-78}|X|$ for all $X\subseteq V(G)$.


title: Deep Learning Improves Template Matching by Normalized Cross Correlation

abstract: Template matching by normalized cross correlation (NCC) is widely used for finding image correspondences. We improve the robustness of this algorithm by preprocessing images with "siamese" convolutional networks trained to maximize the contrast between NCC values of true and false matches. The improvement is quantified using patches of brain images from serial section electron microscopy. Relative to a parameter-tuned bandpass filter, siamese convolutional networks significantly reduce false matches. Furthermore, all false matches can be eliminated by removing a tiny fraction of all matches based on NCC values. The improved accuracy of our method could be essential for connectomics, because emerging petascale datasets may require billions of template matches to assemble 2D images of serial sections into a 3D image stack. Our method is also expected to generalize to many other computer vision applications that use NCC template matching to find image correspondences.


title: Superconformal Blocks for SCFTs with Eight Supercharges

abstract: We show how to treat the superconformal algebras with eight Poincar\'e supercharges in a unified manner for spacetime dimension $2 < d\leq 6$. This formalism is ideally suited for analyzing the quadratic Casimir operator of the superconformal algebra and its use in deriving superconformal blocks. We illustrate this by an explicit construction of the superconformal blocks, for any value of the spacetime dimension, for external protected scalar operators which are the lowest component of flavor current multiplets.


title: Bilinear estimates in Besov spaces generated by the Dirichlet Laplacian

abstract: The purpose of this paper is to establish bilinear estimates in Besov spaces generated by the Dirichlet Laplacian on a domain of Euclidian spaces. These estimates are proved by using the gradient estimates for heat semigroup together with the Bony paraproduct formula and the boundedness of spectral multipliers.


title: Critical phenomena and chemical potential of charged AdS black hole

abstract: We study the thermodynamics and the chemical potential for a five-dimensional charged AdS black hole by treating the cosmological constant as the number of colors $N$ in the boundary gauge theory and its conjugate quantity as the associated chemical potential $\mu$. It is found that there exists a small-large black hole phase transition. The critical phenomena are investigated in the $N^{2}$-$\mu$ chart. In particular, in the reduced parameter space, all the thermodynamic quantities can be rescaled with the black hole charge such that these reduced quantities are charge-independent. Then we obtain the coexistence curve and the phase diagram. The latent heat is also numerically calculated. Moreover, the heat capacity and the thermodynamic scalar are studied. The result indicates that the information of the first-order black hole phase transition is encoded in the heat capacity and scalar. However, the phase transition point cannot be directly calculated with them. Nevertheless, the critical point linked to a second-order phase transition can be determined by either the heat capacity or the scalar. In addition, we calculate the critical exponents of the heat capacity and the scalar for the saturated small and large black holes near the critical point.


title: Modal phenomena of surface and bulk polaritons in magnetic-semiconductor superlattices

abstract: We discuss peculiarities of bulk and surface polaritons propagating in a composite magnetic-semiconductor superlattice influenced by an external static magnetic field. Three particular configurations of magnetization, namely the Voigt, polar and Faraday geometries are considered. In the long-wavelength limit, involving the effective medium theory, the proposed superlattice is described as an anisotropic uniform medium defined by the tensors of effective permittivity and effective permeability. The study is carried out in the frequency band where the characteristic resonant frequencies of the underlying constitutive magnetic and semiconductor materials of the superlattice are different but closely spaced. The effects of mode crossing and anti-crossing in dispersion characteristics of both bulk and surface polaritons are revealed and explained with an assistance of the concept of Morse critical points from the catastrophe theory.


title: Journalists' information needs, seeking behavior, and its determinants on social media

abstract: We describe the results of a qualitative study on journalists' information seeking behavior on social media. Based on interviews with eleven journalists along with a study of a set of university level journalism modules, we determined the categories of information need types that lead journalists to social media. We also determined the ways that social media is exploited as a tool to satisfy information needs and to define influential factors, which impacted on journalists' information seeking behavior. We find that not only is social media used as an information source, but it can also be a supplier of stories found serendipitously. We find seven information need types that expand the types found in previous work. We also find five categories of influential factors that affect the way journalists seek information.


title: X-Ray Amplification by Stimulated Brillouin Scattering

abstract: Plasma-based parametric amplification using stimulated Brillouin scattering offers a route to coherent x-ray pulses orders-of-magnitude more intense than those of the brightest available sources. Brillouin amplification permits amplification of shorter wavelengths with lower pump intensities than Raman amplification, which Landau and collisional damping limit in the x-ray regime. Analytic predictions, numerical solutions of the three-wave coupling equations, and particle-in-cell simulations suggest that Brillouin amplification in solid-density plasmas will allow compression of current x-ray free electron laser pulses to sub-femtosecond durations and unprecedented intensities.


title: The Haldane model under nonuniform strain

abstract: We study the Haldane model under strain using a tight-binding approach, and compare the obtained results with the continuum-limit approximation. As in graphene, nonuniform strain leads to a time-reversal preserving pseudo-magnetic field that induces (pseudo) Landau levels. Unlike a real magnetic field, strain lifts the degeneracy of the zeroth pseudo Landau levels at different valleys. Moreover, for the zigzag edge under uniaxial strain, strain removes the degeneracy within the pseudo-Landau levels by inducing a tilt in their energy dispersion. The latter arises from next-to-leading order corrections to the continuum-limit Hamiltonian, which are absent for a real magnetic field. We show that, for the lowest pseudo-Landau levels in the Haldane model, the dominant contribution to the tilt is different from graphene. In addition, although strain does not strongly modify the dispersion of the edge states, their interplay with the pseudo-Landau levels is different for the armchair and zigzag ribbons. Finally, we study the effect of strain in the band structure of the Haldane model at the critical point of the topological transition, thus shedding light on the interplay between non-trivial topology and strain in quantum anomalous Hall systems.


title: Testing violation of the Leggett-Garg-type inequality in neutrino oscillations of the Daya Bay experiment

abstract: The Leggett-Garg inequality (LGI), derived under the assumption of realism, acts as the temporal Bell's inequality. It is studied in electromagnetic and strong interaction like photonics, superconducting qu-bits and nuclear spin. Until the weak interaction two-state oscillations of neutrinos affirmed the violation of Leggett-Garg-type inequalities (LGtI). We make an empirical test for the deviation of experimental results with the classical limits by analyzing the survival probability data of reactor neutrinos at a distinct range of baseline dividing energies, as an analog to a single neutrino detected at different time. A study of the updated data of Daya-Bay experiment unambiguously depicts an obvious cluster of data over the classical bound of LGtI and shows a $6.1\sigma$ significance of the violation of them.


title: Light bending, static dark energy and related uniqueness of Schwarzschild-de Sitter spacetime

abstract: Since the Schwarzschild-de Sitter spacetime is static inside the cosmological event horizon, if the dark energy state parameter is sufficiently close to $-1$, apparently one could still expect an effectively static geometry, in the attraction dominated region inside the maximum turn around radius, $R_{\rm TA, max}$, of a cosmic structure. We take the first order metric derived recently assuming a static and ideal dark energy fluid with equation of state $P(r)=\alpha\rho(r)$ as a source in Ref. [1], which reproduced the expression for $R_{\rm TA, max}$ found earlier in the cosmological McVittie spacetime. Here we show that the equality originates from the equivalence of geodesic motion in these two backgrounds, in the non-relativistic regime. We extend this metric up to the third order and compute the bending of light using the Rindler-Ishak method. For $ \alpha\neq -1$, a dark energy dependent term appears in the bending equation, unlike the case of the cosmological constant, $\alpha=-1$. Due to this new term in particular, existing data for the light bending at galactic scales yields, $(1+\alpha)\lesssim {\cal O}(10^{-14})$, thereby practically ruling out any such static and inhomogeneous dark energy fluid we started with. Implication of this result pertaining the uniqueness of the Schwarzschild-de Sitter spacetime in such inhomogeneous dark energy background is discussed.


title: Control of single-mode operation in a circular waveguide filled by a longitudinally magnetized gyroelectromagnetic medium

abstract: A substantial control of dispersion features of the hybrid EH01 and HE11 modes of a circular waveguide which is completely filled by a longitudinally magnetized composite finely-stratified ferrite-semiconductor structure is discussed. A relation between the resonant conditions of such a composite gyroelectromagnetic filling of the circular waveguide and dispersion features of the supported modes are studied. Three distinct frequency bands with the single-mode operation under normal as well as anomalous dispersion conditions of the EH01 mode are identified by solving an optimization problem with respect to the filling factors of the composite medium. The possibility of achieving isolated propagation of the HE11 mode is revealed.


title: Stochastic Game in Remote Estimation under DoS Attacks

abstract: This paper studies remote state estimation under denial-of-service (DoS) attacks. A sensor transmits its local estimate of an underlying physical process to a remote estimator via a wireless communication channel. A DoS attacker is capable to interfere the channel and degrades the remote estimation accuracy. Considering the tactical jamming strategies played by the attacker, the sensor adjusts its transmission power. This interactive process between the sensor and the attacker is studied in the framework of a zero-sum stochastic game. To derive their optimal power schemes, we first discuss the existence of stationary Nash equilibrium (SNE) for this game. We then present the monotone structure of the optimal strategies, which helps reduce the computational complexity of the stochastic game algorithm. Numerical examples are provided to illustrate the obtained results.


title: Transiting Exoplanet Monitoring Project (TEMP). II. Refined System Parameters and Transit Timing Analysis of HAT-P-33b

abstract: We present ten $R$-band photometric observations of eight different transits of the hot Jupiter HAT-P-33b, which has been targeted by our Transiting Exoplanet Monitoring Project (TEMP). The data were obtained by two telescopes at the Xinglong Station of National Astronomical Observatories of China (NAOC) from 2013 December through 2016 January, and exhibit photometric scatter of $1.6-3.0\,\rm{mmag}$. After jointly analyzing the previously published photometric data, radial-velocity (RV) measurements, and our new light curves, we revisit the system parameters and orbital ephemeris for the HAT-P-33b system. Our results are consistent with the published values except for the planet-to-star radius ratio ($R_{P}/R_{*}$), the ingress/egress duration ($\tau$) and the total duration ($T_{14}$), which together indicate a slightly shallower and shorter transit shape. Our results are based on more complete light curves, whereas the previously published work had only one complete transit light curve. No significant anomalies in Transit Timing Variations (TTVs) are found, and we place upper mass limits on potential perturbers, largely supplanting the loose constraints provided by the extant RV data. The TTV limits are stronger near mean-motion resonances, especially for the low-order commensurabilities. We can exclude the existence of a perturber with mass larger than 0.6, 0.3, 0.5, 0.5, and $0.3\,{\rm M_\oplus}$ near the 1:3, 1:2, 2:3, 3:2, and 2:1 resonances, respectively.


title: Phase transitions and spin excitations of spin-1 bosons in optical lattice

abstract: We investigate ground state properties of spin-1 bosonic system trapped in optical lattice with extended standard basis operator (SBO) method. For both ferromagnetic ($U_2<0$) and antiferromagnetic ($U_2>0$) systems, we analytically figure out the symmetry properties in Mott-insulator and superfluid phases, which would provide a deeper insight into the MI-SF phase transition process. Then by applying self-consistent approach to the method, we include the effect of quantum and thermal fluctuations and derive the MI-SF transition phase diagram, which is in quantitative agreement with recent Monte-Carlo simulation at zero temperature, and at finite temperature, we find the underestimation of finite-temperature-effect in the mean-field approximation method. If we further consider the spin excitations in the insulating states of spin-1 system in external field, distinct spin phases are expected. Therefore, in the Mott lobes with $n=1$ and $n=2$ atoms per site, we give analytical and numerical boundaries of the singlet, nematic, partially magnetic and ferromagnetic phases in the magnetic phase diagrams.


title: Substitution invariant Sturmian words and binary trees

abstract: We take a global view at substitution invariant Sturmian sequences. We show that homogeneous substitution invariant Sturmian sequences $s_{\alpha,\alpha}$ can be indexed by two binary trees, associated directly to Johannes Kepler's tree of harmonic fractions from 1619. We obtain similar results for the inhomogeneous sequences $s_{\alpha,1-\alpha}$ and $s_{\alpha,0}$.


title: Formation of solar system analogs I: looking for initial conditions through a population synthesis analysis

abstract: Population synthesis models of planetary systems developed during the last $\sim$15 years could reproduce several of the observables of the exoplanet population, and also allowed to constrain planetary formation models. We present our planet formation model, which calculates the evolution of a planetary system during the gaseous phase. The code incorporates relevant physical phenomena for the formation of a planetary system, like photoevaporation, planet migration, gas accretion, water delivery in embryos and planetesimals, a detailed study of the orbital evolution of the planetesimal population, and the treatment of the fusion between embryos, considering their atmospheres. The main goal of this work, unlike other works of planetary population synthesis, is to find suitable scenarios and physical parameters of the disc to form solar system analogs. We are specially interested in the final planet distributions, and in the final surface density, eccentricity and inclination profiles for the planetesimal population. These final distributions will be used as initial conditions for N-body simulations, to study the post-oligarchic formation in a second work. We then consider different formation scenarios, with different planetesimal sizes and different type I migration rates. We find that solar system analogs are favored in massive discs, with low type I migration rates, and small planetesimal sizes. Besides, those rocky planets within their habitables zones are dry when discs dissipate. At last, the final configurations of solar system analogs include information about the mass and semimajor-axis of the planets, water contents, and the properties of the planetesimal remnants.


title: Multisymplecticity of hybridizable discontinuous Galerkin methods

abstract: In this paper, we prove necessary and sufficient conditions for a hybridizable discontinuous Galerkin (HDG) method to satisfy a multisymplectic conservation law, when applied to a canonical Hamiltonian system of partial differential equations. We show that these conditions are satisfied by the "hybridized" versions of several of the most commonly-used finite element methods, including mixed, nonconforming, and discontinuous Galerkin methods. (Interestingly, for the continuous Galerkin method in dimension greater than one, we show that multisymplecticity only holds in a weaker sense.) Consequently, these general-purpose finite element methods may be used for structure-preserving discretization (or semidiscretization) of canonical Hamiltonian systems of ODEs or PDEs. This establishes multisymplecticity for a large class of arbitrarily-high-order methods on unstructured meshes.


title: The Southern HII Region Discovery Survey (SHRDS): Pilot Survey

abstract: The Southern HII Region Discovery Survey (SHRDS) is a survey of the third and fourth quadrants of the Galactic plane that will detect radio recombination line and continuum emission at cm-wavelengths from several hundred HII region candidates using the Australia Telescope Compact Array. The targets for this survey come from the WISE Catalog of Galactic HII Regions, and were identified based on mid-infrared and radio continuum emission. In this pilot project, two different configurations of the Compact Array Broad Band receiver and spectrometer system were used for short test observations. The pilot surveys detected radio recombination line emission from 36 of 53 HII region candidates, as well as seven known \hii regions that were included for calibration. These 36 recombination line detections confirm that the candidates are true HII regions, and allow us to estimate their distances.


title: Boost invariant formulation of the chiral kinetic theory

abstract: We formulate the chiral kinetic equation with the longitudinal boost invariance. We particularly focus on the physical interpretation of the particle number conservation. There appear two terms associated with the expansion, which did not exist in the non-chiral kinetic equation. One is a contribution to the transverse current arising from the side-jump effect, and the other is a change in the density whose flow makes the longitudinal current. We point out a characteristic pattern in the transverse current driven by the expansion, which we call the chiral circular displacement.


title: Age spreads and the temperature dependence of age estimates in Upper Sco

abstract: Past estimates for the age of the Upper Sco Association are typically 11-13 Myr for intermediate-mass stars and 4-5 Myr for low-mass stars. In this study, we simulate populations of young stars to investigate whether this apparent dependence of estimated age on spectral type may be explained by the star formation history of the association. Solar and intermediate mass stars begin their pre-main sequence evolution on the Hayashi track, with fully convective interiors and cool photospheres. Intermediate mass stars quickly heat up and transition onto the radiative Henyey track. As a consequence, for clusters in which star formation occurs on a similar timescale as the transition from a convective to a radiative interior, discrepancies in ages will arise when ages are calculated as a function of temperature instead of mass. Simple simulations of a cluster with constant star formation over several Myr may explain about half of the difference in inferred ages versus photospheric temperature; speculative constructions that consist of a constant star formation followed by a large supernova-driven burst could fully explain the differences, including those between F and G stars where evolutionary tracks may be more accurate. The age spreads of low-mass stars predicted from these prescriptions for star formation are consistent with the observed luminosity spread of Upper Sco. The conclusion that a lengthy star formation history will yield a temperature dependence in ages is expected from the basic physics of pre-main sequence evolution and is qualitatively robust to the large uncertainties in pre-main sequence evolutionary models.


title: Crystal Structure Prediction Supported by Incomplete Experimental Data

abstract: The prediction of material structure from chemical composition has been a long-standing challenge in natural science. Although there have been various methodological developments and successes with computer simulations, the prediction of crystal structures comprising more than several tens of atoms in the unit cell still remains difficult due to the many degrees of freedom, which increase exponentially with the number of atoms. Here we show that when some experimental data is available, even if it is totally insufficient for conventional structure analysis, it can be utilized to support and substantially accelerate structure simulation. In particular, we formulate a cost function based on a weighted sum of interatomic potential energies and a penalty function referred to as "crystallinity", which is defined using limited X-ray diffraction data. This method is applied to well-known polymorphs of $\rm{SiO_2}$ with up to 96 atoms in the simulation cell to find that it reproduces the correct structures efficiently with a very limited number of diffraction peaks. The penalty function is confirmed to destabilize the local minima of the potential energy surface, which facilitates finding the correct structure. This method opens a new avenue for determining and predicting structures that are difficult to determine by conventional methods, such as surface, interface, glass, and amorphous structures.


title: Fully reliable error control for evolutionary problems

abstract: This work is focused on the application of functional-type a posteriori error estimates and corresponding indicators to a class of time-dependent problems. We consider the algorithmic part of their derivation and implementation and also discuss the numerical properties of these bounds that comply with obtained numerical results. This paper examines two different methods of approximate solution reconstruction for evolutionary models, i.e., a time-marching technique and a space-time approach. The first part of the study presents an algorithm for global minimization of the majorant on each of discretization time-cylinders (time-slabs), the effectiveness of this algorithm is confirmed by extensive numerical tests. In the second part of the publication, the application of functional error estimates is discussed with respect to a space-time approach. It is followed by a set of extensive numerical tests that demonstrate the efficiency of proposed error control method.


title: Sharp Threshold of Blow-up and Scattering for the fractional Hartree equation

abstract: We consider the fractional Hartree equation in the $L^2$-supercritical case, and we find a sharp threshold of the scattering versus blow-up dichotomy for radial data: If $ M[u_{0}]^{\frac{s-s_c}{s_c}}E[u_{0}<M[Q]^{\frac{s-s_c}{s_c}}E[Q]$ and $M[u_{0}]^{\frac{s-s_c}{s_c}}\|u_{0}\|^2_{\dot H^s}<M[Q]^{\frac{s-s_c}{s_c}}\| Q\|^2_{\dot H^s}$, then the solution $u(t)$ is globally well-posed and scatters; if $ M[u_{0}]^{\frac{s-s_c}{s_c}}E[u_{0}]<M[Q]^{\frac{s-s_c}{s_c}}E[Q]$ and $M[u_{0}]^{\frac{s-s_c}{s_c}}\|u_{0}\|^2_{\dot H^s}>M[Q]^{\frac{s-s_c}{s_c}}\| Q\|^2_{\dot H^s}$, the solution $u(t)$ blows up in finite time. This condition is sharp in the sense that the solitary wave solution $e^{it}Q(x)$ is global but not scattering, which satisfies the equality in the above conditions. Here, $Q$ is the ground-state solution for the fractional Hartree equation.


title: A New Parallel Message-distribution Technique for Cost-based Steganography

abstract: This paper presents two novel approaches to increase performance bounds of image steganography under the criteria of minimizing distortion. First, in order to efficiently use the images' capacities, we propose using parallel images in the embedding stage. The result is then used to prove sub-optimality of the message distribution technique used by all cost based algorithms including HUGO, S-UNIWARD, and HILL. Second, a new distribution approach is presented to further improve the security of these algorithms. Experiments show that this distribution method avoids embedding in smooth regions and thus achieves a better performance, measured by state-of-the-art steganalysis, when compared with the current used distribution.


title: Which bridge estimator is optimal for variable selection?

abstract: We study the problem of variable selection for linear models under the high-dimensional asymptotic setting, where the number of observations $n$ grows at the same rate as the number of predictors $p$. We consider two-stage variable selection techniques (TVS) in which the first stage uses bridge estimators to obtain an estimate of the regression coefficients, and the second stage simply thresholds this estimate to select the "important" predictors. The asymptotic false discovery proportion (AFDP) and true positive proportion (ATPP) of these TVS are evaluated. We prove that for a fixed ATPP, in order to obtain a smaller AFDP, one should pick a bridge estimator with smaller asymptotic mean square error in the first stage of TVS. Based on such principled discovery, we present a sharp comparison of different TVS, via an in-depth investigation of the estimation properties of bridge estimators. Rather than "order-wise" error bounds with loose constants, our analysis focuses on precise error characterization. Various interesting signal-to-noise ratio and sparsity settings are studied. Our results offer new and thorough insights into high-dimensional variable selection. For instance, we prove that a TVS with Ridge in its first stage outperforms TVS with other bridge estimators in large noise settings; two-stage LASSO becomes inferior when the signal is rare and weak. As a by-product, we show that two-stage methods outperform some standard variable selection techniques, such as LASSO and Sure Independence Screening, under certain conditions.


title: Multi-Task Learning for Contextual Bandits

abstract: Contextual bandits are a form of multi-armed bandit in which the agent has access to predictive side information (known as the context) for each arm at each time step, and have been used to model personalized news recommendation, ad placement, and other applications. In this work, we propose a multi-task learning framework for contextual bandit problems. Like multi-task learning in the batch setting, the goal is to leverage similarities in contexts for different arms so as to improve the agent's ability to predict rewards from contexts. We propose an upper confidence bound-based multi-task learning algorithm for contextual bandits, establish a corresponding regret bound, and interpret this bound to quantify the advantages of learning in the presence of high task (arm) similarity. We also describe an effective scheme for estimating task similarity from data, and demonstrate our algorithm's performance on several data sets.


title: Dictionary-based Monitoring of Premature Ventricular Contractions: An Ultra-Low-Cost Point-of-Care Service

abstract: While cardiovascular diseases (CVDs) are prevalent across economic strata, the economically disadvantaged population is disproportionately affected due to the high cost of traditional CVD management. Accordingly, developing an ultra-low-cost alternative, affordable even to groups at the bottom of the economic pyramid, has emerged as a societal imperative. Against this backdrop, we propose an inexpensive yet accurate home-based electrocardiogram(ECG) monitoring service. Specifically, we seek to provide point-of-care monitoring of premature ventricular contractions (PVCs), high frequency of which could indicate the onset of potentially fatal arrhythmia. Note that a traditional telecardiology system acquires the ECG, transmits it to a professional diagnostic centre without processing, and nearly achieves the diagnostic accuracy of a bedside setup, albeit at high bandwidth cost. In this context, we aim at reducing cost without significantly sacrificing reliability. To this end, we develop a dictionary-based algorithm that detects with high sensitivity the anomalous beats only which are then transmitted. We further compress those transmitted beats using class-specific dictionaries subject to suitable reconstruction/diagnostic fidelity. Such a scheme would not only reduce the overall bandwidth requirement, but also localising anomalous beats, thereby reducing physicians' burden. Finally, using Monte Carlo cross validation on MIT/BIH arrhythmia database, we evaluate the performance of the proposed system. In particular, with a sensitivity target of at most one undetected PVC in one hundred beats, and a percentage root mean squared difference less than 9% (a clinically acceptable level of fidelity), we achieved about 99.15% reduction in bandwidth cost, equivalent to 118-fold savings over traditional telecardiology.


title: Robust Data Geometric Structure Aligned Close yet Discriminative Domain Adaptation

abstract: Domain adaptation (DA) is transfer learning which aims to leverage labeled data in a related source domain to achieve informed knowledge transfer and help the classification of unlabeled data in a target domain. In this paper, we propose a novel DA method, namely Robust Data Geometric Structure Aligned, Close yet Discriminative Domain Adaptation (RSA-CDDA), which brings closer, in a latent joint subspace, both source and target data distributions, and aligns inherent hidden source and target data geometric structures while performing discriminative DA in repulsing both interclass source and target data. The proposed method performs domain adaptation between source and target in solving a unified model, which incorporates data distribution constraints, in particular via a nonparametric distance, i.e., Maximum Mean Discrepancy (MMD), as well as constraints on inherent hidden data geometric structure segmentation and alignment between source and target, through low rank and sparse representation. RSA-CDDA achieves the search of a joint subspace in solving the proposed unified model through iterative optimization, alternating Rayleigh quotient algorithm and inexact augmented Lagrange multiplier algorithm. Extensive experiments carried out on standard DA benchmarks, i.e., 16 cross-domain image classification tasks, verify the effectiveness of the proposed method, which consistently outperforms the state-of-the-art methods.


title: Nonparametric Preference Completion

abstract: We consider the task of collaborative preference completion: given a pool of items, a pool of users and a partially observed item-user rating matrix, the goal is to recover the \emph{personalized ranking} of each user over all of the items. Our approach is nonparametric: we assume that each item $i$ and each user $u$ have unobserved features $x_i$ and $y_u$, and that the associated rating is given by $g_u(f(x_i,y_u))$ where $f$ is Lipschitz and $g_u$ is a monotonic transformation that depends on the user. We propose a $k$-nearest neighbors-like algorithm and prove that it is consistent. To the best of our knowledge, this is the first consistency result for the collaborative preference completion problem in a nonparametric setting. Finally, we demonstrate the performance of our algorithm with experiments on the Netflix and Movielens datasets.


title: Effect of Scattering Efficiency in the Tip Enhanced Raman Spectroscopic Imaging of Nanostructures in the Sub Diffraction Limit

abstract: The experimental limitations in the signal enhancement and spatial resolution in spectroscopic imaging have been always a challenging task in the application of near-field spectroscopy for nanostructured materials in the sub-diffraction limit. In addition, the scattering efficiency also plays an important role in improving signal enhancement and contrast of the spectroscopic imaging of nanostructures by scattering of light. We report the effect of scattering efficiency in the Raman intensity enhancement, and contrast generation in near-field tip-enhanced Raman spectroscopic (TERS) imaging of one dimensional inorganic crystalline nanostructures of Si and AlN having a large variation in polarizability change. The Raman enhancement of pure covalently bonded Si nanowire (NW) is found to be two orders of higher in magnitude for the TERS imaging, as compared to that of AlN nanotip (NT) having a higher degree of ionic bonding, suggesting the importance of scattering efficiency of the materials in TERS imaging. The strong contrast generation due to higher signal enhancement in TERS imaging of Si NW also helped in achieving the better resolved spectroscopic images than that of the AlN NT. The study provides an insight into the role of scattering efficiency in the resolution of near-field spectroscopic images.


title: Deep Rotation Equivariant Network

abstract: Recently, learning equivariant representations has attracted considerable research attention. Dieleman et al. introduce four operations which can be inserted into convolutional neural network to learn deep representations equivariant to rotation. However, feature maps should be copied and rotated four times in each layer in their approach, which causes much running time and memory overhead. In order to address this problem, we propose Deep Rotation Equivariant Network consisting of cycle layers, isotonic layers and decycle layers. Our proposed layers apply rotation transformation on filters rather than feature maps, achieving a speed up of more than 2 times with even less memory overhead. We evaluate DRENs on Rotated MNIST and CIFAR-10 datasets and demonstrate that it can improve the performance of state-of-the-art architectures.


title: VANETs Meet Autonomous Vehicles: A Multimodal 3D Environment Learning Approach

abstract: In this paper, we design a multimodal framework for object detection, recognition and mapping based on the fusion of stereo camera frames, point cloud Velodyne Lidar scans, and Vehicle-to-Vehicle (V2V) Basic Safety Messages (BSMs) exchanged using Dedicated Short Range Communication (DSRC). We merge the key features of rich texture descriptions of objects from 2D images, depth and distance between objects provided by 3D point cloud and awareness of hidden vehicles from BSMs' 3D information. We present a joint pixel to point cloud and pixel to V2V correspondences of objects in frames from the Kitti Vision Benchmark Suite by using a semi-supervised manifold alignment approach to achieve camera-Lidar and camera-V2V mapping of their recognized objects that have the same underlying manifold.


title: Quantum thermodynamic cycle with quantum phase transition

abstract: With the Lipkin-Meshkov-Glick (LMG) model as an illustration, we construct a thermodynamic cycle composed of two isothermal processes and two isomagnetic field processes and study the thermodynamic performance of this cycle accompanied by the quantum phase transition (QPT). We find that for a finite particle system working below the critical temperature, the efficiency of the cycle is capable of approaching the Carnot limit when the external magnetic field \lambda_{1} corresponding to one of the isomagnetic processes reaches the crosspoint of the ground states' energy level, which can become critical point of the QPT in large N limit. Our analysis proves that the system's energy level crossings at low temperature limits can lead to significant efficiency improvement of the quantum heat engine. In the case of the thermodynamics limit, analytical partition function is obtained to study the efficiency of the cycle at high and low temperature limits. At low temperatures, when the magnetic fields of the isothermal processes are located on both sides of the critical point of the QPT, the cycle obtains the maximum efficiency and the Carnot efficiency can be achieved. This observation demonstrate that the QPT of the LMG model below critical temperature is beneficial to the thermodynamic cycle's operation.


title: Dedekind sums take each value infinitely many times

abstract: For $a\in \Bbb Z$ and $b\in\Bbb N$, $(a,b)=1$, let $s(a,b)$ denote the classical Dedekind sum. We show that Dedekind sums take this value infinitely many times in the following sense. There are pairs $(a_i,b_i)$, $i\in\Bbb N$, with $b_i$ tending to infinity as $i$ grows, such that $s(a_i,b_i)=s(a,b)$ for all $i\in \Bbb N$.


title: On Using Time Without Clocks via Zigzag Causality

abstract: Even in the absence of clocks, time bounds on the duration of actions enable the use of time for distributed coordination. This paper initiates an investigation of coordination in such a setting. A new communication structure called a zigzag pattern is introduced, and shown to guarantee bounds on the relative timing of events in this clockless model. Indeed, zigzag patterns are shown to be necessary and sufficient for establishing that events occur in a manner that satisfies prescribed bounds. We capture when a process can know that an appropriate zigzag pattern exists, and use this to provide necessary and sufficient conditions for timed coordination of events using a full-information protocol in the clockless model.


title: On the effect of polydispersity and rotation on the Brinkman force induced by a cloud of particles on a viscous incompressible flow

abstract: In this paper, we are interested in the collective friction of a cloud of particles on the viscous incompressible fluid in which they are moving. The particles velocities are assumed to be given and the fluid is assumed to be driven by the stationary Stokes equations. We consider the limit where the number N of particles goes to infinity with their diameters of order 1/N and their mutual distances of order (1/N)^{1/3}. The rigorous convergence of the fluid velocity to a limit which is solution to a stationary Stokes equation set in the full space but with an extra term, referred to as the Brinkman force, was proven by Desvillettes, Golse and Ricci when the particles are identical spheres in prescribed translations. Our result here is an extension to particles of arbitrary shapes in prescribed translations and rotations. The limit Stokes-Brinkman system involves the particle distribution in position, velocity and shape, through the so-called Stokes' resistance matrices.


title: Technical Design Report (TDR): Searching for a Sterile Neutrino at J-PARC MLF (E56, JSNS2)

abstract: In this document, the technical details of the JSNS$^2$ (J-PARC Sterile Neutrino Search at J-PARC Spallation Neutron Source) experiment are described. The search for sterile neutrinos is currently one of the hottest topics in neutrino physics. The JSNS$^2$ experiment aims to search for the existence of neutrino oscillations with $\Delta m^2$ near 1 eV$^2$ at the J-PARC Materials and Life Science Experimental Facility (MLF). A 1 MW beam of 3 GeV protons incident on a spallation neutron target produces an intense neutrino beam from muon decay at rest. Neutrinos come predominantly from $\mu^+$ decay: $\mu^{+} \to e^{+} + \bar{\nu}_{\mu} + \nu_{e}$. The experiment will search for $\bar{\nu}_{\mu}$ to $\bar{\nu}_{e}$ oscillations which are detected by the inverse beta decay interaction $\bar{\nu}_{e} + p \to e^{+} + n$, followed by gammas from neutron capture on Gd. The detector has a fiducial volume of 17 tons and is located 24 meters away from the mercury target. JSNS$^2$ offers the ultimate direct test of the LSND anomaly. In addition to the sterile neutrino search, the physics program includes cross section measurements with neutrinos with a few 10's of MeV from muon decay at rest and with monochromatic 236 MeV neutrinos from kaon decay at rest. These cross sections are relevant for our understanding of supernova explosions and nuclear physics.


title: Geometric location of periodic points of 2-ramified power series

abstract: In this paper we study the geometric location of periodic points of power series defined over fields of prime characteristic $p$. More specifically, we find a lower bound for the absolute value of all periodic points in the open unit disk of minimal period $p^n$ of 2-ramified power series. We prove that this bound is optimal for a large class of power series. Our main technical result is a computation of the first significant terms of the $p^n$th iterate of 2-ramified power series. As a by-product we obtain a self-contained proof of the characterization of 2-ramified power series.


title: Self-supervised learning of visual features through embedding images into text topic spaces

abstract: End-to-end training from scratch of current deep architectures for new computer vision problems would require Imagenet-scale datasets, and this is not always possible. In this paper we present a method that is able to take advantage of freely available multi-modal content to train computer vision algorithms without human supervision. We put forward the idea of performing self-supervised learning of visual features by mining a large scale corpus of multi-modal (text and image) documents. We show that discriminative visual features can be learnt efficiently by training a CNN to predict the semantic context in which a particular image is more probable to appear as an illustration. For this we leverage the hidden semantic structures discovered in the text corpus with a well-known topic modeling technique. Our experiments demonstrate state of the art performance in image classification, object detection, and multi-modal retrieval compared to recent self-supervised or natural-supervised approaches.


title: Faithful (meta-)encodings of programmable strategies into term rewriting systems

abstract: Rewriting is a formalism widely used in computer science and mathematical logic. When using rewriting as a programming or modeling paradigm, the rewrite rules describe the transformations one wants to operate and rewriting strategies are used to con- trol their application. The operational semantics of these strategies are generally accepted and approaches for analyzing the termination of specific strategies have been studied. We propose in this paper a generic encoding of classic control and traversal strategies used in rewrite based languages such as Maude, Stratego and Tom into a plain term rewriting system. The encoding is proven sound and complete and, as a direct consequence, estab- lished termination methods used for term rewriting systems can be applied to analyze the termination of strategy controlled term rewriting systems. We show that the encoding of strategies into term rewriting systems can be easily adapted to handle many-sorted signa- tures and we use a meta-level representation of terms to reduce the size of the encodings. The corresponding implementation in Tom generates term rewriting systems compatible with the syntax of termination tools such as AProVE and TTT2, tools which turned out to be very effective in (dis)proving the termination of the generated term rewriting systems. The approach can also be seen as a generic strategy compiler which can be integrated into languages providing pattern matching primitives; experiments in Tom show that applying our encoding leads to performances comparable to the native Tom strategies.


title: Discovery of a new branch of the Taurid meteoroid stream as a real source of potentially hazardous bodies

abstract: Taurid meteor shower produces prolonged but usually low activity every October and November. In some years, however, the activity is significantly enhanced. Previous studies based on long-term activity statistics concluded that the enhancement is caused by a swarm of meteoroids locked in 7:2 resonance with Jupiter. Here we present precise data on 144 Taurid fireballs observed by new digital cameras of the European Fireball Network in the enhanced activity year 2015. Orbits of 113 fireballs show common characteristics and form together a well defined orbital structure, which we call new branch. We found that this branch is characterized by longitudes of perihelia lying between 155.9-160o and latitudes of perihelia between 4.2-5.7o. Semimajor axes are between 2.23-2.28 AU and indeed overlap with the 7:2 resonance. Eccentricities are in wide range 0.80-0.90. The orbits form a concentric ring in the inner solar system. The masses of the observed meteoroids were in a wide range from 0.1 g to more than 1000 kg. We found that all meteoroids larger than 300 g were very fragile, while those smaller than 30 g were much more compact. Based on orbital characteristics, we argue that asteroids 2015 TX24 and 2005 UR, both of diameters 200-300 meters, are direct members of the new branch. It is therefore very likely that the new branch contains also numerous still not discovered objects of decameter or even larger size. Since asteroids of sizes of tens to hundreds meters pose a treat to the ground even if they are intrinsically weak, impact hazard increases significantly when the Earth encounters the Taurid new branch every few years. Further studies leading to better description of this real source of potentially hazardous objects, which can be large enough to cause significant regional or even continental damage on the Earth, are therefore extremely important.


title: A $\mathcal C^{2,\alpha}$ estimate of the complex Monge-Amp\`ere equation

abstract: In this paper, we prove a $\mathcal C^{2,\alpha}$-estimate for the solution to the complex Monge-Amp\`ere equation $\det(u_{i\bar{j}})=f$ with $0< f\in \mathcal C^{\alpha}$, under the assumption that $u\in \mathcal C^{1,\beta }$ for some $\beta <1$ which depends on $n$ and $\alpha $.


title: Optical unidirectional amplification in a three-mode optomechanical system

abstract: We study the directional amplification of an optical probe field in a three-mode optomechanical system, where the mechanical resonator interacts with two linearly-coupled optical cavities and the cavities are driven by strong optical pump fields. The optical probe field is injected into one of the cavity modes, and at the same time, the mechanical resonator is subject to a mechanical drive with the driving frequency equal to the frequency difference between the optical probe and pump fields. We show that the transmission of the probe field can be amplified in one direction and de-amplified in the opposite direction. This directional amplification or de-amplification results from the constructive or destruction interference between different transmission paths in this three-mode optomechanical system.


title: Spin-singlet superconductivity in doped topological crystalline insulator Sn_{0.96}In_{0.04}Te

abstract: The In-doped topological crystalline insulator Sn_{1-x}In_xTe is a candidate for a topological superconductor, where pseudo-spin-triplet state has been proposed. To clarify the spin symmetry of Sn_{1-x}In_xTe, we perform ^{125}Te-nuclear magnetic resonance (NMR) measurements in polycrystalline samples with 0< x <0.15. The penetration depth calculated from the NMR line width is T-independent below half the superconducting transition temperature (T_c) in polycrystalline Sn_{0.96}In_{0.04}Te, which indicates a fully-opened superconducting gap. In this sample, the spin susceptibility measured by the spin Knight shift (K_s) at an external magnetic field of H_0 = 0.0872 T decreases below T_c, and K_s(T=0)/K_s(T=T_c) reaches to 0.36 \pm 0.10, which is far below the limiting value 2/3 expected for a spin-triplet state for a cubic crystal structure. Our result indicates that polycrystalline Sn_{0.96}In_{0.04}Te is a spin-singlet superconductor.


title: Bright 5 - 85 MeV Compton gamma-ray pulses from GeV laser-plasma accelerator and plasma mirror

abstract: We convert a GeV laser-plasma electron accelerator into a compact femtosecond-pulsed $\gamma$-ray source by inserting a $100 \mu$m-thick glass plate $\sim3$ cm after the accelerator exit. With near-unity reliability, and requiring only crude alignment, this glass plasma mirror retro-reflected spent drive laser pulses (photon energy $\hbar\omega_L = 1.17$ eV) with $>50\%$ efficiency back onto trailing electrons (peak Lorentz factor $1000 < \gamma_e < 4400$), creating an optical undulator that generated $\sim10^8 \gamma$-ray photons with sub-mrad divergence, estimated peak brilliance $\sim10^{21}$ photons/s/mm$^2$/mrad$^2$/$0.1\%$ bandwidth and negligible bremsstrahlung background. The $\gamma$-ray photon energy $E_\gamma = 4\gamma_e^2 \hbar\omega_L$, inferred from the measured $\gamma_e$ on each shot, peaked from 5 to 85 MeV, spanning a range otherwise available with comparable brilliance only from large-scale GeV-linac-based high-intensity $\gamma$-ray sources.


title: Effect of partially-screened nuclei on fast-electron dynamics

abstract: We analyze the dynamics of fast electrons in plasmas containing partially ionized impurity atoms, where the screening effect of bound electrons must be included. We derive analytical expressions for the deflection and slowing-down frequencies, and show that they are increased significantly compared to the results obtained with complete screening, already at sub-relativistic electron energies. Furthermore, we show that the modifications to the deflection and slowing down frequencies are of equal importance in describing the runaway current evolution. Our results greatly affect fast-electron dynamics and have important implications, e.g. for the efficacy of mitigation strategies for runaway electrons in tokamak devices, and energy loss during relativistic breakdown in atmospheric discharges.


title: Fast-Slow Recurrent Neural Networks

abstract: Processing sequential data of variable length is a major challenge in a wide range of applications, such as speech recognition, language modeling, generative image modeling and machine translation. Here, we address this challenge by proposing a novel recurrent neural network (RNN) architecture, the Fast-Slow RNN (FS-RNN). The FS-RNN incorporates the strengths of both multiscale RNNs and deep transition RNNs as it processes sequential data on different timescales and learns complex transition functions from one time step to the next. We evaluate the FS-RNN on two character level language modeling data sets, Penn Treebank and Hutter Prize Wikipedia, where we improve state of the art results to $1.19$ and $1.25$ bits-per-character (BPC), respectively. In addition, an ensemble of two FS-RNNs achieves $1.20$ BPC on Hutter Prize Wikipedia outperforming the best known compression algorithm with respect to the BPC measure. We also present an empirical investigation of the learning and network dynamics of the FS-RNN, which explains the improved performance compared to other RNN architectures. Our approach is general as any kind of RNN cell is a possible building block for the FS-RNN architecture, and thus can be flexibly applied to different tasks.


title: Representing the suffix tree with the CDAWG

abstract: Given a string $T$, it is known that its suffix tree can be represented using the compact directed acyclic word graph (CDAWG) with $e_T$ arcs, taking overall $O(e_T+e_{{\overline{T}}})$ words of space, where ${\overline{T}}$ is the reverse of $T$, and supporting some key operations in time between $O(1)$ and $O(\log{\log{n}})$ in the worst case. This representation is especially appealing for highly repetitive strings, like collections of similar genomes or of version-controlled documents, in which $e_T$ grows sublinearly in the length of $T$ in practice. In this paper we augment such representation, supporting a number of additional queries in worst-case time between $O(1)$ and $O(\log{n})$ in the RAM model, without increasing space complexity asymptotically. Our technique, based on a heavy path decomposition of the suffix tree, enables also a representation of the suffix array, of the inverse suffix array, and of $T$ itself, that takes $O(e_T)$ words of space, and that supports random access in $O(\log{n})$ time. Furthermore, we establish a connection between the reversed CDAWG of $T$ and a context-free grammar that produces $T$ and only $T$, which might have independent interest.


title: Crust-core transition of a neutron star: effects of the symmetry energy and temperature under strong magnetic fields

abstract: We study the simultaneous effects of the symmetry energy and temperature on the crust-core transition of a magnetar. The dynamical and the thermodynamical spinodals are used to calculate the transition region within a relativistic mean-field approach for the equation of state. Quantizing magnetic fields with intensities in the range of $ 2\times 10^{15}<B<5\times 10 ^{16}$G are considered. Under these strong magnetic fields, the crust extension is very sensitive to the density dependence of the symmetry energy, and the properties that depend on the crust thickness could set a constraint on the equation of state. It is shown that the effect on the extension of the crust-core transition is washed out for temperatures above $10^{9}$ K. However, for temperatures below that value, a noticeable effect exists that grows as the temperature decreases and which should be taken into account when the evolution of magnetars is studied.


title: The supernova-regulated ISM. III. Generation of vorticity, helicity and mean flows

abstract: The forcing of interstellar turbulence, driven mainly by supernova explosions, is irrotational in nature, but the development of significant amounts of vorticity and helicity, accompanied by large-scale dynamo action, has been reported. Several earlier investigations examined vorticity production in simpler systems; here all the relevant processes can be considered simultaneously. We also investigate the mechanisms for the generation of net helicity and large-scale flow in the system. We use a three-dimensional, stratified, rotating and shearing local simulation domain of the size 1x1x2 kpc$^3$, forced with SN explosions occurring at the rate typical of the solar neighbourhood in the Milky Way. In addition to the nominal simulation run with realistic Milky Way parameters, we vary the rotation and shear rates, but keep the absolute value of their ratio fixed. Reversing the sign of shear vs. rotation allows us to separate the rotation- and shear-generated contributions. As in earlier studies, we find the generation of significant amounts of vorticity, with on average 65% of the kinetic energy being in the rotational modes. The vorticity production can be related to the baroclinicity of the flow, especially in the regions of hot, dilute clustered supernova bubbles. In these regions, the vortex stretching acts as a sink of vorticity. The net helicities produced by rotation and shear are of opposite signs for physically motivated rotation laws, with the solar neighbourhood parameters resulting in the near cancellation of the total net helicity. We also find the excitation of oscillatory mean flows, the strength and oscillation period of which depend on the Coriolis and shear parameters; we interpret these as signatures of the anisotropic kinetic (AKA) effect. We use the method of moments to fit for the turbulent transport coeffcients, and find $\alpha_{\rm AKA}$ values of the order 3-5 km/s.


title: Higher order Cheeger inequalities for Steklov eigenvalues

abstract: We prove a lower bound for the $k$-th Steklov eigenvalues in terms of an isoperimetric constant called the $k$-th Cheeger-Steklov constant in three different situations: finite spaces, measurable spaces, and Riemannian manifolds. These lower bounds can be considered as higher order Cheeger type inequalities for the Steklov eigenvalues. In particular it extends the Cheeger type inequality for the first nonzero Steklov eigenvalue previously studied by Escobar in 1997 and by Jammes in 2015 to higher order Steklov eigenvalues. The technique we develop to get this lower bound is based on considering a family of accelerated Markov operators in the finite and mesurable situations and of mass concentration deformations of the Laplace-Beltrami operator in the manifold setting which converges uniformly to the Steklov operator. As an intermediary step in the proof of the higher order Cheeger type inequality, we define the Dirichlet-Steklov connectivity spectrum and show that the Dirichlet connectivity spectra of this family of operators converges to (or bounded by) the Dirichlet-Steklov spectrum uniformly. Moreover, we obtain bounds for the Steklov eigenvalues in terms of its Dirichlet-Steklov connectivity spectrum which is interesting in its own right and is more robust than the higher order Cheeger type inequalities. The Dirichlet-Steklov spectrum is closely related to the Cheeger-Steklov constants.


title: Asymptotic Lipschitz regularity of viscosity solutions of Hamilton-Jacobi equations

abstract: For each continuous initial data $\varphi(x)\in C(M,\mathbb{R})$, we obtain the asymptotic Lipschitz regularity of the viscosity solution of the following evolutionary Hamilton-Jacobi equation with convex and coercive Hamiltonians: \partial_tu(x,t)+H(x,\partial_xu(x,t))=0, u(x,0)=\varphi(x).


title: Phononic Josephson oscillation and self-trapping with two-phonon exchange interaction

abstract: We propose a bosonic Josephson junction (BJJ) in two nonlinear mechanical resonator coupled through two-phonon exchange interaction induced by quadratic optomechanical couplings. The nonlinear dynamic equations and effective Hamiltonian are derived to describe behaviors of the BJJ. We show that the BJJ can work in two different dynamical regimes: Josephson oscillation and macroscopic self-trapping. The system can transfer from one regime to the other one when the self-interaction and asymmetric parameters exceed their critical values. We predict that a transition from Josephson oscillation to macroscopic self-trapping can be induced by the phonon damping in the asymmetric BJJs. Our results opens up a way to demonstrate BJJ with two-phonon exchange interaction and can be applied to other systems, such as the optical and microwave systems.


title: Bloch oscillations sustained by nonlinearity

abstract: We demonstrate that nonlinearity may play a constructive role in supporting Bloch oscillations in a model which is discrete, in one dimension and continuous in the orthogonal one. The model can be experimentally realized in several fields of physics such as optics and Bose-Einstein condensates. We demonstrate that designing an optimal relation between the nonlinearity and the linear gradient strength provides extremely long-lived Bloch oscillations with little degradation. Such robust oscillations can be observed for a broad range of parameters and even for moderate nonlinearities and large enough values of linear potential. We also present an approximate analytical description of the wave packet's evolution featuring a hybrid Bloch oscillating wave-soliton behavior that excellently corresponds to the direct numerical simulations.


title: Magnetopumping current in graphene Corbino pump

abstract: We study conductance and adiabatic pumped charge and spin currents in a graphene quantum pump with Corbino geometry in the presence of an applied perpendicular magnetic field. The pump is driven by the periodic and out of phase modulations of the magnetic field and an electrostatic potential applied to the ring area of the pump. We show that the Zeeman splitting, despite of its smallness, can suppress the conductance oscillations at the zero doping and in a threshold value for the flux piercing the ring area which depends on the inner lead radius and thus on the flux penetrating in it. Moreover, it generates a considerable spin conductance at infinitesimal nonzero doping and at the magnetic flux, that charge conductance starts to suppress. We find that the pumped charge and spin currents increase by the magnetic field with small oscillations until they start to suppress due to the effect of the nonzero doping and the Zeeman splitting. In graphene Corbino pumps with small inner leads the Zeeman splitting shows its effect in a large value of the magnetic field and thus we can get a considerable pumped charge and spin currents at the enough small magnetic fields.


title: Long-lived mesoscopic entanglement between two damped infinite harmonic chains

abstract: We consider two chains, each made of $N$ independent oscillators, immersed in a common thermal bath and study the dynamics of their mutual quantum correlations in the thermodynamic, large-$N$ limit. We show that dissipation and noise due to the presence of the external environment are able to generate collective quantum correlations between the two chains at the mesoscopic level. The created collective quantum entanglement between the two many-body systems turns out to be rather robust, surviving for asymptotically long times even for non vanishing bath temperatures.


title: Maximal quantum Fisher information matrix

abstract: We study the existence of the maximal quantum Fisher information matrix in multi-parameter quantum estimation, which bounds the ultimate precision limit. We show that when the maximal quantum Fisher information matrix exists, it can be directly obtained from the underlying dynamics. Examples are then provided to demonstrate the usefulness of the maximal quantum Fisher information matrix by deriving various tradeoff relations in multi-parameter quantum estimation and obtaining the bounds for the scalings of the precision limit.


title: Experimental generalized quantum suppression law in Sylvester interferometers

abstract: Photonic interference is a key quantum resource for optical quantum computation, and in particular for so-called boson sampling machines. In interferometers with certain symmetries, genuine multiphoton quantum interference effectively suppresses certain sets of events, as in the original Hong-Ou-Mandel effect. Recently, it was shown that some classical and semi-classical models could be ruled out by identifying such suppressions in Fourier interferometers. Here we propose a suppression law suitable for random-input experiments in multimode Sylvester interferometers, and verify it experimentally using 4- and 8-mode integrated interferometers. The observed suppression is stronger than what is observed in Fourier interferometers of the same size, and could be relevant to certification of boson sampling machines and other experiments relying on bosonic interference.


title: Coverings of Spectral Triples

abstract: It is well-known that any covering space of a Riemannian manifold has the natural structure of a Riemannian manifold. This article contains a noncommutative generalization of this fact. Since any Riemannian manifold with a Spin-structure defines a spectral triple, the spectral triple can be regarded as a noncommutative Spin-manifold. Similarly there is an algebraic construction which is a noncommutative generalization of topological covering. This article contains a construction of spectral triple on the "noncommutative covering space".


title: On the Success Probability of Decoding (Partial) Unit Memory Codes

abstract: In this paper, we derive analytic expressions for the success probability of decoding (Partial) Unit Memory codes in memoryless channels. An applications of this result is that these codes outperform individual block codes in certain channels.


title: Quantitative characterization of pore structure of several biochars with 3D imaging

abstract: Pore space characteristics of biochars may vary depending on the used raw material and processing technology. Pore structure has significant effects on the water retention properties of biochar amended soils. In this work, several biochars were characterized with three-dimensional imaging and image analysis. X-ray computed microtomography was used to image biochars at resolution of 1.14 $\mu$m and the obtained images were analysed for porosity, pore-size distribution, specific surface area and structural anisotropy. In addition, random walk simulations were used to relate structural anisotropy to diffusive transport. Image analysis showed that considerable part of the biochar volume consist of pores in size range relevant to hydrological processes and storage of plant available water. Porosity and pore-size distribution were found to depend on the biochar type and the structural anisotopy analysis showed that used raw material considerably affects the pore characteristics at micrometre scale. Therefore attention should be paid to raw material selection and quality in applications requiring optimized pore structure.


title: PET-MRI Joint Reconstruction by Joint Sparsity Based Tight Frame Regularization

abstract: Recent technical advances lead to the coupling of PET and MRI scanners, enabling to acquire functional and anatomical data simultaneously. In this paper, we propose a tight frame based PET-MRI joint reconstruction model via the joint sparsity of tight frame coefficients. In addition, a non-convex balanced approach is adopted to take the different regularities of PET and MRI images into account. To solve the nonconvex and nonsmooth model, a proximal alternating minimization algorithm is proposed, and the global convergence is present based on Kurdyka-Lojasiewicz property. Finally, the numerical experiments show that the our proposed models achieve better performance over the existing PET-MRI joint reconstruction models.


title: Restriction of Odd Degree Characters of $\mathfrak{S}_n$

abstract: Let $n$ and $k$ be natural numbers such that $2^k < n$. We study the restriction to $\mathfrak{S}_{n-2^k}$ of odd-degree irreducible characters of the symmetric group $\mathfrak{S}_n$. This analysis completes the study begun in [Ayyer A., Prasad A., Spallone S., Sem. Lothar. Combin. 75 (2015), Art. B75g, 13 pages] and recently developed in [Isaacs I.M., Navarro G., Olsson J.B., Tiep P.H., J. Algebra 478 (2017), 271-282].


title: Efficient Covariance Approximations for Large Sparse Precision Matrices

abstract: The use of sparse precision (inverse covariance) matrices has become popular because they allow for efficient algorithms for joint inference in high-dimensional models. Many applications require the computation of certain elements of the covariance matrix, such as the marginal variances, which may be non-trivial to obtain when the dimension is large. This paper introduces a fast Rao-Blackwellized Monte Carlo sampling based method for efficiently approximating selected elements of the covariance matrix. The variance and confidence bounds of the approximations can be precisely estimated without additional computational costs. Furthermore, a method that iterates over subdomains is introduced, and is shown to additionally reduce the approximation errors to practically negligible levels in an application on functional magnetic resonance imaging data. Both methods have low memory requirements, which is typically the bottleneck for competing direct methods.


title: Combinatorial n-fold Integer Programming and Applications

abstract: Many fundamental NP-hard problems can be formulated as integer linear programs (ILPs). A famous algorithm by Lenstra solves ILPs in time that is exponential only in the dimension of the program, and polynomial in the size of the ILP. That algorithm became a ubiquitous tool in the design of fixed-parameter algorithms for NP-hard problems, where one wishes to isolate the hardness of a problemby some parameter. However, in many cases using Lenstra's algorithm has two drawbacks: First, the run time of the resulting algorithms is often doubly-exponential in the parameter, and second, an ILP formulation in small dimension cannot easily express problems involving many different costs. Inspired by the work of Hemmecke, Onn and Romanchuk [Math. Prog. 2013], we develop a single-exponential algorithm for so-called combinatorial n-fold integer programs, which are remarkably similar to prior ILP formulations for various problems, but unlike them, also allow variable dimension. We then apply our algorithm to a few representative problems like Closest String, Swap Bribery, Weighted Set Multicover, and obtain exponential speedups in the dependence on the respective parameters, the input size, or both. Unlike Lenstra's algorithm, which is essentially a bounded search tree algorithm, our result uses the technique of augmenting steps. At its heart is a deep result stating that in combinatorial n-fold IPs, existence of an augmenting step implies existence of a \local" augmenting step, which can be found using dynamic programming. Our results provide an important insight into many problems by showing that they exhibit this phenomenon, and highlights the importance of augmentation techniques.


title: Invariance entropy, quasi-stationary measures and control sets

abstract: For control systems in discrete time, this paper discusses measure-theoretic invariance entropy for a subset Q of the state space with respect to a quasi-stationary measure obtained by endowing the control range with a probability measure. The main results show that this entropy is invariant under measurable transformations and that it is already determined by certain subsets of Q which are characterized by controllability properties.


title: Josephson Current in Ballistic Graphene Corbino Disk

abstract: We solve Dirac-Bogoliubov-De-Gennes (DBdG) equation in a superconductor-normal graphene superconductor (SGS) junction with Corbino disk structure to investigate the Josephson current through this junction. We find that the critical current $I_c$ has a nonzero value at Dirac point in which the concentration of the carriers is zero. We show this nonzero critical current depends on the system geometry and it decreases monotonically to zero by increasing the ratio of the outer to inner radii of the Corbino disk ($R_2/R_1$), while in the limit of $R_2/R_1 \rightarrow 1$ it scales like a diffusive Corbino disk. The product of the critical current and the normal-state resistance $I_cR_N$ attains the same value for the planar structure at zero doping. These results reveals the pseudodiffusive behavior of the graphene Corbino Josephson junction similar to the planar structure.


title: Matrix of Polynomials Model based Polynomial Dictionary Learning Method for Acoustic Impulse Response Modeling

abstract: We study the problem of dictionary learning for signals that can be represented as polynomials or polynomial matrices, such as convolutive signals with time delays or acoustic impulse responses. Recently, we developed a method for polynomial dictionary learning based on the fact that a polynomial matrix can be expressed as a polynomial with matrix coefficients, where the coefficient of the polynomial at each time lag is a scalar matrix. However, a polynomial matrix can be also equally represented as a matrix with polynomial elements. In this paper, we develop an alternative method for learning a polynomial dictionary and a sparse representation method for polynomial signal reconstruction based on this model. The proposed methods can be used directly to operate on the polynomial matrix without having to access its coefficients matrices. We demonstrate the performance of the proposed method for acoustic impulse response modeling.


title: Robot Introspection with Bayesian Nonparametric Vector Autoregressive Hidden Markov Models

abstract: Robot introspection, as opposed to anomaly detection typical in process monitoring, helps a robot understand what it is doing at all times. A robot should be able to identify its actions not only when failure or novelty occurs, but also as it executes any number of sub-tasks. As robots continue their quest of functioning in unstructured environments, it is imperative they understand what is it that they are actually doing to render them more robust. This work investigates the modeling ability of Bayesian nonparametric techniques on Markov Switching Process to learn complex dynamics typical in robot contact tasks. We study whether the Markov switching process, together with Bayesian priors can outperform the modeling ability of its counterparts: an HMM with Bayesian priors and without. The work was tested in a snap assembly task characterized by high elastic forces. The task consists of an insertion subtask with very complex dynamics. Our approach showed a stronger ability to generalize and was able to better model the subtask with complex dynamics in a computationally efficient way. The modeling technique is also used to learn a growing library of robot skills, one that when integrated with low-level control allows for robot online decision making.


title: Rapidly star-forming galaxies adjacent to quasars at redshifts exceeding 6

abstract: The existence of massive ($10^{11}$ solar masses) elliptical galaxies by redshift z~4 (when the Universe was 1.5 billion years old) necessitates the presence of galaxies with star-formation rates exceeding 100 solar masses per year at z>6 (corresponding to an age of the Universe of less than 1 billion years). Surveys have discovered hundreds of galaxies at these early cosmic epochs, but their star-formation rates are more than an order of magnitude lower. The only known galaxies with very high star-formation rates at z>6 are, with only one exception, the host galaxies of quasars, but these galaxies also host accreting supermassive (more than $10^9$ solar masses) black holes, which probably affect the properties of the galaxies. Here we report observations of an emission line of singly ionized carbon ([CII] at a wavelength of 158 micrometres) in four galaxies at z>6 that are companions of quasars, with velocity offsets of less than 600 kilometers per second and linear offsets of less than 600 kiloparsecs. The discovery of these four galaxies was serendipitous; they are close to their companion quasars and appear bright in the far-infrared. On the basis of the [CII] measurements, we estimate star-formation rates in the companions of more than 100 solar masses per year. These sources are similar to the host galaxies of the quasars in [CII] brightness, linewidth and implied dynamical masses, but do not show evidence for accreting supermassive black holes. Similar systems have previously been found at lower redshift. We find such close companions in four out of twenty-five z>6 quasars surveyed, a fraction that needs to be accounted for in simulations. If they are representative of the bright end of the [CII] luminosity function, then they can account for the population of massive elliptical galaxies at z~4 in terms of cosmic space density.


title: Improving institutional memory on challenges and methods for estimation of pig herd antimicrobial exposure based on data from the Danish Veterinary Medicines Statistics Program (VetStat)

abstract: With the increasing occurrence of antimicrobial resistance, more attention has been directed towards surveillance of both human and veterinary antimicrobial use. Since the early 2000s, several research papers on Danish pig antimicrobial usage have been published, based on data from the Danish Veterinary Medicines Statistics Program (VetStat). VetStat was established in 2000, as a national database containing detailed information on purchases of veterinary medicine. This paper presents a critical set of challenges originating from static system features, which researchers must address when estimating antimicrobial exposure in Danish pig herds. Most challenges presented are followed by at least one robust solution. A set of challenges requiring awareness from the researcher, but for which no immediate solution was available, were also presented. The selection of challenges and solutions was based on a consensus by a cross-institutional group of researchers working in projects using VetStat data. No quantitative data quality evaluations were performed, as the frequency of errors and inconsistencies in a dataset will vary, depending on the period covered in the data. Instead, this paper focuses on clarifying how VetStat data may be translated to an estimation of the antimicrobial exposure at herd level, by suggesting uniform methods of extracting and editing data, in order to obtain reliable and comparable estimates on pig antimicrobial consumption for research purposes.


title: Towards Understanding the Invertibility of Convolutional Neural Networks

abstract: Several recent works have empirically observed that Convolutional Neural Nets (CNNs) are (approximately) invertible. To understand this approximate invertibility phenomenon and how to leverage it more effectively, we focus on a theoretical explanation and develop a mathematical model of sparse signal recovery that is consistent with CNNs with random weights. We give an exact connection to a particular model of model-based compressive sensing (and its recovery algorithms) and random-weight CNNs. We show empirically that several learned networks are consistent with our mathematical analysis and then demonstrate that with such a simple theoretical framework, we can obtain reasonable re- construction results on real images. We also discuss gaps between our model assumptions and the CNN trained for classification in practical scenarios.


title: Bayesian Compression for Deep Learning

abstract: Compression and computational efficiency in deep learning have become a problem of great significance. In this work, we argue that the most principled and effective way to attack this problem is by adopting a Bayesian point of view, where through sparsity inducing priors we prune large parts of the network. We introduce two novelties in this paper: 1) we use hierarchical priors to prune nodes instead of individual weights, and 2) we use the posterior uncertainties to determine the optimal fixed point precision to encode the weights. Both factors significantly contribute to achieving the state of the art in terms of compression rates, while still staying competitive with methods designed to optimize for speed or energy efficiency.


title: Towards Near Real-Time BGP Deep Analysis: A Big-Data Approach

abstract: BGP (Border Gateway Protocol) serves as the primary routing protocol for the Internet, enabling Autonomous Systems (individual network operators) to exchange network reachability information. Alongside significant on-going research and development efforts, there is a practical need to understand the nature of events that occur on the Internet. Network operators are acutely aware of security-related incidents such as 'Prefix Hijacking' as well as the impact of network instabilities that ripple through the Internet. Recent research focused on the study of BGP anomalies (both network/prefix instability and security-related incidents) has been based on the analysis of historical logs. Further analysis to understand the nature of these anomalous events is not always sufficient to be able to differentiate malicious activities, such as prefix- or sub-prefix- hijacking, from those events caused by inadvertent misconfigurations. In addition, such techniques are challenged by a lack of sufficient resources to store and process data feeds in real-time from multiple BGP Vantage Points (VPs). In this paper, we present a BGP Deep-analysis application developed using the PNDA (Platform for Network Data Analytics) 'Big-Data' platform. PNDA provides a highly scalable environment that enables the ingestion and processing of 'live' BGP feeds from many vantage points in a schema-agnostic manner. The Apache Spark-based application, in conjunction with PNDA's distributed processing capabilities, is able to perform high-level insights as well as near-to-real-time statistical analysis


title: Packing parameters in graphs: New bounds and a solution to an open problem

abstract: In this paper, we investigate the packing parameters in graphs. By applying the Mantel's theorem, We give upper bounds on packing and open packing numbers of triangle-free graphs along with characterizing the graphs for which the equalities hold and exhibit sharp Nordhaus-Gaddum type inequalities for packing numbers. We also solve the open problem of characterizing all connected graphs with $\rho_{o}(G)=n-\omega(G)$ posed in [S. Hamid and S. Saravanakumar, {\em Packing parameters in graphs}, Discuss Math. Graph Theory, 35 (2015), 5--16].


title: Alliance formation with exclusion in the spatial public goods game

abstract: Detecting defection and alarming partners about the possible danger could be essential to avoid being exploited. This act, however, may require a huge individual effort from those who take this job, hence such a strategy seems to be unfavorable. But structured populations can provide an opportunity where a largely unselfish excluder strategy can form an effective alliance with other cooperative strategies, hence they can sweep out defection. Interestingly, this alliance is functioning even at the extremely high cost of exclusion where the sole application of an exclusion strategy would be harmful otherwise. These results may explain why the emergence of extreme selfless behavior is not necessarily against individual selection but could be the result of an evolutionary process.

categories: physics.soc-ph cs.GT nlin.CG physics.bio-ph