title: Model-based ROC (mROC) curve: examining the effect of case-mix and model calibration on the ROC plot

abstract: The performance of risk prediction models is often characterized in terms of discrimination and calibration. The Receiver Operating Characteristic (ROC) curve is widely used for evaluating model discrimination. When evaluating the performance of a risk prediction model in a new sample, the shape of the ROC curve is affected by both case-mix and the postulated model. Further, compared to discrimination, evaluating calibration has not received the same level of attention. Commonly used methods for model calibration involve subjective specification of smoothing or grouping. Leveraging the familiar ROC framework, we introduce the model-based ROC (mROC) curve to assess the calibration of a pre-specified model in a new sample. mROC curve is the ROC curve that should be observed if a pre-specified model is calibrated in the sample. We show the empirical ROC and mROC curves for a sample converge asymptotically if the model is calibrated in that sample. As a consequence, the mROC curve can be used to assess visually the effect of case-mix and model mis-calibration. Further, we propose a novel statistical test for calibration that does not require any smoothing or grouping. Simulations support the adequacy of the test. A case study puts these developments in a practical context. We conclude that mROC can easily be constructed and used to evaluate the effect of case-mix and model calibration on the ROC plot, thus adding to the utility of ROC curve analysis in the evaluation of risk prediction models. R code for the proposed methodology is provided (https://github.com/msadatsafavi/mROC/).


title: A Bayesian approach to tissue-fraction estimation for oncological PET segmentation

abstract: Tumor segmentation in oncological PET is challenging, a major reason being the partial-volume effects that arise due to low system resolution and finite voxel size. The latter results in tissue-fraction effects, i.e. voxels contain a mixture of tissue classes. Conventional segmentation methods are typically designed to assign each voxel in the image as belonging to a certain tissue class. Thus, these methods are inherently limited in modeling tissue-fraction effects. To address the challenge of accounting for partial-volume effects, and in particular, tissue-fraction effects, we propose a Bayesian approach to tissue-fraction estimation for oncological PET segmentation. Specifically, this Bayesian approach estimates the posterior mean of fractional volume that the tumor occupies within each voxel of the image. The proposed method, implemented using a deep-learning-based technique, was first evaluated using clinically realistic 2-D simulation studies with known ground truth, in the context of segmenting the primary tumor in PET images of patients with lung cancer. The evaluation studies demonstrated that the method accurately estimated the tumor-fraction areas and significantly outperformed widely used conventional PET segmentation methods, including a U-net-based method, on the task of segmenting the tumor. In addition, the proposed method was relatively insensitive to partial-volume effects and yielded reliable tumor segmentation for different clinical-scanner configurations. The method was then evaluated using clinical images of patients with stage IIB/III non-small cell lung cancer from ACRIN 6668/RTOG 0235 multi-center clinical trial. Here, the results showed that the proposed method significantly outperformed all other considered methods and yielded accurate tumor segmentation on patient images with Dice similarity coefficient (DSC) of 0.82 (95 % CI: [0.78, 0.86]).


title: From Quantum Foundations of Quantum Field Theory, String Theory and Quantum Gravity to Dark Matter and Dark Energy

abstract: We review our recent work on quantum foundations of quantum mechanics, quantum field theory and quantum gravity (formulated as metastring theory) and various implications for the problems of dark matter and dark energy. The first point concerns the new understanding of quantum theory via the concept of quantum (modular) spacetime endowed with manifest non-locality that is consistent with causality. This view implies the consistency of the fundamental length and Lorentz symmetry, based on the principle of relative (observer dependent) locality. The geometry of such quantum spacetime is encoded in the new concept of Born geometry. This in turn leads to a novel understanding of quantum field theory in a manifestly bi-local representation endowed with metaparticle quanta. A fully dynamical quantum spacetime, with a dynamical Born geometry, leads to quantum gravity (a quantum theory of dynamical geometry of quantum non-locality) in the guise of metastring theory. This generic formulation of string theory implies a radiatively stable positive cosmological constant (viewed as a curvature of the dual spacetime) as a model of dark energy in the observed classical spacetime, as well as metaparticle quanta (the zero modes of the metastring) as the natural quanta of dark matter in this approach.


title: The Kim-Pillay theorem for Abstract Elementary Categories

abstract: We introduce the framework of AECats (abstract elementary categories), generalising both the category of models of some first-order theory and the category of subsets of models. Any AEC and any compact abstract theory ("cat", as introduced by Ben-Yaacov) forms an AECat. In particular, we find applications in positive logic and continuous logic: the category of (subsets of) models of a positive or continuous theory is an AECat. The Kim-Pillay theorem for first-order logic characterises simple theories by the properties dividing independence has. We prove a version of the Kim-Pillay theorem for AECats with the amalgamation property, generalising the first-order version and existing versions for positive logic.


title: Electron confinement by laser-driven azimuthal magnetic fields during direct laser acceleration

abstract: A laser-driven azimuthal plasma magnetic field is known to facilitate electron energy gain from the irradiating laser pulse. The enhancement is due to changes in the orientation between the laser electric field and electron velocity caused by magnetic field deflections. Transverse electron confinement is critical for realizing this concept experimentally. We find that the phase velocity of the laser pulse has a profound impact on the transverse size of electron trajectories. The transverse size remains constant below a threshold energy that depends on the degree of the superluminosity and it increases with the electron energy above the threshold. This increase can cause electron losses in tightly focused laser pulses. We show using 3D particle-in-cell simulations that the electron energy gain can be significantly increased by increasing laser power at fixed intensity due to the increased electron confinement. This finding makes a strong case for designing experiments at multi-PW laser facilities.


title: Learning Cross-domain Generalizable Features by Representation Disentanglement

abstract: Deep learning models exhibit limited generalizability across different domains. Specifically, transferring knowledge from available entangled domain features(source/target domain) and categorical features to new unseen categorical features in a target domain is an interesting and difficult problem that is rarely discussed in the current literature. This problem is essential for many real-world applications such as improving diagnostic classification or prediction in medical imaging. To address this problem, we propose Mutual-Information-based Disentangled Neural Networks (MIDNet) to extract generalizable features that enable transferring knowledge to unseen categorical features in target domains. The proposed MIDNet is developed as a semi-supervised learning paradigm to alleviate the dependency on labeled data. This is important for practical applications where data annotation requires rare expertise as well as intense time and labor. We demonstrate our method on handwritten digits datasets and a fetal ultrasound dataset for image classification tasks. Experiments show that our method outperforms the state-of-the-art and achieve expected performance with sparsely labeled data.


title: Mirror, mirror: Landau-Zener-Stuckelberg-Majorana interferometry of a superconducting qubit in front of a mirror

abstract: We investigate the Landau-Zener-Stuckelberg-Majorana interferometry of a superconducting qubit in a semi-infinite transmission line terminated by a mirror. The transmon-type qubit is at the node of the resonant electromagnetic (EM) field, hiding from the EM field. "Mirror, mirror" briefly describes this system, because the qubit acts as another mirror. We modulate the resonant frequency of the qubit by applying a sinusoidal flux pump. We probe the spectroscopy by measuring the reflection coefficient of a weak probe in the system. Remarkable interference patterns emerge in the spectrum, which can be interpreted as multi-photon resonances in the dressed qubit. Our calculations agree well with the experiments.


title: Velocity and size measurement of droplets from an ultrasonic spray coater using Photon Correlation Spectroscopy and Turbidimetry

abstract: We have developed a combination of light scattering techniques to study and characterize droplets of a ultrasonic spray printer or coater in flight. For this economically relevant printer there is so far no reliable technique to systematically adjust the experimental parameters. We have combined photon correlation spectroscopy and turbidimetry to determine size and speed of the droplets depending on parameters of the printing process as viscosity, concentration and speed of the driving shroud gas. Our method allows to predetermine these parameters to control the properties of the printed films as e.g. thickness from tens of nanometers to micrometers.


title: Complexes of Tournaments, Directionality Filtrations and Persistent Homology

abstract: Complete digraphs are referred to in the combinatorics literature as tournaments. We consider a family of semi-simplicial complexes, that we refer to as "tournaplexes", whose simplices are tournaments. In particular, given a digraph $\mathcal{G}$, we associate with it a "flag tournaplex" which is a tournaplex containing the directed flag complex of $\mathcal{G}$, but also the geometric realisation of cliques that are not directed. We define several types of filtrations on tournaplexes, and exploiting persistent homology, we observe that flag tournaplexes provide finer means of distinguishing graph dynamics than the directed flag complex. We then demonstrate the power of these ideas by applying them to graph data arising from the Blue Brain Project's digital reconstruction of a rat's neocortex.


title: A general variational formulation for relativistic mechanics based on fundamentals of differential geometry

abstract: The first part of this article develops a variational formulation for relativistic mechanics. The results are established through standard tools of variational analysis and differential geometry. The novelty here is that the main motion manifold has a $n+1$ dimensional range. It is worth emphasizing in a first approximation we have neglected the self-interaction energy part. In its second part, this article develops some formalism concerning the causal structure in a general space-time manifold. Finally, the last article section presents a result concerning the existence of a generalized solution for the world sheet manifold variational formulation.


title: Unsafe At Any Level: NHTSA's levels of automation are a liability for autonomous vehicle design and regulation

abstract: Walter Huang, a 38-year-old Apple Inc. engineer, died on March 23, 2018, after his Tesla Model X crashed into a highway barrier in Mountain View, California. Tesla immediately disavowed responsibility for the accident. "The fundamental premise of both moral and legal liability is a broken promise, and there was none here: [Mr. Huang] was well aware that the Autopilot was not perfect [and the] only way for this accident to have occurred is if Mr. Huang was not paying attention to the road, despite the car providing multiple warnings to do so." This is the standard response from Tesla and Uber, the manufacturers of the automated vehicles involved in the six fatal accidents to date: the automated vehicle isn't perfect, the driver knew it wasn't perfect, and if only the driver had been paying attention and heeded the vehicle's warnings, the accident would never have occurred. However, as researchers focused on human-automation interaction in aviation and military operations, we cannot help but wonder if there really are no broken promises and no legal liabilities. Science has a critical role in determining legal liability, and courts appropriately rely on scientists and engineers to determine whether an accident, or harm, was foreseeable. Specifically, a designer could be found liable if, at the time of the accident, scientists knew there was a systematic relationship between the accident and the designer's untaken precaution. Nearly 70 years of research provides an undeniable answer: It is insufficient, inappropriate, and dangerous to automate everything you can and leave the rest to the human. There is a systematic relationship between the design of automated vehicles and the types of accidents that are occurring now and will inevitably continue to occur in the future. These accidents were not unforeseeable and the drivers were not exclusively to blame.


title: Modelling techniques for biomolecular networks

abstract: First we shortly review the different kinds of network modelling methods for systems biology with an emphasis on the different subtypes of logical models, which we review in more detail. Then we show the advantages of Boolean networks models over more mechanistic modelling types like differential equation techniques. Then follows an overlook about connections between different kinds of models and how they can be converted to each other. We also give a short overview about the mathematical frameworks for modelling of logical networks and list available software packages for logical modelling. Then we give an overview about the available standards and ontologies for storing such logical systems biology models and their results. In the end we give a short review about the difference between quantitative and qualitative models and describe the mathematics that specifically deals with qualitative modelling.


title: Mass spectrometry for semi-experimental protein structure determination and modeling

abstract: The structure of proteins is essential for its function. The determination of protein structures is possible by experimental or predicted by computational methods, but also a combination of both approaches is possible. Here, first an overview about experimental structure determination methods with their pros and cons is given. Then we describe how mass spectrometry is useful for semi-experimental integrative protein structure determination. We review the methodology and describe software programs supporting such integrated protein structure prediction approaches, making use of distance constraints got from mass spectrometry cross-linking experiments


title: XLMOD: Cross-linking and chromatography derivatization reagents ontology

abstract: Mass spectrometry has experienced a rapid development since its first application for protein analysis in the 1980s. While the most common use of mass spectrometry for protein analysis is identification and quantification workflows on peptides (digested from their parent protein), there is also a rapidly growing use of mass spectrometry for structural proteomics. One example is the analysis of cross-linked proteins that can give valuable structural information, complementing the information gained by classical protein structure determination methods, useful for integrated methods of structure determination and modeling. For a broad and reproducible application of cross-linking mass spectrometry a standardized representation of cross-linking experimental results is necessary. This paper describes the developing and release of the xlmod ontology from the HUPO-PSI. xlmod contains terms for the description of reagents used in cross-linking experiments and of cross-linker related chemical modifications together with their main properties relevant for planning and performing cross-linking experiments. We also describe how xlmod is used within the new release of HUPO-PSI-s mzIdentML data standard, for reporting the used cross-linking reagents and results in a consistent manner. In addition xlmod contains terms for GC-MS and LC-MS derivatization reagents for specifying them in the upcoming mzTab-M and mzTab-L formats.


title: Graph Neural Networks Meet Neural-Symbolic Computing: A Survey and Perspective

abstract: Neural-symbolic computing has now become the subject of interest of both academic and industry research laboratories. Graph Neural Networks (GNN) have been widely used in relational and symbolic domains, with widespread application of GNNs in combinatorial optimization, constraint satisfaction, relational reasoning and other scientific domains. The need for improved explainability, interpretability and trust of AI systems in general demands principled methodologies, as suggested by neural-symbolic computing. In this paper, we review the state-of-the-art on the use of GNNs as a model of neural-symbolic computing. This includes the application of GNNs in several domains as well as its relationship to current developments in neural-symbolic computing.


title: Bipartisan Paxos: A Modular State Machine Replication Protocol

abstract: There is no shortage of state machine replication protocols. From Generalized Paxos to EPaxos, a huge number of replication protocols have been proposed that achieve high throughput and low latency. However, these protocols all have two problems. First, they do not scale. Many protocols actually slow down when you scale them, instead of speeding up. For example, increasing the number of MultiPaxos acceptors increases quorum sizes and slows down the protocol. Second, they are too complicated. This is not a secret; state machine replication is notoriously difficult to understand. In this paper, we tackle both problems with a single solution: modularity. We present Bipartisan Paxos (BPaxos), a modular state machine replication protocol. Modularity yields high throughput via scaling. We note that while many replication protocol components do not scale, some do. By modularizing BPaxos, we are able to disentangle the two and scale the bottleneck components to increase the protocol's throughput. Modularity also yields simplicity. BPaxos is divided into a number of independent modules that can be understood and proven correct in isolation.


title: A class of equations with three solutions

abstract: Here is one of the results obtained in this paper: Let $\Omega\subset {\bf R}^n$ be a smooth bounded domain, let $q>1$, with $q<{{n+2}\over {n-2}}$ if $n\geq 3$ and let $\lambda_1$ be the first eigenvalue of the problem $$\cases{-\Delta u=\lambda u & in $\Omega$ \cr & \cr u=0 & on $\partial\Omega$\ .\cr}$$ Then, for every $\lambda>\lambda_1$ and for every convex set $S\subseteq L^{\infty}(\Omega)$ dense in $L^2(\Omega)$, there exists $\alpha\in S$ such that the problem $$\cases{-\Delta u=\lambda(u^+-(u^+)^q)+\alpha(x) & in $\Omega$ \cr & \cr u=0 & on $\partial\Omega$\cr}$$ has at least three weak solutions, two of which are global minima in $H^1_0(\Omega)$ of the functional $$u\to {{1}\over {2}}\int_{\Omega}|\nabla u(x)|^2dx-\lambda\int_{\Omega}\left ({{1}\over {2}}|u^+(x)|^2-{{1}\over {q+1}}|u^+(x)|^{q+1}\right )dx-\int_{\Omega}\alpha(x)u(x)dx\ $$ where $u^+=\max\{u,0\}$.


title: Multi-Level Optimal Power Flow Solver in Large Distribution Networks

abstract: Solving optimal power flow (OPF) problems for large distribution networks incurs high computational complexity. We consider a large multi-phase distribution network of tree topology with a deep penetration of active devices. We divide the network into collaborating areas featuring subtree topology and subareas featuring subsubtree topology. We design a multi-level implementation of the primal-dual gradient algorithm to solve the voltage regulation OPF problems while preserving nodal voltage information and topological information within areas and subareas. Numerical results on a 4,521-node system verify that the proposed algorithm can significantly improve the computational speed without compromising any optimality.


title: Asymptotic Smiles for an Affine Jump-Diffusion Model

abstract: In this paper, we study the asymptotic behaviors of implied volatility of an affine jump-diffusion model. Let log stock price under risk-neutral measure follow an affine jump-diffusion model, we show that an explicit form of moment generating function for log stock price can be obtained by solving a set of ordinary differential equations. A large-time large deviation principle for log stock price is derived by applying the G\"{a}rtner-Ellis theorem. We characterize the asymptotic behaviors of the implied volatility in the large-maturity and large-strike regime using rate function in the large deviation principle. The asymptotics of the Black-Scholes implied volatility for fixed-maturity, large-strike and fixed-maturity, small-strike regimes are also studied. Numerical results are provided to validate the theoretical work.


title: Differentiating through the Fr\'echet Mean

abstract: Recent advances in deep representation learning on Riemannian manifolds extend classical deep learning operations to better capture the geometry of the manifold. One possible extension is the Fr\'echet mean, the generalization of the Euclidean mean; however, it has been difficult to apply because it lacks a closed form with an easily computable derivative. In this paper, we show how to differentiate through the Fr\'echet mean for arbitrary Riemannian manifolds. Then, focusing on hyperbolic space, we derive explicit gradient expressions and a fast, accurate, and hyperparameter-free Fr\'echet mean solver. This fully integrates the Fr\'echet mean into the hyperbolic neural network pipeline. To demonstrate this integration, we present two case studies. First, we apply our Fr\'echet mean to the existing Hyperbolic Graph Convolutional Network, replacing its projected aggregation to obtain state-of-the-art results on datasets with high hyperbolicity. Second, to demonstrate the Fr\'echet mean's capacity to generalize Euclidean neural network operations, we develop a hyperbolic batch normalization method that gives an improvement parallel to the one observed in the Euclidean setting.


title: Three-dimensional matching is NP-Hard

abstract: The standard proof of NP-Hardness of 3DM provides a power-$4$ reduction of 3SAT to 3DM. In this note, we provide a linear-time reduction. Under the exponential time hypothesis, this reduction improves the runtime lower bound from $2^{o(\sqrt[4]{m})}$ (under the standard reduction) to $2^{o(m)}$.


title: The Weil-Petersson gradient flow of renormalized volume and 3-dimensional convex cores

abstract: In this paper, we use the Weil-Petersson gradient flow for renormalized volume to study the space $CC(N;S,X)$ of convex cocompact hyperbolic structures on the relatively acylindrical 3-manifold $(N;S)$. Among the cases of interest are the deformation space of an acylindrical manifold and the Bers slice of quasi-Fuchsian space associated to a fixed surface. To treat the possibility of degeneration along flow-lines to peripherally cusped structures, we introduce a surgery procedure to yield a surgered gradient flow that limits to the unique structure $M_{\rm geod} \in CC(N;S,X)$ with totally geodesic convex core boundary facing $S$. Analyzing the geometry of structures along a flow line, we show that if $V_R(M)$ is the renormalized volume of $M$, then $V_R(M)-V_R(M_{\rm geod})$ is bounded below by a linear function of the Weil-Petersson distance $d_{\rm WP}(\partial_c M, \partial_c M_{\rm geod})$, with constants depending only on the topology of $S$. The surgered flow gives a unified approach to a number of problems in the study of hyperbolic 3-manifolds, providing new proofs and generalizations of well-known theorems such as Storm's result that $M_{\rm geod}$ has minimal volume for $N$ acylindrical and the second author's result comparing convex core volume and Weil-Petersson distance for quasifuchsian manifolds.


title: Transversal local rigidity of discrete Abelian actions on Heisenberg nilmanifolds

abstract: In this paper we prove a perturbative result for a class of $\mathbb Z^2$ actions on Heisenberg nilmanifolds, which have Diophantine properties. Along the way we prove cohomological rigidity and obtain a tame splitting for the cohomology with coefficients in smooth vector fields for such actions.


title: Yet another generalization of the Fuglede-Putnam theorem to unbounded operators

abstract: In this note, we give the most natural (perhaps the simplest ever) generalization of the Fuglede-Putnam theorem where all operators involved are unbounded.


title: Stripes on finite domains: Why the zigzag instability is only a partial story

abstract: Stationary periodic patterns are widespread in natural sciences, ranging from nano-scale electrochemical and amphiphilic systems to mesoscale fluid, chemical and biological media and to macro-scale vegetation and cloud patterns. Their formation is usually due to a primary symmetry breaking of a uniform state to stripes, often followed by secondary instabilities to form zigzag and labyrinthine patterns. These secondary instabilities are well studied under idealized conditions of an infinite domain, however, on finite domains, the situation is more subtle since the unstable modes depend also on boundary conditions. Using two prototypical models, the Swift-Hohenberg equation and the forced complex Ginzburg-Landau equation, we consider bounded domains with no flux boundary conditions transversal to the stripes, and reveal a distinct mixed-mode instability that lies in between the classical zigzag and the Eckhaus lines. This explains the stability of stripes in the mildly zigzag unstable regime, and, after crossing the mixed-mode line, the evolution of zigzag stripes in the bulk of the domain and the formation of defects near the boundaries. The results are of particular importance for problems with large time scale separation, such as bulk-heterojunction deformations in organic photovoltaic and vegetation in semi-arid regions, where early temporal transients may play an important role.


title: Finite $N$ unitary matrix model

abstract: We consider one-plaquette unitary matrix model at finite $N$ using exact expression of the partition function for both SU($N$) and U($N$) groups.


title: Robust Robotic Pouring using Audition and Haptics

abstract: Robust and accurate estimation of liquid height lies as an essential part of pouring tasks for service robots. However, vision-based methods often fail in occluded conditions while audio-based methods cannot work well in a noisy environment. We instead propose a multimodal pouring network (MP-Net) that is able to robustly predict liquid height by conditioning on both audition and haptics input. MP-Net is trained on a self-collected multimodal pouring dataset. This dataset contains 300 robot pouring recordings with audio and force/torque measurements for three types of target containers. We also augment the audio data by inserting robot noise. We evaluated MP-Net on our collected dataset and a wide variety of robot experiments. Both network training results and robot experiments demonstrate that MP-Net is robust against noise and changes to the task and environment. Moreover, we further combine the predicted height and force data to estimate the shape of the target container.


title: Calibrated Prediction with Covariate Shift via Unsupervised Domain Adaptation

abstract: Reliable uncertainty estimates are an important tool for helping autonomous agents or human decision makers understand and leverage predictive models. However, existing approaches to estimating uncertainty largely ignore the possibility of covariate shift--i.e., where the real-world data distribution may differ from the training distribution. As a consequence, existing algorithms can overestimate certainty, possibly yielding a false sense of confidence in the predictive model. We propose an algorithm for calibrating predictions that accounts for the possibility of covariate shift, given labeled examples from the training distribution and unlabeled examples from the real-world distribution. Our algorithm uses importance weighting to correct for the shift from the training to the real-world distribution. However, importance weighting relies on the training and real-world distributions to be sufficiently close. Building on ideas from domain adaptation, we additionally learn a feature map that tries to equalize these two distributions. In an empirical evaluation, we show that our proposed approach outperforms existing approaches to calibrated prediction when there is covariate shift.


title: An Evaluation of Knowledge Graph Embeddings for Autonomous Driving Data: Experience and Practice

abstract: The autonomous driving (AD) industry is exploring the use of knowledge graphs (KGs) to manage the vast amount of heterogeneous data generated from vehicular sensors. The various types of equipped sensors include video, LIDAR and RADAR. Scene understanding is an important topic in AD which requires consideration of various aspects of a scene, such as detected objects, events, time and location. Recent work on knowledge graph embeddings (KGEs) - an approach that facilitates neuro-symbolic fusion - has shown to improve the predictive performance of machine learning models. With the expectation that neuro-symbolic fusion through KGEs will improve scene understanding, this research explores the generation and evaluation of KGEs for autonomous driving data. We also present an investigation of the relationship between the level of informational detail in a KG and the quality of its derivative embeddings. By systematically evaluating KGEs along four dimensions -- i.e. quality metrics, KG informational detail, algorithms, and datasets -- we show that (1) higher levels of informational detail in KGs lead to higher quality embeddings, (2) type and relation semantics are better captured by the semantic transitional distance-based TransE algorithm, and (3) some metrics, such as coherence measure, may not be suitable for intrinsically evaluating KGEs in this domain. Additionally, we also present an (early) investigation of the usefulness of KGEs for two use-cases in the AD domain.


title: Robust Model Predictive Control for Nonlinear Systems Using Convex Restriction

abstract: We present an algorithm for robust model predictive control with consideration of uncertainty and safety constraints. Our framework considers a nonlinear dynamical system subject to disturbances from an unknown but bounded uncertainty set. By viewing the system as a fixed point of an operator acting over trajectories, we propose a convex condition on control actions that guarantee safety against the uncertainty set. The proposed condition guarantees that all realizations of the state trajectories satisfy safety constraints. Our algorithm solves a sequence of convex quadratic constrained optimization problems of size n*N, where n is the number of states, and N is the prediction horizon in the model predictive control problem. Compared to existing methods, our approach solves convex problems while guaranteeing that all realizations of uncertainty set do not violate safety constraints. Moreover, we consider the implicit time-discretization of system dynamics to increase the prediction horizon and enhance computational accuracy. Numerical simulations for vehicle navigation demonstrate the effectiveness of our approach.


title: Machine Learning for Achieving Bose-Einstein Condensation of Thulium Atoms

abstract: Bose-Einstein condensation (BEC) is a powerful tool for a wide range of research activities, a large fraction of which are related to quantum simulations. Various problems may benefit from different atomic species, but cooling down novel species interesting for quantum simulations to BEC temperatures requires a substantial amount of optimization and is usually considered as a hard experimental task. In this work, we implemented the Bayesian machine learning technique to optimize the evaporative cooling of thulium atoms and achieved BEC in an optical dipole trap operating near 532 nm. The developed approach could be used to cool down other novel atomic species to quantum degeneracy without additional studies of their properties.


title: Multirate Timestepping for the Incompressible Navier-Stokes Equations in Overlapping Grids

abstract: We develop a multirate timestepper for semi-implicit solutions of the unsteady incompressible Navier-Stokes equations (INSE) based on a recently-developed multidomain spectral element method (SEM). For {\em incompressible} flows, multirate timestepping (MTS) is particularly challenging because of the tight coupling implied by the incompressibility constraint, which manifests as an elliptic subproblem for the pressure at each timestep. The novelty of our approach stems from the development of a stable overlapping Schwarz method applied directly to the Navier-Stokes equations, rather than to the convective, viscous, and pressure substeps that are at the heart of most INSE solvers. Our MTS approach is based on a predictor-corrector (PC) strategy that preserves the temporal convergence of the underlying semi-implicit timestepper. We present numerical results demonstrating that this approach scales to an arbitrary number of overlapping grids, accurately models complex turbulent flow phenomenon, and improves computational efficiency in comparison to singlerate timestepping-based calculations.


title: Information cartography in association rule mining

abstract: Association Rule Mining is a machine learning method for discovering the interesting relations between the attributes in a huge transaction database. Typically, algorithms for Association Rule Mining generate a huge number of association rules, from which it is hard to extract structured knowledge and present this automatically in a form that would be suitable for the user. Recently, an information cartography has been proposed for creating structured summaries of information and visualizing with methodology called "metro maps". This was applied to several problem domains, where pattern mining was necessary. The aim of this study is to develop a method for automatic creation of metro maps of information obtained by Association Rule Mining and, thus, spread its applicability to the other machine learning methods. Although the proposed method consists of multiple steps, its core presents metro map construction that is defined in the study as an optimization problem, which is solved using an evolutionary algorithm. Finally, this was applied to four well-known UCI Machine Learning datasets and one sport dataset. Visualizing the resulted metro maps not only justifies that this is a suitable tool for presenting structured knowledge hidden in data, but also that they can tell stories to users.


title: Self-testing of physical theories, or, is quantum theory optimal with respect to some information-processing task?

abstract: Self-testing usually refers to the task of taking a given set of observed correlations that are assumed to arise via a process that is accurately described by quantum theory, and trying to infer the quantum state and measurements. In other words it is concerned with the question of whether we can tell what quantum black-box devices are doing by looking only at their input-output behaviour and is known to be possible in several cases. Here we introduce a more general question: is it possible to self-test a theory, and, in particular, quantum theory? More precisely, we ask whether within a particular causal structure there are tasks that can only be performed in theories that have the same correlations as quantum mechanics in any scenario. We present a candidate task for such a correlation self-test and analyse it in a range of generalised probabilistic theories (GPTs), showing that none of these perform better than quantum theory. A generalisation of our results showing that all non-quantum GPTs are strictly inferior to quantum mechanics for this task would point to a new way to axiomatise quantum theory, and enable an experimental test that simultaneously rules out such GPTs.


title: Topological transition of a non-Markovian dissipative quantum walk

abstract: We extend non-Hermitian topological quantum walks on a Su-Schrieffer-Heeger (SSH) lattice [M. S. Rudner and L. Levitov, Phys. Rev. Lett. 102, 065703 (2009)] to the case of non-Markovian evolution. This non-Markovian model is established by coupling each unit cell in the SSH lattice to a reservoir formed by a quasi-continuum of levels. We find a topological transition in this model even in the case of non-Markovian evolution, where the walker may visit the reservoir and return to the SSH lattice at a later time. The existence of a topological transition does, however, depend on the low-frequency properties of the reservoir, characterized by a spectral density $J(\epsilon)\propto |\epsilon|^\alpha$. In particular, we find a robust topological transition for a sub-Ohmic ($\alpha<1$) and Ohmic ($\alpha=1$) reservoir, but no topological transition for a super-Ohmic ($\alpha>1$) reservoir. This behavior is directly related to the well-known localization transition for the spin-boson model. We confirm the presence of non-Markovian dynamics by explicitly evaluating a measure of Markovianity for this model.


title: Emotion Recognition System from Speech and Visual Information based on Convolutional Neural Networks

abstract: Emotion recognition has become an important field of research in the human-computer interactions domain. The latest advancements in the field show that combining visual with audio information lead to better results if compared to the case of using a single source of information separately. From a visual point of view, a human emotion can be recognized by analyzing the facial expression of the person. More precisely, the human emotion can be described through a combination of several Facial Action Units. In this paper, we propose a system that is able to recognize emotions with a high accuracy rate and in real time, based on deep Convolutional Neural Networks. In order to increase the accuracy of the recognition system, we analyze also the speech data and fuse the information coming from both sources, i.e., visual and audio. Experimental results show the effectiveness of the proposed scheme for emotion recognition and the importance of combining visual with audio data.


title: Random geometries for optimal control PDE problems based on fictitious domain FEMS and cut elements

abstract: This work investigates an elliptic optimal control problem defined on uncertain domains and discretized by a fictitious domain finite element method and cut elements. Key ingredients of the study are to manage cases considering the usually computationally "forbidden" combination of poorly conditioned equation system matrices due to challenging geometries, optimal control searches with iterative methods, slow convergence to system solutions on deterministic and non--deterministic level, and expensive remeshing due to geometrical changes. We overcome all these difficulties, utilizing the advantages of proper preconditioners adapted to unfitted mesh methods, improved types of Monte Carlo methods, and mainly employing the advantages of embedded FEMs, based on a fixed background mesh computed once even if geometrical changes are taking place. The sensitivity of the control problem is introduced in terms of random domains, employing a Quasi-Monte Carlo method and emphasizing a deterministic target state. The variational discretization concept is adopted, optimal error estimates for the state, adjoint state, and control are derived that confirm the efficiency of the cut finite element method in challenging geometries. The performance of a multigrid scheme especially developed for unfitted finite element discretizations adapted to the optimal control problem is also tested. Some fundamental preconditioners are applied to the arising sparse linear systems coming from the discretization of the state and adjoint state variational forms in the spatial domain. The corresponding convergence rates along with the quality of the prescribed preconditioners are verified by numerical examples.


title: Clinical Text Summarization with Syntax-Based Negation and Semantic Concept Identification

abstract: In the era of clinical information explosion, a good strategy for clinical text summarization is helpful to improve the clinical workflow. The ideal summarization strategy can preserve important information in the informative but less organized, ill-structured clinical narrative texts. Instead of using pure statistical learning approaches, which are difficult to interpret and explain, we utilized knowledge of computational linguistics with human experts-curated biomedical knowledge base to achieve the interpretable and meaningful clinical text summarization. Our research objective is to use the biomedical ontology with semantic information, and take the advantage from the language hierarchical structure, the constituency tree, in order to identify the correct clinical concepts and the corresponding negation information, which is critical for summarizing clinical concepts from narrative text. We achieved the clinically acceptable performance for both negation detection and concept identification, and the clinical concepts with common negated patterns can be identified and negated by the proposed method.


title: A Stochastic Covariance Shrinkage Approach in Ensemble Transform Kalman Filtering

abstract: The Ensemble Kalman Filters (EnKF) employ a Monte-Carlo approach to represent covariance information, and are affected by sampling errors in operational settings where the number of model realizations is much smaller than the model state dimension. To alleviate the effects of these errors EnKF relies on model-specific heuristics such as covariance localization, which takes advantage of the spatial locality of correlations among the model variables. This work proposes an approach to alleviate sampling errors that utilizes a locally averaged-in-time dynamics of the model, described in terms of a climatological covariance of the dynamical system. We use this covariance as the target matrix in covariance shrinkage methods, and develop a stochastic covariance shrinkage approach where synthetic ensemble members are drawn to enrich both the ensemble subspace and the ensemble transformation. We additionally provide for a way in which this methodology can be localized similar to the state-of-the-art LETKF method, and that for a certain model setup, our methodology significantly outperforms it.


title: Survival Cluster Analysis

abstract: Conventional survival analysis approaches estimate risk scores or individualized time-to-event distributions conditioned on covariates. In practice, there is often great population-level phenotypic heterogeneity, resulting from (unknown) subpopulations with diverse risk profiles or survival distributions. As a result, there is an unmet need in survival analysis for identifying subpopulations with distinct risk profiles, while jointly accounting for accurate individualized time-to-event predictions. An approach that addresses this need is likely to improve characterization of individual outcomes by leveraging regularities in subpopulations, thus accounting for population-level heterogeneity. In this paper, we propose a Bayesian nonparametrics approach that represents observations (subjects) in a clustered latent space, and encourages accurate time-to-event predictions and clusters (subpopulations) with distinct risk profiles. Experiments on real-world datasets show consistent improvements in predictive performance and interpretability relative to existing state-of-the-art survival analysis models.


title: Current Oscillations in Quasi-2D Charge-Density-Wave 1T-TaS2 Devices: Revisiting the "Narrow Band Noise" Concept

abstract: We report on current oscillations in quasi two-dimensional (2D) 1T-TaS2 charge-density-wave devices. The MHz-frequency range of the oscillations and the linear dependence of the frequency of the oscillations on the current resemble closely the "narrow band noise," which was often observed in the classical bulk quasi-one-dimensional (1D) trichalcogenide charge-density-wave materials. In bulk quasi-1D materials, the "narrow band noise" was interpreted as a direct evidence of the charge-density-wave sliding. Despite the similarities, we argue that the nature of the MHz oscillations in quasi-2D 1T-TaS2 is different from the "narrow band noise." Analysis of the biasing conditions and current indicate that the observed oscillations are related to the current instabilities due to the voltage-induced transition from the nearly commensurate to incommensurate charge-density-wave phase.


title: Exploring the Bayesian parameter estimation of binary black holes with LISA

abstract: The space-based gravitational wave detector LISA will observe mergers of massive black hole binary systems (MBHBs) to cosmological distances, as well as inspiralling stellar-origin (or stellar-mass) binaries (SBHBs) years before they enter the LIGO/Virgo band. Much remains to be explored for the parameter recovery of both classes of systems. Previous MBHB analyses relied on inspiral-only signals and/or a simplified Fisher matrix analysis, while SBHBs have not yet been extensively analyzed with Bayesian methods. We accelerate likelihood computations by (i) using a Fourier-domain response of the LISA instrument, (ii) using a reduced-order model for non-spinning waveforms that include a merger-ringdown and higher harmonics, (iii) setting the noise realization to zero and computing overlaps in the amplitude/phase representation. We present the first simulations of Bayesian inference for the parameters of massive black hole systems including consistently the merger and ringdown of the signal, as well as higher harmonics. We clarify the roles of LISA response time and frequency dependencies in breaking degeneracies and illustrate how degeneracy breaking unfolds over time. We also find that restricting the merger-dominated signal to its dominant harmonic can make the extrinsic likelihood very degenerate. Including higher harmonics proves to be crucial to break degeneracies and considerably improves the localization of the source, with a surviving bimodality in the sky position. We also present simulations of Bayesian inference for the extrinsic parameters of SBHBs, and show that although unimodal, their posterior distributions can have non-Gaussian features.


title: Multiplicities, pictographs, and volumes

abstract: The present contribution is the written counterpart of a talk given in Yerevan at the SQS'2019 International Workshop "Supersymmetries and Quantum Symmetries" (SQS'2019, 26 August - August 31, 2019). After a short presentation of various pictographs (O-blades, metric honeycombs) that one can use in order to calculate SU(n) multiplicities (Littlewood-Richardson coefficients, Kostka numbers), we briefly discuss the semi-classical limit of these multiplicities in relation with the Horn and Schur volume functions and with the so-called Rn-polynomials that enter the expression of volume functions. For n < 7 the decomposition of the Rn-polynomials on Lie group characters is already known, the case n=7 is obtained here.


title: Contextual-Bandit Based Personalized Recommendation with Time-Varying User Interests

abstract: A contextual bandit problem is studied in a highly non-stationary environment, which is ubiquitous in various recommender systems due to the time-varying interests of users. Two models with disjoint and hybrid payoffs are considered to characterize the phenomenon that users' preferences towards different items vary differently over time. In the disjoint payoff model, the reward of playing an arm is determined by an arm-specific preference vector, which is piecewise-stationary with asynchronous and distinct changes across different arms. An efficient learning algorithm that is adaptive to abrupt reward changes is proposed and theoretical regret analysis is provided to show that a sublinear scaling of regret in the time length $T$ is achieved. The algorithm is further extended to a more general setting with hybrid payoffs where the reward of playing an arm is determined by both an arm-specific preference vector and a joint coefficient vector shared by all arms. Empirical experiments are conducted on real-world datasets to verify the advantages of the proposed learning algorithms against baseline ones in both settings.


title: Decision Trees for Decision-Making under the Predict-then-Optimize Framework

abstract: We consider the use of decision trees for decision-making problems under the predict-then-optimize framework. That is, we would like to first use a decision tree to predict unknown input parameters of an optimization problem, and then make decisions by solving the optimization problem using the predicted parameters. A natural loss function in this framework is to measure the suboptimality of the decisions induced by the predicted input parameters, as opposed to measuring loss using input parameter prediction error. This natural loss function is known in the literature as the Smart Predict-then-Optimize (SPO) loss, and we propose a tractable methodology called SPO Trees (SPOTs) for training decision trees under this loss. SPOTs benefit from the interpretability of decision trees, providing an interpretable segmentation of contextual features into groups with distinct optimal solutions to the optimization problem of interest. We conduct several numerical experiments on synthetic and real data including the prediction of travel times for shortest path problems and predicting click probabilities for news article recommendation. We demonstrate on these datasets that SPOTs simultaneously provide higher quality decisions and significantly lower model complexity than other machine learning approaches (e.g., CART) trained to minimize prediction error.


title: Testing a quantum annealer as a quantum thermal sampler

abstract: Motivated by recent experiments in which specific thermal properties of complex many-body systems were successfully reproduced on a commercially available quantum annealer, we examine the extent to which quantum annealing hardware can reliably sample from the thermal state in a specific basis associated with a target quantum Hamiltonian. We address this question by studying the diagonal thermal properties of the canonical one-dimensional transverse-field Ising model on a D-Wave 2000Q quantum annealing processor. We find that the quantum processor fails to produce the correct expectation values predicted by Quantum Monte Carlo. Comparing to master equation simulations, we find that this discrepancy is best explained by how the measurements at finite transverse fields are enacted on the device. Specifically, measurements at finite transverse field require the system to be quenched from the target Hamiltonian to a Hamiltonian with negligible transverse field, and this quench is too slow. The limitations imposed by such hardware make it an unlikely candidate for thermal sampling, and it remains an open question what thermal expectation values can be robustly estimated in general for arbitrary quantum many-body systems.


title: Recursive max-linear models with propagating noise

abstract: Recursive max-linear vectors model causal dependence between node variables by a structural equation model, expressing each node variable as a max-linear function of its parental nodes in a directed acyclic graph (DAG) and some exogenous innovation. For such a model, there exists a unique minimum DAG, represented by the Kleene star matrix of its edge weight matrix, which identifies the model and can be estimated. For a more realistic statistical modeling we introduce some random observational noise. A probabilistic analysis of this new noisy model reveals that the unique minimum DAG representing the distribution of the non-noisy model remains unchanged and identifiable. Moreover, the distribution of the minimum ratio estimators of the model parameters at their left limits are completely determined by the distribution of the noise variables up to a positive constant. Under a regular variation condition on the noise variables we prove that the estimated Kleene star matrix converges to a matrix of independent Weibull entries after proper centering and scaling.


title: Order-isomorphic twins in permutations

abstract: Let $a_1,\dotsc,a_n$ be a permutation of $[n]$. Two disjoint order-isomorphic subsequences are called \emph{twins}. We show that every permutation of $[n]$ contains twins of length $\Omega(n^{3/5})$ improving the trivial bound of $\Omega(n^{1/2})$. We also show that a random permutation contains twins of length $\Omega(n^{2/3})$, which is sharp.


title: Undergraduate Researcher Personas

abstract: Personalized undergraduate research programs can help increase undergraduate students' participation in research. Personas of undergraduate researchers are a powerful means of encapsulating the richness of students' various goals, motivations, and experiences in research. Student personas can support the design of research programs with student-centered approaches that better engage undergraduate students across multiple disciplines. We used interview data and drew on Self-Determination Theory to identify students' various goals and motivations in order to build three undergraduate researcher personas in physics.


title: Budget-Constrained Bandits over General Cost and Reward Distributions

abstract: We consider a budget-constrained bandit problem where each arm pull incurs a random cost, and yields a random reward in return. The objective is to maximize the total expected reward under a budget constraint on the total cost. The model is general in the sense that it allows correlated and potentially heavy-tailed cost-reward pairs that can take on negative values as required by many applications. We show that if moments of order $(2+\gamma)$ for some $\gamma > 0$ exist for all cost-reward pairs, $O(\log B)$ regret is achievable for a budget $B>0$. In order to achieve tight regret bounds, we propose algorithms that exploit the correlation between the cost and reward of each arm by extracting the common information via linear minimum mean-square error estimation. We prove a regret lower bound for this problem, and show that the proposed algorithms achieve tight problem-dependent regret bounds, which are optimal up to a universal constant factor in the case of jointly Gaussian cost and reward pairs.


title: New rational cubic fourfolds arising from Cremona transformations

abstract: Are Fourier-Mukai equivalent cubic fourfolds birationally equivalent? We obtain an affirmative answer to this question for very general cubic fourfolds of discriminant 20, where we produce birational maps via the Cremona transformation defined by the Veronese surface. By studying how these maps act on the cubics known to be rational, we surprisingly found new rational examples.


title: Imaging Andreev Reflection in Graphene

abstract: Coherent charge transport along ballistic paths can be introduced into graphene by Andreev reflection, for which an electron reflects from a superconducting contact as a hole, while a Cooper pair is transmitted. We use a liquid-helium cooled scanning gate microscope (SGM) to image Andreev reflection in graphene in the magnetic focusing regime, where carriers move along cyclotron orbits between contacts. Images of flow are obtained by deflecting carrier paths and displaying the resulting change in conductance. When electrons enter the the superconductor, Andreev-reflected holes leave for the collecting contact. To test the results, we destroy Andreev reflection with a large current and by heating above the critical temperature. In both cases, the reflected carriers change from holes to electrons.


title: The Marathon 2: A Navigation System

abstract: Developments in mobile robot navigation have enabled robots to operate in warehouses, retail stores, and on sidewalks around pedestrians. Various navigation solutions have been proposed, though few as widely adopted as ROS Navigation. 10 years on, it is still one of the most popular navigation solutions. Yet, ROS Navigation has failed to keep up with modern trends. We propose the new navigation solution, Navigation2, which builds on the successful legacy of ROS Navigation. Navigation2 uses a behavior tree for navigator task orchestration and employs new methods designed for dynamic environments applicable to a wider variety of modern sensors. It is built on top of ROS2, a secure message passing framework suitable for safety critical applications and program lifecycle management. We present experiments in a campus setting utilizing Navigation2 to operate safely alongside students over a marathon as an extension of the experiment proposed in Eppstein et al. The Navigation2 system is freely available at https://github.com/ros-planning/navigation2 with a rich community and instructions.


title: MIndGrasp: A New Training and Testing Framework for Motor Imagery Based 3-Dimensional Assistive Robotic Control

abstract: With increasing global age and disability assistive robots are becoming more necessary, and brain computer interfaces (BCI) are often proposed as a solution to understanding the intent of a disabled person that needs assistance. Most frameworks for electroencephalography (EEG)-based motor imagery (MI) BCI control rely on the direct control of the robot in Cartesian space. However, for 3-dimensional movement, this requires 6 motor imagery classes, which is a difficult distinction even for more experienced BCI users. In this paper, we present a simulated training and testing framework which reduces the number of motor imagery classes to 4 while still grasping objects in three-dimensional space. This is achieved through semi-autonomous eye-in-hand vision-based control of the robotic arm, while the user-controlled BCI achieves movement to the left and right, as well as movement toward and away from the object of interest. Additionally, the framework includes a method of training a BCI directly on the assistive robotic system, which should be more easily transferrable to a real-world assistive robot than using a standard training protocol such as Graz-BCI. Presented results do not consider real human EEG data, but are rather shown as a baseline for comparison with future human data and other improvements on the system.


title: PlaNet of the Bayesians: Reconsidering and Improving Deep Planning Network by Incorporating Bayesian Inference

abstract: In the present paper, we propose an extension of the Deep Planning Network (PlaNet), also referred to as PlaNet of the Bayesians (PlaNet-Bayes). There has been a growing demand in model predictive control (MPC) in partially observable environments in which complete information is unavailable because of, for example, lack of expensive sensors. PlaNet is a promising solution to realize such latent MPC, as it is used to train state-space models via model-based reinforcement learning (MBRL) and to conduct planning in the latent space. However, recent state-of-the-art strategies mentioned in MBRR literature, such as involving uncertainty into training and planning, have not been considered, significantly suppressing the training performance. The proposed extension is to make PlaNet uncertainty-aware on the basis of Bayesian inference, in which both model and action uncertainty are incorporated. Uncertainty in latent models is represented using a neural network ensemble to approximately infer model posteriors. The ensemble of optimal action candidates is also employed to capture multimodal uncertainty in the optimality. The concept of the action ensemble relies on a general variational inference MPC (VI-MPC) framework and its instance, probabilistic action ensemble with trajectory sampling (PaETS). In this paper, we extend VI-MPC and PaETS, which have been originally introduced in previous literature, to address partially observable cases. We experimentally compare the performances on continuous control tasks, and conclude that our method can consistently improve the asymptotic performance compared with PlaNet.


title: Estimating Multiple Precision Matrices with Cluster Fusion Regularization

abstract: We propose a penalized likelihood framework for estimating multiple precision matrices from different classes. Most existing methods either incorporate no information on relationships between the precision matrices, or require this information be known a priori. The framework proposed in this article allows for simultaneous estimation of the precision matrices and relationships between the precision matrices, jointly. Sparse and non-sparse estimators are proposed, both of which require solving a non-convex optimization problem. To compute our proposed estimators, we use an iterative algorithm which alternates between a convex optimization problem solved by blockwise coordinate descent and a k-means clustering problem. Blockwise updates for computing the sparse estimator require solving an elastic net penalized precision matrix estimation problem, which we solve using a proximal gradient descent algorithm. We prove that this subalgorithm has a linear rate of convergence. In simulation studies and two real data applications, we show that our method can outperform competitors that ignore relevant relationships between precision matrices and performs similarly to methods which use prior information often uknown in practice.


title: A Globally Convergent Newton Method for Polynomials

abstract: Newton's method for polynomial root finding is one of mathematics' most well-known algorithms. The method also has its shortcomings: it is undefined at critical points, it could exhibit chaotic behavior and is only guaranteed to converge locally. Based on the {\it Geometric Modulus Principle} for a complex polynomial $p(z)$, together with a {\it Modulus Reduction Theorem} proved here, we develop the {\it Robust Newton's method} (RNM), defined everywhere with a step-size that guarantees an {\it a priori} reduction in polynomial modulus in each iteration. Furthermore, we prove RNM iterates converge globally, either to a root or a critical point. Specifically, given $\varepsilon $ and any seed $z_0$, in $t=O(1/\varepsilon^{2})$ iterations of RNM, independent of degree of $p(z)$, either $|p(z_t)| \leq \varepsilon$ or $|p(z_t) p'(z_t)| \leq \varepsilon$. By adjusting the iterates at {\it near-critical points}, we describe a {\it modified} RNM that necessarily convergence to a root. In combination with Smale's point estimation, RNM results in a globally convergent Newton's method having a locally quadratic rate. We present sample polynomiographs that demonstrate how in contrast with Newton's method RNM smooths out the fractal boundaries of basins of attraction of roots. RNM also finds potentials in computing all roots of arbitrary degree polynomials. A particular consequence of RNM is a simple algorithm for solving cubic equations.


title: Wilsonian Renormalization Group In The Functional Nonperturbative Approach

abstract: We consider a functional relation between a given Wilsonian RG flow, which has to be related to a specific coarse-graining procedure, and an infinite family of (UV cutoff) scale dependent field redefinitions. Within this framework one can define a family of Wilsonian proper-time exact RG equations associated to an arbitrary regulator function. New applications of these RG flow schemes to the Ising Universality class in three dimensions in the derivative expansion are shortly illustrated.


title: Fundamental bounds for scattering from absorptionless electromagnetic structures

abstract: The ability to design the scattering properties of electromagnetic structures is of fundamental interest in optical science and engineering. While there has been great practical success applying local optimization methods to electromagnetic device design, it is unclear whether the performance of resulting designs is close to that of the best possible design. This question remains unsettled for absorptionless electromagnetic devices since the absence of material loss makes it difficult to provide provable bounds on their scattering properties. We resolve this problem by providing non-trivial lower bounds on performance metrics that are convex functions of the scattered fields. Our bounding procedure relies on accounting for a constraint on the electric fields inside the device, which can be provably constructed for devices with small footprints or low dielectric constrast. We illustrate our bounding procedure by studying limits on the scattering cross-sections of dielectric and metallic particles in the absence of material losses.


title: Birth environment of circumbinary planets: are there CBPs on the inclined orbits?

abstract: The inclination distribution of circumbinary planets (CBPs) is an important scientific issue. It is of great significance in estimating the occurrence rate of CBPs and studying their formation and evolution. Although the CBPs currently discovered by the transit method are nearly coplanar, the true inclination distribution of CBPs is still unknown. Previous researches on CBPs mostly regarded them as an isolated binary-planet system, without considering the birth environment of their host binaries. It is generally believed that almost all stars are born in clusters. Therefore, it is necessary to consider the impact of the close encounters of stars on the CBP systems. This article discusses how the close encounters of fly-by stars affect the inclination of CBP. Based on extensive numerical simulations, we found that CBPs in close binary with a spacing of about 0.2 au are almost unaffected by star fly-bys. Their orbits remain coplanar. However, when the spacing of the binary stars is greater than 1 au, the 2-3 fly-bys of the intruding star can excite a considerable inclination even for the CBP near the unstable boundary of the binary. For the planets in the outer region, a single star fly-by can excite an inclination to more than 5 degrees. Especially, CBPs in near polar or retrograde orbits can naturally form through bianry-star encounters. If close binaries are born in open clusters, our simulations suggest that there may be high-inclination CBPs in binaries with spacing > 1 au.


title: AccQOC: Accelerating Quantum Optimal Control Based Pulse Generation

abstract: In the last decades, we have witnessed the rapid growth of Quantum Computing. In the current Noisy Intermediate-Scale Quantum (NISQ) era, the capability of a quantum machine is limited by the decoherence time, gate fidelity and the number of Qubits. Current quantum computing applications are far from the real "quantum supremacy" due to the fragile physical Qubits, which can only be entangled for a few microseconds. Recent works use quantum optimal control to reduce the latency of quantum circuits, thereby effectively increasing quantum volume. However, the key challenge of this technique is the large overhead due to long compilation time. In this paper, we propose AccQOC, a comprehensive static/dynamic hybrid workflow to transform gate groups (equivalent to matrices) to pulses using QOC (Quantum Optimal Control) with a reasonable compilation time budget. AccQOC is composed of static pre-compilation and accelerated dynamic compilation. With the methodology of AccQOC, we reached a balanced point of compilation time and overall latency. The results show that accelerated compilation based on MST achieves 9.88x compilation speedup compared to the standard compilation of each group while maintaining an average 2.43x latency reduction compared with gate-based compilation.


title: Most Probable Dynamics of Stochastic Dynamical Systems with Exponentially Light Jump Fluctuations

abstract: The emergence of the exit events from a bounded domain containing a stable fixed point induced by non-Gaussian L\'evy fluctuations plays a pivotal role in practical physical systems. In the limit of weak noise, we develop a Hamiltonian formalism under the L\'evy fluctuations with exponentially light jumps for one- and two-dimensional stochastic dynamical systems. This formalism is based on a recently proved large deviation principle for dynamical systems under non-Gaussian L\'evy perturbations. We demonstrate how to compute the most probable exit path and the quasi-potential by several examples. Meanwhile, we explore the impacts of the jump measure on the quasi-potential quantitatively and on the most probable exit path qualitatively. Results show that the quasi-potential can be well estimated by an approximate analytical expression. Moreover, we discover that although the most probable exit paths are analogous to the Gaussian case for the isotropic noise, the anisotropic noise leads to significant changes of the structure of the exit paths. These findings shed light on the underlying qualitative mechanism and quantitative feature of the exit phenomenon induced by non-Gaussian noise.


title: Understanding the Intrinsic Robustness of Image Distributions using Conditional Generative Models

abstract: Starting with Gilmer et al. (2018), several works have demonstrated the inevitability of adversarial examples based on different assumptions about the underlying input probability space. It remains unclear, however, whether these results apply to natural image distributions. In this work, we assume the underlying data distribution is captured by some conditional generative model, and prove intrinsic robustness bounds for a general class of classifiers, which solves an open problem in Fawzi et al. (2018). Building upon the state-of-the-art conditional generative models, we study the intrinsic robustness of two common image benchmarks under $\ell_2$ perturbations, and show the existence of a large gap between the robustness limits implied by our theory and the adversarial robustness achieved by current state-of-the-art robust models. Code for all our experiments is available at https://github.com/xiaozhanguva/Intrinsic-Rob.


title: Inverse design of photonic crystals through automatic differentiation

abstract: Gradient-based inverse design in photonics has already achieved remarkable results in designing small-footprint, high-performance optical devices. The adjoint variable method, which allows for the efficient computation of gradients, has played a major role in this success. However, gradient-based optimization has not yet been applied to the mode-expansion methods that are the most common approach to studying periodic optical structures like photonic crystals. This is because, in such simulations, the adjoint variable method cannot be defined as explicitly as in standard finite-difference or finite-element time- or frequency-domain methods. Here, we overcome this through the use of automatic differentiation, which is a generalization of the adjoint variable method to arbitrary computational graphs. We implement the plane-wave expansion and the guided-mode expansion methods using an automatic differentiation library, and show that the gradient of any simulation output can be computed efficiently and in parallel with respect to all input parameters. We then use this implementation to optimize the dispersion of a photonic crystal waveguide, and the quality factor of an ultra-small cavity in a lithium niobate slab. This extends photonic inverse design to a whole new class of simulations, and more broadly highlights the importance that automatic differentiation could play in the future for tracking and optimizing complicated physical models.

categories: physics.optics cs.CE physics.app-ph physics.comp-ph


title: Unblind Your Apps: Predicting Natural-Language Labels for Mobile GUI Components by Deep Learning

abstract: According to the World Health Organization(WHO), it is estimated that approximately 1.3 billion people live with some forms of vision impairment globally, of whom 36 million are blind. Due to their disability, engaging these minority into the society is a challenging problem. The recent rise of smart mobile phones provides a new solution by enabling blind users' convenient access to the information and service for understanding the world. Users with vision impairment can adopt the screen reader embedded in the mobile operating systems to read the content of each screen within the app, and use gestures to interact with the phone. However, the prerequisite of using screen readers is that developers have to add natural-language labels to the image-based components when they are developing the app. Unfortunately, more than 77% apps have issues of missing labels, according to our analysis of 10,408 Android apps. Most of these issues are caused by developers' lack of awareness and knowledge in considering the minority. And even if developers want to add the labels to UI components, they may not come up with concise and clear description as most of them are of no visual issues. To overcome these challenges, we develop a deep-learning based model, called LabelDroid, to automatically predict the labels of image-based buttons by learning from large-scale commercial apps in Google Play. The experimental results show that our model can make accurate predictions and the generated labels are of higher quality than that from real Android developers.


title: Statistical power for cluster analysis

abstract: Cluster algorithms are increasingly popular in biomedical research due to their compelling ability to identify discrete subgroups in data, and their increasing accessibility in mainstream software. While guidelines exist for algorithm selection and outcome evaluation, there are no firmly established ways of computing a priori statistical power for cluster analysis. Here, we estimated power and accuracy for common analysis pipelines through simulation. We varied subgroup size, number, separation (effect size), and covariance structure. We then subjected generated datasets to dimensionality reduction (none, multidimensional scaling, or UMAP) and cluster algorithms (k-means, agglomerative hierarchical clustering with Ward or average linkage and Euclidean or cosine distance, HDBSCAN). Finally, we compared the statistical power of discrete (k-means), "fuzzy" (c-means), and finite mixture modelling approaches (which include latent profile and latent class analysis). We found that outcomes were driven by large effect sizes or the accumulation of many smaller effects across features, and were unaffected by differences in covariance structure. Sufficient statistical power was achieved with relatively small samples (N=20 per subgroup), provided cluster separation is large ({\Delta}=4). Fuzzy clustering provided a more parsimonious and powerful alternative for identifying separable multivariate normal distributions, particularly those with slightly lower centroid separation ({\Delta}=3). Overall, we recommend that researchers 1) only apply cluster analysis when large subgroup separation is expected, 2) aim for sample sizes of N=20 to N=30 per expected subgroup, 3) use multidimensional scaling to improve cluster separation, and 4) use fuzzy clustering or finite mixture modelling approaches that are more powerful and more parsimonious with partially overlapping multivariate normal distributions.


title: Analogue Gravity on a Superconducting Chip

abstract: We describe how analogues of a Hawking evaporating black hole as well as the Unruh effect for an oscillatory, accelerating photodetector in vacuum may be realized using superconducting, microwave circuits that are fashioned out of Josephson tunnel junction and film bulk acoustic resonator elements.


title: AoI and Energy Consumption Oriented Dynamic Status Updating in Caching Enabled IoT Networks

abstract: Caching has been regarded as a promising technique to alleviate energy consumption of sensors in Internet of Things (IoT) networks by responding to users' requests with the data packets stored in the edge caching node (ECN). For real-time applications in caching enabled IoT networks, it is essential to develop dynamic status update strategies to strike a balance between the information freshness experienced by users and energy consumed by the sensor, which, however, is not well addressed. In this paper, we first depict the evolution of information freshness, in terms of age of information (AoI), at each user. Then, we formulate a dynamic status update optimization problem to minimize the expectation of a long term accumulative cost, which jointly considers the users' AoI and sensor's energy consumption. To solve this problem, a Markov Decision Process (MDP) is formulated to cast the status updating procedure, and a model-free reinforcement learning algorithm is proposed, with which the challenge brought by the unknown of the formulated MDP's dynamics can be addressed. Finally, simulations are conducted to validate the convergence of our proposed algorithm and its effectiveness compared with the zero-wait baseline policy.


title: Average Age of Changed Information in the Internet of Things

abstract: The freshness of status updates is imperative in mission-critical Internet of things (IoT) applications. Recently, Age of Information (AoI) has been proposed to measure the freshness of updates at the receiver. However, AoI only characterizes the freshness over time, but ignores the freshness in the content. In this paper, we introduce a new performance metric, Age of Changed Information (AoCI), which captures both the passage of time and the change of information content. Also, we examine the AoCI in a time-slotted status update system, where a sensor samples the physical process and transmits the update packets with a cost. We formulate a Markov Decision Process (MDP) to find the optimal updating policy that minimizes the weighted sum of the AoCI and the update cost. Particularly, in a special case that the physical process is modeled by a two-state discrete time Markov chain with equal transition probability, we show that the optimal policy is of threshold type with respect to the AoCI and derive the closed-form of the threshold. Finally, simulations are conducted to exhibit the performance of the threshold policy and its superiority over the zero-wait baseline policy.


title: Vision-based Robot Manipulation Learning via Human Demonstrations

abstract: Vision-based learning methods provide promise for robots to learn complex manipulation tasks. However, how to generalize the learned manipulation skills to real-world interactions remains an open question. In this work, we study robotic manipulation skill learning from a single third-person view demonstration by using activity recognition and object detection in computer vision. To facilitate generalization across objects and environments, we propose to use a prior knowledge base in the form of a text corpus to infer the object to be interacted with in the context of a robot. We evaluate our approach in a real-world robot, using several simple and complex manipulation tasks commonly performed in daily life. The experimental results show that our approach achieves good generalization performance even from small amounts of training data.


title: Experimental Evaluation of a Pseudo-Doppler Direction-Finding System for Localizing Radio Tags

abstract: We present the design of a radio antenna system for obtaining instantaneous bearing measurements towards a radio emitter. Our work is motivated by applications where robots are used for localizing and tracking radio-tagged wildlife. The traditional method is to use directional antennas that need to be rotated in order find the bearing which is time consuming. Instead, we present a low-cost system capable of finding bearing measurements almost instantaneously using an antenna array. This is particularly appealing for wildlife tracking with Unmanned Aerial Systems (UASs) where remaining stationary can be challenging and energy consuming, in addition to being slow. The proposed system uses existing open source hardware and software systems and leverages principles of pseudo Doppler direction-finding. The resulting system was tested in an anechoic chamber and in outdoor settings. The outdoor tests with particle filtering show that the resulting system is capable of localizing radio tags within 5 meter accuracy starting with an initial estimate of 200m x 200m.


title: Say As You Wish: Fine-grained Control of Image Caption Generation with Abstract Scene Graphs

abstract: Humans are able to describe image contents with coarse to fine details as they wish. However, most image captioning models are intention-agnostic which can not generate diverse descriptions according to different user intentions initiatively. In this work, we propose the Abstract Scene Graph (ASG) structure to represent user intention in fine-grained level and control what and how detailed the generated description should be. The ASG is a directed graph consisting of three types of \textbf{abstract nodes} (object, attribute, relationship) grounded in the image without any concrete semantic labels. Thus it is easy to obtain either manually or automatically. From the ASG, we propose a novel ASG2Caption model, which is able to recognise user intentions and semantics in the graph, and therefore generate desired captions according to the graph structure. Our model achieves better controllability conditioning on ASGs than carefully designed baselines on both VisualGenome and MSCOCO datasets. It also significantly improves the caption diversity via automatically sampling diverse ASGs as control signals.


title: The symplectic group in Polymer Quantum Mechanics

abstract: In this paper we provide the representation of the symplectic group $Sp(2n, \mathbb{R})$ in polymer quantum mechanics. We derive the propagator of the polymer free particle and the polymer harmonic oscillator without considering a polymer scale. The polymer scale is then introduced to reconcile our results with those expressions for the polymer free particle. The propagator for the polymer harmonic oscillator implies non-unitary evolution.


title: Joint Wasserstein Distribution Matching

abstract: Joint distribution matching (JDM) problem, which aims to learn bidirectional mappings to match joint distributions of two domains, occurs in many machine learning and computer vision applications. This problem, however, is very difficult due to two critical challenges: (i) it is often difficult to exploit sufficient information from the joint distribution to conduct the matching; (ii) this problem is hard to formulate and optimize. In this paper, relying on optimal transport theory, we propose to address JDM problem by minimizing the Wasserstein distance of the joint distributions in two domains. However, the resultant optimization problem is still intractable. We then propose an important theorem to reduce the intractable problem into a simple optimization problem, and develop a novel method (called Joint Wasserstein Distribution Matching (JWDM)) to solve it. In the experiments, we apply our method to unsupervised image translation and cross-domain video synthesis. Both qualitative and quantitative comparisons demonstrate the superior performance of our method over several state-of-the-arts.


title: Detective quantum efficiency of photon-counting CdTe and Si detectors for computed tomography: a simulation study

abstract: Purpose: Developing photon-counting CT detectors requires understanding the impact of parameters such as converter material, absorption length and pixel size. We apply a novel linear-systems framework, incorporating spatial and energy resolution, to study realistic silicon (Si) and cadmium telluride (CdTe) detectors at low count rate. Approach: We compared CdTe detector designs with $0.5\times0.5\; \mathrm{mm}^2$ and $0.225\times0.225\; \mathrm{mm}^2$ pixels and Si detector designs with $0.5\times0.5\; \mathrm{mm}^2$ pixels of 30 and 60 mm active absorption length, with and without tungsten scatter blockers. Monte-Carlo simulations of photon transport were used together with Gaussian charge sharing models fitted to published data. Results: For detection in a 300 mm thick object at 120 kVp, the 0.5 mm and 0.225 mm pixel CdTe systems have 28-41 $\%$ and 5-29 $\%$ higher DQE, respectively, than the 60 mm Si system with tungsten, whereas the corresponding numbers for two-material decomposition are 2 $\%$ lower to 11 $\%$ higher DQE and 31-54 $\%$ lower DQE compared to Si. We also show that combining these detectors with dual-spectrum acquisition is beneficial. Conclusions: In the low-count-rate regime, CdTe detector systems outperform the Si systems for detection tasks, while silicon outperforms one or both of the CdTe systems for material decomposition.


title: Deep Reinforcement Learning for Fresh Data Collection in UAV-assisted IoT Networks

abstract: Due to the flexibility and low operational cost, dispatching unmanned aerial vehicles (UAVs) to collect information from distributed sensors is expected to be a promising solution in Internet of Things (IoT), especially for time-critical applications. How to maintain the information freshness is a challenging issue. In this paper, we investigate the fresh data collection problem in UAV-assisted IoT networks. Particularly, the UAV flies towards the sensors to collect status update packets within a given duration while maintaining a non-negative residual energy. We formulate a Markov Decision Process (MDP) to find the optimal flight trajectory of the UAV and transmission scheduling of the sensors that minimizes the weighted sum of the age of information (AoI). A UAV-assisted data collection algorithm based on deep reinforcement learning (DRL) is further proposed to overcome the curse of dimensionality. Extensive simulation results demonstrate that the proposed DRL-based algorithm can significantly reduce the weighted sum of the AoI compared to other baseline algorithms.


title: Fine-grained Video-Text Retrieval with Hierarchical Graph Reasoning

abstract: Cross-modal retrieval between videos and texts has attracted growing attentions due to the rapid emergence of videos on the web. The current dominant approach for this problem is to learn a joint embedding space to measure cross-modal similarities. However, simple joint embeddings are insufficient to represent complicated visual and textual details, such as scenes, objects, actions and their compositions. To improve fine-grained video-text retrieval, we propose a Hierarchical Graph Reasoning (HGR) model, which decomposes video-text matching into global-to-local levels. To be specific, the model disentangles texts into hierarchical semantic graph including three levels of events, actions, entities and relationships across levels. Attention-based graph reasoning is utilized to generate hierarchical textual embeddings, which can guide the learning of diverse and hierarchical video representations. The HGR model aggregates matchings from different video-text levels to capture both global and local details. Experimental results on three video-text datasets demonstrate the advantages of our model. Such hierarchical decomposition also enables better generalization across datasets and improves the ability to distinguish fine-grained semantic differences.


title: Deep Active Learning for Biased Datasets via Fisher Kernel Self-Supervision

abstract: Active learning (AL) aims to minimize labeling efforts for data-demanding deep neural networks (DNNs) by selecting the most representative data points for annotation. However, currently used methods are ill-equipped to deal with biased data. The main motivation of this paper is to consider a realistic setting for pool-based semi-supervised AL, where the unlabeled collection of train data is biased. We theoretically derive an optimal acquisition function for AL in this setting. It can be formulated as distribution shift minimization between unlabeled train data and weakly-labeled validation dataset. To implement such acquisition function, we propose a low-complexity method for feature density matching using self-supervised Fisher kernel (FK) as well as several novel pseudo-label estimators. Our FK-based method outperforms state-of-the-art methods on MNIST, SVHN, and ImageNet classification while requiring only 1/10th of processing. The conducted experiments show at least 40% drop in labeling efforts for the biased class-imbalanced data compared to existing methods.


title: Stable behaviour of infinitely wide deep neural networks

abstract: We consider fully connected feed-forward deep neural networks (NNs) where weights and biases are independent and identically distributed as symmetric centered stable distributions. Then, we show that the infinite wide limit of the NN, under suitable scaling on the weights, is a stochastic process whose finite-dimensional distributions are multivariate stable distributions. The limiting process is referred to as the stable process, and it generalizes the class of Gaussian processes recently obtained as infinite wide limits of NNs (Matthews at al., 2018b). Parameters of the stable process can be computed via an explicit recursion over the layers of the network. Our result contributes to the theory of fully connected feed-forward deep NNs, and it paves the way to expand recent lines of research that rely on Gaussian infinite wide limits.


title: Authentication, Access Control, Privacy, Threats and Trust Management Towards Securing Fog Computing Environments: A Review

abstract: Fog computing is an emerging computing paradigm that has come into consideration for the deployment of IoT applications amongst researchers and technology industries over the last few years. Fog is highly distributed and consists of a wide number of autonomous end devices, which contribute to the processing. However, the variety of devices offered across different users are not audited. Hence, the security of Fog devices is a major concern in the Fog computing environment. Furthermore, mitigating and preventing those security measures is a research issue. Therefore, to provide the necessary security for Fog devices, we need to understand what the security concerns are with regards to Fog. All aspects of Fog security, which have not been covered by other literature works needs to be identified and need to be aggregate all issues in Fog security. It needs to be noted that computation devices consist of many ordinary users, and are not managed by any central entity or managing body. Therefore, trust and privacy is also a key challenge to gain market adoption for Fog. To provide the required trust and privacy, we need to also focus on authentication, threats and access control mechanisms as well as techniques in Fog computing. In this paper, we perform a survey and propose a taxonomy, which presents an overview of existing security concerns in the context of the Fog computing paradigm. We discuss the Blockchain-based solutions towards a secure Fog computing environment and presented various research challenges and directions for future research.


title: Diameter two properties and the Radon-Nikod\'ym property in Orlicz spaces

abstract: Some necessary and sufficient conditions are found for Banach function lattices to have the Radon-Nikod\'ym property. Consequently it is shown that an Orlicz space $L_\varphi$ over a non-atomic $\sigma$-finite measure space $(\Omega, \Sigma,\mu)$, not necessarily separable, has the Radon-Nikod\'ym property if and only if $\varphi$ is an $N$-function at infinity and satisfies the appropriate $\Delta_2$ condition. For an Orlicz sequence space $\ell_\varphi$, it has the Radon-Nikod\'ym property if and only if $\varphi$ satisfies condition $\Delta_2^0$. In the second part the relationships between uniformly $\ell_1^2$ points of the unit sphere of a Banach space and the diameter of the slices are studied. Using these results, a quick proof is given that an Orlicz space $L_\varphi$ has the Daugavet property only if $\varphi$ is linear, so when $L_\varphi$ is isometric to $L_1$. The other consequence is that the Orlicz spaces equipped with the Orlicz norm generated by $N$-functions never have local diameter two property, while it is well-known that when equipped with the Luxemburg norm, it may have that property. Finally, it is shown that the local diameter two property, the diameter two property, the strong diameter two property are equivalent in function and sequence Orlicz spaces with the Luxemburg norm under appropriate conditions on $\varphi$.


title: Intelligent Home 3D: Automatic 3D-House Design from Linguistic Descriptions Only

abstract: Home design is a complex task that normally requires architects to finish with their professional skills and tools. It will be fascinating that if one can produce a house plan intuitively without knowing much knowledge about home design and experience of using complex designing tools, for example, via natural language. In this paper, we formulate it as a language conditioned visual content generation problem that is further divided into a floor plan generation and an interior texture (such as floor and wall) synthesis task. The only control signal of the generation process is the linguistic expression given by users that describe the house details. To this end, we propose a House Plan Generative Model (HPGM) that first translates the language input to a structural graph representation and then predicts the layout of rooms with a Graph Conditioned Layout Prediction Network (GC LPN) and generates the interior texture with a Language Conditioned Texture GAN (LCT-GAN). With some post-processing, the final product of this task is a 3D house model. To train and evaluate our model, we build the first Text-to-3D House Model dataset.


title: On degenerate gamma functions

abstract: Recently, the degenerate gamma functions are introduced as a degenerate version of the usual gamma function by Kim-Kim. In this paper, we investigate several properties of them. Namely, we obtain an analytic continuation as a meromorphic function on the whole complex plane,the difference formula, the values at positive integers, some expressions following from the Weierstrass and Euler formulas for the usual gamma function and an integral representation as the integral along a Hankel contour.


title: The cross cyclomatic complexity: a bi-dimensional measure for program complexity on graphs

abstract: Reduce and control complexity is an essential practice in software design. Cyclomatic complexity (CC) is one of the most popular software metrics, applied for more than 40 years. Despite CC is an interesting metric to highlight the number of branches in a program, it clearly not sufficient to represent the complexity in a piece of software. In this paper, we introduce the cross cyclomatic complexity (CCC), a new bi-dimensional complexity measure on graphs that combines the cyclomatic complexity and the weight of a minimum-weight cycle basis in as pair on the Cartesian plan to characterize program complexity using control flow graphs. Our postulates open a new venue to represent program complexity, and we discuss its implications and opportunities.


title: Learn Task First or Learn Human Partner First: A Hierarchical Task Decomposition Method for Human-Robot Cooperation

abstract: Applying Deep Reinforcement Learning (DRL) to Human-Robot Cooperation (HRC) in dynamic control problems is promising yet challenging as the robot needs to learn the dynamics of the controlled system and dynamics of the human partner. In existing research, the robot powered by DRL adopts coupled observation of the environment and the human partner to learn both dynamics simultaneously. However, such a learning strategy is limited in terms of learning efficiency and team performance. This work proposes a novel task decomposition method with a hierarchical reward mechanism that enables the robot to learn the hierarchical dynamic control task separately from learning the human partner's behavior. The method is validated with a hierarchical control task in a simulated environment with human subject experiments. Our method also provides insight into the design of the learning strategy for HRC. The results show that the robot should learn the task first to achieve higher team performance and learn the human first to achieve higher learning efficiency.


title: A flexible Bayesian framework to estimate age- and cause-specific child mortality over time from sample registration data

abstract: In order to implement disease-specific interventions in young age groups, policy makers in low- and middle-income countries require timely and accurate estimates of age- and cause-specific child mortality. High quality data is not available in settings where these interventions are most needed, but there is a push to create sample registration systems that collect detailed mortality information. Current methods that estimate mortality from this data employ multistage frameworks without rigorous statistical justification that separately estimate all-cause and cause-specific mortality and are not sufficiently adaptable to capture important features of the data. We propose a flexible Bayesian modeling framework to estimate age- and cause-specific child mortality from sample registration data. We provide a theoretical justification for the framework, explore its properties via simulation, and use it to estimate mortality trends using data from the Maternal and Child Health Surveillance System in China.