title: Learned Critical Probabilistic Roadmaps for Robotic Motion Planning

abstract: Sampling-based motion planning techniques have emerged as an efficient algorithmic paradigm for solving complex motion planning problems. These approaches use a set of probing samples to construct an implicit graph representation of the robot's state space, allowing arbitrarily accurate representations as the number of samples increases to infinity. In practice, however, solution trajectories only rely on a few critical states, often defined by structure in the state space (e.g., doorways). In this work we propose a general method to identify these critical states via graph-theoretic techniques (betweenness centrality) and learn to predict criticality from only local environment features. These states are then leveraged more heavily via global connections within a hierarchical graph, termed Critical Probabilistic Roadmaps. Critical PRMs are demonstrated to achieve up to three orders of magnitude improvement over uniform sampling, while preserving the guarantees and complexity of sampling-based motion planning. A video is available at https://youtu.be/AYoD-pGd9ms.


title: A remark on the smallest singular value of powers of Gaussian matrices

abstract: Let $n,k\geq 1$ and let $G$ be the $n\times n$ random matrix with i.i.d. standard real Gaussian entries. We show that there are constants $c_k,C_k>0$ depending only on $k$ such that the smallest singular value of $G^k$ satisfies $$ c_k\,t\leq {\mathbb P}\big\{s_{\min}(G^k)\leq t^k\,n^{-1/2}\big\}\leq C_k\,t,\quad t\in(0,1], $$ and, furthermore, $$ c_k/t\leq {\mathbb P}\big\{\|G^{-k}\|_{HS}\geq t^k\,n^{1/2}\big\}\leq C_k/t,\quad t\in[1,\infty), $$ where $\|\cdot\|_{HS}$ denotes the Hilbert-Schmidt norm.


title: Molecular Beam Epitaxy of lithium niobium oxide multifunctional materials

abstract: The role of stoichiometry and growth temperature in the preferential nucleation of material phases in the Li-Nb-O family are explored yielding an empirical growth phase diagram. It is shown that while single parameter variation often produces multi-phase films, combining substrate temperature control with the previously published lithium flux limited growth allows the repeatable growth of high quality single crystalline films of many different oxide phases. Higher temperatures (800-1050 {\deg}C) than normally used in MBE were necessary to achieve high quality materials. At these temperatures the desorption of surface species is shown to play an important role in film composition. Using this method single phase films of NbO, NbO$_{2}$, LiNbO$_{2}$, Li$_{3}$NbO$_{4}$, LiNbO$_{3}$, and LiNb$_{3}$O$_{8}$ have been achieved in the same growth system, all on c-plane sapphire. Finally, the future of these films in functional oxide heterostructures is briefly discussed.


title: Do People Prefer "Natural" code?

abstract: Natural code is known to be very repetitive (much more so than natural language corpora); furthermore, this repetitiveness persists, even after accounting for the simpler syntax of code. However, programming languages are very expressive, allowing a great many different ways (all clear and unambiguous) to express even very simple computations. So why is natural code repetitive? We hypothesize that the reasons for this lie in fact that code is bimodal: it is executed by machines, but also read by humans. This bimodality, we argue, leads developers to write code in certain preferred ways that would be familiar to code readers. To test this theory, we 1) model familiarity using a language model estimated over a large training corpus and 2) run an experiment applying several meaning preserving transformations to Java and Python expressions in a distinct test corpus to see if forms more familiar to readers (as predicted by the language models) are in fact the ones actually written. We find that these transformations generally produce program structures that are less common in practice, supporting the theory that the high repetitiveness in code is a matter of deliberate preference. Finally, 3) we use a human subject study to show alignment between language model score and human preference for the first time in code, providing support for using this measure to improve code.


title: Molecular Beam Epitaxy Growth of High Crystalline Quality LiNbO$_{3}$

abstract: Lithium niobate is a multi-functional material with wide reaching applications in acoustics, optics, and electronics. Commercial applications for lithium niobate require high crystalline quality currently limited to bulk and ion sliced material. Thin film lithium niobate is an attractive option for a variety of integrated devices, but the research effort has been stagnant due to poor material quality. Both lattice matched and mismatched lithium niobate are grown by molecular beam epitaxy (MBE) and studied to understand the role of substrate and temperature on nucleation conditions and material quality. Growth on sapphire produces partially coalesced columnar grains with atomically flat plateaus and no twin planes. A symmetric rocking curve shows a narrow linewidth with a full width at half-maximum (FWHM) of 8.6 arcsec (0.0024{\deg}) which is comparable to the 5.8 arcsec rocking curve FWHM of the substrate, while the film asymmetric rocking curve is 510 arcsec FWHM. These values indicate that the individual grains are relatively free of long-range disorder detectable by x-ray diffraction (XRD) with minimal measurable tilt and twist and represents the highest structural quality epitaxial material grown on lattice mismatched sapphire without twin planes. Lithium niobate is also grown on lithium tantalate producing high quality coalesced material without twin planes and with a symmetric rocking curve of 193 arcsec, which is nearly equal to the substrate rocking curve of 194 arcsec. The surface morphology of lithium niobate on lithium tantalate is shown to be atomically flat by atomic force microscopy (AFM).


title: Experimental demonstration of an electrostatic orbital angular momentum sorter for electrons

abstract: We report the first experimental demonstration of an electrostatic electron orbital angular momentum (OAM) sorter, which can be used to analyze the OAM states of electrons in a transmission electron microscope. We verify the sorter functionality for several electron beams possessing different superpositions of OAM states, and use it to record the electron beams OAM spectra. Our current electrostatic OAM sorter has an OAM resolution of 2 in the units of h/bar - the reduced Planck constant. It is expected to increase the OAM resolution of the sorter to the optimal resolution of 1 in the future via fine control of the sorting phase elements.


title: Double parton distribution of valence quarks in the pion in chiral quark models

abstract: The valence double parton distribution of the pion is analyzed in the framework of chiral quark models, where in the chiral limit factorization between the longitudinal and transverse degrees of freedom occurs. This feature leads, at the quark-model scale, to a particularly simple distribution of the form $D(x_1,x_2, \vec{q}) = \delta(1-x_1-x_2) F(\vec{q})$, where $x_{1,2}$ are the longitudinal momentum fractions carried the valence quark and antiquark and $\vec{q}$ is their relative transverse momentum. For $\vec{q}=\vec{0}$ this result complies immediately to the Gaunt-Sterling sum rules. The DGLAP evolution to higher scales is carried out in terms of the Mellin moments. We then explore its role on the longitudinal correlation quantified with the ratio of the double distribution to the product of single distributions, $D(x_1,x_2, \vec{q}=\vec{0})/D(x_1)D(x_2)$. We point out that the ratios of moments $\langle x_1^n x_2^m \rangle / \langle x_1^n \rangle \langle x_2^m \rangle $ are independent of the evolution, providing particularly suitable measures to be tested in the upcoming lattice simulations. The transverse form factor $F(\vec{q})$ and its Fourier conjugate in the relative transverse coordinate $\vec{b}$ are obtained in variants of the Nambu--Jona-Lasinio model with the spectral and Pauli-Villars regularizations.The results are valid in the soft-momentum domain. Interestingly, with the spectral regularization of the model, the effective cross section for the double parton scattering of pions is exactly equal to the geometric cross section, $\sigma_{\rm eff}=\pi \langle \vec{b}^2 \rangle$ and yields about 20~mb.


title: Approximation of an algebra by evolution algebras

abstract: It is known that any multiplication of a finite dimensional algebra is determined by a matrix of structural constants. In general, this is a cubic matrix. Difficulty of investigation of an algebra depends on the cubic matrix. Such a cubic matrix defines a quadratic mapping called an evolution operator. In the case of evolution algebras, the cubic matrix consists many zeros allowing to reduce it to a square matrix. In this paper for any finite dimensional algebra we construct a family of evolution algebras corresponding to Jacobian of the evolution operator at a point of the algebra. We obtain some results answering how properties of an algebra depends on the properties of the corresponding family of evolution algebras. Moreover, we consider evolution algebras corresponding to 2 and 3-dimensional nilpotent Leibniz algebras. We prove that such evolution algebras are nilpotent too. Also we classify such evolution algebras.


title: Percentile-Based Residuals for Model Assessment

abstract: Residuals are a key component of diagnosing model fit. The usual practice is to compute standardized residuals using expected values and standard deviations of the observed data, then use these values to detect outliers and assess model fit. Approximate normality of these residuals is key for this process to have good properties, but in many modeling contexts, especially for complex, multi-level models, normality may not hold. In these cases outlier detection and model diagnostics aren't properly calibrated. Alternatively, as we demonstrate, residuals computed from the percentile location of a datum's value in its full predictive distribution lead to well calibrated evaluations of model fit. We generalize an approach described by Dunn and Smyth (1996) and evaluate properties mathematically, via case-studies and by simulation. In addition, we show that the standard residuals can be calibrated to mimic the percentile approach, but that this extra step is avoided by directly using percentile-based residuals. For both the percentile-based residuals and the calibrated standard residuals, the use of full predictive distributions with the appropriate location, spread and shape is necessary for valid assessments.


title: The Exceptional X-ray Evolution of SN1996cr in High Resolution

abstract: We present X-ray spectra spanning 18 years of evolution for SN\,1996cr, one of the five nearest SNe detected in the modern era. {\it Chandra} HETG exposures in 2000, 2004, and 2009 allow us to resolve spectrally the velocity profiles of Ne, Mg, Si, S, and Fe emission lines and monitor their evolution as tracers of the ejecta-circumstellar medium (CSM) interaction. To explain the diversity of X-ray line profiles, we explore several possible geometrical models. Based on the highest signal-to-noise 2009 epoch, we find that a polar geometry with two distinct opening angle configurations and internal obscuration can successfully reproduce all of the observed line profiles. The best fit model consists of two plasma components: (1) a mildly absorbed (2$\times$10$^{21}$\,cm$^{-2}$), cooler ($\approx$2\,keV) with high Ne, Mg, Si, and S abundances associated with a wide polar interaction region (half-opening angle $\approx$58$^{\circ}$); (2) a moderately absorbed (2$\times$10$^{22}$\,cm$^{-2}$), hotter (>~20\,keV) plasma with high Fe abundances and strong internal obscuration associated with a narrow polar interaction region (half-opening angle $\approx$20$^{\circ}$). We extend this model to seven further epochs with lower signal-to-noise ratio and/or lower spectral-resolution between 2000-2018, yielding several interesting trends in absorption, flux, geometry and expansion velocity. We argue that the hotter and colder components are associated with reflected and forward shocks, respectively, at least at later epochs. We discuss the physical implications of our results and plausible explosion scenarios to understand the X-ray data of SN\,1996cr.


title: Breakthrough Listen Follow-up of the Random Transiter (EPIC 249706694/HD 139139) with the Green Bank Telescope

abstract: The star EPIC 249706694 (HD 139139) was found to exhibit 28 transit-like events over an 87 day period during the Kepler mission's K2 Campaign 15 (Rappaport et al. 2019). These events did not fall into an identifiable pattern and could not be explained by a multitude of transit scenarios explored by the authors. We conduct follow-up observations at C-band frequencies with the Green Bank Telescope as part of the ongoing Breakthrough Listen search for technosignatures. We search for narrow band signals above a signal-to-noise threshold of 10 and with Doppler drift rates within +-5 Hz/s. We detect no evidence of technosignatures from EPIC 249706694 and derive an upper limit for the EIRP (Equivalent Isotropic Radiated Power) of putative transmissions to be 10 TW.


title: Political Openness and Armed Conflict: Evidence from Local Councils in Colombia

abstract: In this paper, I empirically investigate how the openness of political institutions to diverse representation can impact conflict-related violence. By exploiting plausibly exogenous variations in the number of councillors in Colombian municipalities, I develop two sets of results. First, regression discontinuity estimates show that larger municipal councils have a considerably greater number of political parties with at least one elected representative. I interpret this result as evidence that larger municipal councils are more open to diverse political participation. The estimates also reveal that non-traditional parties are the main beneficiaries of this greater political openness. Second, regression discontinuity estimates show that political openness substantially decreases conflict-related violence, namely the killing of civilian non-combatants. By exploiting plausibly exogenous variations in local election results, I show that the lower level of political violence stems from greater participation by parties with close links to armed groups. Using data about the types of violence employed by these groups, and representation at higher levels of government, I argue that armed violence has decreased not because of power-sharing arrangements involving armed groups linked to the parties with more political representation, but rather because armed groups with less political power and visibility are deterred from initiating certain types of violence.


title: MelGAN-VC: Voice Conversion and Audio Style Transfer on arbitrarily long samples using Spectrograms

abstract: Traditional voice conversion methods rely on parallel recordings of multiple speakers pronouncing the same sentences. For real-world applications however, parallel data is rarely available. We propose MelGAN-VC, a voice conversion method that relies on non-parallel speech data and is able to convert audio signals of arbitrary length from a source voice to a target voice. We firstly compute spectrograms from waveform data and then perform a domain translation using a Generative Adversarial Network (GAN) architecture. An additional siamese network helps preserving speech information in the translation process, without sacrificing the ability to flexibly model the style of the target speaker. We test our framework with a dataset of clean speech recordings, as well as with a collection of noisy real-world speech examples. Finally, we apply the same method to perform music style transfer, translating arbitrarily long music samples from one genre to another, and showing that our framework is flexible and can be used for audio manipulation applications different from voice conversion.


title: Hilbert Space Fragmentation and Ashkin-Teller Criticality in Fluctuation Coupled Ising Models

abstract: We discuss the effects of exponential fragmentation of the Hilbert space on phase transitions in the context of coupled ferromagnetic Ising models in arbitrary dimension with special emphasis on the one dimensional case. We show that the dynamics generated by quantum fluctuations is bounded within spatial partitions of the system and weak mixing of these partitions caused by global transverse fields leads to a zero temperature phase with ordering in the local product of both Ising copies but no long range order in either species. This leads to a natural connection with the Ashkin-Teller universality class for general lattices. We confirm this for the periodic chain using quantum Monte Carlo simulations. We also point out that our treatment provides an explanation for pseudo-first order behavior seen in the Binder cumulants of the classical frustrated $J_1-J_2$ Ising model and the $q=4$ Potts model in 2D.


title: Stable Intersections of Conformal Cantor Sets

abstract: We investigate stable intersections of conformal Cantor sets and their consequences to dynamical systems. First we define this type of Cantor set and relate it to horseshoes appearing in automorphisms of $\C^2$. Then we study limit geometries, objects related to the asymptotic shape of the Cantor sets, to obtain a criterion that guarantees stable intersection between some configurations. Finally we show that the Buzzard construction of a Newhouse region on $Aut(\C^2)$ can be seem as a case of stable intersection of Cantor sets in our sense and give some (not optimal) estimative on how \say{thick} those sets have to be.


title: Proving Global Optimality of ACOPF Solutions

abstract: We present our latest contributions in terms of mathematical modeling and algorithm development for the global optimization of the ACOPF problem. These contributions allow us to close the optimality gap on a number of open instances in the PGLIB and NESTA benchmark libraries. This is achieved by combining valid cut generation with semidefinite programming-based bound tightening. The mathematical formulations along with the solution algorithms are implemented in the modeling framework Gravity (www.gravityopt.com), an open-source platform for reproducible numerical experiments.


title: Schr\"odinger equations with smooth measure potential and general measure data

abstract: We study equations driven by Schr\"odinger operators consisting of a self-adjoint Dirichlet operator and a singular potential, which belongs to a class of positive Borel measures absolutely continuous with respect to a capacity generated by the operator. In particular, we cover positive potentials exploding on a set of capacity zero. The right-hand side of equations is allowed to be a general bounded Borel measure. The class of self-adjoint Dirichlet operators is quite large. Examples include integro-differential operators with the local part of divergence form. We give a necessary and sufficient condition for the existence of a solution, and prove some regularity and stability results.


title: On Dimension-free Tail Inequalities for Sums of Random Matrices and Applications

abstract: In this paper, we present a new framework to obtain tail inequalities for sums of random matrices. Compared with existing works, our tail inequalities have the following characteristics: 1) high feasibility--they can be used to study the tail behavior of various matrix functions, e.g., arbitrary matrix norms, the absolute value of the sum of the sum of the $j$ largest singular values (resp. eigenvalues) of complex matrices (resp. Hermitian matrices); and 2) independence of matrix dimension --- they do not have the matrix-dimension term as a product factor, and thus are suitable to the scenario of high-dimensional or infinite-dimensional random matrices. The price we pay to obtain these advantages is that the convergence rate of the resulting inequalities will become slow when the number of summand random matrices is large. We also develop the tail inequalities for matrix random series and matrix martingale difference sequence. We also demonstrate usefulness of our tail bounds in several fields. In compressed sensing, we employ the resulted tail inequalities to achieve a proof of the restricted isometry property when the measurement matrix is the sum of random matrices without any assumption on the distributions of matrix entries. In probability theory, we derive a new upper bound to the supreme of stochastic processes. In machine learning, we prove new expectation bounds of sums of random matrices matrix and obtain matrix approximation schemes via random sampling. In quantum information, we show a new analysis relating to the fractional cover number of quantum hypergraphs. In theoretical computer science, we obtain randomness-efficient samplers using matrix expander graphs that can be efficiently implemented in time without dependence on matrix dimensions.

categories: cs.LG math-ph math.MP math.ST quant-ph stat.ML stat.TH


title: Modeling of Dynamical Systems via Successive Graph Approximations

abstract: In this work, we propose a non-parametric technique for online modeling of systems with unknown nonlinear Lipschitz dynamics. The key idea is to successively utilize measurements to approximate the graph of the state-update function using envelopes described by quadratic constraints. The proposed approach is then demonstrated on two control applications: (i) computation of tractable bounds for unmodeled dynamics, and (ii) computation of positive invariant sets. We further highlight the efficacy of the proposed approach via a detailed numerical example.


title: Optimal L-Infinity Frequency Control in Microgrids Considering Actuator Saturation

abstract: Inverter-connected resources can improve transient stability in low-inertia grids by injecting active power to minimize system frequency deviations following disturbances. In practice, most generation and load disturbances are step changes and the engineering figure-of-merit is often the peak overshoot in frequency resulting from these step disturbances. In addition, the inverter-connected resources tend to saturate much more easily than conventional synchronous machines. However, despite these challenges, standard controller designs must deal with averaged quantities through $H_2$ or $H_\infty$ norms and must account for saturation in ad hoc manners. In this paper, we address these challenges by explicitly considering $L_\infty$ control with saturation using a linear matrix inequality-based approach. We show that this approach leads to significant improvements in stability performance.


title: News from Cosmic Ray Air Showers (ICRC 2019 -- Cosmic Ray Indirect Rapport)

abstract: This rapport summarizes the cosmic-ray indirect (CRI) session of the 36th ICRC conference. Updated measurements from several air-shower arrays with higher precision lead to the discovery of new features in the energy spectrum: HAWC measures a softening of the light component (p+He) around $10^{13.5}\,$eV; measurements of the Pierre Auger Observatory show that the second knee is a smooth feature extending at least over the range of $100-200\,$PeV and that the energy spectrum between the ankle and the cut-off cannot be described by a simple broken power law. Measurements of the mass composition confirm that the composition is a varying mixture of protons and nuclei at least up to several $10\,$EeV. A joint effort of several collaborations helps to better assess deficiencies of hadronic interaction models, e.g., by quantifying the muon deficit in the models over the shower energy. Anisotropy measurements with higher precision generally confirm earlier results. A change of the amplitude and phase of the equatorial dipole indicates that in the energy range between the second knee and the ankle there likely is a transition from Galactic to extragalactic sources. However, neither the most energetic Galatic nor the extragalactic sources have been discovered, which remains a primary science goal. Next to more exposure, an increase of measurement accuracy and decrease of systematic uncertainties will provide future progress. It is particularly exciting that new experiments are built and existing experiments upgraded to increase the measurement accuracy of the energy and mass composition, e.g., by combining radio antennas with particle detectors. Last but not least, there is a trend that experiments are designed such that they can target cosmic rays, photons, and neutrinos at the same time, which will facilitate multi-messenger astrophysics at the highest energies.


title: Computational Bounds for Doing Harmonic Analysis on Permutation Modules of Finite Groups

abstract: We develop an approach to finding upper bounds for the number of arithmetic operations necessary for doing harmonic analysis on permutation modules of finite groups. The approach takes advantage of the intrinsic orbital structure of permutation modules, and it uses the multiplicities of irreducible submodules within individual orbital spaces to express the resulting computational bounds. We conclude by showing that these bounds are surprisingly small when dealing with certain permutation modules arising from the action of the symmetric group on tabloids.


title: Knowledge Distillation from Internal Representations

abstract: Knowledge distillation is typically conducted by training a small model (the student) to mimic a large and cumbersome model (the teacher). The idea is to compress the knowledge from the teacher by using its output probabilities as soft-labels to optimize the student. However, when the teacher is considerably large, there is no guarantee that the internal knowledge of the teacher will be transferred into the student; even if the student closely matches the soft-labels, its internal representations may be considerably different. This internal mismatch can undermine the generalization capabilities originally intended to be transferred from the teacher to the student. In this paper, we propose to distill the internal representations of a large model such as BERT into a simplified version of it. We formulate two ways to distill such representations and various algorithms to conduct the distillation. We experiment with datasets from the GLUE benchmark and consistently show that adding knowledge distillation from internal representations is a more powerful method than only using soft-label distillation.


title: Tight estimates of exit and containment probabilities for nonlinear stochastic systems

abstract: Tight estimates of exit/containment probabilities are of particular importance in many control problems. Yet, estimating the exit/containment probabilities is non-trivial: even for linear systems (Ornstein-Uhlenbeck processes), the containment probability can be computed exactly for only some particular values of the system parameters. In this paper, we derive tight bounds on the containment probability for a class of nonlinear stochastic systems. The core idea is to compare the "pull strength" (how hard the deterministic part of the system dynamics pulls towards the origin) experienced by the nonlinear system at hand with that of a well-chosen process for which tight estimates of the containment probability are known or can be numerically obtained (e.g. an Ornstein-Uhlenbeck process). Specifically, the main technical contribution of this paper is to define a suitable dominance relationship between the pull strengths of two systems and to prove that this dominance relationship implies an order relationship between their containment probabilities. We also discuss the link with contraction theory and highlight some examples of applications.


title: Fast approximate simulation of finite long-range spin systems

abstract: Tau leaping is a popular method for performing fast approximate simulation of certain continuous time Markov chain models typically found in chemistry and biochemistry. This method is known to perform well when the transition rates satisfy some form of scaling behaviour. In a similar spirit to tau leaping, we propose a new method for approximate simulation of spin systems which approximates the evolution of spin at each site between sampling epochs as an independent two-state Markov chain. When combined with fast summation methods, our method offers considerable improvement in speed over the standard Doob-Gillespie algorithm. We provide a detailed analysis of the error incurred for both the number of sites incorrectly labelled and for linear functions of the state.


title: Optimizing multigrid reduction-in-time (MGRIT) and Parareal coarse-grid operators for linear advection

abstract: Parallel-in-time methods, such as multigrid reduction-in-time (MGRIT) and Parareal, provide an attractive option for increasing concurrency when simulating time-dependent PDEs in modern high-performance computing environments. While these techniques have been very successful for parabolic equations, it has often been observed that their performance suffers dramatically when applied to advection-dominated problems or purely hyperbolic PDEs using standard rediscretization approaches on coarse grids. In this paper, we apply MGRIT or Parareal to the constant-coefficient linear advection equation, appealing to existing convergence theory to provide insight into the typically non-scalable or even divergent behavior of these solvers for this problem. To overcome these failings, we replace rediscretization on coarse grids with improved coarse-grid operators that are computed by applying optimization techniques to approximately minimize error estimates from the convergence theory. One of our main findings is that, in order to obtain fast convergence as for parabolic problems, coarse-grid operators should take into account the behavior of the hyperbolic problem by tracking the characteristic curves. Our approach is tested for schemes of various orders using explicit or implicit Runge-Kutta methods combined with upwind-finite-difference spatial discretizations. In all cases, we obtain scalable convergence in just a handful of iterations, with parallel tests also showing significant speed-ups over sequential time-stepping. Our insight of tracking characteristics on coarse grids provides a key idea for solving the long-standing problem of efficient parallel-in-time integration for hyperbolic PDEs.


title: Vector bundles with numerically flat reduction on rigid analytic varieties and $p$-adic local systems

abstract: We show how to functorially attach continuous $p$-adic representations of the profinite fundamental group to vector bundles with numerically flat reduction on a proper rigid analytic variety over $\mathbb{C}_p$. This generalizes results by Deninger and Werner for vector bundles on smooth algebraic varieties. Our approach uses fundamental results on the pro-\'etale site of a rigid analytic variety introduced by Scholze. This enables us to get rid of the smoothness condition and to work in the analytic category. Moreover, under some mild conditions, the functor we construct gives a full embedding of the category of vector bundles with numerically flat reduction into the category of continuous $\mathbb{C}_p$-representations. This provides new insights into the $p$-adic Simpson correspondence in the case of a vanishing Higgs-field.


title: Investigation on the generalization of the Sampled Policy Gradient algorithm

abstract: The Sampled Policy Gradient (SPG) algorithm is a new offline actor-critic variant that samples in the action space to approximate the policy gradient. It does so by using the critic to evaluate the sampled actions. SPG offers theoretical promise over similar algorithms such as DPG as it searches the action-Q-value space independently of the local gradient, enabling it to avoid local minima. This paper aims to compare SPG to two similar actor-critic algorithms, CACLA and DPG. The comparison is made across two different environments, two different network architectures, as well as training on on-policy transitions in contrast to using an experience buffer. Results seem to show that although SPG does often not perform the worst, it doesn't always match the performance of the best performing algorithm at a particular task. Further experiments are required to get a better estimate of the qualities of SPG.


title: Large-scale Gastric Cancer Screening and Localization Using Multi-task Deep Neural Network

abstract: Gastric cancer is one of the most common cancers, which ranks third among the leading causes of cancer death. Biopsy of gastric mucosa is a standard procedure in gastric cancer screening test. However, manual pathological inspection is labor-intensive and time-consuming. Besides, it is challenging for an automated algorithm to locate the small lesion regions in the gigapixel whole-slide image and make the decision correctly.To tackle these issues, we collected large-scale whole-slide image dataset with detailed lesion region annotation and designed a whole-slide image analyzing framework consisting of 3 networks which could not only determine the screening result but also present the suspicious areas to the pathologist for reference. Experiments demonstrated that our proposed framework achieves sensitivity of 97.05% and specificity of 92.72% in screening task and Dice coefficient of 0.8331 in segmentation task. Furthermore, we tested our best model in real-world scenario on 10,315 whole-slide images collected from 4 medical centers.


title: Interaction Driven Floquet Engineering of Topological Superconductivity in Rashba Nanowires

abstract: We analyze, analytically and numerically, a periodically driven Rashba nanowire proximity coupled to an $s$-wave superconductor using bosonization and renormalization group analysis in the regime of strong electron-electron interactions. Due to the repulsive interactions, the superconducting gap is suppressed, whereas the Floquet Zeeman gap is enhanced, resulting in a higher effective value of $g$-factor compared to the non-interacting case. The flow equations for different coupling constants, velocities, and Luttinger-liquid parameters explicitly establish that even for small initial values of the Floquet Zeeman gap compared to the superconducting proximity gap, the interactions drive the system into the topological phase and the interband interaction term helps to achieve larger regions of the topological phase in parameter space.


title: ExpertMatcher: Automating ML Model Selection for Clients using Hidden Representations

abstract: Recently, there has been the development of Split Learning, a framework for distributed computation where model components are split between the client and server (Vepakomma et al., 2018b). As Split Learning scales to include many different model components, there needs to be a method of matching client-side model components with the best server-side model components. A solution to this problem was introduced in the ExpertMatcher (Sharma et al., 2019) framework, which uses autoencoders to match raw data to models. In this work, we propose an extension of ExpertMatcher, where matching can be performed without the need to share the client's raw data representation. The technique is applicable to situations where there are local clients and centralized expert ML models, but the sharing of raw data is constrained.


title: Ctrl-Z: Recovering from Instability in Reinforcement Learning

abstract: When learning behavior, training data is often generated by the learner itself; this can result in unstable training dynamics, and this problem has particularly important applications in safety-sensitive real-world control tasks such as robotics. In this work, we propose a principled and model-agnostic approach to mitigate the issue of unstable learning dynamics by maintaining a history of a reinforcement learning agent over the course of training, and reverting to the parameters of a previous agent whenever performance significantly decreases. We develop techniques for evaluating this performance through statistical hypothesis testing of continued improvement, and evaluate them on a standard suite of challenging benchmark tasks involving continuous control of simulated robots. We show improvements over state-of-the-art reinforcement learning algorithms in performance and robustness to hyperparameters, outperforming DDPG in 5 out of 6 evaluation environments and showing no decrease in performance with TD3, which is known to be relatively stable. In this way, our approach takes an important step towards increasing data efficiency and stability in training for real-world robotic applications.


title: Refinements of the trace formula for GL(2)

abstract: We express the discrete noncuspidal terms in the spectral side of the trace formula for GL(2) in terms of orbital integrals, obtaining a geometric expansion for the cuspidal part of the trace formula. Assuming the Ramanujan conjecture for GL(2), this is equal to the tempered part of the trace formula, providing an alternate approach to the method of Frenkel-Langlands-Ngo in the context of Beyond Endoscopy. Using this, we establish a formula which in principle leads to an $r$-trace formula, and conclude with some remarks regarding the primitivization of the trace formula.


title: Comparative study of panel and panelless-based reflectance conversion techniques for agricultural remote sensing

abstract: Small unmanned aircraft systems (sUAS) have allowed for thousands of aerial images to be captured at a moments notice. The simplicity and relative low cost of flying a sUAS has provided remote sensing practitioners, both commercial and academic, with a viable alternative to traditional remote sensing platforms (airplanes and satellites). This paper is an expanded follow-up study to an initial work. Three radiance-to-reflectance conversion methods were tested to determine the optimal technique to use for converting raw digital count imagery to reflectance maps. The first two methods utilized in-scene reflectance conversion panels along with the empirical line method (ELM), while the final method used an upward looking sensor that recorded the band-effective spectral downwelling irradiance. The methods employed were 2-Point ELM, 1-point ELM, and At-altitude Radiance Ratio (AARR). The average signed reflectance factor errors produced by these methods on real sUAS imagery were: -0.005, -0.0028, and -0.0244 respectively. These errors were produced across four altitudes (150, 225, 300 and 375ft), six targets (grass, asphalt, concrete, blue felt, green felt and red felt), five spectral bands (blue, green, red, red edge and near infrared), and three weather conditions (cloudy, partly cloudy and sunny). Finally, data was simulated using the MODTRAN code to generate downwelling irradiance and sensor reaching radiance to compute the theoretical results of the AARR technique. A multitude of variables were varied for these simulations (atmosphere, time, day, target, sensor height, and visibility), which resulted in an overall theoretically achievable signed reflectance factor error of 0.0023.


title: Reversal symmetries for cyclic paths away from thermodynamic equilibrium

abstract: If a system is at thermodynamic equilibrium, an observer cannot tell whether a film of it is being played forward or in reverse: any transition will occur with the same frequency in the forward as in the reverse direction. However, if expenditure of energy changes the rate of even a single transition to yield a non-equilibrium steady state, such time-reversal symmetry undergoes a widespread breakdown, far beyond the point at which the energy is expended. An explosion of interdependency also arises, with steady-state probabilities of system states depending in a complicated manner on the rate of every transition in the system. Nevertheless, in the midst of this global non-equilibrium complexity, we find that cyclic paths have reversibility properties that remain local, and which can exhibit symmetry, no matter how far the system is from thermodynamic equilibrium. Specifically, given any cycle of reversible transitions, the ratio of the frequencies with which the cycle is traversed in one direction versus the other is determined, in the long-time limit, only by the thermodynamic force on the cycle itself, without requiring knowledge of transition rates elsewhere in the system. In particular, if there is no net energy expenditure on the cycle, then, over long times, the cycle traversal frequencies are the same in either direction.


title: BoolSi: a tool for distributed simulations and analysis of Boolean networks

abstract: We present BoolSi, an open-source cross-platform command line tool for distributed simulations of deterministic Boolean networks with synchronous update. It uses MPI standard to support execution on computational clusters, as well as parallel processing on a single computer. BoolSi can be used to model the behavior of complex dynamic networks, such as gene regulatory networks. In particular, it allows for identification and statistical analysis of network attractors. We perform a case study of the activity of a cambium cell to demonstrate the capabilities of the tool.


title: Lone Pair Rotational Dynamics in Solids

abstract: Traditional classifications of crystalline phases focus on nuclear degrees of freedom. Through examination of both electronic and nuclear structure, we introduce the concept of an electronic plastic crystal. Such a material is classified by crystalline nuclear structure, while localized electronic degrees of freedom - here lone pairs - exhibit orientational motion at finite temperatures. This orientational motion is an emergent phenomenon arising from the coupling between electronic structure and polarization fluctuations generated by collective motions, such as phonons. Using ab initio molecular dynamics simulations, we predict the existence of electronic plastic crystal motion in halogen crystals and halide perovskites, and suggest that such motion may be found in a broad range of solids with lone pair electrons. Such fluctuations in the charge density should be observable, in principle via synchrotron scattering.


title: Magnon Blockade in a Hybrid Ferromagnet-Superconductor Quantum System

abstract: The implementation of a single magnon level quantum manipulation is one of the fundamental targets in quantum magnetism with a significant practical relevance for precision metrology, quantum information processing, and quantum simulation. Here, we demonstrate theoretically the feasibility of using a hybrid ferromagnet-superconductor quantum system to prepare a single magnon source based on magnon blockade effects. By numerically solving the quantum master equation, we show that the second-order correlation function of the magnon mode depends crucially on the relation between the qubit-magnon coupling strength and the driving detuning, and simultaneously signatures of the magnon blockade appear only under quite stringent conditions of a cryogenic temperature. In addition to providing perception into the quantum phenomena of magnon, the study of magnon blockade effects will help to develop novel technologies for exploring the undiscovered magnon traits at the quantum level and may find applications in designing single magnon emitters.


title: Towards De-identification of Legal Texts

abstract: In many countries, personal information that can be published or shared between organizations is regulated and, therefore, documents must undergo a process of de-identification to eliminate or obfuscate confidential data. Our work focuses on the de-identification of legal texts, where the goal is to hide the names of the actors involved in a lawsuit without losing the sense of the story. We present a first evaluation on our corpus of NLP tools in tasks such as segmentation, tokenization and recognition of named entities, and we analyze several evaluation measures for our de-identification task. Results are meager: 84% of the documents have at least one name not covered by NER tools, something that might lead to the re-identification of involved names. We conclude that tools must be strongly adapted for processing texts of this particular domain.


title: The Resolution of Keller's Conjecture

abstract: We consider three graphs, $G_{7,3}$, $G_{7,4}$, and $G_{7,6}$, related to Keller's conjecture in dimension 7. The conjecture is false for this dimension if and only if at least one of the graphs contains a clique of size $2^7 = 128$. We present an automated method to solve this conjecture by encoding the existence of such a clique as a propositional formula. We apply satisfiability solving combined with symmetry-breaking techniques to determine that no such clique exists. This result implies that every unit cube tiling of $\mathbb{R}^7$ contains a facesharing pair of cubes. Since a faceshare-free unit cube tiling of $\mathbb{R}^8$ exists (which we also verify), this completely resolves Keller's conjecture.


title: Multiple-objective Reinforcement Learning for Inverse Design and Identification

abstract: The aim of the inverse chemical design is to develop new molecules with given optimized molecular properties or objectives. Recently, generative deep learning (DL) networks are considered as the state-of-the-art in inverse chemical design and have achieved early success in generating molecular structures with desired properties in the pharmaceutical and material chemistry fields. However, satisfying a large number (larger than 10 objectives) of molecular objectives is a limitation of current generative models. To improve the model's ability to handle a large number of molecule design objectives, we developed a Reinforcement Learning (RL) based generative framework to optimize chemical molecule generation. Our use of Curriculum Learning (CL) to fine-tune the pre-trained generative network allowed the model to satisfy up to 21 objectives and increase the generative network's robustness. The experiments show that the proposed multiple-objective RL-based generative model can correctly identify unknown molecules with an 83 to 100 percent success rate, compared to the baseline approach of 0 percent. Additionally, this proposed generative model is not limited to just chemistry research challenges; we anticipate that problems that utilize RL with multiple-objectives will benefit from this framework.


title: Greedy Convex Ensemble

abstract: We consider learning a convex combination of basis models, and present some new theoretical and empirical results that demonstrate the effectiveness of a greedy approach. Theoretically, we first consider whether we can use linear, instead of convex, combinations, and obtain generalization results similar to existing ones for learning from a convex hull. We obtain a negative result that even the linear hull of very simple basis functions can have unbounded capacity, and is thus prone to overfitting; on the other hand, convex hulls are still rich but have bounded capacities. Secondly, we obtain a generalization bound for a general class of Lipschitz loss functions. Empirically, we first discuss how a convex combination can be greedily learned with early stopping, and how a convex combination can be non-greedily learned when the number of basis models is known a priori. Our experiments suggest that the greedy scheme is competitive with or better than several baselines, including boosting and random forests. The greedy algorithm requires little effort in hyper-parameter tuning, and also seems able to adapt to the underlying complexity of the problem. Our code is available at https://github.com/tan1889/gce.


title: Model-based Reinforcement Learning for Predictions and Control for Limit Order Books

abstract: We build a profitable electronic trading agent with Reinforcement Learning that places buy and sell orders in the stock market. An environment model is built only with historical observational data, and the RL agent learns the trading policy by interacting with the environment model instead of with the real-market to minimize the risk and potential monetary loss. Trained in unsupervised and self-supervised fashion, our environment model learned a temporal and causal representation of the market in latent space through deep neural networks. We demonstrate that the trading policy trained entirely within the environment model can be transferred back into the real market and maintain its profitability. We believe that this environment model can serve as a robust simulator that predicts market movement as well as trade impact for further studies.


title: Growth of massive black holes at high-z via accretion predominantly driven by magnetic outflows

abstract: Luminous quasars powered by accreting supermassive black holes (SMBHs) have been found in the early Universe at $z \gtrsim 7.5$, which set a strong constraint on both the seed black hole mass and the rapid growth of the SMBHs. In this work, we explore how the SMBHs are grown through Eddington limited accretion driven predominantly by magnetic outflows. Most angular momentum and the released gravitational energy in the disk can be removed by magnetic outflows, and therefore the mass accretion rate of the black hole (BH) can be high even if the disk is radiating at sub-Eddington luminosity. It is found that the SMBH with several billion solar masses discovered at $z\gtrsim 7$ may probably be grown through chaotic accretion predominantly driven by magnetic outflows from a stellar mass BH, when the disks are radiating at moderate luminosity ($\sim 0.5$ Eddington luminosity) with mild outflows. We find that most SMBHs are spinning at moderate values of spin parameter $a_*$, which implies only a small fraction of quasars may have radio jets.


title: On Odd Rainbow Cycles in Edge-Colored Graphs

abstract: Let $G = (V, E)$ be an $n$-vertex edge-colored graph. In 2013, H. Li proved that if every vertex $v \in V$ is incident to at least $(n+1)/2$ distinctly colored edges, then $G$ admits a rainbow triangle. We prove that the same hypothesis ensures a rainbow $\ell$-cycle $C_{\ell}$ whenever $n \ge 432 \ell$. This result is sharp for all odd integers $\ell \geq 3$, and extends earlier work of the authors for when $\ell$ is even.


title: A cascaded dual-domain deep learning reconstruction method for sparsely spaced multidetector helical CT

abstract: Helical CT has been widely used in clinical diagnosis. Sparsely spaced multidetector in z direction can increase the coverage of the detector provided limited detector rows. It can speed up volumetric CT scan, lower the radiation dose and reduce motion artifacts. However, it leads to insufficient data for reconstruction. That means reconstructions from general analytical methods will have severe artifacts. Iterative reconstruction methods might be able to deal with this situation but with the cost of huge computational load. In this work, we propose a cascaded dual-domain deep learning method that completes both data transformation in projection domain and error reduction in image domain. First, a convolutional neural network (CNN) in projection domain is constructed to estimate missing helical projection data and converting helical projection data to 2D fan-beam projection data. This step is to suppress helical artifacts and reduce the following computational cost. Then, an analytical linear operator is followed to transfer the data from projection domain to image domain. Finally, an image domain CNN is added to improve image quality further. These three steps work as an entirety and can be trained end to end. The overall network is trained using a simulated lung CT dataset with Poisson noise from 25 patients. We evaluate the trained network on another three patients and obtain very encouraging results with both visual examination and quantitative comparison. The resulting RRMSE is 6.56% and the SSIM is 99.60%. In addition, we test the trained network on the lung CT dataset with different noise level and a new dental CT dataset to demonstrate the generalization and robustness of our method.


title: The Daunting Task of Real-World Textual Style Transfer Auto-Evaluation

abstract: The difficulty of textual style transfer lies in the lack of parallel corpora. Numerous advances have been proposed for the unsupervised generation. However, significant problems remain with the auto-evaluation of style transfer tasks. Based on the summary of Pang and Gimpel (2018) and Mir et al. (2019), style transfer evaluations rely on three criteria: style accuracy of transferred sentences, content similarity between original and transferred sentences, and fluency of transferred sentences. We elucidate the problematic current state of style transfer research. Given that current tasks do not represent real use cases of style transfer, current auto-evaluation approach is flawed. This discussion aims to bring researchers to think about the future of style transfer and style transfer evaluation research.


title: Observation of anomalous Ettingshausen effect and large transverse thermoelectric conductivity in permanent magnets

abstract: This study focuses on the potential of permanent magnets as thermoelectric converters. It is found that a SmCo$_5$-type magnet exhibits the large anomalous Ettingshausen effect (AEE) at room temperature and that its charge-to-heat current conversion coefficient is more than one order of magnitude greater than that of typical ferromagnetic metals. The large AEE is an exclusive feature of the SmCo$_5$-type magnet among various permanent magnets in practical use, which is independent of the conventional performance of magnets based on static magnetic properties. The experimental results show that the large AEE originates from the intrinsic transverse thermoelectric conductivity of SmCo$_5$. This finding makes a connection between permanent magnets and thermal energy engineering, providing the basis for creating "thermoelectric permanent magnets."


title: The fastest $\ell_{1,\infty}$ prox in the west

abstract: Proximal operators are of particular interest in optimization problems dealing with non-smooth objectives because in many practical cases they lead to optimization algorithms whose updates can be computed in closed form or very efficiently. A well-known example is the proximal operator of the vector $\ell_1$ norm, which is given by the soft-thresholding operator. In this paper we study the proximal operator of the mixed $\ell_{1,\infty}$ matrix norm and show that it can be computed in closed form by applying the well-known soft-thresholding operator to each column of the matrix. However, unlike the vector $\ell_1$ norm case where the threshold is constant, in the mixed $\ell_{1,\infty}$ norm case each column of the matrix might require a different threshold and all thresholds depend on the given matrix. We propose a general iterative algorithm for computing these thresholds, as well as two efficient implementations that further exploit easy to compute lower bounds for the mixed norm of the optimal solution. Experiments on large-scale synthetic and real data indicate that the proposed methods can be orders of magnitude faster than state-of-the-art methods.


title: Aegis: A Context-aware Security Framework for Smart Home Systems

abstract: Our everyday lives are expanding fast with the introduction of new Smart Home Systems (SHSs). Today, a myriad of SHS devices and applications are widely available to users and have already started to re-define our modern lives. Smart home users utilize the apps to control and automate such devices. Users can develop their own apps or easily download and install them from vendor-specific app markets. App-based SHSs offer many tangible benefits to our lives, but also unfold diverse security risks. Several attacks have already been reported for SHSs. However, current security solutions consider smart home devices and apps individually to detect malicious actions rather than the context of the SHS as a whole. The existing mechanisms cannot capture user activities and sensor-device-user interactions in a holistic fashion. To address these issues, in this paper, we introduce Aegis, a novel context-aware security framework to detect malicious behavior in a SHS. Specifically, Aegis observes the states of the connected smart home entities (sensors and devices) for different user activities and usage patterns in a SHS and builds a contextual model to differentiate between malicious and benign behavior. We evaluated the efficacy and performance of Aegis in multiple smart home settings (i.e., single bedroom, double bedroom, duplex) with real-life users performing day-to-day activities and real SHS devices. We also measured the performance of Aegis against five different malicious behaviors. Our detailed evaluation shows that Aegis can detect malicious behavior in SHS with high accuracy (over 95%) and secure the SHS regardless of the smart home layout, device configuration, installed apps, and enforced user policies. Finally, Aegis achieves minimum overhead in detecting malicious behavior in SHS, ensuring easy deployability in real-life smart environments.


title: An MDL-Based Classifier for Transactional Datasets with Application in Malware Detection

abstract: We design a classifier for transactional datasets with application in malware detection. We build the classifier based on the minimum description length (MDL) principle. This involves selecting a model that best compresses the training dataset for each class considering the MDL criterion. To select a model for a dataset, we first use clustering followed by closed frequent pattern mining to extract a subset of closed frequent patterns (CFPs). We show that this method acts as a pattern summarization method to avoid pattern explosion; this is done by giving priority to longer CFPs, and without requiring to extract all CFPs. We then use the MDL criterion to further summarize extracted patterns, and construct a code table of patterns. This code table is considered as the selected model for the compression of the dataset. We evaluate our classifier for the problem of static malware detection in portable executable (PE) files. We consider API calls of PE files as their distinguishing features. The presence-absence of API calls forms a transactional dataset. Using our proposed method, we construct two code tables, one for the benign training dataset, and one for the malware training dataset. Our dataset consists of 19696 benign, and 19696 malware samples, each a binary sequence of size 22761. We compare our classifier with deep neural networks providing us with the state-of-the-art performance. The comparison shows that our classifier performs very close to deep neural networks. We also discuss that our classifier is an interpretable classifier. This provides the motivation to use this type of classifiers where some degree of explanation is required as to why a sample is classified under one class rather than the other class.


title: Probability, valuations, hyperspace: Three monads on Top and the support as a morphism

abstract: We consider three monads on Top, the category of topological spaces, which formalize topological aspects of probability and possibility in categorical terms. The first one is the Hoare hyperspace monad H, which assigns to every space its space of closed subsets equipped with the lower Vietoris topology. The second is the monad V of continuous valuations, also known as the extended probabilistic powerdomain. We construct both monads in a unified way in terms of double dualization. This reveals a close analogy between them, and allows us to prove that the operation of taking the support of a continuous valuation is a morphism of monads from V to H. In particular, this implies that every H-algebra (topological complete semilattice) is also a V-algebra. Third, we show that V can be restricted to a submonad of tau-smooth probability measures on Top. By composing these two morphisms of monads, we obtain that taking the support of a tau-smooth probability measure is also a morphism of monads.


title: A contextual analysis of the early work of Andrzej Trautman and Ivor Robinson on equations of motion and gravitational radiation

abstract: In a series of papers published in the course of his dissertation work in the mid 1950's, Andrzej Trautman drew upon the slow motion approximation developed by his advisor Infeld, the general covariance based strong conservation laws enunciated by Bergmann and Goldberg, the Riemann tensor attributes explored by Goldberg and related geodesic deviation exploited by Pirani, the permissible metric discontinuities identified by Lichnerowicz, O'Brien and Synge, and finally Petrov's classification of vacuum spacetimes. With several significant additions he produced a comprehensive overview of the state of research in equations of motion and gravitational waves that was presented in a widely cited series of lectures at King's College, London, in 1958. Fundamental new contributions were the formulation of boundary conditions representing outgoing gravitational radiation the deduction of its Petrov type, a covariant expression for null wave fronts, and a derivation of the correct mass loss formula due to radiation emission. Ivor Robinson had already in 1956 developed a bi-vector based technique that had resulted in his rediscovery of exact plane gravitational wave solutions of Einstein's equations. He was the first to characterize shear-free null geodesic congruences. He and Trautman met in London in 1958, and there resulted a long-term collaboration whose initial fruits were the Robinson-Trautman metric, examples of which were exact spherical gravitational waves.


title: Differential graded Lie algebras and Leibniz algebra cohomology

abstract: In this note, we interpret Leibniz algebras as differential graded Lie algebras. Namely, we consider two functors from the category of Leibniz algebras to that of differential graded Lie algebras and show that they naturally give rise to the Leibniz cohomology and the Chevalley-Eilenberg cohomology. As an application, we prove a conjecture stated by Pirashvili in arXiv:1904.00121 [math.KT].


title: A model for the complex reflection coefficient of a collection of parallel layers

abstract: Reflectometry is a technique that uses the light reflected by a sample to determine properties of the sample. Interferometric reflectometry uses interference between two beams, one of which is incident on ---and reflected back by--- a sample and one of which is not, to obtain the complex electric field rather than merely its intensity. Since this interference allows one to retrieve an increased amount of information about the light, it also allows one to obtain more information about the sample, such as a thin layer. We will apply the methods derived here to the case of a planar lipid bilayer.


title: Alternating Recurrent Dialog Model with Large-scale Pre-trained Language Models

abstract: Existing dialog system models require extensive human annotations and are difficult to generalize to different tasks. The recent success of large pre-trained language models such as BERT and GPT-2 (Devlin et al., 2019; Radford et al., 2019) have suggested the effectiveness of incorporating language priors in down-stream NLP tasks. However, how much pre-trained language models can help dialog response generation is still under exploration. In this paper, we propose a simple, general, and effective framework: Alternating Roles Dialog Model (ARDM). ARDM models each speaker separately and takes advantage of the large pre-trained language model. It requires no supervision from human annotations such as belief states or dialog acts to achieve effective conversations. ARDM outperforms or is on par with state-of-the-art methods on two popular task-oriented dialog datasets: CamRest676 and MultiWOZ. Moreover, we can generalize ARDM to more challenging, non-collaborative tasks such as persuasion. In persuasion tasks, ARDM is capable of generating human-like responses to persuade people to donate to a charity.


title: Secret key agreement from correlated data, with no prior information

abstract: A fundamental question that has been studied in cryptography and in information theory is whether two parties can communicate confidentially using exclusively an open channel. We consider the model in which the two parties hold inputs that are correlated in a certain sense. This model has been studied extensively in information theory, and communication protocols have been designed which exploit the correlation to extract from the inputs a shared secret key. However, all the existing protocols are not universal in the sense that they require that the two parties also know some attributes of the correlation. In other words, they require that each party knows something about the other party's input. We present a protocol that does not require any prior additional information. It uses space-bounded Kolmogorov complexity to measure correlation and it allows the two legal parties to obtain a common key that looks random to an eavesdropper that observes the communication and is restricted to use a bounded amount of space for the attack. Thus the protocol achieves complexity-theoretical security, but it does not use any unproven result from computational complexity. On the negative side, the protocol is not efficient in the sense that the computation of the two legal parties uses more space than the space allowed to the adversary.


title: Desingularization of vortices for 2D steady Euler flows via the vorticity method

abstract: In this paper, we consider steady Euler flows in a planar bounded domain in which the vorticity is sharply concentrated in a finite number of disjoint regions of small diameter. Such flows are closely related to the point vortex model and can be regarded as desingularization of point vortices. By an adaption of the vorticity method, we construct a family of steady Euler flows in which the vorticity is concentrated near a global minimum point of the Robin function of the domain, and the corresponding stream function satisfies a semilinear elliptic equation with a given profile function. Furthermore, for any given isolated minimum point $(\bar{x}_1,\cdot\cdot\cdot,\bar{x}_k)$ of the Kirchhoff-Routh function of the domain, we prove that there exists a family of steady Euler flows whose vorticity is supported in $k$ small regions near $\bar{x}_i$, and near each $\bar{x}_i$ the corresponding stream function satisfies a semilinear elliptic equation with a given profile function.


title: An event-triggered transmission scheduling strategy for remote state estimation in the presence of an eavesdropper

abstract: We consider a remote state estimation problem in the presence of an eavesdropper over packet dropping links. A smart sensor transmits its local estimates to a legitimate remote estimator, in the course of which an eavesdropper can randomly overhear the transmission. This problem has been well studied for unstable dynamical systems, but seldom for stable systems. In this paper, we target at stable and marginally stable systems and aim to design an event-triggered scheduling strategy by minimizing the expected error covariance at the remote estimator and keeping that at the eavesdropper above a user-specified lower bound. To this end, we model the evolution of the error covariance as an infinite recurrent Markov chain and develop a recurrence relation to describe the stationary distribution of the state at the eavesdropper. Monotonicity and convergence properties of the expected error covariance are further investigated and employed to solve the optimization problem. Numerical examples are provided to validate the theoretical results.


title: Superconducting proximity effect in a Rashba-type surface state of Pb/Ge(111)

abstract: The Rashba superconductor, in which spin-splitting bands become superconducting, is fascinating as a novel superconducting system in low dimensional systems. Here, we present the results of $\textit{in-situ}$ transport measurements on a Rashba-type surface state of the striped incommensurate (SIC) phase of a Pb atomic layer on Ge(111) surface with additional Pb islands/clusters on it. We found that two-step superconducting transitions at around 7 K and 3 K occurred. The latter superconducting transition is suggested to be induced at the non-superconducting Rashba SIC area because of the lateral proximity effect caused by the superconducting Pb clusters. Our results propose a new type of Rashba superconductor, which is a new platform to understand the Rashba superconducting systems.


title: On the number of limit cycles for polycycles $S^{(2)}$ and $S^{(3)}$ in quadratic Hamilton systems under perturbations of piecewise smooth polynomials

abstract: In this paper, by using Picard-Fuchs equations and Chebyshev criterion, we study the bifurcate of limit cycles for quadratic Hamilton system $S^{(2)}$ and $S^{(3)}$: $\dot{x}= y+2axy+by^2$, $\dot{y}=-x+x^2-ay^2$ with $a\in(-\frac{1}{2},1)$, $b=(1-a)(1+2a)^{1/2}$ and $a=1$, $b=0$ respectively, under perturbations of piecewise smooth polynomials with degree $n$. The discontinuity is the line $y=0$. We bound the number of zeros of first order Melnikov function which controls the number of limit cycles bifurcating from the center. It is proved that the upper bounds of the number of limit cycles for $S^{(2)}$ and $S^{(3)}$ are respectively $25n+161$ $(n\geq3)$ and $24n+126$ $(n\geq3)$ (taking into account the multiplicity).


title: Photon-assisted entanglement and squeezing generation and decoherence suppression via a quadratic optomechanical coupling

abstract: Entanglement and quantum squeezing have wide applications in quantum technologies due to their non-classical characteristics. Here we study entanglement and quantum squeezing in an open spin-optomechanical system, in which a Rabi model (a spin coupled to the mechanical oscillator) is coupled to an ancillary cavity field via a quadratic optomechanical coupling. We find that their performances can be significantly modulated via the photon of the ancillary cavity, which comes from photon-dependent spin-oscillator coupling and detuning. Specifically, a fully switchable spin-oscillator entanglement can be achieved, meanwhile a strong mechanical squeezing is also realized. Moreover, we study the environment-induced decoherence and dissipation, and find that they can be mitigated by increasing the number of photons. This work provides an effective way to manipulate entanglement and quantum squeezing and to suppress decoherence in the cavity quantum electrodynamics with a quadratic optomechanics.


title: Effective Dynamics from Coherent State Path Integral of Full Loop Quantum Gravity

abstract: A new routine is proposed to relate Loop Quant Cosmology (LQC) to Loop Quantum Gravity (LQG) from the perspective of effective dynamics. We derive the big-bang singularity resolution and big bounce from the first principle of full canonical LQG. Our results are obtained in the framework of the reduced phase space quantization of LQG. As a key step in our work, we derive with coherent states a new discrete path integral formula of the transition amplitude generated by the physical Hamiltonian. The semiclassical approximation of the path integral formula gives an interesting set of effective equations of motion (EOMs) for full LQG. When solving the EOMs with homogeneous and isotropic ansatz, we reproduce the LQC effective dynamics in $\mu_0$-scheme. The solution replaces the big-bang singularity by a big bounce. In the end, we comment on the possible relation between the $\bar{\mu}$-scheme of effective dynamics and the continuum limit of the path integral formula.


title: Geck's Conjecture and the Generalized Gelfand-Graev Representations in Bad Characteristic

abstract: For a connected reductive algebraic group $G$ defined over a finite field $\mathbb F_q$, Kawanaka introduced the generalized Gelfand-Graev representations (GGGRs for short) of the finite group $G(\mathbb F_q)$ in the case where $q$ is a power of a good prime for $G$. This representation has been widely studied and used in various contexts. Recently, Geck proposed a conjecture, characterizing Lusztig's special unipotent classes in terms of weighted Dynkin diagrams. Based on this conjecture, he gave a guideline for extending the definition of GGGRs to the case where $q$ is a power of a bad prime for $G$. Here, we will give a proof of Geck's conjecture. Combined with Geck's pioneer work, our proof verifies Geck's conjectural characterization of special unipotent classes, and completes his definition of GGGRs in bad characteristics.


title: The reachable space of the heat equation for a finite rod as a Reproducing Kernel Hilbert Space

abstract: We use some results from the theory of Reproducing Kernel Hilbert Spaces to show that the reachable space of the heat equation for a finite rod with either one or two Dirichlet boundary controls is a RKHS of analytic functions on a square, and we compute its reproducing kernel. We also show that the null reachable space of the heat equation for the half line with Dirichlet boundary data is a RKHS of analytic functions on a sector, whose reproducing kernel is (essentially) the sum of pullbacks of the Bergman and Hardy kernels on the half plane $\mathbb{C}^+$. We also consider the case with Neumann boundary data.


title: A Nonparametric Bayesian Framework for Uncertainty Quantification in Stochastic Simulation

abstract: When we use simulation to assess the performance of stochastic systems, the input models used to drive simulation experiments are often estimated from finite real-world data. There exist both input model and simulation estimation uncertainties in the system performance estimates. Without strong prior information on the input models and the system mean response surface, in this paper, we propose a Bayesian nonparametric framework to quantify the impact from both sources of uncertainty. Specifically, since the real-world data often represent the variability caused by various latent sources of uncertainty, Dirichlet Processes Mixtures (DPM) based nonparametric input models are introduced to model a mixture of heterogeneous distributions, which can faithfully capture the important features of real-world data, such as multi-modality and skewness. Bayesian posteriors of flexible input models characterize the input model estimation uncertainty, which automatically accounts for both model selection and parameter value uncertainty. Then, input model estimation uncertainty is propagated to outputs by using direct simulation. Thus, under very general conditions, our framework delivers an empirical credible interval accounting for both input and simulation uncertainties. A variance decomposition is further developed to quantify the relative contributions from both sources of uncertainty. Our approach is supported by rigorous theoretical and empirical study.


title: Flat band in two-dimensional non-Hermitian optical lattices

abstract: We propose a method to generate a real-energy flat band in a two-dimensional (2D) non-Hermitian Lieb lattice. The coincidence of the flat band eigenstate in both real and momentum spaces is essential for the proposed flat band, which is flexible at the appropriate match between the synthetic magnetic flux and non-Hermiticity. The proposed method is not limited to the 2D non-Hermitian Lieb lattice, and is applied to the 2D non-Hermitian Tasaki's decorated square lattice, dice lattice, and kagome lattice. Our findings make a step forward for the flat band engineering in 2D non-Hermitian optical lattices.


title: Physiological Closed-Loop Control (PCLC) Systems: Review of a Modern Frontier in Automation

abstract: Over the past decade, there has been an unprecedented international focus on improved quality and availability of medical care, which has reignited interest in clinical automation and drawn researchers toward novel solutions in the field of physiological closed-loop control systems (PCLCs). Today, multidisciplinary groups of expert scientists, engineers, clinicians, mathematicians, and policy-makers are combining their knowledge and experience to develop both the next generation of PCLC-based medical equipment and a collaborative commercial/academic infrastructure to support this rapidly expanding frontier. In the following article, we provide a robust introduction to the various aspects of this growing field motivated by the recent and ongoing work supporting two leading technologies: the artificial pancreas (AP) and automated anesthesia. Following a brief high-level overview of the main concepts in automated therapy and some relevant tools from systems and control theory, we explore -- separately -- the developments, challenges, state-of-the-art, and probable directions for AP and automated anesthesia systems. We then close the review with a consideration of the common lessons gleaned from these ventures and the implications they present for future investigations and adjacent research.


title: Bifluxon: Fluxon-Parity-Protected Superconducting Qubit

abstract: We have developed and characterized a symmetry-protected superconducting qubit that offers simultaneous exponential suppression of energy decay from charge and flux noise, and dephasing from flux noise. The qubit consists of a Cooper-pair box (CPB) shunted by a superinductor, thus forming a superconducting loop. Provided the offset charge on the CPB island is an odd number of electrons, the qubit potential corresponds to that of a $\cos \phi / 2$ Josephson element, preserving the parity of fluxons in the loop via Aharonov-Casher interference. In this regime, the logical-state wavefunctions reside in disjoint regions of phase space, thereby ensuring the protection against energy decay. By switching the protection on, we observed a ten-fold increase of the decay time, reaching up to $100 \mu \mathrm{s}$. Though the qubit is sensitive to charge noise, the sensitivity is much reduced in comparison with the charge qubit, and the charge-noise-induced dephasing time of the current device exceeds $1 \mu \mathrm{s}$. Implementation of the full dephasing protection can be achieved in the next-generation devices by combining several $\cos \phi / 2$ Josephson elements in a small array.


title: Ozsvath-Szabo bordered algebras and subquotients of category O

abstract: We show that Ozsv\'ath-Szab\'o's bordered algebra used to efficiently compute knot Floer homology is a graded flat deformation of the regular block of a $\mathfrak{q}$-presentable quotient of parabolic category $\mathcal{O}$. We identify the endomorphism algebra of a minimal projective generator for this block with an explicit quotient of the Ozsv\'ath-Szab\'o algebra using Sartori's diagrammatic formulation of the endomorphism algebra. Both of these algebras give rise to categorifications of tensor products of the vector representation $V^{\otimes n}$ for $U_q(\mathfrak{gl}(1|1))$. Our isomorphism allows us to transport a number of constructions between these two algebras, leading to a new (fully) diagrammatic reinterpretation of Sartori's algebra, new modules over Ozsv\'ath-Szab\'o's algebra lifting various bases of $V^{\otimes n}$, and bimodules over Ozsv\'ath-Szab\'o's algebra categorifying the action of the quantum group element $F$ and its dual on $V^{\otimes n}$.


title: HuggingFace's Transformers: State-of-the-art Natural Language Processing

abstract: Recent progress in natural language processing has been driven by advances in both model architecture and model pretraining. Transformer architectures have facilitated building higher-capacity models and pretraining has made it possible to effectively utilize this capacity for a wide variety of tasks. \textit{Transformers} is an open-source library with the goal of opening up these advances to the wider machine learning community. The library consists of carefully engineered state-of-the art Transformer architectures under a unified API. Backing this library is a curated collection of pretrained models made by and available for the community. \textit{Transformers} is designed to be extensible by researchers, simple for practitioners, and fast and robust in industrial deployments. The library is available at \url{https://github.com/huggingface/transformers}.


title: Semi-parametric Bayes Regression with Network Valued Covariates

abstract: There is an increasing recognition of the role of brain networks as neuroimaging biomarkers in mental health and psychiatric studies. Our focus is posttraumatic stress disorder (PTSD), where the brain network interacts with environmental exposures in complex ways to drive the disease progression. Existing linear models seeking to characterize the relation between the clinical phenotype and the entire edge set in the brain network may be overly simplistic and often involve inflated number of parameters leading to computational burden and inaccurate estimation. In one of the first such efforts, we develop a novel two stage Bayesian framework to find a node-specific lower dimensional representation for the network using a latent scale approach in the first stage, and then use a flexible Gaussian process regression framework for prediction involving the latent scales and other supplementary covariates in the second stage. The proposed approach relaxes linearity assumptions, addresses the curse of dimensionality and is scalable to high dimensional networks while maintaining interpretability at the node level of the network. Extensive simulations and results from our motivating PTSD application show a distinct advantage of the proposed approach over competing linear and non-linear approaches in terms of prediction and coverage.


title: Parametrizing the Detector Response with Neural Networks

abstract: In high energy physics, characterizing the response of a detector to radiation is one of the most important and basic experimental tasks. In many cases, this task is accomplished by parameterizing summary statistics of the full detector response probability density. The parameterized detector response can then be used for calibration as well as for directly improving physics analysis sensitivity. This paper discusses how to parameterize summary statistics of the detector response using neural networks. In particular, neural networks are powerful tools for incorporating multidimensional data and the loss function used during training determines which summary statistic is learned. One common summary statistic that has not been combined with deep learning (as far as the authors are aware) is the mode. A neural network-based approach to mode learning is proposed and empirically demonstrated in the context of high energy jet calibrations. Altogether, the neural network-based toolkit for detector response parameterization can enhance the utility of data collected at high energy physics experiments and beyond.


title: Strong Morita equivalence for inclusions of $C^*$-algebras induced by twisted actions of a countable discrete group

abstract: We consider two twisted actions of a countable discrete group on $\sigma$-unital $C^*$-algebras. Then by taking the reduced crossed products, we get two inclusions of $C^*$-algebras. We suppose that they are strongly Morita equivalent as inclusions of $C^*$-algebras. Also, we suppose that one of the inclusions is irreducible, that is, the relative commutant of one of the $\sigma$-unital $C^*$-algebra in the multiplier $C^*$-algebra of the reduced twisted crossed product is trivial. We show that the two actions are then strongly Morita equivalent up to some automorphism of the group.


title: Approximate semi-amenability of Banach algebras

abstract: In recent work of the authors the notion of a derivation being approximately semi-inner arose as a tool for investigating (approximate) amenability questions for Banach algebras. Here we investigate this property in its own right, together with the consequent one of approximately semi-amenability. Under certain hypotheses regarding approximate identities this new notion is the same as approximate amenability, but more generally it covers some important classes of algebras which are not approximately amenable, in particular Segal algebras on amenable SIN-groups.


title: Absence of replica symmetry breaking in disordered FKG-Ising models under uniform field

abstract: We prove that the variance of spin overlap vanishes in disordered Ising models satisfying the Fortuin-Kasteleyn-Ginibre (FKG) inequality under a uniform field, such as generally distributed random field Ising model, site- and bond-diluted Ising models with the Bernoulli distribution. Chatterjee's proof for the Gaussian random field Ising model is generalized to other independent identically distributed quenched disorder under a uniform field.


title: Design of Phosphorene for Hydrogen Evolution Performance Comparable to Platinum

abstract: Phosphorene, a monolayer of bulk black phosphorus, is promising for light harvest owing to its high charge mobility and tunable direct band gap covering a broad spectral range of light. Here, via atomic-scale first-principles simulations, we report a ultrahigh activity of hydrogen evolution reaction (HER) of phosphorene originated from defective activation. Quantitative evaluation of the Gibbs free energy of the ad desorption of hydrogen (H) to from phosphorene (delta_GH) reveals that atomic vacancies and edges play a dominant role in activating the reaction. We find that the defective states, empty and well-localized around the defect core, are compensated by H species. This induces a proper hydrogen interaction complying with the thermoneutral condition of the free energy (delta_GH=0) comparable to platinum. Our findings of the highly activating defective states suggest the design of non-metal HER catalysts with structural engineering of earth-abundant phosphorus structures.


title: Maturation of siRNA by strand separation: Steered Molecular dynamics study

abstract: RNA interference, particularly siRNA induced gene silencing is becoming an important avenue of modern therapeutics. The siRNA is delivered to the cells as short double helical RNA which becomes single stranded for forming the RISC complex. Significant experimental evidence is available for most of the steps except the process of the separation of the two strands. We have attempted to understand the pathway for double stranded siRNA (dsRNA) to single stranded (ssRNA) molecules using steered molecular dynamics simulations. As the process is completely unexplored we have applied force from all possible directions restraining all possible residues to convert dsRNA to ssRNA. We found pulling one strand along helical axis direction restraining far end of the other strand demands excessive force for ssRNA formation. Pulling a central residue of one strand, in a direction perpendicular to the helix axis, while keeping the base paired residue fixed requires intermediate force for strand separation. Moreover we found that in this process the force requirement is quite high for the first bubble formation (nucleation energy) and the bubble propagation energies are quite small. We hypothesize this could be the mode of action adopted by the proteins in the cells.


title: Forecast Aggregation via Peer Prediction

abstract: Crowdsourcing enables the solicitation of forecasts on a variety of prediction tasks from distributed groups of people. How to aggregate the solicited forecasts, which may vary in quality, into an accurate final prediction remains a challenging yet critical question. Studies have found that weighing expert forecasts more in aggregation can improve the accuracy of the aggregated prediction. However, this approach usually requires access to the historical performance data of the forecasters, which are often not available. In this paper, we study the problem of aggregating forecasts without having historical performance data. We propose using peer prediction methods, a family of mechanisms initially designed to truthfully elicit private information in the absence of ground truth verification, to assess the expertise of forecasters, and then using this assessment to improve forecast aggregation. We evaluate our peer-prediction-aided aggregators on a diverse collection of 14 human forecast datasets. Compared with a variety of existing aggregators, our aggregators achieve a significant and consistent improvement on aggregation accuracy measured by the Brier score and the log score. Our results reveal the effectiveness of identifying experts to improve aggregation even without historical data.


title: Evolution of Intrinsic Vacancies and Prolonged Lifetime of Vacancy Clusters in Black Phosphorene

abstract: Due to the relatively low formation energies and highly mobile characteristics of atomic vacancies in phosphorene, understanding their evolutions becomes crucial for its structural integrity, chemical activities and applications. Herein, by combining first-principles calculations and kinetic Monte Carlo simulation, we investigate the time evolution and formation of atomic vacancy clusters from isolated monovacancies (MVs), aiming to uncover the mechanisms of diffusion, annihilation, reaction of these atomic vacancies. We find that while isolated MVs possess a highly mobile character, they react and form MV pairs which possess much lower mobility and high stability under ambient condition. We also show that the disappearance of MVs at the edge is quite slow due to the relatively high energy barrier, and as a result, around 80% of MVs remains even after two years under ambient condition. Our findings on one hand provide useful information for the structural repairing of phosphorene through chemical functionalization of these vacancy clusters, and on the other hand, suggest that these rather stable vacancy clusters may be used as an activated catalyst.


title: Learning Visual Affordances with Target-Orientated Deep Q-Network to Grasp Objects by Harnessing Environmental Fixtures

abstract: This paper introduces a challenging object grasping task and proposes a self-supervised learning approach. The goal of the task is to grasp an object which is not feasible with a single parallel gripper, but only with harnessing environment fixtures (e.g., walls, furniture, heavy objects). This Slide-to-Wall grasping task assumes no prior knowledge except the partial observation of a target object. Hence the robot should learn an effective policy given a scene observation that may include the target object, environmental fixtures, and any other disturbing objects. We formulate the problem as visual affordances learning for which Target-Oriented Deep Q-Network (TO-DQN) is proposed to efficiently learn visual affordance maps (i.e., Q-maps) to guide robot actions. Since the training necessitates robot's exploration and collision with the fixtures, TO-DQN is first trained safely with a simulated robot manipulator and then applied to a real robot. We empirically show that TO-DQN can learn to solve the task in different environment settings in simulation and outperforms a standard and a variant of Deep Q-Network (DQN) in terms of training efficiency and robustness. The testing performance in both simulation and real-robot experiments shows that the policy trained by TO-DQN achieves comparable performance to humans.


title: Growth of Ultra-high Purity NaI(Tl) Crystal for Dark Matter Searches

abstract: The annual modulation of scintillation event rate observed by the DAMA/LIBRA experiment has been a long-standing controversy in the quest of the direct detection of dark matter. The effort to definitively confirm or refute the annual modulation has turned out to be challenging due to the lack of NaI(Tl) crystals with high enough radio-purity. Most recently, we successfully grew a 6-kg ingot free from contamination during growth, from which a 3.4-kg crystal scintillator was made. The $^{39}$K concentration in the final crystal is estimated to be 4.3$\pm$0.2~ppb, unprecedented for NaI(Tl) crystals. The $^{210}$Pb activity is estimated to be 0.34$\pm$0.04~mBq/kg via $\alpha$ counting of $^{210}$Po, among the lowest of currently-running NaI-based dark matter experiments except DAMA/LIBRA. More importantly, the techniques and protocols we have developed will further contribute to the growth of higher purity NaI(Tl) crystals for dark matter searches.


title: Electric Load and Power Forecasting Using Ensemble Gaussian Process Regression

abstract: We propose a new forecasting method for predicting load demand and generation scheduling. Accurate week-long forecasting of load demand and optimal power generation is critical for efficient operation of power grid systems. In this work, we use a synthetic data set describing a power grid with 700 buses and 134 generators over a 365-days period with data synthetically generated at an hourly rate. The proposed approach for week-long forecasting is based on the Gaussian process regression (GPR) method, with prior covariance matrices of the quantities of interest (QoI) computed from ensembles formed by up to twenty preceding weeks of QoI observations. Then, we use these covariances within the GPR framework to forecast the QoIs for the following week. We demonstrate that the the proposed ensemble GPR (EGPR) method is capable of accurately forecasting weekly total load demand and power generation profiles. The EGPR method is shown to outperform traditional forecasting methods including the standard GPR and autoregressive integrated moving average (ARIMA) methods.


title: Generalized Property-Directed Reachability for Hybrid Systems

abstract: Generalized property-directed reachability (GPDR) belongs to the family of the model-checking techniques called IC3/PDR. It has been successfully applied to software verification; for example, it is the core of Spacer, a state-of-the-art Horn-clause solver bundled with Z3. However, it has yet to be applied to hybrid systems, which involve a continuous evolution of values over time. As the first step towards GPDR- based model checking for hybrid systems, this paper formalizes HGPDR, an adaptation of GPDR to hybrid systems, and proves its soundness. We also implemented a semi-automated proof-of-concept verifier, which allows a user to provide hints to guide verification steps.


title: Direct coupling of ferromagnetic moment and ferroelectric polarization in BiFeO$_3$

abstract: The spin-driven component of electric polarization in a single crystal of multiferroic BiFeO$_{3}$ was experimentally investigated in pulsed high magnetic fields up to 41 T. Sequential measurements of electric polarization for various magnetic field directions provide clear evidence of electric polarization normal to the hexagonal $c$ axis (${\bf P}_{\rm t}$) in not only the cycloidal phase, but also the field-induced canted antiferromagnetic phase. The direction of ${\bf P}_{\rm t}$ is directly coupled with the ferromagnetic moment in the canted antiferromagnetic phase, and thus controlled by changing the direction of the applied magnetic field. This magnetoelectric coupling is reasonably reproduced by the metal-ligand hybridization model.


title: Global Convergence for Replicator Dynamics of Repeated Snowdrift Games

abstract: To understand the emergence and sustainment of cooperative behavior in interacting collectives, we perform global convergence analysis for replicator dynamics of a large, well-mixed population of individuals playing a repeated snowdrift game with four typical strategies, which are always cooperate (ALLC), tit-for-tat (TFT), suspicious tit-for-tat (STFT) and always defect (ALLD). The dynamical model is a three-dimensional ODE system that is parameterized by the payoffs of the base game. Instead of routine searches for evolutionarily stable strategies and sets, we expand our analysis to determining the asymptotic behavior of solution trajectories starting from any initial state, and in particular show that for the full range of payoffs, every trajectory of the system converges to an equilibrium point. The convergence results highlight three findings that are of particular importance for understanding the cooperation mechanisms among self-interested agents playing repeated snowdrift games. First, the inclusion of TFT- and STFT-players, the two types of conditional strategy players in the game, increases the share of cooperators of the overall population compared to the situation when the population consists of only ALLC- and ALLD-players. This confirms findings in biology and sociology that reciprocity may promote cooperation in social collective actions, such as reducing traffic jams and division of labor, where each individual may gain more to play the opposite of what her opponent chooses. Second, surprisingly enough, regardless of the payoffs, there always exists a set of initial conditions under which ALLC players do not vanish in the long run, which does not hold for all the other three types of players. So an ALLC-player, although perceived as the one that can be easily taken advantage of in snowdrift games, has certain endurance in the long run.


title: Supervised feature selection with orthogonal regression and feature weighting

abstract: Effective features can improve the performance of a model, which can thus help us understand the characteristics and underlying structure of complex data. Previous feature selection methods usually cannot keep more local structure information. To address the defects previously mentioned, we propose a novel supervised orthogonal least square regression model with feature weighting for feature selection. The optimization problem of the objection function can be solved by employing generalized power iteration (GPI) and augmented Lagrangian multiplier (ALM) methods. Experimental results show that the proposed method can more effectively reduce the feature dimensionality and obtain better classification results than traditional feature selection methods. The convergence of our iterative method is proved as well. Consequently, the effectiveness and superiority of the proposed method are verified both theoretically and experimentally.