title: Active Magnetoplasmonics with Transparent Conductive Oxide Nanocrystals

abstract: Magnetoplasmonics is highly promising to devise active optical elements: modulating the plasmon resonance condition with magnetic field can boost the performance of refractometric sensors and nanophotonic optical devices. Nevertheless, real life applications are hampered by the magnetoplasmonic trilemma: 1) a good plasmonic metal has sharp optical resonances but low magneto-optical response; 2) a magnetic metal has strong magneto-optical response but a very broad plasmonic resonance; 3) mixing the two components degrades the quality of both features. To overcome the trilemma, we use a different class of materials, transparent conductive oxide nanocrystals (NCs) with plasmonic response in the near infrared. Although non-magnetic, they combine a large cyclotron frequency (due to small electron effective mass) with sharp plasmonic resonances. We benchmark the concept with F- and In- doped CdO (FICO) and Sn-doped In2O3 (ITO) NCs to boost the magneto-optical Faraday rotation and ellipticity, reaching the highest magneto-optical response for a non-magnetic plasmonic material, and exceeding the performance of state-of-the-art ferromagnetic nanostructures. The magnetoplasmonic response of these NCs was rationalized with analytical model based on the excitation of circular magnetoplasmonic modes. Finally, proof of concept experiments demonstrated the superior performance of FICO NCs with respect to current state of the art in magnetoplasmonic refractometric sensing, approaching the sensitivity of leading localized plasmon refractometric methods with the advantage of not requiring complex curve fitting.


title: Jointly Modeling and Clustering Tensors in High Dimensions

abstract: We consider the problem of jointly modeling and clustering populations of tensors by introducing a high-dimensional tensor mixture model with heterogeneous covariances. To effectively tackle the high dimensionality of tensor objects, we employ plausible dimension reduction assumptions that exploit the intrinsic structures of tensors such as low-rankness in the mean and separability in the covariance. In estimation, we develop an efficient high-dimensional expectation-conditional-maximization (HECM) algorithm that breaks the intractable optimization in the M-step into a sequence of much simpler conditional optimization problems, each of which is convex, admits regularization and has closed-form updating formulas. Our theoretical analysis is challenged by both the non-convexity in the EM-type estimation and having access to only the solutions of conditional maximizations in the M-step, leading to the notion of dual non-convexity. We demonstrate that the proposed HECM algorithm, with an appropriate initialization, converges geometrically to a neighborhood that is within statistical precision of the true parameter. The efficacy of our proposed method is demonstrated through comparative numerical experiments and an application to a medical study, where our proposal achieves an improved clustering accuracy over existing benchmarking methods.


title: Algebroids, AKSZ Constructions and Doubled Geometry

abstract: We give a self-contained survey of some approaches aimed at a global description of the geometry underlying double field theory. After reviewing the geometry of Courant algebroids and their incarnations in the AKSZ construction, we develop the theory of metric algebroids including their graded geometry. We use metric algebroids to give a global description of doubled geometry, incorporating the section constraint, as well as an AKSZ-type construction of topological doubled sigma-models. When these notions are combined with ingredients of para-Hermitian geometry, we demonstrate how they reproduce kinematical features of double field theory from a global perspective, including solutions of the section constraint for Riemannian foliated doubled manifolds, as well as a natural notion of generalized T-duality for polarized doubled manifolds. We describe the $L_\infty$-algebras of symmetries of a doubled geometry, and briefly discuss other proposals for global doubled geometry in the literature.


title: Fermi liquids and fractional statistics in one dimension

abstract: Interacting fermion systems in one dimension, which in the low energy approximation are described by Luttinger liquid theory, can be reformulated as systems of weakly interacting particles with fractional exchange statistics. This is shown by use of Landau's Fermi liquid theory, with quasiparticles interpreted as adiabatically dressed fermions. An application of this method is included, where boundary excitations of a two-dimensional quantum Hall electron system are studied.


title: Demystifying Memory Access Patterns of FPGA-Based Graph Processing Accelerators

abstract: Recent advances in reprogrammable hardware (e.g., FPGAs) and memory technology (e.g., DDR4, HBM) promise to solve performance problems inherent to graph processing like irregular memory access patterns on traditional hardware (e.g., CPU). While several of these graph accelerators were proposed in recent years, it remains difficult to assess their performance and compare them on common graph workloads and accelerator platforms, due to few open source implementations and excessive implementation effort. In this work, we build on a simulation environment for graph processing accelerators, to make several existing accelerator approaches comparable. This allows us to study relevant performance dimensions such as partitioning schemes and memory technology, among others. The evaluation yields insights into the strengths and weaknesses of current graph processing accelerators along these dimensions, and features a novel in-depth comparison.


title: Proteno: Text Normalization with Limited Data for Fast Deployment in Text to Speech Systems

abstract: Developing Text Normalization (TN) systems for Text-to-Speech (TTS) on new languages is hard. We propose a novel architecture to facilitate it for multiple languages while using data less than 3% of the size of the data used by the state of the art results on English. We treat TN as a sequence classification problem and propose a granular tokenization mechanism that enables the system to learn majority of the classes and their normalizations from the training data itself. This is further combined with minimal precoded linguistic knowledge for other classes. We publish the first results on TN for TTS in Spanish and Tamil and also demonstrate that the performance of the approach is comparable with the previous work done on English. All annotated datasets used for experimentation will be released at https://github.com/amazon-research/proteno.


title: Recent Advances in Domain Adaptation for the Classification of Remote Sensing Data

abstract: The success of supervised classification of remotely sensed images acquired over large geographical areas or at short time intervals strongly depends on the representativity of the samples used to train the classification algorithm and to define the model. When training samples are collected from an image (or a spatial region) different from the one used for mapping, spectral shifts between the two distributions are likely to make the model fail. Such shifts are generally due to differences in acquisition and atmospheric conditions or to changes in the nature of the object observed. In order to design classification methods that are robust to data-set shifts, recent remote sensing literature has considered solutions based on domain adaptation (DA) approaches. Inspired by machine learning literature, several DA methods have been proposed to solve specific problems in remote sensing data classification. This paper provides a critical review of the recent advances in DA for remote sensing and presents an overview of methods divided into four categories: i) invariant feature selection; ii) representation matching; iii) adaptation of classifiers and iv) selective sampling. We provide an overview of recent methodologies, as well as examples of application of the considered techniques to real remote sensing images characterized by very high spatial and spectral resolution. Finally, we propose guidelines to the selection of the method to use in real application scenarios.


title: An ALMA Gas-dynamical Mass Measurement of the Supermassive Black Hole in the Local Compact Galaxy UGC 2698

abstract: We present 0\farcs{14}-resolution Atacama Large Millimeter/submillimeter Array (ALMA) CO(2$-$1) observations of the circumnuclear gas disk in UGC 2698, a local compact galaxy. The disk exhibits regular rotation with projected velocities rising to 450 km s$^{-1}$ near the galaxy center. We fit gas-dynamical models to the ALMA data cube, assuming the CO emission originates from a dynamically cold, thin disk, and measured the mass of the supermassive black hole (BH) in UGC 2698 to be $M_{\mathrm{BH}} = (2.46 \pm{0.07}$ [$1\sigma$ stat] $^{+0.70}_{-0.78}$ [sys])$\times 10^9$ $M_\odot$. UGC 2698 is part of a sample of nearby early-type galaxies that are plausible $z\sim2$ red nugget relics. Previous stellar-dynamical modeling for three galaxies in the sample found BH masses consistent with the BH mass$-$stellar velocity dispersion ($M_{\mathrm{BH}}-\sigma_\star$) relation but over-massive relative to the BH mass$-$bulge luminosity ($M_{\mathrm{BH}}-L_{\mathrm{bul}}$) correlation, suggesting that BHs may gain the majority of their mass before their host galaxies. However, UGC 2698 is consistent with both $M_{\mathrm{BH}}-\sigma_\star$ and $M_{\mathrm{BH}}-L_{\mathrm{bul}}$. As UGC 2698 has the largest stellar mass and effective radius in the local compact galaxy sample, it may have undergone more recent mergers that brought it in line with the BH scaling relations. Alternatively, given that the three previously-measured compact galaxies are outliers from $M_{\mathrm{BH}}-L_{\mathrm{bul}}$, while UGC 2698 is not, there may be significant scatter at the poorly sampled high-mass end of the relation. Additional gas-dynamical $M_{\mathrm{BH}}$ measurements for the compact galaxy sample will improve our understanding of BH$-$galaxy co-evolution.


title: Impact of gender on the formation and outcome of mentoring relationships in academic research

abstract: Despite increasing representation in graduate training programs, a disproportionate number of women leave academic research before obtaining an independent position. To understand factors underlying this trend, we analyzed a multidisciplinary database of Ph.D. and postdoctoral mentoring relationships covering the years 2000-2020, focusing on data from the life sciences. Student and mentor gender are both associated with differences in rates of student's continuation to independent mentor positions of their own. Although trainees of women mentors are less likely to take on independent positions than trainees of men mentors, this effect is reduced substantially after controlling for several measurements of mentor status. Thus the effect of mentor gender can be explained at least partially by gender disparities in social and financial resources available to mentors. Because trainees and mentors tend to be of the same gender, this association between mentor gender and academic continuation disproportionately impacts women trainees. On average, gender homophily in graduate training is unrelated to mentor status. A notable exception to this trend is the special case of scientists having been granted an outstanding distinction, evidenced by membership in the National Academy of Sciences, being a grantee of the Howard Hughes Medical Institute, or having been awarded the Nobel Prize. This group of mentors trains men graduate students at higher rates than their most successful colleagues. These results suggest that, in addition to other factors that limit career choices for women trainees, gender inequities in mentors' access to resources and prestige contribute to women's attrition from independent research positions.


title: Convergence Rates of Distributed Consensus over Cluster Networks: A Two-Time-Scale Approach

abstract: We study the popular distributed consensus method over networks composed of a number of densely connected clusters with a sparse connection between them. In these cluster networks, the method often constitutes two-time-scale dynamics, where the internal nodes within each cluster reach consensus quickly relative to the aggregate nodes across clusters. Our main contribution is to provide the rate of the distributed consensus method, which characterize explicitly the impacts of the internal and external graphs on the performance of this method. Our main result shows that this rate converges exponentially and only scales with a few number of nodes, which is relatively small to the size of the network. The key technique in our analysis is to consider a Lyapunov function which captures the impacts of different time-scale dynamics on the convergence of the method. Our approach avoids using model reduction, which is the typical way according to singular perturbation theory and relies on relatively simple definitions of the slow and fast variables. In addition, Lyapunov analysis allows us to derive the rate of distributed consensus methods over cluster networks, which is missing from the existing works using singular perturbation theory. We illustrate our theoretical results by a number of numerical simulations over different cluster networks.


title: Sublanguage: A Serious Issue Affects Pretrained Models in Legal Domain

abstract: Legal English is a sublanguage that is important for everyone but not for everyone to understand. Pretrained models have become best practices among current deep learning approaches for different problems. It would be a waste or even a danger if these models were applied in practice without knowledge of the sublanguage of the law. In this paper, we raise the issue and propose a trivial solution by introducing BERTLaw a legal sublanguage pretrained model. The paper's experiments demonstrate the superior effectiveness of the method compared to the baseline pretrained model


title: Poking Holes: Looking for Gaps in LIGO/Virgo's Black Hole Population

abstract: Stellar evolution models predict the existence of a gap in the black hole mass spectrum from $\sim55 M_\odot - 120 M_\odot$ due to pair-instability supernovae (PISNe). We investigate the possible existence of such an "upper" mass gap in the second gravitational wave transient catalog (GWTC-2) by hierarchically modeling the astrophysical distribution of black hole masses. We extend the Truncated and Powerlaw+Peak mass distribution families to allow for an explicit gap in the mass distribution, and apply the extended models to GWTC-2. We find that with the Truncated model there is mild evidence favoring an upper mass gap with log Bayes Factor $\ln \mathcal{B} = 2.79$, inferring the lower and upper bounds at $56.12_{-4.38}^{+7.54} M_\odot$, and $103.74_{-6.32}^{+17.01} M_\odot$ respectively. When using the Powerlaw+Peak model, we find no preference for the gap. When imposing tighter priors on the gap bounds centered on the expected PISN gap bounds, the log Bayes factors in favor of a gap mildly increase. These results are however contingent on the parameter inference for the most massive binary, GW190521, for which follow up analyses showed the source may be an intermediate mass ratio merger that has component masses straddling the gap. Using the GW190521 posterior samples from the analysis in Nitz & Capano (2021), we find an increase in Bayes factors in favor of the gap. However, the overall conclusions are unchanged: There is no preference for a gap when using the Powerlaw+Peak model. This work paves the way for constraining the physics of pair-instability and pulsational pair-instability supernovae and high-mass black hole formation.


title: A survey of active learning algorithms for supervised remote sensing image classification

abstract: Defining an efficient training set is one of the most delicate phases for the success of remote sensing image classification routines. The complexity of the problem, the limited temporal and financial resources, as well as the high intraclass variance can make an algorithm fail if it is trained with a suboptimal dataset. Active learning aims at building efficient training sets by iteratively improving the model performance through sampling. A user-defined heuristic ranks the unlabeled pixels according to a function of the uncertainty of their class membership and then the user is asked to provide labels for the most uncertain pixels. This paper reviews and tests the main families of active learning algorithms: committee, large margin and posterior probability-based. For each of them, the most recent advances in the remote sensing community are discussed and some heuristics are detailed and tested. Several challenging remote sensing scenarios are considered, including very high spatial resolution and hyperspectral image classification. Finally, guidelines for choosing the good architecture are provided for new and/or unexperienced user.


title: Ridge Regression Neural Network for Pediatric Bone Age Assessment

abstract: Bone age is an important measure for assessing the skeletal and biological maturity of children. Delayed or increased bone age is a serious concern for pediatricians, and needs to be accurately assessed in a bid to determine whether bone maturity is occurring at a rate consistent with chronological age. In this paper, we introduce a unified deep learning framework for bone age assessment using instance segmentation and ridge regression. The proposed approach consists of two integrated stages. In the first stage, we employ an image annotation and segmentation model to annotate and segment the hand from the radiographic image, followed by background removal. In the second stage, we design a regression neural network architecture composed of a pre-trained convolutional neural network for learning salient features from the segmented pediatric hand radiographs and a ridge regression output layer for predicting the bone age. Experimental evaluation on a dataset of hand radiographs demonstrates the competitive performance of our approach in comparison with existing deep learning based methods for bone age assessment.


title: Test beam characterization of sensor prototypes for the CMS Barrel MIP Timing Detector

abstract: The MIP Timing Detector will provide additional timing capabilities for detection of minimum ionizing particles (MIPs) at CMS during the High Luminosity LHC era, improving event reconstruction and pileup rejection. The central portion of the detector, the Barrel Timing Layer (BTL), will be instrumented with LYSO:Ce crystals and Silicon Photomultipliers (SiPMs) providing a time resolution of about 30 ps at the beginning of operation, and degrading to 50-60 ps at the end of the detector lifetime as a result of radiation damage. In this work, we present the results obtained using a 120 GeV proton beam at the Fermilab Test Beam Facility to measure the time resolution of unirradiated sensors. A proof-of-concept of the sensor layout proposed for the barrel region of the MTD, consisting of elongated crystal bars with dimensions of about 3 x 3 x 57 mm$^3$ and with double-ended SiPM readout, is demonstrated. This design provides a robust time measurement independent of the impact point of the MIP along the crystal bar. We tested LYSO:Ce bars of different thickness (2, 3, 4 mm) with a geometry close to the reference design and coupled to SiPMs manufactured by Hamamatsu and Fondazione Bruno Kessler. The various aspects influencing the timing performance such as the crystal thickness, properties of the SiPMs (e.g. photon detection efficiency), and impact angle of the MIP are studied. A time resolution of about 28 ps is measured for MIPs crossing a 3 mm thick crystal bar, corresponding to an MPV energy deposition of 2.6 MeV, and of 22 ps for the 4.2 MeV MPV energy deposition expected in the BTL, matching the detector performance target for unirradiated devices.


title: Rethinking Text Line Recognition Models

abstract: In this paper, we study the problem of text line recognition. Unlike most approaches targeting specific domains such as scene-text or handwritten documents, we investigate the general problem of developing a universal architecture that can extract text from any image, regardless of source or input modality. We consider two decoder families (Connectionist Temporal Classification and Transformer) and three encoder modules (Bidirectional LSTMs, Self-Attention, and GRCLs), and conduct extensive experiments to compare their accuracy and performance on widely used public datasets of scene and handwritten text. We find that a combination that so far has received little attention in the literature, namely a Self-Attention encoder coupled with the CTC decoder, when compounded with an external language model and trained on both public and internal data, outperforms all the others in accuracy and computational complexity. Unlike the more common Transformer-based models, this architecture can handle inputs of arbitrary length, a requirement for universal line recognition. Using an internal dataset collected from multiple sources, we also expose the limitations of current public datasets in evaluating the accuracy of line recognizers, as the relatively narrow image width and sequence length distributions do not allow to observe the quality degradation of the Transformer approach when applied to the transcription of long lines.


title: PyTorch Geometric Temporal: Spatiotemporal Signal Processing with Neural Machine Learning Models

abstract: We present PyTorch Geometric Temporal a deep learning framework combining state-of-the-art machine learning algorithms for neural spatiotemporal signal processing. The main goal of the library is to make temporal geometric deep learning available for researchers and machine learning practitioners in a unified easy-to-use framework. PyTorch Geometric Temporal was created with foundations on existing libraries in the PyTorch eco-system, streamlined neural network layer definitions, temporal snapshot generators for batching, and integrated benchmark datasets. These features are illustrated with a tutorial-like case study. Experiments demonstrate the predictive performance of the models implemented in the library on real world problems such as epidemiological forecasting, ridehail demand prediction and web-traffic management. Our sensitivity analysis of runtime shows that the framework can potentially operate on web-scale datasets with rich temporal features and spatial structure.


title: Detect and Classify -- Joint Span Detection and Classification for Health Outcomes

abstract: A health outcome is a measurement or an observation used to capture and assess the effect of a treatment. Automatic detection of health outcomes from text would undoubtedly speed up access to evidence necessary in healthcare decision making. Prior work on outcome detection has modelled this task as either (a) a sequence labelling task, where the goal is to detect which text spans describe health outcomes, or (b) a classification task, where the goal is to classify a text into a pre-defined set of categories depending on an outcome that is mentioned somewhere in that text. However, this decoupling of span detection and classification is problematic from a modelling perspective and ignores global structural correspondences between sentence-level and word-level information present in a given text. To address this, we propose a method that uses both word-level and sentence-level information to simultaneously perform outcome span detection and outcome type classification. In addition to injecting contextual information to hidden vectors, we use label attention to appropriately weight both word and sentence level information. Experimental results on several benchmark datasets for health outcome detection show that our proposed method consistently outperforms decoupled methods, reporting competitive results.


title: Primary Velocity and Orbital Phase Effects on Planetary Detectability from Small Epoch Number Data Sets

abstract: Cross correlation analyses of high resolution spectroscopic data have recently shown great success in directly detecting planetary signals and enabling the characterization of their atmospheres. One such technique aims to observe a system at multiple epochs and combine the measured planetary radial velocities from each epoch into a measurement of the planetary Keplerian orbital velocity $K_p$, constituting a direct detection of the planetary signal. Recent work has shown that in few-epoch ($\sim$5) data sets, unintended structure can arise at a high level, obscuring the planetary detection. In this work, we look to simulations to examine whether there are ways to reduce this structured noise in few-epoch data sets by careful planning of observations. The choice of observation date allows observers to select the primary (stellar) velocity - through a set systemic velocity and chosen barycentric velocity - and the planetary orbital phase, and so we focus on the effects of these two parameters. We find that epochs taken when the primary velocity is near zero, and the stellar lines remain relatively fixed to the telluric rest-frame, greatly reduce the level of structured noise and allow for much stronger planetary detections, on average more than twice the significance of detections made with epochs using randomly selected primary velocities. Following these results, we recommend that observers looking to build up high-resolution multi-epoch data sets target nights when their system has a near-zero primary velocity.


title: Learning User's confidence for active learning

abstract: In this paper, we study the applicability of active learning in operative scenarios: more particularly, we consider the well-known contradiction between the active learning heuristics, which rank the pixels according to their uncertainty, and the user's confidence in labeling, which is related to both the homogeneity of the pixel context and user's knowledge of the scene. We propose a filtering scheme based on a classifier that learns the confidence of the user in labeling, thus minimizing the queries where the user would not be able to provide a class for the pixel. The capacity of a model to learn the user's confidence is studied in detail, also showing the effect of resolution is such a learning task. Experiments on two QuickBird images of different resolutions (with and without pansharpening) and considering committees of users prove the efficiency of the filtering scheme proposed, which maximizes the number of useful queries with respect to traditional active learning.


title: Geometry encoding for numerical simulations

abstract: We present a notion of geometry encoding suitable for machine learning-based numerical simulation. In particular, we delineate how this notion of encoding is different than other encoding algorithms commonly used in other disciplines such as computer vision and computer graphics. We also present a model comprised of multiple neural networks including a processor, a compressor and an evaluator.These parts each satisfy a particular requirement of our encoding. We compare our encoding model with the analogous models in the literature


title: Modelling epidemics: a perspective on mathematical models and their use in the present SARS-CoV-2 epidemic

abstract: In this work we look at several mathematical models that have been constructed during the present pandemic to address dfferent issues of importance to public health policies about epidemic scenarios and thier causes. We start by briefly reviewing the most basic properties of the classic Kermack-McKendrick models and then proceed to look at some generalizations and their applications to contact structures, cocirculation of viral infections, growth patterns of epidemic curves, characterization of probability distributions and passage times necessary to parametrize the compartments that are added to the basic Kermack-McKendrick model. All of these examples revolve around the idea that a system of differential equations constructed for an specifc epidemiological problem, has as central and main theoretical and conceptual support the epidemiological, medical and biological context that motivates its construction and analysis. Fitting differential models to data without a sound perspective on the biological feasibility of the mathematical model will provide an statistical significant fit but will also fall short of providing useful information for the management, mitigation and eventual control of the epidemic of SARS-CoV-2.


title: An $L^2$ Analysis of Reinforcement Learning in High Dimensions with Kernel and Neural Network Approximation

abstract: Reinforcement learning (RL) algorithms based on high-dimensional function approximation have achieved tremendous empirical success in large-scale problems with an enormous number of states. However, most analysis of such algorithms gives rise to error bounds that involve either the number of states or the number of features. This paper considers the situation where the function approximation is made either using the kernel method or the two-layer neural network model, in the context of a fitted Q-iteration algorithm with explicit regularization. We establish an $\tilde{O}(H^3|\mathcal {A}|^{\frac14}n^{-\frac14})$ bound for the optimal policy with $Hn$ samples, where $H$ is the length of each episode and $|\mathcal {A}|$ is the size of action space. Our analysis hinges on analyzing the $L^2$ error of the approximated Q-function using $n$ data points. Even though this result still requires a finite-sized action space, the error bound is independent of the dimensionality of the state space.


title: SALT3: An Improved Type Ia Supernova Model for Measuring Cosmic Distances

abstract: A spectral-energy distribution (SED) model for Type Ia supernovae (SNe Ia) is a critical tool for measuring precise and accurate distances across a large redshift range and constraining cosmological parameters. We present an improved model framework, SALT3, which has several advantages over current models including the leading SALT2 model (SALT2.4). While SALT3 has a similar philosophy, it differs from SALT2 by having improved estimation of uncertainties, better separation of color and light-curve stretch, and a publicly available training code. We present the application of our training method on a cross-calibrated compilation of 1083 SNe with 1207 spectra. Our compilation is $2.5\times$ larger than the SALT2 training sample and has greatly reduced calibration uncertainties. The resulting trained SALT3.K21 model has an extended wavelength range $2000$-$11000$ angstroms (1800 angstroms redder) and reduced uncertainties compared to SALT2, enabling accurate use of low-$z$ $I$ and $iz$ photometric bands. Including these previously discarded bands, SALT3.K21 reduces the Hubble scatter of the low-z Foundation and CfA3 samples by 15% and 10%, respectively. To check for potential systematic uncertainties we compare distances of low ($0.01<z<0.2$) and high ($0.4<z<0.6$) redshift SNe in the training compilation, finding an insignificant $2\pm14$ mmag shift between SALT2.4 and SALT3.K21. While the SALT3.K21 model was trained on optical data, our method can be used to build a model for rest-frame NIR samples from the Roman Space Telescope. Our open-source training code, public training data, model, and documentation are available at https://saltshaker.readthedocs.io/en/latest/, and the model is integrated into the sncosmo and SNANA software packages.


title: Inexact-Proximal Accelerated Gradient Method for Stochastic Nonconvex Constrained Optimization Problems

abstract: Stochastic nonconvex optimization problems with nonlinear constraints have a broad range of applications in intelligent transportation, cyber-security, and smart grids. In this paper, first, we propose an inexact-proximal accelerated gradient method to solve a nonconvex stochastic composite optimization problem where the objective is the sum of smooth and nonsmooth functions, the constraint functions are assumed to be deterministic and the solution to the proximal map of the nonsmooth part is calculated inexactly at each iteration. We demonstrate an asymptotic sublinear rate of convergence for stochastic settings using increasing sample-size considering the error in the proximal operator diminishes at an appropriate rate. Then we customize the proposed method for solving stochastic nonconvex optimization problems with nonlinear constraints and demonstrate a convergence rate guarantee. Numerical results show the effectiveness of the proposed algorithm.


title: Point-to-Point Coulomb Effects in High Brightness Photoelectron Beamlines for Ultrafast Electron Diffraction

abstract: In an effort to increase spatial and temporal resolution of ultrafast electron diffraction and microscopy, ultra-high brightness photocathodes are actively sought to improve electron beam quality. Beam dynamics codes often approximate the Coulomb interaction with mean-field space charge, which is a good approximation in traditional beams. However, point-to-point Coulomb effects, such as disorder induced heating and the Boersch effect, can not be neglected in cold, dense beams produced by such photocathodes. In this paper, we introduce two new numerical methods to calculate the important effects of the photocathode image charge when using a point-to-point interaction model. Equipped with an accurate model of the image charge, we calculate the effects of point-to-point interactions on two high brightness photoemission beamlines for ultrafast diffraction. The first beamline uses a 200 keV gun, whereas the second uses a 5 MeV gun, each operating in the single-shot diffraction regime with $10^5$ electrons/pulse. Assuming a zero photoemission temperature, it is shown that including stochastic Coulomb effects increases the final emittance of these beamlines by over a factor of 2 and decreases the peak transverse phase space density by over a factor of 3 as compared to mean-field simulations. We then introduce a method to compute the energy released by disorder induced heating using the pair correlation function. This disorder induced heating energy was found to scale very near the theoretical result for stationary ultracold plasmas, and it accounts for over half of the emittance growth above mean-field simulations.


title: Memory Order Decomposition of Symbolic Sequences

abstract: We introduce a general method for the study of memory in symbolic sequences based on higher-order Markov analysis. The Markov process that best represents a sequence is expressed as a mixture of matrices of minimal orders, enabling the definition of the so-called memory profile, which unambiguously reflects the true order of correlations. The method is validated by recovering the memory profiles of tunable synthetic sequences. Finally, we scan real data and showcase with practical examples how our protocol can be used to extract relevant stochastic properties of symbolic sequences.

categories: physics.soc-ph cs.IT math.IT physics.data-an


title: Robot to support older people to live independently

abstract: This paper presents an overview on how the PAL Robotics ARI robot is participating in the European SHAPES project to promote healthy and active living among older people, by integrating digital solutions from project partners and adapting the system in order to improve human-robot interaction and user acceptability in a wide range of tasks.


title: Towards Robust Neural Retrieval Models with Synthetic Pre-Training

abstract: Recent work has shown that commonly available machine reading comprehension (MRC) datasets can be used to train high-performance neural information retrieval (IR) systems. However, the evaluation of neural IR has so far been limited to standard supervised learning settings, where they have outperformed traditional term matching baselines. We conduct in-domain and out-of-domain evaluations of neural IR, and seek to improve its robustness across different scenarios, including zero-shot settings. We show that synthetic training examples generated using a sequence-to-sequence generator can be effective towards this goal: in our experiments, pre-training with synthetic examples improves retrieval performance in both in-domain and out-of-domain evaluation on five different test sets.


title: Diffusion in the Anderson model in higher dimensions

abstract: We present an extended microcanonical Lanczos method (MCLM) for a direct evaluation of the diffusion constant and its frequency dependence within the disordered Anderson model of noninteracting particles. The method allows to study systems beyond $10^6$ sites of hypercubic lattices in $ d = 3- 7$ dimensions. Below the transition to localization, where we confirm dynamical scaling behavour, of interest is a wide region of incoherent diffusion, similar to percolating phenomena and to interacting many-body localized systems.


title: Constraining Ultralight Axions with Galaxy Surveys

abstract: Ultralight axions and other bosons are dark matter candidates present in many high energy physics theories beyond the Standard Model. In particular, the string axiverse postulates the existence of up to $\mathcal{O}(100)$ light scalar bosons constituting the dark sector. Considering a mixture of axions and cold dark matter, we obtain upper bounds for the axion relic density $\Omega_a h^2 < 0.004$ for axions of mass $10^{-31}\;\mathrm{eV}\leq m_a \leq 10^{-26}\;\mathrm{eV}$ at 95% confidence. We also improve existing constraints by a factor of over 4.5 and 2.1 for axion masses of $10^{-25}$ eV and $10^{-32}$ eV, respectively. We use the Fourier-space galaxy clustering statistics from the Baryon Oscillation Spectroscopic Survey (BOSS) and demonstrate how galaxy surveys break important degeneracies in the axion parameter space compared to the cosmic microwave background (CMB). We test the validity of the effective field theory of large-scale structure approach to mixed ultralight axion dark matter by making our own mock galaxy catalogs and find an anisotropic ultralight axion signature in the galaxy quadrupole. We also observe an enhancement of the linear galaxy bias from 1.8 to 2.4 when allowing for 5% of the dark matter to be composed of a $10^{-28}$ eV axion in our simulations. Finally, we develop an augmented interpolation scheme allowing a fast computation of the axion contribution to the linear matter power spectrum leading to a 70% reduction of the computational cost for the full Monte Carlo Markov chains analysis.


title: Semisupervised Manifold Alignment of Multimodal Remote Sensing Images

abstract: We introduce a method for manifold alignment of different modalities (or domains) of remote sensing images. The problem is recurrent when a set of multitemporal, multisource, multisensor and multiangular images is available. In these situations, images should ideally be spatially coregistred, corrected and compensated for differences in the image domains. Such procedures require the interaction of the user, involve tuning of many parameters and heuristics, and are usually applied separately. Changes of sensors and acquisition conditions translate into shifts, twists, warps and foldings of the image distributions (or manifolds). The proposed semisupervised manifold alignment (SS-MA) method aligns the images working directly on their manifolds, and is thus not restricted to images of the same resolutions, either spectral or spatial. SS-MA pulls close together samples of the same class while pushing those of different classes apart. At the same time, it preserves the geometry of each manifold along the transformation. The method builds a linear invertible transformation to a latent space where all images are alike, and reduces to solving a generalized eigenproblem of moderate size. We study the performance of SS-MA in toy examples and in real multiangular, multitemporal, and multisource image classification problems. The method performs well for strong deformations and leads to accurate classification for all domains.


title: Defect Engineered Room-Temperature Ferromagnetism in Quasi-Two-Dimensional Nitrided CoTa2O6

abstract: Thermal ammonolysis of quasi-two-dimensional (quasi-2D) CoTa2O6 yields the O2-/N3- and anionic vacancy ordered Co2+Ta5+2O6-xN2x/3$\Box$x/3 (x $\leq$ 0.15) that exhibits a transition from antiferromagnetism to defect engineered above room-temperature ferromagnetism as evidenced by diffraction, spectroscopic and magnetic characterizations. First-principles calculations reveal the origin of ferromagnetism is a particular CoON configuration with N located at Wyckoff position 8j, which breaks mirror symmetry about ab plane. A pressure-induced electronic phase transition is also predicted at around 24.5 GPa, accompanied by insulator-to-metal transition and magnetic moment vanishing.


title: The fast non-ferric kicker system for the Muon $g-2$ Experiment at Fermilab

abstract: We describe the installation, commissioning, and characterization of the new injection kicker system in the Muon $g-2$ Experiment (E989) at Fermilab, which makes a precision measurement of the muon magnetic anomaly. Three Blumlein pulsers drive each of the 1.27-m-long non-ferric kicker magnets, which reside in a storage ring vacuum (SRV) that is subjected to a 1.45 T magnetic field. The new system has been redesigned relative to Muon $g-2$'s predecessor experiment, and we present those details in this manuscript.


title: The SIR-P Model: An Illustration of the Screening Paradox

abstract: In previous work by this author, the screening paradox - the loss of predictive power of screening tests over time $t$ - was mathematically formalized using Bayesian theory. Where $J$ is Youden's statistic, $b$ is the specificity of the screening test and $\phi$ is the prevalence of disease, the ratio of positive predictive values at subsequent time $k$, $\rho(\phi_{k})$, over the original $\rho(\phi_{0})$ at $t_0$ is given by: $\zeta(\phi_{0},k) = \frac{\rho(\phi_{k})}{\rho(\phi_{0})} =\frac{\phi_k(1-b)+J\phi_0\phi_k}{\phi_0(1-b)+J\phi_0\phi_k}$ Herein, we modify the traditional Kermack-McKendrick SIR Model to include the fluctuation of the positive predictive value $\rho(\phi)$ (PPV) of a screening test over time as a function of the prevalence threshold $\phi_e$. We term this modified model the SIR-P model. Where a = sensitivity, b = specificity, $S$ = number susceptible, $I$ = number infected, $R$ = number recovered/dead, $\beta$ = infectious rate, $\gamma$ = recovery rate, and $N$ is the total number in the population, the predictive value $\rho(\phi,t)$ over time $t$ is given by: $\rho(\phi,t) = \frac{a[\frac{\beta IS}{N}-\gamma I]}{ a[\frac{\beta IS}{N}-\gamma I]+(1-b)(1-[\frac{\beta IS}{N}-\gamma I])}$ Otherwise stated: $\rho(\phi,t) = \frac{a\frac{dI}{dt}}{ a\frac{dI}{dt}+(1-b)(1-\frac{dI}{dt})}$ where $\frac{dI}{dt}$ is the fluctuation of infected individuals over time $t$.


title: From Personal Data to Digital Legacy: Exploring Conflicts in the Sharing, Security and Privacy of Post-mortem Data

abstract: As digital technologies become more prevalent there is a growing awareness of the importance of good security and privacy practices. The tools and techniques used to achieve this are typically designed with the living user in mind, with little consideration of how they should or will perform after the user has died. We report on two workshops carried out with users of password managers to explore their views on the post-mortem sharing, security and privacy of a range of common digital assets. We discuss a post-mortem privacy paradox where users recognise value in planning for their digital legacy, yet avoid actively doing so. Importantly, our findings highlight a tension between the use of recommended security tools during life and facilitating appropriate post-mortem access to chosen assets. We offer design recommendations to facilitate and encourage digital legacy planning while promoting good security habits during life.


title: The X17 boson and the $^3$H$(p,e^+ e^-)^4$He and $^3$He$(n,e^+ e^-)^4$He processes: a theoretical analysis

abstract: The present work deals with $e^+$-$e^-$ pair production in the four-nucleon system. We first analyze the process as a purely electromagnetic one in the context of a state-of-the-art approach to nuclear strong-interaction dynamics and nuclear electromagnetic currents, derived from chiral effective field theory ($\chi$EFT). Next, we examine how the exchange of a hypothetical low-mass boson would impact the cross section for such a process. We consider several possibilities, that this boson is either a scalar, pseudoscalar, vector, or axial particle. The ab initio calculations use exact hyperspherical-harmonics methods to describe the bound state and low-energy spectrum of the $A\,$=$\,4$ continuum, and fully account for initial state interaction effects in the $3+1$ clusters. While electromagnetic interactions are treated to high orders in the chiral expansion, the interactions of the hypothetical boson with nucleons are modeled in leading-order $\chi$EFT (albeit, in some instances, selected subleading contributions are also accounted for). We also provide an overview of possible future experiments probing pair production in the $A\,$=$\,4$ system at a number of candidate facilities.facilities.


title: A Novel Hybrid Deep Learning Approach for Non-Intrusive Load Monitoring of Residential Appliance Based on Long Short Term Memory and Convolutional Neural Networks

abstract: Energy disaggregation or nonintrusive load monitoring (NILM), is a single-input blind source discrimination problem, aims to interpret the mains user electricity consumption into appliance level measurement. This article presents a new approach for power disaggregation by using a deep recurrent long short term memory (LSTM) network combined with convolutional neural networks (CNN). Deep neural networks have been shown to be a significant way for these types of problems because of their complexity and huge number of trainable paramters. Hybrid method that proposed in the article could significantly increase the overall accuracy of NILM because it benefits from both network advantages. The proposed method used sequence-to-sequence learning, where the input is a window of the mains and the output is a window of the target appliance. The proposed deep neural network approach has been applied to real-world household energy dataset "REFIT". The REFIT electrical load measurements dataset described in this paper includes whole house aggregate loads and nine individual appliance measurements at 8-second intervals per house, collected continuously over a period of two years from 20 houses around the UK. The proposed method achieve significant performance, improving accuracy and F1-score measures by 95.93% and 80.93% ,respectively which demonstrates the effectiveness and superiority of the proposed approach for home energy monitoring. Comparison of proposed method and other recently published method has been presented and discussed based on accuracy, number of considered appliances and size of the deep neural network trainable parameters. The proposed method shows remarkable performance compare to other previous methods.


title: Skeletal Feature Compensation for Imitation Learning with Embodiment Mismatch

abstract: Learning from demonstrations in the wild (e.g. YouTube videos) is a tantalizing goal in imitation learning. However, for this goal to be achieved, imitation learning algorithms must deal with the fact that the demonstrators and learners may have bodies that differ from one another. This condition -- "embodiment mismatch" -- is ignored by many recent imitation learning algorithms. Our proposed imitation learning technique, SILEM (\textbf{S}keletal feature compensation for \textbf{I}mitation \textbf{L}earning with \textbf{E}mbodiment \textbf{M}ismatch), addresses a particular type of embodiment mismatch by introducing a learned affine transform to compensate for differences in the skeletal features obtained from the learner and expert. We create toy domains based on PyBullet's HalfCheetah and Ant to assess SILEM's benefits for this type of embodiment mismatch. We also provide qualitative and quantitative results on more realistic problems -- teaching simulated humanoid agents, including Atlas from Boston Dynamics, to walk by observing human demonstrations.


title: Ensemble gradient for learning turbulence models from indirect observations

abstract: Training data-driven turbulence models with high fidelity Reynolds stress can be impractical and recently such models have been trained with velocity and pressure measurements. For gradient-based optimization, such as training deep learning models, this requires evaluating the sensitivities of the RANS equations. This paper explores the use of an ensemble approximation of the sensitivities of the RANS equations in training data-driven turbulence models with indirect observations. A deep neural network representing the turbulence model is trained using the network's gradients obtained by backpropagation and the ensemble approximation of the RANS sensitivities. Different ensemble approximations are explored and a method based on explicit projection onto the sample space is presented. As validation, the gradient approximations from the different methods are compared to that from the continuous adjoint equations. The ensemble approximation is then used to learn different turbulence models from velocity observations. In all cases, the learned model predicts improved velocities. However, it was observed that once the sensitivity of the velocity to the underlying model becomes small, the approximate nature of the ensemble gradient hinders further optimization of the underlying model. The benefits and limitations of the ensemble gradient approximation are discussed, in particular as compared to the adjoint equations.


title: Resolving star-forming clumps in a z $\sim$ 2 lensed galaxy: a pixelated Bayesian approach

abstract: We present a pixelized source reconstruction method applied on Integral Field Spectroscopic (IFS) observations of gravitationally lensed galaxies. We demonstrate the effectiveness of this method in a case study on the clumpy morphology of a $z \sim 2$ lensed galaxy behind a group-scale lens. We use a Bayesian forward source modelling approach to reconstruct the surface brightness distribution of the source galaxy on a uniformly pixelized grid while accounting for the image point spread function (PSF). The pixelated approach is sensitive to clump sizes down to 100 pc and resolves smaller clump sizes with an improvement in the signal to noise ratio (SNR) by almost a factor of ten compared with more traditional ray-tracing approaches.


title: Accelerating parameter estimation of gravitational waves from compact binary coalescence using adaptive frequency resolutions

abstract: Bayesian parameter estimation of gravitational waves from compact binary coalescence (CBC) typically requires more than millions of evaluations of computationally expensive template waveforms. We propose a technique to reduce the cost of waveform generation by exploiting the chirping behavior of CBC signal. Our technique does not require waveforms at all frequencies in the frequency range used in the analysis, and does not suffer from the fixed cost due to the upsampling of waveforms. Our technique speeds up the parameter estimation of typical binary neutron star signal by a factor of $\mathcal{O}(10)$ for the low-frequency cutoff of $20\,\mathrm{Hz}$, and $\mathcal{O}(10^2)$ for $5\,\mathrm{Hz}$. It does not require any offline preparations or accurate estimates of source parameters provided by detection pipelines.


title: Detecting Polarized Topics Using Partisanship-aware Contextualized Topic Embeddings

abstract: Growing polarization of the news media has been blamed for fanning disagreement, controversy and even violence. Early identification of polarized topics is thus an urgent matter that can help mitigate conflict. However, accurate measurement of topic-wise polarization is still an open research challenge. To address this gap, we propose Partisanship-aware Contextualized Topic Embeddings (PaCTE), a method to automatically detect polarized topics from partisan news sources. Specifically, utilizing a language model that has been finetuned on recognizing partisanship of the news articles, we represent the ideology of a news corpus on a topic by corpus-contextualized topic embedding and measure the polarization using cosine distance. We apply our method to a dataset of news articles about the COVID-19 pandemic. Extensive experiments on different news sources and topics demonstrate the efficacy of our method to capture topical polarization, as indicated by its effectiveness of retrieving the most polarized topics.


title: A Method to Reveal Speaker Identity in Distributed ASR Training, and How to Counter It

abstract: End-to-end Automatic Speech Recognition (ASR) models are commonly trained over spoken utterances using optimization methods like Stochastic Gradient Descent (SGD). In distributed settings like Federated Learning, model training requires transmission of gradients over a network. In this work, we design the first method for revealing the identity of the speaker of a training utterance with access only to a gradient. We propose Hessian-Free Gradients Matching, an input reconstruction technique that operates without second derivatives of the loss function (required in prior works), which can be expensive to compute. We show the effectiveness of our method using the DeepSpeech model architecture, demonstrating that it is possible to reveal the speaker's identity with 34% top-1 accuracy (51% top-5 accuracy) on the LibriSpeech dataset. Further, we study the effect of two well-known techniques, Differentially Private SGD and Dropout, on the success of our method. We show that a dropout rate of 0.2 can reduce the speaker identity accuracy to 0% top-1 (0.5% top-5).


title: Local Zeta Functions From Calabi-Yau Differential Equations

abstract: The zeta-function of a manifold is closely related to, and sometimes can be calculated completely, in terms of its periods. We report here on a practical and computationally rapid implementation of this procedure for families of Calabi-Yau manifolds with one complex structure parameter phi. Although partly conjectural, it turns out to be possible to compute the matrix of the Frobenius map on the third cohomology group of X(phi) directly from the Picard-Fuchs differential operator of the family. To illustrate our method, we compute tables of the quartic numerators of the zeta-functions for six manifolds of increasing complexity as the parameter phi varies in Fp. For four of these manifolds, we do this for the 500 primes p=5,7,...,3583, while for two manifolds we extend the calculation to 1000 primes. The tables for 5 <= p <= 97 are part of this article while the remaining tables are attached in electronic form. Interest attaches to the cases for which the numerators factorise. Some of these factorisations can be associated with parameter values for which the underlying manifold becomes singular. For the cases we consider here, the singularities are all of conifold or hyperconifold type. In these cases the numerator degenerates to a cubic and this factorises into the product of a linear and a quadratic factor. The quadratic term contains a coefficient that is the p'th coefficient of a modular form. Some of our examples have singularities when the parameter satisfies a polynomial equation that does not factorise over Q. When this happens, the corresponding forms are modular forms with neben type or Hilbert modular forms. The numerator can also factorise into two quadrics. This happens when the Hodge structure of the manifold splits, sometimes this happens for algebraic values of the parameter and we identify, in this way, attractor points of rank two of the parameter space.


title: Quantum advantage and noise reduction in distributed quantum computing

abstract: Distributed quantum computing can give substantial noise reduction due to shallower circuits. An experiment illustrates the advantages in the case of Grover search. This motivates studying the quantum advantage of the distributed version of the Simon and Deutsch-Jozsa algorithm. We show that the distributed Simon algorithm retains the exponential advantage, but the complexity deteriorates from O(n) to O(n^2), where n = log2(N). The distributed Deutsch-Jozsa deteriorates to being probabilistic but retains a quantum advantage over classical random sampling.


title: Surface science motivated by heating of trapped ions from the quantum ground state

abstract: For the past two and a half decades, anomalous heating of trapped ions from nearby electrode surfaces has continued to demonstrate unexpected results. Caused by electric-field noise, this heating of the ions' motional modes remains an obstacle for scalable quantum computation with trapped ions. One of the anomalous features of this electric-field noise is the reported nonmonotonic behavior in the heating rate when a trap is incrementally cleaned by ion bombardment. Motivated by this result, the present work reports on a surface analysis of a sample ion-trap electrode treated similarly with incremental doses of Ar$^+$ ion bombardment. Kelvin probe force microscopy and x-ray photoelectron spectroscopy were used to investigate how the work functions on the electrode surface vary depending on the residual contaminant coverage between each treatment. It is shown that the as-fabricated Au electrode is covered with a hydrocarbon film that is modified after the first treatment, resulting in work functions and core-level binding energies that resemble that of atomic-like carbon on Au. Changes in the spatial distribution of work functions with each treatment, combined with a suggested phenomenological coverage and surface-potential roughness dependence to the heating, appear to be related to the nonmonotonic behavior previously reported.


title: Out-of-Distribution Detection for Dermoscopic Image Classification

abstract: Medical image diagnosis can be achieved by deep neural networks, provided there is enough varied training data for each disease class. However, a hitherto unknown disease class not encountered during training will inevitably be misclassified, even if predicted with low probability. This problem is especially important for medical image diagnosis, when an image of a hitherto unknown disease is presented for diagnosis, especially when the images come from the same image domain, such as dermoscopic skin images. Current out-of-distribution detection algorithms act unfairly when the in-distribution classes are imbalanced, by favouring the most numerous disease in the training sets. This could lead to false diagnoses for rare cases which are often medically important. We developed a novel yet simple method to train neural networks, which enables them to classify in-distribution dermoscopic skin disease images and also detect novel diseases from dermoscopic images at test time. We show that our BinaryHeads model not only does not hurt classification balanced accuracy when the data is imbalanced, but also consistently improves the balanced accuracy. We also introduce an important method to investigate the effectiveness of out-of-distribution detection methods based on presence of varying amounts of out-of-distribution data, which may arise in real-world settings.


title: Machine Learning Approaches for Type 2 Diabetes Prediction and Care Management

abstract: Prediction of diabetes and its various complications has been studied in a number of settings, but a comprehensive overview of problem setting for diabetes prediction and care management has not been addressed in the literature. In this document we seek to remedy this omission in literature with an encompassing overview of diabetes complication prediction as well as situating this problem in the context of real world healthcare management. We illustrate various problems encountered in real world clinical scenarios via our own experience with building and deploying such models. In this manuscript we illustrate a Machine Learning (ML) framework for addressing the problem of predicting Type 2 Diabetes Mellitus (T2DM) together with a solution for risk stratification, intervention and management. These ML models align with how physicians think about disease management and mitigation, which comprises these four steps: Identify, Stratify, Engage, Measure.


title: Investigating the Relative Gas and Small Dust Grain Surface Heights in Protoplanetary Disks

abstract: Dust evolution in protoplanetary disks from small dust grains to pebbles is key to the planet formation process. The gas in protoplanetary disks should influence the vertical distribution of small dust grains ($\sim$1 $\mu m$) in the disk.Utilizing archival near-infrared polarized light and millimeter observations, we can measure the scale height and the flare parameter $\beta$ of the small dust grain scattering surface and $^{12}$CO gas emission surface for three protoplanetary disks IM Lup, HD 163296, and HD 97048 (CU Cha). For two systems, IM Lup and HD 163296, the $^{12}$CO gas and small dust grains at small radii from the star have similar heights but at larger radii ($>$100 au) the dust grain scattering surface height is lower than the $^{12}$CO gas emission surface height. In the case of HD 97048, the small dust grain scattering surface has similar heights to the $^{12}$CO gas emission surface at all radii. We ran a protoplanetary disk radiative transfer model of a generic protoplanetary disk with TORUS and showed that there is no difference between the observed scattering surface and $^{12}$CO emission surface. We also performed analytical modeling of the system and found that gas-to-dust ratios larger than 100 could explain the observed difference in IM Lup and HD 163296. This is the first direct comparison of observations of gas and small dust grain heights distribution in protoplanetary disks. Future observations of gas emission and near-infrared scattered light instruments are needed to look for similar trends in other protoplanetary disks.


title: Estimating and Improving Dynamic Treatment Regimes With a Time-Varying Instrumental Variable

abstract: Estimating dynamic treatment regimes (DTRs) from retrospective observational data is challenging as some degree of unmeasured confounding is often expected. In this work, we develop a framework of estimating properly defined "optimal" DTRs with a time-varying instrumental variable (IV) when unmeasured covariates confound the treatment and outcome, rendering the potential outcome distributions only partially identified. We derive a novel Bellman equation under partial identification, use it to define a generic class of estimands (termed IV-optimal DTRs), and study the associated estimation problem. We then extend the IV-optimality framework to tackle the policy improvement problem, delivering IV-improved DTRs that are guaranteed to perform no worse and potentially better than a pre-specified baseline DTR. Importantly, our IV-improvement framework opens up the possibility of strictly improving upon DTRs that are optimal under the no unmeasured confounding assumption (NUCA). We demonstrate via extensive simulations the superior performance of IV-optimal and IV-improved DTRs over the DTRs that are optimal only under the NUCA. In a real data example, we embed retrospective observational registry data into a natural, two-stage experiment with noncompliance using a time-varying IV and estimate useful IV-optimal DTRs that assign mothers to high-level or low-level neonatal intensive care units based on their prognostic variables.


title: Digital quantum simulation of open quantum systems using quantum imaginary time evolution

abstract: Quantum simulation on emerging quantum hardware is a topic of intense interest. While many studies focus on computing ground state properties or simulating unitary dynamics of closed systems, open quantum systems are an interesting target of study owing to their ubiquity and rich physical behavior. However, their non-unitary dynamics are also not natural to simulate on digital quantum devices. Here, we report algorithms for the digital quantum simulation of the dynamics of open quantum systems governed by a Lindblad equation using adaptations of the quantum imaginary time evolution (QITE) algorithm. We demonstrate the algorithms on IBM Quantum's hardware with simulations of the spontaneous emission of a two level system and the dissipative transverse field Ising model. Our work advances efforts to simulate the dynamics of open quantum systems on quantum hardware.


title: NePTuNe: Neural Powered Tucker Network for Knowledge Graph Completion

abstract: Knowledge graphs link entities through relations to provide a structured representation of real world facts. However, they are often incomplete, because they are based on only a small fraction of all plausible facts. The task of knowledge graph completion via link prediction aims to overcome this challenge by inferring missing facts represented as links between entities. Current approaches to link prediction leverage tensor factorization and/or deep learning. Factorization methods train and deploy rapidly thanks to their small number of parameters but have limited expressiveness due to their underlying linear methodology. Deep learning methods are more expressive but also computationally expensive and prone to overfitting due to their large number of trainable parameters. We propose Neural Powered Tucker Network (NePTuNe), a new hybrid link prediction model that couples the expressiveness of deep models with the speed and size of linear models. We demonstrate that NePTuNe provides state-of-the-art performance on the FB15K-237 dataset and near state-of-the-art performance on the WN18RR dataset.


title: Tidal deformability of strange stars and the GW170817 event

abstract: In this work we consider strange stars formed by quark matter in the color-flavor-locked (CFL) phase of color superconductivity. The CFL phase is described by a Nambu-Jona-Lasinio model with four-fermion vector and diquark interaction channels. The effect of the color superconducting medium on the gluons are incorporated into the model by including the gluon self-energy in the thermodynamic potential. We construct parametrizations of the model by varying the vector coupling $G_V$ and comparing the results to the data on tidal deformability from the GW170817 event, the observational data on maximum masses from massive pulsars such as the MSP J0740+6620, and the mass/radius fits to NICER data for PSR J003+0451. Our results points out to windows for the $G_V$ parameter space of the model, with and without gluon effects included, that are compatible with all these astrophysical constraints, namely, $0.21<G_V/G_S<0.4$, and $0.02<G_V/G_S<0.1$, respectively. We also observe a strong correlation between the tidal deformabilites of the GW170817 event and $G_V$. Our results indicate that strange stars cannot be ruled out in collisions of compact binaries from the structural point of view.


title: A New Scheme for Solving High-Order DG Discretizations of Thermal Radiative Transfer using the Variable Eddington Factor Method

abstract: We present a new approach for solving high-order thermal radiative transfer (TRT) using the Variable Eddington Factor (VEF) method (also known as quasidiffusion). Our approach leverages the VEF equations, which consist of the first and second moments of the $S_N$ transport equation, to more efficiently compute the TRT solution for each time step. The scheme consists of two loops - an outer loop to converge the Eddington tensor and an inner loop to converge the iteration between the temperature equation and the VEF system. By converging the outer iteration, one obtains the fully implicit TRT solution for the given time step with a relatively low number of transport sweeps. However, one could choose to perform exactly one outer iteration (and therefore exactly one sweep) per time step, resulting in a semi-implicit scheme that is both highly efficient and robust. Our results indicate that the error between the one-sweep and fully implicit variants of our scheme may be small enough for consideration in many problems of interest.


title: Optimization of High Entropy Alloy Catalyst for Ammonia Decomposition and Ammonia Synthesis

abstract: The successful synthesis of high entropy alloy (HEA) nanoparticles, a long-sought goal in materials science, opens a new frontier in materials science with applications across catalysis, electronics, structural alloys, and energetic materials. Recently, a Co25Mo45Fe10Ni10Cu10 HEA made of earth-abundant elements was shown to have a high catalytic activity for ammonia decomposition, which rivals that of state-of-the-art, but prohibitively expensive, ruthenium catalyst. Using a computational approach based on first-principles calculations in conjunction with data analytics and machine learning, we build a model to rapidly compute the adsorption energy of H, N, and NHx (x=1,3) species on CoMoFeNiCu alloy surfaces with varied alloy compositions and atomic arrangement. We show that the 25/45 Co/Mo ratio identified experimentally as the most active composition for ammonia decomposition increases the likelihood that the surface adsorbs nitrogen equivalently to that of ruthenium while at the same time interacting moderately strongly with intermediates. Our study underscores the importance of computational modeling and machine learning to identify and optimize HEA alloys across their near-infinite materials design space.


title: On $L_1$-embeddability of unions of $L_1$-embeddable metric spaces and of twisted unions of hypercubes

abstract: We study properties of twisted unions of metric spaces introduced by Johnson, Lindenstrauss, and Schechtman, and by Naor and Rabani. In particular, we prove that under certain natural mild assumptions twisted unions of $L_1$-embeddable metric spaces also embed in $L_1$ with distortions bounded above by constants that do not depend on the metric spaces themselves, or on their size, but only on certain general parameters. This answers a question stated by Naor and by Naor and Rabani. In the second part of the paper we give new simple examples of metric spaces such their every embedding into $L_p$, $1\le p<\infty$, has distortion at least $3$, but which are a union of two subsets, each isometrically embeddable in $L_p$. This extends an analogous result of K.~Makarychev and Y.~Makarychev from Hilbert spaces to $L_p$-spaces, $1\le p<\infty$.


title: A Masked Segmental Language Model for Unsupervised Natural Language Segmentation

abstract: Segmentation remains an important preprocessing step both in languages where "words" or other important syntactic/semantic units (like morphemes) are not clearly delineated by white space, as well as when dealing with continuous speech data, where there is often no meaningful pause between words. Near-perfect supervised methods have been developed for use in resource-rich languages such as Chinese, but many of the world's languages are both morphologically complex, and have no large dataset of "gold" segmentations into meaningful units. To solve this problem, we propose a new type of Segmental Language Model (Sun and Deng, 2018; Kawakami et al., 2019; Wang et al., 2021) for use in both unsupervised and lightly supervised segmentation tasks. We introduce a Masked Segmental Language Model (MSLM) built on a span-masking transformer architecture, harnessing the power of a bi-directional masked modeling context and attention. In a series of experiments, our model consistently outperforms Recurrent SLMs on Chinese (PKU Corpus) in segmentation quality, and performs similarly to the Recurrent model on English (PTB). We conclude by discussing the different challenges posed in segmenting phonemic-type writing systems.


title: Pylot: A Modular Platform for Exploring Latency-Accuracy Tradeoffs in Autonomous Vehicles

abstract: We present Pylot, a platform for autonomous vehicle (AV) research and development, built with the goal to allow researchers to study the effects of the latency and accuracy of their models and algorithms on the end-to-end driving behavior of an AV. This is achieved through a modular structure enabled by our high-performance dataflow system that represents AV software pipeline components (object detectors, motion planners, etc.) as a dataflow graph of operators which communicate on data streams using timestamped messages. Pylot readily interfaces with popular AV simulators like CARLA, and is easily deployable to real-world vehicles with minimal code changes. To reduce the burden of developing an entire pipeline for evaluating a single component, Pylot provides several state-of-the-art reference implementations for the various components of an AV pipeline. Using these reference implementations, a Pylot-based AV pipeline is able to drive a real vehicle, and attains a high score on the CARLA Autonomous Driving Challenge. We also present several case studies enabled by Pylot, including evidence of a need for context-dependent components, and per-component time allocation. Pylot is open source, with the code available at https://github.com/erdos-project/pylot.


title: Human-like informative conversations: Better acknowledgements using conditional mutual information

abstract: This work aims to build a dialogue agent that can weave new factual content into conversations as naturally as humans. We draw insights from linguistic principles of conversational analysis and annotate human-human conversations from the Switchboard Dialog Act Corpus to examine humans strategies for acknowledgement, transition, detail selection and presentation. When current chatbots (explicitly provided with new factual content) introduce facts into a conversation, their generated responses do not acknowledge the prior turns. This is because models trained with two contexts - new factual content and conversational history - generate responses that are non-specific w.r.t. one of the contexts, typically the conversational history. We show that specificity w.r.t. conversational history is better captured by Pointwise Conditional Mutual Information ($\text{pcmi}_h$) than by the established use of Pointwise Mutual Information ($\text{pmi}$). Our proposed method, Fused-PCMI, trades off $\text{pmi}$ for $\text{pcmi}_h$ and is preferred by humans for overall quality over the Max-PMI baseline 60% of the time. Human evaluators also judge responses with higher $\text{pcmi}_h$ better at acknowledgement 74% of the time. The results demonstrate that systems mimicking human conversational traits (in this case acknowledgement) improve overall quality and more broadly illustrate the utility of linguistic principles in improving dialogue agents.


title: CRKSPH-compatible discretization of the SUPG and SAAF transport equations

abstract: The self-adjoint angular flux and streamline-upwind Petrov-Galerkin transport equations are discretized using reproducing kernels with the collocation method to produce a discretization that is compatible with conservative reproducing kernel smoothed particle hydrodynamics. A novel second derivative is derived for the diffusion-like term in the self-adjoint angular flux equation. The resulting equations involve only evaluations of kernels and physical data at the nodal centers.


title: Ab initio determination of ultrahigh thermal conductivity in ternary compounds

abstract: Discovering new materials with ultrahigh thermal conductivity has been a critical research frontier and driven by many important technological applications ranging from thermal management to energy science. Here we have rigorously investigated the fundamental lattice vibrational spectra in ternary compounds and determined the thermal conductivity using a predictive ab initio approach. Phonon transport in B-X-C (X = N, P, As) groups is systematically quantified with different crystal structures and high-order anharmonicity involving a four-phonon process. Our calculation found an ultrahigh room-temperature thermal conductivity through strong carbon-carbon bonding up to 2100 W/mK beyond most common materials and the recently discovered boron arsenide. This study provides fundamental insight into the atomistic design of thermal conductivity and opens up opportunities in new materials searching towards complicated compound structures. DOI: 10.1103/PhysRevB.103.L041203


title: re:Linde et al. (2021): The Bayes factor, HDI-ROPE and frequentist equivalence tests can all be reverse engineered -- almost exactly -- from one another

abstract: Following an extensive simulation study comparing the operating characteristics of three different procedures used for establishing equivalence (the frequentist "TOST", the Bayesian "HDI-ROPE", and the Bayes factor interval null procedure), Linde et al. (2021) conclude with the recommendation that "researchers rely more on the Bayes factor interval null approach for quantifying evidence for equivalence." We redo the simulation study of Linde et al. (2021) in its entirety but with the different procedures calibrated to have the same predetermined maximum type 1 error rate. Our results suggest that, when calibrated in this way, the Bayes Factor, HDI-ROPE, and frequentist equivalence tests all have similar -- almost exactly -- type 2 error rates. In general any advocating for frequentist testing as better or worse than Bayesian testing in terms of empirical findings seems dubious at best. If one decides on which underlying principle to subscribe to in tackling a given problem, then the method follows naturally. Bearing in mind that each procedure can be reverse-engineered from the others (at least approximately), trying to use empirical performance to argue for one approach over another seems like tilting at windmills.


title: Dynamical sweet spot engineering via two-tone flux modulation of superconducting qubits

abstract: Current superconducting quantum processors require strategies for coping with material defects and imperfect parameter targeting in order to scale up while maintaining high performance. To that end, in-situ control of qubit frequencies with magnetic flux can be used to avoid spurious resonances. However, increased dephasing due to 1/f flux noise limits performance at all of these operating points except for noise-protected sweet spots, which are sparse under DC flux bias and monochromatic flux modulation. Here we experimentally demonstrate that two-tone flux modulation can be used to create a continuum of dynamical sweet spots, greatly expanding the range of qubit frequencies achievable while first-order insensitive to slow flux noise. To illustrate some advantages of this flexibility, we use bichromatic flux control to reduce the error rates and gate times of parametric entangling operations between transmons. Independent of gate scheme, the ability to use flux control to freely select qubit frequencies while maintaining qubit coherence represents an important step forward in the robustness and scalability of near-term superconducting qubit devices.


title: Tracing Topic Transitions with Temporal Graph Clusters

abstract: Twitter serves as a data source for many Natural Language Processing (NLP) tasks. It can be challenging to identify topics on Twitter due to continuous updating data stream. In this paper, we present an unsupervised graph based framework to identify the evolution of sub-topics within two weeks of real-world Twitter data. We first employ a Markov Clustering Algorithm (MCL) with a node removal method to identify optimal graph clusters from temporal Graph-of-Words (GoW). Subsequently, we model the clustering transitions between the temporal graphs to identify the topic evolution. Finally, the transition flows generated from both computational approach and human annotations are compared to ensure the validity of our framework.


title: Cross-lingual Entity Alignment with Adversarial Kernel Embedding and Adversarial Knowledge Translation

abstract: Cross-lingual entity alignment, which aims to precisely connect the same entities in different monolingual knowledge bases (KBs) together, often suffers challenges from feature inconsistency to sequence context unawareness. This paper presents a dual adversarial learning framework for cross-lingual entity alignment, DAEA, with two original contributions. First, in order to address the structural and attribute feature inconsistency between entities in two knowledge graphs (KGs), an adversarial kernel embedding technique is proposed to extract graph-invariant information in an unsupervised manner, and project two KGs into the common embedding space. Second, in order to further improve successful rate of entity alignment, we propose to produce multiple random walks through each entity to be aligned and mask these entities in random walks. With the guidance of known aligned entities in the context of multiple random walks, an adversarial knowledge translation model is developed to fill and translate masked entities in pairwise random walks from two KGs. Extensive experiments performed on real-world datasets show that DAEA can well solve the feature inconsistency and sequence context unawareness issues and significantly outperforms thirteen state-of-the-art entity alignment methods.


title: Investigating Failures of Automatic Translation in the Case of Unambiguous Gender

abstract: Transformer based models are the modern work horses for neural machine translation (NMT), reaching state of the art across several benchmarks. Despite their impressive accuracy, we observe a systemic and rudimentary class of errors made by transformer based models with regards to translating from a language that doesn't mark gender on nouns into others that do. We find that even when the surrounding context provides unambiguous evidence of the appropriate grammatical gender marking, no transformer based model we tested was able to accurately gender occupation nouns systematically. We release an evaluation scheme and dataset for measuring the ability of transformer based NMT models to translate gender morphology correctly in unambiguous contexts across syntactically diverse sentences. Our dataset translates from an English source into 20 languages from several different language families. With the availability of this dataset, our hope is that the NMT community can iterate on solutions for this class of especially egregious errors.


title: Removing non-smoothness in solving Black-Scholes equation using a perturbation method

abstract: Black-Scholes equation as one of the most celebrated mathematical models has an explicit analytical solution known as the Black-Scholes formula. Later variations of the equation, such as fractional or nonlinear Black-Scholes equations, do not have a closed form expression for the corresponding formula. In that case, one will need asymptotic expansions, including homotopy perturbation method, to give an approximate analytical solution. However, the solution is non-smooth at a special point. We modify the method by {first} performing variable transformations that push the point to infinity. As a test bed, we apply the method to the solvable Black-Scholes equation, where excellent agreement with the exact solution is obtained. We also extend our study to multi-asset basket and quanto options by reducing the cases to single-asset ones. Additionally we provide a novel analytical solution of the single-asset quanto option that is simple and different from the existing expression.


title: Are Classes Clusters?

abstract: Sentence embedding models aim to provide general purpose embeddings for sentences. Most of the models studied in this paper claim to perform well on STS tasks - but they do not report on their suitability for clustering. This paper looks at four recent sentence embedding models (Universal Sentence Encoder (Cer et al., 2018), Sentence-BERT (Reimers and Gurevych, 2019), LASER (Artetxe and Schwenk, 2019), and DeCLUTR (Giorgi et al., 2020)). It gives a brief overview of the ideas behind their implementations. It then investigates how well topic classes in two text classification datasets (Amazon Reviews (Ni et al., 2019) and News Category Dataset (Misra, 2018)) map to clusters in their corresponding sentence embedding space. While the performance of the resulting classification model is far from perfect, it is better than random. This is interesting because the classification model has been constructed in an unsupervised way. The topic classes in these real life topic classification datasets can be partly reconstructed by clustering the corresponding sentence embeddings.


title: Pareto Self-Supervised Training for Few-Shot Learning

abstract: While few-shot learning (FSL) aims for rapid generalization to new concepts with little supervision, self-supervised learning (SSL) constructs supervisory signals directly computed from unlabeled data. Exploiting the complementarity of these two manners, few-shot auxiliary learning has recently drawn much attention to deal with few labeled data. Previous works benefit from sharing inductive bias between the main task (FSL) and auxiliary tasks (SSL), where the shared parameters of tasks are optimized by minimizing a linear combination of task losses. However, it is challenging to select a proper weight to balance tasks and reduce task conflict. To handle the problem as a whole, we propose a novel approach named as Pareto self-supervised training (PSST) for FSL. PSST explicitly decomposes the few-shot auxiliary problem into multiple constrained multi-objective subproblems with different trade-off preferences, and here a preference region in which the main task achieves the best performance is identified. Then, an effective preferred Pareto exploration is proposed to find a set of optimal solutions in such a preference region. Extensive experiments on several public benchmark datasets validate the effectiveness of our approach by achieving state-of-the-art performance.


title: Counting Phylogenetic Networks with Few Reticulation Vertices: A Second Approach

abstract: Tree-child networks, one of the prominent network classes in phylogenetics, have been introduced for the purpose of modeling reticulate evolution. Recently, the first author together with Gittenberger and Mansouri (2019) showed that the number ${\rm TC}_{\ell,k}$ of tree-child networks with $\ell$ leaves and $k$ reticulation vertices has the first-order asymptotics \[ {\rm TC}_{\ell,k}\sim c_k\left(\frac{2}{e}\right)^{\ell}\ell^{\ell+2k-1},\qquad (\ell\rightarrow\infty). \] Moreover, they also computed $c_k$ for $k=1,2,$ and $3$. In this short note, we give a second approach to the above result which is based on a recent (algorithmic) approach for the counting of tree-child networks due to Cardona and Zhang (2020). This second approach is also capable of giving a simple, closed-form expression for $c_k$, namely, $c_k=2^{k-1}\sqrt{2}/k!$ for all $k\geq 0$.


title: Is there a cap on longevity? A statistical review

abstract: There is sustained and widespread interest in understanding the limit, if any, to the human lifespan. Apart from its intrinsic and biological interest, changes in survival in old age have implications for the sustainability of social security systems. A central question is whether the endpoint of the underlying lifetime distribution is finite. Recent analyses of data on the oldest human lifetimes have led to competing claims about survival and to some controversy, due in part to incorrect statistical analysis. This paper discusses the particularities of such data, outlines correct ways of handling them and presents suitable models and methods for their analysis. We provide a critical assessment of some earlier work and illustrate the ideas through reanalysis of semi-supercentenarian lifetime data. Our analysis suggests that remaining life-length after age 109 is exponentially distributed, and that any upper limit lies well beyond the highest lifetime yet reliably recorded. Lower limits to 95% confidence intervals for the human lifespan are around 130 years, and point estimates typically indicate no upper limit at all.


title: Learning Feature Interactions With and Without Specifications

abstract: Features in product lines and highly configurable systems can interact in ways that are contrary to developers' intent. Current methods to identify such unanticipated feature interactions are costly and inadequate. To address this problem we propose a new approach to learn feature interactions, both in those product lines where constraints on feature combinations are specified and in feature-rich configurable systems where such specifications often are not available. The contribution of the paper is to use program analysis to extract feature-relevant learning models from the source code in order to detect unwanted feature interactions. Where specifications of feature constraints are unavailable, our approach infers the constraints using feature-related data-flow dependency information. Evaluation in experiments on three software product line benchmarks and a highly configurable system shows that this approach is fast and effective. The contribution is to support developers by automatically detecting feature combinations in a new product or version that can interact in unwanted or unrecognized ways. This enables better understanding of latent interactions and identifies software components that should be tested together because their features interact in some configurations.


title: No pure capillary solitary waves exist in 2D finite depth

abstract: We prove that the 2D finite depth capillary water wave equations admit no solitary wave solutions. This closes the existence/non-existence problem for solitary water waves in 2D, under the classical assumptions of incompressibility and irrotationality, and with the physical parameters being gravity, surface tension and the fluid depth.


title: Multivalent Entailment Graphs for Question Answering

abstract: Drawing inferences between open-domain natural language predicates is a necessity for true language understanding. There has been much progress in unsupervised learning of entailment graphs for this purpose. We make three contributions: (1) we reinterpret the Distributional Inclusion Hypothesis to model entailment between predicates of different valencies, like DEFEAT(Biden, Trump) entails WIN(Biden); (2) we actualize this theory by learning unsupervised Multivalent Entailment Graphs of open-domain predicates; and (3) we demonstrate the capabilities of these graphs on a novel question answering task. We show that directional entailment is more helpful for inference than bidirectional similarity on questions of fine-grained semantics. We also show that drawing on evidence across valencies answers more questions than by using only the same valency evidence.


title: Exploring software developers' work practices: Task differences, participation, engagement, and speed of task resolution

abstract: In seeking to understand the processes enacted during software development, an increasing number of studies have mined software repositories. In particular, studies have endeavored to show how teams resolve software defects. Although much of this work has been useful, we contend that large-scale examinations across the range of activities that are commonly performed, beyond defect-related issues alone, would help us to more fully understand the reasons why defects occur as well as their consequences. More generally, these explorations would reveal how team processes occur during all software development efforts. We thus extend such studies by investigating how software practitioners work while undertaking the range of software tasks that are typically performed. Multiple forms of analyses of a longitudinal case study reveal that software practitioners were mostly involved in fixing defects, and that their engagement covaried depending on the nature of the work they were performing. Furthermore, multiple external factors affected speed of task resolution. Our outcomes suggest that behavioral and intrinsic issues may interact with extrinsic factors becoming significant predictors of the speed of software task resolution.


title: Comparison of Grammatical Error Correction Using Back-Translation Models

abstract: Grammatical error correction (GEC) suffers from a lack of sufficient parallel data. Therefore, GEC studies have developed various methods to generate pseudo data, which comprise pairs of grammatical and artificially produced ungrammatical sentences. Currently, a mainstream approach to generate pseudo data is back-translation (BT). Most previous GEC studies using BT have employed the same architecture for both GEC and BT models. However, GEC models have different correction tendencies depending on their architectures. Thus, in this study, we compare the correction tendencies of the GEC models trained on pseudo data generated by different BT models, namely, Transformer, CNN, and LSTM. The results confirm that the correction tendencies for each error type are different for every BT model. Additionally, we examine the correction tendencies when using a combination of pseudo data generated by different BT models. As a result, we find that the combination of different BT models improves or interpolates the F_0.5 scores of each error type compared with that of single BT models with different seeds.


title: Task Space Planning with Complementarity Constraint-based Obstacle Avoidance

abstract: In this paper, we present a task space-based local motion planner that incorporates collision avoidance and constraints on end-effector motion during the execution of a task. Our key technical contribution is the development of a novel kinematic state evolution model of the robot where the collision avoidance is encoded as a complementarity constraint. We show that the kinematic state evolution with collision avoidance can be represented as a Linear Complementarity Problem (LCP). Using the LCP model along with Screw Linear Interpolation (ScLERP) in SE(3), we show that it may be possible to compute a path between two given task space poses by directly moving from the start to the goal pose, even if there are potential collisions with obstacles. The scalability of the planner is demonstrated with experiments using a physical robot. We present simulation and experimental results with both collision avoidance and task constraints to show the efficacy of our approach.


title: Distinct Properties of Vortex Bound States Driven by Temperature

abstract: We investigate the behavior of vortex bound states in the quantum limit by self-consistently solving the Bogoliubov-de Gennes equation. We find that the energies of the vortex bound states deviates from the analytical result $E_\mu=\mu\Delta^2/E_F$ with the half-integer angular momentum $\mu$ in the extreme quantum limit. Specifically, the energy ratio for the first three orders is more close to $1:2:3$ instead of $1:3:5$ at extremely low temperature. The local density of states reveals an Friedel-like behavior associated with that of the pair potential in the extreme quantum limit, which will be smoothed out by thermal effect above a certain temperature even the quantum limit condition, namely $T/T_c<\Delta/E_F$ is still satisfied. Our studies show that the vortex bound states can exhibit very distinct features in different temperature regimes, which provides a comprehensive understanding and should stimulate more experimental efforts for verifications.


title: High-Quality Automated Program Repair

abstract: Automatic program repair (APR) has recently gained attention because it proposes to fix software defects with no human intervention. To automatically fix defects, most APR tools use the developer-written tests to (a) localize the defect, and (b) generate and validate the automatically produced candidate patches based on the constraints imposed by the tests. While APR tools can produce patches that appear to fix the defect for 11-19% of the defects in real-world software, most of the patches produced are not correct or acceptable to developers because they overfit to the tests used during the repair process. This problem is known as the patch overfitting problem. To address this problem, I propose to equip APR tools with additional constraints derived from natural-language software artifacts such as bug reports and requirements specifications that describe the bug and intended software behavior but are not typically used by the APR tools. I hypothesize that patches produced by APR tools while using such additional constraints would be of higher quality. To test this hypothesis, I propose an automated and objective approach to evaluate the quality of patches and propose two novel methods to improve the fault localization and developer-written test suites using natural-language software artifacts. Finally, I propose to use my patch evaluation methodology to analyze the effect of the improved fault localization and test suites on the quality of patches produced by APR tools for real-world defects.


title: Minimal obstructions to $(\infty, k)$-polarity in cographs

abstract: A graph is a cograph if it does not contain a 4-vertex path as an induced subgraph. An $(s, k)$-polar partition of a graph $G$ is a partition $(A, B)$ of its vertex set such that $A$ induces a complete multipartite graph with at most $s$ parts, and $B$ induces the disjoint union of at most $k$ cliques with no other edges. A graph $G$ is said to be $(s, k)$-polar if it admits an $(s, k)$-polar partition. The concepts of $(s, \infty)$-, $(\infty, k)$-, and $(\infty, \infty)$-polar graphs can be analogously defined. Ekim, Mahadev and de Werra pioneered in the research on polar cographs, obtaining forbidden induced subgraph characterizations for $(\infty, \infty)$-polar cographs, as well as for the union of $(\infty, 1)$- and $(1, \infty)$-polar cographs. Recently, a recursive procedure for generating the list of cograph minimal $(s,1)$-polar obstructions for any fixed integer $s$ was found, as well as the complete list of $(\infty, 1)$-polar obstructions. In addition to these results, complete lists of minimal $(s, k)$-polar cograph obstructions are known only for the pair $(2, 2)$. In this work we are concerned with the problem of characterizing $(\infty, k)$-polar cographs for a fixed $k$ through a finite family of forbidden induced subgraphs. As our main result, we provide complete lists of forbidden induced subgraphs for the cases $k=2$ and $k=3$. Additionally, we provide a partial recursive construction for the general case. By considering graph complements, these results extend to $(s, \infty)$-polar cographs.