title: Elastic and Fracture Behavior of Three-Dimensional Ply-to-Ply Angle Interlock Woven Composites: Through-Thickness, Size Effect, and Multiaxial Tests

abstract: This work presents a comprehensive investigation of the elastic and fracture behavior of ply-to-ply angle interlock three-dimensional woven composites. The research investigated novel splitting and wedge-driven out-of-plane fracture tests to shed light on the tensile fracture behavior in the thickness direction and to provide estimates of the out-of-plane tensile strength and fracture energy. In addition, size effect tests on geometrically-scaled Single Edge Notch Tension (SENT) specimens were performed to fully characterize the intra-laminar fracture energy of the material and to study the scaling of structural strength in this type of three-dimensional composites. The results confirmed that size effect in the structural strength of these materials is significant. In fact, even if the range of sizes investigated was broader than in any previous size effect study on traditional laminated composites and two-dimensional textile composites, all the experimental data fell in the transition zone between quasi-ductile and brittle behavior. This implies strong damage tolerance of the investigated three-dimensional composites. The analysis of the data via Bazant's Type II Size Effect Law (SEL) enabled the objective characterization of the intra-laminar fracture energy of three-dimensional composites for the first time. Finally, Arcan rig tests combined with X-ray micro-computed tomography allowed unprecedented insights on the different damage mechanisms under multi-axial nominal loading conditions, particularly tension-dominated and shear-dominated conditions.


title: Polarization Calibration of the BICEP3 CMB polarimeter at the South Pole

abstract: The BICEP3 CMB Polarimeter is a small-aperture refracting telescope located at the South Pole and is specifically designed to search for the possible signature of inflationary gravitational waves in the Cosmic Microwave Background (CMB). The experiment measures polarization on the sky by differencing the signal of co-located, orthogonally polarized antennas coupled to Transition Edge Sensor (TES) detectors. We present precise measurements of the absolute polarization response angles and polarization efficiencies for nearly all of BICEP3s $\sim800$ functioning polarization-sensitive detector pairs from calibration data taken in January 2018. Using a Rotating Polarized Source (RPS), we mapped polarization response for each detector over a full 360 degrees of source rotation and at multiple telescope boresight rotations from which per-pair polarization properties were estimated. In future work, these results will be used to constrain signals predicted by exotic physical models such as Cosmic Birefringence.


title: The Fundamental Plane of Massive Quiescent Galaxies at z~2

abstract: We examine the Fundamental Plane (FP) and mass-to-light ratio ($M/L$) scaling relations using the largest sample of massive quiescent galaxies at $1.5<z<2.5$ to date. The FP ($r_{e}, \sigma_{e}, I_{e}$) is established using $19$ $UVJ$ quiescent galaxies from COSMOS with $Hubble$ $Space$ $Telescope$ $(HST)$ $H_{F160W}$ rest-frame optical sizes and X-shooter absorption line measured stellar velocity dispersions. For a very massive, ${\rm{log}}(M_{\ast}/M_{\odot})>11.26$, subset of 8 quiescent galaxies at $z>2$, from Stockmann et al. (2020), we show that they cannot passively evolve to the local Coma cluster relation alone and must undergo significant structural evolution to mimic the sizes of local massive galaxies. The evolution of the FP and $M/L$ scaling relations, from $z=2$ to present-day, for this subset are consistent with passive aging of the stellar population and minor merger structural evolution into the most massive galaxies in the Coma cluster and other massive elliptical galaxies from the MASSIVE Survey. Modeling the luminosity evolution from minor merger added stellar populations favors a history of merging with "dry" quiescent galaxies.


title: Generalized fiducial factor: an alternative to the Bayes factor for forensic identification of source problems

abstract: One formulation of forensic identification of source problems is to determine the source of trace evidence, for instance, glass fragments found on a suspect for a crime. The current state of the science is to compute a Bayes factor (BF) comparing the marginal distribution of measurements of trace evidence under two competing propositions for whether or not the unknown source evidence originated from a specific source. The obvious problem with such an approach is the ability to tailor the prior distributions (placed on the features/parameters of the statistical model for the measurements of trace evidence) in favor of the defense or prosecution, which is further complicated by the fact that the typical number of measurements of trace evidence is typically sufficiently small that prior choice/specification has a strong influence on the value of the BF. To remedy this problem of prior specification and choice, we develop an alternative to the BF, within the framework of generalized fiducial inference (GFI), that we term a {\em generalized fiducial factor} (GFF). Furthermore, we demonstrate empirically, on the synthetic and real Netherlands Forensic Institute (NFI) casework data, deficiencies in the BF and classical/frequentist likelihood ratio (LR) approaches.


title: A Recursive Method for Real-Time Waveform Fitting with Background Noise Rejection

abstract: We present here a technique for developing a high-throughput algorithm to fit a combination of template pulse shapes while simultaneously subtracting parameterized background noise. By convolving the psuedoinverse of the least-squares fit design matrix along a regularly sampled waveform trace, the time evolution of the fit parameters for each basis function can be determined in real-time. We approximate these sliding linear fit response functions using piecewise polynomials, and develop an FPGA-friendly algorithm to be implemented in high sample-rate data acquisition systems. This is a robust universal filter that compares well to common filters optimized for energy calibration/resolution, as well as filters optimized for timing performance, even when significant noise components are present.


title: Gibbs Point Process Model for Young Star Clusters in M33

abstract: We demonstrate the power of Gibbs point process models from the spatial statistics literature when applied to studies of resolved galaxies. We conduct a rigorous analysis of the spatial distributions of objects in the star formation complexes of M33, including giant molecular clouds (GMCs) and young stellar cluster candidates (YSCCs). We choose a hierarchical model structure from GMCs to YSCCs based on the natural formation hierarchy between them. This approach circumvents the limitations of the empirical two-point correlation function analysis by naturally accounting for the inhomogeneity present in the distribution of YSCCs. We also investigate the effects of GMCs' properties on their spatial distributions. We confirm that the distribution of GMCs and YSCCs are highly correlated. We found that the spatial distributions of YSCCs reaches a peak of clustering pattern at ~250 pc scale compared to a Poisson process. This clustering mainly occurs in regions where the galactocentric distance >~4.5 kpc. Furthermore, the galactocentric distance of GMCs and their mass have strong positive effects on the correlation strength between GMCs and YSCCs. We outline some possible implications of these findings for our understanding of the cluster formation process.


title: Duals of non-zero square

abstract: In this short note, for each non-zero integer n, we construct a 4-manifold containing a smoothly concordant pair of spheres with a common dual of square n but no automorphism carrying one sphere to the other. Our examples, besides showing that the square zero assumption on the dual is necessary in Gabai's and Schneiderman-Teichner's versions of the 4D Light Bulb Theorem, have the interesting feature that both the Freedman-Quinn and Kervaire-Milnor invariant of the pair of spheres vanishes. The proof gives a surprising application of results due to Akbulut-Matveyev and Auckly-Kim-Melvin-Ruberman pertaining to the well-known Mazur cork.


title: A Simplistic Machine Learning Approach to Contact Tracing

abstract: This report is based on the modified NIST challenge, Too Close For Too Long, provided by the SFI Centre for Machine Learning (ML-Labs). The modified challenge excludes the time calculation (too long) aspect. By handcrafting features from phone instrumental data we develop two machine learning models, a GBM and an MLP, to estimate distance between two phones. Our method is able to outperform the leading NIST challenge result by the Hong Kong University of Science and Technology (HKUST) by a significant margin.


title: Too fast to grow: Dynamics of pendant drops sliding on a thin film

abstract: Pendant drops suspended on the underside of a wet substrate are known to accumulate fluid from the surrounding thin liquid film, a process that often results in dripping. The growth of such drops is hastened by their ability to translate over an otherwise uniform horizontal film. Here we show that this scenario is surprisingly reversed when the substrate is slightly tilted ($\approx 2$ deg) ; drops become too fast to grow and shrink over the course of their motion. Combining experiments and numerical simulations, we rationalize the transition between the conventional growth regime and the previously unknown decay regime we report. Using an analytical treatment of the Landau-Levich meniscus that connects the drop to the film, we quantitatively predict the drop dynamics in the two flow regimes and the value of the critical inclination angle where the transition between them occurs.


title: Convex Potential Flows: Universal Probability Distributions with Optimal Transport and Convex Optimization

abstract: Flow-based models are powerful tools for designing probabilistic models with tractable density. This paper introduces Convex Potential Flows (CP-Flow), a natural and efficient parameterization of invertible models inspired by the optimal transport (OT) theory. CP-Flows are the gradient map of a strongly convex neural potential function. The convexity implies invertibility and allows us to resort to convex optimization to solve the convex conjugate for efficient inversion. To enable maximum likelihood training, we derive a new gradient estimator of the log-determinant of the Jacobian, which involves solving an inverse-Hessian vector product using the conjugate gradient method. The gradient estimator has constant-memory cost, and can be made effectively unbiased by reducing the error tolerance level of the convex optimization routine. Theoretically, we prove that CP-Flows are universal density approximators and are optimal in the OT sense. Our empirical results show that CP-Flow performs competitively on standard benchmarks of density estimation and variational inference.


title: Greedy coordinate descent method on non-negative quadratic programming

abstract: The coordinate descent (CD) method has recently become popular for solving very large-scale problems, partly due to its simple update, low memory requirement, and fast convergence. In this paper, we explore the greedy CD on solving non-negative quadratic programming (NQP). The greedy CD generally has much more expensive per-update complexity than its cyclic and randomized counterparts. However, on the NQP, these three CDs have almost the same per-update cost, while the greedy CD can have significantly faster overall convergence speed. We also apply the proposed greedy CD as a subroutine to solve linearly constrained NQP and the non-negative matrix factorization. Promising numerical results on both problems are observed on instances with synthetic data and also image data.


title: On a radical extension of the field of rational functions in several variables

abstract: Let $F$ be a field and let $F(X_1,\dots,X_n)$ be the field of rational functions in $n$ variables $X_1,\dots,X_n$ over $F$. Let $T=X_1+\cdots+X_n\in F(X_1,\dots,X_n)$ and let $m$ be a positive integer such that $\text{char}\,F\nmid m$. Is it possible to express each $X_i$ as a rational function in $X_1^m\dots,X_n^m$ and $T$ over $F$? It is not difficult to prove that this can be done but it is another matter to show how this is done. We answer the above question affirmatively with a nonconstructive proof and a constructive proof.


title: PyLlama: a stable and versatile Python toolkit for the electromagnetic modeling of multilayered anisotropic media

abstract: PyLlama is a handy Python toolkit to compute the electromagnetic reflection and transmission properties of arbitrary multilayered linear media, including the case of anisotropy. Relying on a $4 \times 4$-matrix formalism, PyLlama implements not only the transfer matrix method, that is the most popular choice in existing codes, but also the scattering matrix method, which is numerically stable in all situations (e.g., thick, highly birefringent cholesteric structures at grazing incident angles). PyLlama is also designed to suit the practical needs by allowing the user to create, edit and assemble layers or multilayered domains with great ease. In this article, we present the electromagnetic theory underlying the transfer matrix and scattering matrix methods and outline the architecture and main features of PyLlama. Finally, we validate the code by comparison with available analytical solutions and demonstrate its versatility and numerical stability by modelling cholesteric media of varying complexity. A detailed documentation and tutorial are provided in a separate user manual. Applications of PyLlama range from the design of optical components to the modelling of polaritonic effects in polar crystals, to the study of structurally coloured materials in the living world.


title: Autonomous Cooperative Wall Building by a Team of Unmanned Aerial Vehicles in the MBZIRC 2020 Competition

abstract: This paper presents a system for autonomous cooperative wall building with a team of Unmanned Aerial Vehicles (UAVs). The system was developed for Challenge 2 of the Mohamed Bin Zayed International Robotics Challenge (MBZIRC) 2020. The wall-building scenario of Challenge 2 featured an initial stack of bricks and wall structure where the individual bricks had to be placed by a team of three UAVs. The objective of the task was to maximize collected points for placing the bricks within the restricted construction time while following the prescribed wall pattern. The proposed approach uses initial scanning to find a priori unknown locations of the bricks and the wall structure. Each UAV is then assigned to individual bricks and wall placing locations and further perform grasping and placement using onboard resources only. The developed system consists of methods for scanning a given area, RGB-D detection of bricks and wall placement locations, precise grasping and placing of bricks, and coordination of multiple UAVs. The paper describes the overall system, individual components, experimental verification in demanding outdoor conditions, the achieved results in the competition, and lessons learned. The presented CTU-UPenn-NYU approach achieved the overall best performance among all participants to won the MBZIRC competition by collecting the highest number of points by correct placement of a high number of bricks.


title: Neural network based image reconstruction with astrophysical priors

abstract: With the advent of interferometric instruments with 4 telescopes at the VLTI and 6 telescopes at CHARA, the scientific possibility arose to routinely obtain milli-arcsecond scale images of the observed targets. Such an image reconstruction process is typically performed in a Bayesian framework where the function to minimize is made of two terms: the datalikelihood and the Bayesian prior. This prior should be based on our prior knowledge of the observed source. Up to now,this prior was chosen from a set of generic and arbitrary functions, such as total variation for example. Here, we present an image reconstruction framework using generative adversarial networks where the Bayesian prior is defined using state-of-the-art radiative transfer models of the targeted objects. We validate this new image reconstruction algorithm on synthetic data with added noise. The generated images display a drastic reduction of artefacts and allow a more straight forward astrophysical interpretation. The results can be seen as a first illustration of how neural networks can provide significant improvements to the image reconstruction of a variety of astrophysical sources.


title: GNNUnlock: Graph Neural Networks-based Oracle-less Unlocking Scheme for Provably Secure Logic Locking

abstract: In this paper, we propose GNNUnlock, the first-of-its-kind oracle-less machine learning-based attack on provably secure logic locking that can identify any desired protection logic without focusing on a specific syntactic topology. The key is to leverage a well-trained graph neural network (GNN) to identify all the gates in a given locked netlist that belong to the targeted protection logic, without requiring an oracle. This approach fits perfectly with the targeted problem since a circuit is a graph with an inherent structure and the protection logic is a sub-graph of nodes (gates) with specific and common characteristics. GNNs are powerful in capturing the nodes' neighborhood properties, facilitating the detection of the protection logic. To rectify any misclassifications induced by the GNN, we additionally propose a connectivity analysis-based post-processing algorithm to successfully remove the predicted protection logic, thereby retrieving the original design. Our extensive experimental evaluation demonstrates that GNNUnlock is 99.24%-100% successful in breaking various benchmarks locked using stripped-functionality logic locking, tenacious and traceless logic locking, and Anti-SAT. Our proposed post-processing enhances the detection accuracy, reaching 100% for all of our tested locked benchmarks. Analysis of the results corroborates that GNNUnlock is powerful enough to break the considered schemes under different parameters, synthesis settings, and technology nodes. The evaluation further shows that GNNUnlock successfully breaks corner cases where even the most advanced state-of-the-art attacks fail.


title: Optimal selection of sample-size dependent common subsets of covariates for multi-task regression prediction

abstract: An analyst is given a training set consisting of regression datasets $D_j$ of different sizes, which are distributed according to some $G_j$, $j=1,\ldots,\cal J$, where the distributions $G_j$ are assumed to form a random sample generated by some common source. In particular, the $D_j$'s have a common set of covariates and they are all labeled. The training set is used by the analyst for selection of subsets of covariates denoted by ${P}^*(n)$, whose role is described next. The multi-task problem we consider is as follows: given a number of random labeled datasets (which may be in the training set or not) $D_{J_k}$ of size $n_k$, $k=1,\ldots,K$, estimate separately for each dataset the regression coefficients on the subset of covariates ${P}^*(n_k)$ and then predict future dependent variables given their covariates. Naturally, a large sample size $n_k$ of $D_{J_k}$ allows a larger subset of covariates, and the dependence of the size of the selected covariate subsets on $n_k$ is needed in order to achieve good prediction and avoid overfitting. Subset selection is notoriously difficult and computationally demanding, and requires large samples; using all the regression datasets in the training set together amounts to borrowing strength toward better selection under suitable assumptions. Furthermore, using common subsets for all regressions having a given sample size standardizes and simplifies the data collection and avoids having to select and use a different subset for each prediction task. Our approach is efficient when the relevant covariates for prediction are common to the different regressions, while the models' coefficients may vary between different regressions.


title: A case for robust translation tolerance in humans and CNNs. A commentary on Han et al

abstract: Han et al. (2020) reported a behavioral experiment that assessed the extent to which the human visual system can identify novel images at unseen retinal locations (what the authors call "intrinsic translation invariance") and developed a novel convolutional neural network model (an Eccentricity Dependent Network or ENN) to capture key aspects of the behavioral results. Here we show that their analysis of behavioral data used inappropriate baseline conditions, leading them to underestimate intrinsic translation invariance. When the data are correctly interpreted they show near complete translation tolerance extending to 14{\deg} in some conditions, consistent with earlier work (Bowers et al., 2016) and more recent work Blything et al. (in press). We describe a simpler model that provides a better account of translation invariance.


title: Strictly positive polynomials in the boundary of the SOS cone

abstract: We study the boundary of the cone of real polynomials that can be decomposed as a sum of squares (SOS) of real polynomials. This cone is included in the cone of nonnegative polynomials and both cones share a part of their boundary, which corresponds to polynomials that vanish at at least one point. We focus on the part of the boundary which is not shared, corresponding to strictly positive polynomials. For the cases of polynomials of degree 6 in 3 variables and degree 4 in 4 variables, this boundary has been completely characterized by G. Blekherman. For the cases of more variables or higher degree, results by G. Blekherman, R. Sinn and M. Velasco and other authors based on general conjectures give bounds for the maximum number of polynomials that can appear in a SOS decomposition and the maximum rank of the matrices in the Gram spectrahedron. Combining theoretical results and computational techniques, we compute examples that allow us to prove the optimality of the bounds for all degrees and number of variables. Additionally, we give examples for the following problems: examples in the boundary of the cone that are the sum of less than $n$ squares and have common complex roots, and examples of polynomials in the boundary with length larger than the expected from the dimension.


title: Conformally invariant differential operators on Heisenberg groups and minimal representations

abstract: For a simple real Lie group $G$ with Heisenberg parabolic subgroup $P$, we study the corresponding degenerate principal series representations. For a certain induction parameter the kernel of the conformally invariant system of second order differential operators constructed by Barchini, Kable and Zierau is a subrepresentation which turns out to be the minimal representation. To study this subrepresentation, we take the Heisenberg group Fourier transform in the non-compact picture and show that it yields a new realization of the minimal representation on a space of $L^2$-functions. The Lie algebra action is given by differential operators of order $\leq3$ and we find explicit formulas for the functions constituting the lowest $K$-type. These $L^2$-models were previously known for the groups $\operatorname{SO}(n,n)$, $E_{6(6)}$, $E_{7(7)}$ and $E_{8(8)}$ by Kazhdan and Savin, for the group $G_{2(2)}$ by Gelfand, and for the group $\widetilde{\operatorname{SL}}(3,\mathbb{R})$ by Torasso, using different methods. Our new approach provides a uniform and systematic treatment of these cases and also constructs new $L^2$-models for $E_{6(2)}$, $E_{7(-5)}$ and $E_{8(-24)}$ for which the minimal representation is a continuation of the quaternionic discrete series, and for the groups $\widetilde{\operatorname{SO}}(p,q)$ with either $p\geq q=3$ or $p,q\geq4$ and $p+q$ even. As a byproduct of our construction, we find an explicit formula for the group action of a non-trivial Weyl group element that, together with the simple action of a parabolic subgroup, generates $G$.


title: Electronic and Structural Properties of Janus SMoSe/MoX$_2$ (X=S,Se) In-plane Heterojunctions: A DFT Study

abstract: The electronic and structural properties of Janus MoSSe/MoX$_2$ (X=S,Se) in-plane heterojunctions, endowed with single-atom vacancies, were studied using density functional theory calculations. The stability of these structures was verified from cohesion energy calculations. Results showed that single-atom vacancies induce the appearance of flat midgap states, and a substantial amount of charge is localized in the vicinity of these defects. As a consequence, these heterojunctions presented an intrinsic dipole moment. No bond reconstructions were noted by removing an atom from the lattice, regardless of its chemical species. Our calculations predicted indirect electronic bandgap values between 1.6-1.7 eV.


title: A seamless, extended DG approach for advection-diffusion problems on unbounded domains

abstract: We propose and analyze a seamless extended Discontinuous Galerkin (DG) discretization of advection-diffusion equations on semi-infinite domains. The semi-infinite half line is split into a finite subdomain where the model uses a standard polynomial basis, and a semi-unbounded subdomain where scaled Laguerre functions are employed as basis and test functions. Numerical fluxes enable the coupling at the interface between the two subdomains in the same way as standard single domain DG interelement fluxes. A novel linear analysis on the extended DG model yields unconditional stability with respect to the P\'eclet number. Errors due to the use of different sets of basis functions on different portions of the domain are negligible, as highlighted in numerical experiments with the linear advection-diffusion and viscous Burgers' equations. With an added damping term on the semi-infinite subdomain, the extended framework is able to efficiently simulate absorbing boundary conditions without additional conditions at the interface. A few modes in the semi-infinite subdomain are found to suffice to deal with outgoing single wave and wave train signals more accurately than standard approaches at a given computational cost, thus providing an appealing model for fluid flow simulations in unbounded regions.


title: Contribution of Secondary Neutrinos from Line-of-sight Cosmic Ray Interactions to the IceCube Diffuse Astrophysical Flux

abstract: In ten years of observations, the IceCube neutrino observatory has revealed a neutrino sky in tension with previous expectations for neutrino point source emissions. Astrophysical objects associated with hadronic processes might act as production sites for neutrinos, observed as point sources at Earth. Instead, a nearly isotropic flux of astrophysical neutrinos is observed up to PeV energies, prompting a reassessment of the assumed transport and production physics. This work applies a new physical explanation for neutrino production from populations of active galactic nuclei (AGN) and starburst galaxies to three years of public IceCube point source data. Specifically, cosmic rays (CRs) produced at such sources might interact with extragalactic background light and gas along the line of sight, generating a secondary neutrino flux. This model is tested alongside a number of typical flux weighting schemes, in all cases the all-sky flux contribution being constrained to percent levels of the reported IceCube diffuse astrophysical flux.


title: Advances in secondary ion mass spectroscopy

abstract: Accurate SIMS measurement of nitrogen in niobium relies on the use of closely equivalent standards, made by ion implantation, to convert nitrogen signal intensity to nitrogen content by determination of relative sensitivity factors (RSF). Accurate RSF values for ppm-range nitrogen contents are increasingly critical, as more precision is sought in processes for next-generation superconducting radio-frequency (SRF) accelerator cavities. Factors influencing RSF value measurements were investigated with the aim of reliably attaining better than 10% accuracy in N concentrations at various depths into the bulk. This has been accomplished for materials typical of SRF cavities at the cost of great attention to all aspects.


title: Online Joint Topology Identification and Signal Estimation from Streams with Missing Data

abstract: Identifying the topology underlying a set of time series is useful for tasks such as prediction, denoising, and data completion. Vector autoregressive (VAR) model-based topologies capture dependencies among time series and are often inferred from observed spatio-temporal data. When data are affected by noise and/or missing samples, topology identification and signal recovery (reconstruction) tasks must be performed jointly. Additional challenges arise when i) the underlying topology is time-varying, ii) data become available sequentially, and iii) no delay is tolerated. This study proposes an online algorithm to overcome these challenges in estimating VAR model-based topologies, having constant complexity per iteration, which makes it interesting for big-data scenarios. The inexact proximal online gradient descent framework is used to derive a performance guarantee for the proposed algorithm, in the form of a dynamic regret bound. Numerical tests are also presented, showing the ability of the proposed algorithm to track time-varying topologies with missing data in an online fashion.


title: Multilingual Transfer Learning for QA Using Translation as Data Augmentation

abstract: Prior work on multilingual question answering has mostly focused on using large multilingual pre-trained language models (LM) to perform zero-shot language-wise learning: train a QA model on English and test on other languages. In this work, we explore strategies that improve cross-lingual transfer by bringing the multilingual embeddings closer in the semantic space. Our first strategy augments the original English training data with machine translation-generated data. This results in a corpus of multilingual silver-labeled QA pairs that is 14 times larger than the original training set. In addition, we propose two novel strategies, language adversarial training and language arbitration framework, which significantly improve the (zero-resource) cross-lingual transfer performance and result in LM embeddings that are less language-variant. Empirically, we show that the proposed models outperform the previous zero-shot baseline on the recently introduced multilingual MLQA and TyDiQA datasets.


title: Super-resolution Guided Pore Detection for Fingerprint Recognition

abstract: Performance of fingerprint recognition algorithms substantially rely on fine features extracted from fingerprints. Apart from minutiae and ridge patterns, pore features have proven to be usable for fingerprint recognition. Although features from minutiae and ridge patterns are quite attainable from low-resolution images, using pore features is practical only if the fingerprint image is of high resolution which necessitates a model that enhances the image quality of the conventional 500 ppi legacy fingerprints preserving the fine details. To find a solution for recovering pore information from low-resolution fingerprints, we adopt a joint learning-based approach that combines both super-resolution and pore detection networks. Our modified single image Super-Resolution Generative Adversarial Network (SRGAN) framework helps to reliably reconstruct high-resolution fingerprint samples from low-resolution ones assisting the pore detection network to identify pores with a high accuracy. The network jointly learns a distinctive feature representation from a real low-resolution fingerprint sample and successfully synthesizes a high-resolution sample from it. To add discriminative information and uniqueness for all the subjects, we have integrated features extracted from a deep fingerprint verifier with the SRGAN quality discriminator. We also add ridge reconstruction loss, utilizing ridge patterns to make the best use of extracted features. Our proposed method solves the recognition problem by improving the quality of fingerprint images. High recognition accuracy of the synthesized samples that is close to the accuracy achieved using the original high-resolution images validate the effectiveness of our proposed model.


title: 3D Scattering Tomography by Deep Learning with Architecture Tailored to Cloud Fields

abstract: We present 3DeepCT, a deep neural network for computed tomography, which performs 3D reconstruction of scattering volumes from multi-view images. Our architecture is dictated by the stationary nature of atmospheric cloud fields. The task of volumetric scattering tomography aims at recovering a volume from its 2D projections. This problem has been studied extensively, leading, to diverse inverse methods based on signal processing and physics models. However, such techniques are typically iterative, exhibiting high computational load and long convergence time. We show that 3DeepCT outperforms physics-based inverse scattering methods in term of accuracy as well as offering a significant orders of magnitude improvement in computational time. To further improve the recovery accuracy, we introduce a hybrid model that combines 3DeepCT and physics-based method. The resultant hybrid technique enjoys fast inference time and improved recovery performance.


title: The AFLOW Library of Crystallographic Prototypes: Part 3

abstract: The AFLOW Library of Crystallographic Prototypes has been extended to include a total of 1,100 common crystal structural prototypes (510 new ones with Part 3), comprising all of the inorganic crystal structures defined in the seven-volume Strukturbericht series published in Germany from 1937 through 1943. We cover a history of the Strukturbericht designation system, the evolution of the system over time, and the first comprehensive index of inorganic Strukturbericht designations ever published.


title: Rank-adaptive tensor methods for high-dimensional nonlinear PDEs

abstract: We present a new rank-adaptive tensor method to compute the numerical solution of high-dimensional nonlinear PDEs. The method combines functional tensor train (FTT) series expansions, operator splitting time integration, and a new rank-adaptive algorithm based on a thresholding criterion that limits the component of the PDE velocity vector normal to the FTT tensor manifold. This yields a scheme that can add or remove tensor modes adaptively from the PDE solution as time integration proceeds. The new method is designed to improve computational efficiency, accuracy and robustness in numerical integration of high-dimensional problems. In particular, it overcomes well-known computational challenges associated with dynamic tensor integration, including low-rank modeling errors and the need to invert covariance matrices of tensor cores at each time step. Numerical applications are presented and discussed for linear and nonlinear advection problems in two dimensions, and for a four-dimensional Fokker-Planck equation.


title: Four algorithms to solve symmetric multi-type non-negative matrix tri-factorization problem

abstract: In this paper, we consider the symmetric multi-type non-negative matrix tri-factorization problem (SNMTF), which attempts to factorize several symmetric non-negative matrices simultaneously. This can be considered as a generalization of the classical non-negative matrix tri-factorization problem and includes a non-convex objective function which is a multivariate sixth degree polynomial and a has convex feasibility set. It has a special importance in data science, since it serves as a mathematical model for the fusion of different data sources in data clustering. We develop four methods to solve the SNMTF. They are based on four theoretical approaches known from the literature: the fixed point method (FPM), the block-coordinate descent with projected gradient (BCD), the gradient method with exact line search (GM-ELS) and the adaptive moment estimation method (ADAM). For each of these methods we offer a software implementation: for the former two methods we use Matlab and for the latter Python with the TensorFlow library. We test these methods on three data-sets: the synthetic data-set we generated, while the others represent real-life similarities between different objects. Extensive numerical results show that with sufficient computing time all four methods perform satisfactorily and ADAM most often yields the best mean square error ($\mathrm{MSE}$). However, if the computation time is limited, FPM gives the best $\mathrm{MSE}$ because it shows the fastest convergence at the beginning. All data-sets and codes are publicly available on our GitLab profile.


title: Universal Collective Behavior of the Vibration Spectrum of Highly Connected Disordered Systems

abstract: We study small oscillations of highly connected systems, which represent a limit opposite to the more familiar case of disordered crystals. As a concrete example we analyze the vibrational spectra of composite pendula. Remarkably, these spectra exhibit universal behavior with non-vanishing zero-frequency limit of the density of phonon states. This universality is captured by a natural random matrix model for such systems. We analyze this model using S-transforms of free probability theory and obtain the density of modes explicitly, in the limit of large system size.

categories: math-ph cond-mat.dis-nn cond-mat.stat-mech hep-th math.MP


title: Analog Computation and Representation

abstract: Relative to digital computation, analog computation has been neglected in the philosophical literature. To the extent that attention has been paid to analog computation, it has been misunderstood. The received view -- that analog computation has to do essentially with continuity -- is simply wrong, as shown by careful attention to historical examples of discontinuous, discrete analog computers. Instead of the received view, I develop an account of analog computation in terms of a particular type of analog representation that allows for discontinuity. This account thus characterizes all types of analog computation, whether continuous or discrete. Furthermore, the structure of this account can be generalized to other types of computation: analog computation essentially involves analog representation, whereas digital computation essentially involves digital representation. Besides being a necessary component of a complete philosophical understanding of computation in general, understanding analog computation is important for computational explanation in contemporary neuroscience and cognitive science.


title: A tuning algorithm for a sliding mode controller of buildings with ATMD

abstract: This paper proposes an automatic tuning algorithm for a sliding mode controller (SMC) based on the Ackermann's formula, that attenuates the structural vibrations of a seismically excited building equipped with an Active Tuned Mass Damper (ATMD) mounted on its top floor. The switching gain and sliding surface of the SMC are designed through the proposed tuning algorithm to suppress the structural vibrations by minimizing either the top floor displacement or the control force applied to the ATMD. Moreover, the tuning algorithm selects the SMC parameters to guarantee the following closed-loop characteristics: 1) the transient responses of the structure and the ATMD are sufficiently fast and damped; and 2) the control force, as well as the displacements and velocities of the building and ATMD are within acceptable limits under the frequency band of the earthquake excitation. The proposed SMC shows robustness against the unmodeled dynamics such as the friction of the damper. Experimental results on a reduced scale structure permits demonstrating the efficiency of the tuning algorithm for the SMC, which is compared with the traditional Linear Quadratic Regulator (LQR) and with the Optimal Sliding Mode Controller (OSMC).


title: Bayesian nonstationary and nonparametric covariance estimation for large spatial data

abstract: In spatial statistics, it is often assumed that the spatial field of interest is stationary and its covariance has a simple parametric form, but these assumptions are not appropriate in many applications. Given replicate observations of a Gaussian spatial field, we propose nonstationary and nonparametric Bayesian inference on the spatial dependence. Instead of estimating the quadratic (in the number of spatial locations) entries of the covariance matrix, the idea is to infer a near-linear number of nonzero entries in a sparse Cholesky factor of the precision matrix. Our prior assumptions are motivated by recent results on the exponential decay of the entries of this Cholesky factor for Matern-type covariances under a specific ordering scheme. Our methods are highly scalable and parallelizable. We conduct numerical comparisons and apply our methodology to climate-model output, enabling statistical emulation of an expensive physical model.


title: Power prior models for treatment effect estimation in a small n, sequential, multiple assignment, randomized trial

abstract: A small n, sequential, multiple assignment, randomized trial (snSMART) is a small sample, two-stage design where participants receive up to two treatments sequentially, but the second treatment depends on response to the first treatment. The treatment effect of interest in an snSMART is the first-stage response rate, but outcomes from both stages can be used to obtain more information from a small sample. A novel way to incorporate the outcomes from both stages applies power prior models, in which first stage outcomes from an snSMART are regarded as the primary data and second stage outcomes are regarded as supplemental. We apply existing power prior models to snSMART data, and we also develop new extensions of power prior models. All methods are compared to each other and to the Bayesian joint stage model (BJSM) via simulation studies. By comparing the biases and the efficiency of the response rate estimates among all proposed power prior methods, we suggest application of Fisher's exact test or the Bhattacharyya's overlap measure to an snSMART to estimate the treatment effect in an snSMART, which both have performance mostly as good or better than the BJSM. We describe the situations where each of these suggested approaches is preferred.


title: K3 surfaces and cubic fourfolds with Abelian motive

abstract: Let $\mathcal{C}_d$ denote Hassett's Noether-Lefschetz divisor in the moduli space of cubic fourfolds given by cubics with intersection form of discriminant $d$. We exhibit arithmetic conditions on 20 indexes $d_1,\dots, d_{20}$ that ensure that the divisors $\mathcal{C}_{d_1},\dots,\mathcal{C}_{d_{20}}$ all intersect one another. This allows us to produce examples of rational cubic fourfolds with an associated K3 surface with any prescribed N\'eron-Severi rank. Then we apply this result to show that there exist one dimensional families of cubic fourfolds with finite dimensional Chow motive of Abelian type inside every divisor of special cubic fourfolds. This also implies Abelianity and finite dimensionality of the motive of related Hyperk\"ahler varieties, such as the Fano variety of lines and the LLSvS 8fold. A similar remark, plus some elementary Brill-Noether geometry of the associated K3 surfaces, allows us to show the Abelianity of the motive of an infinity of LSV 10folds, associated to Pfaffian cubic fourfolds. After that, starting from known 4-dimensional families of K3 surfaces, we construct two families $\mathcal{F}$ and $\mathcal{G}$ of cubic fourfolds whose motive is of Abelian type. Cubics from the first family are smooth, and their Chow motive is finite dimensional and Abelian. Those from the second family are singular, and their motives are Schur-finite and Abelian in Voevodsky's triangulated category of motives.


title: ALMA and IRIS Observations of the Solar Chromosphere II: Structure and Dynamics of Chromospheric Plage

abstract: We propose and employ a novel empirical method for determining chromospheric plage regions, which seems to better isolate plage from its surrounding regions compared to other methods commonly used. We caution that isolating plage from its immediate surroundings must be done with care in order to successfully mitigate statistical biases that, for instance, can impact quantitative comparisons between different chromospheric observables. Using this methodology, our analysis suggests that 1.25 mm wavelength free-free emission in plage regions observed with ALMA/Band6 may not form in the low chromosphere as previously thought, but rather in the upper chromospheric parts of dynamic plage features (such as spicules and other bright structures), i.e., near geometric heights of transition region temperatures. We investigate the high degree of similarity between chromospheric plage features observed in ALMA/Band6 (at 1.25 mm wavelength) and IRIS/Si IV at 1393\r{A}. We also show that IRIS/Mg II h and k is not as well correlated with ALMA/Band6 as was previously thought, and we discuss the discrepancies with previous works. Lastly, we report indications for chromospheric heating due to propagating shocks supported by the ALMA/Band6 observations.


title: Long-time behavior of solutions to the generalized Allen-Cahn model with degenerate diffusivity

abstract: The generalized Allen-Cahn equation, \[ u_t=\varepsilon^2(D(u)u_x)_x-\frac{\varepsilon^2}2D'(u)u_x^2-F'(u), \] with nonlinear diffusion, $D = D(u)$, and potential, $F = F(u)$, of the form \[ D(u) = |1-u^2|^{m}, \quad \text{or} \quad D(u) = |1-u|^{m}, \quad m >1, \] and \[ F(u)=\frac{1}{2n}|1-u^2|^{n}, \qquad n\geq2, \] respectively, is studied. These choices correspond to a reaction function that can be derived from a double well potential, and to a generalized degenerate diffusivity coefficient depending on the density $u$ that vanishes at one or at the two wells, $u = \pm 1$. It is shown that interface layer solutions that are equal to $\pm 1$ except at a finite number of thin transitions of width $\varepsilon$ persist for an either exponentially or algebraically long time, depending upon the interplay between the exponents $n$ and $m$. For that purpose, energy bounds for a renormalized effective energy potential of Ginzburg-Landau type are derived.


title: Self-adjoint Laplacians and symmetric diffusions on hyperbolic attractors

abstract: We construct self-adjoint Laplacians and symmetric Markov semigroups on hyperbolic attractors, endowed with Gibbs $u$-measures. If the measure has full support, we can also conclude the existence of an associated symmetric diffusion process. In the special case of partially hyperbolic diffeomorphisms induced by geodesic flows on negatively curved manifolds the Laplacians we consider are self-adjoint extensions of well-known classical leafwise Laplacians. We observe a quasi-invariance property of energy densities in the $u$-conformal case and the existence of nonconstant functions of zero energy.


title: Clustering multivariate functional data using unsupervised binary trees

abstract: We propose a model-based clustering algorithm for a general class of functional data for which the components could be curves or images. The random functional data realizations could be measured with error at discrete, and possibly random, points in the definition domain. The idea is to build a set of binary trees by recursive splitting of the observations. The number of groups are determined in a data-driven way. The new algorithm provides easily interpretable results and fast predictions for online data sets. Results on simulated datasets reveal good performance in various complex settings. The methodology is applied to the analysis of vehicle trajectories on a German roundabout.


title: Path Optimization Sheaves

abstract: Motivated by efforts to incorporate sheaves into networking, we seek to reinterpret pathfinding algorithms in terms of cellular sheaves, using Dijkstra's algorithm as an example. We construct sheaves on a graph with distinguished source and sink vertices, in which paths are represented by sections. The first sheaf is a very general construction that can be applied to other algorithms, while the second is created specifically to capture the decision making of Dijkstra's algorithm. In both cases, Dijkstra's algorithm can be described as a systematic process of extending local sections to global sections. We discuss the relationship between the two sheaves and summarize how other pathfinding algorithms can be interpreted in a similar way. While the sheaves presented here address paths and pathfinding algorithms, we suggest that future work could explore connections to other concepts from graph theory and other networking algorithms. This work was supported by the NASA Internship Project and SCaN Internship Project during the summer of 2020.


title: Image-Graph-Image Translation via Auto-Encoding

abstract: This work presents the first convolutional neural network that learns an image-to-graph translation task without needing external supervision. Obtaining graph representations of image content, where objects are represented as nodes and their relationships as edges, is an important task in scene understanding. Current approaches follow a fully-supervised approach thereby requiring meticulous annotations. To overcome this, we are the first to present a self-supervised approach based on a fully-differentiable auto-encoder in which the bottleneck encodes the graph's nodes and edges. This self-supervised approach can currently encode simple line drawings into graphs and obtains comparable results to a fully-supervised baseline in terms of F1 score on triplet matching. Besides these promising results, we provide several directions for future research on how our approach can be extended to cover more complex imagery.


title: Lowering the scale of Pati-Salam breaking through seesaw mixing

abstract: We analyse the experimental limits on the breaking scale of Pati-Salam extensions of the Standard Model. These arise from the experimental limits on rare-meson decay processes mediated at tree-level by the vector leptoquark in the model. This leptoquark ordinarily couples to to both left- and right-handed SM fermions and therefore the meson decays do not experience a helicity suppression. We find that the current limits vary from $\mathcal{O}(80-2500)$ TeV depending on the choice of matrix structure appearing in the relevant three-generational charged-current interactions. We extensively analyse scenarios where additional fermionic degrees of freedom are introduced, transforming as complete Pati-Salam multiplets. These can lower the scales of Pati-Salam breaking through mass-mixing within the charged-lepton and down-quark sectors, leading to a helicity suppression of the meson decay widths which constrain Pati-Salam breaking. We find four multiplets with varying degrees of viability for this purpose: an $SU(2)_{L/R}$ bidoublet, a pair of $SU(4)$ decuplets and either a $SU(2)_L$ or $SU(2)_R$ triplet all of which contain heavy exotic versions of the SM charged leptons. We find that the Pati-Salam limits can be as low as $\mathcal{O}(5-150)$ TeV with the addition of these four multiplets. We also identify an interesting possible connection between the smallness of the neutrino masses and a helicity suppression of the Pati-Salam limits for three of the four multiplets.


title: Rediscussion of eclipsing binaries. Paper II. The eccentric solar-type system KX Cancri

abstract: KX Cancri is an eclipsing binary containing two G-type stars with an orbital period of 31.2 d and an eccentricity of 0.47. These qualities make it a promising candidate for a benchmark solar-type binary system. We analyse the first light curve of this system to have complete coverage of both primary and secondary eclipses, obtained using the Transiting Exoplanet Survey Satellite (TESS). We augment these data with published radial velocities and measure the masses to be 1.134 +/- 0.003 Msun and 1.124 +/- 0.005 Msun and the radii to be 1.053 +/- 0.006 Rsun and 1.059 +/- 0.005 Rsun. A ratio of the radii near unity is strongly preferred by the TESS data, in contrast to existing ground-based light curves. The distance to the system measured from the radii and Teff values of the stars agrees well with the trigonometric parallax from the Gaia satellite. The properties of the system are consistent with theoretical predictions for a super-solar metallicity and an age of 1.0 to 1.5 Gyr. A detailed analysis of the photospheric properties of the stars based on high-resolution spectra is encouraged.


title: Rediscussion of eclipsing binaries. Paper III. The interferometric, spectroscopic and eclipsing binary V1022 Cassiopeiae

abstract: V1022 Cas has been known as a spectroscopic binary for a century. It was found to be eclipsing based on photometry from the Hipparcos satellite, and an astrometric orbit was recently obtained from near-infrared interferometry. We present the first high-precision measurement of the radii of the stars based on light curves obtained by the TESS satellite. Combined with published radial velocities from high-resolution spectra, we measure the masses of the stars to be 1.626 +/- 0.001 Msun and 1.609 +/- 0.001 Msun, and the radii to be 2.591 +/- 0.026 Rsun and 2.472 +/- 0.027 Rsun. The 12.16-d orbit is eccentric and the stars rotate sub-synchronously, so the system is tidally unevolved. A good match to these masses and radii, and published temperatures of the stars, is found for several sets of theoretical stellar evolutionary models, for a solar metallicity and an age of approximately 2 Gyr. Four separate distance determinations to the system are available, and are in good agreement. The distances are based on surface brightness calibrations, theoretical bolometric corrections, the Gaia parallax, and the angular size of the astrometric orbit. A detailed spectroscopic analysis of the system to measure chemical abundances and more precise temperatures would be helpful.


title: Analyzing Growth Kinematics and Fractal Dimensions of Molybdenum Disulfide Films

abstract: Through the positive role of alkali halides in realizing large area growth of transition metal-di-chalcogenide layers has been validated, the film-growth kinematics has not yet been fully established. This work presents a systematic analysis of the MoS$_2$ morphology for films grown under various pre-treatment conditions of the substrate with sodium chloride (NaCl). At an optimum NaCl concentration, the domain size of the monolayer increased by almost two orders of magnitude compared to alkali-free growth of MoS$_2$. The results show an inverse relationship between fractal dimension and areal coverage of the substrate with monolayers and multi-layers, respectively. Using the Fact-Sage software, the role of NaCl in determining the partial pressures of Mo- and S-based compounds in gaseous phase at the growth temperature is elucidated. The presence of alkali salts is shown to affect the domain size and film morphology by affecting the Mo and S partial pressures. Compared to alkali-free synthesis under the same growth conditions, MoS$_2$ film growth assisted by NaCl results in $\approx$ 81$\%$ of the substrate covered by monolayers. Under ideal growth conditions, at an optimum NaCl concentration, nucleation was suppressed, and domains enlarged, resulting in large area growth of MoS$_2$ monolayers. The monolayers were found to be free of unintentional doping with alkali metal and halogen atoms and exhibit high crystallinity and excellent opto-electronic quality.


title: CommPOOL: An Interpretable Graph Pooling Framework for Hierarchical Graph Representation Learning

abstract: Recent years have witnessed the emergence and flourishing of hierarchical graph pooling neural networks (HGPNNs) which are effective graph representation learning approaches for graph level tasks such as graph classification. However, current HGPNNs do not take full advantage of the graph's intrinsic structures (e.g., community structure). Moreover, the pooling operations in existing HGPNNs are difficult to be interpreted. In this paper, we propose a new interpretable graph pooling framework - CommPOOL, that can capture and preserve the hierarchical community structure of graphs in the graph representation learning process. Specifically, the proposed community pooling mechanism in CommPOOL utilizes an unsupervised approach for capturing the inherent community structure of graphs in an interpretable manner. CommPOOL is a general and flexible framework for hierarchical graph representation learning that can further facilitate various graph-level tasks. Evaluations on five public benchmark datasets and one synthetic dataset demonstrate the superior performance of CommPOOL in graph representation learning for graph classification compared to the state-of-the-art baseline methods, and its effectiveness in capturing and preserving the community structure of graphs.


title: Certifying Incremental Quadratic Constraints for Neural Networks via Convex Optimization

abstract: Abstracting neural networks with constraints they impose on their inputs and outputs can be very useful in the analysis of neural network classifiers and to derive optimization-based algorithms for certification of stability and robustness of feedback systems involving neural networks. In this paper, we propose a convex program, in the form of a Linear Matrix Inequality (LMI), to certify incremental quadratic constraints on the map of neural networks over a region of interest. These certificates can capture several useful properties such as (local) Lipschitz continuity, one-sided Lipschitz continuity, invertibility, and contraction. We illustrate the utility of our approach in two different settings. First, we develop a semidefinite program to compute guaranteed and sharp upper bounds on the local Lipschitz constant of neural networks and illustrate the results on random networks as well as networks trained on MNIST. Second, we consider a linear time-invariant system in feedback with an approximate model predictive controller parameterized by a neural network. We then turn the stability analysis into a semidefinite feasibility program and estimate an ellipsoidal invariant set for the closed-loop system.


title: Making Cross-Domain Recommendations by Associating Disjoint Users and Items Through the Affective Aware Pseudo Association Method

abstract: This paper utilizes an ingenious text-based affective aware pseudo association method (AAPAM) to link disjoint users and items across different information domains and leverage them to make cross-domain content-based and collaborative filtering recommendations. This paper demonstrates that the AAPAM method could seamlessly join different information domain datasets to act as one without any additional cross-domain information retrieval protocols. Besides making cross-domain recommendations, the benefit of joining datasets from different information domains through AAPAM is that it eradicates cold start issues while making serendipitous recommendations.


title: Towards Neural Programming Interfaces

abstract: It is notoriously difficult to control the behavior of artificial neural networks such as generative neural language models. We recast the problem of controlling natural language generation as that of learning to interface with a pretrained language model, just as Application Programming Interfaces (APIs) control the behavior of programs by altering hyperparameters. In this new paradigm, a specialized neural network (called a Neural Programming Interface or NPI) learns to interface with a pretrained language model by manipulating the hidden activations of the pretrained model to produce desired outputs. Importantly, no permanent changes are made to the weights of the original model, allowing us to re-purpose pretrained models for new tasks without overwriting any aspect of the language model. We also contribute a new data set construction algorithm and GAN-inspired loss function that allows us to train NPI models to control outputs of autoregressive transformers. In experiments against other state-of-the-art approaches, we demonstrate the efficacy of our methods using OpenAI's GPT-2 model, successfully controlling noun selection, topic aversion, offensive speech filtering, and other aspects of language while largely maintaining the controlled model's fluency under deterministic settings.


title: Sums of four and more unit fractions and approximate parametrizations

abstract: We prove new upper bounds on the number of representations of rational numbers $\frac{m}{n}$ as a sum of $4$ unit fractions, giving five different regions, depending on the size of $m$ in terms of $n$. In particular, we improve the most relevant cases, when $m$ is small, and when $m$ is close to $n$. The improvements stem from not only studying complete parametrizations of the set of solutions, but simplifying this set appropriately. Certain subsets of all parameters define the set of all solutions, up to applications of divisor functions, which has little impact on the upper bound of the number of solutions. These "approximate parametrizations" were the key point to enable computer programmes to filter through large number of equations and inequalities. Furthermore, this result leads to new upper bounds for the number of representations of rational numbers as sums of more than $4$ unit fractions.


title: A Comment on and Correction to: Opinion dynamics in the presence of increasing agreement pressure

abstract: We identify a counter-example to the consensus result given in [J. Semonsen et al. Opinion dynamics in the presence of increasing agreement pressure. \textit{IEEE Trans. Cyber.}, 49(4): 1270-1278, 2018]. We resolve the counter-example by replacing Lemma 5 in the given reference with a novel variation of the Banach Fixed Point theorem which explains both the numerical results in the reference and the counter-example(s) in this note, and provides a sufficient condition for consensus in systems with increasing peer-pressure. This work is relevant for other papers that have used the proof technique from Semonsen et al. and establishes the veracity of their claims assuming the new sufficient condition.


title: Entanglement of graph states of spin system with Ising interaction and its quantifying on IBM's quantum computer

abstract: We consider graph states generated by operator of evolution with Ising Hamiltonian. The geometric measure of entanglement of a spin with other spins in the graph state is obtained analytically and quantified on IBM's quantum computer, IBM Q Valencia. The results of quantum calculations are in good agreement with the theoretical ones. We conclude that the geometric measure of entanglement of a spin with other spins in the graph state is related with degree of vertex representing the spin in the corresponding graph.


title: Guiding Development Work Across a Software Ecosystem by Visualizing Usage Data

abstract: Software is increasingly produced in the form of ecosystems, collections of interdependent components maintained by a distributed community. These ecosystems act as network organizations, not markets, and thus often lack actionable price-like signals about how the software is used and what impact it has. We introduce a tool, the Scientific Software Network Map, that collects and displays summarized usage data tailored to the needs of actors in software ecosystems. We performed a contextualized walkthrough of the Map with producers and stewards in six scientific software ecosystems that use the R language. We found that they work to maximize diversity rather than quantity of uses, and to minimize coordination costs. We also found that summarized usage data would be useful for justifying ecosystem work to funding agencies; and we discovered a variety of more granular usage needs that would help in adding or maintaining features.


title: The International Virtual Observatory Alliance (IVOA) in 2020

abstract: The International Virtual Observatory Alliance (IVOA) develops the technical standards needed for seamless discovery of and access to astronomy data worldwide, according to the Findable, Accessible, Interoperable and Reusable (FAIR) principles, with the goal of realizing the Virtual Observatory (VO). There are 21 member organizations. The Netherlands VO applied for membership in 2020. Astronomical communities from other nations have shown interest in joining the IVOA. This paper describes the activities of the IVOA in 2020 and summarizes the May and November 2020 virtual "interoperability meetings." The May meeting was the first to be held online and the first to have over 200 registrants.


title: On the modeling of thermal and free carrier nonlinearities in Silicon On Insulator microring resonators

abstract: The temporal dynamics of integrated silicon resonators has been modeled using a set of equations coupling the internal energy, the temperature and the free carrier population. Owing to its simplicity, Newton's law of cooling is the traditional choice for describing the thermal evolution of such systems. In this work, we theoretically and experimentally prove that this can be inadequate in monolithic planar devices, leading to inaccurate predictions. A new equation, that we train to reproduce the correct temperature behaviour, is introduced to fix the discrepancies with the experimental results. We discuss the limitations and the range of validity of our refined model, identifying those cases where Netwon's law provides, nevertheless, accurate solutions. Our modeling describes the phenomena underlying thermal and free carrier instabilities, and is a valuable tool for the engineering of photonic systems which relay on resonator dynamical states, such as all optical spiking neural networks or reservoirs for neuromorphic computing.


title: A Generative Approach for Detection-driven Underwater Image Enhancement

abstract: In this paper, we introduce a generative model for image enhancement specifically for improving diver detection in the underwater domain. In particular, we present a model that integrates generative adversarial network (GAN)-based image enhancement with the diver detection task. Our proposed approach restructures the GAN objective function to include information from a pre-trained diver detector with the goal to generate images which would enhance the accuracy of the detector in adverse visual conditions. By incorporating the detector output into both the generator and discriminator networks, our model is able to focus on enhancing images beyond aesthetic qualities and specifically to improve robotic detection of scuba divers. We train our network on a large dataset of scuba divers, using a state-of-the-art diver detector, and demonstrate its utility on images collected from oceanic explorations of human-robot teams. Experimental evaluations demonstrate that our approach significantly improves diver detection performance over raw, unenhanced images, and even outperforms detection performance on the output of state-of-the-art underwater image enhancement algorithms. Finally, we demonstrate the inference performance of our network on embedded devices to highlight the feasibility of operating on board mobile robotic platforms.


title: A general framework for multimode Gaussian quantum optics and photo-detection: application to Hong-Ou-Mandel interference with filtered heralded single photon sources

abstract: The challenging requirements of large scale quantum information processing using parametric heralded single photon sources involves maximising the interference visibility whilst maintaining an acceptable photon generation rate. By developing a general theoretical framework that allows us to include large numbers of spatial and spectral modes together with linear and non-linear optical elements, we investigate the combined effects of spectral and photon number impurity on the measured Hong--Ou--Mandel interference visibility of parametric photon sources, considering both threshold and number resolving detectors, together with the effects of spectral filtering. We find that for any degree of spectral impurity, increasing the photon generation rate necessarily decreases the interference visibility, even when using number resolving detection. While tight spectral filtering can be used to enforce spectral purity and increased interference visibility at low powers, we find that the induced photon number impurity results in a decreasing interference visibility and heralding efficiency with pump power, while the maximum generation rate is also reduced.


title: Detection of charge states of an InAs nanowire triple quantum dot with an integrated nanowire charge sensor

abstract: A linear triple quantum dot (TQD) integrated with a quantum dot (QD) charge sensor is realized. The TQD and the charge sensor are built from two adjacent InAs nanowires by fine finger gate technique. The charge state configurations of the nanowire TQD are studied by measurements of the direct transport signals of the TQD and by detection of the charge state transitions in the TQD via the nanowire QD sensor. Excellent agreements in the charge stability diagrams of the TQD obtained by the direct transport measurements and by the charge-state transition detection measurements are achieved. It is shown that the charge stability diagrams are featured by three groups of charge state transition lines of different slopes, corresponding to the changes in the electron occupation numbers of the three individual QDs in the TQD. It is also shown that the integrated nanowire QD sensor is highly sensitive and can detect the charge state transitions in the cases where the direct transport signals of the TQD are too weak to be measurable. Tuning to a regime, where all the three QDs in the TQD are close to be on resonance with the Fermi level of the source and drain reservoirs and co-existence of triple and quadruple points becomes possible, has also been demonstrated with the help of the charge sensor in the region where the direct transport signals of the TQD are hardly visible.


title: Three perfect mapping class groups

abstract: We prove that the mapping class group of a surface obtained from removing a Cantor set from either the 2-sphere, the plane, or the interior of the closed 2-disk has no proper countable-index subgroups. The proof is an application of the automatic continuity of these groups, which was established by Mann. As corollaries, we see that these groups do not contain any proper finite-index subgroups and that each of these groups have trivial abelianization.


title: Some remarks on steady solutions to the Euler system in $\mathbb{R}^d$

abstract: We show that the Euler system of gas dynamics in $\mathbb{R}^d$, $d=2,3$, with positive far field density and arbitrary far field entropy, admits infinitely many steady solutions with compactly supported velocity. The same proof yields a similar result for the incompressible Euler system with variable density. In particular, these are examples of global in time smooth (non-trivial) solutions for the corresponding time-dependent systems.


title: Towards a robust out-of-the-box neural network model for genomic data

abstract: The accurate prediction of biological features from genomic data is paramount for precision medicine and sustainable agriculture. For decades, neural network models have been widely popular in fields like computer vision, astrophysics and targeted marketing given their prediction accuracy and their robust performance under big data settings. Yet neural network models have not made a successful transition into the medical and biological world due to the ubiquitous characteristics of biological data such as modest sample sizes, sparsity, and extreme heterogeneity. Here, we investigate the robustness, generalization potential and prediction accuracy of widely used convolutional neural network and natural language processing models with a variety of heterogeneous genomic datasets. Mainly, recurrent neural network models outperform convolutional neural network models in terms of prediction accuracy, overfitting and transferability across the datasets under study. While the perspective of a robust out-of-the-box neural network model is out of reach, we identify certain model characteristics that translate well across datasets and could serve as a baseline model for translational researchers.


title: How to enhance quantum generative adversarial learning of noisy information

abstract: Quantum Machine Learning is where nowadays machine learning meets quantum information science. In order to implement this new paradigm for novel quantum technologies, we still need a much deeper understanding of its underlying mechanisms, before proposing new algorithms to feasibly address real problems. In this context, quantum generative adversarial learning is a promising strategy to use quantum devices for quantum estimation or generative machine learning tasks. However, the convergence behaviours of its training process, which is crucial for its practical implementation on quantum processors, have not been investigated in detail yet. Indeed here we show how different training problems may occur during the optimization process, such as the emergence of limit cycles. The latter may remarkably extend the convergence time in the scenario of mixed quantum states playing a crucial role in the already available noisy intermediate scale quantum devices. Then, we propose new strategies to achieve a faster convergence in any operating regime. Our results pave the way for new experimental demonstrations of such hybrid classical-quantum protocols allowing to evaluate the potential advantages over their classical counterparts.


title: Strong Admissibility for Abstract Dialectical Frameworks

abstract: Abstract dialectical frameworks (ADFs) have been introduced as a formalism for modeling and evaluating argumentation allowing general logical satisfaction conditions. Different criteria used to settle the acceptance of arguments are called semantics. Semantics of ADFs have so far mainly been defined based on the concept of admissibility. However, the notion of strongly admissible semantics studied for abstract argumentation frameworks has not yet been introduced for ADFs. In the current work we present the concept of strong admissibility of interpretations for ADFs. Further, we show that strongly admissible interpretations of ADFs form a lattice with the grounded interpretation as top element.


title: Lifting Lockdown Control Measure Assessment: From Finite to Infinite-dimensional Epidemic Models for COVID-19

abstract: The main focus of this chapter is on public health control strategies which are currently the main way to mitigate COVID-19 pandemic. We introduce and compare compartmental models of increasing complexity for COVID-19 transmission to describe dynamics of the disease spread. We begin by considering an SEAIR model including basic characteristics related to COVID-19. Next, we shall pay attention to age-structure modeling to emphasis the role of age-group individuals on the disease spread. A Model with constant delay is also formulated to show the impact of the latency period on the severity of COVID-19. Since there is evidence that for COVID-19 disease, important relationships exist between what is happening in the host and what is occurring at the population level, we shall link the basic model to in-host dynamics through the so-called threshold-type delay models. Finally, we will include demographic effects to the most complex models and we will conduct rigorous bifurcation analysis to quantify possible factors responsible for disease progression.


title: An IoT Framework for Heart Disease Prediction based on MDCNN Classifier

abstract: Nowadays, heart disease is the leading cause of death worldwide. Predicting heart disease is a complex task since it requires experience along with advanced knowledge. Internet of Things (IoT) technology has lately been adopted in healthcare systems to collect sensor values for heart disease diagnosis and prediction. Many researchers have focused on the diagnosis of heart disease, yet the accuracy of the diagnosis results is low. To address this issue, an IoT framework is proposed to evaluate heart disease more accurately using a Modified Deep Convolutional Neural Network (MDCNN). The smartwatch and heart monitor device that is attached to the patient monitors the blood pressure and electrocardiogram (ECG). The MDCNN is utilized for classifying the received sensor data into normal and abnormal. The performance of the system is analyzed by comparing the proposed MDCNN with existing deep learning neural networks and logistic regression. The results demonstrate that the proposed MDCNN based heart disease prediction system performs better than other methods. The proposed method shows that for the maximum number of records, the MDCNN achieves an accuracy of 98.2 which is better than existing classifiers.


title: The Three Ghosts of Medical AI: Can the Black-Box Present Deliver?

abstract: Our title alludes to the three Christmas ghosts encountered by Ebenezer Scrooge in \textit{A Christmas Carol}, who guide Ebenezer through the past, present, and future of Christmas holiday events. Similarly, our article will take readers through a journey of the past, present, and future of medical AI. In doing so, we focus on the crux of modern machine learning: the reliance on powerful but intrinsically opaque models. When applied to the healthcare domain, these models fail to meet the needs for transparency that their clinician and patient end-users require. We review the implications of this failure, and argue that opaque models (1) lack quality assurance, (2) fail to elicit trust, and (3) restrict physician-patient dialogue. We then discuss how upholding transparency in all aspects of model design and model validation can help ensure the reliability of medical AI.


title: Sketchy With a Chance of Adoption: Can Sketch-Based Telemetry Be Ready for Prime Time?

abstract: Sketching algorithms or sketches have emerged as a promising alternative to the traditional packet sampling-based network telemetry solutions. At a high level, they are attractive because of their high resource efficiency and accuracy guarantees. While there have been significant recent advances in various aspects of sketching for networking tasks, many fundamental challenges remain unsolved that are likely stumbling blocks for adoption. Our contribution in this paper is in identifying and formulating these research challenges across the ecosystem encompassing network operators, platform vendors/developers, and algorithm designers. We hope that these serve as a necessary fillip for the community to enable the broader adoption of sketch-based telemetry.


title: Overcoming Catastrophic Forgetting in Graph Neural Networks

abstract: Catastrophic forgetting refers to the tendency that a neural network "forgets" the previous learned knowledge upon learning new tasks. Prior methods have been focused on overcoming this problem on convolutional neural networks (CNNs), where the input samples like images lie in a grid domain, but have largely overlooked graph neural networks (GNNs) that handle non-grid data. In this paper, we propose a novel scheme dedicated to overcoming catastrophic forgetting problem and hence strengthen continual learning in GNNs. At the heart of our approach is a generic module, termed as topology-aware weight preserving~(TWP), applicable to arbitrary form of GNNs in a plug-and-play fashion. Unlike the main stream of CNN-based continual learning methods that rely on solely slowing down the updates of parameters important to the downstream task, TWP explicitly explores the local structures of the input graph, and attempts to stabilize the parameters playing pivotal roles in the topological aggregation. We evaluate TWP on different GNN backbones over several datasets, and demonstrate that it yields performances superior to the state of the art. Code is publicly available at \url{https://github.com/hhliu79/TWP}.


title: Nonlinear Regression with a Convolutional Encoder-Decoder for Remote Monitoring of Surface Electrocardiograms

abstract: We propose the Nonlinear Regression Convolutional Encoder-Decoder (NRCED), a novel framework for mapping a multivariate input to a multivariate output. In particular, we implement our algorithm within the scope of 12-lead surface electrocardiogram (ECG) reconstruction from intracardiac electrograms (EGM) and vice versa. The goal of performing this task is to allow for improved point-of-care monitoring of patients with an implanted device to treat cardiac pathologies. We will achieve this goal with 12-lead ECG reconstruction and by providing a new diagnostic tool for classifying atypical heartbeats. The algorithm is evaluated on a dataset retroactively collected from 14 patients. Correlation coefficients calculated between the reconstructed and the actual ECG show that the proposed NRCED method represents an efficient, accurate, and superior way to synthesize a 12-lead ECG. We can also achieve the same reconstruction accuracy with only one EGM lead as input. We also tested the model in a non-patient specific way and saw a reasonable correlation coefficient. The model was also executed in the reverse direction to produce EGM signals from a 12-lead ECG and found that the correlation was comparable to the forward direction. Lastly, we analyzed the features learned in the model and determined that the model learns an overcomplete basis of our 12-lead ECG space. We then use this basis of features to create a new diagnostic tool for identifying atypical and diseased heartbeats. This resulted in a ROC curve with an associated area under the curve value of 0.98, demonstrating excellent discrimination between the two classes.


title: Asymptotic hollowness of lattice simplices

abstract: An $(n-1)$-tuple $a = (a(1), \dots, a(n-1))$ consisting of positive integers is said to be asymptotically hollow if there exist infinitely many positive integers $N$ such that the convex hull, $K(a(n))$, in $n$-dimensional Euclidean space of $\{ 0,e_1, \dots, e_{n-1}, \alpha(N)^T\}$ is hollow (has no lattice points in its interior), where $e_i$ run over all but the last standard basis elements, and $\alpha(N) $ is the row $(a(1), \dots, a(N-1), N)$. The tuple is trivial if $\min a(i) = 1$. Nontrivial asymptotically hollow tuples are characterized in terms of modular inequalities, and turn out to be rare. We show that for a tuple $a$, there exists an effectively computable constant $C$ (depending on $a$) such that if for some $N > C$, $K(\alpha(N))$ is (not) hollow, then for all $M > C$, $K(\alpha(M))$ is (not) hollow (respectively). When $n = 4$, the nontrivial asymptotically hollow triples are completely determined; there are eleven of them, together with a one-parameter family.


title: Learning to Resolve Conflicts for Multi-Agent Path Finding with Conflict-Based Search

abstract: Conflict-Based Search (CBS) is a state-of-the-art algorithm for multi-agent path finding. At the high level, CBS repeatedly detects conflicts and resolves one of them by splitting the current problem into two subproblems. Previous work chooses the conflict to resolve by categorizing the conflict into three classes and always picking a conflict from the highest-priority class. In this work, we propose an oracle for conflict selection that results in smaller search tree sizes than the one used in previous work. However, the computation of the oracle is slow. Thus, we propose a machine-learning framework for conflict selection that observes the decisions made by the oracle and learns a conflict-selection strategy represented by a linear ranking function that imitates the oracle's decisions accurately and quickly. Experiments on benchmark maps indicate that our method significantly improves the success rates, the search tree sizes and runtimes over the current state-of-the-art CBS solver.


title: xRAI: Explainable Representations through AI

abstract: We present xRAI an approach for extracting symbolic representations of the mathematical function a neural network was supposed to learn from the trained network. The approach is based on the idea of training a so-called interpretation network that receives the weights and biases of the trained network as input and outputs the numerical representation of the function the network was supposed to learn that can be directly translated into a symbolic representation. We show that interpretation nets for different classes of functions can be trained on synthetic data offline using Boolean functions and low-order polynomials as examples. We show that the training is rather efficient and the quality of the results are promising. Our work aims to provide a contribution to the problem of better understanding neural decision making by making the target function explicit


title: An Atmospheric Dispersion Corrector Design with Milliarcsecond-Level Precision from 1 to 4 microns for High Dispersion Coronagraphy

abstract: Differential atmospheric refraction (DAR) limits the amount of light that can be coupled into a single mode fiber and provides additional complications for any fiber tracking system. We present an atmospheric dispersion corrector (ADC) design based off of two counter-rotating prisms to fit the needs of exoplanet spectroscopy for the Keck Planet Imager and Characterizer (KPIC) from 1.1 to 4.2 microns. Due to strong telluric effects, we find that the default Zemax prescription for DAR between 2 and 4.2 microns to be inaccurate up to 15 mas when comparing against DAR models computed from first principles. Using first-principle models, we developed our own custom ADC optimization solution and achieve less than 4 mas residual dispersion in any individual science band (J, K, L) down to 60 degree zenith angles, while the whole time maintaining less than 3 mas of residual dispersion in the tracking band (H) and less than 2 mas of residual dispersion between the tracking and science bands.


title: Price Suggestion for Online Second-hand Items with Texts and Images

abstract: This paper presents an intelligent price suggestion system for online second-hand listings based on their uploaded images and text descriptions. The goal of price prediction is to help sellers set effective and reasonable prices for their second-hand items with the images and text descriptions uploaded to the online platforms. Specifically, we design a multi-modal price suggestion system which takes as input the extracted visual and textual features along with some statistical item features collected from the second-hand item shopping platform to determine whether the image and text of an uploaded second-hand item are qualified for reasonable price suggestion with a binary classification model, and provide price suggestions for second-hand items with qualified images and text descriptions with a regression model. To satisfy different demands, two different constraints are added into the joint training of the classification model and the regression model. Moreover, a customized loss function is designed for optimizing the regression model to provide price suggestions for second-hand items, which can not only maximize the gain of the sellers but also facilitate the online transaction. We also derive a set of metrics to better evaluate the proposed price suggestion system. Extensive experiments on a large real-world dataset demonstrate the effectiveness of the proposed multi-modal price suggestion system.


title: Vision-based Price Suggestion for Online Second-hand Items

abstract: Different from shopping in physical stores, where people have the opportunity to closely check a product (e.g., touching the surface of a T-shirt or smelling the scent of perfume) before making a purchase decision, online shoppers rely greatly on the uploaded product images to make any purchase decision. The decision-making is challenging when selling or purchasing second-hand items online since estimating the items' prices is not trivial. In this work, we present a vision-based price suggestion system for the online second-hand item shopping platform. The goal of vision-based price suggestion is to help sellers set effective prices for their second-hand listings with the images uploaded to the online platforms. First, we propose to better extract representative visual features from the images with the aid of some other image-based item information (e.g., category, brand). Then, we design a vision-based price suggestion module which takes the extracted visual features along with some statistical item features from the shopping platform as the inputs to determine whether an uploaded item image is qualified for price suggestion by a binary classification model, and provide price suggestions for items with qualified images by a regression model. According to two demands from the platform, two different objective functions are proposed to jointly optimize the classification model and the regression model. For better model training, we also propose a warm-up training strategy for the joint optimization. Extensive experiments on a large real-world dataset demonstrate the effectiveness of our vision-based price prediction system.


title: Simplicial 2-Complex Convolutional Neural Nets

abstract: Recently, neural network architectures have been developed to accommodate when the data has the structure of a graph or, more generally, a hypergraph. While useful, graph structures can be potentially limiting. Hypergraph structures in general do not account for higher order relations between their hyperedges. Simplicial complexes offer a middle ground, with a rich theory to draw on. We develop a convolutional neural network layer on simplicial 2-complexes.


title: A Shapley-Folkman lemma for lattice polytopes

abstract: Some occurrences of $n$ can be replaced by $n-1$ in a special case of the Shapley-Folkman lemma.


title: Efficient Distributed Transposition Of Large-Scale Multigraphs And High-Cardinality Sparse Matrices

abstract: Graph-based representations underlie a wide range of scientific problems. Graph connectivity is typically represented as a sparse matrix in the Compressed Sparse Row format. Large-scale graphs rely on distributed storage, allocating distinct subsets of rows to compute nodes. Efficient matrix transpose is an operation of high importance, providing the reverse graph pathways and a column-ordered matrix view. This operation is well studied for simple graph models. Nevertheless, its resolution for multigraphs and higher-cardinality connectivity matrices is unexistent. We advance state-of-the-art distributed transposition methods by providing a theoretical model, algorithmic details, MPI-based implementation and proof of mathematical soundness for such complex models. Benchmark results demonstrate ideal and almost ideal scaling properties for perfectly- and heterogeneously-balanced datasets, respectively


title: Is seismicity operating at a critical point?

abstract: Seismicity and faulting within the Earth crust are characterized by many scaling laws that are usually interpreted as qualifying the existence of underlying physical mechanisms associated with some kind of criticality in the sense of phase transitions. Using an augmented Epidemic-Type Aftershock Sequence (ETAS) model that accounts for the spatial variability of the background rates $\mu(x,y)$, we present a direct quantitative test of criticality. We calibrate the model to the ANSS catalog of the entire globe, the region around California, and the Geonet catalog for the region around New Zealand using an extended Expectation-Maximization (EM) algorithm including the determination of $\mu(x,y)$. We demonstrate that the criticality reported in previous studies is spurious and can be attributed to a systematic upward bias in the calibration of the branching ratio of the ETAS model, when not accounting correctly for spatial variability. We validate the version of the ETAS model which possesses a space varying background rate $\mu(x,y)$ by performing pseudo prospective forecasting tests. The non-criticality of seismicity has major implications for the prediction of large events.


title: No Krasnoselskii number for general sets in $\mathbb{R}^2$

abstract: For a family $\mathcal{F}$ of sets in $\mathbb{R}^d$, the Krasnoselskii number of $\mathcal{F}$ is the smallest $m$ such that for any $S \in \mathcal{F}$, if every $m$ points of $S$ are visible from a common point in $S$, then any finite subset of $S$ is visible from a single point. More than 35 years ago, Peterson asked whether there exists a Krasnoselskii number for general sets in $\mathbb{R}^d$. Excluding results for special cases of sets with strong topological restrictions, the best known result is due to Breen, who showed that if such a Krasnoselskii number in $\mathbb{R}^2$ exists, then it is larger than $8$. In this paper we answer Peterson's question in the negative by showing that there is no Krasnoselskii number for the family of all sets in $\mathbb{R}^2$. The proof is non-constructive, and uses transfinite induction and the well ordering theorem. In addition, we consider Krasnoselskii numbers with respect to visibility through polygonal paths of length $ \leq n$, for which an analogue of Krasnoselskii's theorem was proved by Magazanik and Perles. We show, by an explicit construction, that for any $n \geq 2$, there is no Krasnoselskii number for the family of general sets in $\mathbb{R}^2$ with respect to visibility through paths of length $\leq n$. (Here the counterexamples are finite unions of line segments.)


title: A parallel cut-cell algorithm for the free-boundary Grad-Shafranov problem

abstract: A parallel cut-cell algorithm is described to solve the free-boundary problem of the Grad-Shafranov equation. The algorithm reformulates the free-boundary problem in an irregular bounded domain and its important aspects include a searching algorithm for the magnetic axis and separatrix, a surface integral along the irregular boundary to determine the boundary values, an approach to optimize the coil current based on a targeting plasma shape, Picard iterations with Aitken's acceleration for the resulting nonlinear problem, and a Cartesian grid embedded boundary method to handle the complex geometry. The algorithm is implemented in parallel using a standard domain-decomposition approach and a good parallel scaling is observed. Numerical results verify the accuracy and efficiency of the free-boundary Grad-Shafranov solver.


title: Performance-Weighed Policy Sampling for Meta-Reinforcement Learning

abstract: This paper discusses an Enhanced Model-Agnostic Meta-Learning (E-MAML) algorithm that generates fast convergence of the policy function from a small number of training examples when applied to new learning tasks. Built on top of Model-Agnostic Meta-Learning (MAML), E-MAML maintains a set of policy parameters learned in the environment for previous tasks. We apply E-MAML to developing reinforcement learning (RL)-based online fault tolerant control schemes for dynamic systems. The enhancement is applied when a new fault occurs, to re-initialize the parameters of a new RL policy that achieves faster adaption with a small number of samples of system behavior with the new fault. This replaces the random task sampling step in MAML. Instead, it exploits the extant previously generated experiences of the controller. The enhancement is sampled to maximally span the parameter space to facilitate adaption to the new fault. We demonstrate the performance of our approach combining E-MAML with proximal policy optimization (PPO) on the well-known cart pole example, and then on the fuel transfer system of an aircraft.


title: On-Request Wireless Charging and Partial Computation Offloading In Multi-Access Edge Computing Systems

abstract: Wireless charging coupled with computation offloading in edge networks offers a promising solution for realizing power-hungry and computation intensive applications on user devices. We consider a multi-access edge computing (MEC) system with collocated MEC server and base-station/access point, each equipped with a massive MIMO antenna array, supporting multiple users requesting data computation and wireless charging. The goal is to minimize the energy consumption for computation offloading and maximize the received energy at the user from wireless charging. The proposed solution is a novel two-stage algorithm employing nested primal-dual and linear programming techniques to perform data partitioning and time allocation for computation offloading and design the optimal energy beamforming for wireless charging, all within MEC-AP transmit power and latency constraints. Algorithm results show that optimal energy beamforming significantly outperforms other schemes such as isotropic or directed charging without beam power allocation. Compared to binary offloading, data partition in partial offloading leads to lower energy consumption and more charging time, and hence offers better wireless charging performance. The charged energy over an extended period of time both with and without computation offloading can be substantial. Opportunistic wireless charging from MEC-AP thus offers a viable untethered approach for supplying energy to user-devices.


title: A MAC-less Neural Inference Processor Supporting Compressed, Variable Precision Weights

abstract: This paper introduces two architectures for the inference of convolutional neural networks (CNNs). Both architectures exploit weight sparsity and compression to reduce computational complexity and bandwidth. The first architecture uses multiply-accumulators (MACs) but avoids unnecessary multiplications by skipping zero weights. The second architecture exploits weight sparsity at the level of their bit representation by substituting resource-intensive MACs with much smaller Bit Layer Multiply Accumulators (BLMACs). The use of BLMACs also allows variable precision weights as variable size integers and even floating points. Some details of an implementation of the second architecture are given. Weight compression with arithmetic coding is also discussed as well as bandwidth implications. Finally, some implementation results for a pathfinder design and various technologies are presented.