title: The residual finiteness of (hyperbolic) automorphism-induced HNN-extensions

abstract: We classify finitely generated, residually finite automorphism-induced HNN-extensions in terms of the residual separability of a single associated subgroup. This classification provides a method to construct automorphism-induced HNN-extensions which are not residually finite. We prove that this method can never yield a "new" counter-example to Gromov's conjecture on the residual finiteness of hyperbolic groups.


title: Ultrathin and Multicolor Optical Cavities with Embedded Metasurfaces

abstract: Over the past years, photonic metasurfaces have demonstrated their remarkable and diverse capabilities for achieving advanced control over light propagation by confining electromagnetic radiation within the deeply subwavelength thickness of these artificial films. Here, we demonstrate that metasurfaces also offer new unparalleled capabilities for decreasing the overall dimensions of integrated optical components and systems. We propose an original approach of embedding a metasurface inside one of the most fundamental optical elements -an optical cavity- to drastically scale-down the thickness of the optical cavity. We apply this methodology to design and implement a metasurfaces-based nano-cavity where the Fabry-Perot interferometric principle has been modified to reduce the cavity thickness below the conventional half-wavelength minimum. In addition, the nano-cavities with embedded metasurfaces can support independently tunable resonances at multiple bands. As a proof-of-concept, using nano-structured metasurfaces within 100-nm nano-cavities, we experimentally demonstrate high spatial resolution color filtering and spectral imaging. The proposed approach can be extrapolated to compact integrated optical systems on-a-chip such as VCSEL, high-resolution spatial light modulators, imaging spectroscopy systems, and bio-sensors.


title: ORBIT: Ordering Based Information Transfer Across Space and Time for Global Surface Water Monitoring

abstract: Many earth science applications require data at both high spatial and temporal resolution for effective monitoring of various ecosystem resources. Due to practical limitations in sensor design, there is often a trade-off in different resolutions of spatio-temporal datasets and hence a single sensor alone cannot provide the required information. Various data fusion methods have been proposed in the literature that mainly rely on individual timesteps when both datasets are available to learn a mapping between features values at different resolutions using local relationships between pixels. Earth observation data is often plagued with spatially and temporally correlated noise, outliers and missing data due to atmospheric disturbances which pose a challenge in learning the mapping from a local neighborhood at individual timesteps. In this paper, we aim to exploit time-independent global relationships between pixels for robust transfer of information across different scales. Specifically, we propose a new framework, ORBIT (Ordering Based Information Transfer) that uses relative ordering constraint among pixels to transfer information across both time and scales. The effectiveness of the framework is demonstrated for global surface water monitoring using both synthetic and real-world datasets.


title: Optical and J,K-photometry of black hole X-ray nova A0620-00 in passive and active stages of quiescence

abstract: Photometric observations of the low-mass X-ray binary system A0620-00=V616 Mon are performed in the optical (unfiltered light, lambda_eff~6400A) and the near-infrared J and K-bands. The mean system flux, the orbital light curve shape and the flickering amplitude dependences on wavelength are examined for two activity stages of the system remaining in quiescence. In 2015-16 A0620-00 was in passive stage (as by Cantrell et al., 2008) exhibiting the regular orbital light curves and low flickering. In less than 230 days in 2016-17 the system switched into active stage: the brightness increased by ~0.2-0.3 mag, the orbital light curve changed while the flickering amplitude increased more than twice. The object regular orbital light curves were fitted by models with "cold" spots on the optical star surface and without those. These models reproduce the observed orbital light curves both in passive and in active stages. The dependence of the mean square flickering amplitude (in fluxes, extinction corrected) on wavelength is computed in the lambda 6400-22000AA range. In active stage, the observed flickering amplitude decreases over the whole studied range and may be represented as Delta F_fl~lambda^{-2} which corresponds to the free-free emission of optically thin high-temperature plasma. In passive stage, flickering obeys Delta F_fl~lambda^{-4} law in the range 6400-12500AA that corresponds to the thermal radiation of optically thick high-temperature plasma. At longer wavelengths the flickering amplitude dependence is flat which may imply existence of a synchrotron component of the relativistic jets emission. These flickering features let us propose that the flickering mechanism includes at least two components: thermal and, apparently, synchrotron, that agrees with the recent discovery of the variable linear polarization of the IR system emission (Russell et al., 2016).


title: Progress Toward A $2\nu\beta\beta$ Measurement For The Majorana Demonstrator

abstract: The MAJORANA DEMONSTRATOR is a $^{76}$Ge-based neutrinoless double-beta decay ($0\nu\beta\beta$) experiment. Staged at the 4850 ft level of the Sanford Underground Research Facility, the DEMONSTRATOR operates an array of high-purity p-type point contact Ge detectors deployed within a graded passive shield and an active muon veto system. The present work concerns the two-neutrino double-beta decay mode ($2\nu\beta\beta$) of $^{76}$Ge. For Ge detectors, having superior energy resolution (0.1%), this mode poses negligible background to the $0\nu\beta\beta$ mode, even for a ton-scale experiment. However, the measurement of the $2\nu\beta\beta$ mode allows for careful systematics checks of active detector mass, enrichment fraction, and pulse shape discrimination cuts related to both the $0\nu\beta\beta$ and $2\nu\beta\beta$ decay modes. A precision measurement of the $2\nu\beta\beta$ shape also allows searches for spectral distortions, possibly indicative of new physics, including $0\nu\beta\beta\chi$. Work is underway to construct a full experimental background model enabling a Bayesian fit to the measured energy spectrum and extraction of a precise $2\nu\beta\beta$ spectrum and half-life.


title: Avalanche precursors of failure in hierarchical fuse networks

abstract: We study precursors of failure in hierarchical random fuse network models which can be considered as idealizations of hierarchical (bio)materials where fibrous assemblies are held together by multi-level (hierarchical) cross-links. When such structures are loaded towards failure, the patterns of precursory avalanche activity exhibit generic scale invariance: Irrespective of load, precursor activity is characterized by power-law avalanche size distributions without apparent cut-off, with power-law exponents that decrease continuously with increasing load. This failure behavior and the ensuing super-rough crack morphology differ significantly from the findings in non-hierarchical structures.


title: Suppressed solitonic cascade in spin-imbalanced superfluid Fermi gas

abstract: Cold atoms experiments offer invaluable information on superfluid dynamics, including decay cascades of topological defects. While the cascade properties are well established for Bose systems, our understanding of their behavior in Fermi counterparts is very limited, in particular in spin-imbalanced systems, where superfluid (paired) and normal (unpaired) particles naturally coexist giving rise to complex spatial structure of the atomic cloud. Here we show, based on a newly developed microscopic approach, that the decay cascades of topological defects are dramatically modified by the spin-polarization. We demonstrate that decay cascades end up at different stages: "dark soliton", "vortex ring" or "vortex line", depending on the polarization. We reveal that it is caused by sucking of unpaired particles into the soliton's internal structure. As a consequence vortex reconnections are hindered and we anticipate that quantum turbulence phenomenon can be significantly affected, indicating new physics induced by polarization effects.


title: Schwartz functions on quasi-Nash varieties

abstract: We introduce a new category called Quasi-Nash, unifying Nash manifolds and algebraic varieties. We define Schwartz functions, tempered functions and tempered distributions in this category. We show that properties that hold on affine spaces, Nash manifolds and algebraic varieties, also hold in this category.


title: Robust and Precise Vehicle Localization based on Multi-sensor Fusion in Diverse City Scenes

abstract: We present a robust and precise localization system that achieves centimeter-level localization accuracy in disparate city scenes. Our system adaptively uses information from complementary sensors such as GNSS, LiDAR, and IMU to achieve high localization accuracy and resilience in challenging scenes, such as urban downtown, highways, and tunnels. Rather than relying only on LiDAR intensity or 3D geometry, we make innovative use of LiDAR intensity and altitude cues to significantly improve localization system accuracy and robustness. Our GNSS RTK module utilizes the help of the multi-sensor fusion framework and achieves a better ambiguity resolution success rate. An error-state Kalman filter is applied to fuse the localization measurements from different sources with novel uncertainty estimation. We validate, in detail, the effectiveness of our approaches, achieving 5-10cm RMS accuracy and outperforming previous state-of-the-art systems. Importantly, our system, while deployed in a large autonomous driving fleet, made our vehicles fully autonomous in crowded city streets despite road construction that occurred from time to time. A dataset including more than 60 km real traffic driving in various urban roads is used to comprehensively test our system.


title: Chiral state-conversion without encircling an exceptional point

abstract: Dynamically varying system parameters along a path enclosing an exceptional point is known to lead to chiral mode conversion. But is it necessary to include this non-Hermitian degeneracy inside the contour for this process to take place? We show that a slow enough variation of parameters, even away from the system's exceptional point, can also lead to a robust asymmetric state exchange. To study this process, we consider a prototypical two-level non-Hermitian Hamiltonian with a constant coupling between elements. Closed form solutions are obtained when the amplification/attenuation coefficients in this arrangement are varied in conjunction with the resonance detuning along a circular contour. Using asymptotic expansions, this input-independent mode conversion is theoretically proven to take place irrespective of whether the exceptional point is enclosed or not upon encirclement. Our results significantly broaden the range of parameter space required for the experimental realization of such chiral mode conversion processes.


title: Set complexity of construction of a regular polygon

abstract: Given a subset of $\mathbb C$ containing $x,y$, one can add $x + y,\,x - y,\,xy$ or (when $y\ne0$) $x/y$ or any $z$ such that $z^2=x$. Let $p$ be a prime Fermat number. We prove that it is possible to obtain from $\{1\}$ a set containing all the $p$-th roots of 1 by $12 p^2$ above operations. This result is different from the standard estimation of complexity of an algorithm computing the $p$-th roots of 1.


title: Octahedral tilt independent magnetism in confined GdTiO$_3$ films

abstract: Polarized neutron reflectometry measurements are presented exploring the evolution of ferrimagnetism in GdTiO$_3$ films as they are confined between SrTiO$_3$ layers of variable thicknesses. As GdTiO$_3$ films approach the thin layer limit and are confined within a substantially thicker SrTiO$_3$ matrix, the TiO$_6$ octahedral tilts endemic to GdTiO$_3$ coherently relax toward the undistorted, cubic phase of SrTiO$_3$. Our measurements reveal that the ferrimagnetic state within the GdTiO$_3$ layers survives as the TiO$_6$ octahedral tilts in the GdTiO$_3$ layers are suppressed. Furthermore, our data suggest that a magnetic dead layer develops within the GdTiO$_3$ layer at each GdTiO$_3$/ SrTiO$_3$ interface. The ferrimagnetic moment inherent to the core GdTiO$_3$ layers is negligibly (in models with dead layers) or only weakly (in models without dead layers) impacted as the octahedral tilt angles are suppressed by more than 50$\%$ and the $t_{2g}$ bandwidth is dramatically renormalized.


title: Hierarchical Modeling of Seed Variety Yields and Decision Making for Future Planting Plans

abstract: Eradicating hunger and malnutrition is a key development goal of the 21st century. We address the problem of optimally identifying seed varieties to reliably increase crop yield within a risk-sensitive decision-making framework. Specifically, we introduce a novel hierarchical machine learning mechanism for predicting crop yield (the yield of different seed varieties of the same crop). We integrate this prediction mechanism with a weather forecasting model, and propose three different approaches for decision making under uncertainty to select seed varieties for planting so as to balance yield maximization and risk.We apply our model to the problem of soybean variety selection given in the 2016 Syngenta Crop Challenge. Our prediction model achieves a median absolute error of 3.74 bushels per acre and thus provides good estimates for input into the decision models.Our decision models identify the selection of soybean varieties that appropriately balance yield and risk as a function of the farmer's risk aversion level. More generally, our models support farmers in decision making about which seed varieties to plant.


title: Direct atomic fabrication and dopant positioning in Si using electron beams with active real time image-based feedback

abstract: Semiconductor fabrication is a mainstay of modern civilization, enabling the myriad applications and technologies that underpin everyday life. However, while sub-10 nanometer devices are already entering the mainstream, the end of the Moore's Law roadmap still lacks tools capable of bulk semiconductor fabrication on sub-nanometer and atomic levels, with probe-based manipulation being explored as the only known pathway. Here we demonstrate that the atomic-sized focused beam of a scanning transmission electron microscope can be used to manipulate semiconductors such as Si on the atomic level, inducing growth of crystalline Si from the amorphous phase, reentrant amorphization, milling, and dopant-front motion. These phenomena are visualized in real time with atomic resolution. We further implement active feedback control based on real-time image analytics to control the e-beam motion, enabling shape control and providing a pathway for atom-by-atom correction of fabricated structures in the near future. These observations open a new epoch for atom-by-atom manufacturing in bulk, the long-held dream of nanotechnology.


title: Completeness of the list of spinor regular ternary quadratic forms

abstract: Extending the notion of regularity introduced by Dickson in 1939, a positive definite ternary integral quadratic form is said to be spinor regular if it represents all the positive integers represented by its spinor genus (that is, all positive integers represented by any form in its spinor genus). Jagy conducted an extensive computer search for primitive ternary quadratic forms that are spinor regular, but not regular, resulting in a list of 29 such forms. In this paper, we will prove that there are no additional forms with this property.


title: Iteration complexity of inexact augmented Lagrangian methods for constrained convex programming

abstract: Augmented Lagrangian method (ALM) has been popularly used for solving constrained optimization problems. Practically, subproblems for updating primal variables in the framework of ALM usually can only be solved inexactly. The convergence and local convergence speed of ALM have been extensively studied. However, the global convergence rate of inexact ALM is still open for problems with nonlinear inequality constraints. In this paper, we work on general convex programs with both equality and inequality constraints. For these problems, we establish the global convergence rate of inexact ALM and estimate its iteration complexity in terms of the number of gradient evaluations to produce a solution with a specified accuracy. We first establish an ergodic convergence rate result of inexact ALM that uses constant penalty parameters or geometrically increasing penalty parameters. Based on the convergence rate result, we apply Nesterov's optimal first-order method on each primal subproblem and estimate the iteration complexity of the inexact ALM. We show that if the objective is convex, then $O(\varepsilon^{-1})$ gradient evaluations are sufficient to guarantee an $\varepsilon$-optimal solution in terms of both primal objective and feasibility violation. If the objective is strongly convex, the result can be improved to $O(\varepsilon^{-\frac{1}{2}}|\log\varepsilon|)$. Finally, by relating to the inexact proximal point algorithm, we establish a nonergodic convergence rate result of inexact ALM that uses geometrically increasing penalty parameters. We show that the nonergodic iteration complexity result is in the same order as that for the ergodic result. Numerical experiments on quadratically constrained quadratic programming are conducted to compare the performance of the inexact ALM with different settings.


title: Galaxy structure from multiple tracers: III. Radial variations in M87's IMF

abstract: We present the first constraints on stellar mass-to-light ratio gradients in an early-type galaxy (ETG) using multiple dynamical tracer populations to model the dark and luminous mass structure simultaneously. We combine the kinematics of the central starlight, two globular cluster populations and satellite galaxies in a Jeans analysis to obtain new constraints on M87's mass structure, employing a flexible mass model which allows for radial gradients in the stellar mass-to-light ratio. We find that, in the context of our model, a radially declining stellar-mass-to-light ratio is strongly favoured. Modelling the stellar mass-to-light ratio as following a power law, $\Upsilon_{\star} \sim R^{-\mu}$, we infer a power-law slope $\mu = -0.54 \pm 0.05$; equally, parameterising the stellar-mass-to-light ratio via a central mismatch parameter relative to a Salpeter IMF, $\alpha$, and scale radius $R_M$, we find $\alpha > 1.48$ at $95\%$ confidence and $R_M = 0.35 \pm 0.04$ kpc. We use stellar population modelling of high-resolution 11-band HST photometry to show that such a steep gradient cannot be achieved by variations in only the metallicity, age, dust extinction and star formation history if the stellar initial mass function (IMF) remains spatially constant. On the other hand, the stellar mass-to-light ratio gradient that we find is consistent with an IMF whose inner slope changes such that it is Salpeter-like in the central $\sim 0.5$ kpc and becomes Chabrier-like within the stellar effective radius. This adds to recent evidence that the non-universality of the IMF in ETGs may be confined to their core regions, and points towards a picture in which the stars in these central regions may have formed in fundamentally different physical conditions.


title: Python Implementation and Construction of Finite Abelian Groups

abstract: Here we present a working framework to establish finite abelian groups in python. The primary aim is to allow new A-level students to work with examples of finite abelian groups using open source software. We include the code used in the implementation of the framework. We also prove some useful results regarding finite abelian groups which are used to establish the functions and help show how number theoretic results can blend with computational power when studying algebra. The groups established are based modular multiplication and addition. We include direct products of cyclic groups meaning the user has access to all finite abelian groups.


title: An energy-based equilibrium contact angle boundary condition on jagged surfaces for phase-field methods

abstract: We consider an energy-based boundary condition to impose an equilibrium wetting angle for the Cahn-Hilliard-Navier-Stokes phase-field model on voxel-set-type computational domains. These domains typically stem from the micro-CT imaging of porous rock and approximate a (on {\mu}m scale) smooth domain with a certain resolution. Planar surfaces that are perpendicular to the main axes are naturally approximated by a layer of voxels. However, planar surfaces in any other directions and curved surfaces yield a jagged/rough surface approximation by voxels. For the standard Cahn-Hilliard formulation, where the contact angle between the diffuse interface and the domain boundary (fluid-solid interface/wall) is 90 degrees, jagged surfaces have no impact on the contact angle. However, a prescribed contact angle smaller or larger than 90 degrees on jagged voxel surfaces is amplified in either direction. As a remedy, we propose the introduction of surface energy correction factors for each fluid-solid voxel face that counterbalance the difference of the voxel-set surface area with the underlying smooth one. The discretization of the model equations is performed with the discontinuous Galerkin method, however, the presented semi-analytical approach of correcting the surface energy is equally applicable to other direct numerical methods such as finite elements, finite volumes, or finite differences, since the correction factors appear in the strong formulation of the model.


title: K3, L3, LP, RM3, A3, FDE: How to Make Many-Valued Logics Work for You

abstract: We investigate some well-known (and a few not-so-well-known) many-valued logics that have a small number (3 or 4) of truth values. For some of them we complain that they do not have any \emph{logical} use (despite their perhaps having some intuitive semantic interest) and we look at ways to add features so as to make them useful, while retaining their intuitive appeal. At the end, we show some surprising results in the system FDE, and its relationships with features of other logics. We close with some new examples of "synonymous logics." An Appendix contains a natural deduction system for our augmented FDE, and proofs of soundness and completeness.


title: Costate-focused models for reinforcement learning

abstract: Many recent algorithms for reinforcement learning are model-free and founded on the Bellman equation. Here we present a method founded on the costate equation and models of the state dynamics. We use the costate -- the gradient of cost with respect to state -- to improve the policy and also to "focus" the model, training it to detect and mimic those features of the environment that are most relevant to its task. We show that this method can handle difficult time-optimal control problems, driving deterministic or stochastic mechanical systems quickly to a target. On these tasks it works well compared to deep deterministic policy gradient, a recent Bellman method. And because it creates a model, the costate method can also learn from mental practice.


title: Fronthaul-Aware Group Sparse Precoding and Signal Splitting in SWIPT C-RAN

abstract: We investigate the precoding, remote radio head (RRH) selection and signal splitting in the simultaneous wireless information and power transferring (SWIPT) cloud radio access networks \mbox{(C-RANs)}. The objective is to minimize the power consumption of the SWIPT C-RAN. Different from the existing literature, we consider the nonlinear fronthaul power consumption and the multiple antenna RRHs. By switching off the unnecessary RRHs, the group sparsity of the precoding coefficients is introduced, which indicates that the precoding process and the RRH selection are coupled. In order to overcome these issues, a group sparse precoding and signal splitting algorithm is proposed based on the majorization-minimization framework, and the convergence behavior is established. Numerical results are used to verify our proposed studies.


title: Probing Split-Ring Resonator Permeabilities with Loop-Gap Resonators

abstract: A method is proposed to experimentally determine the effective complex permeability of split-ring resonator (SRR) arrays used in the design of metamaterials at microwave frequencies. We analyze the microwave response of a loop-gap resonator (LGR) whose bore has been partially loaded with one or more SRRs. Our analysis reveals that the resonance frequency, magnetic plasma frequency, and damping constant of the effective permeability of the SRR array can be extracted from fits to the reflection coefficient (S11) of an inductively-coupled LGR. We propose LGR designs that would allow both a one-dimensional array of SRRs and small three-dimensional arrays of SRRs to be characterized. Finally, we demonstrate the method using a toroidal LGR loaded with a single extended SRR of length z.


title: Zero-Shot Learning via Class-Conditioned Deep Generative Models

abstract: We present a deep generative model for learning to predict classes not seen at training time. Unlike most existing methods for this problem, that represent each class as a point (via a semantic embedding), we represent each seen/unseen class using a class-specific latent-space distribution, conditioned on class attributes. We use these latent-space distributions as a prior for a supervised variational autoencoder (VAE), which also facilitates learning highly discriminative feature representations for the inputs. The entire framework is learned end-to-end using only the seen-class training data. The model infers corresponding attributes of a test image by maximizing the VAE lower bound; the inferred attributes may be linked to labels not seen when training. We further extend our model to a (1) semi-supervised/transductive setting by leveraging unlabeled unseen-class data via an unsupervised learning module, and (2) few-shot learning where we also have a small number of labeled inputs from the unseen classes. We compare our model with several state-of-the-art methods through a comprehensive set of experiments on a variety of benchmark data sets.


title: SV-map between Type I and Heterotic Sigma Models

abstract: The scattering amplitudes of gauge bosons in heterotic and open superstring theories are related by the single-valued projection which yields heterotic amplitudes by selecting a subset of multiple zeta value coefficients in the $\alpha'$ (string tension parameter) expansion of open string amplitudes. In the present work, we argue that this relation holds also at the level of low-energy expansions (or individual Feynman diagrams) of the respective effective actions, by investigating the beta functions of two-dimensional sigma models describing world-sheets of open and heterotic strings. We analyze the sigma model Feynman diagrams generating identical effective action terms in both theories and show that the heterotic coefficients are given by the single-valued projection of the open ones. The single-valued projection appears as a result of summing over all radial orderings of heterotic vertices on the complex plane representing string world-sheet.


title: Understanding the Changing Roles of Scientific Publications via Citation Embeddings

abstract: Researchers may describe different aspects of past scientific publications in their publications and the descriptions may keep changing in the evolution of science. The diverse and changing descriptions (i.e., citation context) on a publication characterize the impact and contributions of the past publication. In this article, we aim to provide an approach to understanding the changing and complex roles of a publication characterized by its citation context. We described a method to represent the publications' dynamic roles in science community in different periods as a sequence of vectors by training temporal embedding models. The temporal representations can be used to quantify how much the roles of publications changed and interpret how they changed. Our study in the biomedical domain shows that our metric on the changes of publications' roles is stable over time at the population level but significantly distinguish individuals. We also show the interpretability of our methods by a concrete example.


title: The holomorphic bosonic string

abstract: We present a holomorphic version of the bosonic string in the formalism of quantum field theory developed by Costello and collaborators. In this paper we focus on the case in which space-time is flat and construct a one-loop exact quantization. Starting from first principles, we arrive at the critical dimension as an obstruction to satisfying the quantum master equation. Moreover, we show how the factorization algebra recovers the BRST cohomology of the string and give another construction of the Gerstenhaber structure. Finally, we show how the factorization homology along closed manifolds encodes the determinant line bundle over the moduli space of Riemann surfaces.


title: Security Issues in Controller Area Networks in Automobiles

abstract: Modern vehicles may contain a considerable number of ECUs (Electronic Control Units) which are connected through various means of communication, with the CAN (Controller Area Network) protocol being the most widely used. However, several vulnerabilities such as the lack of authentication and the lack of data encryption have been pointed out by several authors, which ultimately render vehicles unsafe to their users and surroundings. Moreover, the lack of security in modern automobiles has been studied and analyzed by other researchers as well as several reports about modern car hacking have (already) been published. The contribution of this work aimed to analyze and test the level of security and how resilient is the CAN protocol by taking a BMW E90 (3-series) instrument cluster as a sample for a proof of concept study. This investigation was carried out by building and developing a rogue device using cheap commercially available components while being connected to the same CAN-Bus as a man in the middle device in order to send spoofed messages to the instrument cluster.


title: Bootstrapped synthetic likelihood

abstract: Approximate Bayesian computation (ABC) and synthetic likelihood (SL) techniques have enabled the use of Bayesian inference for models that may be simulated, but for which the likelihood cannot be evaluated pointwise at values of an unknown parameter $\theta$. The main idea in ABC and SL is to, for different values of $\theta$ (usually chosen using a Monte Carlo algorithm), build estimates of the likelihood based on simulations from the model conditional on $\theta$. The quality of these estimates determines the efficiency of an ABC/SL algorithm. In standard ABC/SL, the only means to improve an estimated likelihood at $\theta$ is to simulate more times from the model conditional on $\theta$, which is infeasible in cases where the simulator is computationally expensive. In this paper we describe how to use bootstrapping as a means for improving SL estimates whilst using fewer simulations from the model, and also investigate its use in ABC. Further, we investigate the use of the bag of little bootstraps as a means for applying this approach to large datasets, yielding Monte Carlo algorithms that accurately approximate posterior distributions whilst only simulating subsamples of the full data. Examples of the approach applied to i.i.d., temporal and spatial data are given.

categories: stat.CO cs.AI physics.data-an stat.ME stat.ML


title: Finding Efficient Swimming Strategies in a Three Dimensional Chaotic Flow by Reinforcement Learning

abstract: We apply a reinforcement learning algorithm to show how smart particles can learn approximately optimal strategies to navigate in complex flows. In this paper we consider microswimmers in a paradigmatic three-dimensional case given by a stationary superposition of two Arnold-Beltrami-Childress flows with chaotic advection along streamlines. In such a flow, we study the evolution of point-like particles which can decide in which direction to swim, while keeping the velocity amplitude constant. We show that it is sufficient to endow the swimmers with a very restricted set of actions (six fixed swimming directions in our case) to have enough freedom to find efficient strategies to move upward and escape local fluid traps. The key ingredient is the learning-from-experience structure of the algorithm, which assigns positive or negative rewards depending on whether the taken action is, or is not, profitable for the predetermined goal in the long term horizon. This is another example supporting the efficiency of the reinforcement learning approach to learn how to accomplish difficult tasks in complex fluid environments.

categories: physics.flu-dyn cond-mat.stat-mech cs.LG nlin.CD


title: Quantum Dragon Solutions for Electron Transport through Nanostructures based on Rectangular Graphs

abstract: Electron transport through nanodevices of atoms in a single-layer rectangular arrangement with free (open) boundary conditions parallel to the direction of the current flow is studied within the single-band tight binding model. The Landauer formula gives the electrical conductance to be a function of the electron transmission probability, ${\cal T}(E)$, as a function of the energy $E$ of the incoming electron. A quantum dragon nanodevice is one which has a perfectly conducting channel, namely ${\cal T}(E)=1$ for all energies which are transmitted by the external leads even though there may be arbitrarily strong electron scattering. The rectangular single-layer systems are shown to be able to be quantum dragon devices, both for uniform leads and for dimerized leads. The quantum dragon condition requires appropriate lead-device connections and correlated randomness in the device.


title: BoostJet: Towards Combining Statistical Aggregates with Neural Embeddings for Recommendations

abstract: Recommenders have become widely popular in recent years because of their broader applicability in many e-commerce applications. These applications rely on recommenders for generating advertisements for various offers or providing content recommendations. However, the quality of the generated recommendations depends on user features (like demography, temporality), offer features (like popularity, price), and user-offer features (like implicit or explicit feedback). Current state-of-the-art recommenders do not explore such diverse features concurrently while generating the recommendations. In this paper, we first introduce the notion of Trackers which enables us to capture the above-mentioned features and thus incorporate users' online behaviour through statistical aggregates of different features (demography, temporality, popularity, price). We also show how to capture offer-to-offer relations, based on their consumption sequence, leveraging neural embeddings for offers in our Offer2Vec algorithm. We then introduce BoostJet, a novel recommender which integrates the Trackers along with the neural embeddings using MatrixNet, an efficient distributed implementation of gradient boosted decision tree, to improve the recommendation quality significantly. We provide an in-depth evaluation of BoostJet on Yandex's dataset, collecting online behaviour from tens of millions of online users, to demonstrate the practicality of BoostJet in terms of recommendation quality as well as scalability.


title: The Mpemba index and anomalous relaxation

abstract: The Mpemba effect is a counter-intuitive relaxation phenomenon, where a system prepared at a hot temperature cools down faster than an identical system initiated at a cold temperature when both are quenched to an even colder bath. Such non-monotonic relaxations were observed in various systems, including water, magnetic alloys, polymers, and driven granular gases. We analyze the Mpemba effect in Markovian dynamics and discover that a stronger version of the effect often exists for a carefully chosen set of initial temperatures. In this \emph{strong Mpemba effect}, the relaxation time jumps to a smaller value leading to exponentially faster equilibration dynamics. The number of such special initial temperatures defines the \emph{Mpemba index}, whose parity is a topological property of the system. To demonstrate these concepts, we first analyze the different types of Mpemba relaxations in the mean field anti-ferromagnet Ising model, which demonstrates a surprisingly rich Mpemba phase diagram. Moreover, we show that the strong effect survives the thermodynamic limit and that it is tightly connected with thermal overshoot -- in the relaxation process, the temperature of the relaxing system can decay non-monotonically as a function of time. Using the parity of the Mpemba index, we then study the occurrence of the strong Mpemba effect in a large class of thermal quench processes and show that it happens with non-zero probability even in the thermodynamic limit. This is done by introducing the \emph{isotropic} model for which we obtain analytical lower bound estimates for the probability of the strong Mpemba effects. Consequently, we expect that such exponentially faster relaxations can be observed experimentally in a wide variety of systems.


title: The moment sum-rules for ionic liquids at criticality

abstract: We discuss the first three well known moment charge-charge sum-rules for a general ionic liquid. For the special symmetric case of the Restricted Primitive Model, Das et al. [Phys. Rev. Lett. 107, 215701 (2011)] has recently discovered, through Monte Carlo simulations, that the Stillinger-Lovett or second-moment sum-rule fails at criticality. We critically discuss a possible explanation for this unexpected behavior. On the other hand the fourth-moment sum-rule turns out to be able to account for the results of the simulations at criticality.


title: Anisotropy of the Seebeck and Nernst coefficients in parent compounds of the iron-based superconductors

abstract: In-plane longitudinal and transverse thermoelectric phenomena in two parent compounds of iron-based superconductors are studied. Namely, the Seebeck (S) and Nernst (n) coefficients were measured in the temperature range 10 - 300 K for BaFe2As2 and CaFe2As2 single crystals that were detwinned in-situ. The thermoelectric response shows sizeable anisotropy in the spin density wave state (SDW) for both compounds, while some dissimilarities in the vicinity of the SDW transition can be attributed to the different nature of the phase change in BaFe2As2 and CaFe2As2. Temperature dependences of S and n can be described within a two-band model that contains a contribution from highly mobile, probably Dirac, electrons. The Dirac band seems to be rather isotropic, whereas most of the anisotropy in the transport phenomena could be attributed to "regular" hole-like charge carriers. We also observe that the off-diagonal element of the Peltier tensor axy is not the same for the a and b orthorhombic axes, which indicates that the widely used Mott formula is not applicable to the SDW state of iron-based superconductors.


title: On the Duflot filtration for equivariant cohomology rings and applications to group cohomology

abstract: We study the Duflot filtration on the Borel equivariant cohomology of smooth manifolds with a smooth $p$-torus action. We axiomatize the filtration and prove analog of several structural results about equivariant cohomology rings in this setting. We apply this abstract theory to study the $\mathbb{F}_p$ cohomology rings of classifying spaces of compact Lie groups, and show how to recover geometric results about the cohomology of $BG$ using equivariant cohomology. This includes some results about detection on subgroups and restrictions on associated primes that were previously only known for finite groups. We are particularly interested in the local cohomology modules of equivariant cohomology rings, and we construct a tractable chain complex computing local cohomology. As an application, we study the local cohomology of the group cohomology of the p-Sylow subgroups of $S_{p^n}$ and give vanishing and nonvanishing results for these local cohomology modules that are sharper than those given by the current theory.


title: Robustness of ballistic transport in antidot superlattices

abstract: The magneto-resistance of antidot lattices shows pronounced peaks, which became a hallmark of ballistic electron transport. While most studies agree that they reflect the interplay of regular and chaotic motion in the quasi-classical dynamics, the exact mechanism has been surprisingly controversial. Inspired by recent experiments on graphene antidot lattices showing that the effect survives strong impurity scattering, we give a new explanation of the peaks linked to a fundamental relation between collision times and accessible phase space volumes, accounting for their robustness. Due to the fundamental nature of the mechanism described it will be relevant in many mesoscopic transport phenomena.


title: Longitudinal associations between learning assistants and instructor effectiveness

abstract: A central goal of the Learning Assistant (LA) model is to improve students' learning of science through the transformation of instructor practices. There is minimal existing research on the impact of college physics instructor experiences on their effectiveness. To investigate the association between college introductory physics instructors' experiences with and without LAs and student learning, we drew on data from the Learning About STEM Student Outcomes (LASSO) database. The LASSO database provided us with student-level data (concept inventory scores and demographic data) for 4,365 students and course-level data (instructor experience and course features) for the students' 93 mechanics courses. We performed Hierarchical Multiple Imputation to impute missing data and Hierarchical Linear Modeling to nest students within courses when modeling the associations between instructor experience and student learning. Our models predict that instructors' effectiveness decreases as they gain experience teaching without LAs. However, LA supported environments appear to remediate this decline in effectiveness as instructor effectiveness is maintained while they gain experience teaching with LAs.


title: Optical activity in chiral stacks of 2D semiconductors

abstract: We show that the stacks of two-dimensional semiconductor crystals with the chiral packing exhibit optical activity and circular dichroism. We develop a microscopic theory of these phenomena in the spectral range of exciton transitions which takes into account the spin-dependent hopping of excitons between the layers in the stack and the interlayer coupling of excitons via electromagnetic field. For the stacks of realistic two-dimensional semiconductors such as transition metal dichalcogenides, we calculate the rotation and ellipticity angles of radiation transmitted through such structures. The angles are resonantly enhanced at the frequencies of both bright and dark exciton modes in the stack. We also study the photoluminescence of chiral stacks and show that it is circularly polarized.


title: Systemic inequities in introductory physics courses: the impacts of learning assistants

abstract: Creating equitable performance outcomes among students is a focus of many instructors and researchers. One focus of this effort is examining disparities in physics student performance across genders, which is a well-established problem. Another less common focus is disparities across racial and ethnic groups, which may have received less attention due to low representation rates making it difficult to identify gaps in their performance. In this investigation we examined associations between Learning Assistant (LA) supported courses and improved equity in student performance. We built Hierarchical Linear Models of student performance to investigate how performance differed by gender and by race/ethnicity and how LAs may have moderated those differences. Data for the analysis came from pre-post concept inventories in introductory mechanics courses collected through the Learning About STEM Student Outcomes (LASSO) platform. Our models show that gaps in performance across genders and races/ethnicities were similar in size and increased from pre to post instruction. LA-support is meaningfully and reliably associated with improvement in overall student performance but not with shifts in within-course performance gaps.


title: Crowd Counting Through Walls Using WiFi

abstract: Counting the number of people inside a building, from outside and without entering the building, is crucial for many applications. In this paper, we are interested in counting the total number of people walking inside a building (or in general behind walls), using readily-deployable WiFi transceivers that are installed outside the building, and only based on WiFi RSSI measurements. The key observation of the paper is that the inter-event times, corresponding to the dip events of the received signal, are fairly robust to the attenuation through walls (for instance as compared to the exact dip values). We then propose a methodology that can extract the total number of people from the inter-event times. More specifically, we first show how to characterize the wireless received power measurements as a superposition of renewal-type processes. By borrowing theories from the renewal-process literature, we then show how the probability mass function of the inter-event times carries vital information on the number of people. We validate our framework with 44 experiments in five different areas on our campus (3 classrooms, a conference room, and a hallway), using only one WiFi transmitter and receiver installed outside of the building, and for up to and including 20 people. Our experiments further include areas with different wall materials, such as concrete, plaster, and wood, to validate the robustness of the proposed approach. Overall, our results show that our approach can estimate the total number of people behind the walls with a high accuracy while minimizing the need for prior calibrations.


title: Participation rates of in-class vs. online administration of low-stakes research-based assessments

abstract: This study investigates differences in student participation rates between in-class and online administrations of research-based assessments. A sample of 1,310 students from 25 sections of 3 different introductory physics courses over two semesters were instructed to complete the CLASS attitudinal survey and the concept inventory relevant to their course, either the FCI or the CSEM. Each student was randomly assigned to take one of the surveys in class and the other survey online at home using the Learning About STEM Student Outcomes (LASSO) platform. Results indicate large variations in participation rates across both test conditions (online and in class). A hierarchical generalized linear model (HGLM) of the student data utilizing logistic regression indicates that student grades in the course and faculty assessment administration practices were both significant predictors of student participation. When the recommended online assessments administration practices were implemented, participation rates were similar across test conditions. Implications for student and course assessment methodologies will be discussed.


title: Exact and heuristic algorithms for Cograph Editing

abstract: We present a dynamic programming algorithm for optimally solving the Cograph Editing problem on an $n$-vertex graph that runs in $O(3^n n)$ time and uses $O(2^n)$ space. In this problem, we are given a graph $G = (V, E)$ and the task is to find a smallest possible set $F \subseteq V \times V$ of vertex pairs such that $(V, E \bigtriangleup F)$ is a cograph (or $P_4$-free graph), where $\bigtriangleup$ represents the symmetric difference operator. We also describe a technique for speeding up the performance of the algorithm in practice. Additionally, we present a heuristic for solving the Cograph Editing problem which produces good results on small to medium datasets. In application it is much more important to find the ground truth, not some optimal solution. For the first time, we evaluate whether the cograph property is strict enough to recover the true graph from data to which noise has been added.


title: Least informative distributions in Maximum q-log-likelihood estimation

abstract: We use the Maximum $q$-log-likelihood estimation for Least informative distributions (LID) in order to estimate the parameters in probability density functions (PDFs) efficiently and robustly when data include outlier(s). LIDs are derived by using convex combinations of two PDFs, $f_\epsilon=(1-\epsilon)f_0+\epsilon f_1$. A convex combination of two PDFs is considered as a contamination $f_1$ as outlier(s) to underlying $f_0$ distributions and $f_\epsilon$ is a contaminated distribution. The optimal criterion is obtained by minimizing the change of Maximum q-log-likelihood function when the data have slightly more contamination. In this paper, we make a comparison among ordinary Maximum likelihood, Maximum q-likelihood estimations, LIDs based on $\log_q$ and Huber M-estimation. Akaike and Bayesian information criterions (AIC and BIC) based on $\log_q$ and LID are proposed to assess the fitting performance of functions. Real data sets are applied to test the fitting performance of estimating functions that include shape, scale and location parameters.


title: On monotonicity of some functionals with variable exponent under symmetrization

abstract: We give necessary and sufficient conditions for the P\'olya--Szeg\"o type inequality with variable exponent of summability.


title: Evaluation of Certain Hypergeometric Functions over Finite Fields

abstract: For an odd prime $p$, let $\phi$ denote the quadratic character of the multiplicative group ${\mathbb F}_p^\times$, where ${\mathbb F}_p$ is the finite field of $p$ elements. In this paper, we will obtain evaluations of the hypergeometric functions $ {}_2F_1\left(\begin{matrix} \phi\psi & \psi\\ & \phi \end{matrix};x\right)$, $x\in {\mathbb F}_p$, $x\neq 0, 1$, over ${\mathbb F}_p$ in terms of Hecke character attached to CM elliptic curves for characters $\psi$ of ${\mathbb F}_p^\times$ of order $3$, $4$, $6$, $8$, and $12$.


title: Per-sample immunoglobulin germline inference from B cell receptor deep sequencing data

abstract: The collection of immunoglobulin genes in an individual's germline, which gives rise to B cell receptors via recombination, is known to vary significantly across individuals. In humans, for example, each individual has only a fraction of the several hundred known V alleles. Furthermore, the currently-accepted set of known V alleles is both incomplete (particularly for non-European samples), and contains a significant number of spurious alleles. The resulting uncertainty as to which immunoglobulin alleles are present in any given sample results in inaccurate B cell receptor sequence annotations, and in particular inaccurate inferred naive ancestors. In this paper we first show that the currently widespread practice of aligning each sequence to its closest match in the full set of IMGT alleles results in a very large number of spurious alleles that are not in the sample's true set of germline V alleles. We then describe a new method for inferring each individual's germline gene set from deep sequencing data, and show that it improves upon existing methods by making a detailed comparison on a variety of simulated and real data samples. This new method has been integrated into the partis annotation and clonal family inference package, available at https://github.com/psathyrella/partis, and is run by default without affecting overall run time.


title: Quadric Complexes

abstract: Quadric complexes are square complexes satisfying a certain combinatorial nonpositive curvature condition. These complexes generalize 2-dimensional CAT(0) cube complexes and are a square analog of systolic complexes. We introduce and study the basic properties of these complexes. Using a form of dismantlability for the 1-skeleta of finite quadric complexes we show that every finite group acting on a quadric complex stabilizes a complete bipartite subgraph of its 1-skeleton. Finally, we prove that C(4)-T(4) small cancellation groups act on quadric complexes.


title: Trace spaces of counterexamples to Naimark's Problem

abstract: A counterexample to Naimark's problem is a $C^\ast$-algebra that is not isomorphic to the algebra of compact operators on some Hilbert space, yet still has only one irreducible representation up to unitary equivalence. It is well-known that such algebras must be nonseparable, and in 2004 Akemann and Weaver used the diamond principle (a set theoretic axiom independent from ZFC) to give the first counterexamples. For any such counterexample $A$, the unitary group $U(A)$ acts transitively on the pure states, which are the extreme points of the state space $S(A)$. It is conceivable that this implies (as happens for finite-dimensional simplexes) that the action of $U(A)$ on $S(A)$ has at most one fixed point, i.e. $A$ has at most one trace. We give a strong negative answer here assuming diamond. In particular, we adapt the Akemann-Weaver construction to show that the trace space of a counterexample to Naimark's problem can be affinely homeomorphic to any metrizable Choquet simplex, and can also be nonseparable.


title: Explicit Chabauty-Kim for the Split Cartan Modular Curve of Level 13

abstract: We extend the explicit quadratic Chabauty methods developed in previous work by the first two authors to the case of non-hyperelliptic curves. This results in an algorithm to compute the rational points on a curve of genus $g \ge 2$ over the rationals whose Jacobian has Mordell-Weil rank $g$ and Picard number greater than one, and which satisfies some additional conditions. This algorithm is then applied to the modular curve $X_{s}(13)$, completing the classification of non-CM elliptic curves over $\mathbf{Q}$ with split Cartan level structure due to Bilu-Parent and Bilu-Parent-Rebolledo.


title: AOGNets: Compositional Grammatical Architectures for Deep Learning

abstract: Neural architectures are the foundation for improving performance of deep neural networks (DNNs). This paper presents deep compositional grammatical architectures which harness the best of two worlds: grammar models and DNNs. The proposed architectures integrate compositionality and reconfigurability of the former and the capability of learning rich features of the latter in a principled way. We utilize AND-OR Grammar (AOG) as network generator in this paper and call the resulting networks AOGNets. An AOGNet consists of a number of stages each of which is composed of a number of AOG building blocks. An AOG building block splits its input feature map into N groups along feature channels and then treat it as a sentence of N words. It then jointly realizes a phrase structure grammar and a dependency grammar in bottom-up parsing the "sentence" for better feature exploration and reuse. It provides a unified framework for the best practices developed in state-of-the-art DNNs. In experiments, AOGNet is tested in the CIFAR-10, CIFAR-100 and ImageNet-1K classification benchmark and the MS-COCO object detection and segmentation benchmark. In CIFAR-10, CIFAR-100 and ImageNet-1K, AOGNet obtains better performance than ResNet and most of its variants, ResNeXt and its attention based variants such as SENet, DenseNet and DualPathNet. AOGNet also obtains the best model interpretability score using network dissection. AOGNet further shows better potential in adversarial defense. In MS-COCO, AOGNet obtains better performance than the ResNet and ResNeXt backbones in Mask R-CNN.


title: Knowledge transfer for surgical activity prediction

abstract: Lack of training data hinders automatic recognition and prediction of surgical activities necessary for situation-aware operating rooms. We propose using knowledge transfer to compensate for data deficit and improve prediction. We used two approaches to extract and transfer surgical process knowledge. First, we encoded semantic information about surgical terms using word embedding which boosted learning process. Secondly, we passed knowledge between different clinical datasets of neurosurgical procedures using transfer learning. Transfer learning was shown to be more effective than a simple combination of data, especially for less similar procedures. The combination of two methods provided 22% improvement of activity prediction. We also made several pertinent observations about surgical practices.


title: Quantifying the effect of hydrogen on dislocation dynamics: A three-dimensional discrete dislocation dynamics framework

abstract: We present a new framework to quantify the effect of hydrogen on dislocations using large scale three-dimensional (3D) discrete dislocation dynamics (DDD) simulations. In this model, the first order elastic interaction energy associated with the hydrogen-induced volume change is accounted for. The three-dimensional stress tensor induced by hydrogen concentration, which is in equilibrium with respect to the dislocation stress field, is derived using the Eshelby inclusion model, while the hydrogen diffusion is treated as a continuum process. This newly developed framework is utilized to quantify the effect of different hydrogen concentrations on the dynamics of a glide dislocation in the absence of an applied stress field as well as on the spacing between dislocations in an array of parallel edge dislocations. A shielding effect is observed for materials having a large hydrogen diffusion coefficient, with the shield effect leading to the homogenization of the shrinkage process leading to the glide loop maintaining its circular shape, as well as resulting in a decrease in dislocation separation distances in the array of parallel edge dislocations. On the other hand, for materials having a small hydrogen diffusion coefficient, the high hydrogen concentrations around the edge characters of the dislocations act to pin them. Higher stresses are required to be able to unpin the dislocations from the hydrogen clouds surrounding them. Finally, this new framework can open the door for further large scale studies on the effect of hydrogen on the different aspects of dislocation-mediated plasticity in metals. With minor modifications of the current formulations, the framework can also be extended to account for general inclusion-induced stress field in discrete dislocation dynamics simulations.


title: Local eigenvalue statistics of one-dimensional random non-selfadjoint pseudo-differential operators

abstract: We consider a class of one-dimensional nonselfadjoint semiclassical pseudo-differential operators, subject to small random perturbations, and study the statistical properties of their (discrete) spectra, in the semiclassical limit $h\to 0$. We compare two types of random perturbation: a random potential vs. a random matrix. Hager and Sj\"ostrand had shown that, with high probability, the local spectral density of the perturbed operator follows a semiclassical form of Weyl's law, depending on the value distribution of the principal symbol of our pseudodifferential operator. Beyond the spectral density, we investigate the full local statistics of the perturbed spectrum, and show that it satisfies a form of universality: the statistical only depends on the local spectral density, and of the type of random perturbation, but it is independent of the precise law of the perturbation. This local statistics can be described in terms of the Gaussian Analytic Function, a classical ensemble of random entire functions.


title: Go for a Walk and Arrive at the Answer: Reasoning Over Paths in Knowledge Bases using Reinforcement Learning

abstract: Knowledge bases (KB), both automatically and manually constructed, are often incomplete --- many valid facts can be inferred from the KB by synthesizing existing information. A popular approach to KB completion is to infer new relations by combinatory reasoning over the information found along other paths connecting a pair of entities. Given the enormous size of KBs and the exponential number of paths, previous path-based models have considered only the problem of predicting a missing relation given two entities or evaluating the truth of a proposed triple. Additionally, these methods have traditionally used random paths between fixed entity pairs or more recently learned to pick paths between them. We propose a new algorithm MINERVA, which addresses the much more difficult and practical task of answering questions where the relation is known, but only one entity. Since random walks are impractical in a setting with combinatorially many destinations from a start node, we present a neural reinforcement learning approach which learns how to navigate the graph conditioned on the input query to find predictive paths. Empirically, this approach obtains state-of-the-art results on several datasets, significantly outperforming prior methods.


title: Apprentice: Using Knowledge Distillation Techniques To Improve Low-Precision Network Accuracy

abstract: Deep learning networks have achieved state-of-the-art accuracies on computer vision workloads like image classification and object detection. The performant systems, however, typically involve big models with numerous parameters. Once trained, a challenging aspect for such top performing models is deployment on resource constrained inference systems - the models (often deep networks or wide networks or both) are compute and memory intensive. Low-precision numerics and model compression using knowledge distillation are popular techniques to lower both the compute requirements and memory footprint of these deployed models. In this paper, we study the combination of these two techniques and show that the performance of low-precision networks can be significantly improved by using knowledge distillation techniques. Our approach, Apprentice, achieves state-of-the-art accuracies using ternary precision and 4-bit precision for variants of ResNet architecture on ImageNet dataset. We present three schemes using which one can apply knowledge distillation techniques to various stages of the train-and-deploy pipeline.


title: Smart Inertial Particles

abstract: We performed a numerical study to train smart inertial particles to target specific flow regions with high vorticity through the use of reinforcement learning algorithms. The particles are able to actively change their size to modify their inertia and density. In short, using local measurements of the flow vorticity, the smart particle explores the interplay between its choices of size and its dynamical behaviour in the flow environment. This allows it to accumulate experience and learn approximately optimal strategies of how to modulate its size in order to reach the target high-vorticity regions. We consider flows with different complexities: a two-dimensional stationary Taylor-Green like configuration, a two-dimensional time-dependent flow, and finally a three-dimensional flow given by the stationary Arnold-Beltrami-Childress helical flow. We show that smart particles are able to learn how to reach extremely intense vortical structures in all the tackled cases.


title: Logarithmic Kodaira dimension and zeros of holomorphic log-one-forms

abstract: In this paper, we prove that the zero-locus of any global holomorphic log-one-form on a projective log-smooth pair $\left(X,D\right)$ of log-general type must be non-empty. Applying this result, we give an answer to the algebraic hyperbolicity part of Shafarevich's conjecture, with the generic fiber being Kawamata-log-terminal (klt) and of log-general type.


title: Passive Crowd Speed Estimation in Adjacent Regions With Minimal WiFi Sensing

abstract: In this paper, we propose a methodology for estimating the crowd speed using WiFi devices without relying on people to carry any device. Our approach not only enables speed estimation in the region where WiFi links are, but also in the adjacent possibly WiFi-free regions. More specifically, we use a pair of WiFi links in one region, whose RSSI measurements are then used to estimate the crowd speed, not only in this region, but also in adjacent WiFi-free regions. We first prove how the cross-correlation and the probability of crossing the two links implicitly carry key information about the pedestrian speeds and develop a mathematical model to relate them to pedestrian speeds. We then validate our approach with 108 experiments, in both indoor and outdoor, where up to 10 people walk in two adjacent areas, with variety of speeds per region, showing that our framework can accurately estimate these speeds with only a pair of WiFi links in one region. For instance, the NMSE over all experiments is 0.18. We also evaluate our framework in a museum-type setting and estimate the popularity of different exhibits. We finally run experiments in an aisle in Costco, estimating key attributes of buyers' behaviors.


title: Photogravimagnetic assists of light sails: a mixed blessing for Breakthrough Starshot?

abstract: Upon entering a star system, light sails are subject to both gravitational forces and radiation pressure, and can use both in concert to modify their trajectory. Moreover, stars possess significant magnetic fields, and if the sail is in any way charged, it will feel the Lorentz force also. We investigate the dynamics of so-called "photogravimagnetic assists" of sailcraft around $\alpha$ Centauri A, a potential first destination en route to Proxima Centauri (the goal of the Breakthrough Starshot program). We find that a 10m$^2$ sail with a charge-to-mass-ratio of around 10 $\mu$C/g or higher will need to take account of magnetic field effects during orbital maneouvres. The magnetic field can provide an extra source of deceleration and deflection, and allow capture onto closer orbits around a target star. However, flipping the sign of the sailcraft's charge can radically change resulting trajectories, resulting in complex loop-de-loops around magnetic field lines and essentially random ejection from the star system. Even on well-behaved trajectories, the field can generate off-axis deflections at $\alpha$ Centauri that, while minor, can result in very poor targetting of the final destination (Proxima) post-assist. Fortunately for Breakthrough Starshot, nanosails are less prone to charging en route than their heavier counterparts, but can still accrue relatively high charge at both the origin and destination, when travelling at low speeds. Photogravimagnetic assists are highly non-trivial, and require careful course correction to mitigate against unwanted changes in trajectory.


title: An Optimal and Progressive Approach to Online Search of Top-k Influential Communities

abstract: Community search over large graphs is a fundamental problem in graph analysis. Recent studies propose to compute top-k influential communities, where each reported community not only is a cohesive subgraph but also has a high influence value. The existing approaches to the problem of top-k influential community search can be categorized as index-based algorithms and online search algorithms without indexes. The index-based algorithms, although being very efficient in conducting community searches, need to pre-compute a special-purpose index and only work for one built-in vertex weight vector. In this paper, we investigate on-line search approaches and propose an instance-optimal algorithm LocalSearch whose time complexity is linearly proportional to the size of the smallest subgraph that a correct algorithm needs to access without indexes. In addition, we also propose techniques to make LocalSearch progressively compute and report the communities in decreasing influence value order such that k does not need to be specified. Moreover, we extend our framework to the general case of top-k influential community search regarding other cohesiveness measures. Extensive empirical studies on real graphs demonstrate that our algorithms outperform the existing online search algorithms by several orders of magnitude.


title: End-to-end 3D shape inverse rendering of different classes of objects from a single input image

abstract: In this paper a semi-supervised deep framework is proposed for the problem of 3D shape inverse rendering from a single 2D input image. The main structure of proposed framework consists of unsupervised pre-trained components which significantly reduce the need to labeled data for training the whole framework. using labeled data has the advantage of achieving to accurate results without the need to predefined assumptions about image formation process. Three main components are used in the proposed network: an encoder which maps 2D input image to a representation space, a 3D decoder which decodes a representation to a 3D structure and a mapping component in order to map 2D to 3D representation. The only part that needs label for training is the mapping part with not too many parameters. The other components in the network can be pre-trained unsupervised using only 2D images or 3D data in each case. The way of reconstructing 3D shapes in the decoder component, inspired by the model based methods for 3D reconstruction, maps a low dimensional representation to 3D shape space with the advantage of extracting the basis vectors of shape space from training data itself and is not restricted to a small set of examples as used in predefined models. Therefore, the proposed framework deals directly with coordinate values of the point cloud representation which leads to achieve dense 3D shapes in the output. The experimental results on several benchmark datasets of objects and human faces and comparing with recent similar methods shows the power of proposed network in recovering more details from single 2D images.


title: Hybrid Approach of Relation Network and Localized Graph Convolutional Filtering for Breast Cancer Subtype Classification

abstract: Network biology has been successfully used to help reveal complex mechanisms of disease, especially cancer. On the other hand, network biology requires in-depth knowledge to construct disease-specific networks, but our current knowledge is very limited even with the recent advances in human cancer biology. Deep learning has shown a great potential to address the difficult situation like this. However, deep learning technologies conventionally use grid-like structured data, thus application of deep learning technologies to the classification of human disease subtypes is yet to be explored. Recently, graph based deep learning techniques have emerged, which becomes an opportunity to leverage analyses in network biology. In this paper, we proposed a hybrid model, which integrates two key components 1) graph convolution neural network (graph CNN) and 2) relation network (RN). We utilize graph CNN as a component to learn expression patterns of cooperative gene community, and RN as a component to learn associations between learned patterns. The proposed model is applied to the PAM50 breast cancer subtype classification task, the standard breast cancer subtype classification of clinical utility. In experiments of both subtype classification and patient survival analysis, our proposed method achieved significantly better performances than existing methods. We believe that this work is an important starting point to realize the upcoming personalized medicine.


title: A General Neural Network Hardware Architecture on FPGA

abstract: Field Programmable Gate Arrays (FPGAs) plays an increasingly important role in data sampling and processing industries due to its highly parallel architecture, low power consumption, and flexibility in custom algorithms. Especially, in the artificial intelligence field, for training and implement the neural networks and machine learning algorithms, high energy efficiency hardware implement and massively parallel computing capacity are heavily demanded. Therefore, many global companies have applied FPGAs into AI and Machine learning fields such as autonomous driving and Automatic Spoken Language Recognition (Baidu) [1] [2] and Bing search (Microsoft) [3]. Considering the FPGAs great potential in these fields, we tend to implement a general neural network hardware architecture on XILINX ZU9CG System On Chip (SOC) platform [4], which contains abundant hardware resource and powerful processing capacity. The general neural network architecture on the FPGA SOC platform can perform forward and backward algorithms in deep neural networks (DNN) with high performance and easily be adjusted according to the type and scale of the neural networks.


title: Modal Regression based Atomic Representation for Robust Face Recognition

abstract: Representation based classification (RC) methods such as sparse RC (SRC) have shown great potential in face recognition in recent years. Most previous RC methods are based on the conventional regression models, such as lasso regression, ridge regression or group lasso regression. These regression models essentially impose a predefined assumption on the distribution of the noise variable in the query sample, such as the Gaussian or Laplacian distribution. However, the complicated noises in practice may violate the assumptions and impede the performance of these RC methods. In this paper, we propose a modal regression based atomic representation and classification (MRARC) framework to alleviate such limitation. Unlike previous RC methods, the MRARC framework does not require the noise variable to follow any specific predefined distributions. This gives rise to the capability of MRARC in handling various complex noises in reality. Using MRARC as a general platform, we also develop four novel RC methods for unimodal and multimodal face recognition, respectively. In addition, we devise a general optimization algorithm for the unified MRARC framework based on the alternating direction method of multipliers (ADMM) and half-quadratic theory. The experiments on real-world data validate the efficacy of MRARC for robust face recognition.


title: Real-Time Document Image Classification using Deep CNN and Extreme Learning Machines

abstract: This paper presents an approach for real-time training and testing for document image classification. In production environments, it is crucial to perform accurate and (time-)efficient training. Existing deep learning approaches for classifying documents do not meet these requirements, as they require much time for training and fine-tuning the deep architectures. Motivated from Computer Vision, we propose a two-stage approach. The first stage trains a deep network that works as feature extractor and in the second stage, Extreme Learning Machines (ELMs) are used for classification. The proposed approach outperforms all previously reported structural and deep learning based methods with a final accuracy of 83.24% on Tobacco-3482 dataset, leading to a relative error reduction of 25% when compared to a previous Convolutional Neural Network (CNN) based approach (DeepDocClassifier). More importantly, the training time of the ELM is only 1.176 seconds and the overall prediction time for 2,482 images is 3.066 seconds. As such, this novel approach makes deep learning-based document classification suitable for large-scale real-time applications.


title: Categorical data analysis using a skewed Weibull regression model

abstract: In this paper, we present a Weibull link (skewed) model for categorical response data arising from binomial as well as multinomial model. We show that, for such types of categorical data, the most commonly used models (logit, probit and complementary log-log) can be obtained as limiting cases. We further compare the proposed model with some other asymmetrical models. The Bayesian as well as frequentist estimation procedures for binomial and multinomial data responses are presented in details. The analysis of two data sets to show the efficiency of the proposed model is performed.


title: Hidden Markov Random Field Iterative Closest Point

abstract: When registering point clouds resolved from an underlying 2-D pixel structure, such as those resulting from structured light and flash LiDAR sensors, or stereo reconstruction, it is expected that some points in one cloud do not have corresponding points in the other cloud, and that these would occur together, such as along an edge of the depth map. In this work, a hidden Markov random field model is used to capture this prior within the framework of the iterative closest point algorithm. The EM algorithm is used to estimate the distribution parameters and the hidden component memberships. Experiments are presented demonstrating that this method outperforms several other outlier rejection methods when the point clouds have low or moderate overlap.


title: Pricing Football Players using Neural Networks

abstract: We designed a multilayer perceptron neural network to predict the price of a football (soccer) player using data on more than 15,000 players from the football simulation video game FIFA 2017. The network was optimized by experimenting with different activation functions, number of neurons and layers, learning rate and its decay, Nesterov momentum based stochastic gradient descent, L2 regularization, and early stopping. Simultaneous exploration of various aspects of neural network training is performed and their trade-offs are investigated. Our final model achieves a top-5 accuracy of 87.2% among 119 pricing categories, and places any footballer within 6.32% of his actual price on average.


title: Fast and Efficient Calculations of Structural Invariants of Chirality

abstract: Chirality plays an important role in physics, chemistry, biology, and other fields. It describes an essential symmetry in structure. However, chirality invariants are usually complicated in expression or difficult to evaluate. In this paper, we present five general three-dimensional chirality invariants based on the generating functions. And the five chiral invariants have four characteristics:(1) They play an important role in the detection of symmetry, especially in the treatment of 'false zero' problem. (2) Three of the five chiral invariants decode an universal chirality index. (3) Three of them are proposed for the first time. (4) The five chiral invariants have low order no bigger than 4, brief expression, low time complexity O(n) and can act as descriptors of three-dimensional objects in shape analysis. The five chiral invariants give a geometric view to study the chiral invariants. And the experiments show that the five chirality invariants are effective and efficient, they can be used as a tool for symmetry detection or features in shape analysis.


title: Geburtstage, Wuerfel, Produkte und Karten

abstract: This article, based on a talk, treats some elementary, but not completely simple examples from probability. They concern multiple birthday coincidences, throwing dice, the combinatorics of the German card game "Doppelkopf", and the properties of products of uniformly distributed random numbers. The material, a lot of which was not taken from or found in other sources, should be of interest to all those who lecture or plan to lecture on probability, especially on an elementary or introductory level, and are looking for challenging examples or problems.


title: A complex path around the sign problem

abstract: We review recent attempts at dealing with the sign problem in Monte Carlo calculations by deforming the region of integration in the path integral from real to complex fields. We discuss the theoretical foundations, the algorithmic issues and present some results for low dimensional field theories in both imaginary and real time.


title: Predictive Independence Testing, Predictive Conditional Independence Testing, and Predictive Graphical Modelling

abstract: Testing (conditional) independence of multivariate random variables is a task central to statistical inference and modelling in general - though unfortunately one for which to date there does not exist a practicable workflow. State-of-art workflows suffer from the need for heuristic or subjective manual choices, high computational complexity, or strong parametric assumptions. We address these problems by establishing a theoretical link between multivariate/conditional independence testing, and model comparison in the multivariate predictive modelling aka supervised learning task. This link allows advances in the extensively studied supervised learning workflow to be directly transferred to independence testing workflows - including automated tuning of machine learning type which addresses the need for a heuristic choice, the ability to quantitatively trade-off computational demand with accuracy, and the modern black-box philosophy for checking and interfacing. As a practical implementation of this link between the two workflows, we present a python package 'pcit', which implements our novel multivariate and conditional independence tests, interfacing the supervised learning API of the scikit-learn package. Theory and package also allow for straightforward independence test based learning of graphical model structure. We empirically show that our proposed predictive independence test outperform or are on par to current practice, and the derived graphical model structure learning algorithms asymptotically recover the 'true' graph. This paper, and the 'pcit' package accompanying it, thus provide powerful, scalable, generalizable, and easy-to-use methods for multivariate and conditional independence testing, as well as for graphical model structure learning.


title: Large time behavior of entropy solutions to one-dimensional unipolar hydrodynamic model for semiconductor devices

abstract: We are concerned with the global existence and large time behavior of entropy solutions to the one dimensional unipolar hydrodynamic model for semiconductors in the form of Euler-Poisson equations in a bounded interval. In this paper, we first prove the global existence of entropy solution by vanishing viscosity and compensated compactness framework. In particular, the solutions are uniformly bounded with respect to space and time variables by introducing modified Riemann invariants and the theory of invariant region. Based on the uniform estimates of density, we further show that the entropy solution converges to the corresponding unique stationary solution exponentially in time. No any smallness condition is assumed on the initial data and doping profile. Moreover, the novelty in this paper is about the unform bound with respect to time for the weak solutions of the isentropic Euler-Possion system.


title: Evidence for quantized vortex entry in Nb3Sn under the influence of RF fields

abstract: Evidence for quantized magnetic vortex entry at a defect into a Type-II superconductor operating in the metastable state above Bc1 while under the influence of RF electromagnetic fields has been observed in a polycrystalline coating of Nb3Sn. The measurement was made using a 1.3 GHz microwave cavity for use in particle accelerator R&D and equipped with a temperature mapping system for precision measurements of local heating. At peak RF fields near the quench field of the cavity, discrete jumps in local heating at the quench origin were observed. A number of these jumps were found to be identical in height within measurement error. This suggests entry of a limited number of vortices at a defect. The behavior observed appears to be consistent with recent theoretical predictions regarding vortex entry at grain boundaries in the material.


title: Calibration of atomic trajectories in a large-area dual-atom-interferometer gyroscope

abstract: We propose and demonstrate a method for calibrating atomic trajectories in a large-area dual-atom-interferometer gyroscope. The atom trajectories are monitored by modulating and delaying the Raman transition, and they are precisely calibrated by controlling the laser orientation and the bias magnetic field. To improve the immunity to the gravity effect and the common phase noise, the symmetry and the overlapping of two large-area atomic interference loops are optimized by calibrating the atomic trajectories and by aligning the Raman-laser orientations. The dual-atom-interferometer gyroscope is applied in the measurement of the Earth rotation. The sensitivity is $1.2\times10^{-6}$ rad/s/$\sqrt{Hz}$, and the long-term stability is $6.2\times10^{-8}$ rad/s $@$ $2000$ s.


title: Pure phase BiFeO$_3$ thin films sputtered over Si: A new route towards high magnetization

abstract: We have investigated the structural and magnetic properties of BiFeO$_3$ (BFO) thin films grown over (100)-oriented Si substrates by rf magnetron sputtering in a new route under O$_2$ free low pressure Ar atmosphere. Single-phase BFO films were deposited in a heated substrate and post-annealed in situ. The new routed allows high deposition rate and produce polycrystalline BFO pure phase, confirmed by high resolution X-ray diffraction. Scanning electron and atomic force microscopy reveal very low surface roughness and mean particle size of 33 nm. The BFO phase and composition were confirmed by transmission electron microscopy and line scanning energy-dispersive X-ray spectroscopy in transmission electron microscopy mode. The surface chemistry of the thin film, analyzed by X-ray photoelectron spectroscopy, reveals the presence of Fe$^{3+}$ and Fe$^{2+}$ in a 2:1 ratio, a strong indication that the film contains oxygen vacancies. An hysteretic ferromagnetic behavior with room temperature high saturation magnetization $\sim 165 \times 10^3$ A/m was measured along the film perpendicular and parallel directions. Such high magnetization, deriving from this new route, is explained in the scope of oxygen vacancies, the break of the antiferromagnetic cycloidal order and the increase of spin canting by change in the surface/volume ratio. Understanding the magnetic behavior of a multiferroic thin films is a key for the development of heterogeneous layered structures and multilayered devices and the production of multiferroic materials over Si substrates opens new possibilities in the development of materials that can be directly integrated into the existent semiconductor and spintronic technologies.


title: A new characterization of the dual polar graphs

abstract: In this paper we give a new characterization of the dual polar graphs, extending the work of Brouwer and Wilbrink on regular near polygons. Also as a consequence of our characterization we confirm a conjecture of the authors on non-bipartite distance-regular graphs with smallest eigenvalue at most $-k/2$, where $k$ is the valency of the distance-regular graph, in case of $c_2 \geq3$ and $a_1 =1$.


title: Scaling Trapped Ion Quantum Computers Using Fast Gates and Microtraps

abstract: Most attempts to produce a scalable quantum information processing platform based on ion traps have focused on the shuttling of ions in segmented traps. We show that an architecture based on an array of microtraps with fast gates will outperform architectures based on ion shuttling. This system requires higher power lasers, but does not require the manipulation of potentials or shuttling of ions. This improves optical access, reduces the complexity of the trap, and reduces the number of conductive surfaces close to the ions. The use of fast gates also removes limitations on gate time. The performance of the gates is shown to be robust to the limitations in laser repetition rate and the presence of many ions in the trap array.


title: Extreme Quiescent Variability of the Transient Neutron Star Low-mass X-ray Binary EXO 1745-248 in Terzan 5

abstract: EXO 1745-248 is a transient neutron-star low-mass X-ray binary that resides in the globular cluster Terzan 5. We studied the transient during its quiescent state using 18 Chandra observations of the cluster acquired between 2003 and 2016. We found an extremely variable source, with a luminosity variation in the 0.5-10 keV energy range of $\sim3$ orders of magnitude (between $3\times10^{31}$ erg s$^{-1}$ and $2\times10^{34}$ erg s$^{-1}$) on timescales from years down to only a few days. Using an absorbed power-law model to fit its quiescent spectra, we obtained a typical photon index of $\sim1.4$, indicating that the source is even harder than during outburst and much harder than typical quiescent neutron stars if their quiescent X-ray spectra are also described by a single power-law model. This indicates that EXO 1745-248 is very hard throughout the entire observed X-ray luminosity range. At the highest luminosity, the spectrum fits better when an additional (soft) component is added to the model. All these quiescent properties are likely related to strong variability in the low-level accretion rate in the system. However, its extreme variable behavior is strikingly different from the one observed for other neutron star transients that are thought to still accrete in quiescence. We compare our results to these systems. We also discuss similarities and differences between our target and the transitional millisecond pulsar IGR J18245-2452, which also has hard spectra and strong variability during quiescence.


title: Packing nearly optimal Ramsey R(3,t) graphs

abstract: In 1995 Kim famously proved the Ramsey bound R(3,t) \ge c t^2/\log t by constructing an n-vertex graph that is triangle-free and has independence number at most C \sqrt{n \log n}. We extend this celebrated result, which is best possible up to the value of the constants, by approximately decomposing the complete graph K_n into a packing of such nearly optimal Ramsey R(3,t) graphs. More precisely, for any \epsilon>0 we find an edge-disjoint collection (G_i)_i of n-vertex graphs G_i \subseteq K_n such that (a) each G_i is triangle-free and has independence number at most C_\epsilon \sqrt{n \log n}, and (b) the union of all the G_i contains at least (1-\epsilon)\binom{n}{2} edges. Our algorithmic proof proceeds by sequentially choosing the graphs G_i via a semi-random (i.e., Rodl nibble type) variation of the triangle-free process. As an application, we prove a conjecture in Ramsey theory by Fox, Grinshpun, Liebenau, Person, and Szabo (concerning a Ramsey-type parameter introduced by Burr, Erdos, Lovasz in 1976). Namely, denoting by s_r(H) the smallest minimum degree of r-Ramsey minimal graphs for H, we close the existing logarithmic gap for H=K_3 and establish that s_r(K_3) = \Theta(r^2 \log r).


title: Efficient D-optimal design of experiments for infinite-dimensional Bayesian linear inverse problems

abstract: We develop a computational framework for D-optimal experimental design for PDE-based Bayesian linear inverse problems with infinite-dimensional parameters. We follow a formulation of the experimental design problem that remains valid in the infinite-dimensional limit. The optimal design is obtained by solving an optimization problem that involves repeated evaluation of the log-determinant of high-dimensional operators along with their derivatives. Forming and manipulating these operators is computationally prohibitive for large-scale problems. Our methods exploit the low-rank structure in the inverse problem in three different ways, yielding efficient algorithms. Our main approach is to use randomized estimators for computing the D-optimal criterion, its derivative, as well as the Kullback--Leibler divergence from posterior to prior. Two other alternatives are proposed based on a low-rank approximation of the prior-preconditioned data misfit Hessian, and a fixed low-rank approximation of the prior-preconditioned forward operator. Detailed error analysis is provided for each of the methods, and their effectiveness is demonstrated on a model sensor placement problem for initial state reconstruction in a time-dependent advection-diffusion equation in two space dimensions.


title: (geo)graphs - Complex Networks as a shapefile of nodes and a shapefile of edges for different applications

abstract: Spatial dependency and spatial embedding are basic physical properties of many phenomena modeled by networks. The most indicated computational environment to deal with spatial information is to use Georeferenced Information System (GIS) and Geographical Database Management Systems (GDBMS). Several models have been proposed in this direction, however there is a gap in the literature in generic frameworks for working with Complex Networks in GIS/GDBMS environments. Here we introduce the concept of (geo)graphs: graphs in which the nodes have a known geographical location and the edges have spatial dependence. We present case studies and two open source softwares (GIS4GRAPH and GeoCNet) that indicate how to retrieve networks from GIS data and how to represent networks over GIS data by using (geo)graphs.


title: Convergence of iterative methods based on Neumann series for composite materials: theory and practice

abstract: Iterative Fast Fourier Transform methods are useful for calculating the fields in composite materials and their macroscopic response. By iterating back and forth until convergence, the differential constraints are satisfied in Fourier space, and the constitutive law in real space. The methods correspond to series expansions of appropriate operators and to series expansions for the effective tensor as a function of the component moduli. It is shown that the singularity structure of this function can shed much light on the convergence properties of the iterative Fast Fourier Transform methods. We look at a model example of a square array of conducting square inclusions for which there is an exact formula for the effective conductivity (Obnosov). Theoretically some of the methods converge when the inclusions have zero or even negative conductivity. However, the numerics do not always confirm this extended range of convergence and show that accuracy is lost after relatively few iterations. There is little point in iterating beyond this. Accuracy improves when the grid size is reduced, showing that the discrepancy is linked to the discretization. Finally, it is shown that none of the three iterative schemes investigated over-performs the others for all possible microstructures and all contrasts.


title: Faster Methods for Contracting Infinite 2D Tensor Networks

abstract: We revisit the corner transfer matrix renormalization group (CTMRG) method of Nishino and Okunishi for contracting two-dimensional (2D) tensor networks and demonstrate that its performance can be substantially improved by determining the tensors using an eigenvalue solver as opposed to the power method used in CTMRG. We also generalize the variational uniform matrix product state (VUMPS) ansatz for diagonalizing 1D quantum Hamiltonians to the case of 2D transfer matrices and discuss similarities with the corner methods. These two new algorithms will be crucial to improving the performance of variational infinite projected entangled pair state (PEPS) methods.


title: Solution Uniqueness of Convex Piecewise Affine Functions Based Optimization with Applications to Constrained $\ell_1$ Minimization

abstract: In this paper, we study the solution uniqueness of an individual feasible vector of a class of convex optimization problems involving convex piecewise affine functions and subject to general polyhedral constraints. This class of problems incorporates many important polyhedral constrained $\ell_1$ recovery problems arising from sparse optimization, such as basis pursuit, LASSO, and basis pursuit denoising, as well as polyhedral gauge recovery. By leveraging the max-formulation of convex piecewise affine functions and convex analysis tools, we develop dual variables based necessary and sufficient uniqueness conditions via simple and yet unifying approaches; these conditions are applied to a wide range of $\ell_1$ minimization problems under possible polyhedral constraints. An effective linear program based scheme is proposed to verify solution uniqueness conditions. The results obtained in this paper not only recover the known solution uniqueness conditions in the literature by removing restrictive assumptions but also yield new uniqueness conditions for much broader constrained $\ell_1$-minimization problems.


title: Hybrid Normed Ideal Perturbations of n-tuples of Operators I

abstract: In hybrid normed ideal perturbations of $n$-tuples of operators, the normed ideal is allowed to vary with the component operators. We begin extending to this setting the machinery we developed for normed ideal perturbations based on the modulus of quasicentral approximation and an adaptation of our non-commutative generalization of the Weyl--von~Neumann theorem. For commuting $n$-tuples of hermitian operators, the modulus of quasicentral approximation remains essentially the same when $\cC_n^-$ is replaced by a hybrid $n$-tuple $\cC_{p_1,\dots}^-,\dots,\cC^-_{p_n}$, $p_1^{-1} + \dots + p_n^{-1} = 1$. The proof involves singular integrals of mixed homogeneity.


title: Complex Dynamics of Bus, Tram and Elevator Delays in Transportation System

abstract: It is necessary and important to operate buses and trams on time. The bus schedule is closely related to the dynamic motion of buses. In this part, we introduce the nonlinear maps for describing the dynamics of shuttle buses in the transportation system. The complex motion of the buses is explained by the nonlinear-map models. The transportation system of shuttle buses without passing is similar to that of the trams. The transport of elevators is also similar to that of shuttle buses with freely passing. The complex dynamics of a single bus is described in terms of the piecewise map, the delayed map, the extended circle map and the combined map. The dynamics of a few buses is described by the model of freely passing buses, the model of no passing buses, and the model of increase or decrease of buses. The nonlinear-map models are useful to make an accurate estimate of the arrival time in the bus transportation.