title: A reformulated Krein matrix for star-even polynomial operators with applications

abstract: In its original formulation the Krein matrix was used to locate the spectrum of first-order star-even polynomial operators where both operator coefficients are nonsingular. Such operators naturally arise when considering first-order-in-time Hamiltonian PDEs. Herein the matrix is reformulated to allow for operator coefficients with nontrivial kernel. Moreover, it is extended to allow for the study of the spectral problem associated with quadratic star-even operators, which arise when considering the spectral problem associated with second-order-in-time Hamiltonian PDEs. In conjunction with the Hamiltonian-Krein index (HKI) the Krein matrix is used to study two problems: conditions leading to Hamiltonian-Hopf bifurcations for small spatially periodic waves, and the location and Krein signature of small eigenvalues associated with, e.g., $n$-pulse problems. For the first case we consider in detail a first-order-in-time fifth-order KdV-like equation. In the latter case we use a combination of Lin's method, the HKI, and the Krein matrix to study the spectrum associated with $n$-pulses for a second-order-in-time Hamiltonian system which is used to model the dynamics of a suspension bridge.


title: A time-reversible integrator for the time-dependent Schr\"{o}dinger equation on an adaptive grid

abstract: One of the most accurate methods for solving the time-dependent Schr\"{o}dinger equation uses a combination of the dynamic Fourier method with the split-operator algorithm on a tensor-product grid. To reduce the number of required grid points, we let the grid move together with the wavepacket, but find that the na\"ive algorithm based on an alternate evolution of the wavefunction and grid destroys the time reversibility of the exact evolution. Yet, we show that the time reversibility is recovered if the wavefunction and grid are evolved simultaneously during each kinetic or potential step; this is achieved by using the Ehrenfest theorem together with the splitting method. The proposed algorithm is conditionally stable, symmetric, time-reversible, and conserves the norm of the wavefunction. The preservation of these geometric properties is shown analytically and demonstrated numerically on a three-dimensional harmonic model and collinear model of He-H$_{2}$ scattering. We also show that the proposed algorithm can be symmetrically composed to obtain time-reversible integrators of an arbitrary even order. We observed $10000$-fold speedup by using the tenth- instead of the second- order method to obtain a solution with a time discretization error below $10^{-9}$. Moreover, using the adaptive grid instead of the fixed grid resulted in a 64-fold reduction in the required number of grid points in the harmonic system and made it possible to simulate the He-H$_{2}$ scattering for six times longer, while maintaining reasonable accuracy. Applicability of the algorithm to high-dimensional quantum dynamics is demonstrated using the strongly anharmonic eight-dimensional H\'{e}non--Heiles model.


title: Dynamic Graph Algorithms and Graph Sparsification: New Techniques and Connections

abstract: Graphs naturally appear in several real-world contexts including social networks, the web network, and telecommunication networks. While the analysis and the understanding of graph structures have been a central area of study in algorithm design, the rapid increase of data sets over the last decades has posed new challenges for designing efficient algorithms that process large-scale graphs. These challenges arise from two usual assumptions in classical algorithm design, namely that graphs are static and that they fit into a single machine. However, in many application domains, graphs are subject to frequent changes over time, and their massive size makes them infeasible to be stored in the memory of a single machine. Driven by the need to devise new tools for overcoming such challenges, this thesis focuses on two areas of modern algorithm design that directly deal with processing massive graphs, namely dynamic graph algorithms and graph sparsification. We develop new algorithmic techniques from both dynamic and sparsification perspective for a multitude of graph-based optimization problems which lie at the core of Spectral Graph Theory, Graph Partitioning, and Metric Embeddings. Our algorithms are faster than any previous one and design smaller sparsifiers with better (approximation) quality. More importantly, this work introduces novel reduction techniques that show unexpected connections between seemingly different areas such as dynamic graph algorithms and graph sparsification.


title: Learning Household Task Knowledge from WikiHow Descriptions

abstract: Commonsense procedural knowledge is important for AI agents and robots that operate in a human environment. While previous attempts at constructing procedural knowledge are mostly rule- and template-based, recent advances in deep learning provide the possibility of acquiring such knowledge directly from natural language sources. As a first step in this direction, we propose a model to learn embeddings for tasks, as well as the individual steps that need to be taken to solve them, based on WikiHow articles. We learn these embeddings such that they are predictive of both step relevance and step ordering. We also experiment with the use of integer programming for inferring consistent global step orderings from noisy pairwise predictions.


title: Enabling Intuitive Human-Robot Teaming Using Augmented Reality and Gesture Control

abstract: Human-robot teaming offers great potential because of the opportunities to combine strengths of heterogeneous agents. However, one of the critical challenges in realizing an effective human-robot team is efficient information exchange - both from the human to the robot as well as from the robot to the human. In this work, we present and analyze an augmented reality-enabled, gesture-based system that supports intuitive human-robot teaming through improved information exchange. Our proposed system requires no external instrumentation aside from human-wearable devices and shows promise of real-world applicability for service-oriented missions. Additionally, we present preliminary results from a pilot study with human participants, and highlight lessons learned and open research questions that may help direct future development, fielding, and experimentation of autonomous HRI systems.


title: The bilinear Hilbert transform in UMD spaces

abstract: We prove $L^p$-bounds for the bilinear Hilbert transform acting on functions valued in intermediate UMD spaces. Such bounds were previously unknown for UMD spaces that are not Banach lattices. Our proof relies on bounds on embeddings from Bochner spaces $L^p(\mathbb{R};X)$ into outer Lebesgue spaces on the time-frequency-scale space $\mathbb{R}^3_+$.


title: Multi-Particle Collisions in Microgravity: Coefficient of Restitution and Sticking Threshold for Systems of Mm-Sized Particles

abstract: The current model of planet formation lacks a good understanding of the growth of dust particles inside the protoplanetary disk beyond mm sizes. In order to investigate the low-velocity collisions between this type of particles, the NanoRocks experiment was flown on the International Space Station (ISS) between September 2014 and March 2016. We present the results of this experiment. We quantify the damping of energy in systems of multiple particles in the 0.1 to 1 mm size range while they are in the bouncing regime, and study the formation of clusters through sticking collisions between particles. We developed statistical methods for the analysis of the large quantity of collision data collected by the experiment. We measured the average motion of particles, the moment of clustering, and the cluster size formed. In addition, we ran simple numerical simulations in order to validate our measurements. We computed the average coefficient of restitution (COR) of collisions and find values ranging from 0.55 for systems including a population of fine grains to 0.94 for systems of denser particles. We also measured the sticking threshold velocities and find values around 1 cm/s, consistent with the current dust collision models based on independently collected experimental data. Our findings have the following implications that can be useful for the simulation of particles in PPDs and planetary rings: (1) The average COR of collisions between same-sized free-floating particles at low speeds (< 2 cm/s) is not dependent on the collision velocity; (2) The simplified approach of using a constant COR value will accurately reproduce the average behavior of a particle system during collisional cooling; (3) At speeds below 5 mm/s, the influence of particle rotation becomes apparent on the collision behavior; (4) Current dust collision models predicting sticking thresholds are robust.


title: Towards A Robot Explanation System: A Survey and Our Approach to State Summarization, Storage and Querying, and Human Interface

abstract: As robot systems become more ubiquitous, developing understandable robot systems becomes increasingly important in order to build trust. In this paper, we present an approach to developing a holistic robot explanation system, which consists of three interconnected components: state summarization, storage and querying, and human interface. To find trends towards and gaps in the development of such an integrated system, a literature review was performed and categorized around those three components, with a focus on robotics applications. After the review of each component, we discuss our proposed approach for robot explanation. Finally, we summarize the system as a whole and review its functionality.


title: Null hypersurfaces in Kerr-(A)dS spacetimes

abstract: For many purposes, a three-dimensional foliation of spacetime is more advantageous to understanding its light cone structure. We derive the equations describing such foliations for the Kerr geometry with non-zero cosmological constant, and show that they reduce to null hypersurfaces in vacuum (anti-)de Sitter spacetime in the limit of zero mass. Furthermore, we find that these null hypersurfaces are free of caustics everywhere for $r>0$. Our construction has applications in numerical studies of rotating black holes, and in defining Kruskal coordinates for rotating black holes with non-zero cosmological constant.


title: Controlling a Random Population is EXPTIME-hard

abstract: Bertrand et al. [1] (LMCS 2019) describe two-player zero-sum games in which one player tries to achieve a reachability objective in $n$ games (on the same finite arena) simultaneously by broadcasting actions, and where the opponent has full control of resolving non-deterministic choices. They show EXPTIME completeness for the question if such games can be won for every number $n$ of games. We consider the almost-sure variant in which the opponent randomizes their actions, and where the player tries to achieve the reachability objective eventually with probability one. The lower bound construction in [1] does not directly carry over to this randomized setting. In this note we show EXPTIME hardness for the almost-sure problem by reduction from Countdown Games.


title: Degenerate elastic networks

abstract: We minimize a linear combination of the Willmore and the length functional among networks in $\mathbb{R}^d$ belonging to a given class determined by the number of curves, the order of the junctions and the angles between curves at the junctions. Since this class lacks compactness, we characterize the set of limits of sequences of networks bounded in energy, providing an explicit representation of the relaxed problem. This is expressed in terms of the new notion of degenerate elastic networks that, rather surprisingly, involves only the properties of the given class, without reference to the curvature. In the case of $d=2$ we also give an equivalent description of degenerate elastic networks by means of a combinatorial definition easy to validate by a finite algorithm. Moreover we provide examples, counterexamples, and additional results that motivate our study and show the sharpness of our characterization.


title: Uniqueness and nonuniqueness of limits of Teichmueller harmonic map flow

abstract: The harmonic map energy of a map from a closed, constant-curvature surface to a closed target manifold can be seen as a functional on the space of maps and domain metrics. We consider the gradient flow for this energy. In the absence of singularities, previous theory established that the flow converges to a branched minimal immersion, but only at a sequence of times converging to infinity, and only after pulling back by a sequence of diffeomorphisms. In this paper we investigate whether it is necessary to pull back by these diffeomorphisms, and whether the convergence is uniform as time tends to infinity.


title: Zero-Shot Action Recognition in Videos: A Survey

abstract: Zero-Shot Action Recognition has attracted attention in the last years and many approaches have been proposed for recognition of objects, events and actions in images and videos. There is a demand for methods that can classify instances from classes that are not present in the training of models, especially in the complex problem of automatic video understanding, since collecting, annotating and labeling videos are difficult and laborious tasks. We have identified that there are many methods available in the literature, however, it is difficult to categorize which techniques can be considered state of the art. Despite the existence of some surveys about zero-shot action recognition in still images and experimental protocol, there is no work focused on videos. Therefore, we present a survey of the methods that comprise techniques to perform visual feature extraction and semantic feature extraction as well to learn the mapping between these features considering specifically zero-shot action recognition in videos. We also provide a complete description of datasets, experiments and protocols, presenting open issues and directions for future work, essential for the development of the computer vision research field.


title: Trispectrum reconstruction of non-Gaussian noise

abstract: Using a qubit to probe non-Gaussian noise environments is theoretically studied in the context of classical random telegraph processes. Protocols for control pulses are developed to effectively scan higher noise correlations, offering valuable information on the charge environment of the qubit. Specifically, the noise power spectrum and trispectrum are reconstructed simultaneously for a wide range of qubit-fluctuator coupling strengths, demonstrating the method's robustness. These protocols are readily testable in various qubit systems with well-developed quantum control, including quantum dot spins, superconducting qubits and NV centers in diamond.


title: Compositional Synthesis of Decentralized Robust Set-Invariance Controllers for Large-scale Linear Systems

abstract: Ensuring constraint satisfaction in large-scale systems with hard constraints is vital in many safety critical systems. The challenge is to design controllers that are efficiently synthesized offline, easily implementable online, and provide formal correctness guarantees. In this paper, we provide a method to compute correct-by-construction controllers for a network of coupled linear systems with additive bounded disturbances such that i) the design of the controllers is fully compositional - we use an optimization-based approach that iteratively computes subsystem-level assume-guarantee contracts in the form of robust control invariant sets; and ii) the controllers are decentralized hence online implementation requires only the local state information. We present illustrative examples, including a case study on a system with 1000 dimensions.


title: Bannai-Ito algebras and the universal R-matrix of osp(1|2)

abstract: The Bannai-Ito algebra $BI(n)$ is viewed as the centralizer of the action of $\mathfrak{osp}(1|2)$ in the $n$-fold tensor product of the universal algebra of this Lie superalgebra. The generators of this centralizer are constructed with the help of the universal $R$-matrix of $\mathfrak{osp}(1|2)$. The specific structure of the $\mathfrak{osp}(1|2)$ embeddings to which the centralizing elements are attached as Casimir elements is explained. With the generators defined, the structure relations of $BI(n)$ are derived from those of $BI(3)$ by repeated action of the coproduct and using properties of the $R$-matrix and of the generators of the symmetric group $\mathfrak S_n$.


title: Responsive Planning and Recognition for Closed-Loop Interaction

abstract: Many intelligent systems currently interact with others using at least one of fixed communication inputs or preset responses, resulting in rigid interaction experiences and extensive efforts developing a variety of scenarios for the system. Fixed inputs limit the natural behavior of the user in order to effectively communicate, and preset responses prevent the system from adapting to the current situation unless it was specifically implemented. Closed-loop interaction instead focuses on dynamic responses that account for what the user is currently doing based on interpretations of their perceived activity. Agents employing closed-loop interaction can also monitor their interactions to ensure that the user responds as expected. We introduce a closed-loop interactive agent framework that integrates planning and recognition to predict what the user is trying to accomplish and autonomously decide on actions to take in response to these predictions. Based on a recent demonstration of such an assistive interactive agent in a turn-based simulated game, we also discuss new research challenges that are not present in the areas of artificial intelligence planning or recognition alone.


title: Coproducts of proximity spaces

abstract: In this paper, we introduce coproducts of proximity spaces. After exploring several of their basic properties, we show that given a collection of proximity spaces, the coproduct of their Smirnov compactifications proximally and densely embeds in the Smirnov compactification of the coproduct of the original proximity spaces. We also show that the dense proximity embedding is a proximity isomorphism if and only if the index set is finite. After constructing a number of examples of coproducts and their Smirnov compactifications, we explore several properties of the Smirnov compactification of the coproduct, including its metrizability, connectedness of the boundary, dimension, and its relation to the Stone-Cech compactification. In particular, we show that the Smirnov compactification of the infinite coproduct is never metrizable and that its boundary is highly disconnected. We also show that the proximity dimension of the Smirnov compactification of the coproduct equals the supremum of the covering dimensions of the individual Smirnov compactifications and that the Smirnov compactification of the coproduct is homeomorphic to the Stone-Cech compactification if and only if each individual proximity space is equipped with the Stone-Cech proximity. We finish with an example of a coproduct with the covering dimension $0$ but the proximity dimension $\infty.$


title: Recommendation or Discrimination?: Quantifying Distribution Parity in Information Retrieval Systems

abstract: Information retrieval (IR) systems often leverage query data to suggest relevant items to users. This introduces the possibility of unfairness if the query (i.e., input) and the resulting recommendations unintentionally correlate with latent factors that are protected variables (e.g., race, gender, and age). For instance, a visual search system for fashion recommendations may pick up on features of the human models rather than fashion garments when generating recommendations. In this work, we introduce a statistical test for "distribution parity" in the top-K IR results, which assesses whether a given set of recommendations is fair with respect to a specific protected variable. We evaluate our test using both simulated and empirical results. First, using artificially biased recommendations, we demonstrate the trade-off between statistically detectable bias and the size of the search catalog. Second, we apply our test to a visual search system for fashion garments, specifically testing for recommendation bias based on the skin tone of fashion models. Our distribution parity test can help ensure that IR systems' results are fair and produce a good experience for all users.


title: LDPC Codes Achieve List Decoding Capacity

abstract: We show that Gallager's ensemble of Low-Density Parity Check (LDPC) codes achieves list-decoding capacity with high probability. These are the first graph-based codes shown to have this property. This result opens up a potential avenue towards truly linear-time list-decodable codes that achieve list-decoding capacity. Our result on list decoding follows from a much more general result: any $\textit{local}$ property satisfied with high probability by a random linear code is also satisfied with high probability by a random LDPC code from Gallager's distribution. Local properties are properties characterized by the exclusion of small sets of codewords, and include list-decodability, list-recoverability and average-radius list-decodability. In order to prove our results on LDPC codes, we establish sharp thresholds for when local properties are satisfied by a random linear code. More precisely, we show that for any local property $\mathcal{P}$, there is some $R^*$ so that random linear codes of rate slightly less than $R^*$ satisfy $\mathcal{P}$ with high probability, while random linear codes of rate slightly more than $R^*$, with high probability, do not. We also give a characterization of the threshold rate $R^*$.


title: Lubricated friction of textured soft substrates

abstract: The understanding of sliding friction for wet, patterned surfaces from first principles is challenging. While emerging applications have sought design principles from biology, a general framework is lacking because soft interfaces experience a multiphysics coupling between solid deformation and fluid dissipation. We investigate the elastohydrodynamic sliding of >50 patterned sliding pairs comprising elastomers, thermosets, and hydrogels, and discover that texturing induces a critical transition in the macroscopic friction coefficient. This critical friction scales universally, without any fitting parameters, with the reduced elastic modulus and the pattern geometry. To capture the frictional dissipation, we separate the flow curve into two regimes and account for the contributions of shear and normal forces applied by the fluid on the patterns. Our model combines Reynolds' equations and elastic deformation to provide physical insights that allow engineering of the elastohydrodynamic friction in a class of soft tribopairs using pattern geometry, material elasticity, and fluid properties.


title: Bayesian analysis of longitudinal studies with treatment by indication

abstract: It is often of interest in observational studies to measure the causal effect of a treatment on time-to-event outcomes. In a medical setting, observational studies commonly involve patients who initiate medication therapy and others who do not, and the goal is to infer the effect of medication therapy on time until recovery, a pre-defined level of improvement, or some other time-to-event outcome. A difficulty with such studies is that the notion of a medication initiation time does not exist in the control group. We propose an approach to infer causal effects of an intervention in longitudinal observational studies when the time of treatment assignment is only observed for treated units and where treatment is given by indication. We present a framework for conceptualizing an underlying randomized experiment in this setting based on separating the process that governs the time of study arm assignment from the mechanism that determines the assignment. Our approach involves inferring the missing times of assignment followed by estimating treatment effects. This approach allows us to incorporate uncertainty about the missing times of study arm assignment, which induces uncertainty in both the selection of the control group and the measurement of time-to-event outcomes for these controls. We demonstrate our approach to study the effects on mortality of inappropriately prescribing phosphodiesterase type 5 inhibitors (PDE5Is), a medication contraindicated for groups 2 and 3 pulmonary hypertension, using administrative data from the Veterans Affairs (VA) health care system.


title: Strong Gravitational Lensing for Photon Coupled to Weyl Tensor in Kiselev Black Hole

abstract: The ambition of the present work is to highlight the phenomena of strong gravitational lensing and deflection angle for the photons coupling with Weyl tensor in a Kiselev black hole. Here, we have extended the prior work of Chen and Jing \cite{1} for Schwarzschild black hole to Kiselev black hole. For this purpose, the equation of motion for the photons coupled to Weyl tensor, null geodesic and equation of photon sphere in a Kiselev black hole spacetime have been formulated. It is found that the equation of motion of the photons depends not only on the coupling between photon and Weyl tensor, but also on the polarization direction of the photons. There is a critical value of the coupling parameter $\alpha$ for existence of the marginally circular photon orbit outside the event horizon, which depends on the parameters of black hole and the polarization direction of photons. Further, the polarization directions of coupled photon and the coupling parameter $\alpha$, both modify the features of the photon sphere, the angle of deflection and the functions $(\bar{a}$ and $\bar{b})$ for the strong gravitational lensing in Kiselev black hole spacetime. In addition to this, the observable gravitational lensing quantities and the shadows of the Kiselev black hole spacetime are presented in detail.


title: Adaptive Scheduling for Multi-Task Learning

abstract: To train neural machine translation models simultaneously on multiple tasks (languages), it is common to sample each task uniformly or in proportion to dataset sizes. As these methods offer little control over performance trade-offs, we explore different task scheduling approaches. We first consider existing non-adaptive techniques, then move on to adaptive schedules that over-sample tasks with poorer results compared to their respective baseline. As explicit schedules can be inefficient, especially if one task is highly over-sampled, we also consider implicit schedules, learning to scale learning rates or gradients of individual tasks instead. These techniques allow training multilingual models that perform better for low-resource language pairs (tasks with small amount of data), while minimizing negative effects on high-resource tasks.


title: A Random Network Model for the Analysis of Blockchain Designs with Communication Delay

abstract: This paper proposes a random network model for blockchains, a distributed hierarchical data structure of blocks that has found several applications in various industries. The model is parametric on two probability distribution functions governing block production and communication delay, which are key to capture the complexity of the mechanism used to synchronize the many distributed local copies of a blockchain. The proposed model is equipped with simulation algorithms for both bounded and unbounded number of distributed copies of the blockchain. They are used to study fast blockchain systems, i.e., blockchains in which the average time of block production can match the average time of message broadcasting used for blockchain synchronization. In particular, the model and the algorithms are useful to understand efficiency criteria associated with fast blockchains for identifying, e.g., when increasing the block production will have negative impact on the stability of the distributed data structure given the network's broadcast delay.


title: Coupling Rendering and Generative Adversarial Networks for Artificial SAS Image Generation

abstract: Acquisition of Synthetic Aperture Sonar (SAS) datasets is bottlenecked by the costly deployment of SAS imaging systems, and even when data acquisition is possible,the data is often skewed towards containing barren seafloor rather than objects of interest. We present a novel pipeline, called SAS GAN, which couples an optical renderer with a generative adversarial network (GAN) to synthesize realistic SAS images of targets on the seafloor. This coupling enables high levels of SAS image realism while enabling control over image geometry and parameters. We demonstrate qualitative results by presenting examples of images created with our pipeline. We also present quantitative results through the use of t-SNE and the Fr\'echet Inception Distance to argue that our generated SAS imagery potentially augments SAS datasets more effectively than an off-the-shelf GAN.


title: The Computational Complexity of Finding Temporal Paths under Waiting Time Constraints

abstract: Computing a (short) path between two vertices is one of the most fundamental primitives in graph algorithmics. In recent years, the study of paths in temporal graphs, that is, graphs where the vertex set is fixed but the edge set changes over time, gained more and more attention. A path is time-respecting, or temporal, if it uses edges with non-decreasing time stamps. We investigate a basic constraint for temporal paths, where the time spent at each vertex must not exceed a given duration $\Delta$, referred to as $\Delta$-restless temporal paths. This constraint arises naturally in the modeling of real-world processes like packet routing in communication networks and infection transmission routes of diseases where recovery confers lasting resistance. While finding temporal paths without waiting time restrictions is known to be doable in polynomial time, we show that the "restless variant" of this problem becomes computationally hard even in very restrictive settings. For example, it is W[1]-hard when parameterized by the distance to disjoint path of the underlying graph, which implies W[1]-hardness for many other parameters like feedback vertex number and pathwidth. A natural question is thus whether the problem becomes tractable in some natural settings. We explore several natural parameterizations, presenting FPT algorithms for three kinds of parameters: (1) output-related parameters (here, the maximum length of the path), (2) classical parameters applied to the underlying graph (e.g., feedback edge number), and (3) a new parameter called timed feedback vertex number, which captures finer-grained temporal features of the input temporal graph, and which may be of interest beyond this work.


title: Influence of nuclear spin polarization on spin echo signal of NV center qubit

abstract: We consider the spin echo dynamics of a nitrogen-vacancy center qubit based the $S\!= \! 1$ ground state spin manifold, caused by a dynamically polarized nuclear environment. We show that the echo signal acquires then a nontrivially time-dependent phase shift. This effect should be observable for polarization $\approx \! 0.5$ of nuclei within $\sim \! 1$ nm from the qubit, and for the NV center initialized in a superposition of $m\! = \! 0$ and either $m\! =\! 1$ or $m\! =\! -1$ states. This phase shift is much smaller when the NV center is prepared in a superposition of $m\! = \! 1$ and $m\! =\! -1$ states, i.e. when the qubit couples to the spin environment in a way analogous to that of spin-$1/2$. For nuclear environment devoid of spins strongly coupled to the qubit, the phase shift is well described within Gaussian approximation, which provides an explanation for the dependence of the shift magnitude on the choice of states on which the qubit is based, and makes it clear that its presence is related to the linear response of the environment perturbed by an evolving qubit. Consequently, its observation signifies the presence environment-mediated self-interaction of the qubit, and hence, it invalidates the notion that the nuclear environment acts as a source of external noise driving the qubit. We also show how a careful comparison of the echo signal from qubits based on $m\! = \! 0,1$ and $m\! =\! \pm 1$ manifolds, can distinguish between effectively Gaussian and non-Gaussian environment.


title: SuRF: a New Method for Sparse Variable Selection, with Application in Microbiome Data Analysis

abstract: In this paper, we present a new variable selection method for regression and classification purposes. Our method, called Subsampling Ranking Forward selection (SuRF), is based on LASSO penalised regression, subsampling and forward-selection methods. SuRF offers major advantages over existing variable selection methods in terms of both sparsity of selected models and model inference. We provide an R package that can implement our method for generalized linear models. We apply our method to classification problems from microbiome data, using a novel agglomeration approach to deal with the special tree-like correlation structure of the variables. Existing methods arbitrarily choose a taxonomic level a priori before performing the analysis, whereas by combining SuRF with these aggregated variables, we are able to identify the key biomarkers at the appropriate taxonomic level, as suggested by the data. We present simulations in multiple sparse settings to demonstrate that our approach performs better than several other popularly used existing approaches in recovering the true variables. We apply SuRF to two microbiome data sets: one about prediction of pouchitis and another for identifying samples from two healthy individuals. We find that SuRF can provide a better or comparable prediction with other methods while controlling the false positive rate of variable selection.


title: Sub-Kelvin cooling for two kilopixel bolometer arrays in the PIPER receiver

abstract: The Primordial Inflation Polarization Explorer (PIPER) is a balloon-borne telescope mission to search for inflationary gravitational waves from the early universe. PIPER employs two 32x40 arrays of superconducting transition-edge sensors, which operate at 100 mK. An open bucket dewar of liquid helium maintains the receiver and telescope optics at 1.7 K. We describe the thermal design of the receiver and sub-kelvin cooling with a continuous adiabatic demagnetization refrigerator (CADR). The CADR operates between 70-130 mK and provides ~10 uW cooling power at 100 mK, nearly five times the loading of the two detector assemblies. We describe electronics and software to robustly control the CADR, overall CADR performance in flight-like integrated receiver testing, and practical considerations for implementation in the balloon float environment.


title: MinneApple: A Benchmark Dataset for Apple Detection and Segmentation

abstract: In this work, we present a new dataset to advance the state-of-the-art in fruit detection, segmentation, and counting in orchard environments. While there has been significant recent interest in solving these problems, the lack of a unified dataset has made it difficult to compare results. We hope to enable direct comparisons by providing a large variety of high-resolution images acquired in orchards, together with human annotations of the fruit on trees. The fruits are labeled using polygonal masks for each object instance to aid in precise object detection, localization, and segmentation. Additionally, we provide data for patch-based counting of clustered fruits. Our dataset contains over 41, 000 annotated object instances in 1000 images. We present a detailed overview of the dataset together with baseline performance analysis for bounding box detection, segmentation, and fruit counting as well as representative results for yield estimation. We make this dataset publicly available and host a CodaLab challenge to encourage comparison of results on a common dataset. To download the data and learn more about MinneApple please see the project website: http://rsn.cs.umn.edu/index.php/MinneApple. Up to date information is available online.


title: Co-Attentive Cross-Modal Deep Learning for Medical Evidence Synthesis and Decision Making

abstract: Modern medicine requires generalised approaches to the synthesis and integration of multimodal data, often at different biological scales, that can be applied to a variety of evidence structures, such as complex disease analyses and epidemiological models. However, current methods are either slow and expensive, or ineffective due to the inability to model the complex relationships between data modes which differ in scale and format. We address these issues by proposing a cross-modal deep learning architecture and co-attention mechanism to accurately model the relationships between the different data modes, while further reducing patient diagnosis time. Differentiating Parkinson's Disease (PD) patients from healthy patients forms the basis of the evaluation. The model outperforms the previous state-of-the-art unimodal analysis by 2.35%, while also being 53% more parameter efficient than the industry standard cross-modal model. Furthermore, the evaluation of the attention coefficients allows for qualitative insights to be obtained. Through the coupling with bioinformatics, a novel link between the interferon-gamma-mediated pathway, DNA methylation and PD was identified. We believe that our approach is general and could optimise the process of medical evidence synthesis and decision making in an actionable way.


title: Electric field and current induced electroforming modes in NbOx

abstract: Electroforming is used to initiate the memristive response in metal/oxide/metal devices by creating a filamentary conduction path in the oxide film. Here we use a simple photoresist-based detection technique to map the spatial distribution of conductive filaments formed in Nb/NbOx/Pt devices, and correlate these with current-voltage characteristics and in-situ thermoreflectance measurements to identify distinct modes of electroforming in low and high conductivity NbOx films. In low conductivity films the filaments are randomly distributed within the oxide film, consistent with a field-induced weakest-link mechanism, while in high conductivity films they are concentrated in the center of the film. In the latter case the current-voltage characteristics and in-situ thermoreflectance imaging show that electroforming is associated with current bifurcation into regions of low and high current density. This is supported by finite element modelling of the current distribution and shown to be consistent with predictions of a simple core-shell model of the current distribution. These results clearly demonstrate two distinct modes of electroforming in the same materials system and show that the dominant mode depends on the conductivity of the film, with field-induced electroforming dominant in low conductivity films and current-bifurcation induced electroforming dominant in high conductivity films. Finally, we demonstrate S-type and snap-back negative differential resistance in the high conductivity films and explain this behavior in terms of two-zone model.


title: Local Decode and Update for Big Data Compression

abstract: This paper investigates data compression that simultaneously allows local decoding and local update. The main result is a universal compression scheme for memoryless sources with the following features. The rate can be made arbitrarily close to the entropy of the underlying source, contiguous fragments of the source can be recovered or updated by probing or modifying a number of codeword bits that is on average linear in the size of the fragment, and the overall encoding and decoding complexity is quasilinear in the blocklength of the source. In particular, the local decoding or update of a single message symbol can be performed by probing or modifying a constant number of codeword bits. This latter part improves over previous best known results for which local decodability or update efficiency grows logarithmically with blocklength.


title: Pump-triple sum-frequency-probe spectroscopy of transition metal dichalcogenides

abstract: Triple sum-frequency (TSF) spectroscopy measures multidimensional spectra by resonantly exciting multiple quantum coherences of vibrational and electronic states. In this work we demonstrate pump-TSF-probe spectroscopy in which a pump excites a sample and some time later three additional electric fields generate a probe field which is measured. We demonstrate pump-TSF-probe spectroscopy on polycrystalline, smooth, thin films and spiral nanostructures of both MoS2 and WS2. The pump-TSF-probe spectra are qualitatively similar to the more conventional transient-reflectance spectra. While transient-reflectance sensitivity suffers under low surface coverage, pump-TSF-probe sensitivity is independent of the sample coverage and nanostructure morphologies. Our results demonstrate that pump-TSF-probe is a valuable methodology for studying microscopic material systems.


title: Spatio-spectral networks for color-texture analysis

abstract: Texture is one of the most-studied visual attribute for image characterization since the 1960s. However, most hand-crafted descriptors are monochromatic, focusing on the gray scale images and discarding the color information. In this context, this work focus on a new method for color texture analysis considering all color channels in a more intrinsic approach. Our proposal consists of modeling color images as directed complex networks that we named Spatio-Spectral Network (SSN). Its topology includes within-channel edges that cover spatial patterns throughout individual image color channels, while between-channel edges tackle spectral properties of channel pairs in an opponent fashion. Image descriptors are obtained through a concise topological characterization of the modeled network in a multiscale approach with radially symmetric neighborhoods. Experiments with four datasets cover several aspects of color-texture analysis, and results demonstrate that SSN overcomes all the compared literature methods, including known deep convolutional networks, and also has the most stable performance between datasets, achieving $98.5(\pm1.1)$ of average accuracy against $97.1(\pm1.3)$ of MCND and $96.8(\pm3.2)$ of AlexNet. Additionally, an experiment verifies the performance of the methods under different color spaces, where results show that SSN also has higher performance and robustness.


title: Application of Genetic Algorithm for More Efficient Multi-Layer Thickness Optimization in Solar Cells

abstract: Thin-film solar cells are predominately designed similar to a stacked structure. Optimizing the layer thicknesses in this stack structure is crucial to extract the best efficiency of the solar cell. The commonplace method used in optimization simulations, such as for optimizing the optical spacer layers' thicknesses, is the parameter sweep. Our simulation study shows that the implementation of a meta-heuristic method like the genetic algorithm results in a significantly faster and accurate search method when compared to the brute-force parameter sweep method in both single and multi-layer optimization. While other sweep methods can also outperform the brute-force method, they do not consistently exhibit $100\%$ accuracy in the optimized results like our genetic algorithm. We have used a well-studied P3HT-based structure to test our algorithm. Our best-case scenario was observed to use $60.84\%$ fewer simulations than the brute-force method.


title: A look at the way we look at Complex Networks Robustness and Resilience

abstract: This article offers a brief overview of the current research topics concerning strategies to mitigate the adverse effects of perturbations in complex networks. It addresses the issue of an unclear use of Robustness and Resilience terminology and proposes a common interdisciplinary framework for comparing strategies across different fields. It concludes with a high-level discussion of existing challenges and suggestions for future research.


title: Effect of the electron-lattice coupling on the charge and magnetic order in rare-earth nickelates

abstract: We investigate the impact of electron-lattice coupling on the stability of various magnetic orders in rare-earth nickelates. We use the Hartree-Fock approximation, at zero temperature, to study an effective, two-band model with correlations characterized by a Hubbard $U$ and a Hund's $J$. This is coupled to breathing-mode distortions of the octahedral oxygen cages, described semi-classically, with a Holstein term. We analyze the effect of the various parameters on the resulting phase diagram, in particular on the charge disproportionation and on the magnetic order. We confirm that the coupling to the lattice cooperates with Hund's coupling and thus encourages charge disproportionation. We also find that it favors the fully disproportionated, 4-site periodic magnetic order of type $\Uparrow 0 \Downarrow 0$. Other convergent magnetic phases, such as the collinear $\uparrow\uparrow\downarrow\downarrow$ and non-collinear $\uparrow\rightarrow\downarrow\leftarrow$ states, do not couple to the lattice because of their lack of charge disproportionation. Novel phases, e.g. with charge disproportionation but no magnetic order, are also found to be stabilized in specific conditions.


title: Magnetospheric Accretion as a Source of H$\alpha$ Emission from Proto-planets around PDS 70

abstract: Advances in high-resolution imaging have revealed H$\alpha$ emission separated from the host star. It is generally believed that the emission is associated with forming planets in protoplanetary disks. However, the nature of this emission is still not fully understood. Here we report a modeling effort of H$\alpha$ emission from the planets around the young star PDS 70. Using standard magnetospheric accretion models previously applied to accreting young stars, we find that the observed line fluxes can be reproduced using a range of parameters relevant to PDS 70b and c, with the mean mass accretion rate of log(${\rm \dot{M}}$) = $-8.0\pm0.6$ M$_{\rm Jup}$ yr$^{-1}$ and $-8.1\pm0.6$ M$_{\rm Jup}$ yr$^{-1}$ for PDS 70b and PDS 70c, respectively. Our results suggest that H$\alpha$ emission from young planets can originate in the magnetospheric accretion of mass from the circumplanetary disk. We find that empirical relationships between mass accretion rate and H$\alpha$ line properties frequently used in T Tauri stars are not applicable in the planetary mass regime. In particular, the correlations between line flux and mass accretion rate underpredict the accretion rate by about an order of magnitude, and the width at the 10% height of the line is insensitive to the accretion rate at ${\rm \dot{M}}$ $< 10^{-8}$ M$_{\rm Jup}$ yr$^{-1}$.


title: Distributed Focus and Digital Zoom

abstract: We explore integrated microcamera focus systems for array cameras. We propose a new model for system camera integration relying on fast action focus mechanisms with >10mm aperture. Rather than reducing resolution or expanding aperture size, such systems can be used in arrays to enable digital zoom. We show that a common mechanism supports camera modules with focal lengths ranging from 25 to 60 mm. Designs for each focal length include a fixed objective lens group and an adjustable back focus group. Increasing the focal power of the front focal group enables the travel range of available microcamera modules to accommodate long focal length systems. We present design examples both discrete and multiscale array camera systems.


title: Direct transformation from Cartesian into geodetic coordinates on a triaxial ellipsoid

abstract: This paper presents two new direct symbolic-numerical algorithms for the transformation of Cartesian coordinates into geodetic coordinates considering the general case of a triaxial reference ellipsoid. The problem in both algorithms is reduced to finding a real positive root of a sixth degree polynomial. The first approach consists of algebraic manipulations of the equations describing the geometry of the problem and the second one uses Gr\"obner bases. In order to perform numerical tests and accurately compare efficiency and reliability, our algorithms together with the iterative methods presented by M. Ligas (2012) and J. Feltens (2009) have been implemented in C++. The numerical tests have been accomplished by considering 10 celestial bodies, referenced in the available literature. The obtained results clearly show that our algorithms improve the aforementioned iterative methods, in terms of both efficiency and accuracy.


title: SUSY S-matrix Bootstrap and Friends

abstract: We consider the 2D S-matrix bootstrap in the presence of supersymmetry, $\mathbb{Z}_2$ and $\mathbb{Z}_4$ symmetry. At the boundary of the allowed S-matrix space we encounter well known integrable models such as the supersymmetric sine-Gordon and restricted sine-Gordon models, novel elliptic deformations thereof, as well as a two parameter family of $\mathbb{Z}_4$ elliptic S-matrices previously proposed by Zamolodchikov. We highlight an intricate web of relations between these various S-matrices.


title: Collective Multi-mode Vacuum Rabi Splitting

abstract: We report the experimental observation of collective multi-mode vacuum Rabi splitting in free space. In contrast to optical cavities, the atoms couple to a continuum of modes, and the optical thickness of the cloud provides a measure of this coupling. The splitting, also referred as normal mode splitting, is monitored through the Rabi oscillations in the scattered intensity, and the results are fully explained by a linear-dispersion theory.


title: A data-driven method for quantifying the impact of a genetic circuit on its host

abstract: Genetic circuits are designed to implement certain logic in living cells, keeping burden on the host cell minimal. However, manipulating the genome often will have a significant impact for various reasons (usage of the cell machinery to express new genes, toxicity of genes, interactions with native genes, etc.). In this work we utilize Koopman operator theory to construct data-driven models of transcriptomic-level dynamics from noisy and temporally sparse RNAseq measurements. We show how Koopman models can be used to quantify impact on genetic circuits. We consider an experimental example, using high-throughput RNAseq measurements collected from wild-type E. coli, single gate components transformed in E. coli, and a NAND circuit composed from individual gates in E. coli, to explore how Koopman subspace functions encode increasing circuit interference on E. coli chassis dynamics. The algorithm provides a novel method for quantifying the impact of synthetic biological circuits on host-chassis dynamics.


title: Scalable approximation of Green's function for estimation of anharmonic energy corrections

abstract: A method based on separated integration to estimate anharmonic corrections to energy and vibration of molecules in a second-order diagrammatic vibrational many-body Green's function formalism has already been presented. A severe bottleneck in extending this approach to bigger molecules is that the storage of the Green's function scales exponentially with the number of atoms in the molecule. In this article, we present a method that overcomes this limitation by approximating the Green's function in the Hierarchical Tucker tensor format. We illustrate that the storage cost is linear in dimension and hence one can obtain accurate representations of the Green's function for a molecule of any size. Application of this method to estimate the second-order correction to energy of molecules illustrates the advantage of this approach.


title: Linear Size Planar Manhattan Network for Convex Point Sets

abstract: Let $G = (V, E)$ be an edge-weighted geometric graph such that every edge is horizontal or vertical. The weight of an edge $uv \in E$ is its length. Let $ W_G (u,v)$ denote the length of a shortest path between a pair of vertices $u$ and $v$ in $G$. The graph $G$ is said to be a Manhattan network for a given point set $ P $ in the plane if $P \subseteq V$ and $\forall p,q \in P$, $ W_G (p,q)=|pq|_1$. In addition to $ P$, graph $G$ may also include a set $T$ of Steiner points in its vertex set $V$. In the Manhattan network problem, the objective is to construct a Manhattan network of small size for a set of $ n $ points. This problem was first considered by Gudmundsson et al.\cite{gudmundsson2007small}. They give a construction of a Manhattan network of size $\Theta(n \log n)$ for general point set in the plane. We say a Manhattan network is planar if it can be embedded in the plane without any edge crossings. In this paper, we construct a linear size planar Manhattan network for convex point set in linear time using $\mathcal{ O}(n)$ Steiner points. We also show that, even for convex point set, the construction in Gudmundsson et al. \cite{gudmundsson2007small} needs $\Omega (n \log n)$ Steiner points and the network may not be planar.


title: Diamond Mirrors for High-Power Lasers

abstract: High-power lasers have numerous scientific and industrial applications. Some key areas include laser cutting and welding in manufacturing, directed energy in fusion reactors or defense applications, laser surgery in medicine, and advanced photolithography in the semiconductor industry. These applications require optical components, in particular mirrors, that withstand high optical powers for directing light from the laser to the target. Ordinarily, mirrors are comprised of multilayer coatings of different refractive index and thickness. At high powers, imperfections in these layers lead to absorption of light, resulting in thermal stress and permanent damage to the mirror. Here we design, simulate, fabricate, and demonstrate monolithic and highly reflective dielectric mirrors which operate under high laser powers without damage. The mirrors are realized by etching nanostructures into the surface of single-crystal diamond, a material with exceptional optical and thermal properties. We measure reflectivities of greater than 98% and demonstrate damage-free operation using 10 kW of continuous-wave laser light at 1070 nm, with intensities up to 4.6 MW/cm2. In contrast, at these laser powers, we observe damage to a standard dielectric mirror based on optical coatings. Our results initiate a new category of broadband optics that operate in extreme conditions.


title: F-Cooper: Feature based Cooperative Perception for Autonomous Vehicle Edge Computing System Using 3D Point Clouds

abstract: Autonomous vehicles are heavily reliant upon their sensors to perfect the perception of surrounding environments, however, with the current state of technology, the data which a vehicle uses is confined to that from its own sensors. Data sharing between vehicles and/or edge servers is limited by the available network bandwidth and the stringent real-time constraints of autonomous driving applications. To address these issues, we propose a point cloud feature based cooperative perception framework (F-Cooper) for connected autonomous vehicles to achieve a better object detection precision. Not only will feature based data be sufficient for the training process, we also use the features' intrinsically small size to achieve real-time edge computing, without running the risk of congesting the network. Our experiment results show that by fusing features, we are able to achieve a better object detection result, around 10% improvement for detection within 20 meters and 30% for further distances, as well as achieve faster edge computing with a low communication delay, requiring 71 milliseconds in certain feature selections. To the best of our knowledge, we are the first to introduce feature-level data fusion to connected autonomous vehicles for the purpose of enhancing object detection and making real-time edge computing on inter-vehicle data feasible for autonomous vehicles.


title: Reduced order models for spectral domain inversion: Embedding into the continuous problem and generation of internal data

abstract: We generate data-driven reduced order models (ROMs) for inversion of the one and two dimensional Schr\"odinger equation in the spectral domain given boundary data at a few frequencies. The ROM is the Galerkin projection of the Schr\"odinger operator onto the space spanned by solutions at these sample frequencies. The ROM matrix is in general full, and not good for extracting the potential. However, using an orthogonal change of basis via Lanczos iteration, we can transform the ROM to a block triadiagonal form from which it is easier to extract $q$. In one dimension, the tridiagonal matrix corresponds to a three-point staggered finite-difference system for the Schr\"odinger operator discretized on a so-called spectrally matched grid which is almost independent of the medium. In higher dimensions, the orthogonalized basis functions play the role of the grid steps. The orthogonalized basis functions are localized and also depend only very weakly on the medium, and thus by embedding into the continuous problem, the reduced order model yields highly accurate internal solutions. That is to say, we can obtain, just from boundary data, very good approximations of the solution of the Schr\"odinger equation in the whole domain for a spectral interval that includes the sample frequencies. We present inversion experiments based on the internal solutions in one and two dimensions.


title: Spontaneous time-reversal symmetry breaking without magnetism in a $S=1$ chain

abstract: States of matter that break time-reversal symmetry are invariably associated with magnetism or circulating currents. Recently, one of us proposed a phase, the directional scalar spin chiral order (DSSCO), as an exception: it breaks time-reversal symmetry via chiral ordering of spins along a particular direction, but is spin-rotation symmetric. In this work, we prove the existence of this state via state-of-the-art density matrix renormalization group (DMRG) analysis on a spin-1 chain with nearest-neighbor bilinear-biquadratic interactions and additional third-neighbor ferromagnetic Heisenberg exchange. Despite the large entanglement introduced by the third-neighbor coupling, we are able to access system sizes up to $L=918$ sites. We find first order phase transitions from the DSSCO into the famous Haldane phase as well as a spin-quadrupolar phase where spin nematic correlations dominate. In the Haldane phase, we propose and demonstrate a method for detecting the topological edge states using DMRG that could be useful for other topological phases too.


title: A transparent referendum protocol with immutable proceedings and verifiable outcome for trustless networks

abstract: High voter turnout in elections and referendums is very desirable in order to ensure a robust democracy. Secure electronic voting is a vision for the future of elections and referendums. Such a system can counteract factors that hinder strong voter turnout such as the requirement of physical presence during limited hours at polling stations. However, this vision brings transparency and confidentiality requirements that render the design of such solutions challenging. Specifically, the counting must be implemented in a reproducible way and the ballots of individual voters must remain concealed. In this paper, we propose and evaluate a referendum protocol that ensures transparency, confidentiality, and integrity, in trustless networks. The protocol is built by combining Secure Multi-Party Computation (SMPC) and Distributed Ledger or Blockchain technology. The persistence and immutability of the protocol communication allows verifiability of the referendum outcome on the client side. Voters therefore do not need to trust in third parties. We provide a formal description and conduct a thorough security evaluation of our proposal.


title: Optimization on the Surface of the (Hyper)-Sphere

abstract: Thomson problem is a classical problem in physics to study how $n$ number of charged particles distribute themselves on the surface of a sphere of $k$ dimensions. When $k=2$, i.e. a 2-sphere (a circle), the particles appear at equally spaced points. Such a configuration can be computed analytically. However, for higher dimensions such as $k \ge 3$, i.e. the case of 3-sphere (standard sphere), there is not much that is understood analytically. Finding global minimum of the problem under these settings is particularly tough since the optimization problem becomes increasingly computationally intensive with larger values of $k$ and $n$. In this work, we explore a wide variety of numerical optimization methods to solve the Thomson problem. In our empirical study, we find stochastic gradient based methods (SGD) to be a compelling choice for this problem as it scales well with the number of points.


title: Finite element modeling of micropolar-based phononic crystals

abstract: The performance of a Cosserat/micropolar solid as a numerical vehicle to represent dispersive media is explored. The study is conducted using the finite element method with emphasis on Hermiticity, positive definiteness, principle of virtual work and Bloch-Floquet boundary conditions. The periodic boundary conditions are given for both translational and rotational degrees of freedom and for the associated force- and couple-traction vectors. Results in terms of band structures for different material cells and mechanical parameters are provided.


title: Nonlocal interferences induced by the phase of the wavefunction for a particle in a cavity with moving boundaries

abstract: We investigate the dynamics of a particle in a confined periodic system---a time-dependent oscillator confined by infinitely high and moving walls---and focus on the evolution of the phase of the wavefunction. It is shown that, for some specific initial states in this potential, the phase of the wavefunction throughout the cavity depends on the walls motion. We further elaborate a thought experiment based on interferences devised to detect this form of single-particle nonlocality from a relative phase. We point out that, within the non-relativistic formalism based on the Schr\"odinger equation (SE), detecting this form of nonlocality can give rise to signaling. We believe this effect is an artifact, but the standard relativistic corrections to the SE do not appear to fix it. Specific illustrations are given, with analytical results in the adiabatic approximation, and numerical computations to show that contributions from high-energy states (corresponding to superluminal velocities) are negligible.


title: Thermal effects of Very Special Relativity Quantum Electrodynamics

abstract: In this paper we compute the high-temperature effective Lagrangian for the quantum electrodynamics defined in the framework of Very Special Relativity (VSR). The main aspect of the VSR setting is that it modifies the gauge invariance, admitting now different types of interactions appearing in a nonlocal form. In order to explore the richness of these new couplings, we employ the usual Matsubara imaginary-time formalism to evaluate the effective Lagrangian at one and two-loop order. We illustrate the leading VSR Lorentz violation modifications by computing high-temperature internal energy density and establishing a comparison with the expected contribution.


title: Cleaning our own Dust: Simulating and Separating Galactic Dust Foregrounds with Neural Networks

abstract: Separating galactic foreground emission from maps of the cosmic microwave background (CMB), and quantifying the uncertainty in the CMB maps due to errors in foreground separation are important for avoiding biases in scientific conclusions. Our ability to quantify such uncertainty is limited by our lack of a model for the statistical distribution of the foreground emission. Here we use a Deep Convolutional Generative Adversarial Network (DCGAN) to create an effective non-Gaussian statistical model for intensity of emission by interstellar dust. For training data we use a set of dust maps inferred from observations by the Planck satellite. A DCGAN is uniquely suited for such unsupervised learning tasks as it can learn to model a complex non-Gaussian distribution directly from examples. We then use these simulations to train a second neural network to estimate the underlying CMB signal from dust-contaminated maps. We discuss other potential uses for the trained DCGAN, and the generalization to polarized emission from both dust and synchrotron.


title: Continuous Abstraction of Nonlinear Systems using Sum-of-Squares Programming

abstract: We present a control design procedure for nonlinear control systems in which we represent a potentially high dimensional system with a low dimensional continuous-state abstraction. The abstraction generates a reference which the original system follows with a low level controller. We propose sum-of-squares programming as a tool to design this controller and to provide an upper bound on the relative error between the system and its abstraction. We compute the low level controller simultaneously with a simulation function that gives the boundedness guarantee for the relative error.


title: Distribution of laser shot noise energy delivered to a levitated nanoparticle

abstract: This paper quantifies the rate at which laser shot noise energy is delivered to a nanoparticle for the various scenarios commonly encountered in levitated optomechanics. While previous articles have the same form and dependencies, the proportionality constants often differ in the literature. This paper resolves these discrepancies. The rate at which energy is delivered to an optically trapped particle's respective degrees of freedom depends on the radiation pattern of scattered light as well as the direction of laser propagation. For a traveling plane wave with linearly polarized light, in the Rayleigh regime this leads the translational shot noise heating rate to be proportional to $1/10$ of the total rate in the laser polarization direction, $7/10$ in the laser propagation direction, and $2/10$ in the direction perpendicular to both. Analytical expressions for the shot noise heating rate are provided in the Rayleigh limit as well as numerical calculations for particles in the Mie regime for silica and diamond. For completeness, numerical calculations of the shot noise heating for silica Mie particles at the focus of a strongly focused laser beam are calculated for varying numerical aperture and common laser wavelengths. Both numerical calculations show that the Rayleigh expression generally gives an overestimate of the shot noise heating especially for larger radii, but is still a good approximation even for incident focal fields. The exception to the relative decrease is when a Mie resonance is reached which was found for diamond. Lastly, Rayleigh expressions for the rotational shot noise heating for a symmetric top-like particle for linear, elliptically, and unpolarized light are also provided.


title: Fuzzy Logic Control for Mixed conventional/braking Actuation Mobile Robots

abstract: The use of conventional actuators in robotic systems (electric motors in particular), while often offering advantages in terms of flexibility and controllability, suffer from primary actuator failure,due to unexpected complexities in their environment, which can lead to loss of controllability. Conventional actuators can impose disadvantages on mechanical complexity, weight, and cost. Here,the Mixed conventional/braking Actuation Mobile Robot (MAMR), a new mobile robot platform,is proposed to tackle such drawbacks in actuation and explore the use and control of braking actuation. This platform substitutes the drive motors used in Ackermann steering with brakes that have only two states, ON and OFF. Additionally, the conventional drive wheels are replaced by a single,omni-directional wheel that only supports a driving force in the robots longitudinal direction. The ability of braking actuators in providing controllability under actuator failure is one of the primary motivations of this work. To validate the reliability and accuracy of MAMR approach, this paper studies the design of such robotic systems, the design and synthesis of fuzzy logic controllers along with the experimental assessments of these controllers in real-time. The experimental tests point out the controller performance enhancement using fuzzy logic controllers and MAMR.


title: Theory of anisotropic elastoresistivity of two-dimensional extremely strongly correlated metals

abstract: There is considerable recent interest in the phenomenon of anisotropic electroresistivity of correlated metals. While some interesting work has been done on the iron-based superconducting systems, not much is known for the cuprate materials. Here we study the anisotropy of elastoresistivity for cuprates in the normal state. We present theoretical results for the effect of strain on resistivity, and additionally on the optical weight and local density of states. We use the recently developed extremely strongly correlated Fermi liquid theory in two dimensions, which accounts quantitatively for the unstrained resistivities for three families of single-layer cuprates. The strained hoppings of a tight-binding model are roughly modeled analogously to strained transition metals. The strained resistivity for a two-dimensional $t$-$t'$-$J$ model are then obtained, using the equations developed in recent work. Our quantitative predictions for these quantities have the prospect of experimental tests in the near future, for strongly correlated materials such as the hole-doped and electron-doped high-$T_c$ materials.


title: P$^2$IM: Scalable and Hardware-independent Firmware Testing via Automatic Peripheral Interface Modeling (extended version)

abstract: Dynamic testing or fuzzing of embedded firmware is severely limited by hardware-dependence and poor scalability, partly contributing to the widespread vulnerable IoT devices. We propose a software framework that continuously executes a given firmware binary while channeling inputs from an off-the-shelf fuzzer, enabling hardware-independent and scalable firmware testing. Our framework, using a novel technique called P$^2$IM, abstracts diverse peripherals and handles firmware I/O on the fly based on automatically generated models. P$^2$IM is oblivious to peripheral designs and generic to firmware implementations, and therefore, applicable to a wide range of embedded devices. We evaluated our framework using 70 sample firmware and 10 firmware from real devices, including a drone, a robot, and a PLC. It successfully executed 79% of the sample firmware without any manual assistance. We also performed a limited fuzzing test on the real firmware, which unveiled 7 unique unknown bugs.


title: Learned imaging with constraints and uncertainty quantification

abstract: We outline new approaches to incorporate ideas from deep learning into wave-based least-squares imaging. The aim, and main contribution of this work, is the combination of handcrafted constraints with deep convolutional neural networks, as a way to harness their remarkable ease of generating natural images. The mathematical basis underlying our method is the expectation-maximization framework, where data are divided in batches and coupled to additional "latent" unknowns. These unknowns are pairs of elements from the original unknown space (but now coupled to a specific data batch) and network inputs. In this setting, the neural network controls the similarity between these additional parameters, acting as a "center" variable. The resulting problem amounts to a maximum-likelihood estimation of the network parameters when the augmented data model is marginalized over the latent variables.


title: Rethinking the Micro-Foundation of Opinion Dynamics: Rich Consequences of the Weighted-Median Mechanism

abstract: To identify the main mechanisms underlying complex opinion formation processes in social systems, researchers have long been exploring simple mechanistic mathematical models. Most existing opinion dynamics models are built on a common micro-foundation, i.e., the weighted-averaging opinion update. However, we argue that this universally-adopted mechanism features a non-negligible unrealistic feature, which brings unnecessary difficulties in seeking a proper balance between model complexity and predictive power. In this paper, we propose the weighted-median mechanism as a new micro-foundation of opinion dynamics, which, with minimal assumptions, fundamentally resolves the inherent unrealistic feature of the weighted-averaging mechanism. Derived from the cognitive dissonance theory in psychology, the weighted-median mechanism is supported by online experiment data and broadens the applicability of opinion dynamics models to multiple-choice issues with ordered discrete options. Moreover, the weighted-median mechanism, despite being the simplest in form, captures various non-trivial real-world features of opinion evolution, while some widely-studied averaging-based models fail to.


title: Good sequencings for small Mendelsohn triple systems

abstract: A Mendelsohn triple system of order $v$ (or MTS$(v)$) is a decomposition of the complete graph into directed 3-cyles. We denote the directed 3-cycle with edges $(x,y)$, $(y,z)$ and $(z,x)$ by $(x,y,z)$, $(y,z,x)$ or $(z,x,y)$. An $\ell$-good sequencing of a MTS$(v)$ is a permutation of the points of the design, say $[x_1 \; \cdots \; x_v]$, such that, for every triple $(x,y,z)$ in the design, it is not the case that $x = x_i$, $y = x_j$ and $z = x_k$ with $i < j < k$ and $k-i+1 \leq \ell$; or with $j < k < i$ and $i-j+1 \leq \ell$; or with $k < i < j$ and $j-k+1 \leq \ell$.


title: Some improvement on non-parametric estimation of income distribution and poverty index

abstract: In this paper, we propose an estimator of Foster, Greer and Thorbecke class of measures $\displaystyle P(z,\alpha) = \int_0^{z}\Big(\frac{z-x}{z}\Big)^{\alpha}f(x)\, dx$, where $z>0$ is the poverty line, $f$ is the probabily density function of the income distribution and $\alpha$ is the so-called poverty aversion. The estimator is constructed with a bias reduced kernel estimator. Uniform almost sure consistency and uniform mean square consistenty are established. A simulation study indicates that our new estimator performs well.


title: Combating Conservativeness in Data-Driven Optimization under Uncertainty: A Solution Path Approach

abstract: In data-driven optimization, solution feasibility is often ensured through a "safe" reformulation of the uncertain constraints, such that an obtained data-driven solution is guaranteed to be feasible for the oracle formulation with high statistical confidence. Such approaches generally involve an implicit estimation of the whole feasible set that can scale rapidly with the problem dimension, in turn leading to over-conservative solutions. In this paper, we investigate a validation-based strategy to avoid set estimation by exploiting the intrinsic low dimensionality among all possible solutions output from a given reformulation. We demonstrate how our obtained solutions satisfy statistical feasibility guarantees with light dimension dependence, and how they are asymptotically optimal and thus regarded as the least conservative with respect to the considered reformulation classes. We apply this strategy to several data-driven optimization paradigms including (distributionally) robust optimization, sample average approximation and scenario optimization. Numerical experiments show encouraging performances of our strategy compared to established benchmarks.


title: Sliding Mode Control for Mixed Conventional/Braking Actuation Mobile Robots

abstract: The Mixed convention/braking Actuation Mobile Robot (MAMR) was designed to tackle some of the drawbacks of conventional mobile robots such as losing controllability due to primary actuator failures, mechanical complexity, weight, and cost. It replaces conventional steering wheels with braking actuators and conventional drive wheels with a single omni-directional wheel. This makes it fall under the category of under-actuated mobile robots. The brakes have only two states, ON and OFF, resulting in discontinuous dynamics. This inspires the use of a discontinuous control law to control the system. This work presents a Sliding Mode Controller (SMC) design to park the MAMR system from a given initial configuration to a desired final configuration. Experimental results are presented to validate the parking control of the MAMR.


title: Decentralized Proximal Gradient Algorithms with Linear Convergence Rates

abstract: This work studies a class of non-smooth decentralized multi-agent optimization problems where the agents aim at minimizing a sum of local strongly-convex smooth components plus a common non-smooth term. We propose a general primal-dual algorithmic framework that unifies many existing state-of-the-art algorithms. We establish linear convergence of the proposed method to the exact solution in the presence of the non-smooth term. Moreover, for the more general class of problems with agent specific non-smooth terms, we show that linear convergence cannot be achieved (in the worst case) for the class of algorithms that uses the gradients and the proximal mappings of the smooth and non-smooth parts, respectively. We further provide a numerical counterexample that shows how some state-of-the-art algorithms fail to converge linearly for strongly-convex objectives and different local non-smooth terms.


title: An Alert-Generation Framework for Improving Resiliency in Human-Supervised, Multi-Agent Teams

abstract: Human-supervision in multi-agent teams is a critical requirement to ensure that the decision-maker's risk preferences are utilized to assign tasks to robots. In stressful complex missions that pose risk to human health and life, such as humanitarian-assistance and disaster-relief missions, human mistakes or delays in tasking robots can adversely affect the mission. To assist human decision making in such missions, we present an alert-generation framework capable of detecting various modes of potential failure or performance degradation. We demonstrate that our framework, based on state machine simulation and formal methods, offers probabilistic modeling to estimate the likelihood of unfavorable events. We introduce smart simulation that offers a computationally-efficient way of detecting low-probability situations compared to standard Monte-Carlo simulations. Moreover, for certain class of problems, our inference-based method can provide guarantees on correctly detecting task failures.


title: Modeling the Saturation of the Bell Instability Using Hybrid Simulations

abstract: The nonresonant cosmic ray instability, predicted by Bell (2004), is thought to play an important role in the acceleration and confinement of cosmic rays (CRs) close to supernova remnants. Despite its importance, the exact mechanism responsible for the saturation of the instability has not been determined, and there is no first-principle prediction for the amplitude of the saturated magnetic field. Using a survey of self-consistent kinetic hybrid simulations (with kinetic ions and fluid electrons), we study the saturation of the non-resonant streaming instability as a function of the parameters of both the thermal background plasma and the CR population. The strength of the saturated magnetic field has important implications for both CR acceleration in supernova remnants and CR diffusion in the Galaxy.


title: Ptychographic phase-retrieval by proximal algorithms

abstract: We derive a set of ptychography phase-retrieval iterative engines based on proximal algorithms originally developed in convex optimization theory, and discuss their connections with existing ones. The use of proximal operator creates a simple frame work that allows us to incorporate the effect of noise from a maximum-likelihood principle. We focus on three particular algorithms, namely proximal minimization, alternating direction method of multiplier and accelerated proximal gradient, and benckmark their performance with numerical simulations and experimental x-ray data. Among them, accelerated proximal gradient shows superior performance in terms of both accuracy and convergence rate for a noisy dataset.


title: Repulsion in controversial debate drives public opinion into fifty-fifty stalemate

abstract: Opinion formation is a process with strong implications for public policy. In controversial debates with large consequences, the public opinion is often trapped in a fifty-fifty stalemate, jeopardizing broadly accepted political decisions. Emergent effects from millions of private discussions make it hard to understand or influence this kind of opinion dynamics. Here we demonstrate that repulsion from opinions favors fifty-fifty stalemates. We study a voter model where agents can have two opinions or an undecided state in-between. In pairwise discussions, undecided agents can be convinced or repelled from the opinion expressed by another agent. If repulsion happens in at least one of four cases, as in controversial debates, the frequencies of both opinions equalize. Further we include transitions of decided agents to the undecided state. If that happens often, the share of undecided agents becomes large, as can be measured with the share of undecided answers in polls.


title: The scattering matrix for 0th order pseudodifferential operators

abstract: We use microlocal radial estimates to prove the full limiting absorption principle for $P$, a self-adjoint 0th order pseudodifferential operator satisfying hyperbolic dynamical assumptions as of Colin de Verdi\`ere and Saint-Raymond. We define the scattering matrix for $P-\omega$ with generic $\omega \in \mathbb R$ and show that the scattering matrix extends to a unitary operator on appropriate $L^2$ spaces. After conjugation with natural reference operators, the scattering matrix becomes a $0$th order Fourier integral operator with a canonical relation associated to the bicharacteristics of $P-\omega$. The operator $P$ gives a microlocal model of internal waves in stratified fluids as illustrated in the paper of Colin de Verdi\`ere and Saint-Raymond.


title: Multi-Perspective, Simultaneous Embedding

abstract: We describe MPSE: a Multi-Perspective Simultaneous Embedding method for visualizing high-dimensional data, based on multiple pairwise distances between the data points. Specifically, MPSE computes positions for the points in 3D and provides different views into the data by means of 2D projections (planes) that preserve each of the given distance matrices. We consider two versions of the problem: fixed projections and variable projections. MPSE with fixed projections takes as input a set of pairwise distance matrices defined on the data points, along with the same number of projections and embeds the points in 3D so that the pairwise distances are preserved in the given projections. MPSE with variable projections takes as input a set of pairwise distance matrices and embeds the points in 3D while also computing the appropriate projections that preserve the pairwise distances. The proposed approach can be useful in multiple scenarios: from creating simultaneous embedding of multiple graphs on the same set of vertices, to reconstructing a 3D object from multiple 2D snapshots, to analyzing data from multiple points of view. We provide a functional prototype of MPSE that is based on an adaptive and stochastic generalization of multi-dimensional scaling to multiple distances and multiple variable projections. We provide an extensive quantitative evaluation with datasets of different sizes and using different number of projections, as well as several examples that illustrate the quality of the resulting solutions.


title: Deep Learning Realm for Geophysics: Seismic Acquisition, Processing, Interpretation, and Inversion

abstract: Applying deep-learning models to geophysical applications has attracted special attentions during the past a couple of years. There are several papers published in this domain involving with different topics primarily focusing on synthetic data. Based on the presented results, reaching the satisfaction from geophysicists is largely far away the performance of the conventional algorithms. Moreover, it is extremely hard to find documents with detailed illustrations about what are the best practices that researchers should follow regarding how to design appropriate deep-learning models to effectively and precisely tackle problems relevant to the field data. This paper serves as a summarization to demonstrate successful stories and share with extensive experiences we have gained during past several years in the process of designing and deploying of deep-learning models to the geophysical projects in a large scale. Four different disciplines are discussed individually with seismic acquisition, processing, interpretation, and inversion. Finally, special attentions about designing an effective deep-neural-networks especially for geophysics are discussed.


title: Status of the Keck Planet Imager and Characterizer Phase II Development

abstract: The Keck Planet Imager and Characterizer comprises of a series of upgrades to the Keck II adaptive optics system and instrument suite to improve the direct imaging and high resolution spectroscopy capabilities of the facility instruments NIRC2 and NIRSPEC, respectively. Phase I of KPIC includes a NIR pyramid wavefront sensor and a Fiber Injection Unit (FIU) to feed NIRSPEC with a single mode fiber, which have already been installed and are currently undergoing commissioning. KPIC will enable High Dispersion Coronagraphy (HDC) of directly imaged exoplanets for the first time, providing potentially improved detection significance and spectral characterization capabilities compared to direct imaging. In favorable cases, Doppler imaging, spin measurements, and molecule mapping are also possible. This science goal drives the development of phase II of KPIC, which is scheduled to be deployed in early 2020. Phase II optimizes the system throughput and contrast using a variety of additional submodules, including a 952 element deformable mirror, phase induced amplitude apodization lenses, an atmospheric dispersion compensator, multiple coronagraphs, a Zernike wavefront sensor, and multiple science ports. A testbed is being built in the Exoplanet Technology Lab at Caltech to characterize and test the design of each of these submodules before KPIC phase II is deployed to Keck. This paper presents an overview of the design of phase II and report on results from laboratory testing.


title: Topological charge pumping in spin-dependent superlattices with glide symmetry

abstract: Topological charge pumping represents an important quantum phenomenon that shows the fundamental connection to the topological properties of dynamical systems. Here, we introduce a pumping process in a spin-dependent double-well optical lattice with glide symmetry. In the dynamic process, the glide symmetry protects the band touching points and topological properties of the system are characterised by the non-Abelian Berry curvature. By engineering suitable form of interaction between different spin components, the model not only demonstrates topological phase transition, but also shows hybridisation between the spatial and temporal domain with novel topological features captured by the Wilson line along the synthetic directions. Our work provides a new model based on ultracold atoms towards the implementation of versatile topological matters and topological phenomena.


title: Aero-Resonant Migration

abstract: The process of planet conglomeration, which primarily unfolds in a geometrically thin disk of gas and dust, is often accompanied by dynamical excitation of the forming planets and planetesimals. The ensuing orbital crossing can lead to large-scale collisional fragmentation, populating the system with icy and rocky debris. In a gaseous nebula, such leftover solid matter tends to spiral down towards the host star due to aerodynamic drag. Along the way, the inward drifting debris can encounter planets and gravitationally couple to them via mean-motion resonances, sapping them of their orbital energy and causing them to migrate. Here, we develop a simple theory for this migration mechanism, which we call "Aero-Resonant Migration" (ARM), in which small planetesimals ($10\,$m $\lesssim s \lesssim 10\,$km) undergo orbital decay due to aerodynamic drag and resonantly shepherd planets ahead of them. Using a combination of analytical calculations and numerical experiments, we show that ARM is a robust migration mechanism, able to significantly transport planets on timescales $\lesssim1$ Myr, and present simple formulae for the ARM rate.


title: Relativistic, Viscous, Radiation Hydrodynamic Simulations of Geometrically Thin Disks. II. Disk Variability

abstract: We perform detailed variability analysis of two-dimensional viscous, radiation hydrodynamic numerical simulations of Shakura-Sunyaev thin disks around a stellar mass black hole. Disk models are initialized on both the gas-, as well as radiation-, pressure-dominated branches of the thermal equilibrium curve, with mass accretion rates spanning the range from $\dot{M} = 0.01 L_\mathrm{Edd}/c^2$ to $10 L_\mathrm{Edd}/c^2$. An analysis of temporal variations of the numerically simulated disk reveals multiple robust, coherent oscillations. Considering the local mass flux variability, we find an oscillation occurring at the maximum radial epicyclic frequency, $3.5\times 10^{-3}\,t_\mathrm{g}^{-1}$, a possible signature of a trapped fundamental ${\it g}$-mode. Although present in each of our simulated models, the trapped ${\it g}$-mode feature is most prominent in the gas-pressure-dominated case. The total pressure fluctuations in the disk suggest strong evidence for standing-wave ${\it p}$-modes, some trapped in the inner disk close to the ISCO, others present in the middle/outer parts of the disk. Knowing that the trapped ${\it g}$-mode frequency and maximum radial epicyclic frequency differ by only $0.01\%$ in the case of a non-rotating black hole, we simulated an additional initially gas-pressure-dominated disk with a dimensionless black hole spin parameter $a_* = 0.5$. The oscillation frequency in the spinning black hole case confirms that this oscillation is indeed a trapped ${\it g}$-mode. All the numerical models we report here also show a set of high frequency oscillations at the vertical epicyclic and breathing mode frequencies. The vertical oscillations show a 3:2 frequency ratio with oscillations occurring approximately at the radial epicyclic frequency, which could be of astrophysical importance in observed twin peak, high-frequency quasi-periodic oscillations.


title: Selfie Drone Stick: A Natural Interface for Quadcopter Photography

abstract: A physical selfie stick extends the user's reach, enabling the acquisition of personal photos that include more of the background scene. Similarly, a quadcopter can capture photos from vantage points unattainable by the user; but teleoperating a quadcopter to good viewpoints is a difficult task. This paper presents a natural interface for quadcopter photography, the SelfieDroneStick that allows the user to guide the quadcopter to the optimal vantage point based on the phone's sensors. Users specify the composition of their desired long-range selfies using their smartphone, and the quadcopter autonomously flies to a sequence of vantage points from where the desired shots can be taken. The robot controller is trained from a combination of real-world images and simulated flight data. This paper describes two key innovations required to deploy deep reinforcement learning models on a real robot: 1) an abstract state representation for transferring learning from simulation to the hardware platform, and 2) reward shaping and staging paradigms for training the controller. Both of these improvements were found to be essential in learning a robot controller from simulation that transfers successfully to the real robot.


title: Learning to Communicate and Energize: Modulation, Coding and Multiple Access Designs for Wireless Information-Power Transmission

abstract: The explosion of the number of low-power devices in the next decades calls for a re-thinking of wireless network design, namely, unifying wireless transmission of information and power so as to make the best use of the RF spectrum, radiation, and infrastructure for the dual purpose of communicating and energizing. This paper provides a novel learning-based approach towards such wireless network design. To that end, a parametric model of a practical energy harvester, accounting for various sources of nonlinearities, is proposed using a nonlinear regression algorithm applied over collected real data. Relying on the proposed model, the learning problem of modulation design for Simultaneous Wireless Information-Power Transmission (SWIPT) over a point-to-point link is studied. Joint optimization of the transmitter and the receiver is implemented using Neural Network (NN)-based autoencoders. The results reveal that by increasing the receiver power demand, the baseband transmit modulation constellation converges to an On-Off keying signalling. Utilizing the observations obtained via learning, an algorithmic SWIPT modulation design is proposed. It is observed via numerical results that the performance loss of the proposed modulations are negligible compared to the ones obtained from learning. Extension of the studied problem to learning modulation design for multi-user SWIPT scenarios and coded modulation design for point-to-point SWIPT are considered. The major conclusion of this work is to utilize learning-based results to design non learning-based algorithms, which perform as well. In particular, inspired by the results obtained via learning, an algorithmic approach for coded modulation design is proposed, which performs very close to its learning counterparts, and is significantly superior due to its high real-time adaptability to new system design parameters.


title: Flight Controller Synthesis Via Deep Reinforcement Learning

abstract: Traditional control methods are inadequate in many deployment settings involving control of Cyber-Physical Systems (CPS). In such settings, CPS controllers must operate and respond to unpredictable interactions, conditions, or failure modes. Dealing with such unpredictability requires the use of executive and cognitive control functions that allow for planning and reasoning. Motivated by the sport of drone racing, this dissertation addresses these concerns for state-of-the-art flight control by investigating the use of deep neural networks to bring essential elements of higher-level cognition for constructing low level flight controllers. This thesis reports on the development and release of an open source, full solution stack for building neuro-flight controllers. This stack consists of the methodology for constructing a multicopter digital twin for synthesize the flight controller unique to a specific aircraft, a tuning framework for implementing training environments (GymFC), and a firmware for the world's first neural network supported flight controller (Neuroflight). GymFC's novel approach fuses together the digital twinning paradigm for flight control training to provide seamless transfer to hardware. Additionally, this thesis examines alternative reward system functions as well as changes to the software environment to bridge the gap between the simulation and real world deployment environments. Work summarized in this thesis demonstrates that reinforcement learning is able to be leveraged for training neural network controllers capable, not only of maintaining stable flight, but also precision aerobatic maneuvers in real world settings. As such, this work provides a foundation for developing the next generation of flight control systems.


title: Transactional Smart Contracts in Blockchain Systems

abstract: This paper presents TXSC, a framework that provides smart contract developers with transaction primitives. These primitives allow developers to write smart contracts without the need to reason about the anomalies that can arise due to concurrent smart contract function executions.


title: The O(N) S-matrix Monolith

abstract: We consider the scattering matrices of massive quantum field theories with no bound states and a global $O(N)$ symmetry in two spacetime dimensions. In particular we explore the space of two-to-two S-matrices of particles of mass $m$ transforming in the vector representation as restricted by the general conditions of unitarity, crossing, analyticity and $O(N)$ symmetry. We found a rich structure in that space by using convex maximization and in particular its convex dual minimization problem. At the boundary of the allowed space special geometric points such as vertices were found to correspond to integrable models. The dual convex minimization problem provides a novel and useful approach to the problem allowing, for example, to prove that generically the S-matrices so obtained saturate unitarity and, in some cases, that they are at vertices of the allowed space.