title: Hamiltonians for the zeros of a general family of zeta functions

abstract: Towards the Hilbert-P\'olya conjecture, in this paper, we present a general construction of Hamiltonian $\hat{H}_{f},$ which leads a general family of Hurwitz zeta functions $(-1)^{z_{n}-1}L(f,z_{n},x+1)$ defined by Mellin transform becomes their eigenstates under a suitable boundary condition, and the eigenvalues $E_{n}$ have the property that $z_{n}=\frac{1}{2}(1-iE_{n})$ are the zeros of a general family of zeta functions $L(f,z)\equiv L(f,z,1)$.


title: Suspicious News Detection Using Micro Blog Text

abstract: We present a new task, suspicious news detection using micro blog text. This task aims to support human experts to detect suspicious news articles to be verified, which is costly but a crucial step before verifying the truthfulness of the articles. Specifically, in this task, given a set of posts on SNS referring to a news article, the goal is to judge whether the article is to be verified or not. For this task, we create a publicly available dataset in Japanese and provide benchmark results by using several basic machine learning techniques. Experimental results show that our models can reduce the cost of manual fact-checking process.


title: Calibration of imperfect geophysical models by multiple satellite interferograms with measurement bias

abstract: Model calibration consists of using experimental or field data to estimate the unknown parameters of a mathematical model. The presence of model discrepancy and measurement bias in the data complicates this task. Satellite interferograms, for instance, are widely used for calibrating geophysical models in geological hazard quantification. In this work, we used satellite interferograms to relate ground deformation observations to the properties of the magma chamber at K\={\i}lauea Volcano in Hawai`i. We derived closed-form marginal likelihoods and implemented posterior sampling procedures that simultaneously estimate the model discrepancy of physical models, and the measurement bias from the atmospheric error in satellite interferograms. We found that model calibration by aggregating multiple interferograms and downsampling the pixels in the interferograms can reduce the computation complexity compared to calibration approaches based on multiple data sets. The conditions that lead to no loss of information from data aggregation and downsampling are studied. Simulation illustrates that both discrepancy and measurement bias can be estimated, and real applications demonstrate that modeling both effects helps obtain a reliable estimation of a physical model's unobserved parameters and enhance its predictive accuracy. We implement the computational tools in the RobustCalibration package available on CRAN.


title: Towards Smart City Innovation Under the Perspective of Software-Defined Networking, Artificial Intelligence and Big Data

abstract: Smart city projects address many of the current problems afflicting high populated areas and cities and, as such, are a target for government, institutions and private organizations that plan to explore its foreseen advantages. In technical terms, smart city projects present a complex set of requirements including a large number users with highly different and heterogeneous requirements. In this scenario, this paper proposes and analyses the impact and perspectives on adopting software-defined networking and artificial intelligence as innovative approaches for smart city project development and deployment. Big data is also considered as an inherent element of most smart city project that must be tackled. A framework layered view is proposed with a discussion about software-defined networking and machine learning impacts on innovation followed by a use case that demonstrates the potential benefits of cognitive learning for smart cities. It is argued that the complexity of smart city projects do require new innovative approaches that potentially result in more efficient and intelligent systems.


title: Viscous Self Interacting Dark Matter Cosmology For Small Redshift

abstract: The viscosity of dark matter in cosmological models may cause an accelerated expansion and when this effect is sufficiently large, it can explain the dark energy. In this work, attributing the origin of viscosity to self-interaction of dark matter, we study the viscous cosmology at small redshift $(0\leq z\leq2.5)$. Assuming the cluster scale to be virialized and by modeling a power law behavior of velocity gradients, we calculate the Hubble expansion rate, $H(z)$ and the deceleration parameter, $q(z)$. We then perform a $\chi^{2}$ analysis to estimate the best fit model parameters. By using the best fit values, we explain the cosmic chronometer and type Ia supernova data. We conclude that if the dissipative effects become prominent only at the late time of cosmic evolution and are smaller at higher redshift, we can explain the observational data without requiring any dark energy component. Our analysis is independent of any specific model of self interacting dark matter.


title: Microwave control of thermal magnon spin transport

abstract: We observe that an rf microwave field strongly influences the transport of incoherent thermal magnons in yttrium iron garnet. Ferromagnetic resonance in the nonlinear regime suppresses thermal magnon transport by 95%. The transport is also modulated at non-resonant conditions in two cases, both related to the magnon band minimum. Firstly, a strong enhancement of the nonlocal signal appears at a static magnetic field below the resonance condition. This increase only occurs at one field polarity and can be as large as 800%. We attribute this effect to magnon kinetic processes, which give rise to band-minimum magnons and high-energy chiral surface modes. Secondly, the signal increases at a static field above the resonance condition, where the rf frequency coincides with the magnon band minimum. Our study gives insight into the interplay between coherent and incoherent spin dynamics: The rf field modifies the occupation of relevant magnon states and, via kinetic processes, the magnon spin transport.


title: Extreme asteroids in the Pan-STARRS 1 Survey

abstract: Using the first 18 months of the Pan-STARRS 1 survey we have identified 33 candidate high-amplitude objects for follow-up observations and carried out observations of 22 asteroids. 4 of the observed objects were found to have observed amplitude $A_{obs}\geq 1.0$ mag. We find that these high amplitude objects are most simply explained by single rubble pile objects with some density-dependent internal strength, allowing them to resist mass shedding even at their highly elongated shapes. 3 further objects although below the cut-off for 'high-amplitude' had a combination of elongation and rotation period which also may require internal cohesive strength, depending on the density of the body. We find that none of the 'high-amplitude asteroids' identified here require any unusual cohesive strengths to resist rotational fission. 3 asteroids were sufficiently observed to allow for shape and spin pole models to be determined through light curve inversion. 45864 was determined to have retrograde rotation with spin pole axes $\lambda=218\pm 10^{\circ}, \beta=-82\pm 5^{\circ}$ and asteroid 206167 was found to have best fit spin pole axes $\lambda= 57 \pm 5^{\circ}$, $\beta=-67 \pm 5^{\circ}$. An additional object not initially measured with $A_{obs}>1.0$ mag, 49257, was determined to have a shape model which does suggest a high-amplitude object. Its spin pole axes were best fit for values $\lambda=112\pm 6^{\circ}, \beta=6\pm 5^{\circ}$. In the course of this project to date no large super-fast rotators ($P_{rot} < 2.2$ h) have been identified.


title: Merging the A- and Q-spectral theories for digraphs

abstract: Let $G$ be a digraph and $A(G)$ be the adjacency matrix of $G$. Let $D(G)$ be the diagonal matrix with outdegrees of vertices of $G$. For any real $\alpha\in[0,1]$, Liu et al. \cite{LWCL} defined the matrix $A_\alpha(G)$ as $$A_\alpha(G)=\alpha D(G)+(1-\alpha)A(G).$$ The largest modulus of the eigenvalues of $A_\alpha(G)$ is called the $A_\alpha$ spectral radius of $G$. In this paper, we determine the digraphs which attain the maximum (or minimum) $A_\alpha$ spectral radius among all strongly connected digraphs with given parameters such as girth, clique number, vertex connectivity or arc connectivity. We also discuss a number of open problems.


title: Thermo-micro-mechanical simulation of bulk metal forming processes

abstract: The newly proposed microstructural constitutive model for polycrystal viscoplasticity in cold and warm regimes (Motaman and Prahl, 2019), is implemented as a microstructural solver via user-defined material subroutine in a finite element (FE) software. Addition of the microstructural solver to the default thermal and mechanical solvers of a standard FE package enabled coupled thermo-micro-mechanical or thermal-microstructural-mechanical (TMM) simulation of cold and warm bulk metal forming processes. The microstructural solver, which incrementally calculates the evolution of microstructural state variables (MSVs) and their correlation to the thermal and mechanical variables, is implemented based on the constitutive theory of isotropic hypoelasto-viscoplastic (HEVP) finite (large) strain/deformation. The numerical integration and algorithmic procedure of the FE implementation are explained in detail. Then, the viability of this approach is shown for (TMM-) FE simulation of an industrial multistep warm forging.


title: A no-regret generalization of hierarchical softmax to extreme multi-label classification

abstract: Extreme multi-label classification (XMLC) is a problem of tagging an instance with a small subset of relevant labels chosen from an extremely large pool of possible labels. Large label spaces can be efficiently handled by organizing labels as a tree, like in the hierarchical softmax (HSM) approach commonly used for multi-class problems. In this paper, we investigate probabilistic label trees (PLTs) that have been recently devised for tackling XMLC problems. We show that PLTs are a no-regret multi-label generalization of HSM when precision@k is used as a model evaluation metric. Critically, we prove that pick-one-label heuristic - a reduction technique from multi-label to multi-class that is routinely used along with HSM - is not consistent in general. We also show that our implementation of PLTs, referred to as extremeText (XT), obtains significantly better results than HSM with the pick-one-label heuristic and XML-CNN, a deep network specifically designed for XMLC problems. Moreover, XT is competitive to many state-of-the-art approaches in terms of statistical performance, model size and prediction time which makes it amenable to deploy in an online system.


title: Average Convergence Rate of Evolutionary Algorithms II: Continuous Optimization

abstract: The average convergence rate (ACR) measures how fast the approximation error of an evolutionary algorithm converges to zero per generation. It is defined as the geometric average of the reduction rate of the approximation error over consecutive generations. This paper makes a theoretical analysis of the ACR in continuous optimization. The obtained results are summarized as follows. According to the limit property, the ACR is classified into two categories: (1) linear ACR whose limit inferior value is larger than a positive and (2) sublinear ACR whose value converges to zero. Then, it is proven that the ACR is linear for evolutionary programming using positive landscape-adaptive mutation, but sublinear for that using landscape-invariant or zero landscape-adaptive mutation. The relationship between the ACR and the decision space dimension is also classified into two categories: (1) polynomial ACR whose value is larger than the reciprocal of a polynomial function of the dimension for any generation, and (2) exponential ACR whose value is less than the reciprocal of an exponential function of the dimension for an exponential long period. It is proven that for easy problems such as linear functions, the ACR of the (1+1) adaptive random univariate search is polynomial. But for hard functions such as the deceptive function, the ACR of both the (1+1) adaptive random univariate search and evolutionary programming is exponential.


title: Beyond Structured Programming

abstract: The correctness of a structured program is, at best, plausible. Though this is a step forward compared to what came before, it falls short of verified correctness. To verify a structured program according to Hoare's method one is faced with the problem of finding assertions to fit existing code. In 1971 this mode of verification was declared by Dijkstra as too hard to be of practical use---he advised that proof and code were to grow together. A method for doing this was independently published by Reynolds in 1978 and by van Emden in 1979. The latter was further developed to attain the form of matrix code. This form of code not only obviates the need of fitting assertions to existing code, but helps in discovering an algorithm that reaches a given postcondition from a fixed precondition. In this paper a keyboard-editable version of matrix code is presented that uses E.W. Dijkstra's guarded commands as starting point. The result is reached by using Floyd's method rather than Hoare's as starting point.


title: Serendipitous X-ray Sources in the Chandra HRC Field around Alpha Centauri

abstract: For more than a decade, Alpha Centauri AB (G2V+K1V) has been observed by Chandra, in a long-term program to follow coronal (T~1MK) activity cycles of the two sunlike stars. Over 2008.4-2017.8, nineteen HRC-I exposures were taken, each about 10 ks in duration, and spaced about six months apart. Beyond monitoring the AB X-ray luminosities, the HRC-I sequence represents a unique decadal record of the dozen, or so, serendipitous X-ray sources in the Alpha Cen field, which is at low Galactic latitude and thus dominated by nearby stars. For the present study, the ten brightest candidates were considered. Only a handful of these were persistent; most were variable, some highly so, flaring in a few epochs, weak or absent in the others. All ten X-ray sources have Gaia objects within about 2 arcseconds; mostly late-type dwarfs, but a few giants. However, two of the proposed optical counterparts have statistically significant offsets, and possible conflicts between the X-ray and optical properties. Another of the candidates brightened a factor of 100 in X-rays during a single exposure, briefly attaining super-flare status. The Gaia counterpart is anomalously blue for its absolute G-magnitude, and likely is a WD+dM pair. To the extent that the low Galactic latitude field is representative, the Chandra time-domain view emphasizes that the high-energy stellar sky is biased toward transient sources, so any snapshot survey surely will miss many of the most interesting objects.


title: Resource Control in P2P Cryptocurrency Networks

abstract: For decentralised P2P networks, it is very important to have a mechanism in place that allows the nodes to control resource usage and prevent flooding and denial-of-service attacks with spam. In this paper, we discuss and compare the different approaches to fully decentralised resource control that are used by projects in the cryptocurrency space. The introduced methods are then applied to design a decentralised exchange for Namecoin names (or more generally, crypto assets) as an example.


title: Periodicity of Multidimensional Continued Fractions

abstract: It is known that the continued fraction expansion of a real number is periodic if and only if the number is a quadratic irrational. In an attempt to generalize this phenomenon to other settings, Jun-Ichi Tamura and Shin-Ichi Yasutomi have developed a new algorithm for multidimensional continued fractions (Algebraic Jacobi-Perron algorithm) that involves cubic irrationals, and proved periodicity in some cubic number fields, such as $\mathbb{Q}(\sqrt[3]{m^3+1})$ where $m\in\mathbb{Z}$, and $\mathbb{Q}(\delta_m)$ where $\delta_m$ is a root of $x^3-mx+1=0,\,\,m\in\mathbb{Z},\,\, m\geq3$ with the algorithm. In this paper, we study some other types of number fields that give rise to periodic continued fractions using the Algebraic Jacobi-Perron algorithm obtaining results for $\mathbb{Q}(\sqrt[l]{m^l+1})$ for any positive integer $l$. Furthermore, we find that some families of cubic equations, such as $x^3+3ax^2+bx+ab-2a^3+1=0,\,b\leq3a^2-3,\,a,b\in\mathbb{Z}$, have roots that have periodic multidimensional continued fractions.


title: The Variational Deficiency Bottleneck

abstract: We introduce a bottleneck method for learning data representations based on information deficiency, rather than the more traditional information sufficiency. A variational upper bound allows us to implement this method efficiently. The bound itself is bounded above by the variational information bottleneck objective, and the two methods coincide in the regime of single-shot Monte Carlo approximations. The notion of deficiency provides a principled way of approximating complicated channels by relatively simpler ones. We show that the deficiency of one channel with respect to another has an operational interpretation in terms of the optimal risk gap of decision problems, capturing classification as a special case. Experiments demonstrate that the deficiency bottleneck can provide advantages in terms of minimal sufficiency as measured by information bottleneck curves, while retaining robust test performance in classification tasks.


title: Applications of Envelopes

abstract: Intuitively, an envelope of a family of curves is a curve that is tangent to a member of the family at each point. Here we use envelopes of families of circles to study objects from matrix theory and hyperbolic geometry. First we explore relationships between numerical ranges of $2\times 2$ matrices and families of circles to study the elliptical range theorem. Then we deduce a relationship between envelopes and the boundaries of families of intersecting circles and use it to find the boundaries of various families of pseudohyperbolic disks.


title: Saddle-Node Bifurcation of Periodic Orbits for a Delay Differential Equation

abstract: We consider the scalar delay differential equation $$ \dot{x}(t)=-x(t)+f_{K}(x(t-1)) $$ with a nondecreasing feedback function $f_{K}$ depending on a parameter $K$, and we verify that a saddle-node bifurcation of periodic orbits takes place as $K$ varies. The nonlinearity $f_{K}$ is chosen so that it has two unstable fixed points (hence the dynamical system has two unstable equilibria), and these fixed points remain bounded away from each other as $K$ changes. The generated periodic orbits are of large amplitude in the sense that they oscillate about both unstable fixed points of $f_{K}$.


title: Numerical Range and Compressions of the Shift

abstract: The numerical range of a bounded, linear operator on a Hilbert space is a set in $\mathbb{C}$ that encodes important information about the operator. In this survey paper, we first consider numerical ranges of matrices and discuss several connections with envelopes of families of curves. We then turn to the shift operator, perhaps the most important operator on the Hardy space $H^2(\mathbb{D})$, and compressions of the shift operator to model spaces, i.e.~spaces of the form $H^2 \ominus \theta H^2$ where $\theta$ is inner. For these compressions of the shift operator, we provide a survey of results on the connection between their numerical ranges and the numerical ranges of their unitary dilations. We also discuss related results for compressed shift operators on the bidisk associated to rational inner functions and conclude the paper with a brief discussion of the Crouzeix conjecture.


title: The C-Band All-Sky Survey (C-BASS): Constraining diffuse Galactic radio emission in the North Celestial Pole region

abstract: The C-Band All-Sky Survey C-BASS is a high-sensitivity all-sky radio survey at an angular resolution of 45 arcmin and a frequency of 4.7 GHz. We present a total intensity 4.7 GHz map of the North Celestial Pole (NCP) region of sky, above declination +80 deg, which is limited by source confusion at a level of ~0.6 mK rms. We apply the template-fitting (cross-correlation) technique to WMAP and Planck data, using the C-BASS map as the synchrotron template, to investigate the contribution of diffuse foreground emission at frequencies ~20-40 GHz. We quantify the anomalous microwave emission (AME) that is correlated with far-infrared dust emission. The AME amplitude does not change significantly (<10%) when using the higher frequency C-BASS 4.7 GHz template instead of the traditional Haslam 408 MHz map as a tracer of synchrotron radiation. We measure template coefficients of $9.93\pm0.35$ and $9.52\pm0.34$ K per unit $\tau_{353}$ when using the Haslam and C-BASS synchrotron templates, respectively. The AME contributes $55\pm2\,\mu$K rms at 22.8 GHz and accounts for ~60% of the total foreground emission. Our results suggest that a harder (flatter spectrum) component of synchrotron emission is not dominant at frequencies >5 GHz; the best-fitting synchrotron temperature spectral index is $\beta=-2.91\pm0.04$ from 4.7 to 22.8 GHz and $\beta=-2.85\pm0.14$ from 22.8 to 44.1 GHz. Free-free emission is weak, contributing ~$7\,\mu$K rms (~7%) at 22.8 GHz. The best explanation for the AME is still electric dipole emission from small spinning dust grains.


title: ERS approximation for solving Schr\"odinger's equation and applications

abstract: A new technique was recently developed to approximate the solution of the Schroedinger equation. This approximation (dubbed ERS) is shown to yield a better accuracy than the WKB-approximation. Here, we review the ERS approximation and its application to one and three-dimensional systems. In particular, we treat bound state solutions. We further focus on random potentials in a quantum wire and discuss the solution in the context of Anderson localization.


title: Inhomogeneous Reionization Models in Cosmological Hydrodynamical Simulations

abstract: In this work we present a new hybrid method to simulate the thermal effects of the reionization in cosmological hydrodynamical simulations. The method improves upon the standard approach used in simulations of the intergalactic medium (IGM) and galaxy formation without a significant increase of the computational cost allowing for efficient exploration of the parameter space. The method uses a small set of phenomenological input parameters and combines a semi-numerical reionization model to solve for the topology of reionization and an approximate model of how reionization heats the IGM, with the massively parallel \texttt{Nyx} hydrodynamics code, specifically designed to solve for the structure of diffuse IGM gas. We have produced several large-scale high resolution cosmological hydrodynamical simulations ($2048^3$, $L_{\rm box} = 40$ Mpc/h) with different instantaneous and inhomogeneous HI reionization models that use this new methodology. We study the IGM thermal properties of these models and find that large scale temperature fluctuations extend well beyond the end of reionization. Analyzing the 1D flux power spectrum of these models, we find up to $\sim 50\%$ differences in the large scale properties (low modes, $k\lesssim0.01$ s/km) of the post-reionization power spectrum due to the thermal fluctuations. We show that these differences could allow one to distinguish between different reionization scenarios already with existing Ly$\alpha$ forest measurements. Finally, we explore the differences in the small-scale cutoff of the power spectrum and we find that, for the same heat input, models show very good agreement provided that the reionization redshift of the instantaneous reionization model happens at the midpoint of the inhomogeneous model.


title: Element sets for high-order Poincar\'e mapping of perturbed Keplerian motion

abstract: The propagation and Poincar\'e mapping of perturbed Keplerian motion is a key topic in celestial mechanics and astrodynamics, e.g. to study the stability of orbits or design bounded relative trajectories. The high-order transfer map (HOTM) method enables efficient mapping of perturbed Keplerian orbits over many revolutions. For this, the method uses the high-order Taylor expansion of a Poincar\'e or stroboscopic map, which is accurate close to the expansion point. In this paper, we investigate the performance of the HOTM method using different element sets for building the high-order map. The element sets investigated are the classical orbital elements, modified equinoctial elements, Hill variables, cylindrical coordinates and Deprit's ideal elements. The performances of the different coordinate sets are tested by comparing the accuracy and efficiency of mapping low-Earth and highly-elliptical orbits perturbed by $J_2$ with numerical propagation. The accuracy of HOTM depends strongly on the choice of elements and type of orbit. A new set of elements is introduced that enables extremely accurate mapping of the state, even for high eccentricities and higher-order zonal perturbations. Finally, the high-order map is shown to be very useful for the determination and study of fixed points and centre manifolds of Poincar\'e maps.


title: Direct Quantitative Photoacoustic Tomography for realistic acoustic media

abstract: Quantitative photo-acoustic tomography (QPAT) seeks to reconstruct a distribution of optical attenuation coefficients inside a sample from a set of time series of pressure data that is measured outside the sample. The associated inverse problems involve two steps, namely acoustic and optical, which can be solved separately or as a direct composite problem. We adopt the latter approach for realistic acoustic media that possess heterogeneous and often not accurately known distributions for sound speed and ambient density, as well as an attenuation following a frequency power law that is evident in tissue media. We use a Diffusion Approximation (DA) model for the optical portion of the problem. We solve the corresponding composite inverse problem using three total variation (TV) regularised optimisation approaches. Accordingly, we develop two Krylov-subspace inexact-Newton algorithms that utilise the Jacobian matrix in a matrix-free manner in order to handle the computational cost. Additionally, we use a gradient-based algorithm that computes a search direction using the L-BFGS method, and applies a TV regularisation based on the Alternating Direction Method of Multipliers (ADMM) as a benchmark, because this method is popular for QPAT and direct QPAT. The results indicate the superiority of the developed inexact Newton algorithms over gradient-based Quasi-Newton approaches for a comparable computational complexity.


title: Evolution of ultra-relativistic hollow-electron-beam wakefield drivers during their propagation in plasmas

abstract: Ultra-relativistic hollow electron beams can drive plasma wakefields ($\sim$ GV/m) suitable for positron acceleration. Stable propagation of hollow electron beams for long distances in plasmas is required to accelerate positrons to high energies by these plasma wakefields. In this work, we show by quasi-static kinetic simulations using the code WAKE that an ultra-relativistic azimuthally-symmetric hollow electron beam propagates in a plasma by developing fish-bone like structure and shifting its bulk, differentially along its length (rear part fastest), towards its axis due to the decrease in the betatron time period of beam electrons from the beam-front to beam-rear. Hollow electron beams with small radius collapse into their axis due to the pull by the secondary wakefields generated by some of the beam electrons reaching the axis. Hollow beams with sufficiently large radius, however, can propagate stably in plasmas for several meters and be used for positron acceleration.


title: Constraining minimal Type-III Seesaw Model from naturalness, Lepton Flavor Violation and Electroweak Vacuum stability

abstract: We study the minimal type-III seesaw model in which we extend the SM by adding two $SU(2)_L$ triplet fermions with zero hypercharge to explain the origin of the non-zero neutrino masses. We show that the naturalness conditions and the limits from lepton flavor violating decays provide very stringent bounds on the model parameters along with the constraints from the stability/metastability of the electroweak vacuum. We perform a detailed analysis of the model parameter space including all the constraints for both normal as well as inverted hierarchies of the light neutrino masses. We find that most of the region that are allowed by lepton flavor violating decays and naturalness fall in the stable/metastable region depending on the values of the standard model parameters.


title: Marine crustaceans with hairy appendages: role of hydrodynamic boundary layers in sensing and feeding

abstract: Decapod crustaceans have appendages with an array of rigid hairs covered in chemo-receptors, used to sense and track food. Crustaceans directly influence the flow behavior by changing the speed of flow past the hairy surface, thereby manipulating the Reynolds number (Re). Hairs act either as a rake - diverting flow around the hair array, or as a sieve - filtering flow through the hairs. In our experiments, we uncover a third transitional phase: deflection - where the flow partially penetrates the hair array and is deflected laterally. We develop a reduced order model that predicts the flow phase based on the depth of the boundary layer on a single hair. This model with no fitting parameters agrees very well with our experimental data. Additionally, our model agrees well with measurements of both chemo-sensing and suspension-feeding crustaceans, and can be generalized for many different geometries.


title: Accelerated Inference in Markov Random Fields via Smooth Riemannian Optimization

abstract: Markov Random Fields (MRFs) are a popular model for several pattern recognition and reconstruction problems in robotics and computer vision. Inference in MRFs is intractable in general and related work resorts to approximation algorithms. Among those techniques, semidefinite programming (SDP) relaxations have been shown to provide accurate estimates while scaling poorly with the problem size and being typically slow for practical applications. Our first contribution is to design a dual ascent method to solve standard SDP relaxations that takes advantage of the geometric structure of the problem to speed up computation. This technique, named Dual Ascent Riemannian Staircase (DARS), is able to solve large problem instances in seconds. Our second contribution is to develop a second and faster approach. The backbone of this second approach is a novel SDP relaxation combined with a fast and scalable solver based on smooth Riemannian optimization. We show that this approach, named Fast Unconstrained SEmidefinite Solver (FUSES), can solve large problems in milliseconds. Contrarily to local MRF solvers, e.g., loopy belief propagation, our approaches do not require an initial guess. Moreover, we leverage recent results from optimization theory to provide per-instance sub-optimality guarantees. We demonstrate the proposed approaches in multi-class image segmentation problems. Extensive experimental evidence shows that (i) FUSES and DARS produce near-optimal solutions, attaining an objective within 0.1% of the optimum, (ii) FUSES and DARS are remarkably faster than general-purpose SDP solvers, and FUSES is more than two orders of magnitude faster than DARS while attaining similar solution quality, (iii) FUSES is faster than local search methods while being a global solver.


title: 3D Terrain Segmentation in the SWIR Spectrum

abstract: We focus on the automatic 3D terrain segmentation problem using hyperspectral shortwave IR (HS-SWIR) imagery and 3D Digital Elevation Models (DEM). The datasets were independently collected, and metadata for the HS-SWIR dataset are unavailable. We explore an overall slope of the SWIR spectrum that correlates with the presence of moisture in soil to propose a band ratio test to be used as a proxy for soil moisture content to distinguish two broad classes of objects: live vegetation from impermeable manmade surface. We show that image based localization techniques combined with the Optimal Randomized RANdom Sample Consensus (RANSAC) algorithm achieve precise spatial matches between HS-SWIR data of a portion of downtown Los Angeles (LA (USA)) and the Visible image of a geo-registered 3D DEM, covering a wider-area of LA. Our spectral-elevation rule based approach yields an overall accuracy of 97.7%, segmenting the object classes into buildings, houses, trees, grass, and roads/parking lots.


title: An inverse technique for reconstructing ocean's density stratification from surface data

abstract: In this article, we propose an inverse technique that accurately reconstructs the ocean's density stratification profile simply from free surface elevation data. Satellite observations suggest that ocean surface contains the signature of internal tides, which are internal gravity waves generated by the barotropic tides. Since internal tides contain the information of ocean's density stratification, the latter can in principle be reconstructed from the free surface signature. We consider a simple theoretical model that approximates a continuously stratified ocean as discrete layers of constant buoyancy frequency; this facilitates the derivation of a closed-form dispersion relation. First, we numerically simulate internal tide generation for toy ocean scenarios and subsequently perform Space-Time Fourier Transform (STFT) of the free surface, which yields internal tide spectra with wavenumbers corresponding to the tidal frequency. The density profile is reconstructed by substituting these wavenumbers into the dispersion relation. Finally, we consider a more realistic situation with rotation, bottom topography, shear and density profiles representative of the Strait of Gibraltar. Density reconstruction in the presence and absence of shear are respectively found to be $90.2\%$ and $94.2\%$ accurate.


title: Modeling Perceptual Aliasing in SLAM via Discrete-Continuous Graphical Models

abstract: Perceptual aliasing is one of the main causes of failure for Simultaneous Localization and Mapping (SLAM) systems operating in the wild. Perceptual aliasing is the phenomenon where different places generate a similar visual (or, in general, perceptual) footprint. This causes spurious measurements to be fed to the SLAM estimator, which typically results in incorrect localization and mapping results. The problem is exacerbated by the fact that those outliers are highly correlated, in the sense that perceptual aliasing creates a large number of mutually-consistent outliers. Another issue stems from the fact that most state-of-the-art techniques rely on a given trajectory guess (e.g., from odometry) to discern between inliers and outliers and this makes the resulting pipeline brittle, since the accumulation of error may result in incorrect choices and recovery from failures is far from trivial. This work provides a unified framework to model perceptual aliasing in SLAM and provides practical algorithms that can cope with outliers without relying on any initial guess. We present two main contributions. The first is a Discrete-Continuous Graphical Model (DC-GM) for SLAM: the continuous portion of the DC-GM captures the standard SLAM problem, while the discrete portion describes the selection of the outliers and models their correlation. The second contribution is a semidefinite relaxation to perform inference in the DC-GM that returns estimates with provable sub-optimality guarantees. Experimental results on standard benchmarking datasets show that the proposed technique compares favorably with state-of-the-art methods while not relying on an initial guess for optimization.


title: Stein Variational Gradient Descent as Moment Matching

abstract: Stein variational gradient descent (SVGD) is a non-parametric inference algorithm that evolves a set of particles to fit a given distribution of interest. We analyze the non-asymptotic properties of SVGD, showing that there exists a set of functions, which we call the Stein matching set, whose expectations are exactly estimated by any set of particles that satisfies the fixed point equation of SVGD. This set is the image of Stein operator applied on the feature maps of the positive definite kernel used in SVGD. Our results provide a theoretical framework for analyzing the properties of SVGD with different kernels, shedding insight into optimal kernel choice. In particular, we show that SVGD with linear kernels yields exact estimation of means and variances on Gaussian distributions, while random Fourier features enable probabilistic bounds for distributional approximation. Our results offer a refreshing view of the classical inference problem as fitting Stein's identity or solving the Stein equation, which may motivate more efficient algorithms.


title: Resolutions of ideals associated to subspace arrangements

abstract: Given a collection of $t$ subspaces in an $n$-dimensional $\mathbb{K} $-vector space $W$ we can associate to them $t$ vanishing ideals in the symmetric algebra $\mathcal{S}(W^*) = \mathbb{K}[x_1,x_2,\dots,x_n]$. As a subspace is defined by a set of linear equations, its vanishing ideal is generated by linear forms so it is a linear ideal. Conca and Herzog showed that the Castelnuovo-Mumford regularity of the product of $t$ linear ideals is equal to $t$. Derksen and Sidman showed that the Castelnuovo-Mumford regularity of the intersection of $t$ linear ideals is at most $t$ and they also showed that similar results hold for a more general class of ideals constructed from linear ideals. In this paper we show that analogous results hold when we replace the symmetric algebra $\mathcal{S}(W^*)$ with the exterior algebra $ \bigwedge(W^*)$ and work over a field of characteristic 0. To prove these results we rely on the functoriality of free resolutions and construct a functor $\Omega$ from the category of polynomial functors to itself. The functor $\Omega$ transforms resolutions of ideals in the symmetric algebra to resolutions of ideals in the exterior algebra.


title: Designing Refund Bonus Schemes for Provision Point Mechanism in Civic Crowdfunding

abstract: Civic crowdfunding (CC) is a popular medium for raising funds for civic projects from interested agents. With Blockchains gaining traction, we can implement CC in a reliable, transparent, and secure manner with smart contracts (SCs). The fundamental challenge in CC is free-riding. PPR, the proposal by Zubrickas [23] of giving refund bonus to the contributors, in the case of the project not getting provisioned, has attractive properties. However, as observed by Chandra et al. [7], PPR faces a challenge wherein the agents defer their contribution until the deadline. We define this delaying of contributions as a race condition. To address this, their proposal, PPS, considers the temporal aspects of a contribution. However, PPS is computationally complex, expensive to implement as an SC, and it being sophisticated, it is difficult to explain to a layperson. In this work, our goal is to identify all essential properties a refund bonus scheme must satisfy in order to curb free-riding while avoiding the race condition. We prove Contribution Monotonicity and Time Monotonicity are sufficient conditions for this. We propose three elegant refund bonus schemes satisfying these two conditions leading to three novel mechanisms for CC - PPRG, PPRE, and PPRP. We show that PPRG is the most cost-effective mechanism when deployed as an SC. We show that under certain modest assumptions on valuations of the agents, in PPRG, the project is funded at equilibrium.


title: In-situ liquid SAXS studies on the early stage of calcium carbonate formation

abstract: Calcium carbonate is a model system to investigate the mechanism of solid formation by precipitation from solutions, and it is often considered in the debated classical and non-classical nucleation mechanism. Despite the great scientific relevance of calcium carbonate in different areas of science, little is known about the early stage of its formation. We, therefore, designed contactless devices capable to provide informative investigations on the early stages of the precipitation pathway of calcium carbonate in supersaturated solutions using classical scattering methods such as Wide-Angle X-ray Scattering (WAXS) and Small-Angle X-ray Scattering (SAXS) techniques. In particular, SAXS was exploited for investigating the size of entities formed from supersaturated solutions before the critical conditions for amorphous calcium carbonate (ACC) nucleation are attained. The saturation level was controlled by mixing four diluted solutions (i.e., NaOH, CaCl2, NaHCO3, H2O) at constant T and pH. The scattering data were collected on a liquid jet generated about 75 sec after the mixing point. The data were modeled using parametric statistical models providing insight about the size distribution of denser matter in the liquid jet. Theoretical implications on the early stage of solid formation pathway are inferred.


title: On a categoriacal theory for emergence

abstract: It is well-known that biological phenomena are emergent. Emergent phenomena are quite interesting and amazing. However, they are difficult to be understood. Due to this difficulty, we propose a theory to describe emergence based on a powerful mathematical tool, namely, Theory of Categories. In order to do this, we first utilize constructs (categories whose objects are structured sets), their operations and their corresponding generalized underlying functor (which are not necessary faithful) to characterize emergence. After this, we introduce and show several results concerning homomorphism (isomorphism) between emergences, representability, pullback, pushout, equalizer, product and co-product of emergences among other concepts. Finally, we explain how our theory fits in studies involving biological systems.


title: Uncertain Trees: Dealing with Uncertain Inputs in Regression Trees

abstract: Tree-based ensemble methods, as Random Forests and Gradient Boosted Trees, have been successfully used for regression in many applications and research studies. Furthermore, these methods have been extended in order to deal with uncertainty in the output variable, using for example a quantile loss in Random Forests (Meinshausen, 2006). To the best of our knowledge, no extension has been provided yet for dealing with uncertainties in the input variables, even though such uncertainties are common in practical situations. We propose here such an extension by showing how standard regression trees optimizing a quadratic loss can be adapted and learned while taking into account the uncertainties in the inputs. By doing so, one no longer assumes that an observation lies into a single region of the regression tree, but rather that it belongs to each region with a certain probability. Experiments conducted on several data sets illustrate the good behavior of the proposed extension.


title: The Swampland Spectrum Conjecture in Inflation

abstract: The quantum gravity conjectures that aim to separate the landscape from the swampland among the low energy theories were originally formulated in the context of scalar field spaces spanned by moduli. Because these conjectures have implications for cosmology they have recently been considered in a more general context for scalar field theories with potentials, in particular inflation. From an effective field theory perspective the presence of a potential induces a natural metric that makes the distance measure $D_V$ along scalar field trajectories dependent on the potential. The present paper proposes a modified formulation in terms of $D_V$ of those conjectures that involve trajectory distances.


title: Minimum Reload Cost Graph Factors

abstract: The concept of Reload cost in a graph refers to the cost that occurs while traversing a vertex via two of its incident edges. This cost is uniquely determined by the colors of the two edges. This concept has various applications in transportation networks, communication networks, and energy distribution networks. Various problems using this model are defined and studied in the literature. The problem of finding a spanning tree whose diameter with respect to the reload costs is the smallest possible, the problems of finding a path, trail or walk with minimum total reload cost between two given vertices, problems about finding a proper edge coloring of a graph such that the total reload cost is minimized, the problem of finding a spanning tree such that the sum of the reload costs of all paths between all pairs of vertices is minimized, and the problem of finding a set of cycles of minimum reload cost, that cover all the vertices of a graph, are examples of such problems. % In this work we focus on the last problem. Noting that a cycle cover of a graph is a 2-factor of it, we generalize the problem to that of finding an $r$-factor of minimum reload cost of an edge colored graph. We prove several NP-hardness results for special cases of the problem. Namely, bounded degree graphs, planar graphs, bounded total cost, and bounded number of distinct costs. For the special case of $r=2$, our results imply an improved NP-hardness result. On the positive side, we present a polynomial-time solvable special case which provides a tight boundary between the polynomial and hard cases in terms of $r$ and the maximum degree of the graph. We then investigate the parameterized complexity of the problem, prove W[1]-hardness results and present an FPT algorithm.


title: Hull Form Optimization with Principal Component Analysis and Deep Neural Network

abstract: Designing and modifying complex hull forms for optimal vessel performances have been a major challenge for naval architects. In the present study, Principal Component Analysis (PCA) is introduced to compress the geometric representation of a group of existing vessels, and the resulting principal scores are manipulated to generate a large number of derived hull forms, which are evaluated computationally for their calm-water performances. The results are subsequently used to train a Deep Neural Network (DNN) to accurately establish the relation between different hull forms and their associated performances. Then, based on the fast, parallel DNN-based hull-form evaluation, the large-scale search for optimal hull forms is performed.


title: Multi-Agent Common Knowledge Reinforcement Learning

abstract: Cooperative multi-agent reinforcement learning often requires decentralised policies, which severely limit the agents' ability to coordinate their behaviour. In this paper, we show that common knowledge between agents allows for complex decentralised coordination. Common knowledge arises naturally in a large number of decentralised cooperative multi-agent tasks, for example, when agents can reconstruct parts of each others' observations. Since agents an independently agree on their common knowledge, they can execute complex coordinated policies that condition on this knowledge in a fully decentralised fashion. We propose multi-agent common knowledge reinforcement learning (MACKRL), a novel stochastic actor-critic algorithm that learns a hierarchical policy tree. Higher levels in the hierarchy coordinate groups of agents by conditioning on their common knowledge, or delegate to lower levels with smaller subgroups but potentially richer common knowledge. The entire policy tree can be executed in a fully decentralised fashion. As the lowest policy tree level consists of independent policies for each agent, MACKRL reduces to independently learnt decentralised policies as a special case. We demonstrate that our method can exploit common knowledge for superior performance on complex decentralised coordination tasks, including a stochastic matrix game and challenging problems in StarCraft II unit micromanagement.


title: High-resolution abundance analysis of four red giants in the globular cluster NGC 6558

abstract: NGC 6558 is a bulge globular cluster with a blue horizontal branch (BHB), combined with a metallicity of [Fe/H]~-1.0. It is similar to HP 1 and NGC 6522, which could be among the oldest objects in the Galaxy. Element abundances in these clusters could reveal the nature of the first supernovae. We aim to carry out detailed spectroscopic analysis for four red giants of NGC 6558, in order to derive the abundances of the light elements C, N, O, Na, Al, the alpha-elements Mg, Si, Ca, Ti, and the heavy elements Y, Ba, and Eu. High-resolution spectra of four stars with FLAMES-UVES@VLT UT2-Kueyen were analysed. Spectroscopic parameter-derivation was based on excitation and ionization equilibrium of FeI and FeII. This analysis results in a metallicity of [Fe/H] = -1.17+-0.10 for NGC 6558. We find the expected alpha-element enhancements in O and Mg with [O/Fe]=+0.40, [Mg/Fe]=+0.33, and low enhancements in Si and Ca. Ti has a moderate enhancement of [Ti/Fe]=+0.22. The r-element Eu appears very enhanced with a mean value of [Eu/Fe]=+0.63. Ba appears to have a solar abundance ratio relative to Fe. NGC 6558 shows an abundance pattern that could be typical of the oldest inner bulge globular clusters, together with the pattern in the similar clusters NGC 6522 and HP 1. They show low abundances of the odd-Z elements Na and Al, and of the explosive nucleosynthesis alpha-elements Si, Ca, and Ti. The hydrostatic burning alpha-elements O and Mg are normally enhanced as expected in old stars enriched with yields from core-collapse supernovae, and the iron-peak elements Mn, Cu, Zn show low abundances, which is expected for Mn and Cu, but not for Zn. Finally, the cluster trio NGC 6558, NGC 6522, and HP 1 show similar abundance patterns.


title: Computational geometry and the U.S. Supreme Court

abstract: We use the United States Supreme Court as an illuminative context in which to discuss three different spatial voting preference models: an instance of the widely used single-peaked preferences, and two models that are more novel in which vote outcomes have a strength in addition to a location. We introduce each model from a formal axiomatic perspective, briefly discuss practical motivation for each in terms of judicial behavior, prove mathematical relationships among the voting coalitions compatible with each model, and then study the two-dimensional setting by presenting computational tools for working with the models and by exploring these with judicial voting data from the Supreme Court.


title: Wi-Motion: A Robust Human Activity Recognition Using WiFi Signals

abstract: Recent research has shown that human motions and positions can be recognized through WiFi signals. The key intuition is that different motions and positions introduce different multi-path distortions in WiFi signals and generate different patterns in the time-series of channel state information (CSI). In this paper, we propose Wi-Motion, a WiFi-based human activities recognition system. Unlike existing systems, Wi-Motion adopts the amplitude and phase information extracted from the CSI sequence to construct the classifiers respectively, and combines the results using a combination strategy based on posterior probability. As the simulation results shows, Wi-Motion can recognize six human activities with the mean accuracy of 98:4%.


title: On Measurement of the Spatio-Frequency Property of OFDM Backscattering

abstract: Orthogonal frequency-division multiplexing (OFDM) backscatter system, such as Wi-Fi backscatter, has recently been recognized as a promising technique for the IoT connectivity, due to its ubiquitous and low-cost property. This paper investigates the spatial-frequency property of the OFDM backscatter which takes the distance and the angle into account in different frequency bands. We deploy three typical scenarios for performing measurements to evaluate the received signal strength from the backscatter link. The impact of the distances among the transmitter, the tag and the receiver, as well as the angle between the transmitter and the tag is observed through the obtained measurement data. From the evaluation results, it is found that the best location of tag is either close to the receiver or the transmitter which depends on the frequency band, and the best angle is 90 degrees between the transmitter and the receiver. This work opens the shed light on the spatial deployment of the backscatter tag in different frequency band with the aim of improving the performance and reducing the interference.


title: From Communication to Sensing : Recognizing and Counting Repetitive Motions with Wireless Backscattering

abstract: Recently several ground-breaking RF-based motion recognition systems were proposed to detect and/or recognize macro/micro human movements. These systems often suffer from various interferences caused by multiple-users moving simultaneously, resulting in extremely low recognition accuracy. To tackle this challenge, we propose a novel system, called Motion-Fi, which marries battery free wireless backscattering and device-free sensing. Motion-Fi is an accurate, interference tolerable motion-recognition system, which counts repetitive motions without using scenario-dependent templates or profiles and enables multi-users performing certain motions simultaneously because of the relatively short transmission range of backscattered signals. Although the repetitive motions are fairly well detectable through the backscattering signals in theory, in reality they get blended into various other system noises during the motion. Moreover, irregular motion patterns among users will lead to expensive computation cost for motion recognition. We build a backscattering wireless platform to validate our design in various scenarios for over 6 months when different persons, distances and orientations are incorporated. In our experiments, the periodicity in motions could be recognized without any learning or training process, and the accuracy of counting such motions can be achieved within 5% count error. With little efforts in learning the patterns, our method could achieve 93.1% motion-recognition accuracy for a variety of motions. Moreover, by leveraging the periodicity of motions, the recognition accuracy could be further improved to nearly 100% with only 3 repetitions. Our experiments also show that the motions of multiple persons separating by around 2 meters cause little accuracy reduction in the counting process.


title: GD358: three decades of observations for the in-depth asteroseismology of a DBV star

abstract: We report on the analysis of 34 years of photometric observations of the pulsating helium atmosphere white dwarf GD358. The complete data set includes archival data from 1982-2006, and 1195.2 hours of new observations from 2007- 2016. From this data set, we extract 15 frequencies representing g-mode pulsation modes, adding 4 modes to the 11 modes known previously. We present evidence that these 15 modes are ell = 1 modes, 13 of which belong to a consecutive sequence in radial overtone k. We perform a detailed asteroseismic analysis using models that include parameterized, complex carbon and oxygen core composition profiles to fit the periods. Recent spectroscopic analyses place GD358 near the red edge of the DBV instability strip, at 24,000 plus or minus 500 K and a log g of 7.8 plus or minus 0.08 dex. The surface gravity translates to a mass range of 0.455 to 0.540 solar masses. Our best fit model has a temperature of 23,650 K and a mass of 0.5706 solar masses. That is slightly more massive than suggested by most the recent spectroscopy. We find a pure helium layer mass of 10^-5.50, consistent with the result of previous studies and the outward diffusion of helium over time.


title: Intelligent Nanophotonics: Merging Photonics and Artificial Intelligence at the Nanoscale

abstract: Nanophotonics has been an active research field over the past two decades, triggered by the rising interests in exploring new physics and technologies with light at the nanoscale. As the demands of performance and integration level keep increasing, the design and optimization of nanophotonic devices become computationally expensive and time-inefficient. Advanced computational methods and artificial intelligence, especially its subfield of machine learning, have led to revolutionary development in many applications, such as web searches, computer vision, and speech/image recognition. The complex models and algorithms help to exploit the enormous parameter space in a highly efficient way. In this review, we summarize the recent advances on the emerging field where nanophotonics and machine learning blend. We provide an overview of different computational methods, with the focus on deep learning, for the nanophotonic inverse design. The implementation of deep neural networks with photonic platforms is also discussed. This review aims at sketching an illustration of the nanophotonic design with machine learning and giving a perspective on the future tasks.


title: Flash Photography for Data-Driven Hidden Scene Recovery

abstract: Vehicles, search and rescue personnel, and endoscopes use flash lights to locate, identify, and view objects in their surroundings. Here we show the first steps of how all these tasks can be done around corners with consumer cameras. Recent techniques for NLOS imaging using consumer cameras have not been able to both localize and identify the hidden object. We introduce a method that couples traditional geometric understanding and data-driven techniques. To avoid the limitation of large dataset gathering, we train the data-driven models on rendered samples to computationally recover the hidden scene on real data. The method has three independent operating modes: 1) a regression output to localize a hidden object in 2D, 2) an identification output to identify the object type or pose, and 3) a generative network to reconstruct the hidden scene from a new viewpoint. The method is able to localize 12cm wide hidden objects in 2D with 1.7cm accuracy. The method also identifies the hidden object class with 87.7% accuracy (compared to 33.3% random accuracy). This paper also provides an analysis on the distribution of information that encodes the occluded object in the accessible scene. We show that, unlike previously thought, the area that extends beyond the corner is essential for accurate object localization and identification.


title: Regularization Effect of Fast Gradient Sign Method and its Generalization

abstract: Fast Gradient Sign Method (FGSM) is a popular method to generate adversarial examples that make neural network models robust against perturbations. Despite its empirical success, its theoretical property is not well understood. This paper develops theory to explain the regularization effect of Generalized FGSM, a class of methods to generate adversarial examples. Motivated from the relationship between FGSM and LASSO penalty, the asymptotic properties of Generalized FGSM are derived in the Generalized Linear Model setting, which is essentially the 1-layer neural network setting with certain activation functions. In such simple neural network models, I prove that Generalized FGSM estimation is root n-consistent and weakly oracle under proper conditions. The asymptotic results are also highly similar to penalized likelihood estimation. Nevertheless, Generalized FGSM introduces additional bias when data sampling is not sign neutral, a concept I introduce to describe the balance-ness of the noise signs. Although the theory in this paper is developed under simple neural network settings, I argue that it may give insights and justification for FGSM in deep neural network settings as well.


title: Normal real affine varieties with circle actions

abstract: We provide a complete description of normal affine algebraic varieties over the real numbers endowed with an effective action of the real circle, that is, the real form of the complex multiplicative group whose real locus consists of the unitary circle in the real plane. Our approach builds on the geometrico-combinatorial description of normal affine varieties with effective actions of split tori in terms of proper polyhedral divisors on semiprojective varieties due to Altmann and Hausen


title: Covariant Equations of Motion of Extended Bodies with Arbitrary Mass and Spin Multipoles

abstract: This paper employs the post-Newtonian approximations of scalar-tensor theory of gravity along with the Cartesian STF tensors and the Blanchet-Damour multipole formalism to derive translational and rotational equations of motion of N extended bodies with arbitrary distribution of mass and velocity. We assume that spacetime can be covered by a global coordinate chart which is Minkowskian at infinity. We also introduce N local coordinate charts adapted to each body and covering a finite domain of space around the body. Gravitational field of each body is parametrized by an infinite set of the body's mass and spin multipoles and tidal multipoles of external N-1 bodies. The origin of the local coordinates is set moving along accelerated worldline of the center of mass of the body by an appropriate choice of the internal and external dipole moments of its gravitational field. Translational and rotational equations of motion are derived by integrating microscopic equations of matter and applying the method of asymptotic matching. The matching is also used for separating the post-Newtonian self-field effects from the external gravitational environment and for constructing the effective background spacetime manifold. It allows us to present the equations of translational and rotational motion of each body in covariant form by making use of the Einstein principle of equivalence. Our approach significantly generalizes the Mathisson-Papapetrou-Dixon covariant equations of motion with regard to the number of body's multipoles and the post-Newtonian terms taken into account. The equations of translational and rotational motion derived in the present paper include the infinite set of mass and spin multipoles of the bodies and can be used for much more accurate prediction of orbital dynamics of extended bodies in inspiraling binary systems before they merge.


title: The CoSTAR Block Stacking Dataset: Learning with Workspace Constraints

abstract: A robot can now grasp an object more effectively than ever before, but once it has the object what happens next? We show that a mild relaxation of the task and workspace constraints implicit in existing object grasping datasets can cause neural network based grasping algorithms to fail on even a simple block stacking task when executed under more realistic circumstances. To address this, we introduce the JHU CoSTAR Block Stacking Dataset (BSD), where a robot interacts with 5.1 cm colored blocks to complete an order-fulfillment style block stacking task. It contains dynamic scenes and real time-series data in a less constrained environment than comparable datasets. There are nearly 12,000 stacking attempts and over 2 million frames of real data. We discuss the ways in which this dataset provides a valuable resource for a broad range of other topics of investigation. We find that hand-designed neural networks that work on prior datasets do not generalize to this task. Thus, to establish a baseline for this dataset, we demonstrate an automated search of neural network based models using a novel multiple-input HyperTree MetaModel, and find a final model which makes reasonable 3D pose predictions for grasping and stacking on our dataset. The CoSTAR BSD, code, and instructions are available at https://sites.google.com/site/costardataset.


title: Limit Cycles in a Model of Olfactory Sensory Neurons

abstract: We propose an approach to study small limit cycle bifurcations on a center manifold in analytic or smooth systems depending on parameters. We then apply it to the investigation of limit cycle bifurcations in a model of calcium oscillations in the cilia of olfactory sensory neurons and show that it can have two limit cycles: a stable cycle appearing after a Bautin (generalized Hopf) bifurcation and an unstable cycle appearing after a subcritical Hopf bifurcation.


title: Revealing the Broad Line Region of NGC 1275: The Relationship to Jet Power

abstract: NGC 1275 is one of the most conspicuous active galactic nuclei (AGN) in the local Universe. The radio jet currently emits a flux density of $\sim 10$ Jy at $\sim 1$ mm wavelengths, down from the historic high of $\sim 65$ Jy in 1980. Yet, the nature of the AGN in NGC 1275 is still controversial. It has been debated whether this is a broad emission line (BEL) Seyfert galaxy, an obscured Seyfert galaxy, a narrow line radio galaxy or a BL-Lac object. We clearly demonstrate a persistent H$\beta$ BEL over the last 35 years with a full width half maximum (FWHM) of 4150 - 6000 km/s. We also find a prominent P$\alpha$ BEL (FWHM $\approx 4770 $ km/s) and a weak CIV BEL (FWHM $\approx 4000 $ km/s), H$\beta$/CIV $\approx 2$. A far UV HST observation during suppressed jet activity reveals a low luminosity continuum. The H$\beta$ BEL luminosity is typical of broad line Seyfert galaxies with similar far UV luminosity. X-ray observations indicate a softer ionizing continuum than expected for a broad line Seyfert galaxy with similar far UV luminosity. This is opposite of the expectation of advection dominated accretion. The AGN continuum appears to be thermal emission from a low luminosity, optically thick, accretion flow with a low Eddington ratio, $\sim 0.0001$. The soft, weak ionizing continuum is consistent with the relatively weak CIV BEL. Evidence that the BEL luminosity is correlated with the jet mm wave luminosity is presented. Apparently, the accretion rate regulates jet power.


title: FBMC Prototype Filter Design via Convex Optimization

abstract: In this work, we propose a prototype filter design for Filter Bank MultiCarrier (FBMC) systems based on convex optimization, aiming superior spectrum features while maintaining a high symbol reconstruction quality. Initially, the proposed design is written as a non-convex Quadratically Constrained Quadratic Programming (QCQP), which is relaxed into a convex QCQP guided by a line search. Through the resulting problem, we design three prototype filters: Type-I, II and III. In particular, the Type-II filter shows a slightly better performance than the classical Mirabasi-Martin design, while Type-I and III filters offer a much more contained spectrum than most of the prototype filters suitable for FBMC applications. Furthermore, numerical results corroborate the effectiveness of the designed filters as the proposed filters offer fast decay and contained spectrum while not jeopardizing symbol reconstruction in practice.


title: On the role of ML estimation and Bregman divergences in sparse representation of covariance and precision matrices

abstract: Sparse representation of structured signals requires modelling strategies that maintain specific signal properties, in addition to preserving original information content and achieving simpler signal representation. Therefore, the major design challenge is to introduce adequate problem formulations and offer solutions that will efficiently lead to desired representations. In this context, sparse representation of covariance and precision matrices, which appear as feature descriptors or mixture model parameters, respectively, will be in the main focus of this paper.


title: Algorithmic information distortions and incompressibility in uniform multidimensional networks

abstract: This article presents a theoretical investigation of generalized encoded forms of networks in a uniform multidimensional space. First, we study encoded networks with (finite) arbitrary node dimensions (or aspects), such as time instants or layers. In particular, we study these networks that are formalized in the form of multiaspect graphs. In the context of node-aligned non-uniform (or node-unaligned non-uniform and uniform) multidimensional spaces, previous results has shown that, unlike classical graphs, the algorithmic information of a multidimensional network is not in general dominated by the algorithmic information of the binary sequence that determines the presence or absence of edges. In the present work, first we demonstrate the existence of such algorithmic information distortions for node-aligned uniform multidimensional networks. Secondly, we show that there are particular cases of infinite nesting families of finite uniform multidimensional networks such that each member of these families is incompressible. From these results, we also recover the network topological properties and equivalences in irreducible information content of multidimensional networks in comparison to their isomorphic classical graph counterpart in the previous literature. These results together establish a universal algorithmic approach and set limitations and conditions for irreducible information content analysis in comparing arbitrary networks with a large number of dimensions, such as multilayer networks.


title: Regularized integral representations of the reciprocal Gamma function

abstract: This paper establishes a real integral representation of the reciprocal $\Gamma$ function in terms of a regularized hypersingular integral. The equivalence with the usual complex representation is demonstrated. A regularized complex representation along the usual Hankel path is derived.


title: The B-Exponential Divergence and its Generalizations with Applications to Parametric Estimation

abstract: In this paper a new family of minimum divergence estimators based on the Bregman divergence is proposed, where the defining convex function has an exponential nature. These estimators avoid the necessity of using an intermediate kernel density and many of them also have strong robustness properties. It is further demonstrated that the proposed approach can be extended to construct a class of generalized estimating equations, where the pool of the resultant estimators encompass a large variety of minimum divergence estimators and range from highly robust to fully efficient based on the choice of the tuning parameters. All of the resultant estimators are M-estimators, where the defining functions make explicit use of the form of the parametric model. The properties of these estimators are discussed in detail; the theoretical results are substantiated by simulation and real data examples. It is observed that in many cases, certain robust estimators from the above generalized class provide better compromises between robustness and efficiency compared to the existing standards.


title: Generic Expansions by a Reduct

abstract: Consider the expansion $T_S$ of a theory $T$ by a predicate for a submodel of a reduct $T_0$ of $T$. We present a setup in which this expansion admits a model companion $TS$. We show that the nice features of the theory $T$ transfer to $TS$. In particular, we study conditions for which this expansion preserves the $\NSOP{1}$-ness, the simplicity or the stability of the starting theory $T$. We give concrete examples of new $\NSOP{1}$ not simple theories obtained by this process, among them the expansion of a perfect $\omega$-free $\PAC$ field of positive characteristic by generic additive subgroups, and the expansion of an algebraically closed field of \emph{any} characteristic by a generic multiplicative subgroup.


title: Nearly subadditive sequences

abstract: We show that the de Bruijn-Erd\H{o}s condition for the error term in their improvement of Fekete's Lemma is not only sufficient but also necessary in the following strong sense. Suppose that given a sequence $0\leq f(1)\leq f(2)\leq f(3)\leq \dots $ such that \begin{equation}\sum_{ n=1}^{\infty} f(n)/n^2 = \infty. \end{equation} Then, there exists a sequence $\{b(n)\}_{n=1,2,\dots}$ satisfying \begin{equation}\label{eq1} b(n+m) \leq b(n) + b(m) + f(n+m) \end{equation} such that the sequence of slopes $\{ b(n)/n\}_{n=1,2,\dots}$ takes every rational number. When the series is bounded we improve their result as follows. If there exist $N$ and real $\mu >1$ such that near $f$-subadditivity holds for all pairs $(n,m)$ with $N\leq n\leq m \leq \mu n$, then $\lim_n b(n)/n $ exists.


title: Skin optical attenuation coefficient modulation with arterial pulsation

abstract: Video plethysmographic (vPPG) signals arise from subtle color modulations in reflected light. vPPG signals can be used to monitor vital signs remotely in a noninvasive manner. We have used optical coherence tomography (OCT) imaging to find whether vPPG signals may result from skin optical properties being modulated by arterial transmural pressure propagation. Our results confirm this hypothesis and show that the attenuation coefficient of even the most external, non-vascularized skin, changes at the same frequency as the arterial pulsation.


title: Total Power Minimization: Joint Antenna Selection and Beamforming Design

abstract: In this paper, we consider the total power minimization problem when we have signal-to-interference-plus-noise ratio (SINR) constraints. The consumed power in the circuits depends on the number of active antennas, which can be modeled using zero-norm. Due to the difficulty of dealing with the non-convex zero-norm, we used the standard alternate weighted one-norm approach. We addressed the total power minimization for a narrowband system with and without per-antenna power constraints (PAPCs). We derived iterative closed-form expressions in both cases. Then we analysed the case when we have multiple bands operating at the same time. Analogous closed-form expressions are provided. Our simulation results show that significant gains can be obtained in terms of the total power required compared to standard methods that do not take into account the circuit power.


title: Towards Robust Deep Neural Networks

abstract: We investigate the topics of sensitivity and robustness in feedforward and convolutional neural networks. Combining energy landscape techniques developed in computational chemistry with tools drawn from formal methods, we produce empirical evidence indicating that networks corresponding to lower-lying minima in the optimization landscape of the learning objective tend to be more robust. The robustness estimate used is the inverse of a proposed sensitivity measure, which we define as the volume of an over-approximation of the reachable set of network outputs under all additive $l_{\infty}$-bounded perturbations on the input data. We present a novel loss function which includes a sensitivity term in addition to the traditional task-oriented and regularization terms. In our experiments on standard machine learning and computer vision datasets, we show that the proposed loss function leads to networks which reliably optimize the robustness measure as well as other related metrics of adversarial robustness without significant degradation in the classification error. Experimental results indicate that the proposed method outperforms state-of-the-art sensitivity-based learning approaches with regards to robustness to adversarial attacks. We also show that although the introduced framework does not explicitly enforce an adversarial loss, it achieves competitive overall performance relative to methods that do.

categories: stat.ML cond-mat.stat-mech cs.LG math.OC


title: Co-Toeplitz Quantization: A Simple Case

abstract: The author has introduced in a recent paper a new class of operators, called co-Toeplitz operators, with symbols in a co-algebra. This is the categorical dual to Toeplitz operators which have symbols in an algebra. The mapping from a symbol to its co-Toeplitz operator gives a quantization scheme, called co-Toeplitz quantization. A new, quite simple particular case of co-Toeplitz quantization is introduced in this note. Examples are given in order to show some of its properties.


title: Maximal contact and symplectic structures

abstract: We introduce a procedure for gluing Weinstein domains along Weinstein subdomains. By gluing along flexible subdomains, we show that any finite collection of high-dimensional Weinstein domains with the same topology are Weinstein subdomains of a `maximal' Weinstein domain also with the same topology. As an application, we produce exotic cotangent bundles containing many closed regular Lagrangians that are formally Lagrangian isotopic but not Hamiltonian isotopic and also give a new construction of exotic Weinstein structures on Euclidean space. We describe a similar construction in the contact setting which we use to produce `maximal' contact structures and extend several existing results in low-dimensional contact geometry to high-dimensions. We prove that all contact manifolds have symplectic caps, introduce a general procedure for producing contact manifolds with many Weinstein fillings, and give a new proof of the existence of codimension two contact embeddings.


title: Cooperative Deep Reinforcement Learning for Multiple-Group NB-IoT Networks Optimization

abstract: NarrowBand-Internet of Things (NB-IoT) is an emerging cellular-based technology that offers a range of flexible configurations for massive IoT radio access from groups of devices with heterogeneous requirements. A configuration specifies the amount of radio resources allocated to each group of devices for random access and for data transmission. Assuming no knowledge of the traffic statistics, the problem is to determine, in an online fashion at each Transmission Time Interval (TTI), the configurations that maximizes the long-term average number of IoT devices that are able to both access and deliver data. Given the complexity of optimal algorithms, a Cooperative Multi-Agent Deep Neural Network based Q-learning (CMA-DQN) approach is developed, whereby each DQN agent independently control a configuration variable for each group. The DQN agents are cooperatively trained in the same environment based on feedback regarding transmission outcomes. CMA-DQN is seen to considerably outperform conventional heuristic approaches based on load estimation.


title: Low-shot Learning via Covariance-Preserving Adversarial Augmentation Networks

abstract: Deep neural networks suffer from over-fitting and catastrophic forgetting when trained with small data. One natural remedy for this problem is data augmentation, which has been recently shown to be effective. However, previous works either assume that intra-class variances can always be generalized to new classes, or employ naive generation methods to hallucinate finite examples without modeling their latent distributions. In this work, we propose Covariance-Preserving Adversarial Augmentation Networks to overcome existing limits of low-shot learning. Specifically, a novel Generative Adversarial Network is designed to model the latent distribution of each novel class given its related base counterparts. Since direct estimation of novel classes can be inductively biased, we explicitly preserve covariance information as the `variability' of base examples during the generation process. Empirical results show that our model can generate realistic yet diverse examples, leading to substantial improvements on the ImageNet benchmark over the state of the art.


title: Real-time Action Recognition with Dissimilarity-based Training of Specialized Module Networks

abstract: This paper addresses the problem of real-time action recognition in trimmed videos, for which deep neural networks have defined the state-of-the-art performance in the recent literature. For attaining higher recognition accuracies with efficient computations, researchers have addressed the various aspects of limitations in the recognition pipeline. This includes network architecture, the number of input streams (where additional streams augment the color information), the cost function to be optimized, in addition to others. The literature has always aimed, though, at assigning the adopted network (or networks, in case of multiple streams) the task of recognizing the whole number of action classes in the dataset at hand. We propose to train multiple specialized module networks instead. Each module is trained to recognize a subset of the action classes. Towards this goal, we present a dissimilarity-based optimized procedure for distributing the action classes over the modules, which can be trained simultaneously offline. On two standard datasets--UCF-101 and HMDB-51--the proposed method demonstrates a comparable performance, that is superior in some aspects, to the state-of-the-art, and that satisfies the real-time constraint. We achieved 72.5\% accuracy on the challenging HMDB-51 dataset. By assigning fewer and unalike classes to each module network, this research paves the way to benefit from light-weight architectures without compromising recognition accuracy.


title: Post-prognostics decision in Cyber-Physical Systems

abstract: Prognostics and Health Management (PHM) offers several benefits for predictive maintenance. It predicts the future behavior of a system as well as its Remaining Useful Life (RUL). This RUL is used to planned the maintenance operation to avoid the failure, the stop time and optimize the cost of the maintenance and failure. However, with the development of the industry the assets are nowadays distributed this is why the PHM needs to be developed using the new IT. In our work we propose a PHM solution based on Cyber physical system where the physical side is connected to the analyze process of the PHM which are developed in the cloud to be shared and to benefit of the cloud characteristics


title: Evolution of Photospheric Vector Magnetic Field Associated with Moving Flare Ribbons As Seen By GST

abstract: The photospheric response to solar flares, also known as coronal back reaction, is often observed as sudden flare-induced changes in vector magnetic field and sunspot motions. However, it remains obscure whether evolving flare ribbons, the flare signature closest to the photosphere, are accompanied by changes in vector magnetic field therein. Here we explore the relationship between the dynamics of flare ribbons in the chromosphere and variations of magnetic fields in the underlying photosphere, using high-resolution off-band H-alpha images and near-infrared vector magnetograms of the M6.5 flare on 2015 June 22 observed with the 1.6 m Goode Solar Telescope. We find that changes of photospheric fields occur at the arrival of the flare ribbon front, thus propagating analogously to flare ribbons. In general, the horizontal field increases and the field lines become more inclined to the surface. When ribbons sweep through regions that undergo a rotational motion, the fields transiently turn more vertical with decreased horizontal field and inclination angle, and then restore and/or become more horizontal than before the ribbon arrival. The ribbon propagation decelerates near the sunspot rotation center, where the vertical field becomes permanently enhanced. Similar magnetic field changes are discernible in magnetograms from the Helioseismic and Magnetic Imager (HMI), and an inward collapse of coronal magnetic fields is inferred from the time sequence of non-linear force-free field models extrapolated from HMI magnetograms. We conclude that photospheric fields respond nearly instantaneously to magnetic reconnection in the corona.


title: Upper bounds on the second largest prime factor of an odd perfect number

abstract: Acquaah and Konyagin showed that if $N$ is an odd perfect number where $N= p_1^{a_1}p_2^{a_2} \cdots p_k^{a_k}$ where $p_1 < p_2 \cdots < p_k$ then one must have $p_k < 3^{1/3}N^{1/3}$. Using methods similar to theirs, we show that $p_{k-1}< (2N)^{1/5}$ and that $p_{k-1}p_k < 6^{1/4}N^{1/2}.$ We also show that if $p_k$ and $p_{k-1}$ are close to each other than these bounds can be further strengthened.


title: Middle-Out Decoding

abstract: Despite being virtually ubiquitous, sequence-to-sequence models are challenged by their lack of diversity and inability to be externally controlled. In this paper, we speculate that a fundamental shortcoming of sequence generation models is that the decoding is done strictly from left-to-right, meaning that outputs values generated earlier have a profound effect on those generated later. To address this issue, we propose a novel middle-out decoder architecture that begins from an initial middle-word and simultaneously expands the sequence in both directions. To facilitate information flow and maintain consistent decoding, we introduce a dual self-attention mechanism that allows us to model complex dependencies between the outputs. We illustrate the performance of our model on the task of video captioning, as well as a synthetic sequence de-noising task. Our middle-out decoder achieves significant improvements on de-noising and competitive performance in the task of video captioning, while quantifiably improving the caption diversity. Furthermore, we perform a qualitative analysis that demonstrates our ability to effectively control the generation process of our decoder.


title: Temperature dependence of the solid-liquid interface free energy of Ni and Al from molecular dynamics simulation of nucleation

abstract: The temperature dependence of the solid-liquid interfacial free energy, {\gamma}, is investigated for Al and Ni at the undercooled temperature regime based on a recently developed persistent-embryo method. The atomistic description of the nucleus shape is obtained from molecular dynamics simulations. The computed {\gamma} shows a linear dependence on the temperature. The values of {\gamma} extrapolated to the melting temperature agree well with previous data obtained by the capillary fluctuation method. Using the temperature dependence of {\gamma}, we estimate the nucleation free energy barrier in a wide temperature range from the classical nucleation theory. The obtained data agree very well with the results from the brute-force molecular dynamics simulations.


title: Life at the edge: complexity and criticality in biological function

abstract: Why life is complex and --most importantly-- what is the origin of the over abundance of complexity in nature? This is a fundamental scientific question which, paraphrasing the late Per Bak, "is screaming to be answered but seldom is even being asked". In these lectures we review recent attempts across several scales to understand the origins of complex biological problems from the perspective of critical phenomena. To illustrate the approach three cases are discussed, namely the large scale brain dynamics, the characterisation of spontaneous fluctuations of proteins and the physiological complexity of the cell mitochondria network.


title: Gaussian Process Prior Variational Autoencoders

abstract: Variational autoencoders (VAE) are a powerful and widely-used class of models to learn complex data distributions in an unsupervised fashion. One important limitation of VAEs is the prior assumption that latent sample representations are independent and identically distributed. However, for many important datasets, such as time-series of images, this assumption is too strong: accounting for covariances between samples, such as those in time, can yield to a more appropriate model specification and improve performance in downstream tasks. In this work, we introduce a new model, the Gaussian Process (GP) Prior Variational Autoencoder (GPPVAE), to specifically address this issue. The GPPVAE aims to combine the power of VAEs with the ability to model correlations afforded by GP priors. To achieve efficient inference in this new class of models, we leverage structure in the covariance matrix, and introduce a new stochastic backpropagation strategy that allows for computing stochastic gradients in a distributed and low-memory fashion. We show that our method outperforms conditional VAEs (CVAEs) and an adaptation of standard VAEs in two image data applications.


title: Large triangle packings and Tuza's conjecture in sparse random graphs

abstract: The triangle packing number $\nu(G)$ of a graph $G$ is the maximum size of a set of edge-disjoint triangles in $G$. Tuza conjectured that in any graph $G$ there exists a set of at most $2\nu(G)$ edges intersecting every triangle in $G$. We show that Tuza's conjecture holds in the random graph $G=G(n,m)$, when $m \le 0.2403n^{3/2}$ or $m\ge 2.1243n^{3/2}$. This is done by analyzing a greedy algorithm for finding large triangle packings in random graphs.


title: A Convex Duality Framework for GANs

abstract: Generative adversarial network (GAN) is a minimax game between a generator mimicking the true model and a discriminator distinguishing the samples produced by the generator from the real training samples. Given an unconstrained discriminator able to approximate any function, this game reduces to finding the generative model minimizing a divergence measure, e.g. the Jensen-Shannon (JS) divergence, to the data distribution. However, in practice the discriminator is constrained to be in a smaller class $\mathcal{F}$ such as neural nets. Then, a natural question is how the divergence minimization interpretation changes as we constrain $\mathcal{F}$. In this work, we address this question by developing a convex duality framework for analyzing GANs. For a convex set $\mathcal{F}$, this duality framework interprets the original GAN formulation as finding the generative model with minimum JS-divergence to the distributions penalized to match the moments of the data distribution, with the moments specified by the discriminators in $\mathcal{F}$. We show that this interpretation more generally holds for f-GAN and Wasserstein GAN. As a byproduct, we apply the duality framework to a hybrid of f-divergence and Wasserstein distance. Unlike the f-divergence, we prove that the proposed hybrid divergence changes continuously with the generative model, which suggests regularizing the discriminator's Lipschitz constant in f-GAN and vanilla GAN. We numerically evaluate the power of the suggested regularization schemes for improving GAN's training performance.


title: Deep Limits of Residual Neural Networks

abstract: Neural networks have been very successful in many applications; we often, however, lack a theoretical understanding of what the neural networks are actually learning. This problem emerges when trying to generalise to new data sets. The contribution of this paper is to show that, for the residual neural network model, the deep layer limit coincides with a parameter estimation problem for a nonlinear ordinary differential equation. In particular, whilst it is known that the residual neural network model is a discretisation of an ordinary differential equation, we show convergence in a variational sense. This implies that optimal parameters converge in the deep layer limit. This is a stronger statement than saying for a fixed parameter the residual neural network model converges (the latter does not in general imply the former). Our variational analysis provides a discrete-to-continuum $\Gamma$-convergence result for the objective function of the residual neural network training step to a variational problem constrained by a system of ordinary differential equations; this rigorously connects the discrete setting to a continuum problem.


title: Non-Commutative Integrability of the Grassmann Pentagram Map

abstract: The pentagram map is a discrete integrable system first introduced by Schwartz in 1992. It was proved to be intregable by Schwartz, Ovsienko, and Tabachnikov in 2010. Gekhtman, Shapiro, and Vainshtein studied Poisson geometry associated to certain networks embedded in a disc or annulus, and its relation to cluster algebras. Later, Gekhtman et al. and Tabachnikov reinterpreted the pentagram map in terms of these networks, and used the associated Poisson structures to give a new proof of integrability. In 2011, Mari Beffa and Felipe introduced a generalization of the pentagram map to certain Grassmannians, and proved it was integrable. We reinterpret this Grassmann pentagram map in terms of noncommutative algebra, in particular the double brackets of Van den bergh, and generalize the approach of Gekhtman et al. to establish a noncommutative version of integrability.


title: Transmission across a bilayer graphene region

abstract: The transmission across a graphene bilayer region is calculated for two different types of connections to monolayer leads. A transfer matrix algorithm based on a tight binding model is developed to obtain the ballistic transmission beyond linear response. The two configurations are found to behave similarly when no gate voltage is applied. For a finite gate voltage, both develop a conductance gap characteristic of a biased bilayer, but only one shows a pronounced conductance step at the gap edge. A gate voltage domain wall applied to the bilayer region renders the conductance of the two configurations similar. For a microstructure consisting of equally spaced domain walls, we find a high sensitivity to the domain size. This is attributed to the presence of topologically protected in-gap states localized at domain walls, which hybridize as the domain size becomes of the order of their confining scale. Our results show that transmission through a bilayer region can be manipulated by a gate voltage in ways not previously anticipated.


title: Another Enumeration of Caterpillar Trees

abstract: A caterpillar tree is a connected, acyclic, graph in which all vertices are either a member of a central path, or joined to that central path by a single edge. In other words, caterpillar trees are the class of trees which become path graphs after removing all leaves. In 1973, F. Harary and A.J. Schwenk provided two proofs found in [1] which show that the number of non-isomorphic caterpillars with N vertices is given by the formula $2^{N-4} + \ 2^{ \lfloor \frac{N - 4}{2}\rfloor}$, where $\lfloor \ \rfloor$ denotes the floor function. The first proof follows from a special case of an application of P\'{o}lya's Enumeration theorem on graphs with integer-weighted vertices. The second proof proceeds through an appropriate edge labelling of the caterpillars. The proof presented here owes much of its insight to the first two, but has the benefit of utilizing a natural labelling for the caterpillars. We will proceed by labelling the vertices of the caterpillars with integer-weights, followed by an application of the orbit-counting theorem.


title: Testing non-classicality with exact Wigner currents for an anharmonic quantum system

abstract: Phase-space features of the Wigner flow for an anharmonic quantum system driven by the harmonic oscillator potential modified by the addition of an inverse square (one-dimension Coulomb-like) contribution are analytically described in terms of Wigner functions and Wigner currents. Reporting about three correlated continuity equations which quantify the flux of quantum information in the phase-space, the non-classicality profile of such an anharmonic system can be consistently obtained in terms of the fluxes of {\em probability}, {\em purity} and {\em von Neumann-like entropy}. Considering that quantum fluctuations can be identified from distortions over the classical regime, they can be quantified through the above-mentioned information fluxes whenever some {\em classically bounded} volume of the phase-space is selected. Our results suggest that the Wigner flow approach works as a probe of quantumness and classicality for a large set of anharmonic quantum systems driven by quantum wells.


title: On buffered double autoregressive time series models

abstract: A buffered double autoregressive (BDAR) time series model is proposed in this paper to depict the buffering phenomenon of conditional mean and conditional variance in time series. To build this model, a novel flexible regime switching mechanism is introduced to modify the classical threshold time series model by capturing the stickiness of signal. Besides, considering the inadequacy of traditional models under the lack of information, a signal retrospection is run in this model to provide a more accurate judgment. Moreover, formal proofs suggest strict stationarity and geometric ergodicity of BDAR model under several sufficient conditions. A Gaussian quasi-maximum likelihood estimation (QMLE) is employed and the asymptotic distributions of its estimators are derived. It has been demonstrated that the estimated thresholds of the BDAR model are $n$-consistent, each of which converges weakly to a functional of a two-sided compound Poisson process. The remaining parameters are $\sqrt{n}$-consistent and asymptotically normal. Furthermore, a model selection criteria and its asymptotic property have been established. Simulation studies are constructed to evaluate the finite sample performance of QMLE and model selection criteria. Finally, an empirical analysis of Hang Seng Index (HSI) using BDAR model reveals the asymmetry of investors' preference over losses and gains as well as the asymmetry of volatility structure.


title: Sample Complexity for Nonlinear Dynamics

abstract: We consider the identification problems for nonlinear dynamical systems. An explicit sample complexity bound in terms of the number of data points required to recover the models accurately is derived. Our results extend recent sample complexity results for linear dynamics. Our approach for obtaining sample complexity bounds for nonlinear dynamics relies on a linear, albeit infinite dimensional, representation of nonlinear dynamics provided by Koopman and Perron-Frobenius operator. We exploit the linear property of these operators to derive the sample complexity bounds. Such complexity bounds will play a significant role in data-driven learning and control of nonlinear dynamics. Several numerical examples are provided to highlight our theory.