title: A Kernelized Manifold Mapping to Diminish the Effect of Adversarial Perturbations

abstract: The linear and non-flexible nature of deep convolutional models makes them vulnerable to carefully crafted adversarial perturbations. To tackle this problem, we propose a non-linear radial basis convolutional feature mapping by learning a Mahalanobis-like distance function. Our method then maps the convolutional features onto a linearly well-separated manifold, which prevents small adversarial perturbations from forcing a sample to cross the decision boundary. We test the proposed method on three publicly available image classification and segmentation datasets namely, MNIST, ISBI ISIC 2017 skin lesion segmentation, and NIH Chest X-Ray-14. We evaluate the robustness of our method to different gradient (targeted and untargeted) and non-gradient based attacks and compare it to several non-gradient masking defense strategies. Our results demonstrate that the proposed method can increase the resilience of deep convolutional neural networks to adversarial perturbations without accuracy drop on clean data.


title: Designing reactive power control rules for smart inverters using support vector machines

abstract: Smart inverters have been advocated as a fast-responding mechanism for voltage regulation in distribution grids. Nevertheless, optimal inverter coordination can be computationally demanding, and preset local control rules are known to be subpar. Leveraging tools from machine learning, the design of customized inverter control rules is posed here as a multi-task learning problem. Each inverter control rule is modeled as a possibly nonlinear function of local and/or remote control inputs. Given the electric coupling, the function outputs interact to yield the feeder voltage profile. Using an approximate grid model, inverter rules are designed jointly to minimize a voltage deviation objective based on anticipated load and solar generation scenarios. Each control rule is described by a set of coefficients, one for each training scenario. To reduce the communication overhead between the grid operator and the inverters, we devise a voltage regulation objective that is shown to promote parsimonious descriptions for inverter control rules. Numerical tests using real-world data on a benchmark feeder demonstrate the advantages of the novel nonlinear rules and explore the trade-off between voltage regulation and sparsity in rule descriptions.


title: Non-Hermitian dynamics without dissipation in quantum systems

abstract: Models based on non-Hermitian Hamiltonians can exhibit a range of surprising and potentially useful phenomena. Physical realizations typically involve couplings to sources of incoherent gain and loss; this is problematic in quantum settings, because of the unavoidable fluctuations associated with this dissipation. Here, we present several routes for obtaining unconditional non-Hermitian dynamics in non-dissipative quantum systems. We exploit the fact that quadratic bosonic Hamiltonians that do not conserve particle number give rise to non-Hermitian dynamical matrices. We discuss the nature of these mappings from non-Hermitian to Hermitian Hamiltonians, and explore applications to quantum sensing, entanglement dynamics and topological band theory. The systems we discuss could be realized in a variety of photonic and phononic platforms using the ubiquitous resource of parametric driving.


title: Separated Edge-Soliton Meditated Dynamic Switching of Vortex Chirality and Polarity

abstract: Magnetic vortices are characterized by the senses of in-plane magnetization chirality and by the polarity of the vortex core. The electrical control of vortex polarity and chirality is highly demanded not only for fundamental understanding on spin dynamics in nano-disks under different circumstances, but also for technological applications, such as magnetic non-volatile memories and spin torque oscillators for neuromorphic computing. Here we report a novel approach that enables one to electrically control both the vortex chirality and polarity with low energy consumption. Thorough micromagnetic simulations, we show that in thin nano-disks of diameter larger than 160 nm, with the presence of current-induced Oersted field, the dynamic transformation of the edge solitons is able to efficiently switch both vortex chirality and polarity with low current under certain circumstances. We then developed an approach to directly write any of the four vortex states by electrical current pulses from a random state. We further investigated the switching phase diagram as a function of disk diameters. The results show that the switching process is highly size-dependent. As disk diameter is smaller than 160 nm, the switch of VC chirality and polarity always takes place at the same time, resulting in an unchanged handedness before and after switch. Furthermore, the critical switch current can be as low as 3 x 10^6 A/cm^2 , indicating a possible way for low current switch of vortex chirality in small disks.


title: Theoretical model of HD 163296 presently forming in-situ planets and comparison with the models of AS 209, HL Tau, and TW Hya

abstract: We fit an isothermal oscillatory density model to the disk of HD 163296 in which planets have presumably already formed and they are orbiting at least within the four observed dark gaps. This 156 AU large axisymmetric disk shows various physical properties comparable to those of AS 209, HL Tau, and TW Hya that we have modeled previously; but it compares best to AS 209. The disks of HD 163296 and AS 209 are comparable in size and they share similar values of the power-law index $k\approx 0$ (a radial density profile of the form $\rho(R)\propto R^{-1}$), the rotational parameter $\beta_0$ (to within a factor of 3); a relatively small inner core radius (although this parameter for HD 163296 is exceptionally small, $R_1\simeq 0.15$ AU, presumably due to unresolved planets in the inner 50 AU); the scale length $R_0$ and the Jeans gravitational frequency $\Omega_J$ (to within factors of 1.4); the equation of state ($c_0^2/\rho_0$) and the central density $\rho_0$ (to within factors of 2); and the core angular velocity $\Omega_0$ (to within a factor of 4.5). In the end, we compare all six nebular disks that we have modeled so far.


title: Kantorovich Duality and Optimal Transport Problems on Magnetic Graphs

abstract: We consider Lipschitz- and Arens-Eells-type function spaces constructed for magnetic graphs, which are adapted to this setting from the area of optimal transport on discrete spaces. After establishing the duality between these spaces, we prove a characterization of the extreme points of the unit ball in the magnetic Lipschitz space as well as a result identifying the magnetic Arens-Eells space as a quotient of the classical Arens-Eells space of an associated classical graph called the lift graph.


title: A Strongly Asymptotically Optimal Agent in General Environments

abstract: Reinforcement Learning agents are expected to eventually perform well. Typically, this takes the form of a guarantee about the asymptotic behavior of an algorithm given some assumptions about the environment. We present an algorithm for a policy whose value approaches the optimal value with probability 1 in all computable probabilistic environments, provided the agent has a bounded horizon. This is known as strong asymptotic optimality, and it was previously unknown whether it was possible for a policy to be strongly asymptotically optimal in the class of all computable probabilistic environments. Our agent, Inquisitive Reinforcement Learner (Inq), is more likely to explore the more it expects an exploratory action to reduce its uncertainty about which environment it is in, hence the term inquisitive. Exploring inquisitively is a strategy that can be applied generally; for more manageable environment classes, inquisitiveness is tractable. We conducted experiments in "grid-worlds" to compare the Inquisitive Reinforcement Learner to other weakly asymptotically optimal agents.


title: Structural and magnetic properties in sputtered iron oxide epitaxial thin films -- Magnetite Fe$_3$O$_4$ and epsilon ferrite e-Fe$_2$O$_3$

abstract: Epitaxial thin film fabrication of iron oxides including magnetite Fe3O4 and epsilon-ferrite epsilon-Fe2O3 with the potential for advancing electromagnetic devices has been investigated, which led to the first ever epsilon-ferrite epitaxial layer being synthesized in the conventional sputtering process. Concerning Fe3O4 (100) / MgO (100) films, a cube-on-cube epitaxial relationship and sharp rocking curves with FWHM of 50 - 350 arcsec were confirmed regardless of the small amount of Ge additions. Sputtering Ar gas pressure PAr heavily influenced their magnetic and transport properties. High PAr = 15 mTorr caused a high magnetization of 6.52 kG for the Ge added sample and the clear Verwey transition at 122 K for the non Ge addition case. Conversion electron Mossbauer spectroscopy (CEMS) measurements revealed that low PAr < 10 mTorr causes Fe/O off-stoichiometry on the oxidizing side for the non Ge addition case and the reductive side for the Ge addition case, respectively. Regarding the epsilon-Fe2O3 (001) / SrTiO3(111) epilayer synthesis, bilayer microstructure composed of an approximately 5nm thick initially grown epsilon-Fe2O3 (001) epilayer and subsequently grown e-Fe2O3 (001) epilayer was confirmed from cross-sectional TEM observations. The coexistence of magnetically hard and soft phases was confirmed from the magnetization measurements. As a possible application of the single nm thick epsilon-Fe2O3 layer, 4-resistive-state multiferroic tunnel junction (MFTJ) is considered.


title: Systematic Design of Decentralized Algorithms for Consensus Optimization

abstract: We propose a separation principle that enables a systematic way of designing decentralized algorithms used in consensus optimization. Specifically, we show that a decentralized optimization algorithm can be constructed by combining a non-decentralized base optimization algorithm and decentralized consensus tracking. The separation principle provides modularity in both the design and analysis of algorithms under an automated convergence analysis framework using integral quadratic constraints (IQCs). We show that consensus tracking can be incorporated into the IQC-based analysis. The workflow is illustrated through the design and analysis of a decentralized algorithm based on the alternating direction method of multipliers.


title: Mixed-Triggered Reliable Control for Singular Networked Cascade Control Systems with Randomly Occurring Cyber Attack

abstract: In this paper, the issue of mixed-triggered reliable dissipative control is investigated for singular networked cascade control systems (NCCSs) with actuator saturation and randomly occurring cyber attacks. In order to utilize the limited communication resources effectively, a more general mixed-triggered scheme is established which includes both schemes namely time-triggered and event-triggered in a single framework. In particular, two main factors are incorporated to the proposed singular NCCS model namely, actuator saturation and randomly occurring cyber attack, which is an important role to damage the overall network security. By employing Lyapunov-Krasovskii stability theory, a new set of sufficient conditions in terms of linear matrix inequalities (LMIs) is derived to guarantee the singular NCCSs to be admissible and strictly (Q,S,R)-dissipative. Subsequently, a power plant boiler-turbine system based on a numerical example is provided to demonstrate the effectiveness of the proposed control scheme.


title: Thermal and Electrical Properties of Hybrid Composites with Graphene and Boron Nitride Fillers

abstract: We report on the thermal and electrical properties of hybrid epoxy composites with graphene and boron nitride fillers. The thicknesses, lateral dimensions, and aspect ratios of each filler material were intentionally selected for geometric similarity to one another, in contrast to prior studies that utilized dissimilar filler geometries to achieve a "synergistic" effect. We demonstrate that the electrically-conductive graphene and electrically-insulating boron nitride fillers allow one to effectively engineer the thermal and electrical conductivities of their resulting composites. By varying the constituent fraction of boron nitride to graphene in a composite with ~44% total filler loading, one can tune the thermal conductivity enhancement from a factor of x15 to x35 and increase the electrical conductivity by many orders of magnitude. The obtained results are important for the development of next-generation thermal interface materials with controllable electrical properties necessary for applications requiring either electrical grounding or insulation.


title: Learning Modular Safe Policies in the Bandit Setting with Application to Adaptive Clinical Trials

abstract: The stochastic multi-armed bandit problem is a well-known model for studying the exploration-exploitation trade-off. It has significant possible applications in adaptive clinical trials, which allow for dynamic changes in the treatment allocation probabilities of patients. However, most bandit learning algorithms are designed with the goal of minimizing the expected regret. While this approach is useful in many areas, in clinical trials, it can be sensitive to outlier data, especially when the sample size is small. In this paper, we define and study a new robustness criterion for bandit problems. Specifically, we consider optimizing a function of the distribution of returns as a regret measure. This provides practitioners more flexibility to define an appropriate regret measure. The learning algorithm we propose to solve this type of problem is a modification of the BESA algorithm [Baransi et al., 2014], which considers a more general version of regret. We present a regret bound for our approach and evaluate it empirically both on synthetic problems as well as on a dataset from the clinical trial literature. Our approach compares favorably to a suite of standard bandit algorithms.


title: Prediction of Human Trajectory Following a Haptic Robotic Guide Using Recurrent Neural Networks

abstract: Social intelligence is an important requirement for enabling robots to collaborate with people. In particular, human path prediction is an essential capability for robots in that it prevents potential collision with a human and allows the robot to safely make larger movements. In this paper, we present a method for predicting the trajectory of a human who follows a haptic robotic guide without using sight, which is valuable for assistive robots that aid the visually impaired. We apply a deep learning method based on recurrent neural networks using multimodal data: (1) human trajectory, (2) movement of the robotic guide, (3) haptic input data measured from the physical interaction between the human and the robot, (4) human depth data. We collected actual human trajectory and multimodal response data through indoor experiments. Our model outperformed the baseline result while using only the robot data with the observed human trajectory, and it shows even better results when using additional haptic and depth data.


title: IVOA: Introspective Vision for Obstacle Avoidance

abstract: Vision, as an inexpensive yet information rich sensor, is commonly used for perception on autonomous mobile robots. Unfortunately, accurate vision-based perception requires a number of assumptions about the environment to hold -- some examples of such assumptions, depending on the perception algorithm at hand, include purely lambertian surfaces, texture-rich scenes, absence of aliasing features, and refractive surfaces. In this paper, we present an approach for introspective vision for obstacle avoidance (IVOA) -- by leveraging a supervisory sensor that is occasionally available, we detect failures of stereo vision-based perception from divergence in plans generated by vision and the supervisory sensor. By projecting the 3D coordinates where the plans agree and disagree onto the images used for vision-based perception, IVOA generates a training set of reliable and unreliable image patches for perception. We then use this training dataset to learn a model of which image patches are likely to cause failures of the vision-based perception algorithm. Using this model, IVOA is then able to predict whether the relevant image patches in the observed images are likely to cause failures due to vision (both false positives and false negatives). We empirically demonstrate with extensive real-world data from both indoor and outdoor environments, the ability of IVOA to accurately predict the failures of two distinct vision algorithms.


title: Similarity-based Random Survival Forest

abstract: Predicting time-to-event outcomes in large databases can be a challenging but important task. One example of this is in predicting the time to a clinical outcome for patients in intensive care units (ICUs), which helps to support critical medical treatment decisions. In this context, the time to an event of interest could be, for example, survival time or time to recovery from a disease/ailment observed within the ICU. The massive health datasets generated from the uptake of Electronic Health Records (EHRs) are quite heterogeneous as patients can be quite dissimilar in their relationship between the feature vector and the outcome, adding more noise than information to prediction. In this paper, we propose a modified random forest method for survival data that identifies similar cases in an attempt to improve accuracy for predicting time-to-event outcomes; this methodology can be applied in various settings, including with ICU databases. We also introduce an adaptation of our methodology in the case of dependent censoring. Our proposed method is demonstrated in the Medical Information Mart for Intensive Care (MIMIC-III) database, and, in addition, we present properties of our methodology through a comprehensive simulation study. Introducing similarity to the random survival forest method indeed provides improved predictive accuracy compared to random survival forest alone across the various analyses we undertook.


title: The generalized TAP free energy II

abstract: In a recent paper [14], we developed the generalized TAP approach for mixed $p$-spin models with Ising spins at positive temperature. Here we extend these results in two directions. We find a simplified representation for the energy of the generalized TAP states in terms of the Parisi measure of the model and, in particular, show that the energy of all states at a given distance from the origin is the same. Furthermore, we prove the analogues of the positive temperature results at zero temperature, which concern the ground-state energy and the organization of ground-state configurations in space.


title: Active Authentication using an Autoencoder regularized CNN-based One-Class Classifier

abstract: Active authentication refers to the process in which users are unobtrusively monitored and authenticated continuously throughout their interactions with mobile devices. Generally, an active authentication problem is modelled as a one class classification problem due to the unavailability of data from the impostor users. Normally, the enrolled user is considered as the target class (genuine) and the unauthorized users are considered as unknown classes (impostor). We propose a convolutional neural network (CNN) based approach for one class classification in which a zero centered Gaussian noise and an autoencoder are used to model the pseudo-negative class and to regularize the network to learn meaningful feature representations for one class data, respectively. The overall network is trained using a combination of the cross-entropy and the reconstruction error losses. A key feature of the proposed approach is that any pre-trained CNN can be used as the base network for one class classification. Effectiveness of the proposed framework is demonstrated using three publically available face-based active authentication datasets and it is shown that the proposed method achieves superior performance compared to the traditional one class classification methods. The source code is available at: github.com/otkupjnoz/oc-acnn.


title: A Fundamental Performance Limitation for Adversarial Classification

abstract: Despite the widespread use of machine learning algorithms to solve problems of technological, economic, and social relevance, provable guarantees on the performance of these data-driven algorithms are critically lacking, especially when the data originates from unreliable sources and is transmitted over unprotected and easily accessible channels. In this paper we take an important step to bridge this gap and formally show that, in a quest to optimize their accuracy, binary classification algorithms -- including those based on machine-learning techniques -- inevitably become more sensitive to adversarial manipulation of the data. Further, for a given class of algorithms with the same complexity (i.e., number of classification boundaries), the fundamental tradeoff curve between accuracy and sensitivity depends solely on the statistics of the data, and cannot be improved by tuning the algorithm.


title: An optimal result for global classical and bounded solutions in a two-dimensional Keller-Segel-Navier-Stokes system with sensitivity

abstract: This paper deals with a boundary-value problem for a coupled chemotaxis-Navier-Stokes system involving tensor-valued sensitivity with saturation $$\left\{ \begin{array}{l} n_t+u\cdot\nabla n=\Delta n-\nabla\cdot(nS(x,n,c)\nabla c),\quad x\in \Omega, t>0, c_t+u\cdot\nabla c=\Delta c-c+n,\quad x\in \Omega, t>0,\\ u_t+\kappa(u \cdot \nabla)u+\nabla P=\Delta u+n\nabla \phi,\quad x\in \Omega, t>0,\\ \nabla\cdot u=0,\quad x\in \Omega, t>0, \end{array}\right.$$ which describes chemotaxis-fluid interaction in cases when the evolution of the chemoattractant is essentially dominated by production through cells, where $\kappa\in \mathbb{R},\phi\in W^{2,\infty}(\Omega)$ and $S$ is a given function with values in $\mathbb{R}^{2\times2}$ which fulfills $$|S(x,n,c)| \leq C_S (1 + n)^{-\alpha}$$ with some $C _S > 0$ and $\alpha \geq 0.$ If $\alpha>0$ and $\Omega\subseteq \mathbb{R}^2$ is a {\bf bounded} domain with smooth boundary, then for all reasonably regular initial data, a corresponding initial-boundary value problem for $(KSNF)$ possesses a global classical solution which is bounded on $\Omega\times(0,\infty)$. This extends a recent result by Wang-Winkler-Xiang (Annali della Scuola Normale Superiore di Pisa-Classe di Scienze. XVIII, (2018), 2036--2145) which asserts global existence of bounded solutions under the constraint $\Omega\subseteq \mathbb{R}^2$ is a bounded {\bf convex domain} with smooth boundary. Moreover, we shall improve the result of Wang-Xiang (J. Diff. Eqns., 259(2015), 7578--7609), who proved the possibility of global and bounded, in the case that ${\bf\kappa\equiv0}$ and $\alpha>0$. In comparison to the result for the corresponding fluid-free system, the {\bf optimal condition} on the parameter $\alpha$ for both {\bf global existence} and {\bf boundedness} are obtained.


title: Observational Constraints on the Rastall gravity from Rotation Curves of Low Surface Brightness Galaxies

abstract: The Rastall gravity is a modification of Einstein's general relativity, in which the energy-momentum conservation is not satisfied and depends on the gradient of the Ricci curvature. It is in dispute whether the Rastall gravity is equivalent to the general relativity (GR). In this work, we constrain the theory using the rotation curves of Low Surface Brightness (LSB) spiral galaxies. Through fitting the rotation curves of LSB galaxies, we obtain the parameter $\beta$ of the Rastall gravity. The $\beta$ values of LSB galaxies satisfy Weak Energy Condition (WEC) and Strong Energy Condition(SEC). Combining the $\beta$ values of type Ia supernovae and gravitational lensing of elliptical galaxies on the Rastall gravity, we conclude that the Rastall gravity is equivalent to the general relativity.


title: Detection of latent heteroscedasticity and group-based regression effects in linear models via Bayesian model selection

abstract: Standard linear modeling approaches make potentially simplistic assumptions regarding the structure of categorical effects that may obfuscate more complex relationships governing data. For example, recent work focused on the two-way unreplicated layout has shown that hidden groupings among the levels of one categorical predictor frequently interact with the ungrouped factor. We extend the notion of a "latent grouping factor" to linear models in general. The proposed work allows researchers to determine whether an apparent grouping of the levels of a categorical predictor reveals a plausible hidden structure given the observed data. Specifically, we offer Bayesian model selection-based approaches to reveal latent group-based heteroscedasticity, regression effects, and/or interactions. Failure to account for such structures can produce misleading conclusions. Since the presence of latent group structures is frequently unknown a priori to the researcher, we use fractional Bayes factor methods and mixture $g$-priors to overcome lack of prior information.


title: Initial data for general relativistic simulations of multiple electrically charged black holes with linear and angular momenta

abstract: A general relativistic, stationary and axisymmetric black hole in a four-dimensional asymptotically-flat spacetime is fully determined by its mass, angular momentum and electric charge. The expectation that astrophysically relevant black holes do not posses charge has resulted in a limited number of investigations of moving and charged black holes in the dynamical, strong-field gravitational (and electromagnetic) regime, where numerical studies are necessary. Apart from having a theoretical interest, the advent of multimessenger astronomy with gravitational waves offers new ways to think about charged black holes. In this work, we initiate an exploration of charged binary black holes by generating valid initial data for general relativistic simulations of black hole systems that have generic electric charge, linear and angular momenta. We develop our initial data formalism within the framework of the conformal transverse-traceless (Bowen-York) technique using the puncture approach, and apply the theory of isolated horizons to attribute physical parameters (mass, charge and angular momentum) to each hole. We implemented our formalism in the case of a binary system by modifying the publicly available TwoPunctures and QuasiLocalMeasures codes. We demonstrate that our code can recover existing solutions and that it has excellent self-convergence properties for a generic configuration of two black holes.


title: Thermoelectric microscopy of magnetic skyrmions

abstract: The magnetic skyrmion is a nanoscale topological object characterized by the winding of the magnetic moments, appearing in magnetic materials with broken inversion symmetry. Because of its low current threshold for driving, the skyrmions have been intensely studied toward novel storage applications by using electron-beam, X-ray, and visible light microscopies. Here, we show that the skyrmions can be imaged via thermoelectric signals of spin-caloritronic phenomena in combination with focused heating. We applied a local temperature gradient to a CoFeB/Ta/W multilayer film with scanning a heating position and, by exploiting a Hall bar structure, mapped the magnitude as well as the direction of the resultant thermoelectric current distributions, which link to skyrmions' inner magnetic textures. This method also enables the observation of skyrmion dynamics under driving current pulses, being a useful imaging technique for the development of skyrmion devices.


title: Spatiotemporal Pyramid Network for Video Action Recognition

abstract: Two-stream convolutional networks have shown strong performance in video action recognition tasks. The key idea is to learn spatiotemporal features by fusing convolutional networks spatially and temporally. However, it remains unclear how to model the correlations between the spatial and temporal structures at multiple abstraction levels. First, the spatial stream tends to fail if two videos share similar backgrounds. Second, the temporal stream may be fooled if two actions resemble in short snippets, though appear to be distinct in the long term. We propose a novel spatiotemporal pyramid network to fuse the spatial and temporal features in a pyramid structure such that they can reinforce each other. From the architecture perspective, our network constitutes hierarchical fusion strategies which can be trained as a whole using a unified spatiotemporal loss. A series of ablation experiments support the importance of each fusion strategy. From the technical perspective, we introduce the spatiotemporal compact bilinear operator into video analysis tasks. This operator enables efficient training of bilinear fusion operations which can capture full interactions between the spatial and temporal features. Our final network achieves state-of-the-art results on standard video datasets.


title: SECNLP: A Survey of Embeddings in Clinical Natural Language Processing

abstract: Traditional representations like Bag of words are high dimensional, sparse and ignore the order as well as syntactic and semantic information. Distributed vector representations or embeddings map variable length text to dense fixed length vectors as well as capture the prior knowledge which can transferred to downstream tasks. Even though embedding has become de facto standard for representations in deep learning based NLP tasks in both general and clinical domains, there is no survey paper which presents a detailed review of embeddings in Clinical Natural Language Processing. In this survey paper, we discuss various medical corpora and their characteristics, medical codes and present a brief overview as well as comparison of popular embeddings models. We classify clinical embeddings into nine types and discuss each embedding type in detail. We discuss various evaluation methods followed by possible solutions to various challenges in clinical embeddings. Finally, we conclude with some of the future directions which will advance the research in clinical embeddings.


title: Wall effects of eccentric spheres machine learning for convenient computation

abstract: In confined systems, such as the inside of a biological cell, the outer boundary or wall can affect the dynamics of internal particles. In many cases of interest both the internal particle and outer wall are approximately spherical. Therefore, quantifying the wall effects from an outer spherical boundary on the motion of an internal eccentric sphere is very useful. However, when the two spheres are not concentric, the problem becomes non-trivial. In this paper we improve existing analytical methods to evaluate these wall effects and then train a feed-forward artificial neural network within a broader model. The final model generally performed with $\sim 0.001\%$ error within the training domain and $\sim 0.05\%$ when the outer spherical wall was extrapolated to an infinite plane. Through this model, the wall effects of an outer spherical boundary on the arbitrary motion of an internal sphere for all experimentally achievable configurations can now be conveniently and efficiently determined.


title: Cell Density Quantification with TurboSPI: R2* Mapping with Compensation for Off-Resonance Fat Modulation

abstract: Tracking the migration of superparamagnetic iron oxide (SPIO) labeled immune cells in vivo is valuable for understanding the immunogenic response to cancer and therapies. Quantitative cell tracking using compressed sensing TurboSPI-based R2* mapping is a promising development to improve accuracy in longitudinal studies on immune recruitment. The phase-encoded TurboSPI sequence provides high fidelity relaxation data in the form of signal time-courses with high temporal resolution. However, early in vivo applications of this method revealed that simple mono-exponential R2* fitting performs poorly due to the contaminant fat signal in voxels surrounding regions of interest, such as flank tumors and lymph nodes adjacent to adipose tissue. This is especially problematic if there is poor infiltration to the tumor such that immune cells remain near the periphery. The presence of an off-resonance fat isochromat results in modulations in the signal time-course can be erroneously fit as R2* signal decay, thereby overestimating the density of SPIO labeled cells. Simply excluding any voxel with fat-typical modulations results in underestimates in voxels that have mixed content. We propose using a more comprehensive dual-decay (R2f* and R2w*) Dixon-based signal model that accounts for the potential presence of fat in a voxel to better estimate SPIO induced de-phasing. In silico single voxel simulations illustrate how the proposed signal model provides stable R2w* estimates that are invariant to fat content. The proposed dual-decay model outperforms previous methods when applied to in vitro samples of SPIO labeled cells and oil prepared with oil content >15%. Preliminary in vivo results show that, compared to previous methods, the dual-decay Dixon model improves the balance of R2* specificity versus sensitivity, which in turn will result in more reliable analysis in future cell tracking studies.


title: CodeNet: Training Large Scale Neural Networks in Presence of Soft-Errors

abstract: This work proposes the first strategy to make distributed training of neural networks resilient to computing errors, a problem that has remained unsolved despite being first posed in 1956 by von Neumann. He also speculated that the efficiency and reliability of the human brain is obtained by allowing for low power but error-prone components with redundancy for error-resilience. It is surprising that this problem remains open, even as massive artificial neural networks are being trained on increasingly low-cost and unreliable processing units. Our coding-theory-inspired strategy, "CodeNet," solves this problem by addressing three challenges in the science of reliable computing: (i) Providing the first strategy for error-resilient neural network training by encoding each layer separately; (ii) Keeping the overheads of coding (encoding/error-detection/decoding) low by obviating the need to re-encode the updated parameter matrices after each iteration from scratch. (iii) Providing a completely decentralized implementation with no central node (which is a single point of failure), allowing all primary computational steps to be error-prone. We theoretically demonstrate that CodeNet has higher error tolerance than replication, which we leverage to speed up computation time. Simultaneously, CodeNet requires lower redundancy than replication, and equal computational and communication costs in scaling sense. We first demonstrate the benefits of CodeNet in reducing expected computation time over replication when accounting for checkpointing. Our experiments show that CodeNet achieves the best accuracy-runtime tradeoff compared to both replication and uncoded strategies. CodeNet is a significant step towards biologically plausible neural network training, that could hold the key to orders of magnitude efficiency improvements.


title: Non-Gaussian noise spectroscopy with a superconducting qubit sensor

abstract: Accurate characterization of the noise influencing a quantum system of interest has far-reaching implications across quantum science, ranging from microscopic modeling of decoherence dynamics to noise-optimized quantum control. While the assumption that noise obeys Gaussian statistics is commonly employed, noise is generically non-Gaussian in nature. In particular, the Gaussian approximation breaks down whenever a qubit is strongly coupled to discrete noise sources or has a non-linear response to the environmental degrees of freedom. Thus, in order to both scrutinize the applicability of the Gaussian assumption and capture distinctive non-Gaussian signatures, a tool for characterizing non-Gaussian noise is essential. Here, we experimentally validate a quantum control protocol which, in addition to the spectrum, reconstructs the leading higher-order spectrum of engineered non-Gaussian dephasing noise using a superconducting qubit as a sensor. This first experimental demonstration of non-Gaussian noise spectroscopy represents a major step toward demonstrating a complete spectral estimation toolbox for quantum devices.


title: Beyond the linear analysis of stability in higher derivative gravity with the Bianchi-I metric

abstract: The study of stability of gravitational perturbations in higher derivative gravity has shown that at the linear level the massive unphysical ghost is not generated from vacuum if the initial seed of metric perturbation has frequency essentially below the Planck threshold. The mathematical knowledge indicated that the linear stability is supposed to hold even at the nonperturbative level, but in such a complicated case it is important to perform a verification of this statement. We compare the asymptotic stability solutions at the linear and full nonperturbative levels for the Bianchi-I metric with small anisotropies, which can be regarded as an extreme, zero frequency limit of a gravitational wave. As one should expect from the combination of previous analysis and general mathematical theorems, there is a good correspondence between linear stability and the nonperturbative asymptotic behavior.


title: Robust commuter movement inference from connected mobile devices

abstract: The preponderance of connected devices provides unprecedented opportunities for fine-grained monitoring of the public infrastructure. However while classical models expect high quality application-specific data streams, the promise of the Internet of Things (IoT) is that of an abundance of disparate and noisy datasets from connected devices. In this context, we consider the problem of estimation of the level of service of a city-wide public transport network. We first propose a robust unsupervised model for train movement inference from wifi traces, via the application of robust clustering methods to a one dimensional spatio-temporal setting. We then explore the extent to which the demand-supply gap can be estimated from connected devices. We propose a classification model of real-time commuter patterns, including both a batch training phase and an online learning component. We describe our deployment architecture and assess our system accuracy on a large-scale anonymized dataset comprising more than 10 billion records.


title: Efficient quantum error correction of dephasing induced by a common fluctuator

abstract: Quantum error correction is expected to be essential in large-scale quantum technologies. However, the substantial overhead of qubits it requires is thought to greatly limit its utility in smaller, near-term devices. Here we introduce a new family of special-purpose quantum error-correcting codes that offer an exponential reduction in overhead compared to the usual repetition code. They are tailored for a common and important source of decoherence in current experiments, whereby a register of qubits is subject to phase noise through coupling to a common fluctuator, such as a resonator or a spin defect. The smallest instance encodes one logical qubit into two physical qubits, and corrects decoherence to leading-order using a constant number of one- and two-qubit operations. More generally, while the repetition code on $n$ qubits corrects errors to order $t^{O(n)}$, with $t$ the time between recoveries, our codes correct to order $t^{O(2^n)}$. Moreover, they are robust to model imperfections in small- and intermediate-scale devices, where they already provide substantial gains in error suppression. As a result, these hardware-efficient codes open a potential avenue for useful quantum error correction in near-term, pre-fault tolerant devices.


title: An FPT Algorithm for Minimum Additive Spanner Problem

abstract: For a positive integer $t$ and a graph $G$, an additive $t$-spanner of $G$ is a spanning subgraph in which the distance between every pair of vertices is at most the original distance plus $t$. Minimum Additive $t$-Spanner Problem is to find an additive $t$-spanner with the minimum number of edges in a given graph, which is known to be NP-hard. Since we need to care about global properties of graphs when we deal with additive $t$-spanners, Minimum Additive $t$-Spanner Problem is hard to handle, and hence only few results are known for it. In this paper, we study Minimum Additive $t$-Spanner Problem from the viewpoint of parameterized complexity. We formulate a parameterized version of the problem in which the number of removed edges is regarded as a parameter, and give a fixed-parameter algorithm for it. We also extend our result to $(\alpha, \beta)$-spanners.


title: Early Detection of Influenza outbreaks in the United States

abstract: Public health surveillance systems often fail to detect emerging infectious diseases, particularly in resource limited settings. By integrating relevant clinical and internet-source data, we can close critical gaps in coverage and accelerate outbreak detection. Here, we present a multivariate algorithm that uses freely available online data to provide early warning of emerging influenza epidemics in the US. We evaluated 240 candidate predictors and found that the most predictive combination does \textit{not} include surveillance or electronic health records data, but instead consists of eight Google search and Wikipedia pageview time series reflecting changing levels of interest in influenza-related topics. In cross validation on 2010-2016 data, this algorithm sounds alarms an average of 16.4 weeks prior to influenza activity reaching the Center for Disease Control and Prevention (CDC) threshold for declaring the start of the season. In an out-of-sample test on data from the rapidly-emerging fall wave of the 2009 H1N1 pandemic, it recognized the threat five weeks in advance of this surveillance threshold. Simpler algorithms, including fixed week-of-the-year triggers, lag the optimized alarms by only a few weeks when detecting seasonal influenza, but fail to provide early warning in the 2009 pandemic scenario. This demonstrates a robust method for designing next generation outbreak detection algorithms. By combining scan statistics with machine learning, it identifies tractable combinations of data sources (from among thousands of candidates) that can provide early warning of emerging infectious disease threats worldwide.


title: Charge and spin-specific local integrals of motion in a disordered Hubbard model

abstract: While many-body localization has primarily been studied in systems with a single local degree of freedom, experimental studies of many-body localization in cold atom systems motivate exploration of the disordered Hubbard model. With two coupled local degrees of freedom it is natural to ask how localization in spin relates to disorder in charge and vice versa. Most prior work has addressed disorder in only one of these sectors and often has not used measures of localization which distinguish between charge and spin. Here we explore localization in the Hubbard model with a wide range of independent values of charge and spin disorder, using measures of localization based on charge and spin-specific integrals of motion. Our results demonstrate symmetry between the response of the spin to charge disorder and vice versa, and we find very weak disorder in one channel, so long as the disorder in the other channel is sufficiently strong, results in localization in both channels. The strength of disorder required in the less disordered channel declines as the system size increases. Further, the weaker the disorder in the less-disordered channel, the longer the time scale at which localization appears in the dynamics of this degree of freedom.


title: Regression models for compositional data: General log-contrast formulations, proximal optimization, and microbiome data applications

abstract: Compositional data sets are ubiquitous in science, including geology, ecology, and microbiology. In microbiome research, compositional data primarily arise from high-throughput sequence-based profiling experiments. These data comprise microbial compositions in their natural habitat and are often paired with covariate measurements that characterize physicochemical habitat properties or the physiology of the host. Inferring parsimonious statistical associations between microbial compositions and habitat- or host-specific covariate data is an important step in exploratory data analysis. A standard statistical model linking compositional covariates to continuous outcomes is the linear log-contrast model. This model describes the response as a linear combination of log-ratios of the original compositions and has been extended to the high-dimensional setting via regularization. In this contribution, we propose a general convex optimization model for linear log-contrast regression which includes many previous proposals as special cases. We introduce a proximal algorithm that solves the resulting constrained optimization problem exactly with rigorous convergence guarantees. We illustrate the versatility of our approach by investigating the performance of several model instances on soil and gut microbiome data analysis tasks.


title: Same But Different: Distance Correlations Between Topological Summaries

abstract: Persistent homology allows us to create topological summaries of complex data. In order to analyse these statistically, we need to choose a topological summary and a relevant metric space in which this topological summary exists. While different summaries may contain the same information (as they come from the same persistence module), they can lead to different statistical conclusions since they lie in different metric spaces. The best choice of metric will often be application-specific. In this paper we discuss distance correlation, which is a non-parametric tool for comparing data sets that can lie in completely different metric spaces. In particular we calculate the distance correlation between different choices of topological summaries. We compare some different topological summaries for a variety of random models of underlying data via the distance correlation between the samples. We also give examples of performing distance correlation between topological summaries and other scalar measures of interest, such as a paired random variable or a parameter of the random model used to generate the underlying data. This article is meant to be expository in style, and will include the definitions of standard statistical quantities in order to be accessible to non-statisticians.


title: Primitive prime divisors in the critical orbits of one-parameter families of rational polynomials

abstract: For a rational polynomial $f$ and rational numbers $c, u$, we put $f_c(x):=f(x)+c$, and consider the Zsigmondy set $\mathcal{Z}(f_c,u)$ associated to the sequence $\{f_c^n(u)-u\}_{n\geq 0}$, where $f_c^n$ is the $n$-st iteration of $f_c$. In this paper, we prove that if $u$ is a rational critical point of $f$, then there exists an $\mathbf M_f>0$ such that $\mathbf M_f\geq \max_{c\in \mathbb{Q}}\{\#\mathcal{Z}(f_c,u)\}$.


title: Low-rank matrix recovery via regularized nuclear norm minimization

abstract: In this paper, we theoretically investigate the low-rank matrix recovery problem in the context of the unconstrained regularized nuclear norm minimization (RNNM) framework. Our theoretical findings show that, the RNNM method is able to provide a robust recovery of any matrix $X$ (not necessary to be exactly low-rank) from its few noisy measurements $\textbf{b}=\mathcal{A}(X)+\textbf{n}$ with a bounded constraint $\|\textbf{n}\|_{2}\leq\epsilon$, provided that the $tk$-order restricted isometry constant (RIC) of $\mathcal{A}$ satisfies a certain constraint related to $t>0$. Specifically, the obtained recovery condition in the case of $t>4/3$ is found to be same with the sharp condition established previously by Cai and Zhang (2014) to guarantee the exact recovery of any rank-$k$ matrix via the constrained nuclear norm minimization method. More importantly, to the best of our knowledge, we are the first to establish the $tk$-order RIC based coefficient estimate of the robust null space property in the case of $0<t\leq1$.


title: The Application of Multi-block ADMM on Isotonic Regression Problems

abstract: The multi-block ADMM has received much attention from optimization researchers due to its excellent scalability. In this paper, the multi-block ADMM is applied to solve two large-scale problems related to isotonic regression. Numerical experiments show that the multi-block ADMM is convergent when the chosen parameter is small enough and the multi-block ADMM scales well compared with baselines.


title: A Parallel Decomposition Scheme for Solving Long-Horizon Optimal Control Problems

abstract: We present a temporal decomposition scheme for solving long-horizon optimal control problems. In the proposed scheme, the time domain is decomposed into a set of subdomains with partially overlapping regions. Subproblems associated with the subdomains are solved in parallel to obtain local primal-dual trajectories that are assembled to obtain the global trajectories. We provide a sufficient condition that guarantees convergence of the proposed scheme. This condition states that the effect of perturbations on the boundary conditions (i.e., initial state and terminal dual/adjoint variable) should decay asymptotically as one moves away from the boundaries. This condition also reveals that the scheme converges if the size of the overlap is sufficiently large and that the convergence rate improves with the size of the overlap. We prove that linear quadratic problems satisfy the asymptotic decay condition, and we discuss numerical strategies to determine if the condition holds in more general cases. We draw upon a non-convex optimal control problem to illustrate the performance of the proposed scheme.


title: The flare catalog and the flare activity in the Kepler mission

abstract: We present a flare catalog of the {\it Kepler} mission using the long-cadence data of Data Release 25. This catalog comprises 3420 flare stars and 162,262 flare events. A comparison shows that the flare catalogs of previous studies are seriously polluted by various false positive signals and artifacts. The incidence of flare stars rises with decreasing temperature, which accords with the theoretical analysis. The flare frequency distributions (FFDs) from F-type stars to M-type stars obey a power-law relation with $\alpha \sim 2$, indicating that they have the same mechanism on generating flares. The remarkable incidence and the deviation of FFDs on A-type flare stars imply that they generate flares in a different way. The activity--rotation relation is consistent with previous studies at low temperature band, whereas it becomes dispersive with increasing temperature. Combined with the Gyrochronology, we find that the mixing of stars of two different dynamos gives rise to the dispersion. We thereby propose a scenario on understanding the activity--rotation relation across the H-R diagram. Based on the scenario and the correspondence of dynamo with regard to activity and rotation, we suggest a new expression on the activity--rotation relation, in which the segmentation is on the basis of the dynamo rather than the rotation period. The rotation distribution of flare stars shows that about 70\% of flare stars rotate faster than 10 days and the rate approaches 95\% at 30 days. Based on the incidence and the rotation distribution of flare stars, we estimate that a superflare with energy $\sim 10^{34}$ erg occurs on the Sun at least once in 5500 years.


title: Ultra-Scalable Spectral Clustering and Ensemble Clustering

abstract: This paper focuses on scalability and robustness of spectral clustering for extremely large-scale datasets with limited resources. Two novel algorithms are proposed, namely, ultra-scalable spectral clustering (U-SPEC) and ultra-scalable ensemble clustering (U-SENC). In U-SPEC, a hybrid representative selection strategy and a fast approximation method for K-nearest representatives are proposed for the construction of a sparse affinity sub-matrix. By interpreting the sparse sub-matrix as a bipartite graph, the transfer cut is then utilized to efficiently partition the graph and obtain the clustering result. In U-SENC, multiple U-SPEC clusterers are further integrated into an ensemble clustering framework to enhance the robustness of U-SPEC while maintaining high efficiency. Based on the ensemble generation via multiple U-SEPC's, a new bipartite graph is constructed between objects and base clusters and then efficiently partitioned to achieve the consensus clustering result. It is noteworthy that both U-SPEC and U-SENC have nearly linear time and space complexity, and are capable of robustly and efficiently partitioning ten-million-level nonlinearly-separable datasets on a PC with 64GB memory. Experiments on various large-scale datasets have demonstrated the scalability and robustness of our algorithms. The MATLAB code and experimental data are available at https://www.researchgate.net/publication/330760669.


title: Unifying quantum heat transfer and superradiant signature in a nonequilibrium collective-qubit system: a polaron-transformed Redfield approach

abstract: We investigate full counting statistics of quantum heat transfer in a collective-qubit system, constructed by multi-qubits interacting with two thermal baths. The nonequilibrium polaron-transformed Redfield approach embedded with an auxiliary counting field is applied to obtain the steady state heat current and fluctuations, which enables us to study the impact of the qubit-bath interaction in a wide regime. The heat current, current noise and skewness are all found to clearly unify the limiting results in the weak and strong couplings, respectively. Moreover, the superradiant heat transfer is clarified as a system-size-dependent effect, and large number of qubits dramatically suppresses the nonequilibrium superradiant signature.


title: Limit Theorems for Network Dependent Random Variables

abstract: This paper is concerned with cross-sectional dependence arising because observations are interconnected through an observed network. Following Doukhan and Louhichi (1999), we measure the strength of dependence by covariances of nonlinearly transformed variables. We provide a law of large numbers and central limit theorem for network dependent variables. We also provide a method of calculating standard errors robust to general forms of network dependence. For that purpose, we rely on a network heteroskedasticity and autocorrelation consistent (HAC) variance estimator, and show its consistency. The results rely on conditions characterized by tradeoffs between the rate of decay of dependence across a network and network's denseness. Our approach can accommodate data generated by network formation models, random fields on graphs, conditional dependency graphs, and large functional-causal systems of equations.


title: Indigenous use of stellar scintillation to predict weather and seasonal change

abstract: Indigenous peoples across the world observe the motions and positions of stars to develop seasonal calendars. Additionally, changing properties of stars, such as their brightness and colour, are also used for predicting weather. Combining archival studies with ethnographic fieldwork in Australia's Torres Strait, we explore the various ways Indigenous peoples utilise stellar scintillation (twinkling) as an indicator for predicting weather and seasonal change, discussing the scientific underpinnings of this knowledge. By observing subtle changes in the ways the stars twinkle, Meriam people gauge changing trade winds, approaching wet weather, and temperature changes. We then explore how the Northern Dene of Arctic North America utilise stellar scintillation to forecast weather.


title: Learning low-precision neural networks without Straight-Through Estimator(STE)

abstract: The Straight-Through Estimator (STE) is widely used for back-propagating gradients through the quantization function, but the STE technique lacks a complete theoretical understanding. We propose an alternative methodology called alpha-blending (AB), which quantizes neural networks to low-precision using stochastic gradient descent (SGD). Our method (AB) avoids STE approximation by replacing the quantized weight in the loss function by an affine combination of the quantized weight w_q and the corresponding full-precision weight w with non-trainable scalar coefficient $\alpha$ and $1-\alpha$. During training, $\alpha$ is gradually increased from 0 to 1; the gradient updates to the weights are through the full-precision term, $(1-\alpha)w$, of the affine combination; the model is converted from full-precision to low-precision progressively. To evaluate the method, a 1-bit BinaryNet on CIFAR10 dataset and 8-bits, 4-bits MobileNet v1, ResNet_50 v1/2 on ImageNet dataset are trained using the alpha-blending approach, and the evaluation indicates that AB improves top-1 accuracy by 0.9%, 0.82% and 2.93% respectively compared to the results of STE based quantization.


title: A caloritronics-based Mott neuristor

abstract: Machine learning imitates the basic features of biological neural networks to efficiently perform tasks such as pattern recognition. This has been mostly achieved at a software level, and a strong effort is currently being made to mimic neurons and synapses with hardware components, an approach known as neuromorphic computing. CMOS-based circuits have been used for this purpose, but they are non-scalable, limiting the device density and motivating the search for neuromorphic materials. While recent advances in resistive switching have provided a path to emulate synapses at the 10 nm scale, a scalable neuron analogue is yet to be found. Here, we show how heat transfer can be utilized to mimic neuron functionalities in Mott nanodevices. We use the Joule heating created by current spikes to trigger the insulator-to-metal transition in a biased VO2 nanogap. We show that thermal dynamics allow the implementation of the basic neuron functionalities: activity, leaky integrate-and-fire, volatility and rate coding. By using local temperature as the internal variable, we avoid the need of external capacitors, which reduces neuristor size by several orders of magnitude. This approach could enable neuromorphic hardware to take full advantage of the rapid advances in memristive synapses, allowing for much denser and complex neural networks. More generally, we show that heat dissipation is not always an undesirable effect: it can perform computing tasks if properly engineered.


title: NoRML: No-Reward Meta Learning

abstract: Efficiently adapting to new environments and changes in dynamics is critical for agents to successfully operate in the real world. Reinforcement learning (RL) based approaches typically rely on external reward feedback for adaptation. However, in many scenarios this reward signal might not be readily available for the target task, or the difference between the environments can be implicit and only observable from the dynamics. To this end, we introduce a method that allows for self-adaptation of learned policies: No-Reward Meta Learning (NoRML). NoRML extends Model Agnostic Meta Learning (MAML) for RL and uses observable dynamics of the environment instead of an explicit reward function in MAML's finetune step. Our method has a more expressive update step than MAML, while maintaining MAML's gradient based foundation. Additionally, in order to allow more targeted exploration, we implement an extension to MAML that effectively disconnects the meta-policy parameters from the fine-tuned policies' parameters. We first study our method on a number of synthetic control problems and then validate our method on common benchmark environments, showing that NoRML outperforms MAML when the dynamics change between tasks.


title: The Duality in Quantum Work

abstract: An open question of fundamental importance in quantum thermodynamics is how to describe the statistics of work for initial state with quantum coherence. In this paper, work statistics is considered from a fully new perspective of "wave-particle" duality. Based on the generalized quantum work measurement, predictability of energy levels $\mathcal{D}_W$ and effectiveness of coherence $\mathcal{V}_W$ are defined, and they obey inequality $\mathcal{D}_W^2+\mathcal{V}_W^2\leq1$, which is the fundamental tradeoff relations between the contributions of population and coherence to quantum work distribution. As an application, we consider a driven two-level system and discuss the condition of the bound of above tradeoff relation. These results shed light on the effects of quantum coherence in quantum thermodynamics.


title: Distributed Optimal Voltage Control with Asynchronous and Delayed Communication

abstract: The increased penetration of volatile renewable energy into distribution networks necessities more efficient distributed voltage control. In this paper, we design distributed feedback control algorithms where each bus can inject \emph{both active and reactive} power into the grid to regulate the voltages. The control law on each bus is only based on local voltage measurements and communication to its physical neighbors. Moreover, the buses can perform their updates \emph{asynchronously} without receiving information from their neighbors for periods of time. The algorithm enforces \emph{hard upper and lower limits} on the active and reactive powers at every iteration. We prove that the algorithm converges to the optimal feasible voltage profile, assuming linear power flows. This provable convergence is maintained under bounded communication delays and asynchronous communications. We further numerically test the performance of the algorithm using the full \emph{nonlinear AC power flow} model. Our simulations show the effectiveness of our algorithm on realistic networks with both static and fluctuating loads, even in the presence of communication delays.


title: Reinforcement Learning on Variable Impedance Controller for High-Precision Robotic Assembly

abstract: Precise robotic manipulation skills are desirable in many industrial settings, reinforcement learning (RL) methods hold the promise of acquiring these skills autonomously. In this paper, we explicitly consider incorporating operational space force/torque information into reinforcement learning; this is motivated by humans heuristically mapping perceived forces to control actions, which results in completing high-precision tasks in a fairly easy manner. Our approach combines RL with force/torque information by incorporating a proper operational space force controller; where we also exploit different ablations on processing this information. Moreover, we propose a neural network architecture that generalizes to reasonable variations of the environment. We evaluate our method on the open-source Siemens Robot Learning Challenge, which requires precise and delicate force-controlled behavior to assemble a tight-fit gear wheel set.


title: Incremental Visual-Inertial 3D Mesh Generation with Structural Regularities

abstract: Visual-Inertial Odometry (VIO) algorithms typically rely on a point cloud representation of the scene that does not model the topology of the environment. A 3D mesh instead offers a richer, yet lightweight, model. Nevertheless, building a 3D mesh out of the sparse and noisy 3D landmarks triangulated by a VIO algorithm often results in a mesh that does not fit the real scene. In order to regularize the mesh, previous approaches decouple state estimation from the 3D mesh regularization step, and either limit the 3D mesh to the current frame or let the mesh grow indefinitely. We propose instead to tightly couple mesh regularization and state estimation by detecting and enforcing structural regularities in a novel factor-graph formulation. We also propose to incrementally build the mesh by restricting its extent to the time-horizon of the VIO optimization; the resulting 3D mesh covers a larger portion of the scene than a per-frame approach while its memory usage and computational complexity remain bounded. We show that our approach successfully regularizes the mesh, while improving localization accuracy, when structural regularities are present, and remains operational in scenes without regularities.


title: Radon numbers and the fractional Helly theorem

abstract: A basic measure of the combinatorial complexity of a convexity space is its Radon number. In this paper we show a fractional Helly theorem for convexity spaces with a bounded Radon number, answering a question of Kalai. As a consequence we also get a weak epsilon-net theorem for convexity spaces with a bounded Radon number. This answers a question of Bukh and extends a recent result of Moran and Yehudayoff.


title: Neural Networks Trained on Natural Scenes Exhibit Gestalt Closure

abstract: The Gestalt laws of perceptual organization, which describe how visual elements in an image are grouped and interpreted, have traditionally been thought of as innate despite their ecological validity. We use deep-learning methods to investigate whether natural scene statistics might be sufficient to derive the Gestalt laws. We examine the law of closure, which asserts that human visual perception tends to "close the gap" by assembling elements that can jointly be interpreted as a complete figure or object. We demonstrate that a state-of-the-art convolutional neural network, trained to classify natural images, exhibits closure on synthetic displays of edge fragments, as assessed by similarity of internal representations. This finding provides support for the hypothesis that the human perceptual system is even more elegant than the Gestaltists imagined: a single law---adaptation to the statistical structure of the environment---might suffice as fundamental.


title: On the mass accretion rate and infrared excess in Herbig Ae/Be Stars

abstract: The present study makes use of the unprecedented capability of the Gaia mission to obtain the stellar parameters such as distance, age, and mass of HAeBe stars. The accuracy of Gaia DR2 astrometry is demonstrated from the comparison of the Gaia DR2 distances of 131 HAeBe stars with the previously estimated values from the literature. This is one of the initial studies to estimate the age and mass of a confirmed sample of HAeBe stars using both the photometry and distance from the Gaia mission. Mass accretion rates are calculated from $H\alpha$ line flux measurements of 106 HAeBe stars. Since we used distances and the stellar masses derived from the Gaia DR2 data in the calculation of mass accretion rate, our estimates are more accurate than previous studies. The mass accretion rate is found to decay exponentially with age, from which we estimated a disk dissipation timescale of $1.9\pm 0.1$ Myr. Mass accretion rate and stellar mass exhibits a power law relation of the form, $\dot{M}_{acc}$ $\propto$ $M_{*}^{2.8\pm0.2}$. From the distinct distribution in the values of the infrared spectral index, $n_{2-4.6}$, we suggest the possibility of difference in the disk structure between Herbig Be and Herbig Ae stars.


title: Real-Time Source Independent Quantum Random Number Generator with Squeezed States

abstract: Random numbers are a fundamental ingredient for many applications including simulation, modelling and cryptography. Sound random numbers should be independent and uniformly distributed. Moreover, for cryptographic applications they should also be unpredictable. We demonstrate a real-time self-testing source independent quantum random number generator (QRNG) that uses squeezed light as source. We generate secure random numbers by measuring the quadratures of the electromagnetic field without making any assumptions on the source; only the detection device is trusted. We use a homodyne detection to alternatively measure the Q and P conjugate quadratures of our source. Using the entropic uncertainty relation, measurements on P allow us to estimate a bound on the min-entropy of Q conditioned on any classical or quantum side information that a malicious eavesdropper may detain. This bound gives the minimum number of secure bits we can extract from the Q measurement. We discuss the performance of different estimators for this bound. We operate this QRNG with a squeezed state and we compare its performance with a QRNG using thermal states. The real-time bit rate was 8.2 kb/s when using the squeezed source and between 5.2-7.2 kb/s when the thermal state source was used.


title: COMIC: Towards A Compact Image Captioning Model with Attention

abstract: Recent works in image captioning have shown very promising raw performance. However, we realize that most of these encoder-decoder style networks with attention do not scale naturally to large vocabulary size, making them difficult to be deployed on embedded system with limited hardware resources. This is because the size of word and output embedding matrices grow proportionally with the size of vocabulary, adversely affecting the compactness of these networks. To address this limitation, this paper introduces a brand new idea in the domain of image captioning. That is, we tackle the problem of compactness of image captioning models which is hitherto unexplored. We showed that, our proposed model, named COMIC for COMpact Image Captioning, achieves comparable results in five common evaluation metrics with state-of-the-art approaches on both MS-COCO and InstaPIC-1.1M datasets despite having an embedding vocabulary size that is 39x - 99x smaller. The source code and models are available at: https://github.com/jiahuei/COMIC-Compact-Image-Captioning-with-Attention


title: Designing Optimal Multiplex Networks for Certain Laplacian Spectral Properties

abstract: We discuss the design of interlayer edges in a multiplex network, under a limited budget, with the goal of improving its overall performance. We analyze the following three problems separately; first, we maximize the smallest nonzero eigenvalue, also known as the algebraic connectivity; secondly, we minimize the largest eigenvalue, also known as the spectral radius; and finally, we minimize the spectral width. Maximizing the algebraic connectivity requires identical weights on the interlayer edges for budgets less than a threshold value. However, for larger budgets, the optimal weights are generally non-uniform. The dual formulation transforms the problem into a graph realization (embedding) problem that allows us to give a fuller picture. Namely, before the threshold budget, the optimal realization is one-dimensional with nodes in the same layer embedded to a single point; while, beyond the threshold, the optimal embeddings generally unfold into spaces with dimension bounded by the multiplicity of the algebraic connectivity. Finally, for extremely large budgets the embeddings revert again to lower dimensions. Minimizing the largest eigenvalue is driven by the spectral radius of the individual networks and its corresponding eigenvector. Before a threshold, the total budget is distributed among interlayer edges corresponding to the nodal lines of this eigenvector, and the optimal largest eigenvalue of the Laplacian remains constant. For larger budgets, the weight distribution tends to be almost uniform. In the dual picture, the optimal graph embedding is one-dimensional and non-homogeneous at first and beyond this threshold, the optimal embedding expands to be multi-dimensional, and for larger values of the budget, the two layers fill the embedding space. Finally, we show how these two problems are connected to minimizing the spectral width.


title: Reflection of a Point Object in an Arbitrary Curved Mirror

abstract: In this work, I have derived the equation of the curve obtained on reflection of a point object in an arbitrary curved mirror if the object and the mirror are placed on the 2D Cartesian plane. I have used only the basic laws of reflection of classical geometric optics and elementary coordinate geometry. Several examples are provided and compared with Gaussian optics. We also see how the equations reduce to the standard mirror formula under the paraxial approximation.


title: Isoscalar dipole excitations in $^{16}$O

abstract: Isoscalar (IS) monopole and dipole excitations in $^{16}$O were investigated by the method of shifted basis antisymmetrized molecular dynamics combined with the generator coordinate method. Significant strengths of the IS monopole and dipole transitions were obtained in the low-energy region below the giant resonances. In addition to the compressive mode, which mainly contributes to the high-energy strengths for the IS dipole giant resonance, we obtained a variety of low-energy dipole modes such as the vortical dipole mode in the $1^-_1$ state of the vibrating tetrahedral $4\alpha$ and the $^{12}$C+$\alpha$ cluster structure in the $1^-_2$ state. The $1^-_1$ state contributes to the significant low-energy strength of the IS dipole transition as 5\% of the energy-weighted sum rule, which describes well the experimental data observed by the $\alpha$ inelastic scattering.


title: On the finiteness of solutions for polynomial-factorial Diophantine equations

abstract: We study the Diophantine equations obtained by equating a polynomial and the factorial function, and prove the finiteness of integer solutions under certain conditions. For example, we show that there exists only finitely many $l$ such that $l!$ is represented {by} $N_A(x)$, where $N_A$ is a norm form constructed from the field norm of a field extension $K/\mathbf Q$. We also deal with the equation $N_A(x)=l!_S$, where $l!_S$ is the Bhargava factorial. In this paper, we also show that the Oesterl\'e-Masser conjecture implies that for any infinite subset $S$ of $\mathbf Z$ and for any polynomial $P(x)\in\mathbf Z[x]$ of degree $2$ or more the equation $P(x)=l!_S$ has only finitely many solutions $(x,l)$. For some special infinite subsets $S$ of $\mathbf Z$, we can show the finiteness of solutions for the equation $P(x)=l!_S$ unconditionally.


title: Flux density measurements for 32 pulsars in the 20 cm band

abstract: Flux density measurements provide fundamental observational parameters that describe a pulsar. In the current pulsar catalogue, 27% of radio pulsars have no flux density measurement in the 20 cm observing band. Here, we present the first measurements of the flux densities in this band for 32 pulsars observed using the Parkes radio telescope and provide updated pulse profiles for these pulsars. We have used both archival and new observations to make these measurements. Various schemes exist for measuring flux densities. We show how the flux densities measured vary between these methods and how the presence of radio-frequency-interference will bias flux density measurements


title: Unsupervised Cross-spectral Stereo Matching by Learning to Synthesize

abstract: Unsupervised cross-spectral stereo matching aims at recovering disparity given cross-spectral image pairs without any supervision in the form of ground truth disparity or depth. The estimated depth provides additional information complementary to individual semantic features, which can be helpful for other vision tasks such as tracking, recognition and detection. However, there are large appearance variations between images from different spectral bands, which is a challenge for cross-spectral stereo matching. Existing deep unsupervised stereo matching methods are sensitive to the appearance variations and do not perform well on cross-spectral data. We propose a novel unsupervised cross-spectral stereo matching framework based on image-to-image translation. First, a style adaptation network transforms images across different spectral bands by cycle consistency and adversarial learning, during which appearance variations are minimized. Then, a stereo matching network is trained with image pairs from the same spectra using view reconstruction loss. At last, the estimated disparity is utilized to supervise the spectral-translation network in an end-to-end way. Moreover, a novel style adaptation network F-cycleGAN is proposed to improve the robustness of spectral translation. Our method can tackle appearance variations and enhance the robustness of unsupervised cross-spectral stereo matching. Experimental results show that our method achieves good performance without using depth supervision or explicit semantic information.


title: Topological conjugacy between induced non-autonomous set-valued systems and subshifts of finite type

abstract: This paper establishes topological (equi-)semiconjugacy and (equi-)conjugacy between induced non-autonomous set-valued systems and subshifts of finite type. First, some necessary and sufficient conditions are given for a non-autonomous discrete system to be topologically semiconjugate or conjugate to a subshift of finite type. Further, several sufficient conditions for it to be topologically equi-semiconjugate or equi-conjugate to a subshift of finite type are obtained. Consequently, estimations of topological entropy and several criteria of Li-Yorke chaos and distributional chaos in a sequence are derived. Second, the relationships of several related dynamical behaviors between the non-autonomous discrete system and its induced set-valued system are investigated. Based on these results, the paper furthermore establishes the topological (equi-)semiconjugacy and (equi-)conjugacy between induced set-valued systems and subshifts of finite type. Consequently, estimations of the topological entropy for the induced set-valued system are obtained, and several criteria of Li-Yorke chaos and distributional chaos in a sequence are established. Some of these results not only extend the existing related results for autonomous discrete systems to non-autonomous discrete systems, but also relax the assumptions of the counterparts in the literature. Two examples are finally provided for illustration.


title: From Knowledge Map to Mind Map: Artificial Imagination

abstract: Imagination is one of the most important factors which makes an artistic painting unique and impressive. With the rapid development of Artificial Intelligence, more and more researchers try to create painting with AI technology automatically. However, lacking of imagination is still a main problem for AI painting. In this paper, we propose a novel approach to inject rich imagination into a special painting art Mind Map creation. We firstly consider lexical and phonological similarities of seed word, then learn and inherit original painting style of the author, and finally apply Dadaism and impossibility of improvisation principles into painting process. We also design several metrics for imagination evaluation. Experimental results show that our proposed method can increase imagination of painting and also improve its overall quality.


title: CloudPSS: A High-Performance Power System Simulator Based on Cloud Computing

abstract: With the increasing computations in power system simulations, high-performance and cost-effective power system simulator is highly required. In this paper, a cloud-computing based power system simulator, namely CloudPSS, is designed. Based on an open service integrating framework, a self-developed electromagnetic transients (EMT) simulator with an automatic code generator is provided to accelerate EMT simulations using heterogeneous devices in the cloud, such as CPU and GPU. Test results show that CloudPSS can achieve significant speedups for both large-scale simulation tasks and multi-scenario tasks. Moreover, benefit from the time-sharing of cloud computing resources, the computational cost can be greatly reduced.


title: Exact Solution for the Portfolio Diversification Problem Based on Maximizing the Risk Adjusted Return

abstract: The potential benefits of portfolio diversification have been known to investors for a long time. Markowitz (1952) suggested the seminal approach for optimizing the portfolio problem based on finding the weights as budget shares that minimize the variance of the underlying portfolio. Hatemi-J and El-Khatib (2015) suggested finding the weights that will result in maximizing the risk adjusted return of the portfolio. This approach seems to be preferred by the rational investors since it combines risk and return when the optimal budget shares are sought for. The current paper provides a general solution for this risk adjusted return problem that can be utilized for any potential number of assets that are included in the portfolio.


title: Stochastic Online Learning with Probabilistic Graph Feedback

abstract: We consider a problem of stochastic online learning with general probabilistic graph feedback, where each directed edge in the feedback graph has probability $p_{ij}$. Two cases are covered. (a) The one-step case, where after playing arm $i$ the learner observes a sample reward feedback of arm $j$ with independent probability $p_{ij}$. (b) The cascade case where after playing arm $i$ the learner observes feedback of all arms $j$ in a probabilistic cascade starting from $i$ -- for each $(i,j)$ with probability $p_{ij}$, if arm $i$ is played or observed, then a reward sample of arm $j$ would be observed with independent probability $p_{ij}$. Previous works mainly focus on deterministic graphs which corresponds to one-step case with $p_{ij} \in \{0,1\}$, an adversarial sequence of graphs with certain topology guarantees, or a specific type of random graphs. We analyze the asymptotic lower bounds and design algorithms in both cases. The regret upper bounds of the algorithms match the lower bounds with high probability.


title: Automatic microscopic cell counting by use of deeply-supervised density regression model

abstract: Accurately counting cells in microscopic images is important for medical diagnoses and biological studies, but manual cell counting is very tedious, time-consuming, and prone to subjective errors, and automatic counting can be less accurate than desired. To improve the accuracy of automatic cell counting, we propose here a novel method that employs deeply-supervised density regression. A fully convolutional neural network (FCNN) serves as the primary FCNN for density map regression. Innovatively, a set of auxiliary FCNNs are employed to provide additional supervision for learning the intermediate layers of the primary CNN to improve network performance. In addition, the primary CNN is designed as a concatenating framework to integrate multi-scale features through shortcut connections in the network, which improves the granularity of the features extracted from the intermediate CNN layers and further supports the final density map estimation.


title: Orbital Torque: Torque Generation by Orbital Current Injection

abstract: We propose a mechanism of torque generation by injection of an orbital current, which we call $\textit{orbital torque}$. In a magnetic bilayer consisting of a nonmagnet (NM) and a ferromagnet (FM), we consider a situation where the spin-orbit coupling (SOC) is present only in the FM. Although the SOC is absent in the NM, the orbital Hall effect can arise in the NM. When the resulting orbital Hall current is injected to the FM, the SOC of the FM converts the orbital angular momentum into spin, which exerts torque to the magnetization of the FM. Remarkably, even for small SOC strength comparable to that of $3d$ FMs, the orbital torque can be comparable to the spin torque induced by the spin Hall effect of the NM with strong SOC. This provides a way to experimentally probe the OHE and opens a new venue to achieving spin-torque devices based on light elements that exhibit gigantic orbital response. Experimental implications are discussed.


title: Role of Faddeevian anomaly in the s-wave scattering of chiral fermion off dilaton black holes towards preservation of information

abstract: It was found that s-wave scattering of chiral fermion off dilaton black-hole when studied with a model generated from the chiral Schwinger model with standard Jackiw-Rajaraman type of anomaly provided information non preserving result. However this scattering problem when studied with a model generated from chiral Schwinger model with generalized Faddeevian type of anomaly. rendered information preserving result and it had well agreement with Hawking's revised proposal related to information loss. A minute and equitable investigation in detail has been carried out here to show how Faddeevian type of anomaly scores over the Jackiw-Rajaraman type of anomaly in connection with the s-wave scattering of chiral fermion of dilaton black-hole related to preservation of information.


title: Borcherds' method for Enriques surfaces

abstract: We classify all primitive embeddings of the lattice of numerical equivalence classes of divisors of an Enriques surface with the intersection form multiplied by 2 into an even unimodular hyperbolic lattice of rank 26. These embeddings have a property that facilitates the computation of the automorphism group of an Enriques surface by Borcherds' method.


title: Wasserstein Hamiltonian flows

abstract: We establish kinetic Hamiltonian flows in density space embedded with the $L^2$-Wasserstein metric tensor. We derive the Euler-Lagrange equation in density space, which introduces the associated Hamiltonian flows. We demonstrate that many classical equations, such as Vlasov equation, Schr{\"o}dinger equation and Schr{\"o}dinger bridge problem, can be rewritten as the formalism of Hamiltonian flows in density space.


title: On the volume of sectional-hyperbolic sets

abstract: For a transitive sectional-hypebolic set $\Lambda$ with positive volume on a $d$-dimensional manifold $M$($d\ge3$), we show that $\Lambda=M$ and $\Lambda$ is a uniformly hyperbolic set without singularities


title: Development of an alpha-particle imaging detector based on a low radioactivity micro-time-projection chamber

abstract: An important issue for rare-event-search experiments, such as the search for dark matter or neutrinoless double beta decay, is to reduce radioactivity of the detector materials and the experimental environment. The selection of materials with low radioactive impurities, such as isotopes of the uranium and thorium chains, requires a precise measurement of surface and bulk radioactivity. Focused on the first one, an alpha-particle detector has been developed based on a gaseous micro-time-projection chamber. A low-alpha mu-PIC with reduced alpha-emission background was installed in the detector. The detector offers the advantage of position sensitivity, which allows the alpha-particle contamination of the sample to be imaged and the background to be measured at the same time. The detector performance was measured by using an alpha-particle source. The measurement with a sample was also demonstrated and the sensitivity is discussed.


title: A systematic study of soft X-ray pulse profiles of magnetars in quiescence

abstract: Magnetars are neutron stars with extremely high surface magnetic fields. They show diverse X-ray pulse profiles in the quiescent state. We perform a systematic Fourier analysis of their soft X-ray pulse profiles. We find that most magnetars have a single-peaked profile and hence have low amplitudes of the second Fourier harmonics (A2). On the other hand, the pulsed fraction (PF) spreads over a wide range. We compared the results with theoretical profiles assuming various surface hotspot asymmetries, viewing geometries, and beaming functions. We found that a single value of the intensity ratio r between two antipodal hotspots is unable to reproduce the observed distribution of A2 and PF for all magnetars. The inferred r is probably anticorrelated with the thermal luminosity, implying that high-luminosity magnetars tend to have two symmetric hotspots. Our results are consistent with theoretical predictions, for which the existence of an evolving toroidal magnetic field breaks the symmetry of the surface temperature.


title: Zero-Shot Task Transfer

abstract: In this work, we present a novel meta-learning algorithm, i.e. TTNet, that regresses model parameters for novel tasks for which no ground truth is available (zero-shot tasks). In order to adapt to novel zero-shot tasks, our meta-learner learns from the model parameters of known tasks (with ground truth) and the correlation of known tasks to zero-shot tasks. Such intuition finds its foothold in cognitive science, where a subject (human baby) can adapt to a novel-concept (depth understanding) by correlating it with old concepts (hand movement or self-motion), without receiving explicit supervision. We evaluated our model on the Taskonomy dataset, with four tasks as zero-shot: surface-normal, room layout, depth, and camera pose estimation. These tasks were chosen based on the data acquisition complexity and the complexity associated with the learning process using a deep network. Our proposed methodology out-performs state-of-the-art models (which use ground truth)on each of our zero-shot tasks, showing promise on zero-shot task transfer. We also conducted extensive experiments to study the various choices of our methodology, as well as showed how the proposed method can also be used in transfer learning. To the best of our knowledge, this is the firstsuch effort on zero-shot learning in the task space.


title: Differentiable Causal Computations via Delayed Trace

abstract: We investigate causal computations taking sequences of inputs to sequences of outputs where the $n$th output depends on the first $n$ inputs only. We model these in category theory via a construction taking a Cartesian category $C$ to another category $St(C)$ with a novel trace-like operation called "delayed trace", which misses yanking and dinaturality axioms of the usual trace. The delayed trace operation provides a feedback mechanism in $St(C)$ with an implicit guardedness guarantee. When $C$ is equipped with a Cartesian differential operator, we construct a differential operator for $St(C)$ using an abstract version of backpropagation through time, a technique from machine learning based on unrolling of functions. This obtains a swath of properties for backpropagation through time, including a chain rule and Schwartz theorem. Our differential operator is also able to compute the derivative of a stateful network without requiring the network to be unrolled.


title: The generalized Auslander-Reiten duality on a module category

abstract: We characterize the generalized Auslander--Reiten duality on the category of finitely presented modules over some certain Hom-finite category. Examples include the category FI of finite sets with injections, and the one VI of finite dimensional vector spaces with linear injections over a finite field.


title: The Number of Convex Polyominoes with Given Height and Width

abstract: We give a new combinatorial proof for the number of convex polyominoes whose minimum enclosing rectangle has given dimensions. We also count the subclass of these polyominoes that contain the lower left corner of the enclosing rectangle (directed polyominoes). We indicate how to sample random polyominoes in these classes. As a side result, we calculate the first and second moments of the number of common points of two monotone lattice paths between two given points.


title: Background modelling for $\gamma$-ray spectroscopy with INTEGRAL/SPI

abstract: The coded-mask spectrometer-telescope SPI on board the INTEGRAL observatory records photons in the energy range between 20 and 8000 keV. A robust and versatile method to model the dominating instrumental background (BG) radiation is difficult to establish for such a telescope in the rapidly changing space environment. From long-term monitoring of SPI's Germanium detectors, we built up a spectral parameter data base, which characterises the instrument response as well as the BG behaviour. We aim to build a self-consistent and broadly applicable BG model for typical science cases of INTEGRAL/SPI, based on this data base. The general analysis method for SPI relies on distinguishing between illumination patterns on the 19-element Germanium detector array from BG and sky in a maximum likelihood framework. We illustrate how the complete set of measurements, even including the exposures of the sources of interest, can be used to define a BG model. We apply our method to different science cases, including point-like, diffuse, continuum, and line emission, and evaluate the adequacy in each case. From likelihood values and the number of fitted parameters, we determine how strong the impact of the unknown BG variability is. We find that the number of fitted parameters, i.e. how often the BG has to be re-normalised, depends on the emission type (diffuse with many observations over a large sky region, or point-like with concentrated exposure around one source), and the spectral energy range and bandwidth. A unique time scale, valid for all analysis issues, is not applicable for INTEGRAL/SPI, but must and can be inferred from the chosen data set. We conclude that our BG modelling method is usable in a large variety of INTEGRAL/SPI science cases, and provides nearly systematics-free and robust results.


title: Quasi-Continuous Metasurfaces for Orbital Angular Momentum Generation

abstract: A quasi-continuous composite perfect electric conductor-perfect magnetic conductor metasurface and a systematic metasurface design process are proposed for the orbital angular momentum (OAM) generation. The metasurfaces reflect the incident left circularly polarized (LCP)/right circularly polarized (RCP) plane wave to RCP/LCP vortex beams carrying OAM at normal or oblique direction. Unlike conventional metasurfaces that are composed of discrete scatterers, the scatterers on the proposed metasurface form a quasi-continuous pattern. The patterning of the metasurface is calculated through grating vectors, and no optimization of single scatterer is required. Furthermore, the distortions from local-response discontinuity of discrete scatterers are avoided. This letter provides great convenience to high-quality OAM generation.


title: Quadratic residues and related permutations

abstract: Let $p$ be an odd prime. For any $p$-adic integer $a$ we let $\overline{a}$ denote the unique integer $x$ with $-p/2<x<p/2$ and $x-a$ divisible by $p$. In this paper we study some permutations involving quadratic residues modulo $p$. For instance, we consider the following three sequences. \begin{align*} &A_0: \overline{1^2},\ \overline{2^2},\ \cdots,\ \overline{((p-1)/2)^2},\\ &A_1: \overline{a_1},\ \overline{a_2},\ \cdots,\ \overline{a_{(p-1)/2}},\\ &A_2: \overline{g^2},\ \overline{g^4},\ \cdots,\ \overline{g^{p-1}}, \end{align*} where $g\in\Z$ is a primitive root modulo $p$ and $1\le a_1<a_2<\cdots<a_{(p-1)/2}\le p-1$ are all quadratic residues modulo $p$. Obviously $A_i$ is a permutation of $A_j$ and we call this permutation $\sigma_{i,j}$. Sun obtained the sign of $\sigma_{0,1}$ when $p\equiv 3\pmod4$. In this paper we give the sign of $\sigma_{0,1}$ and determine the sign $\sigma_{0,2}$ when $p\equiv 1\pmod 4$.


title: A graphical categorification of the two-variable Chebyshev polynomials of the second kind

abstract: We show that the $A_2$ clasps in the Karoubi envelope of $A_2$ spider satisfy the recursive formula of the two-variable Chebyshev polynomials of the second kind associated with a root system of type $A_2$. The $A_2$ spider is a diagrammatic description of the representation category for $U_q(\mathfrak{sl}_3)$ and the $A_2$ clasps are projectors. Our categorification also gives a natural definition of a $q$-deformation of the two-variable Chebyshev polynomials. This paper is constructed based only on the linear skein theory and graphical calculus.


title: Trace operator on von Koch's snowflake

abstract: We study properties of the boundary trace operator on the Sobolev space $W^1_1(\Omega)$. Using the density result by Koskela and Zhang, we define a surjective operator \mbox{$Tr: W^1_1(\Omega_K)\rightarrow X(\Omega_K)$}, where $\Omega_K$ is von Koch's snowflake and $X(\Omega_K)$ is a trace space with the quotient norm. Since $\Omega_K$ is a uniform domain whose boundary is Ahlfors-regular with an exponent strictly bigger than one, it was shown by L. Mal\'y that there exists a right inverse to $Tr$, i.e. a linear operator $S: X(\Omega_K) \rightarrow W^1_1(\Omega_K)$ such that $Tr \circ S= Id_{X(\Omega_K)}$. In this paper we provide a different, purely combinatorial proof based on geometrical structure of von Koch's snowflake. Moreover we identify the isomorphism class of the trace space as $\ell_1$. As an additional consequence of our approach we obtain a simple proof of the Peetre's theorem about non-existence of the right inverse for domain $\Omega$ with regular boundary, which explains Banach space geometry cause for this phenomenon.


title: A difference-of-convex approach for split feasibility with applications to matrix factorizations and outlier detection

abstract: The split feasibility problem is to find an element in the intersection of a closed set $C$ and the linear preimage of another closed set $D$, assuming the projections onto $C$ and $D$ are easy to compute. This class of problems arises naturally in many contemporary applications such as compressed sensing. While the sets $C$ and $D$ are typically assumed to be convex in the literature, in this paper, we allow both sets to be possibly nonconvex. We observe that, in this setting, the split feasibility problem can be formulated as an optimization problem with a difference-of-convex objective so that standard majorization-minimization type algorithms can be applied. Here we focus on the nonmonotone proximal gradient algorithm with majorization studied in [15, Appendix A]. We show that, when this algorithm is applied to a split feasibility problem, the sequence generated clusters at a stationary point of the problem under mild assumptions. We also study local convergence property of the sequence under suitable assumptions on the closed sets involved. Finally, we perform numerical experiments to illustrate the efficiency of our approach on solving split feasibility problems that arise in completely positive matrix factorization, (uniformly) sparse matrix factorization, and outlier detection.