title: A Fuzzy-Rough based Binary Shuffled Frog Leaping Algorithm for Feature Selection

abstract: Feature selection and attribute reduction are crucial problems, and widely used techniques in the field of machine learning, data mining and pattern recognition to overcome the well-known phenomenon of the Curse of Dimensionality, by either selecting a subset of features or removing unrelated ones. This paper presents a new feature selection method that efficiently carries out attribute reduction, thereby selecting the most informative features of a dataset. It consists of two components: 1) a measure for feature subset evaluation, and 2) a search strategy. For the evaluation measure, we have employed the fuzzy-rough dependency degree (FRFDD) in the lower approximation-based fuzzy-rough feature selection (L-FRFS) due to its effectiveness in feature selection. As for the search strategy, a new version of a binary shuffled frog leaping algorithm is proposed (B-SFLA). The new feature selection method is obtained by hybridizing the B-SFLA with the FRDD. Non-parametric statistical tests are conducted to compare the proposed approach with several existing methods over twenty two datasets, including nine high dimensional and large ones, from the UCI repository. The experimental results demonstrate that the B-SFLA approach significantly outperforms other metaheuristic methods in terms of the number of selected features and the classification accuracy.


title: Coaction for Feynman integrals and diagrams

abstract: We propose a general coaction for families of integrals appearing in the evaluation of Feynman diagrams, such as multiple polylogarithms and generalized hypergeometric functions. We further conjecture a link between this coaction and graphical operations on Feynman diagrams. At one-loop order, there is a basis of integrals for which this correspondence is fully explicit. We discuss features and present examples of the diagrammatic coaction on two-loop integrals. We also present the coaction for the functions ${}_{p+1}F_p$ and Appell $F_1$.


title: Efficient closed domination in digraph products

abstract: A digraph $D$ is an efficient closed domination digraph if there exists a subset $S$ of $V(D)$ for which the closed out-neighborhoods centered in vertices of $S$ form a partition of $V(D)$. In this work we deal with efficient closed domination digraphs among several product of digraphs. We completely describe the efficient closed domination digraphs among lexicographic and strong products of digraphs. We characterize those direct products of digraphs that are efficient closed domination digraphs, where factors are either two cycles or two paths. Among Cartesian product of digraphs, we describe all such efficient closed domination digraphs such that they are a Cartesian product digraph either with a cycle or with a star.


title: Observables of low-lying supersymmetric vectorlike leptonic generations via loop corrections

abstract: A correlated analysis of observables arising from loop induced effects from a vectorlike generation is given. The observables include flavor changing radiative decays $\mu\to e \gamma, \tau\to \mu \gamma, \tau\to e \gamma$, electric dipole moments of the charged leptons $e,\mu, \tau$, and corrections to magnetic dipole moments of $g_\mu-2$ and $g_e-2$. In this work we give a full analysis of the corrections to these observables by taking into account both the supersymmetry loops as well as the exchange of a vectorlike leptonic generation. Thus the fermion mass matrix involves a $5\times 5$ mixing matrix while the scalar sector involves a $10\times 10$ mixing matrix including the CP violating phases from the vectorlike sector. The analysis is done under the constraint of the Higgs boson mass at the experimentally measured value. The loops considered include the exchange of $W$ and $Z$ bosons and of leptons and a mirror lepton, and the exchange of charginos and neutralinos, sleptons and mirror sleptons. The correction to the diphoton decay of the Higgs $h\to \gamma\gamma$ including the exchange of the vectorlike leptonic multiplet is also computed.


title: Annihilating-Ideal Graph of $C(X)$

abstract: In this article we study the annihilating-ideal graph of the ring $C(X)$. We have tried to associate the graph properties of $\mathbb{AG}(X)$, the ring properties of $C(X)$ and the topological properties of $X$. We have shown that $ X $ has an isolated point \ff $ \mathbb{R} $ is a direct summand of $ C(X) $ if and only if $ \mathbb{AG}(X) $ is not triangulated. Radius, girth, dominating number and clique number of the $\mathbb{AG}(X)$ are investigated. We have proved that $ c(X) \leqslant \mathrm{dt}(\mathbb{AG}(X)) \leqslant w(X) $ and $ \mathrm{clique} \mathbb{AG}(X) = \chi \mathbb{AG}(X) = c(X) $.


title: Markov branching processes with disasters: extinction, survival and duality to p-jump processes

abstract: A $p$-jump process is a piecewise deterministic Markov process with jumps by a factor of $p$. We prove a limit theorem for such processes on the unit interval. Via duality with respect to probability generating functions, we deduce limiting results for the survival probabilities of time-homogeneous branching processes with arbitrary offspring distributions, underlying binomial disasters. Extending this method, we obtain corresponding results for time-inhomogeneous birth-death processes underlying time-dependent binomial disasters and continuous state branching processes with $p$-jumps.


title: Ulrich bundles on three dimensional scrolls

abstract: In this paper we construct Ulrich bundles of low rank on three-dimensional scrolls (with respect to the tautological line bundle). We pay special attention to the four types of threefold scrolls in $\mathbb{P}^5$ which were classified in [Ott92].


title: $\Omega$ versus Graviphoton

abstract: I study the deformation of the topological string by $\bar\Omega$, the complex conjugate of the $\Omega$-deformation. Namely, I identify $\bar\Omega$ in terms of a physical state in the string spectrum and verify that the deformed Yang-Mills and ADHM actions are reproduced. This completes the study initiated recently [1] where we show that $\bar\Omega$ decouples from the one-loop topological amplitudes in heterotic string theory. Similarly to the N=2* deformation, I show that the quadratic terms in the effective action play a crucial role in obtaining the correct realisation of the full $\Omega$-deformation. Finally, I comment on the differences between the graviphoton and the $\Omega$-deformation in general and discuss possible $\bar\Omega$ remnants at the boundary of the string moduli space.


title: News Session-Based Recommendations using Deep Neural Networks

abstract: News recommender systems are aimed to personalize users experiences and help them to discover relevant articles from a large and dynamic search space. Therefore, news domain is a challenging scenario for recommendations, due to its sparse user profiling, fast growing number of items, accelerated item's value decay, and users preferences dynamic shift. Some promising results have been recently achieved by the usage of Deep Learning techniques on Recommender Systems, specially for item's feature extraction and for session-based recommendations with Recurrent Neural Networks. In this paper, it is proposed an instantiation of the CHAMELEON -- a Deep Learning Meta-Architecture for News Recommender Systems. This architecture is composed of two modules, the first responsible to learn news articles representations, based on their text and metadata, and the second module aimed to provide session-based recommendations using Recurrent Neural Networks. The recommendation task addressed in this work is next-item prediction for users sessions: "what is the next most likely article a user might read in a session?" Users sessions context is leveraged by the architecture to provide additional information in such extreme cold-start scenario of news recommendation. Users' behavior and item features are both merged in an hybrid recommendation approach. A temporal offline evaluation method is also proposed as a complementary contribution, for a more realistic evaluation of such task, considering dynamic factors that affect global readership interests like popularity, recency, and seasonality. Experiments with an extensive number of session-based recommendation methods were performed and the proposed instantiation of CHAMELEON meta-architecture obtained a significant relative improvement in top-n accuracy and ranking metrics (10% on Hit Rate and 13% on MRR) over the best benchmark methods.


title: Multiparty Dependent Session Types (Extended Abstract)

abstract: Programs are more distributed and concurrent today than ever before, and structural communications are at the core. Constructing and debugging such programs are hard due to the lack of formal specification/verification of concurrency. This work formalizes the first multiparty dependent session types as an expressive and practical type discipline for enforcing communication protocols. The type system is formulated in the setting of multi-threaded $\lambda$-calculus with inspirations from multirole logic, a generalization of classical logic we discovered earlier. We prove its soundness by a novel technique called deadlock-freeness reducibility. The soundness of the type system implies communication fidelity and absence of deadlock.


title: On the shortest distance between orbits and the longest common substring problem

abstract: In this paper, we study the behaviour of the shortest distance between orbits and show that under some rapidly mixing conditions, the decay of the shortest distance depends on the correlation dimension. For irrational rotations, we prove a different behaviour depending on the irrational exponent of the angle of the rotation. For random processes, this problem corresponds to the longest common substring problem. We extend the result of Arratia and Waterman on sequence matching to $\alpha$-mixing processes with exponential decay.


title: Optimal Gradient Checkpoint Search for Arbitrary Computation Graphs

abstract: Deep Neural Networks(DNNs) require huge GPU memory when training on modern image/video databases. Unfortunately, the GPU memory is physically finite, which limits the image resolutions and batch sizes that could be used in training for better DNN performance. Unlike solutions that require physically upgrade GPUs, the Gradient CheckPointing(GCP) training trades computation for more memory beyond existing GPU hardware. GCP only stores a subset of intermediate tensors, called Gradient Checkpoints (GCs), during forward. Then during backward, extra local forwards are conducted to compute the missing tensors. The total training memory cost becomes the sum of (1) the memory cost of the gradient checkpoints and (2) the maximum memory cost of local forwards. To achieve maximal memory cut-offs, one needs optimal algorithms to select GCs. Existing GCP approaches rely on either manual input of GCs or heuristics-based GC search on Linear Computation Graphs (LCGs), and cannot apply to Arbitrary Computation Graphs(ACGs). In this paper, we present theories and optimal algorithms on GC selection that, for the first time, are applicable to ACGs and achieve the maximal memory cut-offs. Extensive experiments show that our approach not only outperforms existing approaches (only applicable on LCGs), and is applicable to a vast family of LCG and ACG networks, such as Alexnet, VGG, ResNet, Densenet, Inception Net and highly complicated DNNs by Network Architecture Search. Our work enables GCP training on ACGs, and cuts off up-to 80% of training memory with a moderate time overhead (~30%-50%). Codes are available


title: Testing General Relativity with Black Hole-Pulsar Binaries

abstract: Binary pulsars allow us to carry out precision tests of gravity and have placed stringent bounds on a broad class of theories beyond general relativity. Current and future radio telescopes, such as FAST, SKA, and MeerKAT, may find a new astrophysical system, a pulsar orbiting around a black hole, which will provide us a new source for probing gravity. In this paper, we systematically study the prospects of testing general relativity with such black hole-pulsar binaries. We begin by finding a mapping between generic non-Einsteinian parameters in the orbital decay rate and theoretical constants in various modified theories of gravity and then summarize this mapping with a ready-to-use list. Theories we study here include scalar-tensor theories, varying $G$ theories, massive gravity theories, generic screening gravity and quadratic curvature-corrected theories. We next use simulated measurement accuracy of the orbital decay rate with FAST/SKA to derive projected upper bounds on the generic non-Einsteinian parameters. We find that such bounds from black hole-pulsars can be stronger than those from neutron star-pulsar and neutron star-white dwarf binaries by a few orders of magnitude when the correction enters at negative post-Newtonian orders. By mapping such bounds on generic parameters to those on various modified theories of gravity, we find that one can constrain the time variation in Newton's constant $G$ to be comparable and complementary to the current strongest bound from solar system experiments. We also study how well one can probe quadratic gravity from black hole quadrupole moment measurements of black hole-pulsars. We find that bounds on the parity-violating sector of quadratic gravity can be stronger than current bounds by six orders of magnitude. These results suggest that a new discovery of black hole-pulsars in the future will provide powerful ways to probe gravity further.


title: Solving NP-hard problems with bistable polaritonic networks

abstract: A lattice of locally bistable driven-dissipative cavity polaritons is found theoretically to effectively simulate the Ising model, also enabling an effective transverse field. We benchmark the system performance for spin glass problems, and study the scaling of the ground state energy deviation and success probability as a function of system size. As particular examples we consider NP-hard problems embedded in the Ising model, namely graph partitioning and the knapsack problem. We find that locally bistable polariton networks act as classical simulators for solving optimization problems, which can potentially present an improvement within the exponential complexity class.


title: Delay-Performance Tradeoffs in Causal Microphone Array Processing

abstract: In real-time listening enhancement applications, such as hearing aid signal processing, sounds must be processed with no more than a few milliseconds of delay to sound natural to the listener. Listening devices can achieve better performance with lower delay by using microphone arrays to filter acoustic signals in both space and time. Here, we analyze the tradeoff between delay and squared-error performance of causal multichannel Wiener filters for microphone array noise reduction. We compute exact expressions for the delay-error curves in two special cases and present experimental results from real-world microphone array recordings. We find that delay-performance characteristics are determined by both the spatial and temporal correlation structures of the signals.


title: Unsupervised machine learning for detection of phase transitions in off-lattice systems II. Applications

abstract: We outline how principal component analysis (PCA) can be applied to particle configuration data to detect a variety of phase transitions in off-lattice systems, both in and out of equilibrium. Specifically, we discuss its application to study 1) the nonequilibrium random organization (RandOrg) model that exhibits a phase transition from quiescent to steady-state behavior as a function of density, 2) orientationally and positionally driven equilibrium phase transitions for hard ellipses, and 3) compositionally driven demixing transitions in the non-additive binary Widom-Rowlinson mixture.


title: Unsupervised machine learning for detection of phase transitions in off-lattice systems I. Foundations

abstract: We demonstrate the utility of an unsupervised machine learning tool for the detection of phase transitions in off-lattice systems. We focus on the application of principal component analysis (PCA) to detect the freezing transitions of two-dimensional hard-disk and three-dimensional hard-sphere systems as well as liquid-gas phase separation in a patchy colloid model. As we demonstrate, PCA autonomously discovers order-parameter-like quantities that report on phase transitions, mitigating the need for a priori construction or identification of a suitable order parameter--thus streamlining the routine analysis of phase behavior. In a companion paper, we further develop the method established here to explore the detection of phase transitions in various model systems controlled by compositional demixing, liquid crystalline ordering, and non-equilibrium active forces.


title: Asymptotics in Spin-Boson type models

abstract: In this paper we investigate a family of models for a qubit interacting with a bosonic field. More precisely, we find asymptotic limits of the Hamiltonian as the strength of the interaction tends to infinity. The main result has two applications. First of all, we show that self-energy renormalisation schemes similar to that of the Nelson model will never give a physically interesting result. This is because any limit obtained through such a scheme would be independent of the qubit. Secondly, we find that exited states exist in the massive Spin-Boson models for sufficiently large interaction strengths. We are also able to compute the asymptotic limit of many physical quantities.


title: Experimental validation of the hybrid scattering model of installed jet noise

abstract: Jet installation causes jet noise to be amplified significantly at low frequencies and its physical mechanism must be understood to develop effective aircraft noise reduction strategies. A hybrid semi-empirical prediction model has recently been developed based on the instability-wave-scattering mechanism. However, its validity and accuracy remain to be tested. To do so, in this paper we carry out a systematic installed jet-noise experiment in the laboratory using a flat plate instead of an aircraft wing. We show that reducing $H$ (the separation distance between the flat plate and jet centreline) causes stronger low-frequency noise enhancement while resulting in little change to the noise shielding and enhancement at high frequencies. Decreasing $L$ (the axial distance between the jet exit plane and the trailing edge of the plate) results in reduced noise amplification at low frequencies and also weakens both the shielding and enhancement at high frequencies. Increasing the jet Mach number abates the installation effects. It is shown that the hybrid model developed in the earlier work agrees with experimental measurements and can capture the effects of varying $H$, $L$ and jet Mach number extremely well. It is concluded that the model captures the correct physics and can serve as an accurate and robust prediction tool. This new physical understanding provides insights into innovative strategies for suppressing installed jet noise.


title: Subsampled R\'enyi Differential Privacy and Analytical Moments Accountant

abstract: We study the problem of subsampling in differential privacy (DP), a question that is the centerpiece behind many successful differentially private machine learning algorithms. Specifically, we provide a tight upper bound on the R\'enyi Differential Privacy (RDP) (Mironov, 2017) parameters for algorithms that: (1) subsample the dataset, and then (2) applies a randomized mechanism M to the subsample, in terms of the RDP parameters of M and the subsampling probability parameter. Our results generalize the moments accounting technique, developed by Abadi et al. (2016) for the Gaussian mechanism, to any subsampled RDP mechanism.


title: Non-existence of ground states in the translation invariant Nelson model

abstract: In this paper we consider the massless translation invariant Nelson model with ultraviolet cutoff. It is proven that the fiber operators have no ground state if there is no infrared cutoff.


title: Towards Composable Bias Rating of AI Services

abstract: A new wave of decision-support systems are being built today using AI services that draw insights from data (like text and video) and incorporate them in human-in-the-loop assistance. However, just as we expect humans to be ethical, the same expectation needs to be met by automated systems that increasingly get delegated to act on their behalf. A very important aspect of an ethical behavior is to avoid (intended, perceived, or accidental) bias. Bias occurs when the data distribution is not representative enough of the natural phenomenon one wants to model and reason about. The possibly biased behavior of a service is hard to detect and handle if the AI service is merely being used and not developed from scratch, since the training data set is not available. In this situation, we envisage a 3rd party rating agency that is independent of the API producer or consumer and has its own set of biased and unbiased data, with customizable distributions. We propose a 2-step rating approach that generates bias ratings signifying whether the AI service is unbiased compensating, data-sensitive biased, or biased. The approach also works on composite services. We implement it in the context of text translation and report interesting results.


title: Combined cluster and atomic displacement expansion for solid solutions and magnetism

abstract: Finite temperature disordered solid solutions and magnetic materials are difficult to study directly using first principles calculations, due to the large unit cells and many independent samples that are required. In this work, we develop a combined cluster expansion and atomic displacement expansion, which we fit to first principles energies, forces, and stresses. We then use the expansion to calculate thermodynamic quantities at nearly first principles levels of accuracy. We demonstrate that by treating all the relevant degrees of freedom explicitly, we can achieve improved convergence properties as compared to a simple cluster expansion, and our model naturally includes both configurational and vibrational entropy. In addition, we can treat coupling between structural and chemical or magnetic degrees of freedom. As examples, we use our expansion to calculate properties of Si$_{1-x}$Ge$_x$, magnetic MnO, Al with vacancies, and Ba$_x$Sr$_{1-x}$TiO$_3$.


title: Correlations of multiplexed quantum ghost images and improvement of the quality of restored image

abstract: The currently used ghost image schemes traditionally involve two-mode entangled light states or incoherent radiation. Here, application of four-mode entangled light states is considered. It is shown that multiplexed ghost images (MGI) formed by four-mode entangled quantum light states have mutual spatial correlations determined by the 8th order field correlation functions. A special algorithm to calculate high-order correlations of Bose operators was developed. We also demonstrate that the accounting of MGI correlations allows us to improve the quality of the restored image of an object when processing MGI by measurement reduction method. Computer modelling of recovery of the image from MGI was carried out. It is established that in the considered example the signal-to-noise ratio of the reduced ghost image is $4.6$ times higher than the best signal-to-noise ratio for the ghost images themselves.


title: Multimodal Sensing and Interaction for a Robotic Hand Orthosis

abstract: Wearable robotic hand rehabilitation devices can allow greater freedom and flexibility than their workstation-like counterparts. However, the field is generally lacking effective methods by which the user can operate the device: such controls must be effective, intuitive, and robust to the wide range of possible impairment patterns. Even when focusing on a specific condition, such as stroke, the variety of encountered upper limb impairment patterns means that a single sensing modality, such as electromyography (EMG), might not be sufficient to enable controls for a broad range of users. To address this significant gap, we introduce a multimodal sensing and interaction paradigm for an active hand orthosis. In our proof-of-concept implementation, EMG is complemented by other sensing modalities, such as finger bend and contact pressure sensors. We propose multimodal interaction methods that utilize this sensory data as input, and show they can enable tasks for stroke survivors who exhibit different impairment patterns. We believe that robotic hand orthoses developed as multimodal sensory platforms with help address some of the key challenges in physical interaction with the user.


title: Implementation of Smart Contracts Using Hybrid Architectures with On- and Off-Blockchain Components

abstract: Recently, decentralised (on-blockchain) platforms have emerged to complement centralised (off-blockchain) platforms for the implementation of automated, digital (smart) contracts. However, neither alternative can individually satisfy the requirements of a large class of applications. On-blockchain platforms suffer from scalability, performance, transaction costs and other limitations. Off-blockchain platforms are afflicted by drawbacks due to their dependence on single trusted third parties. We argue that in several application areas, hybrid platforms composed from the integration of on- and off-blockchain platforms are more able to support smart contracts that deliver the desired quality of service (QoS). Hybrid architectures are largely unexplored. To help cover the gap, in this paper we discuss the implementation of smart contracts on hybrid architectures. As a proof of concept, we show how a smart contract can be split and executed partially on an off-blockchain contract compliance checker and partially on the Rinkeby Ethereum network. To test the solution, we expose it to sequences of contractual operations generated mechanically by a contract validator tool.


title: Monte Carlo simulations of a disordered superconductor-metal quantum phase transition

abstract: We investigate the quantum phase transitions of a disordered nanowire from superconducting to metallic behavior by employing extensive Monte Carlo simulations. To this end, we map the quantum action onto a (1+1)-dimensional classical XY model with long-range interactions in imaginary time. We then analyze the finite-size scaling behavior of the order parameter susceptibility, the correlation time, the superfluid density, and the compressibility. We find strong numerical evidence for the critical behavior to be of infinite-randomness type and to belong to the random transverse-field Ising universality class, as predicted by a recent strong disorder renormalization group calculation.


title: Battery Life-Cycle Optimization and Runtime Control for Commercial Buildings Demand Side Management: A New York City Case Study

abstract: In metropolitan areas populated with commercial buildings, electric power supply is stringent especially during business hours. Demand side management using battery is a promising solution to mitigate peak demands, however long payback time creates barriers for large scale adoption. In this paper, we have developed a design phase battery life-cycle cost assessment tool and a runtime controller for the building owners, taking into account the degradation of battery. In the design phase, perfect knowledge on building load profile is assumed to estimate ideal payback time. In runtime, stochastic programming and load predictions are applied to address the uncertainties in loads for producing optimal battery operation. For validation, we have performed numerical experiments using the real-life tariff model serves New York City, Zn/MnO2 battery, and state-of-the-art building simulation tool. Experimental results shows a small gap between design phase assessment and runtime control. To further examine the proposed methods, we have applied the same tariff model and performed numerical experiments on nine weather zones and three types of commercial buildings. On contrary to the common practice of shallow discharging battery for preventing phenomenal degradation, experimental results show promising payback time achieved by optimally deep discharge a battery.


title: Speech Separation Using Partially Asynchronous Microphone Arrays Without Resampling

abstract: We consider the problem of separating speech sources captured by multiple spatially separated devices, each of which has multiple microphones and samples its signals at a slightly different rate. Most asynchronous array processing methods rely on sample rate offset estimation and resampling, but these offsets can be difficult to estimate if the sources or microphones are moving. We propose a source separation method that does not require offset estimation or signal resampling. Instead, we divide the distributed array into several synchronous subarrays. All arrays are used jointly to estimate the time-varying signal statistics, and those statistics are used to design separate time-varying spatial filters in each array. We demonstrate the method for speech mixtures recorded on both stationary and moving microphone arrays.


title: The PROSPECT Reactor Antineutrino Experiment

abstract: The Precision Reactor Oscillation and Spectrum Experiment, PROSPECT, is designed to make both a precise measurement of the antineutrino spectrum from a highly-enriched uranium reactor and to probe eV-scale sterile neutrinos by searching for neutrino oscillations over meter-long baselines. PROSPECT utilizes a segmented $^6$Li-doped liquid scintillator detector for both efficient detection of reactor antineutrinos through the inverse beta decay reaction and excellent background discrimination. PROSPECT is a movable 4-ton antineutrino detector covering distances of 7m to 13m from the High Flux Isotope Reactor core. It will probe the best-fit point of the $\bar\nu_e$ disappearance experiments at 4$\sigma$ in 1 year and the favored regions of the sterile neutrino parameter space at more than 3$\sigma$ in 3 years. PROSPECT will test the origin of spectral deviations observed in recent $\theta_{13}$ experiments, search for sterile neutrinos, and address the hypothesis of sterile neutrinos as an explanation of the reactor anomaly. This paper describes the design, construction, and commissioning of PROSPECT and reports first data characterizing the performance of the PROSPECT antineutrino detector.


title: Universal Approximation with Quadratic Deep Networks

abstract: Recently, deep learning has achieved huge successes in many important applications. In our previous studies, we proposed quadratic/second-order neurons and deep quadratic neural networks. In a quadratic neuron, the inner product of a vector of data and the corresponding weights in a conventional neuron is replaced with a quadratic function. The resultant quadratic neuron enjoys an enhanced expressive capability over the conventional neuron. However, how quadratic neurons improve the expressing capability of a deep quadratic network has not been studied up to now, preferably in relation to that of a conventional neural network. Regarding this, we ask four basic questions in this paper: (1) for the one-hidden-layer network structure, is there any function that a quadratic network can approximate much more efficiently than a conventional network? (2) for the same multi-layer network structure, is there any function that can be expressed by a quadratic network but cannot be expressed with conventional neurons in the same structure? (3) Does a quadratic network give a new insight into universal approximation? (4) To approximate the same class of functions with the same error bound, is a quantized quadratic network able to enjoy a lower number of weights than a quantized conventional network? Our main contributions are the four interconnected theorems shedding light upon these four questions and demonstrating the merits of a quadratic network in terms of expressive efficiency, unique capability, compact architecture and computational capacity respectively.


title: Unconventional transport in low-density two-dimensional Rashba systems

abstract: Rashba spin-orbit coupling appears in 2D systems lacking inversion symmetry, and causes the spin-splitting of otherwise degenerate energy bands into an upper and lower helicity band. In this paper, we explore how impurity scattering affects transport in the ultra-low-density regime where electrons are confined to the lower helicity band. A previous study has investigated the conductivity in this regime using a treatment in the first Born approximation. In this work, we use the full T-matrix to uncover new features of the conductivity. We first compute the conductivity within a semiclassical Boltzmann framework and show that it exhibits an unconventional density dependence due to the unusual features of the group velocity in the single particle dispersion, as well as quantized plateaus as a function of the logarithm of the electron density. We support this with a calculation using the Kubo formula and find that these plateaus persist in the full quantum theory. We suggest that this quantization may be seen in a pump-probe experiment.


title: WeedMap: A large-scale semantic weed mapping framework using aerial multispectral imaging and deep neural network for precision farming

abstract: We present a novel weed segmentation and mapping framework that processes multispectral images obtained from an unmanned aerial vehicle (UAV) using a deep neural network (DNN). Most studies on crop/weed semantic segmentation only consider single images for processing and classification. Images taken by UAVs often cover only a few hundred square meters with either color only or color and near-infrared (NIR) channels. Computing a single large and accurate vegetation map (e.g., crop/weed) using a DNN is non-trivial due to difficulties arising from: (1) limited ground sample distances (GSDs) in high-altitude datasets, (2) sacrificed resolution resulting from downsampling high-fidelity images, and (3) multispectral image alignment. To address these issues, we adopt a stand sliding window approach that operates on only small portions of multispectral orthomosaic maps (tiles), which are channel-wise aligned and calibrated radiometrically across the entire map. We define the tile size to be the same as that of the DNN input to avoid resolution loss. Compared to our baseline model (i.e., SegNet with 3 channel RGB inputs) yielding an area under the curve (AUC) of [background=0.607, crop=0.681, weed=0.576], our proposed model with 9 input channels achieves [0.839, 0.863, 0.782]. Additionally, we provide an extensive analysis of 20 trained models, both qualitatively and quantitatively, in order to evaluate the effects of varying input channels and tunable network hyperparameters. Furthermore, we release a large sugar beet/weed aerial dataset with expertly guided annotations for further research in the fields of remote sensing, precision agriculture, and agricultural robotics.


title: Optimal 3D-Trajectory Design and Resource Allocation for Solar-Powered UAV Communication Systems

abstract: In this paper, we investigate the resource allocation algorithm design for multicarrier solar-powered unmanned aerial vehicle (UAV) communication systems. In particular, the UAV is powered by solar energy enabling sustainable communication services to multiple ground users. We study the joint design of the three-dimensional (3D) aerial trajectory and the wireless resource allocation for maximization of the system sum throughput over a given time period. As a performance benchmark, we first consider an offline resource allocation design assuming non-causal knowledge of the channel gains. The algorithm design is formulated as a mixed-integer non-convex optimization problem taking into account the aerodynamic power consumption, solar energy harvesting, a finite energy storage capacity, and the quality-of-service (QoS) requirements of the users. Despite the non-convexity of the optimization problem, we solve it optimally by applying monotonic optimization to obtain the optimal 3D-trajectory and the optimal power and subcarrier allocation policy. Subsequently, we focus on online algorithm design which only requires real-time and statistical knowledge of the channel gains. The optimal online resource allocation algorithm is motivated by the offline scheme and entails a high computational complexity. Hence, we also propose a low-complexity iterative suboptimal online scheme based on successive convex approximation. Our results unveil the tradeoff between solar energy harvesting and power-efficient communication. In particular, the solar-powered UAV first climbs up to a high altitude to harvest a sufficient amount of solar energy and then descents again to a lower altitude to reduce the path loss of the communication links to the users it serves.


title: First-principle prediction of the existence of C64-graphyne and its nitrogen and boron substitutions

abstract: By using of the first-principles calculations based on density functional theory, a novel monolayer planar structure named C64-graphyne is predicted. Tetratomic and hexatomic rings, as well as C-C triple bonds exist in this new stable structure with the lattice parameter of 9.291 {\AA}. The carbon hexatomic ring in C64-graphyne contains two quite distinct C-C bonds, which is known as cyclohexatriene. Its electronic band structures show a semiconductor nature with a narrow direct band gap of 0.35 eV. Interestingly, by substituting one nitrogen atom for a carbon atom in the hexatomic ring or a sp hybridized carbon atom on the carbon chain of C64-graphyne, two stable planar structures are obtained. Such doping with nitrogen atom induced metal properties in this monolayer structure. Further investigation shows that the alternating substitution of boron and nitrogen atoms for all carbon atoms in C64-graphyne also induced a new stable structure, the (BN)64 structure. The electronic studies proved its insulator property with the band gap of 4.08 eV.


title: An Ontology-Based Recommender System with an Application to the Star Trek Television Franchise

abstract: Collaborative filtering based recommender systems have proven to be extremely successful in settings where user preference data on items is abundant. However, collaborative filtering algorithms are hindered by their weakness against the item cold-start problem and general lack of interpretability. Ontology-based recommender systems exploit hierarchical organizations of users and items to enhance browsing, recommendation, and profile construction. While ontology-based approaches address the shortcomings of their collaborative filtering counterparts, ontological organizations of items can be difficult to obtain for items that mostly belong to the same category (e.g., television series episodes). In this paper, we present an ontology-based recommender system that integrates the knowledge represented in a large ontology of literary themes to produce fiction content recommendations. The main novelty of this work is an ontology-based method for computing similarities between items and its integration with the classical Item-KNN (K-nearest neighbors) algorithm. As a study case, we evaluated the proposed method against other approaches by performing the classical rating prediction task on a collection of Star Trek television series episodes in an item cold-start scenario. This transverse evaluation provides insights into the utility of different information resources and methods for the initial stages of recommender system development. We found our proposed method to be a convenient alternative to collaborative filtering approaches for collections of mostly similar items, particularly when other content-based approaches are not applicable or otherwise unavailable. Aside from the new methods, this paper contributes a testbed for future research and an online framework to collaboratively extend the ontology of literary themes to cover other narrative content.


title: Pressure Induced Compression of Flatbands in Twisted Bilayer Graphene

abstract: We investigate the bandwidth compression due to out of plane pressure of the moire flatbands near charge neutrality in twisted bilayer graphene for a continuous range of small rotation angles of up to $\sim2.5^{\circ}$. The flatband bandwidth minima angles are found to grow linearly with interlayer coupling {\omega} and decrease with Fermi velocity. Application of moderate pressure values of up to 2.5 GPa achievable through a hydraulic press should allow accessing a flatband for angles as large as $\sim 1.5$^{\circ}$ instead of $\sim 1 \circ$ at zero pressure. This reduction of the moir\'e pattern length for larger twist angle implies an increase of the effective Coulomb interaction scale per moire cell by about 50% and enhance roughly by a factor of $\sim 2$ the elastic energy that resists the commensuration strains due to the moire pattern. Our results suggest that application of pressure on twisted bilayer graphene nanodevices through a hydraulic press will notably facilitate the device preparation efforts required for exploring the ordered phases near magic angle flatbands.


title: Estimating the Maximum Intensities of Soft X-ray Flares Using Extreme Value Theory

abstract: Solar flares are one of the most energetic events in the solar system, their impact on Earth at ground level and its atmosphere remains under study. The repercussions of this phenomenon in our technological infrastructure includes radio blackouts and errors in geopositional and navigation systems that are considered natural hazards in ever more countries. Occurrence frequency and intensity of the most energetic solar flares are been taken into account in national programs for civil protection in order to reduce the risk and increase the resilience from Space Weather events. In this work we use the statistical theory of extreme values as well as other statistical methods in order to asses the magnitudes of the most extreme solar flare events expected to occur in a given period of time. We found that the data set under study presents a dual tail behaviour. Our results show that on average we can expect one solar flare greater than X23 each 25 years, that is to say, one such event each two solar cycles.


title: Sourcerer's Apprentice and the study of code snippet migration

abstract: On the worldwide web, not only are webpages connected but source code is too. Software development is becoming more accessible to everyone and the licensing for software remains complicated. We need to know if software licenses are being maintained properly throughout their reuse and evolution. This motivated the development of the Sourcerer's Apprentice, a webservice that helps track clone relicensing, because software typically employ software licenses to describe how their software may be used and adapted. But most developers do not have the legal expertise to sort out license conflicts. In this paper we put the Apprentice to work on empirical studies that demonstrate there is much sharing between StackOverflow code and Python modules and Python documentation that violates the licensing of the original Python modules and documentation: software snippets shared through StackOverflow are often being relicensed improperly to CC-BY-SA 3.0 without maintaining the appropriate attribution. We show that many snippets on StackOverflow are inappropriately relicensed by StackOverflow users, jeopardizing the status of the software built by companies and developers who reuse StackOverflow snippets.


title: Cosmic ray acceleration in accretion flows of galaxy clusters

abstract: We investigate acceleration of cosmic rays by shocks and accretion flows in galaxy clusters. Numerical results for spectra of accelerated particles and nonthermal emission are presented. It is shown that the acceleration of protons and nuclei in the nearby galaxy cluster Virgo can explain the observed spectra of ultra high energy cosmic rays.


title: The growth and development of living organisms from the thermodynamic point of view

abstract: The living organism is considered as an open system, whereas Prigogine's approach to the thermodynamics of such systems is used. The approach allows one to formulate the law of individual growth and development (ontogenesis) of the living organism, whereas it has taken into account that the development and functioning of the system are occurring under the special internal program. The thermodynamic equation of growth is followed a method of estimation of the specific entropy of organism. The theory is compared with experimental data, whereas estimates of the specific entropy of some species were done; it shows a reduction of specific entropy in the evolutionary row: yeast - insects - reptiles - birds.


title: Can graphene bilayers be the membrane mimetic materials?

abstract: Since the pioneering works of the founder of membrane mimetic chemistry Janos H. Fendler it is known that a number of atomic or molecular clusters and films (including nanoscale ones) are capable of mimicking the membrane functions. Membrane mimetic materials can be either soft matter or solid state materials. Conducting films (including those with magnetic properties) and semiconductors are also known to possess membrane mimetic properties. If we consider the agent exchange through the membrane in the operator form, the chemical composition of the membranes and their models, as well as the difference between the atomic and molecular clusters or layers become not so essential, and hence, membrane mimetic chemistry of nano- and mesostructures do not differ significantly within the agent-based approach. This invited review containing several parts reflects the main aspects of the author's report at the conference "Graphene: a molecule and 2D-crystal" (September 8-12, 2015, Novosibirsk, Russia) and considers various aspects of the similarity between the graphene nanostructures, membranes and bionic membrane-like nanomaterials.

categories: physics.bio-ph physics.atm-clus physics.chem-ph q-bio.MN


title: Methods of electron microscopy of biological and abiogenic structures in artificial gas atmospheres

abstract: This paper reviews opportunities for using electron microscopy in various gas atmospheres for the analysis and morpho-physiological modification of biological structures. The approaches that allow varying the gaseous phase content, as well as temperature, humidity, and pressure, are considered. The applicability of both kinetic and dynamic approaches to the tissue and bioinorganic structure manipulations is pointed out. The possibility of simulation of the beam-induced formation and disintegration of abiogenetic molecular structures is also mentioned as a particular case of the electron beam influence and treatment of the precursor medium in an artificial atmosphere.


title: Probability Calibration Trees

abstract: Obtaining accurate and well calibrated probability estimates from classifiers is useful in many applications, for example, when minimising the expected cost of classifications. Existing methods of calibrating probability estimates are applied globally, ignoring the potential for improvements by applying a more fine-grained model. We propose probability calibration trees, a modification of logistic model trees that identifies regions of the input space in which different probability calibration models are learned to improve performance. We compare probability calibration trees to two widely used calibration methods---isotonic regression and Platt scaling---and show that our method results in lower root mean squared error on average than both methods, for estimates produced by a variety of base learners.


title: Layer k-projection and unfolding electronic bands at interfaces

abstract: The k-projection method provides an approach to separate the contributions from different constituents in heterostructure systems and can act as an aid to connect the results of experiments and calculations. We show that the technique can be used to "unfold" the calculated electronic bands of interfaces and supercells, and provide local band structure by integrating the projected states over specified regions of space, a step that can be implemented efficiently using fast Fourier transforms. We apply the method to investigate the effects of interfaces in heterostructures consisting of a graphene bilayer on H-saturated SiC(0001), BAs monolayer on the ferromagnetic semiconductor CrI3, silicene on Ag(111), and to the Bi2Se3 surface. Our results reveal that the band structure of the graphene bilayer around the Dirac point is strongly dependent on the termination of SiC(0001): on the C face, the graphene is n doped and a gap of ~0.13 eV is opened, whereas, on the Si face, the graphene is essentially unchanged and neutral. We show that for BAs/CrI3, the magnetic proximity effect can effectively induce a spin splitting up to about 50 meV in BAs. For silicene/Ag(111), our calculations reproduce the angle-resolved photoemission spectroscopy results, including linearly dispersing bands at the edge of the first Brillouin zone of Ag(111); although these states result from the interaction between the silicene overlayer and the substrate, we demonstrate that they are not Dirac states.


title: Learning Stabilizable Dynamical Systems via Control Contraction Metrics

abstract: We propose a novel framework for learning stabilizable nonlinear dynamical systems for continuous control tasks in robotics. The key idea is to develop a new control-theoretic regularizer for dynamics fitting rooted in the notion of stabilizability, which guarantees that the learned system can be accompanied by a robust controller capable of stabilizing any open-loop trajectory that the system may generate. By leveraging tools from contraction theory, statistical learning, and convex optimization, we provide a general and tractable semi-supervised algorithm to learn stabilizable dynamics, which can be applied to complex underactuated systems. We validated the proposed algorithm on a simulated planar quadrotor system and observed notably improved trajectory generation and tracking performance with the control-theoretic regularized model over models learned using traditional regression techniques, especially when using a small number of demonstration examples. The results presented illustrate the need to infuse standard model-based reinforcement learning algorithms with concepts drawn from nonlinear control theory for improved reliability.


title: Automatic Detection and Diagnosis of Biased Online Experiments

abstract: We have seen a massive growth of online experiments at LinkedIn, and in industry at large. It is now more important than ever to create an intelligent A/B platform that can truly democratize A/B testing by allowing everyone to make quality decisions, regardless of their skillset. With the tremendous knowledge base created around experimentation, we are able to mine through historical data, and discover the most common causes for biased experiments. In this paper, we share four of such common causes, and how we build into our A/B testing platform the automatic detection and diagnosis of such root causes. These root causes range from design-imposed bias, self-selection bias, novelty effect and trigger-day effect. We will discuss in detail what each bias is and the scalable algorithm we developed to detect the bias. Surfacing up the existence and root cause of bias automatically for every experiment is an important milestone towards intelligent A/B testing.


title: Bragg Coherent Modulation Imaging: Strain- and Defect- Sensitive Single Views of Extended Samples

abstract: Nanoscale heterogeneity (including size, shape, strain, and defects) significantly impacts material properties and how they function. Bragg coherent x-ray imaging methods have emerged as a powerful tool to investigate, in three-dimensional detail, the local material response to external stimuli in reactive environments, thereby enabling explorations of the structure-defect-function relationship at the nanoscale. Although progress has been made in understanding this relationship, coherent imaging of extended samples is relatively slow (typically requiring many minutes) due to the experimental constraints required to solve the phase problem. Here, we develop Bragg coherent modulation imaging (BCMI), which uses a modulator to solve the phase problem thereby enabling fast, local imaging of an extended sample. Because a known modulator is essential to the technique, we first demonstrate experimentally that an unknown modulator structure can be recovered by using the exit wave diversity that exists in a standard Bragg coherent diffraction imaging (BCDI) experiment. We then show with simulations that a known modulator constraint is sufficient to solve the phase problem and enable a single view of an extended sample that is sensitive to defects and dislocations. Our results pave the way for BCMI investigations of strain and defect dynamics in complex extended crystals with temporal resolution limited by the exposure time of a single diffraction pattern.


title: Cognitive Techniques for Early Detection of Cybersecurity Events

abstract: The early detection of cybersecurity events such as attacks is challenging given the constantly evolving threat landscape. Even with advanced monitoring, sophisticated attackers can spend as many as 146 days in a system before being detected. This paper describes a novel, cognitive framework that assists a security analyst by exploiting the power of semantically rich knowledge representation and reasoning with machine learning techniques. Our Cognitive Cybersecurity system ingests information from textual sources, and various agents representing host and network-based sensors, and represents this information in a knowledge graph. This graph uses terms from an extended version of the Unified Cybersecurity Ontology. The system reasons over the knowledge graph to derive better actionable intelligence to security administrators, thus decreasing their cognitive load and increasing their confidence in the system. We have developed a proof of concept framework for our approach and demonstrate its capabilities using a custom-built ransomware instance that is similar to WannaCry.


title: CRUM: Checkpoint-Restart Support for CUDA's Unified Memory

abstract: Unified Virtual Memory (UVM) was recently introduced on recent NVIDIA GPUs. Through software and hardware support, UVM provides a coherent shared memory across the entire heterogeneous node, migrating data as appropriate. The older CUDA programming style is akin to older large-memory UNIX applications which used to directly load and unload memory segments. Newer CUDA programs have started taking advantage of UVM for the same reasons of superior programmability that UNIX applications long ago switched to assuming the presence of virtual memory. Therefore, checkpointing of UVM will become increasingly important, especially as NVIDIA CUDA continues to gain wider popularity: 87 of the top 500 supercomputers in the latest listings are GPU-accelerated, with a current trend of ten additional GPU-based supercomputers each year. A new scalable checkpointing mechanism, CRUM (Checkpoint-Restart for Unified Memory), is demonstrated for hybrid CUDA/MPI computations across multiple computer nodes. CRUM supports a fast, forked checkpointing, which mostly overlaps the CUDA computation with storage of the checkpoint image in stable storage. The runtime overhead of using CRUM is 6% on average, and the time for forked checkpointing is seen to be a factor of up to 40 times less than traditional, synchronous checkpointing.


title: Toward Multimodal Interaction in Scalable Visual Digital Evidence Visualization Using Computer Vision Techniques and ISS

abstract: Visualization requirements in Forensic Lucid have to do with different levels of case knowledge abstraction, representation, aggregation, as well as the operational aspects as the final long-term goal of this proposal. It encompasses anything from the finer detailed representation of hierarchical contexts to Forensic Lucid programs, to the documented evidence and its management, its linkage to programs, to evaluation, and to the management of GIPSY software networks. This includes an ability to arbitrarily switch between those views combined with usable multimodal interaction. The purpose is to determine how the findings can be applied to Forensic Lucid and investigation case management. It is also natural to want a convenient and usable evidence visualization, its semantic linkage and the reasoning machinery for it. Thus, we propose a scalable management, visualization, and evaluation of digital evidence using the modified interactive 3D documentary system - Illimitable Space System - (ISS) to represent, semantically link, and provide a usable interface to digital investigators that is navigable via different multimodal interaction techniques using Computer Vision techniques including gestures, as well as eye-gaze and audio.


title: The equivalence of several conjectures on independence of $\ell$

abstract: We consider several conjectures on the independence of $\ell$ of the \'etale cohomology of (singular, open) varieties over $\bar{\mathbf F}_p$. The main result is that independence of $\ell$ of the Betti numbers $h^i_{\text{c}}(X,\mathbf Q_\ell)$ for arbitrary varieties is equivalent to independence of $\ell$ of homological equivalence $\sim_{\text{hom},\ell}$ for cycles on smooth projective varieties. We give several other equivalent statements. As a surprising consequence, we prove that independence of $\ell$ of Betti numbers for smooth quasi-projective varieties implies the same result for arbitrary separated finite type $k$-schemes.


title: Dynamical Privacy in Distributed Computing -- Part II: PPSC Gossip Algorithms

abstract: In the first part of the paper, we have studied the computational privacy risks in distributed computing protocols against local or global dynamics eavesdroppers, and proposed a Privacy-Preserving-Summation-Consistent (PPSC) mechanism as a generic privacy encryption subroutine for consensus-based distributed computations. In this part of this paper, we show that the conventional deterministic and random gossip algorithms can be used to realize the PPSC mechanism over a given network. At each time step, a node is selected to interact with one of its neighbors via deterministic or random gossiping. Such node generates a random number as its new state, and sends the subtraction between its current state and that random number to the neighbor; then the neighbor updates its state by adding the received value to its current state. We establish concrete privacy-preservation conditions by proving the impossibility for the reconstruction of the network input from the output of the gossip-based PPSC mechanism against eavesdroppers with full network knowledge, and by showing that the PPSC mechanism can achieve differential privacy at arbitrary privacy levels. The convergence is characterized explicitly and analytically for both deterministic and randomized gossiping, which is essentially achieved in a finite number of steps. Additionally, we illustrate that the proposed algorithms can be easily made adaptive in real-world applications by making realtime trade-offs between resilience against node dropout or communication failure and privacy preservation capabilities.


title: Social Robots for People with Developmental Disabilities: A User Study on Design Features of a Graphical User Interface

abstract: Social robots, also known as service or assistant robots, have been developed to improve the quality of human life in recent years. The design of socially capable and intelligent robots can vary, depending on the target user groups. In this work, we assess the effect of social robots' roles, functions, and communication approaches in the context of a social agent providing service or entertainment to users with developmental disabilities. In this paper, we describe an exploratory study of interface design for a social robot that assists people suffering from developmental disabilities. We developed series of prototypes and tested one in a user study that included three residents with various function levels. This entire study had been recorded for the following qualitative data analysis. Results show that each design factor played a different role in delivering information and in increasing engagement. We also note that some of the fundamental design principles that would work for ordinary users did not apply to our target user group. We conclude that social robots could benefit our target users, and acknowledge that these robots were not suitable for certain scenarios based on the feedback from our users.


title: An inhomogeneous most likely path formalism for proton computed tomography

abstract: Purpose: Multiple Coulomb scattering (MCS) poses a challenge in proton CT (pCT) image reconstruction. The assumption of straight paths is replaced with Bayesian models of the most likely path (MLP). Current MLP-based pCT reconstruction approaches assume a water scattering environment. We propose an MLP formalism based on accurate determination of scattering moments in inhomogeneous media. Methods: Scattering power relative to water (RScP) was calculated for a range of human tissues and investigated against relative stopping power (RStP). Monte Carlo simulation was used to compare the new inhomogeneous MLP formalism to the water approach in a slab geometry and a human head phantom. An MLP-Spline-Hybrid method was investigated for improved computational efficiency. Results: A piecewise-linear correlation between RStP and RScP was shown, which may assist in iterative pCT reconstruction. The inhomogeneous formalism predicted Monte Carlo proton paths through a water cube with thick bone inserts to within 1.0 mm for beams ranging from 210 to 230 MeV incident energy. Improvement in accuracy over the conventional MLP ranged from 5% for a 230 MeV beam to 17% for 210 MeV. There was no noticeable gain in accuracy when predicting 200 MeV proton paths through a clinically relevant human head phantom. The MLP-Spline-Hybrid method reduced computation time by half while suffering negligible loss of accuracy. Conclusions: We have presented an MLP formalism that accounts for material composition. In most clinical cases a water scattering environment can be assumed, however in certain cases of significant heterogeneity the proposed algorithm may improve proton path estimation.


title: EagleEye: Attack-Agnostic Defense against Adversarial Inputs (Technical Report)

abstract: Deep neural networks (DNNs) are inherently vulnerable to adversarial inputs: such maliciously crafted samples trigger DNNs to misbehave, leading to detrimental consequences for DNN-powered systems. The fundamental challenges of mitigating adversarial inputs stem from their adaptive and variable nature. Existing solutions attempt to improve DNN resilience against specific attacks; yet, such static defenses can often be circumvented by adaptively engineered inputs or by new attack variants. Here, we present EagleEye, an attack-agnostic adversarial tampering analysis engine for DNN-powered systems. Our design exploits the {\em minimality principle} underlying many attacks: to maximize the attack's evasiveness, the adversary often seeks the minimum possible distortion to convert genuine inputs to adversarial ones. We show that this practice entails the distinct distributional properties of adversarial inputs in the input space. By leveraging such properties in a principled manner, EagleEye effectively discriminates adversarial inputs and even uncovers their correct classification outputs. Through extensive empirical evaluation using a range of benchmark datasets and DNN models, we validate EagleEye's efficacy. We further investigate the adversary's possible countermeasures, which implies a difficult dilemma for her: to evade EagleEye's detection, excessive distortion is necessary, thereby significantly reducing the attack's evasiveness regarding other detection mechanisms.


title: Finiteness of trivial solutions of factorial products yielding a factorial over number fields

abstract: We consider a Bertrand type estimate for primes splitting completely. As one of its applications, we show the finiteness of trivial solutions of Diophantine equation about the factorial function over number fields except for the case the rational number field.


title: Free boundary problems arising in the theory of maximal solutions of equations with exponential nonlinearities

abstract: We consider equations of the form $\Delta u +\lambda^2 V(x)e^{\,u}=\rho$ in various two dimensional settings. We assume that $V>0$ is a given function, $\lambda>0$ is a small parameter and $\rho=\mathcal O(1)$ or $\rho\to +\infty$ as $\lambda\to 0$. In a recent paper we prove the existence of the maximal solutions for a particular choice $V\equiv 1$, $\rho=0$ when the problem is posed in doubly connected domains under Dirichlet boundary conditions. We related the maximal solutions with a novel free boundary problem. The purpose of this note is to derive the corresponding free boundary problems in other settings. Solvability of such problems is, viewed formally, the necessary condition for the existence of the maximal solution.


title: 3BrY nucleation: from dimers to needle-like clusters

abstract: The nucleation of 1,3,5-Tris(4-bromophenyl)- benzene (3BrY) occurs through the formation of amorphous intermediates. However, a detailed experimental characterization of the structural and dynamic features of such early- stage precursors is prevented by their transient nature. In this work, we apply molecular mod- eling to investigate both nucleation and growth of 3BrY clusters from solution. Consistently with experimental findings, we observe that dis- ordered clusters consisting of 10-15 monomers can spontaneously emerge from solution. Such clusters are poorly ordered and can easily fluctuate in size and shape. When clusters grow to a larger size (up to 200 monomer units) they develop a markedly elongated morphology, reminiscent of the needle-like 3BrY crystals observed experimentally. The growth process is characterized by a continuous rearrangement of ordered and disordered domains accompa- nied by a dynamical exchange of molecules and oligomers with the solution, in a mechanism resembling the self-assembly of non-covalent supra-molecular fibers.


title: Maximal solution of the Liouville equation in doubly connected domains

abstract: In this paper we consider the Liouville equation $\Delta u +\lambda^2 e^{\,u}=0$ with Dirichlet boundary conditions in a two dimensional, doubly connected domain $\Omega$. We show that there exists a simple, closed curve $\gamma\subset \Omega$ such that for a sequence $\lambda_n\to 0$ and a sequence of solutions $u_{n}$ it holds $\frac{u_{n}}{\log\frac{1}{\lambda_n}}\to H$, where $H$ is a harmonic function in $\Omega\setminus\gamma$ and $\frac{\lambda_n^2}{\log\frac{1}{\lambda_n}}\int_\Omega e^{\,u_n}\,dx\to 8\pi c_\Omega$, where $c_\Omega$ is a constant depending on the conformal class of $\Omega$ only.


title: Simulation of quantum circuits by low-rank stabilizer decompositions

abstract: Recent work has explored using the stabilizer formalism to classically simulate quantum circuits containing a few non-Clifford gates. The computational cost of such methods is directly related to the notion of stabilizer rank, which for a pure state $\psi$ is defined to be the smallest integer $\chi$ such that $\psi$ is a superposition of $\chi$ stabilizer states. Here we develop a comprehensive mathematical theory of the stabilizer rank and the related approximate stabilizer rank. We also present a suite of classical simulation algorithms with broader applicability and significantly improved performance over the previous state-of-the-art. A new feature is the capability to simulate circuits composed of Clifford gates and arbitrary diagonal gates, extending the reach of a previous algorithm specialized to the Clifford+T gate set. We implemented the new simulation methods and used them to simulate quantum algorithms with 40-50 qubits and over 60 non-Clifford gates, without resorting to high-performance computers. We report a simulation of the Quantum Approximate Optimization Algorithm in which we process superpositions of $\chi\sim10^6$ stabilizer states and sample from the full n-bit output distribution, improving on previous simulations which used $\sim 10^3$ stabilizer states and sampled only from single-qubit marginals. We also simulated instances of the Hidden Shift algorithm with circuits including up to 64 T gates or 16 CCZ gates; these simulations showcase the performance gains available by optimizing the decomposition of a circuit's non-Clifford components.


title: Recurrent extensions of real-valued self-similar Markov processes

abstract: Let $X=(X_t, t\geq 0)$ be a self-similar Markov process taking values in $\mathbb{R}$ such that the state 0 is a trap. In this paper, we present a necessary and sufficient condition for the existence of a self-similar recurrent extension of $X$ that leaves 0 continuously. The condition is expressed in terms of the associated Markov additive process via the Lamperti-Kiu representation. Our results extend those of Fitzsimmons (2006) and Rivero (2005, 2007) where the existence and uniqueness of a recurrent extension for positive self similar Markov processes were treated. In particular, we describe the recurrent extension of a stable L\'evy process which to the best of our knowledge has not been studied before.


title: FMCode: A 3D In-the-Air Finger Motion Based User Login Framework for Gesture Interface

abstract: Applications using gesture-based human-computer interface require a new user login method with gestures because it does not have a traditional input method to type a password. However, due to various challenges, existing gesture-based authentication systems are generally considered too weak to be useful in practice. In this paper, we propose a unified user login framework using 3D in-air-handwriting, called FMCode. We define new types of features critical to distinguish legitimate users from attackers and utilize Support Vector Machine (SVM) for user authentication. The features and data-driven models are specially designed to accommodate minor behavior variations that existing gesture authentication methods neglect. In addition, we use deep neural network approaches to efficiently identify the user based on his or her in-air-handwriting, which avoids expansive account database search methods employed by existing work. On a dataset collected by us with over 100 users, our prototype system achieves 0.1% and 0.5% best Equal Error Rate (EER) for user authentication, as well as 96.7% and 94.3% accuracy for user identification, using two types of gesture input devices. Compared to existing behavioral biometric systems using gesture and in-air-handwriting, our framework achieves the state-of-the-art performance. In addition, our experimental results show that FMCode is capable to defend against client-side spoofing attacks, and it performs persistently in the long run. These results and discoveries pave the way to practical usage of gesture-based user login over the gesture interface.


title: A Theory of Dichotomous Valuation with Applications to Variable Selection

abstract: An econometric or statistical model may undergo a marginal gain if we admit a new variable to the model, and a marginal loss if we remove an existing variable from the model. Assuming equality of opportunity among all candidate variables, we derive a valuation framework by the expected marginal gain and marginal loss in all potential modeling scenarios. However, marginal gain and loss are not symmetric; thus, we introduce three unbiased solutions. When used in variable selection, our new approaches significantly outperform several popular methods used in practice. The results also explore some novel traits of the Shapley value.

categories: stat.ML cs.LG econ.EM math.OC math.ST stat.TH


title: Interaction of domain walls and vortices in the two-dimensional O(2) and O(3) principal chiral models

abstract: Using extensive Monte Carlo simulations, we investigate the critical properties of domain walls, vortices and $\mathbb{Z}_2$ vortices in the Ising-$O(2)$ and Ising-$O(3)\otimes O(2)$ models. We have consider the nontrivial case when disorder in the Ising order parameter induces disorder in the continuous parameter. Such a situation arises when a domain wall becomes opaque for continuous parameter correlations. We find that in this case the vortex density at the BKT transition (or crossover) point turns out to be non-universal, while the wall density at the Ising transition remains universal, i.e. in agreement with the Ising model. An important part of this study is the numerical measurement of defect-defect correlators. We find that the wall-vortex correlator tends to zero in the thermodynamic limit at the Ising point, which explains the universality of the wall density. A possible multicritical behavior of the models is also discussed.


title: Re-examining the premise of isobaric collisions and a novel method to measure the chiral magnetic effect

abstract: In this proceeding we will show that the expectations of the isobaric $^{96}_{44}\mathrm{Ru}+^{96}_{44}\mathrm{Ru}$ and $^{96}_{40}\mathrm{Zr}+^{96}_{40}\mathrm{Zr}$ collisions on chiral magnetic effect (CME) search may not hold as originally anticipated due to large uncertainties in the isobaric nuclear structures. We demonstrate this using Woods-Saxon densities and the proton and neutron densities calculated by the density functional theory. Furthermore, a novel method is proposed to gauge background and possible CME contributions in the same system, intrinsically better than the isobaric collisions of two different systems. We illustrate the method with Monte Carlo Glauber and AMPT (A Multi-Phase Transport) simulations.


title: A time-space Hausdorff fractal model for non-Fickian transport in porous media

abstract: This paper presents a time-space Hausdorff derivative model for depicting solute transport in aquifers or water flow in heterogeneous porous media. In this model, the time and space Hausdorff derivatives are defined on non-Euclidean fractal metrics with power law scaling transform which, respectively, connect the temporal and spatial complexity during transport. The Hausdorff derivative model can be transformed to an advection-dispersion equation with time- and space-dependent dispersion and convection coefficients. This model is a fractal partial differential equation (PDE) defined on a fractal space and differs from the fractional PDE which is derived for non-local transport of particles on a non-fractal Euclidean space. As an example of applications of this model, an explicit solution with a constant diffusion coefficient and flow velocity subject to an instantaneous source is derived and fitted to the breakthrough curves of tritium as a tracer in porous media. These results are compared with those of a scale-dependent dispersion model and a time-scale dependent dispersion model. Overall, it is found that the fractal PDE based on the Hausdorff derivatives better captures the early arrival and heavy tail in the scaled breakthrough curves for variable transport distances. The estimated parameters in the fractal Hausrdorff model represent clear physical mechanisms such as linear relationships between the orders of Hausdorff derivatives and the transport distance. The mathematical formulation is applicable to both solute transport and water flow in porous media.


title: Conditional Decoupling of Quantum Information

abstract: Insights from quantum information theory show that correlation measures based on quantum entropy are fundamental tools that reveal the entanglement structure of multipartite states. In that spirit, Groisman, Popescu, and Winter [Physical Review A 72, 032317 (2005)] showed that the quantum mutual information $I(A;B)$ quantifies the minimal rate of noise needed to erase the correlations in a bipartite state of quantum systems $AB$. Here, we investigate correlations in tripartite systems $ABE$. In particular, we are interested in the minimal rate of noise needed to apply to the systems $AE$ in order to erase the correlations between $A$ and $B$ given the information in system $E$, in such a way that there is only negligible disturbance on the marginal $BE$. We present two such models of conditional decoupling, called deconstruction and conditional erasure cost of tripartite states $ABE$. Our main result is that both are equal to the conditional quantum mutual information $I(A;B|E)$ -- establishing it as an operational measure for tripartite quantum correlations.

categories: quant-ph cond-mat.stat-mech cs.IT hep-th math.IT


title: Multi-modal Cycle-consistent Generalized Zero-Shot Learning

abstract: In generalized zero shot learning (GZSL), the set of classes are split into seen and unseen classes, where training relies on the semantic features of the seen and unseen classes and the visual representations of only the seen classes, while testing uses the visual representations of the seen and unseen classes. Current methods address GZSL by learning a transformation from the visual to the semantic space, exploring the assumption that the distribution of classes in the semantic and visual spaces is relatively similar. Such methods tend to transform unseen testing visual representations into one of the seen classes' semantic features instead of the semantic features of the correct unseen class, resulting in low accuracy GZSL classification. Recently, generative adversarial networks (GAN) have been explored to synthesize visual representations of the unseen classes from their semantic features - the synthesized representations of the seen and unseen classes are then used to train the GZSL classifier. This approach has been shown to boost GZSL classification accuracy, however, there is no guarantee that synthetic visual representations can generate back their semantic feature in a multi-modal cycle-consistent manner. This constraint can result in synthetic visual representations that do not represent well their semantic features. In this paper, we propose the use of such constraint based on a new regularization for the GAN training that forces the generated visual features to reconstruct their original semantic features. Once our model is trained with this multi-modal cycle-consistent semantic compatibility, we can then synthesize more representative visual representations for the seen and, more importantly, for the unseen classes. Our proposed approach shows the best GZSL classification results in the field in several publicly available datasets.


title: Relativistic non-dipole effects in strong-field atomic ionization at moderate intensities

abstract: We present a detailed experimental and theoretical study on the relativistic non-dipole effects in strong-field atomic ionisation by near-infrared linearly-polarised few-cycle laser pulses in the intensity range 1014 -1015 W/cm2. We record high-resolution photoelectron momentum distributions of argon using a reaction microscope and compare our measurements with a truly ab-initio fully relativistic 3D model based on the time-dependent Dirac equation. We observe counter-intuitive peak shifts of the transverse electron momentum distribution in the direction opposite to that of laser propagation as a function of laser intensity and demonstrate an excellent agreement between experimental results and theoretical predictions.


title: Mini-Tracker concepts for the SALT transient follow-up program

abstract: The MeerKAT radio telescope array, the Large Synoptic Survey Telescope (LSST), and eventually the Square Kilometer Array (SKA) will usher in a remarkable new era in astronomy, with thousands of transients being discovered and transmitted to the astronomical community in near-real-time each night. Immediate spectroscopic follow-up will be critical to understanding their early-time physics - a task to which the Southern African Large Telescope (SALT) is uniquely suited, given its southerly latitude and the 14-degree-diameter uncorrected field (patrol area) of its 10-m spherical primary mirror. A new telescope configuration is envisioned, incorporating multiple mini-trackers that range around a much larger patrol area of 35 degrees in diameter. Each mini-tracker is equipped with a small spherical aberration corrector feeding an efficient, low resolution spectrograph to perform contemporaneous follow-up observations.


title: Collimated GeV attosecond electron-positron bunches from a plasma channel driven by 10 PW lasers

abstract: High-energy positrons and bright {\gamma}-ray sources are unique both for fundamental research and practical applications. However, GeV electron-positron pair jets and {\gamma}-ray flashes are still hardly produced in laboratories. Here we demonstrate that, by irradiating two 10 PW-scale laser pulses onto a near-critical density plasma channel, highly-directional GeV electron-positron pairs and bright {\gamma}-ray beams can be efficiently generated. Three-dimensional particle-in-cell simulations show that GeV positron jets show high density (8*10^21/cm^3), attosecond duration (400 as) and a divergence angle of 14{\deg}. Additionally, ultrabright (2*10^25 photons/s/mm2/mrad2/0.1%BW) collimated attosecond (370 as) {\gamma}-ray flashes with a laser energy conversion efficiency of 5.6% are emitted. Once realized in experiment, it may open up new possibilities for a wide variety of applications.


title: Quark Wigner distribution of the pion meson in light-cone quark model

abstract: We study the quark Wigner distributions of the pion meson which are the phase-space distributions containing the most general one-parton information.Using the pion wavefunctions deduced from a light-cone quark model, we calculate the Wigner distributions of the unpolarized, longitudinally polarized and transversely polarized quarks inside the pion meson via the Fock-state overlap representation, respectively.We present the numerical results of the transverse Wigner distributions in which the longitudinal momentum fraction is integrated out. The mixed Wigner distributions as functions of $b_y$ and $k_x$ are also presented.As a intermediate step, we provide the analytic results for the generalized transverse momentum dependent distributions of the pion meson in the same model.


title: Action Anticipation By Predicting Future Dynamic Images

abstract: Human action-anticipation methods predict what is the future action by observing only a few portion of an action in progress. This is critical for applications where computers have to react to human actions as early as possible such as autonomous driving, human-robotic interaction, assistive robotics among others. In this paper, we present a method for human action anticipation by predicting the most plausible future human motion. We represent human motion using Dynamic Images and make use of tailored loss functions to encourage a generative model to produce accurate future motion prediction. Our method outperforms the currently best performing action-anticipation methods by 4% on JHMDB-21, 5.2% on UT-Interaction and 5.1% on UCF 101-24 benchmarks.


title: Sleep-wake classification via quantifying heart rate variability by convolutional neural network

abstract: Fluctuations in heart rate are intimately tied to changes in the physiological state of the organism. We examine and exploit this relationship by classifying a human subject's wake/sleep status using his instantaneous heart rate (IHR) series. We use a convolutional neural network (CNN) to build features from the IHR series extracted from a whole-night electrocardiogram (ECG) and predict every 30 seconds whether the subject is awake or asleep. Our training database consists of 56 normal subjects, and we consider three different databases for validation; one is private, and two are public with different races and apnea severities. On our private database of 27 subjects, our accuracy, sensitivity, specificity, and AUC values for predicting the wake stage are 83.1%, 52.4%, 89.4%, and 0.83, respectively. Validation performance is similar on our two public databases. When we use the photoplethysmography instead of the ECG to obtain the IHR series, the performance is also comparable. A robustness check is carried out to confirm the obtained performance statistics. This result advocates for an effective and scalable method for recognizing changes in physiological state using non-invasive heart rate monitoring. The CNN model adaptively quantifies IHR fluctuation as well as its location in time and is suitable for differentiating between the wake and sleep stages.


title: Study of Polarization-Driven Shortening for Polar Codes Designed with the Gaussian Approximation

abstract: This paper presents a polarization-driven (PD) shortening technique for the design of rate-compatible polar codes. The proposed shortening strategy consists of reducing the generator matrix by relating its row index with the channel polarization index. We assume that the shortened bits are known by both the encoder and the decoder and employ successive cancellation (SC) for decoding the shortened codes constructed by the proposed PD technique. A performance analysis is then carried out based on the Spectrum Distance (SD). Simulations show that the proposed PD-based shortened polar codes outperform existing shortened polar codes.


title: Connectivity and Blockage Effects in Millimeter-Wave Air-To-Everything Networks

abstract: Millimeter-wave (mmWave) offers high data rate and bandwidth for air-to-everything (A2X) communications including air-to-air, air-to-ground, and air-to-tower. MmWave communication in the A2X network is sensitive to buildings blockage effects. In this paper, we propose an analytical framework to define and characterise the connectivity for an aerial access point (AAP) by jointly using stochastic geometry and random shape theory. The buildings are modelled as a Boolean line-segment process with fixed height. The blocking area for an arbitrary building is derived and minimized by optimizing the altitude of AAP. A lower bound on the connectivity probability is derived as a function of the altitude of AAP and different parameters of users and buildings including their densities, sizes, and heights. Our study yields guidelines on practical mmWave A2X networks deployment.


title: Observation of a Fine Structure in $e^+e^- \to$ hadrons Production at the Nucleon-antinucleon Threshold

abstract: A study of hadron production at the nucleon-antinucleon threshold has been performed with the CMD-3 detector at the VEPP-2000 $e^+e^-$ collider. A very fast rise with about 1 MeV width has been observed in the $e^+e^- \to p\bar p$ cross section. A sharp drop in the $e^+e^- \to 3(\pi^+\pi^-)$ cross section has been confirmed and found to have a less than 2 MeV width, in agreement with the observed fast rise of the $e^+e^- \to p\bar p$ cross section. For the first time a similar sharp drop is demonstrated in the $e^+e^- \to K^+K^-\pi^+\pi^-$ cross section. The behavior of the $e^+e^- \to 3(\pi^+\pi^-),~ K^+K^-\pi^+\pi^-$ cross sections cannot be explained by an interference of any resonance amplitude with continuum, therefore this phenomenon cannot be due to a narrow near-threshold resonance. No such structure has been observed in the $e^+e^- \to 2(\pi^+\pi^-)$ cross section.


title: New irrational polygons with Ehrhart-theoretic period collapse

abstract: In a recent paper, Cristofaro-Gardiner--Li--Stanley [CGLS15] constructed examples of irrational triangles whose Ehrhart functions (i.e. lattice-point count) are polynomials when restricted to positive integer dilation factors. This is very surprising because the Ehrhart functions of rational polygons are usually only quasi-polynomials. We demonstrate that most of their triangles can also be obtained by a simple cut-and-paste procedure that allows us to build new examples with more sides. Our examples might potentially have applications in the theory of symplectic embeddings.


title: Mitigation of Side-Effect Modulation in Optical OFDM VLC Systems

abstract: Side-effect modulation (SEM) has the potential to be a significant source of interference in future visible light communication (VLC) systems. SEM is a variation in the intensity of the light emitted by a luminaire and is usually a side-effect caused by the power supply used to drive the luminaires. For LED luminaires powered by a switched mode power supply, the SEM can be at much higher frequencies than that emitted by conventional incandescent or fluorescent lighting. It has been shown that the SEM caused by commercially available LED luminaires is often periodic and of low power. In this paper, we investigate the impact of typical forms of SEM on the performance of optical OFDM VLC systems; both ACO-OFDM and DCO-OFDM are considered. Our results show that even low levels of SEM power can significantly degrade the bit-error-rate performance. To solve this problem, an SEM mitigation scheme is described. The mitigation scheme is decision-directed and is based on estimating and subtracting the fundamental component of the SEM from the received signal. We describe two forms of the algorithm; one uses blind estimation while the other uses pilot-assisted estimation based on a training sequence. Decision errors, resulting in decision noise, limit the performance of the blind estimator even when estimation is based on very long signals. However, the pilot system can achieve more accurate estimations, thus better performance. Results are first presented for typical SEM waveforms for the case where the fundamental frequency of the SEM is known. The algorithms are then extended to include a frequency estimation step and the mitigation algorithm is shown also to be effective in this case.


title: Explicit computations of Fourier transforms of polyhedral cones

abstract: The Fourier transforms of polyhedral cones can be used, via Brion's theorem, to compute various geometric quantities of polytopes, such as volumes, moments, and lattice-point counts. We present a novel method of computing these conic Fourier transforms by polynomial interpolation. Given the fact that computing volumes of polytopes is #P-hard (Dyer--Frieze [DF88]), we cannot hope for fast algorithms in the general case. However, with extra assumptions on the combinatorics of the cone, we demonstrate it is possible to compute its Fourier transform efficiently.


title: Duality of (2,3,5)-distributions and Lagrangian cone structures

abstract: As was shown by a part of the authors, for a given $(2, 3, 5)$-distribution $D$ on a $5$-dimensional manifold $Y$, there is, locally, a Lagrangian cone structure $C$ on another $5$-dimensional manifold $X$ which consists of abnormal or singular paths of $(Y, D)$. We give a characterization of the class of Lagrangian cone structures corresponding to $(2, 3, 5)$-distributions. Thus we complete the duality between $(2, 3, 5)$-distributions and Lagrangian cone structures via pseudo-product structures of type $G_2$. A local example of non-flat perturbations of the global model of flat Lagrangian cone structure which corresponds to $(2,3,5)$-distributions is given.


title: Depth Estimation via Affinity Learned with Convolutional Spatial Propagation Network

abstract: Depth estimation from a single image is a fundamental problem in computer vision. In this paper, we propose a simple yet effective convolutional spatial propagation network (CSPN) to learn the affinity matrix for depth prediction. Specifically, we adopt an efficient linear propagation model, where the propagation is performed with a manner of recurrent convolutional operation, and the affinity among neighboring pixels is learned through a deep convolutional neural network (CNN). We apply the designed CSPN to two depth estimation tasks given a single image: (1) To refine the depth output from state-of-the-art (SOTA) existing methods; and (2) to convert sparse depth samples to a dense depth map by embedding the depth samples within the propagation procedure. The second task is inspired by the availability of LIDARs that provides sparse but accurate depth measurements. We experimented the proposed CSPN over two popular benchmarks for depth estimation, i.e. NYU v2 and KITTI, where we show that our proposed approach improves in not only quality (e.g., 30% more reduction in depth error), but also speed (e.g., 2 to 5 times faster) than prior SOTA methods.


title: Identifying Financial Institutions by Transaction Signatures

abstract: Financial data aggregators and Personal Financial Management (PFM) services are software products that help individuals manage personal finances by collecting information from multiple accounts at various Financial Institutes (FIs), presenting data in a coherent and concentrated way, and highlighting insights and suggestions. Money transfers consist of two sides and a direction. From the perspective of a financial data aggregator, an incoming transaction consists of a date, an amount, and a description string, but not the explicit identity of the sending FI. In this paper we investigate supervised learning based methods to infer the identity of the sending FI from the description string of a money transfer transaction, using a blend of traditional and RNN based NLP methods. Our approach is based on the observation that the textual description field associated with a transactions is subjected to various types of normalizations and standardizations, resulting in unique patterns that identify the issuer. We compare multiple methods using a large real-word dataset of over $10$ million transactions.


title: Analysis of Dynamic Pull-in Voltage of a Graphene MEMS Model

abstract: Bifurcation analysis of dynamic pull-in for a lumped mass model is presented. The restoring force of the spring is derived based on the nonlinear constitutive stress-strain law and the driving force of the mass attached to the spring is based on the electrostatic Coulomb force, respectively. The analysis is performed on the resulting nonlinear spring-mass equation with initial conditions. The necessary and sufficient conditions for the existence of periodic solutions are derived analytically and illustrated numerically. The conditions for bifurcation points on the parameters associated with the second-order elastic stiffness constant and the voltage are determined.


title: The Heun operator of Hahn type

abstract: The Heun-Hahn operator on the uniform grid is defined. This operator is shown to map polynomials of degree $n$ to polynomials of degree $n+1$, to be tridiagonal in bases made out of either Pochhammer or Hahn polynomials and to be bilinear in the operators of the Hahn algebra. The extension of this algebra that includes the Heun-Hahn operator as generator is described. Biorthogonal rational functions on uniform grids are shown to be related to this framework.


title: Knots Connected by Wide Ribbons

abstract: A ribbon is, intuitively, a smooth mapping of an annulus $S^1 \times I$ in 3-space having constant width $\varepsilon$. This can be formalized as a triple $(x,\varepsilon, \mathbf{u})$ where $x$ is smooth curve in 3-space and $\mathbf{u}$ is a unit vector field based along $x$. In the 1960s and 1970s, G. Calugareanu, G. H. White, and F. B. Fuller proved relationships between the geometry and topology of thin ribbons, in particular the "Link = Twist + Writhe" theorem that has been applied to help understand properties of double-stranded DNA. Although ribbons of small width have been studied extensively, it appears that less is known about ribbons of large width whose images (even via a smooth map) can be singular or self-intersecting. Suppose $K$ is a smoothly embedded knot in $\mathbb{R}^3$. Given a regular parameterization $\mathbf{x}(s)$, and a smooth unit vector field $\mathbf{u}(s)$ based along $K$, we may define a ribbon of width $R$ associated to $\mathbf{x}$ and $\mathbf{u}$ as the set of all points $\mathbf{x}(s) + r\mathbf{u}(s)$, $r \in [0,R]$. For large $R$, these wide ribbons typically have self-intersections. In this paper, we analyze how the knot type of the outer ribbon edge $\mathbf{x}(s) + R\mathbf{u}(s)$ relates to that of the original knot $K$. We show that, generically, there is an eventual limiting knot type of the outer ribbon edge as $R$ gets arbitrary large. We prove that this eventual knot type is one of only finitely many possibilities which depend just on the vector field $\mathbf{u}$. However, the particular knot type within the finite set depends on the parameterized curves $\mathbf{x}(s)$, $\mathbf{u}(s)$, and their interactions. Finally, we show how to control the curves and their parameterizations so that given two knot types $K_1$ and $K_2$, we can find a smooth ribbon of constant width connecting curves of these two knot types.


title: Nucleon gravitational form factors from instantons: forces between quark and gluon subsystems

abstract: Using the instanton picture of the QCD vacuum we compute the nucleon $\bar c^Q(t)$ form factor of the quark part of the energy momentum tensor (EMT). This form factor describes the non-conservation of the quark part of EMT and contributes to the quark pressure distribution inside the nucleon. Also it can be interpreted in terms of forces between quark and gluon subsystems inside the nucleon. We show that this form factor is parametrically small in the instanton packing fraction. Numerically we obtain for the nucleon EMT a small value of $\bar c^Q(0)\simeq 1.4\cdot 10^{-2}$ at the low normalisation point of $\sim 0.4$ GeV$^2$. This smallness implies interesting physics picture - the forces between quark and gluon mechanical subsystems are smaller than the forces inside each subsystem. The forces from side of gluon subsystem squeeze the quark subsystem - they are compression forces. Additionally, the smallness of $\bar c^Q(t)$ might justify Teryaev's equipartition conjecture. We estimate that the contribution of $\bar c^Q (t)$ to the pressure distribution inside the nucleon is in the range of 1 -20 % relative to the contribution of the quark $D$-term.