title: Spiral arms in CALIFA galaxies traced by non-circular velocities, abundances, and extinctions

abstract: We derive maps of the observed velocity of ionized gas, the oxygen abundance, and the extinction (Balmer decrement) across the area of the four spiral galaxies NGC36, NGC180, NGC6063, and NGC7653 from integral field spectroscopy obtained by the Calar Alto Legacy Integral Field Area (CALIFA) survey. We searched for spiral arms through Fourier analysis of the spatial distribution of three tracers (non-circular motion, enhancement of the oxygen abundance, and of the extinction) in the discs of our target galaxies. The spiral arms (two-armed logarithmic spirals in the deprojected map) are shown in each target galaxy for each tracer considered. The pitch angles of the spiral arms in a given galaxy obtained with the three different tracers are close to each other. The enhancement of the oxygen abundance in the spiral arms as compared to the abundance in the interarm regions at a given galactocentric distance is small; within a few per cent. We identified a metallicity gradient in our target galaxies. Both barred galaxies in our sample show flatter gradients than the two galaxies without bars. Galactic inclination, position angle of the major axis, and the rotation curve were also obtained for each target galaxy using the Fourier analysis of the two-dimensional velocity map.


title: Learning Structural Node Embeddings Via Diffusion Wavelets

abstract: Nodes residing in different parts of a graph can have similar structural roles within their local network topology. The identification of such roles provides key insight into the organization of networks and can be used for a variety of machine learning tasks. However, learning structural representations of nodes is a challenging problem, and it has typically involved manually specifying and tailoring topological features for each node. In this paper, we develop GraphWave, a method that represents each node's network neighborhood via a low-dimensional embedding by leveraging heat wavelet diffusion patterns. Instead of training on hand-selected features, GraphWave learns these embeddings in an unsupervised way. We mathematically prove that nodes with similar network neighborhoods will have similar GraphWave embeddings even though these nodes may reside in very different parts of the network, and our method scales linearly with the number of edges. Experiments in a variety of different settings demonstrate GraphWave's real-world potential for capturing structural roles in networks, and our approach outperforms existing state-of-the-art baselines in every experiment, by as much as 137%.


title: Maximally Recoverable LRCs: A field size lower bound and constructions for few heavy parities

abstract: The explosion in the volumes of data being stored online has resulted in distributed storage systems transitioning to erasure coding based schemes. Local Reconstruction Codes (LRCs) have emerged as the codes of choice for these applications. These codes can correct a small number of erasures by accessing only a small number of remaining coordinates. An $(n,r,h,a,q)$-LRC is a linear code over $\mathbb{F}_q$ of length $n$, whose codeword symbols are partitioned into $g=n/r$ local groups each of size $r$. It has $h$ global parity checks and each local group has $a$ local parity checks. Such an LRC is Maximally Recoverable (MR), if it corrects all erasure patterns which are information-theoretically correctable under the stipulated structure of local and global parity checks. We show the first non-trivial lower bounds on the field size required for MR LRCs. When $a,h$ are constant and the number of local groups $g \ge h$, while $r$ may grow with $n$, our lower bound simplifies to $q\ge \Omega_{a,h}\left(n\cdot r^{\min\{a,h-2\}}\right).$ No superlinear (in $n$) lower bounds were known prior to this work for any setting of parameters. MR LRCs deployed in practice have a small number of global parities, typically $h=2,3$. We complement our lower bounds by giving constructions with small field size for $h\le 3$. For $h=2$, we give a linear field size construction. We also show a surprising application of elliptic curves and arithmetic progression free sets in the construction of MR LRCs.


title: Simulated Encounters of Parker Solar Probe with A Coronal-Hole Jet

abstract: Solar coronal jets are small, transient, collimated ejections most easily observed in coronal holes (CHs). The upcoming Parker Solar Probe (PSP) mission provides the first opportunity to encounter CH jets in situ near the Sun and examine their internal structure and dynamics. Using projected mission orbital parameters, we have simulated PSP encounters with a fully three-dimensional magnetohydrodynamic (MHD) model of a CH jet. We find that three internal jet regions, featuring different wave modes and levels of compressibility, have distinct identifying signatures detectable by PSP. The leading Alfv\'{e}n wave front and its immediate wake are characterized by trans-Alfv\'{e}nic plasma flows with mild density enhancements. This front exhibits characteristics of a fast switch-on MHD shock, whose arrival is signaled by the sudden onset of large-amplitude transverse velocity and magnetic-field oscillations highly correlated in space and time. The trailing portion is characterized by supersonic but sub-Alfv\'{e}nic outflows of dense plasma with uncorrelated velocity and magnetic-field oscillations. This compressible region contains most of the jet's mass. The volume between the immediate wake and dense jet, the remote wake, mixes and transitions the characteristics of the two other regions. In addition to probing each region separately, we also simulate a co-rotational PSP-jet encounter. In this scenario, the simulated spacecraft hovers over the jet-producing CH, as may occur during the mission's co-rotational phases, sampling each jet region in turn. We estimate that PSP will encounter numerous CH jets over the lifetime of the mission.


title: Crystal Graph Convolutional Neural Networks for an Accurate and Interpretable Prediction of Material Properties

abstract: The use of machine learning methods for accelerating the design of crystalline materials usually requires manually constructed feature vectors or complex transformation of atom coordinates to input the crystal structure, which either constrains the model to certain crystal types or makes it difficult to provide chemical insights. Here, we develop a crystal graph convolutional neural networks framework to directly learn material properties from the connection of atoms in the crystal, providing a universal and interpretable representation of crystalline materials. Our method provides a highly accurate prediction of density functional theory calculated properties for eight different properties of crystals with various structure types and compositions after being trained with $10^4$ data points. Further, our framework is interpretable because one can extract the contributions from local chemical environments to global properties. Using an example of perovskites, we show how this information can be utilized to discover empirical rules for materials design.


title: Power Modelling for Heterogeneous Cloud-Edge Data Centers

abstract: Existing power modelling research focuses not on the method used for developing models but rather on the model itself. This paper aims to develop a method for deploying power models on emerging processors that will be used, for example, in cloud-edge data centers. Our research first develops a hardware counter selection method that appropriately selects counters most correlated to power on ARM and Intel processors. Then, we propose a two stage power model that works across multiple architectures. The key results are: (i) the automated hardware performance counter selection method achieves comparable selection to the manual selection methods reported in literature, and (ii) the two stage power model can predict dynamic power more accurately on both ARM and Intel processors when compared to classic power models.


title: Phylogenetic Derivative: A Tool for Assessing Local Tree Reconstruction in the Presence of Recombination

abstract: Recently, much attention has been given to understanding recombination events along a chromosome in a variety of field. For instance, many population genetics problems are limited by the inaccuracy of inferred evolutionary histories of chromosomes sampled randomly from a population. This evolutionary history differs among genomic locations as an artifact of recombination events along a chromosome. Thus, much recent attention has been focused on identifying these recombination points. However, many proposed methods either make simplifying, but unrealistic, assumptions about recombination along a chromosome, or are unable to scale to large genome-wide data like what has become commonplace in statistical genetics. Here, we introduce a \emph{phylogenetic derivative} to describe the relatedness of neighboring trees along a chromosome. This phylogenetic derivative is a computationally efficient, flexible metric that can be also be used assess the prevalence of recombination across a chromosome. These proposed methods are tested and perform well in analyzing both simulated data and a real mouse data set.


title: Allocating and splitting free energy to maximize molecular machine flux

abstract: Biomolecular machines transduce between different forms of energy. These machines make directed progress and increase their speed by consuming free energy, typically in the form of nonequilibrium chemical concentrations. Machine dynamics are often modeled by transitions between a set of discrete metastable conformational states. In general, the free energy change associated with each transition can increase the forward rate constant, decrease the reverse rate constant, or both. In contrast to previous optimizations, we find that in general flux is neither maximized by devoting all free energy changes to increasing forward rate constants nor by solely decreasing reverse rate constants. Instead the optimal free energy splitting depends on the detailed dynamics. Extending our analysis to machines with vulnerable states (from which they can break down), in the strong driving corresponding to in vivo cellular conditions, processivity is maximized by reducing the occupation of the vulnerable state.


title: Revisit Fuzzy Neural Network: Demystifying Batch Normalization and ReLU with Generalized Hamming Network

abstract: We revisit fuzzy neural network with a cornerstone notion of generalized hamming distance, which provides a novel and theoretically justified framework to re-interpret many useful neural network techniques in terms of fuzzy logic. In particular, we conjecture and empirically illustrate that, the celebrated batch normalization (BN) technique actually adapts the normalized bias such that it approximates the rightful bias induced by the generalized hamming distance. Once the due bias is enforced analytically, neither the optimization of bias terms nor the sophisticated batch normalization is needed. Also in the light of generalized hamming distance, the popular rectified linear units (ReLU) can be treated as setting a minimal hamming distance threshold between network inputs and weights. This thresholding scheme, on the one hand, can be improved by introducing double thresholding on both extremes of neuron outputs. On the other hand, ReLUs turn out to be non-essential and can be removed from networks trained for simple tasks like MNIST classification. The proposed generalized hamming network (GHN) as such not only lends itself to rigorous analysis and interpretation within the fuzzy logic theory but also demonstrates fast learning speed, well-controlled behaviour and state-of-the-art performances on a variety of learning tasks.


title: Lower Bounds for Higher-Order Convex Optimization

abstract: State-of-the-art methods in convex and non-convex optimization employ higher-order derivative information, either implicitly or explicitly. We explore the limitations of higher-order optimization and prove that even for convex optimization, a polynomial dependence on the approximation guarantee and higher-order smoothness parameters is necessary. As a special case, we show Nesterov's accelerated cubic regularization method to be nearly tight.


title: Multi-modal Aggregation for Video Classification

abstract: In this paper, we present a solution to Large-Scale Video Classification Challenge (LSVC2017) [1] that ranked the 1st place. We focused on a variety of modalities that cover visual, motion and audio. Also, we visualized the aggregation process to better understand how each modality takes effect. Among the extracted modalities, we found Temporal-Spatial features calculated by 3D convolution quite promising that greatly improved the performance. We attained the official metric mAP 0.8741 on the testing set with the ensemble model.


title: Experimental characterization of the transition to coherence collapse in a semiconductor laser with optical feedback

abstract: Semiconductor lasers with time-delayed optical feedback display a wide range of dynamical regimes, which have found various practical applications. They also provide excellent testbeds for data analysis tools for characterizing complex signals. Recently, several of us have analyzed experimental intensity time-traces and quantitatively identified the onset of different dynamical regimes, as the laser's current increases. Specifically, we identified the onset of low-frequency fluctuations (LFFs), where the laser intensity displays abrupt dropouts, and the onset of coherence collapse (CC), where the intensity fluctuations are highly irregular. Here we map these regimes when both, the laser current and the feedback strength vary. We show that the shape of the distribution of intensity fluctuations (characterized by the standard deviation, the skewness, and the kurtosis) allows to distinguish between noise, LFFs, and CC, and to quantitatively determine (in spite of the gradual nature of the transitions) the boundaries of the three regimes. Ordinal analysis of the inter-dropout time intervals consistently identifies the three regimes occurring in the same parameter regions as the analysis of the intensity distribution. Simulations of the well-known time-delayed Lang-Kobayashi model are in good qualitative agreement with the observations.


title: VO 2 -based radiative thermal transistor in the static regime

abstract: We study a near-field radiative thermal transistor analogous to an electronic one made of a VO 2 base placed between two silica semi-infinite plates playing the roles of the transistor collector and emitter. The fact that VO 2 exhibits an insulator to metal transition is exploited to modulate and/or amplify heat fluxes between the emitter and the collector, by applying a thermal current on the VO 2 base. We study the transistor behavior in 4 typical regimes where the emitter-base and base-collector separation distances can be larger or smaller than the thermal wavelength, and in which the VO 2 layer can be opaque or transparent. Thermal currents variations with the base temperatures are calculated and analyzed. An optimum configuration for base thickness and separation distance maximizing the thermal transistor modulation factor is found.


title: Dissociative recombination by frame transformation to Siegert pseudostates: A~comparison with a numerically solvable model

abstract: We present a simple two-dimensional model of the indirect dissociative recombination process. The model has one electronic and one nuclear degree of freedom and it can be solved to high precision, without making any physically motivated approximations, by employing the exterior complex scaling method together with the finite-elements method and discrete variable representation. The approach is applied to solve a model for dissociative recombination of H$_2^+$ and the results serve as a benchmark to test validity of several physical approximations commonly used in the computational modeling of dissociative recombination for real molecular targets. The second, approximate, set of calculations employs a combination of multi-channel quantum defect theory and frame transformation into a basis of Siegert pseudostates. The cross sections computed with the two methods are compared in detail for collision energies from 0 to 2 eV.


title: Automatic estimation of attractor invariants

abstract: The invariants of an attractor have been the most used resource to characterize a nonlinear dynamics. Their estimation is a challenging endeavor in short-time series and/or in presence of noise. In this article we present two new coarse-grained estimators for the correlation dimension and for the correlation entropy. They can be easily estimated from the calculation of two U-correlation integrals. We have also developed an algorithm that is able to automatically obtain these invariants and the noise level in order to process large data sets. This method has been statistically tested through simulations in low-dimensional systems. The results show that it is robust in presence of noise and short data lengths. In comparison with similar approaches, our algorithm outperforms the estimation of the correlation entropy.


title: Similarity-based Multi-label Learning

abstract: Multi-label classification is an important learning problem with many applications. In this work, we propose a principled similarity-based approach for multi-label learning called SML. We also introduce a similarity-based approach for predicting the label set size. The experimental results demonstrate the effectiveness of SML for multi-label classification where it is shown to compare favorably with a wide variety of existing algorithms across a range of evaluation criterion.


title: Modeling and Real-Time Scheduling of DC Platform Supply Vessel for Fuel Efficient Operation

abstract: DC marine architecture integrated with variable speed diesel generators (DGs) has garnered the attention of the researchers primarily because of its ability to deliver fuel efficient operation. This paper aims in modeling and to autonomously perform real-time load scheduling of dc platform supply vessel (PSV) with an objective to minimize specific fuel oil consumption (SFOC) for better fuel efficiency. Focus has been on the modeling of various components and control routines, which are envisaged to be an integral part of dc PSVs. Integration with photovoltaic-based energy storage system (ESS) has been considered as an option to cater for the short time load transients. In this context, this paper proposes a real-time transient simulation scheme, which comprises of optimized generation scheduling of generators and ESS using dc optimal power flow algorithm. This framework considers real dynamics of dc PSV during various marine operations with possible contingency scenarios, such as outage of generation systems, abrupt load changes, and unavailability of ESS. The proposed modeling and control routines with real-time transient simulation scheme have been validated utilizing the real-time marine simulation platform. The results indicate that the coordinated treatment of renewable based ESS with DGs operating with optimized speed yields better fuel savings. This has been observed in improved SFOC operating trajectory for critical marine missions. Furthermore, SFOC minimization at multiple suboptimal points with its treatment in the real-time marine system is also highlighted.


title: Evidence for the Confinement of Magnetic Monopoles in Quantum Spin Ice

abstract: Magnetic monopoles are hypothesised elementary particles connected by Dirac strings that behave like infinitely thin solenoids. Despite decades of searches, free magnetic monopoles and their Dirac strings have eluded experimental detection, although there is substantial evidence for deconfined magnetic monopole quasiparticles in spin ice materials. Here we report the detection of a hierarchy of unequally-spaced magnetic excitations \emph{via} high resolution inelastic neutron spectroscopic measurements on the quantum spin ice candidate Pr$_{2}$Sn$_{2}$O$_{7}$. These excitations are well-described by a simple model of monopole pairs bound by a linear potential with an effective tension of 0.642(8) K~$\cdot$\AA$^{-1}$ at 1.65~K. The success of the linear potential model suggests that these low energy magnetic excitations are direct spectroscopic evidence for the confinement of magnetic monopole quasiparticles in the quantum spin ice candidate Pr$_{2}$Sn$_{2}$O$_{7}$.


title: E-beam manipulation of Si atoms on graphene edges with aberration-corrected STEM

abstract: The burgeoning field of atomic level material control holds great promise for future breakthroughs in quantum and memristive device manufacture and fundamental studies of atomic-scale chemistry. Realization of atom-by atom control of matter represents a complex and ongoing challenge. Here, we explore the feasibility of controllable motion of dopant Si atoms at the edges of graphene via the sub-atomically focused electron beam in a scanning transmission electron microscope (STEM). We demonstrate that the graphene edges can be cleaned of Si atoms and then subsequently replenished from nearby source material. It is also shown how Si edge atoms may be pushed from the edge of a small hole into the bulk of the graphene lattice and from the bulk of the lattice back to the edge. This is accomplished through sputtering of the edge of the graphene lattice to bury or uncover Si dopant atoms. These experiments form an initial step toward general atomic scale material control.


title: Improved approximation of layout problems on random graphs

abstract: Inspired by previous work of Diaz, Petit, Serna, and Trevisan (Approximating layout problems on random graphs Discrete Mathematics, 235, 2001, 245--253), we show that several well-known graph layout problems are approximable to within a factor arbitrarily close to 1 of the optimal with high probability for random graphs drawn from an Erd\"os-Renyi distribution with appropriate sparsity conditions. Moreover, we show that the same results hold for the analogous problems on directed acyclic graphs.


title: Photobleaching of randomly-rotating fluorescently-decorated particles

abstract: Randomly rotating particles that have been isotropically labeled with rigidly linked fluorophores will undergo non-isotropic (patchy) photobleaching under illumination due to the dipole coupling of fluorophores with light. For a rotational diffusion rate $D$ of the particle and a photobleaching timescale $\tau$ of the fluorophores, the dynamics of this process are characterized by the dimensionless combination $D \tau$. We find significant interparticle fluctuations at intermediate $D \tau$. These fluctuations vanish at both large and small $D \tau$, or at small or large elapsed times $t$. Associated with these fluctuations between particles, we also observe transient non-monotonicities of the brightness of individual particles. These non-monotonicities can be as much 20\% of the original brightness. We show that these novel photobleach-fluctuations dominate over variability of single-fluorophore orientation when there are at least $10^3$ fluorophores on individual particles.


title: Single shot all optical switching of intrinsic micron size magnetic domains of a Pt/Co/Pt ferromagnetic stack

abstract: We demonstrate that the magnetization reversal in a ferromagnetic Pt/Co/Pt stack can be induced by a single femtosecond laser pulse. We find that the size of the switched spot is comparable to the size of the intrinsic magnetic domains. It requires an absorbed energy density of $\sim$4 mJ.cm$^{-2}$, beyond which the excited spot fragments into a multidomain structure. The switching process occurs back and forth with subsequent laser pulses and it is helicity-independent. Furthermore, the dynamics of the magnetization reversal occurs in a timescale less than one microsecond. These results suggest that all optical switching in ferromagnetic films requires to match the laser spot with the specific domain sizes.


title: Insights on Variance Estimation for Blocked and Matched Pairs Designs

abstract: Evaluating blocked randomized experiments from a potential outcomes perspective has two primary branches of work. The first focuses on larger blocks, with multiple treatment and control units in each block. The second focuses on matched pairs, with a single treatment and control unit in each block. These literatures not only provide different estimators for the standard errors of the estimated average impact, but they are also built on different sets of assumptions. Neither literature handles cases with blocks of varying size that contain singleton treatment or control units, a case which can occur in a variety of contexts, such as with different forms of matching or post-stratification. In this paper, we reconcile the literatures by carefully examining the performance of variance estimators under several different frameworks. We then use these insights to derive novel variance estimators for experiments containing blocks of different sizes.


title: Hasse diagrams of posets with up to 7 elements, and the number of posets with 10 elements, without the use of computer programs

abstract: Let $P(n)$ be the set of all posets with $n$ elements. Let $P^{(j)}(n)$, $1\leq j\leq 2^n,$ be the number of all posets with $n$ elements possessing exactly $j$ antichains. We have determined the numbers $P^{(j)}(7),$ $1\leq j\leq 128$, and using a result of M.~Ern\'e [Ern\'e, M., On the cardinalities of finite topologies and the number of antichains in partially ordered sets, Discrete Mathematics 35 (1981), 119-133.], we compute $|P(10)|$ without the aid of any computer program. We include the Hasse diagrams of all the non-isomorphic posets of $P(7)$. We also present formulas for the number of connected posets of certain forms, and use them to compute $|P(n)|$ with $1\le n\le 8$ by a different method.


title: A Treatise on Sucker's Bets

abstract: In 1970, Statistics giant, Bradley Efron, amazed the world by coming up with a set of four dice, let's call them A,B,C,D, whose faces are marked with [0,0,4,4,4,4], [3,3,3,3,3,3],[2,2,2,2,6,6],[1,1,1,5,5,5] respectively, where die A beats die B, die B beats die C, die C beats die D, but, surprise surprise, die D beats die A! This was an amazing demonstration that "being more likely to win" is not a transitive relation. But that was only one example, and of course, instead of dice, we can use decks of cards, where they are called (by Martin Gardner, who popularized this way back in 1970) , "sucker's bets". Can you find all such examples, with a specified number of decks, and sizes? If you have a computer algebra system (in our case Maple), you sure can! Not only that, we can figure out how likely such sucker bets are, and derive, fully automatically, statistical information!


title: The Implicit Bias of Gradient Descent on Separable Data

abstract: We examine gradient descent on unregularized logistic regression problems, with homogeneous linear predictors on linearly separable datasets. We show the predictor converges to the direction of the max-margin (hard margin SVM) solution. The result also generalizes to other monotone decreasing loss functions with an infimum at infinity, to multi-class problems, and to training a weight layer in a deep network in a certain restricted setting. Furthermore, we show this convergence is very slow, and only logarithmic in the convergence of the loss itself. This can help explain the benefit of continuing to optimize the logistic or cross-entropy loss even after the training error is zero and the training loss is extremely small, and, as we show, even if the validation loss increases. Our methodology can also aid in understanding implicit regularization n more complex models and with other optimization methods.


title: Inference for stochastic kinetic models from multiple data sources for joint estimation of infection dynamics from aggregate reports and virological data

abstract: Before the current pandemic, influenza and respiratory syncytial virus (RSV) were the leading etiological agents of seasonal acute respiratory infections (ARI) around the world. In this setting, medical doctors typically based the diagnosis of ARI on patients' symptoms alone and did not routinely conduct virological tests necessary to identify individual viruses, limiting the ability to study the interaction between multiple pathogens and to make public health recommendations. We consider a stochastic kinetic model (SKM) for two interacting ARI pathogens circulating in a large population and an empirically-motivated background process for infections with other pathogens causing similar symptoms. An extended marginal sampling approach, based on the linear noise approximation to the SKM, integrates multiple data sources and additional model components. We infer the parameters defining the pathogens' dynamics and interaction within a Bayesian model and explore the posterior trajectories of infections for each illness based on aggregate infection reports from six epidemic seasons collected by the state health department and a subset of virological tests from a sentinel program at a general hospital in San Luis Potos\'{i}, M\'{e}xico. We interpret the results and make recommendations for future data collection strategies.


title: Diameter and curvature control under mean curvature flow

abstract: We prove that for the mean curvature flow of two-convex hypersurfaces the intrinsic diameter stays uniformly controlled as one approaches the first singular time. We also derive sharp $L^{n-1}$-estimates for the regularity scale of the level set flow with two-convex initial data. Our proof relies on a detailed analysis of cylindrical regions ($\varepsilon$-tubes) under mean curvature flow. The results are new even in the most classical case of mean convex surfaces evolving by mean curvature flow in $\mathbb{R}^3$.


title: Multi-level Residual Networks from Dynamical Systems View

abstract: Deep residual networks (ResNets) and their variants are widely used in many computer vision applications and natural language processing tasks. However, the theoretical principles for designing and training ResNets are still not fully understood. Recently, several points of view have emerged to try to interpret ResNet theoretically, such as unraveled view, unrolled iterative estimation and dynamical systems view. In this paper, we adopt the dynamical systems point of view, and analyze the lesioning properties of ResNet both theoretically and experimentally. Based on these analyses, we additionally propose a novel method for accelerating ResNet training. We apply the proposed method to train ResNets and Wide ResNets for three image classification benchmarks, reducing training time by more than 40% with superior or on-par accuracy.


title: Sharp estimates for oscillatory integral operators via polynomial partitioning

abstract: The sharp range of $L^p$-estimates for the class of H\"ormander-type oscillatory integral operators is established in all dimensions under a positive-definite assumption on the phase. This is achieved by generalising a recent approach of the first author for studying the Fourier extension operator, which utilises polynomial partitioning arguments.


title: Real-Time Robust Finger Gaits Planning under Object Shape and Dynamics Uncertainties

abstract: Dexterous manipulation has broad applications in assembly lines, warehouses and agriculture. To perform large-scale manipulation tasks for various objects, a multi-fingered robotic hand sometimes has to sequentially adjust its grasping gestures, i.e. the finger gaits, to address the workspace limits and guarantee the object stability. However, realizing finger gaits planning in dexterous manipulation is challenging due to the complicated grasp quality metrics, uncertainties on object shapes and dynamics (mass and moment of inertia), and unexpected slippage under uncertain contact dynamics. In this paper, a dual-stage optimization based planner is proposed to handle these challenges. In the first stage, a velocity-level finger gaits planner is introduced by combining object grasp quality with hand manipulability. The proposed finger gaits planner is computationally efficient and realizes finger gaiting without 3D model of the object. In the second stage, a robust manipulation controller using robust control and force optimization is proposed to address object dynamics uncertainties and external disturbances. The dual-stage planner is able to guarantee stability under unexpected slippage caused by uncertain contact dynamics. Moreover, it does not require velocity measurement or expensive 3D/6D tactile sensors. The proposed dual-stage optimization based planner is verified by simulations on Mujoco.


title: Bayesian Spatial Binary Regression for Label Fusion in Structural Neuroimaging

abstract: Alzheimer's disease is a neurodegenerative condition that accelerates cognitive decline relative to normal aging. It is of critical scientific importance to gain a better understanding of early disease mechanisms in the brain to facilitate effective, targeted therapies. The volume of the hippocampus is often used in diagnosis and monitoring of the disease. Measuring this volume via neuroimaging is difficult since each hippocampus must either be manually identified or automatically delineated, a task referred to as segmentation. Automatic hippocampal segmentation often involves mapping a previously manually segmented image to a new brain image and propagating the labels to obtain an estimate of where each hippocampus is located in the new image. A more recent approach to this problem is to propagate labels from multiple manually segmented atlases and combine the results using a process known as label fusion. To date, most label fusion algorithms employ voting procedures with voting weights assigned directly or estimated via optimization. We propose using a fully Bayesian spatial regression model for label fusion that facilitates direct incorporation of covariate information while making accessible the entire posterior distribution. Our results suggest that incorporating tissue classification (e.g, gray matter) into the label fusion procedure can greatly improve segmentation when relatively homogeneous, healthy brains are used as atlases for diseased brains. The fully Bayesian approach also produces meaningful uncertainty measures about hippocampal volumes, information which can be leveraged to detect significant, scientifically meaningful differences between healthy and diseased populations, improving the potential for early detection and tracking of the disease.


title: Automated Design using Neural Networks and Gradient Descent

abstract: We propose a novel method that makes use of deep neural networks and gradient decent to perform automated design on complex real world engineering tasks. Our approach works by training a neural network to mimic the fitness function of a design optimization task and then, using the differential nature of the neural network, perform gradient decent to maximize the fitness. We demonstrate this methods effectiveness by designing an optimized heat sink and both 2D and 3D airfoils that maximize the lift drag ratio under steady state flow conditions. We highlight that our method has two distinct benefits over other automated design approaches. First, evaluating the neural networks prediction of fitness can be orders of magnitude faster then simulating the system of interest. Second, using gradient decent allows the design space to be searched much more efficiently then other gradient free methods. These two strengths work together to overcome some of the current shortcomings of automated design.


title: A Novikov fundamental group

abstract: Given a $1$-cohomology class $u$ on a closed manifold $M$, we define a Novikov fundamental group associated to $u$, generalizing the usual fundamental group in the same spirit as Novikov homology generalizes Morse homology to the case of non exact $1$-forms. As an application, lower bounds for the minimal number of index $1$ and $2$ critical points of Morse closed $1$-forms are obtained, that are different in nature from those derived from the Novikov homology.


title: Anderson localization and Brewster anomaly of electromagnetic waves in randomly-stratified anisotropic media

abstract: Anderson localization of $p$-polarized waves and the Brewster anomaly phenomenon, which is the delocalization of $p$-polarized waves at a special incident angle, in randomly-stratified anisotropic media are studied theoretically for two different random models. In the first model, the random parts of the transverse and longitudinal components of the dielectric tensor, between which the longitudinal component is the one in the stratification direction, are assumed to be uncorrelated, while, in the second model, they are proportional to each other. We calculate the localization length in a precise way using the invariant imbedding method. From analytical considerations, we provide an interpretation of the Brewster anomaly as a phenomenon arising when the wave impedance is effectively uniform. Similarly, the ordinary Brewster effect is interpreted as an impedance matching phenomenon. We derive the existence condition for the Brewster anomaly and concise analytical expressions for the localization length, which are accurate in the weak disorder regime. We find that the Brewster anomaly can arise only when disorder is sufficiently weak and only in the second model with a positive ratio of the random parts. The incident angle at which the anomaly occurs depends sensitively on the ratio of the random parts and the average values of the tensor components. In the cases where the critical angle of total reflection exists, the angle at which the anomaly occurs can be either bigger or smaller than the critical angle. When the transverse and longitudinal components are uncorrelated, localization is dominated by the the transverse component at small incident angles. When only the longitudinal component is random, the localization length diverges as $\theta^{-4}$ as the incident angle $\theta$ goes to zero and is also argued to diverge for all $\theta$ in the strong disorder limit.


title: Convolutional Neural Networks Via Node-Varying Graph Filters

abstract: Convolutional neural networks (CNNs) are being applied to an increasing number of problems and fields due to their superior performance in classification and regression tasks. Since two of the key operations that CNNs implement are convolution and pooling, this type of networks is implicitly designed to act on data described by regular structures such as images. Motivated by the recent interest in processing signals defined in irregular domains, we advocate a CNN architecture that operates on signals supported on graphs. The proposed design replaces the classical convolution not with a node-invariant graph filter (GF), which is the natural generalization of convolution to graph domains, but with a node-varying GF. This filter extracts different local features without increasing the output dimension of each layer and, as a result, bypasses the need for a pooling stage while involving only local operations. A second contribution is to replace the node-varying GF with a hybrid node-varying GF, which is a new type of GF introduced in this paper. While the alternative architecture can still be run locally without requiring a pooling stage, the number of trainable parameters is smaller and can be rendered independent of the data dimension. Tests are run on a synthetic source localization problem and on the 20NEWS dataset.


title: Optimal Control of Wireless Computing Networks

abstract: Augmented information (AgI) services allow users to consume information that results from the execution of a chain of service functions that process source information to create real-time augmented value. Applications include real-time analysis of remote sensing data, real-time computer vision, personalized video streaming, and augmented reality, among others. We consider the problem of optimal distribution of AgI services over a wireless computing network, in which nodes are equipped with both communication and computing resources. We characterize the wireless computing network capacity region and design a joint flow scheduling and resource allocation algorithm that stabilizes the underlying queuing system while achieving a network cost arbitrarily close to the minimum, with a tradeoff in network delay. Our solution captures the unique chaining and flow scaling aspects of AgI services, while exploiting the use of the broadcast approach coding scheme over the wireless channel.


title: Quantum effects of Aharonov-Bohm type and noncommutative quantum mechanics

abstract: Quantum mechanics in noncommutative space modifies the standard result of the Aharonov-Bohm effect for electrons and other recent quantum effects. Here we obtain the phase in noncommutative space for the Spavieri effect, a generalization of Aharonov-bohm effect which involving a coherent superposition of particles with opposite charges moving along single open interferometric path. By means of the experimental considerations a limit $\sqrt{\theta}\simeq\left( 0,13T\text{eV}\right) ^{-1}$ is achieved, improving by 10 orders of magnitude the derived by Chaichian et. al. for the Aharonov-Bohm effect. It is also shown that the noncommutative phases of the Aharonov-Casher and He-McKellar-Willkens effects are nullified in the current experimental tests.


title: Prediction of Topological Crystalline Insulator and Topological Phase Transitions in Two-dimensional PbTe Films

abstract: Topological phases, especially topological crystalline insulators (TCIs), have been intensively explored observed experimentally in three-dimensional (3D) materials. However, the two-dimensional (2D) films are explored much less than 3D TCI, and even 2D topological insulators. Based on ab initio calculations, here we investigate the electronic and topological properties of 2D PbTe(001) few-layers. The monolayer and trilayer PbTe are both intrinsic 2D TCIs with a large band gap reaching 0.27 eV, indicating a high possibility for room-temperature observation of quantized conductance. The origin of TCI phase can be attributed to the p band inversion,which is determined by the competitions of orbital hybridization and quantum confinement. We also observe a semimetal-TCI-normal insulator transition under biaxial strains, whereas a uniaxial strains lead to Z2 nontrivial states. Especially, the TCI phase of PbTe monolayer remains when epitaxial grow on NaI semiconductor substrate. Our findings on the controllable quantum states with sizable band gaps present an ideal platform for realizing future topological quantum devices with ultralow dissipation.


title: Another look at the impact of an eV-mass sterile neutrino on the effective neutrino mass of neutrinoless double-beta decays

abstract: The possible existence of an eV-mass sterile neutrino, slightly mixing with ordinary active neutrinos, is not yet excluded by neutrino oscillation experiments. Assuming neutrinos to be Majorana particles, we explore the impact of such a sterile neutrino on the effective neutrino mass of neutrinoless double-beta decays $\langle m \rangle^\prime_{ee} \equiv m^{}_1 |V^{}_{e1}|^2 e^{{\rm i}\rho} + m^{}_2 |V^{}_{e2}|^2 + m^{}_3 |V^{}_{e3}|^2 e^{{\rm i}\sigma} + m^{}_4 |V^{}_{e4}|^2 e^{{\rm i}\omega}$, where $m^{}_i$ and $V^{}_{ei}$ (for $i = 1, 2, 3, 4$) denote respectively the absolute masses and the first-row elements of the 4$\times$4 neutrino flavor mixing matrix $V$, for which a full parametrization involves three Majorana-type CP-violating phases $\{\rho, \sigma, \omega\}$. A zero effective neutrino mass $|\langle m \rangle^\prime_{ee}| = 0$ is possible no matter whether three active neutrinos take the normal or inverted mass ordering, and its implications for the parameter space are examined in great detail. In particular, given the best-fit values of $m^{}_4 \approx 1.3~{\rm eV}$ and $|V^{}_{e4}|^2 \approx 0.019$ from the latest global analysis of neutrino oscillation data, a three-dimensional view of $|\langle m \rangle^\prime_{ee}|$ in the $(m^{}_1, \rho)$-plane is presented and further compared with that of the counterpart $|\langle m \rangle^{}_{ee}|$ in the absence of any sterile neutrino.


title: Rotation of inertial frames by angular momentum of matter and waves

abstract: We elucidate the dynamics of a thin spherical material shell with a tangential pressure, using a new approach. This is both simpler than the traditional method of extrinsic curvature junction conditions (which we also employ), and suggests an expression for a `gravitational potential energy' of the shell. Such a shell, if slowly spinning, can rotationally drag the inertial frames within it through a finite angle as it collapses and rebounds from a minimum radius. Rebounding `spherical' and cylindrical pulses of rotating gravitational waves were studied previously. Here we calculate their angular momentum and show that their rotational frame dragging is in agreement with that of the rotating spherical shell and a rotating cylindrical dust shell. This shows that Machian effects occur equally for material and analogous `immaterial' sources.


title: Deep Residual Learning for Small-Footprint Keyword Spotting

abstract: We explore the application of deep residual learning and dilated convolutions to the keyword spotting task, using the recently-released Google Speech Commands Dataset as our benchmark. Our best residual network (ResNet) implementation significantly outperforms Google's previous convolutional neural networks in terms of accuracy. By varying model depth and width, we can achieve compact models that also outperform previous small-footprint variants. To our knowledge, we are the first to examine these approaches for keyword spotting, and our results establish an open-source state-of-the-art reference to support the development of future speech-based interfaces.


title: Bandlimited approximations and estimates for the Riemann zeta-function

abstract: In this paper, we provide explicit upper and lower bounds for the argument of the Riemann zeta-function and its antiderivatives in the critical strip under the assumption of the Riemann hypothesis. This extends the previously known bounds for these quantities on the critical line (and sharpens the error terms in such estimates). Our tools come not only from number theory, but also from Fourier analysis and approximation theory. An important element in our strategy is the ability to solve a Fourier optimization problem with constraints, namely, the problem of majorizing certain real-valued even functions by bandlimited functions, optimizing the $L^1(\mathbb{R})-$error. Deriving explicit formulae for the Fourier transforms of such optimal approximations plays a crucial role in our approach.


title: Diff-DAC: Distributed Actor-Critic for Average Multitask Deep Reinforcement Learning

abstract: We propose a fully distributed actor-critic algorithm approximated by deep neural networks, named \textit{Diff-DAC}, with application to single-task and to average multitask reinforcement learning (MRL). Each agent has access to data from its local task only, but it aims to learn a policy that performs well on average for the whole set of tasks. During the learning process, agents communicate their value-policy parameters to their neighbors, diffusing the information across the network, so that they converge to a common policy, with no need for a central node. The method is scalable, since the computational and communication costs per agent grow with its number of neighbors. We derive Diff-DAC's from duality theory and provide novel insights into the standard actor-critic framework, showing that it is actually an instance of the dual ascent method that approximates the solution of a linear program. Experiments suggest that Diff-DAC can outperform the single previous distributed MRL approach (i.e., Dist-MTLPS) and even the centralized architecture.


title: Consistency of Lipschitz learning with infinite unlabeled data and finite labeled data

abstract: We study the consistency of Lipschitz learning on graphs in the limit of infinite unlabeled data and finite labeled data. Previous work has conjectured that Lipschitz learning is well-posed in this limit, but is insensitive to the distribution of the unlabeled data, which is undesirable for semi-supervised learning. We first prove that this conjecture is true in the special case of a random geometric graph model with kernel-based weights. Then we go on to show that on a random geometric graph with self-tuning weights, Lipschitz learning is in fact highly sensitive to the distribution of the unlabeled data, and we show how the degree of sensitivity can be adjusted by tuning the weights. In both cases, our results follow from showing that the sequence of learned functions converges to the viscosity solution of an $\infty$-Laplace type equation, and studying the structure of the limiting equation.


title: Sharp mixed norm spherical restriction

abstract: Let $d\geq 2$ be an integer and let $2d/(d-1) < q \leq \infty$. In this paper we investigate the sharp form of the mixed norm Fourier extension inequality \begin{equation*} \big\|\widehat{f\sigma}\big\|_{L^q_{{\rm rad}}L^2_{{\rm ang}}(\mathbb{R}^d)} \leq {\bf C}_{d,q}\, \|f\|_{L^2(\mathbb{S}^{d-1},{\rm d}\sigma)}, \end{equation*} established by L. Vega in 1988. Letting $\mathcal{A}_d \subset (2d/(d-1), \infty]$ be the set of exponents for which the constant functions on $\mathbb{S}^{d-1}$ are the unique extremizers of this inequality, we show that: (i) $\mathcal{A}_d$ contains the even integers and $\infty$; (ii) $\mathcal{A}_d$ is an open set in the extended topology; (iii) $\mathcal{A}_d$ contains a neighborhood of infinity $(q_0(d), \infty]$ with $q_0(d) \leq \left(\tfrac{1}{2} + o(1)\right) d\log d$. In low dimensions we show that $q_0(2) \leq 6.76\,;\,q_0(3) \leq 5.45 \,;\, q_0(4) \leq 5.53 \,;\, q_0(5) \leq 6.07$. In particular, this breaks for the first time the even exponent barrier in sharp Fourier restriction theory. The crux of the matter in our approach is to establish a hierarchy between certain weighted norms of Bessel functions, a nontrivial question of independent interest within the theory of special functions.


title: Lower Bounds for Two-Sample Structural Change Detection in Ising and Gaussian Models

abstract: The change detection problem is to determine if the Markov network structures of two Markov random fields differ from one another given two sets of samples drawn from the respective underlying distributions. We study the trade-off between the sample sizes and the reliability of change detection, measured as a minimax risk, for the important cases of the Ising models and the Gaussian Markov random fields restricted to the models which have network structures with $p$ nodes and degree at most $d$, and obtain information-theoretic lower bounds for reliable change detection over these models. We show that for the Ising model, $\Omega\left(\frac{d^2}{(\log d)^2}\log p\right)$ samples are required from each dataset to detect even the sparsest possible changes, and that for the Gaussian, $\Omega\left( \gamma^{-2} \log(p)\right)$ samples are required from each dataset to detect change, where $\gamma$ is the smallest ratio of off-diagonal to diagonal terms in the precision matrices of the distributions. These bounds are compared to the corresponding results in structure learning, and closely match them under mild conditions on the model parameters. Thus, our change detection bounds inherit partial tightness from the structure learning schemes in previous literature, demonstrating that in certain parameter regimes, the naive structure learning based approach to change detection is minimax optimal up to constant factors.


title: Deep crustal heating by neutrinos from the surface of accreting neutron stars

abstract: We present a new mechanism for deep crustal heating in accreting neutron stars. Charged pions ($\pi^+$) are produced in nuclear collisions on the neutron star surface during active accretion and upon decay they provide a flux of neutrinos into the neutron star crust. For massive and/or compact neutron stars, neutrinos deposit $\approx 1\textrm{--} 2 \, \mathrm{MeV}$ of heat per accreted nucleon into the inner crust. The strength of neutrino heating is comparable to the previously known sources of deep crustal heating, such as from pycnonuclear fusion reactions, and is relevant for studies of cooling neutron stars. We model the thermal evolution of a transient neutron star in a low-mass X-ray binary, and in the particular case of the neutron star MXB~1659-29 we show that additional deep crustal heating requires a higher thermal conductivity for the neutron star inner crust. A better knowledge of pion production cross sections near threshold would improve the accuracy of our predictions.


title: Deep Generative Dual Memory Network for Continual Learning

abstract: Despite advances in deep learning, neural networks can only learn multiple tasks when trained on them jointly. When tasks arrive sequentially, they lose performance on previously learnt tasks. This phenomenon called catastrophic forgetting is a fundamental challenge to overcome before neural networks can learn continually from incoming data. In this work, we derive inspiration from human memory to develop an architecture capable of learning continuously from sequentially incoming tasks, while averting catastrophic forgetting. Specifically, our contributions are: (i) a dual memory architecture emulating the complementary learning systems (hippocampus and the neocortex) in the human brain, (ii) memory consolidation via generative replay of past experiences, (iii) demonstrating advantages of generative replay and dual memories via experiments, and (iv) improved performance retention on challenging tasks even for low capacity models. Our architecture displays many characteristics of the mammalian memory and provides insights on the connection between sleep and learning.


title: Coupled breathing modes in one-dimensional Skyrmion lattices

abstract: We explored strong coupling of dynamic breathing modes in one-dimensional (1D) skyrmion lattices periodically arranged in thin-film nanostrips. The coupled breathing modes exhibit characteristic concave-down dispersions that represent the in-phase high-energy mode at zero wavenumber (k=0) and the anti-phase low-energy mode at the Brillouin zone boundary (k=kBZ). The band width of the allowed modes increases with decreasing inter-distance between nearest-neighboring skyrmions. Furthermore, the collective breathing modes propagate very well through the thin-film nanostrips, as fast as 200 ~ 700 m/s, which propagation is controllable by the strength of magnetic fields applied perpendicularly to the film plane. The breathing modes in 1D skyrmion lattices potentially formed in such nanostrips possibly can be used as information carriers in information processing devices.


title: Topology Adaptive Graph Convolutional Networks

abstract: Spectral graph convolutional neural networks (CNNs) require approximation to the convolution to alleviate the computational complexity, resulting in performance loss. This paper proposes the topology adaptive graph convolutional network (TAGCN), a novel graph convolutional network defined in the vertex domain. We provide a systematic way to design a set of fixed-size learnable filters to perform convolutions on graphs. The topologies of these filters are adaptive to the topology of the graph when they scan the graph to perform convolution. The TAGCN not only inherits the properties of convolutions in CNN for grid-structured data, but it is also consistent with convolution as defined in graph signal processing. Since no approximation to the convolution is needed, TAGCN exhibits better performance than existing spectral CNNs on a number of data sets and is also computationally simpler than other recent methods.


title: Renormalization Group Invariance of the Pole Mass in the Multi-Higgs system

abstract: We have investigated the renormalization group running of the pole mass in the multi-Higgs theory in two different types of the gauge fixing conditions. It turns out that the pole mass when expressed in terms of the Lagrangian parameters, is invariant under the renormalization group with the beta and gamma functions of the symmetric phase.


title: Intermediate-Mass-Elements in Young Supernova Remnants Reveal Neutron Star Kicks by Asymmetric Explosions

abstract: The birth properties of neutron stars yield important information on the still debated physical processes that trigger the explosion and on intrinsic neutron-star physics. These properties include the high space velocities of young neutron stars with average values of several 100 km/s, whose underlying "kick" mechanism is not finally clarified. There are two competing possibilities that could accelerate NSs during their birth: anisotropic ejection of either stellar debris or neutrinos. We here present new evidence from X-ray measurements that chemical elements between silicon and calcium in six young gaseous supernova remnants are preferentially expelled opposite to the direction of neutron star motion. There is no correlation between the kick velocities and magnetic field strengths of these neutron stars. Our results support a hydrodynamic origin of neutron-star kicks connected to asymmetric explosive mass ejection, and they conflict with neutron-star acceleration scenarios that invoke anisotropic neutrino emission caused by particle and nuclear physics in combination with very strong neutron-star magnetic fields.


title: Combinatorial proof of an identity of Andrews and Yee

abstract: Recently, Andrews and Yee studied two-variable generalizations of two identities involving partition functions $p_\omega(n)$ and $p_\nu(n)$ introduced by Andrews, Dixit and Yee. In this paper, we present a combinatorial proof of an interesting identity in their work.


title: Monotone Covering Properties defined by Closure-Preserving Operators

abstract: We continue Gartside, Moody, and Stares' study of versions of monotone paracompactness. We show that the class of spaces with a monotone closure-preserving open operator is strictly larger than those with a monotone open locally-finite operator. We prove that monotonically metacompact GO-spaces have a monotone open locally-finite operator, and so do GO-spaces with a monotone (open or not) closure-preserving operator, whose underlying LOTS has a $\sigma$-closed-discrete dense subset. A GO-space with a $\sigma$-closed-discrete dense subset and a monotone closure-preserving operator is metrizable. A compact LOTS with a monotone open closure-preserving operator is metrizable.


title: The $q$-Schur algebras and $q$-Schur dualities of finite type

abstract: We formulate a $q$-Schur algebra associated to an arbitrary $W$-invariant finite set $X_{\texttt f}$ of integral weights for a complex simple Lie algebra with Weyl group $W$. We establish a $q$-Schur duality between the $q$-Schur algebra and Hecke algebra associated to $W$. We then realize geometrically the $q$-Schur algebra and duality, and construct a canonical basis for the $q$-Schur algebra with positivity. With suitable choices of $X_{\texttt f}$ in classical types, we recover the $q$-Schur algebras in the literature. Our $q$-Schur algebras are closely related to the category $\mathcal O$, where the type $G_2$ is studied in detail.


title: Searching for Axion Dark Matter using Radio Telescopes

abstract: We investigate the use of next generation radio telescopes such as the Square Kilometre Array (SKA) to detect axion two-photon coupling in the astrophysical environment. The uncertainty surrounding astrophysical magnetic fields presents new challenges, but with a frequency range corresponding to axions of mass $1.7-57\mu$eV and a spectral profile with a number of distinguishing features, SKA-mid offers a tantalising opportunity to constrain axion dark matter properties. To determine the sensitivity of SKA-mid to an axion signal, we consider observations of the Galactic centre and interstellar medium, and find that this new telescope could allow us to probe axion couplings $\gtrsim10^{-16}$GeV$^{-1}$.


title: Quantum attacks on Bitcoin, and how to protect against them

abstract: The key cryptographic protocols used to secure the internet and financial transactions of today are all susceptible to attack by the development of a sufficiently large quantum computer. One particular area at risk are cryptocurrencies, a market currently worth over 150 billion USD. We investigate the risk of Bitcoin, and other cryptocurrencies, to attacks by quantum computers. We find that the proof-of-work used by Bitcoin is relatively resistant to substantial speedup by quantum computers in the next 10 years, mainly because specialized ASIC miners are extremely fast compared to the estimated clock speed of near-term quantum computers. On the other hand, the elliptic curve signature scheme used by Bitcoin is much more at risk, and could be completely broken by a quantum computer as early as 2027, by the most optimistic estimates. We analyze an alternative proof-of-work called Momentum, based on finding collisions in a hash function, that is even more resistant to speedup by a quantum computer. We also review the available post-quantum signature schemes to see which one would best meet the security and efficiency requirements of blockchain applications.


title: Distributed Change Detection via Average Consensus over Networks

abstract: Distributed change-point detection has been a fundamental problem when performing real-time monitoring using sensor-networks. We propose a distributed detection algorithm, where each sensor only exchanges CUSUM statistic with their neighbors based on the average consensus scheme, and an alarm is raised when local consensus statistic exceeds a pre-specified global threshold. We provide theoretical performance bounds showing that the performance of the fully distributed scheme can match the centralized algorithms under some mild conditions. Numerical experiments demonstrate the good performance of the algorithm especially in detecting asynchronous changes.


title: Almost-global tracking of the unactuated joint in a pendubot

abstract: Tracking the unactuated configuration variable in an underactuated system, in a global sense, has not received much attention. Here we present a scheme to do so for a pendubot - a two link robot actuated only at the first link. We propose a control law for almost-global asymptotic tracking (AGAT) of a smooth reference trajectory for the unactuated second joint of the pendubot. The control law achieves almost-global tracking for any smooth reference trajectory specified for the unactuated joint. Further, we generalize the proposed scheme to an n-link system with as many (or more) degrees of actuation than unactuation, and show that the result holds.


title: Speeding up Context-based Sentence Representation Learning with Non-autoregressive Convolutional Decoding

abstract: Context plays an important role in human language understanding, thus it may also be useful for machines learning vector representations of language. In this paper, we explore an asymmetric encoder-decoder structure for unsupervised context-based sentence representation learning. We carefully designed experiments to show that neither an autoregressive decoder nor an RNN decoder is required. After that, we designed a model which still keeps an RNN as the encoder, while using a non-autoregressive convolutional decoder. We further combine a suite of effective designs to significantly improve model efficiency while also achieving better performance. Our model is trained on two different large unlabelled corpora, and in both cases the transferability is evaluated on a set of downstream NLP tasks. We empirically show that our model is simple and fast while producing rich sentence representations that excel in downstream tasks.


title: Partitioning Relational Matrices of Similarities or Dissimilarities using the Value of Information

abstract: In this paper, we provide an approach to clustering relational matrices whose entries correspond to either similarities or dissimilarities between objects. Our approach is based on the value of information, a parameterized, information-theoretic criterion that measures the change in costs associated with changes in information. Optimizing the value of information yields a deterministic annealing style of clustering with many benefits. For instance, investigators avoid needing to a priori specify the number of clusters, as the partitions naturally undergo phase changes, during the annealing process, whereby the number of clusters changes in a data-driven fashion. The global-best partition can also often be identified.


title: Information-Based Optimal Subdata Selection for Big Data Linear Regression

abstract: Extraordinary amounts of data are being produced in many branches of science. Proven statistical methods are no longer applicable with extraordinary large data sets due to computational limitations. A critical step in big data analysis is data reduction. Existing investigations in the context of linear regression focus on subsampling-based methods. However, not only is this approach prone to sampling errors, it also leads to a covariance matrix of the estimators that is typically bounded from below by a term that is of the order of the inverse of the subdata size. We propose a novel approach, termed information-based optimal subdata selection (IBOSS). Compared to leading existing subdata methods, the IBOSS approach has the following advantages: (i) it is significantly faster; (ii) it is suitable for distributed parallel computing; (iii) the variances of the slope parameter estimators converge to 0 as the full data size increases even if the subdata size is fixed, i.e., the convergence rate depends on the full data size; (iv) data analysis for IBOSS subdata is straightforward and the sampling distribution of an IBOSS estimator is easy to assess. Theoretical results and extensive simulations demonstrate that the IBOSS approach is superior to subsampling-based methods, sometimes by orders of magnitude. The advantages of the new approach are also illustrated through analysis of real data.


title: A Geometric Perspective on the Power of Principal Component Association Tests in Multiple Phenotype Studies

abstract: Joint analysis of multiple phenotypes can increase statistical power in genetic association studies. Principal component analysis, as a popular dimension reduction method, especially when the number of phenotypes is high-dimensional, has been proposed to analyze multiple correlated phenotypes. It has been empirically observed that the first PC, which summarizes the largest amount of variance, can be less powerful than higher order PCs and other commonly used methods in detecting genetic association signals. In this paper, we investigate the properties of PCA-based multiple phenotype analysis from a geometric perspective by introducing a novel concept called principal angle. A particular PC is powerful if its principal angle is $0^o$ and is powerless if its principal angle is $90^o$. Without prior knowledge about the true principal angle, each PC can be powerless. We propose linear, non-linear and data-adaptive omnibus tests by combining PCs. We show that the omnibus PC test is robust and powerful in a wide range of scenarios. We study the properties of the proposed methods using power analysis and eigen-analysis. The subtle differences and close connections between these combined PC methods are illustrated graphically in terms of their rejection boundaries. Our proposed tests have convex acceptance regions and hence are admissible. The $p$-values for the proposed tests can be efficiently calculated analytically and the proposed tests have been implemented in a publicly available R package {\it MPAT}. We conduct simulation studies in both low and high dimensional settings with various signal vectors and correlation structures. We apply the proposed tests to the joint analysis of metabolic syndrome related phenotypes with data sets collected from four international consortia to demonstrate the effectiveness of the proposed combined PC testing procedures.


title: Single wavelength 480 Gb/s direct detection over 80km SSMF enabled by Stokes Vector Kramers Kronig transceiver

abstract: We propose 4D modulation with directed detection employing a novel Stokes-Vector Kramers-Kronig transceiver. It shows that employing Stokes vector receiver, transmitted digital carrier and Kramers-Kronig detection offers an effective way to de-rotate polarization multiplexed complex double side band signal without using a local oscillator at receiver. The impact of system parameters and configurations including carrier-to-signal-power ratio, guard band of the digital carrier, oversampling ratio and real MIMO is experimentally investigated. Finally, a record 480 Gb/s data rate over 80 km SSMF is achieved in a 60 Gbaud PDM-16QAM single carrier experiment with a BER below the threshold of 2.0x10-2


title: Capturing the Future by Replaying the Past

abstract: Delimited continuations are the mother of all monads! So goes the slogan inspired by Filinski's 1994 paper, which showed that delimited continuations can implement any monadic effect, letting the programmer use an effect as easily as if it was built into the language. It's a shame that not many languages have delimited continuations. Luckily, exceptions and state are also the mother of all monads! In this Pearl, we show how to implement delimited continuations in terms of exceptions and state, a construction we call $\textit{thermometer continuations}$. While traditional implementations of delimited continuations require some way of "capturing" an intermediate state of the computation, the insight of thermometer continuations is to reach this intermediate state by replaying the entire computation from the start, guiding it using a recording it so that the same thing happens until the captured point. Along the way, we explain delimited continuations and monadic reflection, show how the Filinski construction lets thermometer continuations express any monadic effect, share an elegant special-case for nondeterminism, and discuss why our construction is not prevented by theoretical results that exceptions and state cannot macro-express continuations.


title: Dual Skipping Networks

abstract: Inspired by the recent neuroscience studies on the left-right asymmetry of the human brain in processing low and high spatial frequency information, this paper introduces a dual skipping network which carries out coarse-to-fine object categorization. Such a network has two branches to simultaneously deal with both coarse and fine-grained classification tasks. Specifically, we propose a layer-skipping mechanism that learns a gating network to predict which layers to skip in the testing stage. This layer-skipping mechanism endows the network with good flexibility and capability in practice. Evaluations are conducted on several widely used coarse-to-fine object categorization benchmarks, and promising results are achieved by our proposed network model.


title: A Range-Doppler-Angle Estimation Method for Passive Bistatic Radar

abstract: In this paper, an effective target detection and localization method is proposed for a passive bistatic radar (PBR) system. The PBR system consists of a commercial FM radio station, which is a non-cooperative illuminator of opportunity (IO), referred to as the transmitter antenna and multiple surveillance antennas that form an antenna array, e.g., uniform linear array (ULA). Unlike other literatures where the reference signal is received by a directional antenna, here, the reference signal (direct path) is estimated by beamforming method. Then a modified extensive cancellation algorithm (MECA) based on (least squares) LS method is proposed to solve the disturbance cancellation. After cancelling the disturbance, the matched filter (MF) and LS methods are used for range-Doppler estimation of targets, and then the angles of targets are estimated based on beamforming method. The proposed method is suitable for an antenna array. Simulation results are presented to illustrate the superiority of the proposed MECA disturbance cancellation method and parameter estimation method.


title: Minimax Rates and Efficient Algorithms for Noisy Sorting

abstract: There has been a recent surge of interest in studying permutation-based models for ranking from pairwise comparison data. Despite being structurally richer and more robust than parametric ranking models, permutation-based models are less well understood statistically and generally lack efficient learning algorithms. In this work, we study a prototype of permutation-based ranking models, namely, the noisy sorting model. We establish the optimal rates of learning the model under two sampling procedures. Furthermore, we provide a fast algorithm to achieve near-optimal rates if the observations are sampled independently. Along the way, we discover properties of the symmetric group which are of theoretical interest.


title: Blocking Probability and Spatial Throughput Characterization for Cellular-Enabled UAV Network with Directional Antenna

abstract: The past few years have witnessed a tremendous increase on the use of unmanned aerial vehicles (UAVs) in civilian applications, which call for high-performance communication between UAVs and their ground clients, especially when they are densely deployed. To achieve this goal, cellular base stations (BSs) can be leveraged to provide a new and promising solution to support massive UAV communications simultaneously in a cost-effective way. However, different from terrestrial communication channels, UAV-to-BS channels are usually dominated by the light-of-sight (LoS) link, which aggravates the co-channel interference and renders the spatial frequency reuse in existing cellular networks ineffective. In this paper, we consider the use of a directional antenna at each UAV to confine the interference to/from other UAV users within a limited region and hence improve the spatial reuse of the spectrum. Under this model, a UAV user may be temporarily blocked from communication if it cannot find any BS in its antenna main-lobe, or it finds that all BSs under its main-lobe are simultaneously covered by those of some other UAVs and hence suffer from strong co-channel interference. Assuming independent homogeneous Poisson point processes (HPPPs) for the UAVs' and ground BSs' locations respectively, we first analytically derive a closed-form upper bound for the UAV blocking probability and then characterize the achievable average spatial throughput of the cellular-enabled UAV communication network, in terms of various key parameters including the BS/UAV densities as well as the UAV's flying altitude and antenna beamwidth. Simulation results verify that the derived bound is practically tight, and further show that adaptively adjusting the UAV altitude and/or beamwidth with different BS/UAV densities can significantly reduce the UAV blocking probability and hence improve the network spatial throughput.


title: Static Black Hole and Vacuum Energy: Thin Shell and Incompressible Fluid

abstract: With the back reaction of the vacuum energy-momentum tensor consistently taken into account, we study static spherically symmetric black-hole-like solutions to the semi-classical Einstein equation. The vacuum energy is assumed to be given by that of 2-dimensional massless scalar fields, as a widely used model in the literature for black holes. The solutions have no horizon. Instead, there is a local minimum in the radius. We consider thin shells as well as incompressible fluid as the matter content of the black-hole-like geometry. The geometry has several interesting features due to the back reaction of vacuum energy. In particular, Buchdahl's inequality can be violated without divergence in pressure, even if the surface is below the Schwarzschild radius. At the same time, the surface of the star can not be far below the Schwarzschild radius for a density not much higher than the Planck scale, and the proper distance from its surface to the origin can be very short even for very large Schwarzschild radius. The results also imply that, contrary to the folklore, in principle the Boulware vacuum can be physical for black holes.


title: Introduction and reconciliation of the ROS and aging paradoxes

abstract: This paper suggests that aging is influenced synthetically by pro-aging factors such as ROS and anti-aging factors such as protective responses. The anti-aging effect may be side effects of retrograde responses motivated against adverse circumstances. ROS may be more closely correlated with metabolism rather than aging.


title: Summability methods and Fourier analysis on $\mathbb{R}^{\times}$

abstract: We introduce a certain class of summability methods on $L^{\infty}([1, \infty))$ which is defined by convolution in the group algebra $L^1(\mathbb{R}^{\times})$. This class contains an integral version of Ces\`{a}ro summability method and we particularly give a necessary and sufficient condition for a summability method in the class to equivalent to this one. Also, we consider another class of summability methods concerned with convolution in the group algebra $L^1(\mathbb{R})$ and give similar results.


title: Label Embedding Network: Learning Label Representation for Soft Training of Deep Networks

abstract: We propose a method, called Label Embedding Network, which can learn label representation (label embedding) during the training process of deep networks. With the proposed method, the label embedding is adaptively and automatically learned through back propagation. The original one-hot represented loss function is converted into a new loss function with soft distributions, such that the originally unrelated labels have continuous interactions with each other during the training process. As a result, the trained model can achieve substantially higher accuracy and with faster convergence speed. Experimental results based on competitive tasks demonstrate the effectiveness of the proposed method, and the learned label embedding is reasonable and interpretable. The proposed method achieves comparable or even better results than the state-of-the-art systems. The source code is available at \url{https://github.com/lancopku/LabelEmb}.


title: M-RWTL: Learning Signal-Matched Rational Wavelet Transform in Lifting Framework

abstract: Transform learning is being extensively applied in several applications because of its ability to adapt to a class of signals of interest. Often, a transform is learned using a large amount of training data, while only limited data may be available in many applications. Motivated with this, we propose wavelet transform learning in the lifting framework for a given signal. Significant contributions of this work are: 1) the existing theory of lifting framework of the dyadic wavelet is extended to more generic rational wavelet design, where dyadic is a special case and 2) the proposed work allows to learn rational wavelet transform from a given signal and does not require large training data. Since it is a signal-matched design, the proposed methodology is called Signal-Matched Rational Wavelet Transform Learning in the Lifting Framework (M-RWTL). The proposed M-RWTL method inherits all the advantages of lifting, i.e., the learned rational wavelet transform is always invertible, method is modular, and the corresponding M-RWTL system can also incorporate nonlinear filters, if required. This may enhance the use of RWT in applications which is so far restricted. M-RWTL is observed to perform better compared to standard wavelet transforms in the applications of compressed sensing based signal reconstruction.


title: Local approximation of a metapopulation's equilibrium

abstract: We consider the approximation of the equilibrium of a metapopulation model, in which a finite number of patches are randomly distributed over a bounded subset $\Omega$ of Euclidean space. The approximation is good when a large number of patches contribute to the colonization pressure on any given unoccupied patch, and when the quality of the patches varies little over the length scale determined by the colonization radius. If this is the case, the equilibrium probability of a patch at $z$ being occupied is shown to be close to $q_1(z)$, the equilibrium occupation probability in Levins's model, at any point $z \in \Omega$ not too close to the boundary, if the local colonization pressure and extinction rates appropriate to $z$ are assumed. The approximation is justified by giving explicit upper and lower bounds for the occupation probabilities, expressed in terms of the model parameters. Since the patches are distributed randomly, the occupation probabilities are also random, and we complement our bounds with explicit bounds on the probability that they are satisfied at all patches simultaneously.


title: Role of intervening Mg II absorbers on the rotation measure and fractional polarisation of the background quasars

abstract: We probed the magnetic fields in high-redshift galaxies using excess extragalactic contribution to residual rotation measure (RRM) for quasar sightlines with intervening Mg II absorbers. Based on a large sample of 1132 quasars, we have computed RRM distributions broadening using median absolute deviation from mean ($\sigma^{md}_{rrm}$), and found it to be 17.1$\pm0.7$ rad m$^{-2}$ for 352 sightlines having Mg II intervening absorbers in comparison to its value of $15.1\pm0.6$ rad m$^{-2}$ for 780 sightlines without such absorbers, resulting in an excess broadening ($\sigma_{rrm}^{ex}$) of $8.0\pm1.9$ rad m$^{-2}$ among these two subsamples. This value of $\sigma_{rrm}^{ex}$, has allowed us to constrain the average strength of magnetic field (rest frame) in high redshift galaxies responsible for these Mg II absorbers, to be $\sim 1.3\pm0.3 \mu G$ at a median redshift of 0.92. This estimate of magnetic field is consistent with the reported estimate in earlier studies based on radio-infrared correlation and energy equipartition for galaxies in local universe. A similar analysis on subsample split based on the radio spectral index, $\alpha$, (with $F_{\nu}\propto \nu^{\alpha}$) for flat ($\alpha$ $\geq -0.3$; 315 sources) and steep ($\alpha$ $\leq -0.7$; 476 sources) spectrum sources shows a significant $\sigma_{rrm}^{ex}$ (at 3.5$\sigma$ level) for the former and absent in latter. An anti-correlation found between the $\sigma^{md}_{rrm}$ and percentage polarisation ($p$) with similar Pearson correlation of $-0.62$ and $-0.87$ for subsample with and without Mg II, respectively, suggests main contribution for decrements in the $p$ value to be intrinsic to the local environment of quasars.


title: Geometric Decomposition-Based Formulation for Time Derivatives of Instantaneous Impact Point

abstract: A new analytic formulation to express the time derivatives of the instantaneous impact point (IIP) of a rocket is proposed. The geometric relationship on a plane tangential to the IIP is utilized to decompose the inertial IIP rate vector into the downrange and crossrange components, and a systematic procedure to determine the component values is presented. The new formulation shows significant advantages over the existing formulation such that the procedure and final expressions for the IIP derivatives are easy to understand and more compact. The validity of the proposed formulation was demonstrated through numerical simulation.


title: A Study of All-Convolutional Encoders for Connectionist Temporal Classification

abstract: Connectionist temporal classification (CTC) is a popular sequence prediction approach for automatic speech recognition that is typically used with models based on recurrent neural networks (RNNs). We explore whether deep convolutional neural networks (CNNs) can be used effectively instead of RNNs as the "encoder" in CTC. CNNs lack an explicit representation of the entire sequence, but have the advantage that they are much faster to train. We present an exploration of CNNs as encoders for CTC models, in the context of character-based (lexicon-free) automatic speech recognition. In particular, we explore a range of one-dimensional convolutional layers, which are particularly efficient. We compare the performance of our CNN-based models against typical RNNbased models in terms of training time, decoding time, model size and word error rate (WER) on the Switchboard Eval2000 corpus. We find that our CNN-based models are close in performance to LSTMs, while not matching them, and are much faster to train and decode.


title: Dynamics of a dry-rebounding drop: observations, simulations, and modeling

abstract: Dynamics of a dry-rebounding drop was studied experimentally, numerically, and theoretically. Experimental results were reproduced by our computational fluid dynamics simulations, from which time series of kinetic energy, potential energy, and surface energy were obtained. The time series of these energies quantitatively clarified the energy conversion and loss during the dry-rebound. These results were interpreted by using an imaginary spring model and a spherical harmonic analysis. The spring model explained the vertical deformation of the drop, however, could not completely explain the energy loss, the timings of the energy loss did not match. From a viewpoint of the spherical harmonic deformation of a drop, the deformation of the drop after the impact was found to be a combination of two vibrational motions. One of the two vibrational motions is an inertial motion derived from the free-fall and the another is a pressure-induced motion derived from a pressure surge due to the sudden stop of the bottom part of the drop at the impact. The existence of the pressure surge at the impact was confirmed in the simulated results. The pressure-induced motion resists the inertial motion and consequently dumps the kinetic energy of the drop.


title: Total-Text: A Comprehensive Dataset for Scene Text Detection and Recognition

abstract: Text in curve orientation, despite being one of the common text orientations in real world environment, has close to zero existence in well received scene text datasets such as ICDAR2013 and MSRA-TD500. The main motivation of Total-Text is to fill this gap and facilitate a new research direction for the scene text community. On top of the conventional horizontal and multi-oriented texts, it features curved-oriented text. Total-Text is highly diversified in orientations, more than half of its images have a combination of more than two orientations. Recently, a new breed of solutions that casted text detection as a segmentation problem has demonstrated their effectiveness against multi-oriented text. In order to evaluate its robustness against curved text, we fine-tuned DeconvNet and benchmark it on Total-Text. Total-Text with its annotation is available at https://github.com/cs-chan/Total-Text-Dataset


title: The pion vector form factor from Lattice QCD at the physical point

abstract: We present an investigation of the electromagnetic pion form factor, $F_\pi(Q^2)$, at small values of the four-momentum transfer $Q^2$ ($\lesssim 0.25$ GeV$^2$), based on the gauge configurations generated by European Twisted Mass Collaboration with $N_f = 2$ twisted-mass quarks at maximal twist including a clover term. Momentum is injected using non-periodic boundary conditions and the calculations are carried out at a fixed lattice spacing ($a \simeq 0.09$ fm) and with pion masses equal to its physical value, 240 MeV and 340 MeV. Our data are successfully analyzed using Chiral Perturbation Theory at next-to-leading order in the light-quark mass. For each pion mass two different lattice volumes are used to take care of finite size effects. Our final result for the squared charge radius is $\langle r^2 \rangle_\pi = 0.443~(29)$ fm$^2$, where the error includes several sources of systematic errors except the uncertainty related to discretization effects. The corresponding value of the SU(2) chiral low-energy constant $\overline{\ell}_6$ is equal to $\overline{\ell}_6 = 16.2 ~ (1.0)$.


title: Termination in Convex Sets of Distributions

abstract: Convex algebras, also called (semi)convex sets, are at the heart of modelling probabilistic systems including probabilistic automata. Abstractly, they are the Eilenberg-Moore algebras of the finitely supported distribution monad. Concretely, they have been studied for decades within algebra and convex geometry. In this paper we study the problem of extending a convex algebra by a single point. Such extensions enable the modelling of termination in probabilistic systems. We provide a full description of all possible extensions for a particular class of convex algebras: For a fixed convex subset $D$ of a vector space satisfying additional technical condition, we consider the algebra of convex subsets of $D$. This class contains the convex algebras of convex subsets of distributions, modelling (nondeterministic) probabilistic automata. We also provide a full description of all possible extensions for the class of free convex algebras, modelling fully probabilistic systems. Finally, we show that there is a unique functorial extension, the so-called black-hole extension.


title: Trainable back-propagated functional transfer matrices

abstract: Connections between nodes of fully connected neural networks are usually represented by weight matrices. In this article, functional transfer matrices are introduced as alternatives to the weight matrices: Instead of using real weights, a functional transfer matrix uses real functions with trainable parameters to represent connections between nodes. Multiple functional transfer matrices are then stacked together with bias vectors and activations to form deep functional transfer neural networks. These neural networks can be trained within the framework of back-propagation, based on a revision of the delta rules and the error transmission rule for functional connections. In experiments, it is demonstrated that the revised rules can be used to train a range of functional connections: 20 different functions are applied to neural networks with up to 10 hidden layers, and most of them gain high test accuracies on the MNIST database. It is also demonstrated that a functional transfer matrix with a memory function can roughly memorise a non-cyclical sequence of 400 digits.


title: Efficient Localized Inference for Large Graphical Models

abstract: We propose a new localized inference algorithm for answering marginalization queries in large graphical models with the correlation decay property. Given a query variable and a large graphical model, we define a much smaller model in a local region around the query variable in the target model so that the marginal distribution of the query variable can be accurately approximated. We introduce two approximation error bounds based on the Dobrushin's comparison theorem and apply our bounds to derive a greedy expansion algorithm that efficiently guides the selection of neighbor nodes for localized inference. We verify our theoretical bounds on various datasets and demonstrate that our localized inference algorithm can provide fast and accurate approximation for large graphical models.


title: Overcoming the SQL in gravitational wave detectors using spin systems with negative effective mass

abstract: Quantum back action (QBA) of a measurement limits the precision of observation of the motion of a free mass. This profound effect dabbed the "Heisenberg microscope" in the early days of quantum mechanics, leads to the standard quantum limit (SQL) stemming from the balance between the measurement sensitivity and the QBA. Here we consider the measurement of motion of a free mass performed in a quantum reference frame with an effective negative mass which is not limited by QBA. As a result, the disturbance on the motion of a free mass can be measured beyond SQL. QBA-limited detection of motion for a free mass is extremely challenging, but there are devices where this effect is expected to play an essential role, namely, gravitational wave detectors (GWD) such as LIGO and VIRGO. Recent reports on observation of gravitational waves have opened new horizons in cosmology and astrophysics. Here we present a general idea and a detailed numerical analysis for QBA-evading measurement of the gravitational wave effect on the GWD mirrors which can be considered free masses under relevant conditions. The measurement is performed by two entangled beams of light probing the GWD and an auxiliary atomic spin ensemble, respectively. The latter plays a role of a free negative mass. We show that under realistic conditions the sensitivity of the GWD can be significantly increased over the entire frequency band of interest.


title: An adaptive variational procedure for the conservative and positivity preserving Allen-Cahn phase-field model

abstract: We present an adaptive variational procedure for unstructured meshes to capture fluid-fluid interfaces in two-phase flows. The two phases are modeled by the phase-field finite element formulation, which involves the conservative Allen-Cahn equation coupled with the incompressible Navier-Stokes equations. The positivity preserving variational formulation is designed to maintain the bounded and stable solution of the Allen-Cahn equation. For the adaptivity procedure, we consider the residual-based error estimates for the underlying differential equations of the two-phase system. In particular, the adaptive refinement/coarsening is carried out by the newest vertex bisection algorithm by evaluating the residual error indicators based on the error estimates of the Allen-Cahn equation. The coarsening algorithm avoids the storage of the tree data structures for the hierarchical mesh, thus providing the ease of numerical implementation. The proposed nonlinear adaptive partitioned procedure aims at reducing the amount of coarsening while maintaining the convergence properties of the underlying nonlinear coupled differential equations. We investigate the adaptive phase-field finite element scheme through the spinodal decomposition in a complex curved geometry and the volume-conserved interface motion driven by the mean curvature flow for two circles in a square domain. We then assess the accuracy and efficiency of the proposed procedure by modeling the free-surface motion in a sloshing tank. The mesh adaptivity remarkably reduces the degrees of freedom and the computational cost by nearly half for similar accuracy and improves the mass conservation. Finally, we apply the adaptive numerical framework to solve the application of a dam-breaking problem with topological changes.


title: Doppelgangers: the Ur-Operation and Posets of Bounded Height

abstract: In the early 1970's, Richard Stanley and Kenneth Johnson introduced and laid the groundwork for studying the order polynomial of partially ordered sets (posets). Decades later, Hamaker, Patrias, Pechenik, and Williams introduced the term "doppelgangers": equivalence classes of posets under the equivalence relation given by equality of the order polynomial. We provide necessary and sufficient conditions on doppelgangers through application of both old and novel tools, including new recurrences and the Ur-operation: a new generalized poset operation. In addition, we prove that the doppelgangers of posets P of bounded height $|P|-k$ may be classified up to systems of $k$ diophantine equations in $2^{O(k^2)}$ time, and similarly that the order polynomial of such posets may be computed in $O(|P|)$ time.


title: Symmetric Implication Zroupoids and Weak Associative Laws

abstract: An algebra $\mathbf A = \langle A, \to, 0 \rangle$, where $\to$ is binary and $0$ is a constant, is called an implication zroupoid ($\mathcal I$-zroupoid, for short) if $\mathbf A$ satisfies the identities: $(x \to y) \to z \approx ((z' \to x) \to (y \to z)')'$ and $0'' \approx 0$, where $x' : = x \to 0$. An implication zroupoid is symmetric if it satisfies $x'' \approx x$ and $(x \to y')' \approx (y \to x')'$. The variety of symmetric $\mathcal I$-zroupoids is denoted by $\mathcal S$. We began a systematic analysis of weak associative laws of length $\leq 4$ in [CS16e], by examining the identities of Bol-Moufang type in the context of the variety $\mathcal S$. In this paper we complete the analysis by investigating the rest of the weak associative laws of length $\leq 4$ relative to $\mathcal S$. We show that, of the 155 subvarieties of $\mathcal S$ defined by the weak associative laws of size $\leq 4$, there are exactly $6$ distinct ones. We also give an explicit description of the poset of the (distinct) subvarieties of $\mathcal S$ defined by weak associative laws of length $\leq 4$.


title: On the geometry of the theory space in the ERG formalism

abstract: We consider the theory space as a manifold whose coordinates are given by the couplings appearing in the Wilson action. We discuss how to introduce connections on this theory space. A particularly intriguing connection can be defined directly from the solution of the exact renormalization group (ERG) equation. We advocate a geometric viewpoint that lets us define straightforwardly physically relevant quantities invariant under the changes of a renormalization scheme.