title: Improving Pre-Trained Multilingual Models with Vocabulary Expansion

abstract: Recently, pre-trained language models have achieved remarkable success in a broad range of natural language processing tasks. However, in multilingual setting, it is extremely resource-consuming to pre-train a deep language model over large-scale corpora for each language. Instead of exhaustively pre-training monolingual language models independently, an alternative solution is to pre-train a powerful multilingual deep language model over large-scale corpora in hundreds of languages. However, the vocabulary size for each language in such a model is relatively small, especially for low-resource languages. This limitation inevitably hinders the performance of these multilingual models on tasks such as sequence labeling, wherein in-depth token-level or sentence-level understanding is essential. In this paper, inspired by previous methods designed for monolingual settings, we investigate two approaches (i.e., joint mapping and mixture mapping) based on a pre-trained multilingual model BERT for addressing the out-of-vocabulary (OOV) problem on a variety of tasks, including part-of-speech tagging, named entity recognition, machine translation quality estimation, and machine reading comprehension. Experimental results show that using mixture mapping is more promising. To the best of our knowledge, this is the first work that attempts to address and discuss the OOV issue in multilingual settings.


title: Total Least Squares Regression in Input Sparsity Time

abstract: In the total least squares problem, one is given an $m \times n$ matrix $A$, and an $m \times d$ matrix $B$, and one seeks to "correct" both $A$ and $B$, obtaining matrices $\hat{A}$ and $\hat{B}$, so that there exists an $X$ satisfying the equation $\hat{A}X = \hat{B}$. Typically the problem is overconstrained, meaning that $m \gg \max(n,d)$. The cost of the solution $\hat{A}, \hat{B}$ is given by $\|A-\hat{A}\|_F^2 + \|B - \hat{B}\|_F^2$. We give an algorithm for finding a solution $X$ to the linear system $\hat{A}X=\hat{B}$ for which the cost $\|A-\hat{A}\|_F^2 + \|B-\hat{B}\|_F^2$ is at most a multiplicative $(1+\epsilon)$ factor times the optimal cost, up to an additive error $\eta$ that may be an arbitrarily small function of $n$. Importantly, our running time is $\tilde{O}( \mathrm{nnz}(A) + \mathrm{nnz}(B) ) + \mathrm{poly}(n/\epsilon) \cdot d$, where for a matrix $C$, $\mathrm{nnz}(C)$ denotes its number of non-zero entries. Importantly, our running time does not directly depend on the large parameter $m$. As total least squares regression is known to be solvable via low rank approximation, a natural approach is to invoke fast algorithms for approximate low rank approximation, obtaining matrices $\hat{A}$ and $\hat{B}$ from this low rank approximation, and then solving for $X$ so that $\hat{A}X = \hat{B}$. However, existing algorithms do not apply since in total least squares the rank of the low rank approximation needs to be $n$, and so the running time of known methods would be at least $mn^2$. In contrast, we are able to achieve a much faster running time for finding $X$ by never explicitly forming the equation $\hat{A} X = \hat{B}$, but instead solving for an $X$ which is a solution to an implicit such equation. Finally, we generalize our algorithm to the total least squares problem with regularization.


title: Optimal Precoding for Multiuser MIMO Systems With Phase Quantization and PSK Modulation via Branch-and-Bound

abstract: MIMO systems are considered as most promising for wireless communications. However, with an increasing number of radio front ends the corresponding energy consumption and costs become an issue, which can be relieved by the utilization of low-resolution quantizers. In this study we propose an optimal precoding algorithm constrained to constant envelope signals and phase quantization that maximizes the minimum distance to the decision threshold at the receivers using a branch-and-bound strategy. The proposed algorithm is superior to the existing methods in terms of bit error rate. Numerical results show that the proposed approach has significantly lower complexity than exhaustive search.


title: The Infall Motion in the Low-Mass Protostellar Binary NGC 1333 IRAS 4A1/4A2

abstract: We report ALMA observations of NGC 1333 IRAS 4A, a young low-mass protostellar binary, referred as 4A1 and 4A2. With multiple H$_2$CO transitions and HNC (4$-$3) observed at a resolution of 0.25" ($\sim$70 au), we investigate the gas kinematics of 4A1 and 4A2. Our results show that, on the large angular scale ($\sim$10"), 4A1 and 4A2 each display a well-collimated outflow along the N-S direction, and an S-shaped morphology is discerned in the outflow powered by 4A2. On the small scale ($\sim$0.3"), 4A1 and 4A2 exhibit distinct spectral features toward the continuum centroid, with 4A1 showing simple symmetric profiles predominantly in absorption and 4A2 demonstrating rather complicated profiles in emission as well as in absorption. Based on radiative transfer modeling exercises, we find that the physical parameters inferred from earlier low-resolution observations cannot be directly extrapolated down to the 4A1 inner region. Possible reasons for the discrepancies between the observed and modelled profiles are discussed. We constrain the mass infall rate in 4A1 to be at most around 3$\times$10$^{-5}$ M$_{\odot}$ year$^{-1}$ at the layer of 75 au. For the kinematics of the 4A2 inner envelope, the absorbing dips in the H$_2$CO spectra are skewed toward the redshifted side and likely signatures of inward motion. These absorbing dips are relatively narrow. This is, like the case for 4A1, significantly slower than the anticipated inflow speed. We estimate a mass infall rate of 3.1$-$6.2 $\times$ 10$^{-5}$ M$_\odot$ year$^{-1}$ at the layer of 100 au in 4A2.


title: A Momentum-Based Foot Placement Strategy for Stable Postural Control of Robotic Spring-Mass Running with Point Feet

abstract: A long-standing argument in model-based control of locomotion is about the level of complexity that a model should have to define a behavior such as running. Even though goldilocks model based on biomechanical evidence is often sought, it is unclear what complexity level qualifies to be such a model. This dilemma deepens further for bipedal robotic running with point feet, since these robots are underactuated, while tracking center-of-mass (COM) trajectories defined by the spring-loaded inverted pendulum (SLIP) model of running allocates all control inputs, leaving angular coordinates of the robot's trunk uncontrolled. Existing work in the literature approach this problem either by trading off COM trajectories against upright trunk posture during stance or by adopting more detailed models that include effects of trunk angular dynamics. In this paper, we present a new approach based on modifying foot placement targets of the SLIP model. Theoretical analysis and numerical results show that the proposed approach outperforms these traditional strategies.


title: Thermally driven two-magnet nano-oscillator with large spin-charge conversion

abstract: Next-generation spintronic applications require material properties that can be hardly met by one material candidate. Here we demonstrate that by combining insulating and metallic magnets, enhanced spin-charge conversion and energy-efficient thermal spin currents can be realized. We develop a nanowire device consisting of an yttrium iron garnet and permalloy bi-layer. An interfacial temperature gradient drives the nanowire magnetization into auto-oscillations at gigahertz frequencies. Interfacial spin coupling and magnetoresistance of the permalloy layer translate spin dynamics into sizable microwave signals. The results show prospect for energy-efficient spintronic devices and present an experimental realization of magnon condensation in a heterogeneous magnetic system.


title: Interpreting Undesirable Pixels for Image Classification on Black-Box Models

abstract: In an effort to interpret black-box models, researches for developing explanation methods have proceeded in recent years. Most studies have tried to identify input pixels that are crucial to the prediction of a classifier. While this approach is meaningful to analyse the characteristic of blackbox models, it is also important to investigate pixels that interfere with the prediction. To tackle this issue, in this paper, we propose an explanation method that visualizes undesirable regions to classify an image as a target class. To be specific, we divide the concept of undesirable regions into two terms: (1) factors for a target class, which hinder that black-box models identify intrinsic characteristics of a target class and (2) factors for non-target classes that are important regions for an image to be classified as other classes. We visualize such undesirable regions on heatmaps to qualitatively validate the proposed method. Furthermore, we present an evaluation metric to provide quantitative results on ImageNet.


title: Co-axial Helicity Injection on the STOR-M Tokamak

abstract: Injection of relatively high density spheromaks with significant helicity-content into a tokamak has been proposed as a means for fueling and current drive. The CHI (Co-axial Helicity Injection) device was devised to inject current to the STOR-M tokamak. Various circuit modifications were made to the CHI controls, enabling testing of various injection configurations. The charge/discharge circuits for CT formation/acceleration and stuffing field were modified, and the power supplies and power converters were replaced. Various modifications were implemented to solve the original slow bank triggering problems. The CHI device was mounted on STOR-M for radial and vertical CT injection at various times. Spheromak injection into STOR-M usually resulted in disruption of the tokamak discharge. After modifying the CHI device to operate at increased power, it looked like tokamak current was increased by a few kiloamps just prior to disruption, but careful testing proved that the signals indicating a current drive were actually spurious, caused by inductive pickup. The CHI device was attached to a portable vacuum chamber that was constructed from spare parts, to characterise the CTs produced. Magnetic probes were constructed to measure poloidal and toroidal field near the CT edge. Langmuir probes were made and returned reasonable estimates for edge CT density and temperature.


title: Combined Energy and Comfort Optimization of Air Conditioning System in Connected and Automated Vehicles

abstract: In this paper, we propose a combined energy and comfort optimization (CECO) strategy for the air conditioning (A/C) system of the connected and automated vehicles (CAVs). By leveraging the weather and traffic predictions enabled by the emerging CAV technologies, the proposed strategy is able to minimize the A/C system energy consumption while maintaining the occupant thermal comfort (OTC) within the comfort constraints, where the comfort is quantified by a modified predictive mean vote (PMV) model adapted for an automotive application. A general CECO problem is formulated and addressed using model predictive control (MPC) and weather/traffic previews. Depending on the ways of exploiting the preview information and enforcing the OTC constraint, different MPCs are developed based on solving different variations of the general CECO problem. The CECO-based MPCs are then tested in simulation using an automotive A/C system simulation model (CoolSim) as the virtual testbed. The simulation results show that, over SC03 driving cycle, the proposed CECO-based MPCs outperform the baseline cabin temperature tracking controller, reducing the A/C system energy consumption by up to 7.6%, while achieving better OTC according to the PMV-based metrics. This energy saving in A/C system translates to 3.1% vehicle fuel economy improvement. The trade-off between energy efficiency and OTC for different control scenarios is also highlighted.


title: Electrodiffusion with calcium-activated potassium channels in dendritic spine

abstract: We investigate calcium signaling feedback through calcium-activated potassium channels of a dendritic spine by applying the immersed boundary method with electrodiffusion. We simulate the stochastic gating of such ion channels and the resulting spatial distribution of concentration, current, and membrane voltage within the dendritic spine. In this simulation, the permeability to ionic flow across the membrane is regulated by the amplitude of chemical potential barriers. With spatially localized ion channels, chemical potential barriers are locally and stochastically regulated. This regulation represents the ion channel gating with multiple subunits, the open and closed states governed by a continuous-time Markov process. The model simulation recapitulates an inhibitory action on voltage-sensitive calcium channels by the calcium-activated potassium channels in a stochastic manner as a \emph{non-local} feedback loop. The model predicts amplified calcium influx with more closely placed channel complexes, proposing a potential mechanism of differential calcium handling by channel distributions. This work provides a foundation for future computer simulation studies of dendritic spine motility and structural plasticity.


title: Quasi-periodic dynamics and a Neimark-Sacker bifurcation in nonlinear random walks on complex networks

abstract: We study the dynamics of nonlinear random walks on complex networks. We investigate the role and effect of directed network topologies on long-term dynamics. While a period-doubling bifurcation to alternating patterns occurs at a critical bias parameter value, we find that some directed structures give rise to a different kind of bifurcation that gives rise to quasi-periodic dynamics. This does not occur for all directed network structure, but only when the network structure is sufficiently directed. We find that the onset of quasi-periodic dynamics is the result of a Neimark-Sacker bifurcation, where a pair of complex-conjugate eigenvalues of the system Jacobian passes through the unit circle, destabilizing the stationary distribution with high-dimensional rotations. We investigate the nature of these bifurcations, study the onset of quasi-periodic dynamics as network structure is tuned to be more directed, and present an analytically tractable case of a four-neighbor ring.


title: Stray light analysis of the Xinglong 2.16m telescope

abstract: An effort towards understanding of the stray light problems for the Xinglong 2.16-m telescope was presented to estimate the stray light performance of the telescope itself and provide a method for improving the stray light suppression. The stray light analysis for 2.16-m telescope model, which consists the onion shaped dome, telescope structure, equatorial mount and telescope optics, has been performed with two cases (1) point to 60{\deg} and (2) point to zenith, in both azimuth and elevation direction. The Point Source Normalized Irradiance Transmittance (PSNIT), which is generally used for assessing stray light and uncorrelated to entrance aperture, was calculated with a series of off-axis angles. It shows that the PSNIT values are less than 10-7 when off-axis angles are larger than 20{\deg}. The dominant contributors of stray light (primary and secondary mirror, telescope structure and dome) were identified to give advice for performance improvement. The analyses indicate that significant benefit can be realized with adding only 5 vanes inside the bottom portion of the secondary baffle. In the case of point to zenith, the PSNIT values will decrease about 40% at average.


title: Co-design for Security and Performance: LMI Tools

abstract: We present a convex optimization to reduce the impact of sensor falsification attacks in linear time invariant systems controlled by observer-based feedback. We accomplish this by finding optimal observer and controller gain matrices that minimize the size of the reachable set of attack-induced states. To avoid trivial solutions, we integrate a covariance-based $\|H\|_2$ closed-loop performance constraint, for which we develop a novel linearization for this typically nonlinear, non-convex problem. We demonstrate the effectiveness of this linear matrix inequality framework through a numerical case study.


title: Accurate atomic correlation and total energies for correlation consistent effective core potentials

abstract: Very recently, we introduced a set of correlation consistent effective core potentials (ccECPs) constructed within full many-body approaches. By employing significantly more accurate correlated approaches we were able to reach a new level of accuracy for the resulting effective core Hamiltonians. We also strived for simplicity of use and easy transferability into a variety of electronic structure methods in quantum chemistry and condensed matter physics. Here, as a reference for future use, we present exact or nearly-exact total energy calculations for these ccECPs. The calculations cover H-Kr elements and are based on the state-of-the-art configuration interaction (CI), coupled-cluster (CC), and quantum Monte Carlo (QMC) calculations with systematically eliminated/improved errors. In particular, we carry out full CI/CCSD(T)/CCSDT(Q) calculations with cc-pVnZ with up to n=6 basis sets and we estimate the complete basis set limits. Using combinations of these approaches, we achieved an accuracy of $\approx$ 1-10 mHa for K-Zn atoms and $\approx$ 0.1-0.3 mHa for all other elements $-$ within about 1% or better of the ccECP total correlation energies. We also estimate the corresponding kinetic energies within the feasible limit of full CI calculations. In order to provide data for QMC calculations, we include fixed-node diffusion Monte Carlo energies for each element that give quantitative insights into the fixed-node biases for single-reference trial wave functions. The results offer a clear benchmark for future high accuracy calculations in a broad variety of correlated wave function methods such as CI and CC as well is in stochastic approaches such as real space sampling QMC.


title: SoK: Blockchain Technology and Its Potential Use Cases

abstract: Bitcoin's success has led to significant interest in its underlying components, particularly Blockchain technology. Over 10 years after Bitcoin's initial release, the community still suffers from a lack of clarity regarding what properties defines Blockchain technology, its relationship to similar technologies, and which of its proposed use-cases are tenable and which are little more than hype. In this paper we answer four common questions regarding Blockchain technology: (1) what exactly is Blockchain technology, (2) what capabilities does it provide, and (3) what are good applications for Blockchain technology, and (4) how does it relate to other approache distributed technologies (e.g., distributed databases). We accomplish this goal by using grounded theory (a structured approach to gathering and analyzing qualitative data) to thoroughly analyze a large corpus of literature on Blockchain technology. This method enables us to answer the above questions while limiting researcher bias, separating thought leadership from peddled hype and identifying open research questions related to Blockchain technology. The audience for this paper is broad as it aims to help researchers in a variety of areas come to a better understanding of Blockchain technology and identify whether it may be of use in their own research.


title: Completely Positive Quantum Trajectories with Applications to Quantum State Smoothing

abstract: Here, we are concerned with comparing estimation schemes for the quantum state under continuous measurement (quantum trajectories), namely quantum state filtering and, as introduced by us [Phys. Rev. Lett. 115, 180407 (2015)], quantum state smoothing. Unfortunately, the cumulative errors in the most typical simulations of quantum trajectories with a total time of simulation $T$ can reach orders of $T \Delta t$. Moreover, these errors may correspond to deviations from valid quantum evolution as described by a completely positive map. Here we introduce a higher-order method that reduces the cumulative errors in the complete positivity of the evolution to of order $T\Delta t^2$, whether for linear (unnormalised) or nonlinear (normalised) quantum trajectories. Our method also guarantees that the discrepancy in the average evolution between different detection methods (different `unravellings', such as quantum jumps or quantum diffusion) is similarly small. This equivalence is essential for comparing quantum state filtering to quantum state smoothing, as the latter assumes that all irreversible evolution is unravelled, although the estimator only has direct knowledge of some records. In particular, here we compare, for the first time, the average difference between filtering and smoothing conditioned on an event of which the estimator lacks direct knowledge: a photon detection within a certain time window. We find that the smoothed state is actually {\em less pure}, both before and after the time of the jump. Similarly, the fidelity of the smoothed state with the `true' (maximal knowledge) state is also lower than that of the filtered state before the jump. However, after the jump, the fidelity of the smoothed state is higher.


title: ASSD: Attentive Single Shot Multibox Detector

abstract: This paper proposes a new deep neural network for object detection. The proposed network, termed ASSD, builds feature relations in the spatial space of the feature map. With the global relation information, ASSD learns to highlight useful regions on the feature maps while suppressing the irrelevant information, thereby providing reliable guidance for object detection. Compared to methods that rely on complicated CNN layers to refine the feature maps, ASSD is simple in design and is computationally efficient. Experimental results show that ASSD competes favorably with the state-of-the-arts, including SSD, DSSD, FSSD and RetinaNet.


title: Nonlinear Dynamic Models of Conflict via Multiplexed Interaction Networks

abstract: The risk of conflict is exasperated by a multitude of internal and external factors. Current multivariate analysis paints diverse causal risk profiles that vary with time. However, these profiles evolve and a universal model to understand that evolution remains absent. Most of the current conflict analysis is data-driven and conducted at the individual country or region level, often in isolation. Consistent consideration of multi-scale interactions and their non-linear dynamics is missing. Here, we develop a multiplexed network model, where each city is modelled as a non-linear bi-stable system with stable states in either war or peace. The causal factor categories which exasperate the risk of conflict are each modelled as a network layer. We consider 3 layers: (1) core geospatial network of interacting cities reflecting ground level interactions, (2) cultural network of interacting countries reflecting cultural groupings, and (3) political network of interacting countries reflecting alliances. Together, they act as drivers to push cities towards or pull cities away from war. Using a variety of data sources relative to 2002-2016, we show, that our model correctly predicts the transitions from war to peace and peace to war with F1 score of 0.78 to 0.92 worldwide at the city scale resolution. As many conflicts during this period are auto-regressive (e.g. the War on Terror in Afghanistan and Iraq, the Narco War across the Americas), we can predict the emergence of new war or new peace. We demonstrate successful predictions across a wide range of conflict genres and we perform causal discovery by identifying which model component led to the correct prediction. In the cases of Somalia (2008-13), Myanmar (2013-15), Colombia (2011-14), Libya (2014-16), and Yemen (2011-13) we identify the set of most likely causal factors and how it may differ across a country and change over time.


title: Open-source neuronavigation for multimodal non-invasive brain stimulation using 3D Slicer

abstract: In recent years, non-invasive neuro-modulation methods such as Focused Ultrasound (FUS) have gained popularity. The aim of this work is to introduce the use of existing open-source technology for surgical navigation to the field of multimodal non-invasive brain stimulation. Unlike homegrown and commercial systems, the use of well-documented, well maintained, and freely available open-source components minimizes the learning curve, maximizes technology transfer outcome, and fosters reproducible science for complex, guided neuromodulation systems. The described system significantly lowers the entry bar to clinical research and experimentation in the field of non-invasive brain stimulation. Our contribution is two-fold. First, a high-level overview of the components of the descried system is given in this manuscript. Second, all files are made available online, with a comprehensive step-by-step manual, quickly allowing researchers to build a custom system. A spatial accuracy of 0.93 mm was found through validation using a robotic positioning system.


title: Matrix Completion Using Alternating Minimization for Distribution System State Estimation

abstract: This paper examines the problem of state estimation in power distribution systems under low-observability conditions. The recently proposed constrained matrix completion method which combines the standard matrix completion method and power flow constraints has been shown to be effective in estimating voltage phasors under low-observability conditions using single-snapshot information. However, the method requires solving a semidefinite programming (SDP) problem, which becomes computationally infeasible for large systems and if multiple-snapshot (time-series) information is used. This paper proposes an efficient algorithm to solve the constrained matrix completion problem with time-series data. This algorithm is based on reformulating the matrix completion problem as a bilinear (non-convex) optimization problem, and applying the alternating minimization algorithm to solve this problem. This paper proves the summable convergence of the proposed algorithm, and demonstrates its efficacy and scalability via IEEE 123-bus system and a real utility feeder system. This paper also explores the value of adding more data from the history in terms of computation time and estimation accuracy.


title: Leveraging Multimodal Haptic Sensory Data for Robust Cutting

abstract: Cutting is a common form of manipulation when working with divisible objects such as food, rope, or clay. Cooking in particular relies heavily on cutting to divide food items into desired shapes. However, cutting food is a challenging task due to the wide range of material properties exhibited by food items. Due to this variability, the same cutting motions cannot be used for all food items. Sensations from contact events, e.g., when placing the knife on the food item, will also vary depending on the material properties, and the robot will need to adapt accordingly. In this paper, we propose using vibrations and force-torque feedback from the interactions to adapt the slicing motions and monitor for contact events. The robot learns neural networks for performing each of these tasks and generalizing across different material properties. By adapting and monitoring the skill executions, the robot is able to reliably cut through more than 20 different types of food items and even detect whether certain food items are fresh or old.


title: Constraint Energy Minimizing Generalized Multiscale Discontinuous Galerkin Method

abstract: Numerical simulation of flow problems and wave propagation in heterogeneous media has important applications in many engineering areas. However, numerical solutions on the fine grid are often prohibitively expensive, and multiscale model reduction techniques are introduced to efficiently solve for an accurate approximation on the coarse grid. In this paper, we propose an energy minimization based multiscale model reduction approach in the discontinuous Galerkin discretization setting. The main idea of the method is to extract the non-decaying component in the high conductivity regions by identifying dominant modes with small eigenvalues of local spectral problems, and define multiscale basis functions in coarse oversampled regions by constraint energy minimization problems. The multiscale basis functions are in general discontinuous on the coarse grid and coupled by interior penalty discontinuous Galerkin formulation. The minimal degree of freedom in representing high-contrast features is achieved through the design of local spectral problems, which provides the most compressed local multiscale space. We analyze the method for solving Darcy flow problem and show that the convergence is linear in coarse mesh size and independent of the contrast, provided that the oversampling size is appropriately chosen. Numerical results are presented to show the performance of the method for simulation on flow problem and wave propagation in high-contrast heterogeneous media.


title: $p$-torsion \'etale sheaves on the Jacobian of a curve

abstract: Suppose $X$ is a smooth, proper, geometrically connected curve over $\mathbb F_q$ with an $\mathbb F_q$-rational point $x_0$. For any $\mathbb F_q^{\times}$-character $\sigma$ of $\pi_1(X)$ trivial on $x_0$, we construct a functor $\mathbb L_n^{\sigma}$ from the derived category of coherent sheaves on the moduli space of deformations of $\sigma$ over the Witt ring $W_n(\mathbb F_q)$ to the derived category of constructible $W_n(\mathbb F_q)$-sheaves on the Jacobian of $X$. The functors $\mathbb L_n^{\sigma}$ categorify the Artin reciprocity map for geometric class field theory with $p$-torsion coefficients. We then give a criterion for the fully faithfulness of (an enhanced version of) $\mathbb L_n^{\sigma}$ in terms of the Hasse-Witt matrix of $X$.


title: Color dependence of clustering of massive galaxies at 0.5$\le z \le$2.5: similar spatial distributions between green valley galaxies and AGNs

abstract: We present a measurement of the spatial clustering of rest-frame UV-selected massive galaxies at $0.5\le z \le 2.5$ in the COSMOS/UltraVISTA field. Considering four separate redshift bins with $\Delta z=0.5$, we construct three galaxy populations, i.e., red sequence (RS), blue cloud (BC), and green valley (GV) galaxies, according to their rest-frame extinction-corrected UV colors. The correlation lengths of these populations are confirmed to be dependent on their rest-frame UV color and redshift: UV redder galaxies are found to be more clustered. In all redshift bins, the GV galaxies generally have medium clustering amplitudes and are hosted within dark matter halos whose masses are more or less between those of RS and BC galaxies; and the clustering amplitude of GV galaxies is close to that of AGNs in the same redshift bin, suggesting that AGN activity may be responsible for transforming galaxy colors. After carefully examining their stellar masses, we find that the clustering amplitudes of galaxy samples with different colors are all similar once they have a similar median stellar mass and that the median stellar mass alone may be a good predictor of galaxy clustering.


title: Comment on "$\Phi$ memristor: Real memristor found" by F. Z. Wang, L. Li, L. Shi, H. Wu, and L. O. Chua [J. Appl. Phys. 125, 054504 (2019)]

abstract: Wang et al. claim [J. Appl. Phys. 125, 054504 (2019)] that a current-carrying wire interacting with a magnetic core represents a memristor. Here, we demonstrate that this claim is false. We first show that such memristor "discovery" is based on incorrect physics, which does not even capture basic properties of magnetic core materials, such as their magnetic hysteresis. Moreover, the predictions of Wang et al.'s model contradict the experimental curves presented in their paper. Additionally, the theoretical pinched hysteresis loops presented by Wang et al. can not be reproduced if their model is used, and there are serious flaws in their "negative memristor" emulator design. Finally, a simple gedanken experiment shows that the proposed $\Phi$-memristor would fail the memristor test we recently suggested in J. Phys. D: Appl. Phys. 52, 01LT01 (2019). The device "discovered" by Wang et al. is just an inductor with memory.


title: Playing Atari Ball Games with Hierarchical Reinforcement Learning

abstract: Human beings are particularly good at reasoning and inference from just a few examples. When facing new tasks, humans will leverage knowledge and skills learned before, and quickly integrate them with the new task. In addition to learning by experimentation, human also learn socio-culturally through instructions and learning by example. In this way humans can learn much faster compared with most current artificial intelligence algorithms in many tasks. In this paper, we test the idea of speeding up machine learning through social learning. We argue that in solving real-world problems, especially when the task is designed by humans, and/or for humans, there are typically instructions from user manuals and/or human experts which give guidelines on how to better accomplish the tasks. We argue that these instructions have tremendous value in designing a reinforcement learning system which can learn in human fashion, and we test the idea by playing the Atari games Tennis and Pong. We experimentally demonstrate that the instructions provide key information about the task, which can be used to decompose the learning task into sub-systems and construct options for the temporally extended planning, and dramatically accelerate the learning process.


title: Creating locally interacting Hamiltonians in the synthetic frequency dimension for photons

abstract: The recent emerging field of synthetic dimension in photonics offers a variety of opportunities for manipulating different internal degrees of freedom of photons such as the spectrum of light. While nonlinear optical effects can be incorporated into these photonic systems with synthetic dimensions, these nonlinear effects typically result in long-range interactions along the frequency axis. Thus it has been difficult to use of the synthetic dimension concept to study a large class of Hamiltonians that involves local interactions. Here we show that a Hamiltonian that is locally interacting along the synthetic dimension can be achieved in a dynamically-modulated ring resonator incorporating $\chi^{(3)}$ nonlinearity, provided that the group velocity dispersion of the waveguide forming the ring is specifically designed. As a demonstration we numerically implement a Bose-Hubbard model, and explore photon blockade effect, in the synthetic frequency space. Our work opens new possiblities for studying fundamental many-body physics in the synthetic space in photonics, with potential applications in optical quantum communication and quantum computation.


title: A novel metric for community detection

abstract: Research into detection of dense communities has recently attracted increasing attention within network science, various metrics for detection of such communities have been proposed. The most popular metric -- Modularity -- is based on the so-called rule that the links within communities are denser than external links among communities, has become the default. However, this default metric suffers from ambiguity, and worse, all augmentations of modularity and based on a narrow intuition of what it means to form a "community". We argue that in specific, but quite common systems, links within a community are not necessarily more common than links between communities. Instead we propose that the defining characteristic of a community is that links are more predictable within a community rather than between communities. In this paper, based on the effect of communities on link prediction, we propose a novel metric for the community detection based directly on this feature. We find that our metric is more robustness than traditional modularity. Consequently, we can achieve an evaluation of algorithm stability for the same detection algorithm in different networks. Our metric also can directly uncover the false community detection, and infer more statistical characteristics for detection algorithms.


title: On the Separation Structure of Lalley-Gatzouras Fractals

abstract: We obtain a necessary and sufficient condition for Lalley-Gatzouras sets to be uniform disconnected. This enable us to find all Lalley-Gatzouras sets which are quasisymmetrically equivalent to the Cantor ternary set. As another application, we also study the limit behavior of the gap sequence of Lalley-Gatzouras sets.


title: Telescope: an interactive tool for managing large scale analysis from mobile devices

abstract: In today's world of big data, computational analysis has become a key driver of biomedical research. Recent exponential growth in the volume of available omics data has reshaped the landscape of contemporary biology, creating demand for a continuous feedback loop that seamlessly integrates experimental biology techniques and bioinformatics tools. High-performance computational facilities are capable of processing considerable volumes of data, yet often lack an easy-to-use interface to guide the user in supervising and adjusting bioinformatics analysis in real-time. Here we report the development of Telescope, a novel interactive tool that interfaces with high-performance computational clusters to deliver an intuitive user interface for controlling and monitoring bioinformatics analyses in real-time. Telescope was designed to natively operate with a simple and straightforward interface using Web 2.0 technology compatible with most modern devices (e.g., tablets and personal smartphones). Telescope provides a modern and elegant solution to integrate computational analyses into the experimental environment of biomedical research. Additionally, it allows biomedical researchers to leverage the power of large computational facilities in a user-friendly manner. Telescope is freely available at https://github.com/Mangul-Lab-USC/telescope.


title: On a Class of Sums with Unexpectedly High Cancellation, and its Applications

abstract: Following attempts at an analytic proof of the Pentagonal Number Theorem, we report on the discovery of a general principle leading to an unexpected cancellation of oscillating sums. After stating the motivation, and our theorem, we apply it to prove several results on the Prouhet-Tarry-Escott Problem, integer partitions, and the distribution of prime numbers. Regarding the Prouhet-Tarry-Escott problem, we show that \begin{align*} \sum_{|\ell|\leq x}(4x^2-4\ell^2)^{2r}-\sum_{|\ell|<x}(4x^2-(2\ell+1)^2)^{2r}=\text{polynomial w.r.t. } x \text{ with degree }2r-1. \end{align*} This can perhaps be proved using properties of Bernoulli polynomials, but the claim fell out of our method in a more natural and motivated way. Using this result, we solve an approximate version of the PTE Problem, and in doing so our work in the approximate case exceeds the bounds one can prove using a pigeonhole argument, which seems remarkable. Also, we prove that $$ \sum_{\ell^2 < n} (-1)^\ell p(n-\ell^2)\ \sim\ (-1)^n 2^{-3/4} n^{-1/4} \sqrt{p(n)}, $$ where $p(n)$ is the usual partition function. We get the following "Weak pentagonal number theorem", in which we can replace the partition function $p(n)$ with Chebyshev $\Psi$ function: $$ \sum_{0 < \ell < \sqrt{xT}/2} \Psi([e^{\sqrt{x - \frac{(2\ell)^2}{T}}},\ e^{\sqrt{x - \frac{(2\ell-1)^2}{T}}}])\ =\Psi(e^{\sqrt{x}})\left(\frac{1}{2} + O\left (e^{-0.196\sqrt{x}}\right)\right), $$ where $T=e^{0.786\sqrt{x}}$, where $\Psi([a,b]) := \sum_{n\in [a,b]} \Lambda(n)$ and $\Psi(x) = \Psi([1,x])$, where $\Lambda$ is the von Mangoldt function. Note that this last equation (sum over $\ell$) is stronger than one would get using a strong form of the Prime Number Theorem and also a naive use of the Riemann Hypothesis in each interval, since the widths of the intervals are smaller than $e^{\frac{1}{2} \sqrt{x}}$, making the RH estimate ``trivial".


title: DMM-Net: Differentiable Mask-Matching Network for Video Object Segmentation

abstract: In this paper, we propose the differentiable mask-matching network (DMM-Net) for solving the video object segmentation problem where the initial object masks are provided. Relying on the Mask R-CNN backbone, we extract mask proposals per frame and formulate the matching between object templates and proposals at one time step as a linear assignment problem where the cost matrix is predicted by a CNN. We propose a differentiable matching layer by unrolling a projected gradient descent algorithm in which the projection exploits the Dykstra's algorithm. We prove that under mild conditions, the matching is guaranteed to converge to the optimum. In practice, it performs similarly to the Hungarian algorithm during inference. Meanwhile, we can back-propagate through it to learn the cost matrix. After matching, a refinement head is leveraged to improve the quality of the matched mask. Our DMM-Net achieves competitive results on the largest video object segmentation dataset YouTube-VOS. On DAVIS 2017, DMM-Net achieves the best performance without online learning on the first frames. Without any fine-tuning, DMM-Net performs comparably to state-of-the-art methods on SegTrack v2 dataset. At last, our matching layer is very simple to implement; we attach the PyTorch code ($<50$ lines) in the supplementary material. Our code is released at https://github.com/ZENGXH/DMM_Net.


title: A Radio Signal Modulation Recognition Algorithm Based on Residual Networks and Attention Mechanisms

abstract: To solve the problem of inaccurate recognition of types of communication signal modulation, a RNN neural network recognition algorithm combining residual block network with attention mechanism is proposed. In this method, 10 kinds of communication signals with Gaussian white noise are generated from standard data sets, such as MASK, MPSK, MFSK, OFDM, 16QAM, AM and FM. Based on the original RNN neural network, residual block network is added to solve the problem of gradient disappearance caused by deep network layers. Attention mechanism is added to the network to accelerate the gradient descent. In the experiment, 16QAM, 2FSK and 4FSK are used as actual samples, IQ data frames of signals are used as input, and the RNN neural network combined with residual block network and attention mechanism is trained. The final recognition results show that the average recognition rate of real-time signals is over 93%. The network has high robustness and good use value.


title: Noisy Batch Active Learning with Deterministic Annealing

abstract: We study the problem of training machine learning models incrementally with batches of samples annotated with noisy oracles. We select each batch of samples that are important and also diverse via clustering and importance sampling. More importantly, we incorporate model uncertainty into the sampling probability to compensate for poor estimation of the importance scores when the training data is too small to build a meaningful model. Experiments on benchmark image classification datasets (MNIST, SVHN, CIFAR10, and EMNIST) show improvement over existing active learning strategies. We introduce an extra denoising layer to deep networks to make active learning robust to label noises and show significant improvements.


title: Stratified Space Learning: Reconstructing Embedded Graphs

abstract: Many data-rich industries are interested in the efficient discovery and modelling of structures underlying large data sets, as it allows for the fast triage and dimension reduction of large volumes of data embedded in high dimensional spaces. The modelling of these underlying structures is also beneficial for the creation of simulated data that better represents real data. In particular, for systems testing in cases where the use of real data streams might prove impractical or otherwise undesirable. We seek to discover and model the structure by combining methods from topological data analysis with numerical modelling. As a first step in combining these two areas, we examine the recovery of the abstract graph $G$ structure, and model a linear embedding $|G|$ given only a noisy point cloud sample $X$ of $|G|$.


title: Hidden Stratification Causes Clinically Meaningful Failures in Machine Learning for Medical Imaging

abstract: Machine learning models for medical image analysis often suffer from poor performance on important subsets of a population that are not identified during training or testing. For example, overall performance of a cancer detection model may be high, but the model still consistently misses a rare but aggressive cancer subtype. We refer to this problem as hidden stratification, and observe that it results from incompletely describing the meaningful variation in a dataset. While hidden stratification can substantially reduce the clinical efficacy of machine learning models, its effects remain difficult to measure. In this work, we assess the utility of several possible techniques for measuring and describing hidden stratification effects, and characterize these effects on multiple medical imaging datasets. We find evidence that hidden stratification can occur in unidentified imaging subsets with low prevalence, low label quality, subtle distinguishing features, or spurious correlates, and that it can result in relative performance differences of over 20% on clinically important subsets. Finally, we explore the clinical implications of our findings, and suggest that evaluation of hidden stratification should be a critical component of any machine learning deployment in medical imaging.


title: Distinctive Thermoelectric Properties of Supersaturated Si-Ge-P Compounds: Achieving Figure of Merit ZT > 3.6

abstract: The efficiency of energy conversion in thermoelectric generators (TEGs) is directly proportional to electrical conductivity and Seebeck coefficient while inversely to thermal conductivity. The challenge is to optimize these interdependent parameters simultaneously. In this work, the problem is addressed with a novel approach of nanostructuring and constructive electronic structure modification to achieve a very high value of dimensionless figure of merit ZT greater than 3.6 at 1000 K with negative Seebeck coefficient. Supersaturated solid-solutions of Si-Ge containing 1 atomic percent Fe and 10 atomic percent P are prepared by high-energy ball milling. The bulk samples consisting of ultra-fine nano-crystallites 9.7 nm are obtained by the sophisticated low-temperature & high-pressure sintering process. Despite that the electrical resistivity is slightly high due to the localization of electrons is associated with the highly disordered structure and low electrical density of states near the chemical potential, a very low thermal conductivity \k{appa} less than 1 W m-1K-1 and very large magnitude of Seebeck coefficient exceeding 470 uV K-1 are achieved in association with the nanostructuring and the Fe 3d impurity states, respectively, to realize a very large magnitude of ZT.


title: A right inverse of differential operator $\Delta+a$ in weighted Hilbert space $L^2(\mathbb{R}^n,e^{-|x|^2})$

abstract: In this note, we prove the existence of weak solutions of a Poisson type equation in the weighted Hilbert space $L^2(\mathbb{R}^n,e^{-|x|^2})$.


title: Study of $e^{+}e^{-} \to D^{+} D^{-} \pi^{+} \pi^{-} $ at center-of-mass energies from 4.36 to 4.60 GeV

abstract: We report a study of the $e^{+}e^{-} \to D^{+} D^{-} \pi^{+} \pi^{-}$ process using $e^{+}e^{-}$ collision data samples with an integrated luminosity of $2.5\,\rm{fb}^{-1}$ at center-of-mass energies from 4.36 to $4.60 \rm{GeV}$, collected with the BESIII detector at the BEPCII storage ring. The $D_{1}(2420)^+$ is observed in the $D^{+} \pi^{+} \pi^{-}$ mass spectrum. The mass and width of the $D_{1}(2420)^+$ are measured to be $(2427.2\pm 1.0_{\rm stat.}\pm 1.2_{\rm syst.}) \rm{MeV}/c^2$ and $(23.2\pm 2.3_{\rm stat.} \pm2.3_{\rm syst.}) \rm{MeV}$, respectively. The first errors are statistical and the second ones are systematic. In addition, the Born cross sections of the $e^{+}e^{-} \to D_{1}(2420)^+D^- + c.c. \to D^{+} D^{-} \pi^{+} \pi^{-}$ and $e^{+}e^{-} \to \psi(3770) \pi^{+} \pi^{-} \to D^{+} D^{-} \pi^{+} \pi^{-}$ processes are measured as a function of the center-of-mass energy.


title: Valley separation via trigonal warping

abstract: Monolayer Graphene contains two inequivalent local minimum, valleys, located at $K$ and $K'$ in the Brillouin zone. There has been considerable interest in the use of these two valleys as a doublet for information processing. Herein I propose a method to resolve valley currents spatially, using only a weak magnetic field. Due to the trigonal warping of the valleys, a spatial offset appears in the guiding centre co-ordinate, and is strongly enhanced due to collimation. This can be exploited to spatially separate valley states. Based on current experimental devices, spatial separation is possible for densities well within current experimental limits. Using numerical simulations, I demonstrate the spatial separation of the valley states.


title: Dynamics of time-periodic reaction-diffusion equations with front-like initial data on $\mathbb{R}$

abstract: This paper is concerned with the Cauchy problem $$u_t=u_{xx} +f(t,u), \,\,\, x\in\mathbb{R},\,t>0, $$ $$u(0,x)= u_0(x), \,\,\, x\in\mathbb{R},$$ where $f$ is a rather general nonlinearity that is periodic in $t$, and satisfies $f(\cdot,0)\equiv 0$ and that the corresponding ODE has a positive periodic solution $p(t)$. Assuming that $u_0$ is front-like, that is, $u_0(x)$ is close to $p(0)$ for $x\approx -\infty$ and close to $0$ for $x\approx \infty$, we aim to determine the long-time dynamical behavior of the solution $u(t,x)$ by using the notion of propagation terrace introduced by Ducrot, Giletti and Matano (2014). We establish the existence and uniqueness of propagating terrace for a very large class of nonlinearities, and show the convergence of the solution $u(t,x)$ to the terrace as $t\to\infty$ under various conditions on $f$ or $u_0$. We first consider the special case where $u_0$ is a Heaviside type function, and prove the converge result without requiring any non-degeneracy on $f$. Furthermore, if $u_0$ is more general such that it can be trapped between two Heaviside type functions, but not necessarily monotone, we show that the convergence result remains valid under a rather mild non-degeneracy assumption on $f$. Lastly, in the case where $f$ is a non-degenerate multistable nonlinearity, we show the global and exponential convergence for a much larger class of front-like initial data.


title: Coherent control of nitrogen-vacancy center spins in silicon carbide at room temperature

abstract: Solid-state color centers with manipulatable spin qubits and telecom-ranged fluorescence are ideal platforms for quantum communications and distributed quantum computations. In this work, we coherently control the nitrogen-vacancy (NV) center spins in silicon carbide at room temperature, in which telecom-wavelength emission is detected. We increase the NV concentration six-fold through optimization of implantation conditions. Hence, coherent control of NV center spins is achieved at room temperature and the coherence time T2 can be reached to around 17.1 {\mu}s. Furthermore, investigation of fluorescence properties of single NV centers shows that they are room temperature photostable single photon sources at telecom range. Taking advantages of technologically mature materials, the experiment demonstrates that the NV centers in silicon carbide are promising platforms for large-scale integrated quantum photonics and long-distance quantum networks.


title: Adaptive ROI Generation for Video Object Segmentation Using Reinforcement Learning

abstract: In this paper, we aim to tackle the task of semi-supervised video object segmentation across a sequence of frames where only the ground-truth segmentation of the first frame is provided. The challenges lie in how to online update the segmentation model initialized from the first frame adaptively and accurately, even in presence of multiple confusing instances or large object motion. The existing approaches rely on selecting the region of interest for model update, which however, is rough and inflexible, leading to performance degradation. To overcome this limitation, we propose a novel approach which utilizes reinforcement learning to select optimal adaptation areas for each frame, based on the historical segmentation information. The RL model learns to take optimal actions to adjust the region of interest inferred from the previous frame for online model updating. To speed up the model adaption, we further design a novel multi-branch tree based exploration method to fast select the best state action pairs. Our experiments show that our work improves the state-of-the-art of the mean region similarity on DAVIS 2016 dataset to 87.1%.


title: Mapping with Reflection -- Detection and Utilization of Reflection in 3D Lidar Scans

abstract: This paper presents a method to detect reflection of 3D light detection and ranging (Lidar) scans and uses it to classify the points and also map objects outside the line of sight. Our software uses several approaches to analyze the point cloud, including intensity peak detection, dual return detection, plane fitting, and finding the boundaries. These approaches can classify the point cloud and detect the reflection in it. By mirroring the reflection points on the detected window pane and adding classification labels on the points, we can improve the map quality in a Simultaneous Localization and Mapping (SLAM) framework. Experiments using real scan data and ground truth data showcase the effectiveness of our method.


title: Menger Convexity and Hausdorff Metric

abstract: Shimizu and Takahashi have shown that every decreasing sequence of nonempty, bounded, closed, convex subsets of a complete, uniformly Takahashi convex metric space has nonempty intersection. It is well known that the Menger convexity is a generalization of the Takahashi convexity. In this article, we acquire a nonempty intersection property, in terms of the Hausdorff metric, for Menger convex metric spaces, that also provides a class of reflexive Menger convex spaces. We apply this result to discuss the fixed point results for generalized hybrid mappings defined on Menger convex spaces.


title: Vortex sheets in ideal 3D fluids, coadjoint orbits, and characters

abstract: We describe the coadjoint orbits of the group of volume preserving diffeomorphisms of $\mathbb{R}^3$ associated to the motion of closed vortex sheets in ideal 3D fluids. We show that these coadjoint orbits can be identified with nonlinear Grassmannians of compact surfaces enclosing a given volume and endowed with a closed 1-form describing the vorticity density. If the vorticity density has discrete period group and is nonvanishing, the vortex sheet is given by a surface of genus one fibered by its vortex lines over a circle. We determine the Hamilton equations for such vortex sheets relative to the Hamiltonian function suggested in Khesin (2012) and prove that there are no stationary solutions having rotational symmetries. These coadjoint orbits are shown to be prequantizable if the period group of the 1-form and the volume enclosed by the surface satisfy an Onsager-Feynman relation, as argued in Goldin et al. (1991) for the case of open vortex sheets (tubes/ribbons). We find a character for the prequantizable coadjoint orbits, as well as a polarization group on which the character extends, which is a first step beyond prequantization.


title: Reweighted Proximal Pruning for Large-Scale Language Representation

abstract: Recently, pre-trained language representation flourishes as the mainstay of the natural language understanding community, e.g., BERT. These pre-trained language representations can create state-of-the-art results on a wide range of downstream tasks. Along with continuous significant performance improvement, the size and complexity of these pre-trained neural models continue to increase rapidly. Is it possible to compress these large-scale language representation models? How will the pruned language representation affect the downstream multi-task transfer learning objectives? In this paper, we propose Reweighted Proximal Pruning (RPP), a new pruning method specifically designed for a large-scale language representation model. Through experiments on SQuAD and the GLUE benchmark suite, we show that proximal pruned BERT keeps high accuracy for both the pre-training task and the downstream multiple fine-tuning tasks at high prune ratio. RPP provides a new perspective to help us analyze what large-scale language representation might learn. Additionally, RPP makes it possible to deploy a large state-of-the-art language representation model such as BERT on a series of distinct devices (e.g., online servers, mobile phones, and edge devices).


title: Surface Reconstructions of Heusler Compounds in the Ni-Ti-Sn (001) System

abstract: As progress is made on thin-film synthesis of Heusler compounds, a more complete understanding of the surface will be required to control their properties, especially as functional heterostructures are explored. Here, the surface reconstructions of semiconducting half-Heusler NiTiSn(001), and Ni1+xTiSn(001) (x=0.0-1.0) are explored as a way to optimize growth conditions during molecular beam epitaxy. Density functional theory (DFT) calculations were carried out to guide the interpretation of the experimental results. For NiTiSn(001) a c(2x2) surface reconstruction was observed for Sn rich samples, while a (1x1) unreconstructed surface was observed for Ti-rich samples. A narrow range around 1:1:1 stoichiometry exhibited a (2x1) surface reconstruction. Electrical transport is used to relate the observed reflection high energy electron diffraction (RHEED) pattern during and after growth with carrier concentration and stoichiometry. Scanning tunneling microscopy and RHEED were used to examine surface reconstructions, the results of which are in good agreement with density functional calculations. X-ray photoelectron spectroscopy was used to determine surface termination and stoichiometry. Atomic surface models are proposed, which suggest Sn-dimers form in reconstructed Ni1+xTiSn(001) half-Heusler surfaces (x<0.25) with a transition to Ni terminated surfaces for x > 0.25.


title: Improving Federated Learning Personalization via Model Agnostic Meta Learning

abstract: Federated Learning (FL) refers to learning a high quality global model based on decentralized data storage, without ever copying the raw data. A natural scenario arises with data created on mobile phones by the activity of their users. Given the typical data heterogeneity in such situations, it is natural to ask how can the global model be personalized for every such device, individually. In this work, we point out that the setting of Model Agnostic Meta Learning (MAML), where one optimizes for a fast, gradient-based, few-shot adaptation to a heterogeneous distribution of tasks, has a number of similarities with the objective of personalization for FL. We present FL as a natural source of practical applications for MAML algorithms, and make the following observations. 1) The popular FL algorithm, Federated Averaging, can be interpreted as a meta learning algorithm. 2) Careful fine-tuning can yield a global model with higher accuracy, which is at the same time easier to personalize. However, solely optimizing for the global model accuracy yields a weaker personalization result. 3) A model trained using a standard datacenter optimization method is much harder to personalize, compared to one trained using Federated Averaging, supporting the first claim. These results raise new questions for FL, MAML, and broader ML research.


title: Classification of time-domain waveforms using a speckle-based optical reservoir computer

abstract: Reservoir computing is a recurrent machine learning framework that expands the dimensionality of a problem by mapping an input signal into a higher-dimension reservoir space that can capture and predict features of complex, non-linear temporal dynamics. Here, we report on a bulk optical demonstration of an analog reservoir computer using speckles generated by propagating a laser beam modulated with a spatial light modulator through a multimode waveguide. We demonstrate that the hardware can successfully perform a multivariate audio classification task performed using the Japanese vowel speakers public data set. We perform full wave optical calculations of this architecture implemented in a chip-scale platform using an SiO2 waveguide and demonstrate that it performs as well as a fully numerical implementation of reservoir computing. As all the optical components used in the experiment can be fabricated using a commercial photonic integrated circuit foundry, our result demonstrates a framework for building a scalable, chip-scale, reservoir computer capable of performing optical signal processing.


title: High-order strong methods for stochastic differential equations with colored noises

abstract: The key difficulty to develop efficient high-order methods for integrating stochastic differential equations lies in the calculations of the multiple stochastic integrals. This letter suggests a scheme to compute the stochastic integrals for the colored noises based on the white noise representation. The multiple stochastic integrals involving one and two stationary noises can be conveniently generated together with noises using the discrete Fourier transformation. Based on the calculated stochastic integrals, we obtain simple fourth-order and third-order strong methods for equations with a single and multiple noises, respectively. Numerical tests verify the accuracy of the suggested methods.


title: The status of 21cm interferometric experiments

abstract: Interferometric experiments of the reionization era offer the advantages of measuring power in spatial modes with increased sensitivity afforded by multiple independent sky measurements. Here we review early work to measure this signal, current experiments, and future opportunities, highlighting the lessons learned along the way that have shaped the research field and experimental design. In particular, this chapter discusses the history, progress, challenges and forecasts for detection and exploration of the spatial structure of the 21~cm brightness temperature signal in the Epoch of Reionisation using interferometric experiments. We discuss GMRT, PAPER, LOFAR, MWA, and the future HERA and SKA.


title: Radiation properties of an oscillating atom in the presence of external fields

abstract: In the first part of the present work, the correction to photon emission rate of an oscillating two-level atom in the presence of electromagnetic quantum vacuum field has been investigated for two different configurations: (i) Atom is trapped in the vicinity of a perfect conductor (ii) Atom is trapped between two perfect conductors. In the second part, the correction to the decay rate of an initially excited oscillating two-level atom due to the presence of a perfect conducting surface is found.


title: Invisible Marker: Automatic Annotation of Segmentation Masks for Object Manipulation

abstract: We propose a method to annotate segmentation masks accurately and automatically using invisible marker for object manipulation. Invisible marker is invisible under visible (regular) light conditions, but becomes visible under invisible light, such as ultraviolet (UV) light. By painting objects with the invisible marker, and by capturing images while alternately switching between regular and UV light at high speed, massive annotated datasets are created quickly and inexpensively. We show a comparison between our proposed method and manual annotations. We demonstrate semantic segmentation for deformable objects including clothes, liquids, and powders under controlled environmental light conditions. In addition, we show demonstrations of liquid pouring tasks under uncontrolled environmental light conditions in complex environments such as inside the office, house, and outdoors. Furthermore, it is possible to capture data while the camera is in motion so it becomes easier to capture large datasets, as shown in our demonstration.


title: Deep Gated Multi-modal Learning: In-hand Object Pose Changes Estimation using Tactile and Image Data

abstract: For in-hand manipulation, estimation of the object pose inside the hand is one of the important functions to manipulate objects to the target pose. Since in-hand manipulation tends to cause occlusions by the hand or the object itself, image information only is not sufficient for in-hand object pose estimation. Multiple modalities can be used in this case, the advantage is that other modalities can compensate for occlusion, noise, and sensor malfunctions. Even though deciding the utilization rate of a modality (referred to as reliability value) corresponding to the situations is important, the manual design of such models is difficult, especially for various situations. In this paper, we propose deep gated multi-modal learning, which self-determines the reliability value of each modality through end-to-end deep learning. For the experiments, an RGB camera and a GelSight tactile sensor were attached to the parallel gripper of the Sawyer robot, and the object pose changes were estimated during grasping. A total of 15 objects were used in the experiments. In the proposed model, the reliability values of the modalities were determined according to the noise level and failure of each modality, and it was confirmed that the pose change was estimated even for unknown objects.


title: Valley filters, accumulators, and switches induced in graphene quantum dots by lines of adsorbed hydrogen atoms

abstract: We present electronic structure and quantum transport calculations that predict conducting channels induced in graphene quantum dots by lines of adsorbed hydrogen atoms to function as highly efficient, experimentally realizable valley filters, accumulators and switches. The underlying physics is a novel property of graphene Dirac point resonances (DPRs) that is revealed here, namely, that an electric current passing through a DPR-mediated conducting channel in a given direction is carried by electrons of em only one of the two graphene valleys. Our predictions apply to lines of hydrogen atoms adsorbed on graphene quantum dots that are either free standing or supported on a hexagonal boron nitride substrate.


title: Novel Reconciliation Protocol Based on Spinal Code for Continuous-variable Quantum Key Distribution

abstract: Reconciliation is a crucial procedure in post-processing of continuous variable quantum key distribution (CV-QKD) system, which is used to make two distant legitimate parties share identical corrected keys. The adaptive reconciliation is necessary and important for practical systems to cope with the variable channel. Many researchers adopt the punctured LDPC codes to implement adaptive reconciliation. In this paper, a novel rateless reconciliation protocol based on spinal code is proposed, which can achieve a high-efficiency and adaptive reconciliation in a larger range of SNRs. Due to the short codes length and simple tructure, our protocol is easy to implement without the complex codes designs of fixed rate codes, e.g., LDPC codes. Meanwhile, the structure of our protocol is highly parallel, which is suitable for hardware implementation, thus it also has the potential of high-speed hardware implementation. Besides, the security of proposed protocol is proved in theory. Experiment results show that the reconciliation efficiency maintains around 95% for ranging SNRs in a larger range (0,0.5), even exceeds 96.5% at extremely low SNR (<= 0.03) by using this novel scheme. The proposed protocol makes the long-distance CV-QKD systems much easier and stable to perform a high-performance and adaptive reconciliation.


title: Edge Expansion and Spectral Gap of Nonnegative Matrices

abstract: The classic graphical Cheeger inequalities state that if $M$ is an $n\times n$ symmetric doubly stochastic matrix, then \[ \frac{1-\lambda_{2}(M)}{2}\leq\phi(M)\leq\sqrt{2\cdot(1-\lambda_{2}(M))} \] where $\phi(M)=\min_{S\subseteq[n],|S|\leq n/2}\left(\frac{1}{|S|}\sum_{i\in S,j\not\in S}M_{i,j}\right)$ is the edge expansion of $M$, and $\lambda_{2}(M)$ is the second largest eigenvalue of $M$. We study the relationship between $\phi(A)$ and the spectral gap $1-\text{Re}\lambda_{2}(A)$ for any doubly stochastic matrix $A$ (not necessarily symmetric), where $\lambda_{2}(A)$ is a nontrivial eigenvalue of $A$ with maximum real part. Fiedler showed that the upper bound on $\phi(A)$ is unaffected, i.e., $\phi(A)\leq\sqrt{2\cdot(1-\text{Re}\lambda_{2}(A))}$. With regards to the lower bound on $\phi(A)$, there are known constructions with \[ \phi(A)\in\Theta\left(\frac{1-\text{Re}\lambda_{2}(A)}{\log n}\right), \] indicating that at least a mild dependence on $n$ is necessary to lower bound $\phi(A)$. In our first result, we provide an exponentially better construction of $n\times n$ doubly stochastic matrices $A_{n}$, for which \[\phi(A_{n})\leq\frac{1-\text{Re}\lambda_{2}(A_{n})}{\sqrt{n}}.\] In fact, all nontrivial eigenvalues of our matrices are $0$, even though the matrices are highly nonexpanding. We further show that this bound is in the correct range (up to the exponent of $n$), by showing that for any doubly stochastic matrix $A$, \[\phi(A)\geq\frac{1-\text{Re}\lambda_{2}(A)}{35\cdot n}.\] Our second result extends these bounds to general nonnegative matrices $R$, obtaining a two-sided quantitative refinement of the Perron-Frobenius theorem in which the edge expansion $\phi(R)$ (appropriately defined), a quantitative measure of the irreducibility of $R$, controls the gap between the Perron-Frobenius eigenvalue and the next-largest real part of any eigenvalue.


title: The Convex Geometry of Integrator Reach Sets

abstract: We study the convex geometry of the forward reach sets for integrator dynamics in finite dimensions with bounded control. We derive closed-form expressions for the volume and the diameter (i.e., maximal width) of these sets in terms of the state space dimension, control bound, and time. These results are novel, and use convex analysis to give an analytical handle on the "size" of the integrator reach set. Several concrete examples are provided to illustrate our results. We envision that the ideas presented here will motivate further theoretical and algorithmic development in reach set computation.


title: Risk-Averse Planning Under Uncertainty

abstract: We consider the problem of designing policies for partially observable Markov decision processes (POMDPs) with dynamic coherent risk objectives. Synthesizing risk-averse optimal policies for POMDPs requires infinite memory and thus undecidable. To overcome this difficulty, we propose a method based on bounded policy iteration for designing stochastic but finite state (memory) controllers, which takes advantage of standard convex optimization methods. Given a memory budget and optimality criterion, the proposed method modifies the stochastic finite state controller leading to sub-optimal solutions with lower coherent risk.


title: Research into the sampling methods of digital beam position measurement

abstract: BPM (Beam Position Measurement) system is one of the most important beam diagnostic instruments in accelerators. A fully digital BPM (DBPM) has been designed for SSRF (Shanghai Synchrotron Radiation Facility). As Analog-to-Digital Converter (ADC) is one crucial part in the DBPM system, the sampling methods should be studied to achieve optimum performance. We implemented different sampling modes and compared them through tests. Meanwhile, long term variation among four sampling channels is another concern, which would introduce errors in beam position measurement. We designed an interleaved distribution scheme to address this issue. To evaluate these sampling methods, we conducted commissioning tests with the beam in SSRF. Test results indicate that with proper sampling methods, a turn-by-turn (TBT) position resolution better than 1 um is achieved, and the slow-acquisition (SA) position resolution is enhanced from 4.28 {\mu}m to 0.17 {\mu}m successfully with interleaved distribution applied, both beyond requirement.


title: Dynamics in a time-discrete food-chain model with strong pressure on preys

abstract: Ecological systems are complex dynamical systems. Modelling efforts on ecosystems' dynamical stability have revealed that population dynamics, being highly nonlinear, can be governed by complex fluctuations. Indeed, experimental and field research has provided mounting evidence of chaos in species' abundances, especially for discrete-time systems. Discrete-time dynamics, mainly arising in boreal and temperate ecosystems for species with non-overlapping generations, have been largely studied to understand the dynamical outcomes due to changes in relevant ecological parameters. The local and global dynamical behaviour of many of these models is difficult to investigate analytically in the parameter space and, typically, numerical approaches are employed when the dimension of the phase space is large. In this article we provide topological and dynamical results for a map modelling a discrete-time, three-species food chain with two predator species interacting on the same prey population. The domain where dynamics live is characterized, as well as the so-called escaping regions, for which the species go rapidly to extinction after surpassing the carrying capacity. We also provide a full description of the local stability of equilibria within a volume of the parameter space given by the prey's growth rate and the predation rates. We have found that the increase of the pressure of predators on the prey results in chaos. The entry into chaos is achieved via a supercritical Neimarck-Sacker bifurcation followed by period-doubling bifurcations of invariant curves. Interestingly, an increasing predation directly on preys can shift the extinction of top predators to their survival, allowing an unstable persistence of the three species by means of periodic and strange chaotic attractors.


title: Bootstrap Cross-validation Improves Model Selection in Pharmacometrics

abstract: Cross-validation assesses the predictive ability of a model, allowing one to rank models accordingly. Although the nonparametric bootstrap is almost always used to assess the variability of a parameter, it can be used as the basis for cross-validation if one keeps track of which items were not selected in a given bootstrap iteration. The items which were selected constitute the training data and the omitted items constitute the testing data. This bootstrap cross-validation (BS-CV) allows model selection to be made on the basis of predictive ability by comparing the median values of ensembles of summary statistics of testing data. BS-CV is herein demonstrated using several summary statistics, including a new one termed the simple metric for prediction quality (SMPQ), and using the warfarin data included in the Monolix distribution with 13 pharmacokinetics (PK) models and 12 pharmacodynamics (PD) models. Of note the two best PK models by AIC had the worst predictive ability, underscoring the danger of using single realizations of a random variable (such as AIC) as the basis for model selection. Using these data BS-CV was able to discriminate between similar indirect response models (inhibition of input versus stimulation of output). This could be useful in situations in which the mechanism of action is unknown (unlike warfarin).


title: x_opt and x_c in Cuprate Superconductors

abstract: The characteristic hole concentrations x_opt = 0.16 and x_c = 0.19 in the under- to overdoped transition regime of cuprate superconductors are shown to be intimately interrelated by mesoscale organization of bond-like RVB arrays extending over 3 x 4 copper sites and comprising 2 holes spaced by 3a_0. Columnar organization maximizes T_c but reduces the superfluid stiffness, while zigzag (tweedy) organization maximizes the superfluid stiffness at reduced T_c. The latter goes hand in hand with the closure of the pseudogap. Columnar and zigzag mesoscale organizations of the bond-like RVB arrays are displacively transformed in each other between x_opt and x_c by weak overdoping.


title: Gallium-Boron-Phosphide (GaBP$_2$): A New III-V Semiconductor for photovoltaics

abstract: Using machine learning (ML) approach, we unearthed a new III-V semiconducting material having an optimal bandgap for high efficient photovoltaics with the chemical composition of Gallium-Boron-Phosphide(GaBP$_2$, space group: Pna2$_1$). ML predictions are further validated by state of the art ab-initio density functional theory (DFT) simulations. The stoichiometric Heyd-Scuseria-Ernzerhof (HSE) bandgap of GaBP$_2$ is noted to 1.65 eV, a close ideal value (1.4-1.5 eV) to reach the theoretical Queisser-Shockley limit. The calculated electron mobility is similar to that of silicon. Unlike perovskites, the newly discovered material is thermally, dynamically and mechanically stable. Above all the chemical composition of GaBP$_2$ are non-toxic and relatively earth-abundant, making it a new generation of PV material. Using ML, we show that with a minimal set of features the bandgap of III-III-V and II-IV-V semiconductor can be predicted up to an RMSE of less than 0.4 eV. We presented a set of scaling laws, which can be used to estimate the bandgap of new III-III-V and II-IV-V semiconductor, with three different crystal phases, within an RMSE of approx. 0.5 eV.


title: Anderson localization induced by complex potential

abstract: Uncorrelated disorder potential in one-dimensional lattice definitely induces Anderson localization, while quasiperiodic potential can lead to both localized and extended phases, depending on the potential strength. We investigate the Anderson localization in one-dimensional lattice with nonHermitian complex disorder and quasiperiodic potential. We present a non-Hermitian SSH chain to demonstrate that a nonzero imaginary disorder on-site potential can induce the standard Anderson localization. We also show that the non-Hermitian Aubry-Andre (AA) model exhibits a transition in parameter space, which separates the localization and delocalization phases and is determined by the self-duality of the model. It indicates that a pure imaginary quasiperiodic potential takes the same role as a real quasiperiodic potential in the transition point between localization and delocalization. Remarkably, a system with complex quasiperiodic potential exhibits interference-like pattern on the transition points.


title: Distributionally Robust Tuning of Anomaly Detectors in Cyber-Physical Systems with Stealthy Attacks

abstract: Designing resilient control strategies for mitigating stealthy attacks is a crucial task in emerging cyber-physical systems. In the design of anomaly detectors, it is common to assume Gaussian noise models to maintain tractability; however, this assumption can lead the actual false alarm rate to be significantly higher than expected. We propose a distributionally robust anomaly detector for noise distributions in moment-based ambiguity sets. We design a detection threshold that guarantees that the actual false alarm rate is upper bounded by the desired one by using generalized Chebyshev inequalities. Furthermore, we highlight an important trade-off between the worst-case false alarm rate and the potential impact of a stealthy attacker by efficiently computing an outer ellipsoidal bound for the attack-reachable states corresponding to the distributionally robust detector threshold. We illustrate this trade-off with a numerical example and compare the proposed approach with a traditional chi-squared detector.


title: Region-wise Generative Adversarial ImageInpainting for Large Missing Areas

abstract: Recently deep neutral networks have achieved promising performance for filling large missing regions in image inpainting tasks. They usually adopted the standard convolutional architecture over the corrupted image, leading to meaningless contents, such as color discrepancy, blur and artifacts. Moreover, most inpainting approaches cannot well handle the large continuous missing area cases. To address these problems, we propose a generic inpainting framework capable of handling with incomplete images on both continuous and discontinuous large missing areas, in an adversarial manner. From which, region-wise convolution is deployed in both generator and discriminator to separately handle with the different regions, namely existing regions and missing ones. Moreover, a correlation loss is introduced to capture the non-local correlations between different patches, and thus guides the generator to obtain more information during inference. With the help of our proposed framework, we can restore semantically reasonable and visually realistic images. Extensive experiments on three widely-used datasets for image inpainting tasks have been conducted, and both qualitative and quantitative experimental results demonstrate that the proposed model significantly outperforms the state-of-the-art approaches, both on the large continuous and discontinuous missing areas.


title: Investigations of Interacting Quantum Systems in a Novel General Approximation Scheme (NGAS)

abstract: A novel general approximation scheme (NGAS) is proposed, which is potentially applicable for arbitrary strength of interaction in quantum theory. The scheme utilizes an input Hamiltonian, which is exactly solvable. The effects of interaction are incorporated into the input Hamiltonian through a non-linear feedback mechanism by self-consistency requirements. The method is nonperturbative, yet improvable in a new formulation of perturbation theory, designated as: mean field perturbation theory (MFPT). We apply the scheme to anharmonic interactions in one dimension and obtain uniformly accurate results by Borel summation, for the quartic-, sextic-, octic-anharmonic oscillators and to the quartic double-well oscillator (QDWO) for arbitrary strength of coupling. The results for the QDWO may be contrasted with those from the standard formulation of perturbation theory , where Borel-summation fails even for small values of the coupling strength.


title: Study on possible molecular states composed of $\Lambda_c\bar D$ ($\Lambda_b B$) and $\Sigma_c\bar D$ ($\Sigma_b B$) within the Bethe-Salpeter framework

abstract: $P_c(4312)$ observed by the LHCb collaboration is confirmed as a pentaquark and its structure, production, and decay behaviors attract great attention from theorists and experimentalists. Since its mass is very close to sum of $\Sigma_c$ and $\bar D$ masses, it is naturally tempted to be considered as a molecular state composed of $\Sigma_c$ and $\bar D$. Moreover, $P_c(4312)$ is observed in the channel with $J/\psi p$ final state, requiring that isospin conservation $P_c(4312)$ is an isospin-1/2 eigenstate. In literature, several groups used various models to estimate its spectrum. We systematically study the pentaquarks within the framework of the Bethe-Salpeter equation; thus $P_c(4312)$ is an excellent target because of the available data. We calculate the spectrum of $P_c(4312)$ in terms of the Bethe-Salpter equations and further study its decay modes. Some predictions on other possible pentaquark states that can be tested in future experiments are made.


title: The Two-sided Jet Structures of NGC 1052 at Scales from 300 to $4 \times 10^7$ Schwarzschild Radii

abstract: We investigated the jet width profile with distance along the jet in the nearby radio galaxy NGC 1052 at radial distances between $\sim300$ to $4 \times 10^7$ Schwarzschild Radii($R_{\rm S}$) from the central engine on both their approaching and receding jet sides. The width of jets was measured in images obtained with the VLBI Space Observatory Programme (VSOP), the Very Long Baseline Array (VLBA), and the Very Large Array (VLA). The jet-width profile of receding jets are apparently consistent with that of approaching jets throughout the measuring distance ranges, indicating symmetry at least up to the sphere of gravitational influence of the central black hole. The power-law index $a$ of the jet-width profile ($w_{\rm{jet}} \propto r^{a}$, where $w_{\rm jet}$ is the jet width, $r$ is the distance from the central engine in the unit of $R_{\rm S}$) apparently shows a transition from $a \sim 0$ to $a \sim 1$, i.e., the cylindrical-to-conical jet structures, at a distance of $\sim1\times10^{4} \ R_{\mathrm{S}}$. The cylindrical jet shape at the small distances is reminiscent of the innermost jets in 3C 84. Both the central engines of NGC 1052 and 3C 84 are surrounded by dense material, part of which is ionized and causes heavy free-free absorption.


title: Finite Horizon Density Steering for Multi-input State Feedback Linearizable Systems

abstract: In this paper, we study the feedback synthesis problem for steering the joint state density or ensemble subject to multi-input state feedback linearizable dynamics. This problem is of interest to many practical applications including that of dynamically shaping a robotic swarm. Our results here show that it is possible to exploit the structural nonlinearities to derive the feedback controllers steering the joint density from a prescribed shape to another while minimizing the expected control effort to do so. The developments herein build on our previous work, and extend the theory of the Schr\"{o}dinger bridge problem subject to feedback linearizable dynamics.


title: On families of optimal Hardy-weights for linear second-order elliptic operators

abstract: We construct families of optimal Hardy-weights for a subcritical linear second-order elliptic operator using a one-dimensional reduction. More precisely, we first characterise all optimal Hardy-weights with respect to one-dimensional subcritical Sturm-Liouville operators on a given interval, and then apply this result to obtain families of optimal Hardy inequalities for general linear second-order elliptic operators in higher dimensions. As an application, we prove a new Rellich inequality.


title: Learnable Tree Filter for Structure-preserving Feature Transform

abstract: Learning discriminative global features plays a vital role in semantic segmentation. And most of the existing methods adopt stacks of local convolutions or non-local blocks to capture long-range context. However, due to the absence of spatial structure preservation, these operators ignore the object details when enlarging receptive fields. In this paper, we propose the learnable tree filter to form a generic tree filtering module that leverages the structural property of minimal spanning tree to model long-range dependencies while preserving the details. Furthermore, we propose a highly efficient linear-time algorithm to reduce resource consumption. Thus, the designed modules can be plugged into existing deep neural networks conveniently. To this end, tree filtering modules are embedded to formulate a unified framework for semantic segmentation. We conduct extensive ablation studies to elaborate on the effectiveness and efficiency of the proposed method. Specifically, it attains better performance with much less overhead compared with the classic PSP block and Non-local operation under the same backbone. Our approach is proved to achieve consistent improvements on several benchmarks without bells-and-whistles. Code and models are available at https://github.com/StevenGrove/TreeFilter-Torch.


title: Clustering Uncertain Data via Representative Possible Worlds with Consistency Learning

abstract: Clustering uncertain data is an essential task in data mining for the internet of things. Possible world based algorithms seem promising for clustering uncertain data. However, there are two issues in existing possible world based algorithms: (1) They rely on all the possible worlds and treat them equally, but some marginal possible worlds may cause negative effects. (2) They do not well utilize the consistency among possible worlds, since they conduct clustering or construct the affinity matrix on each possible world independently. In this paper, we propose a representative possible world based consistent clustering (RPC) algorithm for uncertain data. First, by introducing representative loss and using Jensen-Shannon divergence as the distribution measure, we design a heuristic strategy for the selection of representative possible worlds, thus avoiding the negative effects caused by marginal possible worlds. Second, we integrate a consistency learning procedure into spectral clustering to deal with the representative possible worlds synergistically, thus utilizing the consistency to achieve better performance. Experimental results show that our proposed algorithm performs better than the state-of-the-art algorithms.


title: Actions of homeomorphism groups of manifolds admitting a nontrivial finite free action

abstract: In this paper, we study the action of $\text{Homeo}_0(M)$, the identity component of the group of homeomorphisms of an $n$-dimensional manifold $M$ with an $\mathbb{F}_p$-free action, on another manifold $N$ of dimension $n+k<2n$. We prove that if $M$ is not an $\mathbb{F}_p$-homology sphere, then $N\cong M\times K$ for a homology manifold $K$ such that the action of $\text{Homeo}_0(M)$ on $M$ is standard and on $K$ is trivial. In particular, for $M=S^n$ a sphere, any nontrivial action is a generalization of the "coning-off" construction.


title: Non-Hermitian Linear Response Theory

abstract: Linear response theory lies at the heart of quantum many-body physics because it builds up connections between the dynamical response to an external probe and correlation functions at equilibrium. Here we consider the dynamical response of a Hermitian system to a non-Hermitian probe, and we develop a non-Hermitian linear response theory that can also relate this dynamical response to equilibrium properties. As an application of our theory, we consider the real-time dynamics of momentum distribution induced by one-body and two-body dissipations. We find that, for many cases, the dynamics of momentum occupation and the width of momentum distribution obey the same universal function, governed by the single-particle spectral function. We also find that, for critical state with no well-defined quasi-particles, the dynamics are slower than normal state and our theory provides a model independent way to extract the critical exponent. We apply our results to analyze recent experiment on the Bose-Hubbard model and find surprising good agreement between theory and experiment. We also propose to further verify our theory by carrying out a similar experiment on a one-dimensional Luttinger liquid.


title: TORM: Fast and Accurate Trajectory Optimization of Redundant Manipulator given an End-Effector Path

abstract: A redundant manipulator has multiple inverse kinematics solutions per end-effector pose. Accordingly, there can be many trajectories for joints that follow a given endeffector path in the Cartesian space. In this paper, we present a trajectory optimization of a redundant manipulator (TORM) to synthesize a trajectory that follows a given end-effector path accurately, while achieving smoothness and collisionfree manipulation. Our method holistically incorporates three desired properties into the trajectory optimization process by integrating the Jacobian-based inverse kinematics solving method and an optimization-based motion planning approach. Specifically, we optimize a trajectory using two-stage gradient descent to reduce potential competition between different properties during the update. To avoid falling into local minima, we iteratively explore different candidate trajectories with our local update. We compare our method with state-of-the-art methods in test scenes including external obstacles and two non-obstacle problems. Our method robustly minimizes the pose error in a progressive manner while satisfying various desirable properties.


title: Margin-Based Generalization Lower Bounds for Boosted Classifiers

abstract: Boosting is one of the most successful ideas in machine learning. The most well-accepted explanations for the low generalization error of boosting algorithms such as AdaBoost stem from margin theory. The study of margins in the context of boosting algorithms was initiated by Schapire, Freund, Bartlett and Lee (1998) and has inspired numerous boosting algorithms and generalization bounds. To date, the strongest known generalization (upper bound) is the $k$th margin bound of Gao and Zhou (2013). Despite the numerous generalization upper bounds that have been proved over the last two decades, nothing is known about the tightness of these bounds. In this paper, we give the first margin-based lower bounds on the generalization error of boosted classifiers. Our lower bounds nearly match the $k$th margin bound and thus almost settle the generalization performance of boosted classifiers in terms of margins.


title: Extreme Spin Squeezing from Deep Reinforcement Learning

abstract: Spin squeezing (SS) is a recognized resource for realizing measurement precision beyond the standard quantum limit $\propto 1/\sqrt{N}$. The rudimentary one-axis twisting (OAT) interaction can facilitate SS and has been realized in diverse experiments, but it cannot achieve extreme SS for precision at Heisenberg limit $\propto 1/{N}$. Aided by deep reinforcement learning (DRL), we discover size-independent universal rules for realizing nearly extreme SS with OAT interaction using merely a handful of rotation pulses. More specifically, only 6 pairs of pulses are required for up to $10^4$ particles, while the time taken to reach extreme SS remains on the same order of the optimal OAT squeezing time, which makes our scheme viable for experiments that reported OAT squeezing. This study highlights the potential of DRL for controlled quantum dynamics.


title: Online Learning of Dynamical Systems: An Operator Theoretic Approach

abstract: In this paper, we provide an algorithm for online computation of Koopman operator in real-time using streaming data. In recent years, there has been an increased interest in data-driven analysis of dynamical systems, with operator theoretic techniques being the most popular. Existing algorithms, like Dynamic Mode Decomposition (DMD) and Extended Dynamic Mode Decomposition (EDMD), use the entire data set for computation of the Koopman operator. However, many real life applications like power system analysis, biological systems, building systems etc. requires the real-time computation and updating of the Koopman operator, as new data streams in. In this paper, we propose an iterative algorithm for online computation of Koopman operator such that at each time step the Koopman operator is updated incrementally. In particular, we propose a Recursive Extended Dynamic Decomposition (rEDMD) algorithm for computation of Koopman operator from streaming data. Further, we test the algorithm in three different dynamical systems, namely, a linear system, a nonlinear system and a system governed by a Partial Differential Equation (PDE) and illustrate the computational efficiency of the iterative algorithm over the existing DMD and EDMD algorithms.


title: Comparing Static and Dynamic Weighted Software Coupling Metrics

abstract: Coupling metrics are an established way to measure software architecture quality with respect to modularity. Static coupling metrics are obtained from the source or compiled code of a program, while dynamic metrics use runtime data gathered e.g., by monitoring a system in production. We study \emph{weighted} dynamic coupling that takes into account how often a connection is executed during a system's run. We investigate the correlation between dynamic weighted metrics and their static counterparts. We use data collected from four different experiments, each monitoring production use of a commercial software system over a period of four weeks. We observe an unexpected level of correlation between the static and the weighted dynamic case as well as revealing differences between class- and package-level analyses.


title: Strongly constrained stochastic processes: the multi-ends Brownian bridge

abstract: In a recent article, Krapivsky and Redner (J. Stat. Mech. 093208 (2018)) established that the distribution of the first hitting times for a diffusing particle subject to hitting an absorber is independent of the direction of the external flow field. In the present paper, we build upon this observation and investigate when the conditioning on the diffusion leads to a process that is totally independent of the flow field. For this purpose, we adopt the Langevin approach, or more formally the theory of conditioned stochastic differential equations. This technique allows us to derive a large variety of stochastic processes: in particular, we introduce a new kind of Brownian bridge ending at two different final points and calculate its fundamental probabilities. This method is also very well suited for generating statistically independent paths. Numerical simulations illustrate our findings.


title: Semi-classical electronic transport properties of ternary compound AlGaAs$_2$: Role of different scattering mechanisms

abstract: We present a comprehensive investigation of semi-classical transport properties of n-type ternary compound AlGaAs2, using Rode's iterative method. Four scattering mechanisms, have been included in our transport calculation, namely, ionized impurity, piezoelectric, acoustic deformation and polar optical phonon (POP). The scattering rates have been calculated in terms of ab-initio parameters. We consider AlGaAs2 to have two distinct crystal geometries, one in tetragonal phase (space group: ), while the other one having body centered tetragonal crystal structure (space group:). We have observed higher electron mobility in the body centered tetragonal phase, thereby making it more suitable for high mobility device application, over the tetragonal phase. In order to understand the differences in electron moblities for these two phases, curvatures of the E-k graph of the conduction bands for these phases have been compared. At room temperature, the dominant contribution in electron mobility was found to be provided by inelastic POP scattering. We have also noted that mobility is underestimated in relaxation time approximation as compared with the Rode's iterative approach.