title: The Quadrupole Moment of Compact Binaries to the Fourth post-Newtonian Order: II. Dimensional Regularization and Renormalization

abstract: The regularization and renormalization of the radiative mass-type quadrupole moment of inspiralling compact binaries (without spins) is investigated at the fourth post-Newtonian (4PN) approximation of general relativity. As clear from the conservative 4PN equations of motion, a dimensional regularization has to be implemented in order to properly treat the non-linear interactions experienced by gravitational waves during their propagation toward future null infinity. By implementing such procedure, we show that the poles coming from the source moment (computed in a companion paper) are exactly cancelled in the radiative moment, as expected for a physical quantity. We thus define and obtain a "renormalized" source quadrupole, three-dimensional by nature, which is an important step towards the computation of the gravitational-wave flux with 4PN accuracy. Furthermore, we explicitly prove the equivalence between the dimensional regularization and the previously used Hadamard partie finie scheme up to the 3PN order.


title: Fractional integral inequalities via Atangana-Baleanu operators for convex and concave functions

abstract: Recently, many fractional integral operators were introduced by different mathematicians. One of these fractional operators, Atangana-Baleanu fractional integral operator, was defined by Atangana and Baleanu in [2]. In this study, firstly, a new identity by using Atangana-Baleanu fractional integral operators are proved. Then, new fractional integral inequalities have been obtained for convex and concave functions with the help of this identity and some certain integral inequalities


title: A universal H\"older estimate up to dimension 4 for stable solutions to half-Laplacian semilinear equations

abstract: We study stable solutions to the equation $(-\Delta)^{1/2} u = f(u)$, posed in a bounded domain of $\mathbb{R}^n$. For nonnegative convex nonlinearities, we prove that stable solutions are smooth in dimensions $n\leq 4$. This result, which was known only for $n=1$, follows from a new interior H\"older estimate that is completely independent of the nonlinearity $f$. A main ingredient in our proof is a new geometric form of the stability condition. It is still unknown for other fractions of the Laplacian and, surprisingly, it requires convexity of the nonlinearity. From it, we deduce higher order Sobolev estimates that allow us to extend the techniques developed by Cabr\'e, Figalli, Ros-Oton, and Serra for the Laplacian. In this way we obtain, besides the H\"older bound for $n\leq 4$, a universal $H^{1/2}$ estimate in all dimensions. Our $L^\infty$ bound is expected to hold for $n\leq 8$, but this has been settled only in the radial case or when $f(u) = \lambda e^u$. For other fractions of the Laplacian, the expected optimal dimension for boundedness of stable solutions has been reached only when $f(u) = \lambda e^u$, even in the radial case.


title: Photon-photon scattering from a UV-complete gravity QFT

abstract: Quantum quadratic gravity (QQG) produces a tree-level differential cross section for $\gamma\gamma\to\gamma\gamma$ that is well-behaved at all energies. From this we can study how the corrections to low energy scattering amplitudes are related to the UV physics, in particular to the exchange of the massive graviparticles. An effective forward scattering amplitude is obtained by separating out the effects of the $t$-channel graviton pole. This is possible due to the UV-completeness, and even though the Froissart bound is not satisfied. We then consider photon-photon scattering to two graviparticles and a further imaginary contribution to the $\gamma\gamma\to\gamma\gamma$ forward scattering amplitude. Unitarity without positivity is a key property of QQG and it impacts all our results.


title: Retrieving Exoplanet Atmospheres using Planetary Infrared Excess: Prospects for the Nightside of WASP-43 b and other Hot Jupiters

abstract: To increase the sample size of future atmospheric characterization efforts, we build on the planetary infrared excess (PIE) technique that has been proposed as a means to detect and characterize the thermal spectra of transiting and non-transiting exoplanets using sufficiently broad wavelength coverage to uniquely constrain the stellar and planetary spectral components from spatially unresolved observations. We performed simultaneous retrievals of stellar and planetary spectra for the archetypal planet WASP-43b in its original configuration and a non-transiting configuration to determine the efficacy of the PIE technique for characterizing the planet's nightside atmospheric thermal structure and composition using typical out-of-transit JWST observations. We found that using PIE with JWST should enable the stellar and planetary spectra to be disentangled with no degeneracies seen between the two flux sources, thus allowing robust constraints on the planet's nightside thermal structure and water abundance to be retrieved. The broad wavelength coverage achieved by combining spectra from NIRISS, NIRSpec, and MIRI enables PIE retrievals that are within 10% of the precision attained using traditional secondary eclipse measurements, although mid-IR observations with MIRI alone may face up to 3.5 times lower precision on the planet's irradiation temperature. For non-transiting planets with unconstrained radius priors, we were able to identify and break the degeneracy between planet radius and irradiation temperature using data that resolved the peak of both the stellar and planetary spectra, thus potentially increasing the number of planets amenable to atmospheric characterization with JWST and other future mission concepts.


title: Contextual Combinatorial Bandits with Changing Action Sets via Gaussian Processes

abstract: We consider a contextual bandit problem with a combinatorial action set and time-varying base arm availability. At the beginning of each round, the agent observes the set of available base arms and their contexts and then selects an action that is a feasible subset of the set of available base arms to maximize its cumulative reward in the long run. We assume that the mean outcomes of base arms are samples from a Gaussian Process (GP) indexed by the context set ${\cal X}$, and the expected reward is Lipschitz continuous in expected base arm outcomes. For this setup, we propose an algorithm called Optimistic Combinatorial Learning and Optimization with Kernel Upper Confidence Bounds (O'CLOK-UCB) and prove that it incurs $\tilde{O}(\sqrt{\lambda^*(K)KT\overline{\gamma}_{T}} )$ regret with high probability, where $\overline{\gamma}_{T}$ is the maximum information gain associated with the set of base arm contexts that appeared in the first $T$ rounds, $K$ is the maximum cardinality of any feasible action over all rounds and $\lambda^*(K)$ is the maximum eigenvalue of all covariance matrices of selected actions up to time $T$, which is a function of $K$. To dramatically speed up the algorithm, we also propose a variant of O'CLOK-UCB that uses sparse GPs. Finally, we experimentally show that both algorithms exploit inter-base arm outcome correlation and vastly outperform the previous state-of-the-art UCB-based algorithms in realistic setups.


title: Mapped discontinuous Galerkin interpolations and sheared boundary conditions

abstract: Translations or, more generally, coordinate transformations of scalar fields arise in several applications, such as weather, accretion disk and magnetized plasma turbulence modeling. In local studies of accretion disks and magnetized plasmas these coordinate transformations consist of an analytical mapping and enter via sheared-shift boundary conditions. This work introduces a discontinuous Galerkin algorithm to compute these coordinate transformations or boundary conditions based on projections and quadrature-free integrals. The procedure is high-order accurate, preserves certain moments exactly and works in multiple dimensions. Tests of the proposed approach with increasing complexity are presented, beginning with translations of one and two dimensional fields, followed by 3D and 5D simulations with sheared (twist-shift) boundary conditions. The results show that the algorithm is (p+1)-order accurate in the DG representation and (p+2)-order accurate in the cell averages, with p being the order of the polynomial basis functions. Quantification of the algorithm's diffusion and, for shearing boundary conditions, discussion of aliasing errors are provided.


title: Measuring chemical likeness of stars with RSCA

abstract: Identification of chemically similar stars using elemental abundances is core to many pursuits within Galactic archaeology. However, measuring the chemical likeness of stars using abundances directly is limited by systematic imprints of imperfect synthetic spectra in abundance derivation. We present a novel data-driven model that is capable of identifying chemically similar stars from spectra alone. We call this Relevant Scaled Component Analysis (RSCA). RSCA finds a mapping from stellar spectra to a representation that optimizes recovery of known open clusters. By design, RSCA amplifies factors of chemical abundance variation and minimizes those of non-chemical parameters, such as instrument systematics. The resultant representation of stellar spectra can therefore be used for precise measurements of chemical similarity between stars. We validate RSCA using 185 cluster stars in 22 open clusters in the APOGEE survey. We quantify our performance in measuring chemical similarity using a reference set of 151,145 field stars. We find that our representation identifies known stellar siblings more effectively than stellar abundance measurements. Using RSCA, 1.8% of pairs of field stars are as similar as birth siblings, compared to 2.3% when using stellar abundance labels. We find that almost all of the information within spectra leveraged by RSCA fits into a two-dimensional basis, which we link to [Fe/H] and alpha-element abundances. We conclude that chemical tagging of stars to their birth clusters remains prohibitive. However, using the spectra has noticeable gain, and our approach is poised to benefit from larger datasets and improved algorithm designs.


title: An Exploration of the Mentorship Needs of Research Software Engineers

abstract: As a newly designated professional title, research software engineers (RSEs) link the two worlds of software engineering and research science. They lack clear development and training opportunities, particularly in the realm of mentoring. In this paper, we discuss mentorship as it pertains to the unique needs of RSEs and propose ways in which organizations and institutions can support mentor/mentee relationships for RSEs


title: Less than 1% of Core-Collapse Supernovae in the local universe occur in elliptical galaxies

abstract: We present observations of three Core-collapse supernovae (CCSNe) in elliptical hosts, detected by the Zwicky Transient Facility Bright Transient Survey (BTS). SN 2019ape is a SN Ic that exploded in the main body of a typical elliptical galaxy. Its properties are consistent with an explosion of a regular SN Ic progenitor. A secondary g-band light curve peak could indicate interaction of the ejecta with circumstellar material (CSM). An H$\alpha$-emitting source at the explosion site suggests a residual local star formation origin. SN 2018fsh and SN 2020uik are SNe II which exploded in the outskirts of elliptical galaxies. SN 2020uik shows typical spectra for SNe II, while SN 2018fsh shows a boxy nebular H$\alpha$ profile, a signature of CSM interaction. We combine these 3 SNe with 7 events from the literature and analyze their hosts as a sample. We present multi-wavelength photometry of the hosts, and compare this to archival photometry of all BTS hosts. Using the spectroscopically complete BTS we conclude that $0.3\%^{+0.3}_{-0.1}$ of all CCSNe occur in elliptical galaxies. We derive star-formation rates and stellar masses for the host-galaxies and compare them to the properties of other SN hosts. We show that CCSNe in ellipticals have larger physical separations from their hosts compared to SNe Ia in elliptical galaxies, and discuss implications for star-forming activity in elliptical galaxies.


title: Geometry and organization of coherent structures in stably stratified atmospheric boundary layers

abstract: Global intermittency is observed in the stably stratified Atmospheric Boundary Layer (ABL) and corresponds to having large nonturbulent flow regions to develop in an otherwise turbulent flow. In this paper, the differences between continuous and intermittent turbulence are quantified with the help of coherent structures. Eight classes of coherent structures are identified from literature, most of which are indicated by scalar criteria derived from velocity fields. A method is developed to geometrically classify structures into three categories: blob-like, tube-like or sheet-like. An alternate definition of the intermittency factor $\gamma$ based on coherent structures is introduced to separate turbulent and nonturbulent parts of a flow. Applying this conditioning technique and the geometrical characterization on direct numerical simulations (DNS) of an Ekman flow, we find the following: (i) structures with similar geometries (either tube-like or sheet-like) are found regardless of the strength of stratification; (ii) global intermittency affects all regions of the ABL - viscous sublayer, buffer layer, inner, and outer layer; (iii) for the highly stratified case, sweep/ejection pairs form well-separated clusters within the viscous sublayer which can possibly explain the abundance of hairpin-like vortices with a particular orientation; (iv) nonturbulent regions are occupied with streamwise velocity fluctuations and there is a switch between high- and low-speed streaks at a particular height for all stratified cases.


title: On the fundamental group of non-collapsed ancient Ricci flows

abstract: We show that any manifold admitting a non-collapsed, ancient Ricci flow must have finite fundamental group. This generalizes what was known for $\kappa$-solutions in dimensions 2, 3. We furthermore show that this fundamental group must be a quotient of the fundamental group of the regular part of any tangent flow at infinity.


title: Vertex Algebra constructions for (analytic) Geometric Langlands in genus zero

abstract: We employ two-dimensional chiral algebra techniques to produce solutions of certain differential and integral equations which occur in the context of the Analytic Geometric Langlands Program.


title: Planar solutions of higher-spin theory. Nonlinear corrections

abstract: Leading order higher-spin corrections to the linearized higher-spin black brane are analyzed in four dimensions. It is shown that the static solution that respects planar symmetry exists in the bosonic case at given order. Its higher-spin Weyl tensors are found in a closed form and are shown to have the double copy origin. The effect of higher-spin fields to form a strictly positive scalar condensate for any values of higher-spin charges is observed.


title: Ab-initio investigation of finite size effects in rutile titania nanoparticles with semilocal and nonlocal density functionals

abstract: In this work, we employ hybrid and generalized gradient approximation (GGA) level density functional theory (DFT) calculations to investigate the convergence of surface properties and band structure of rutile titania (TiO$_2$) nanoparticles with particle size. The surface energies and band structures are calculated for cuboidal particles with minimum dimension ranging from 3.7 \r{A} (24 atoms) to 10.3 \r{A} (384 atoms) using a highly-parallel real-space DFT code to enable hybrid level DFT calculations of larger nanoparticles than are typically practical. We deconvolute the geometric and electronic finite size effects in surface energy, and evaluate the influence of defects on band structure and density of states (DOS). The electronic finite size effects in surface energy vanish when the minimum length scale of the nanoparticles becomes greater than 10 \r{A}. We show that this length scale is consistent with a computationally efficient numerical analysis of the characteristic length scale of electronic interactions. The surface energy of nanoparticles having minimum dimension beyond this characteristic length can be approximated using slab calculations that account for the geometric defects. In contrast, the finite size effects on the band structure is highly dependent on the shape and size of these particles. The DOS for cuboidal particles and more realistic particles constructed using the Wulff algorithm reveal that defect states within the bandgap play a key role in determining the band structure of nanoparticles and the bandgap does not converge to the bulk limit for the particle sizes investigated.


title: Thumb Stabilization and Assistance in a Robotic Hand Orthosis for Post-Stroke Hemiparesis

abstract: We propose a dual-cable method of stabilizing the thumb in the context of a hand orthosis designed for individuals with upper extremity hemiparesis after stroke. This cable network adds opposition/reposition capabilities to the thumb, and increases the likelihood of forming a hand pose that can successfully manipulate objects. In addition to a passive-thumb version (where both cables are of fixed length), our approach also allows for a single-actuator active-thumb version (where the extension cable is actuated while the abductor remains passive), which allows a range of motion intended to facilitate creating and maintaining grasps. We performed experiments with five chronic stroke survivors consisting of unimanual resistive-pull tasks and bimanual twisting tasks with simulated real-world objects; these explored the effects of thumb assistance on grasp stability and functional range of motion. Our results show that both active- and passive-thumb versions achieved similar performance in terms of improving grasp force generation over a no-device baseline, but active thumb stabilization enabled users to maintain grasps for longer durations.


title: Multi-Modal Attack Detection for Cyber-Physical Additive Manufacturing

abstract: Cyber-Physical Additive Manufacturing (AM) constructs a physical 3D object layer-by-layer according to its digital representation and has been vastly applied to fast prototyping and the manufacturing of functional end-products across fields. The computerization of traditional production processes propels these technological advancements; however, this also introduces new vulnerabilities, necessitating the study of cyberattacks on these systems. The AM Sabotage Attack is one kind of kinetic cyberattack that originates from the cyber domain and can eventually lead to physical damage, injury, or even death. By introducing inconspicuous yet damaging alterations in any specific process of the AM digital process chain, the attackers can compromise the structural integrity of a manufactured component in a manner that is invisible to a human observer. If the manufactured objects are critical for their system, those attacks can even compromise the whole system's structural integrity and pose a severe safety risk to its users. For example, an inconspicuous void (less than 1 mm in dimension) placed in the 3D design of a tensile test specimen can reduce its yield load by 14%. However, security studies primarily focus on securing digital assets, overlooking the fact that AM systems are CPSs.


title: An Overview of the Drone Open-Source Ecosystem

abstract: Unmanned aerial systems capable of beyond visual line of sight operation can be organized into a top-down hierarchy of layers including flight supervision, command and control, simulation of systems, operating systems, and physical hardware. Flight supervision includes unmanned air traffic management, flight planning, authorization, and remote identification. Command and control ensure drones can be piloted safely. Simulation of systems concerns how drones may react to different environments and how changing conditions and information provide input to a piloting system. Electronic hardware controlling drone operation is typically accessed using an operating system. Each layer in the hierarchy has an ecosystem of open-source solutions. In this brief survey we describe representative open-source examples for each level of the hierarchy.


title: Two Sites of r-Process Production Assessed on the Basis of the Age-tagged Abundances of Solar Twins

abstract: Solar twins, i.e., stars that are nearly identical to the Sun, including their metallicities, in the solar vicinity show ages widely distributed from 0-10 Gyr. This fact matches the orbital history of solar twins in the new paradigm of galactic dynamics, in which stars radially move on the disk when they encounter transient spiral arms. This finding suggests that older twins were born closer to the Galactic center and traveled a longer distance to reach their present location, according to the hypothesis that chemical enrichment occurs more quickly and that solar metallicity is attained on a shorter timescale with a decreasing Galactocentric distance (R_GC). We show that abundance patterns covering a wide range of heavy elements for solar twins sharing similar ages are identical and that their variation among different age groups can be understood on the basis of the age-R_GC connection within the framework of Galactic chemical evolution. This study identifies the Galactic bulge as the birthplace of the oldest solar twins. Based on this scheme, we find that the relation between [r-process/Fe] and R_GC for the inner Galactic region is incompatible with the hypothesis of a sole site for r-process production, that is, neutron star mergers, whose delay time distribution could be approximated by the power-law form (proportional to t^n). Alternatively, this relation suggests the presence of two distinct sites for r-process production: short-lived massive stars, ending with specific core-collapse supernovae, and neutron star mergers that are heavily inclined to emerge with longer delay times, as represented by n=0-0.5.


title: Floquet stability of oscillating boundary layers on adiabatic slopes

abstract: The presence of a no-slip, impermeable, adiabatic, sloped boundary in an otherwise quiescent, stably stratified, Boussinesq flow generates baroclinic vorticity within a diffusive boundary layer. Such conditions are typical of the oscillating boundary layers on adiabatic abyssal slopes, sloped lake bathymetry, and sloped coastal bathymetry in the absence of high-wavenumber internal waves, mean flows, far-field turbulence on larger scales, and resonant tidal-bathymetric interaction. We investigate the linear stability of the oscillating flow within non-dimensional parameter space typical of the $M_2$ tide and hydraulically smooth, mid-latitude abyssal slopes through Floquet linear stability analysis. The flow dynamics depend on three non-dimensional variables: the Reynolds number for Stokes' second problem (Re), the Prandtl number, and a frequency ratio that accounts for the resonance conditions (C, criticality) of the buoyant restoring force and the tidal forcing. The Floquet analysis results suggest that oscillating laminar boundary layers on adiabatic abyssal slopes are increasingly unstable as Reynolds number, criticality parameter, and/or spanwise disturbance wavenumber are increased. We also show that the two-dimensional Floquet linear instability necessarily generates three-dimensional baroclinic vorticity, which suggests that the evolution of the gravitational instabilities may be nonlinear as $t\rightarrow\infty$.


title: Probing modified gravity with magnetically levitated resonators

abstract: We present an experimental procedure, based on Meissner effect levitation of neodymium ferromagnets, as a method of measuring the gravitational interactions between mg masses. The scheme consists of two superconducting lead traps, with a magnet levitating in each trap. The levitating magnets behave as harmonic oscillators, and by carefully driving the motion of one magnet on resonance with the other, we find that it should be easily possible to measure the gravitational field produced by a 4~mg sphere, with the gravitational attraction from masses as small as 30~$\mu$g predicted to be measurable within realistic a realistic measurement time frame. We apply this acceleration sensitivity to one concrete example and show the ability of testing models of modified Newtonian dynamics.


title: Generic Generalized Diagonal Matrices

abstract: Generalized diagonal matrices are matrices that have two ladders of entries that are zero in the upper right and bottom left corners. The minors of generic generalized diagonal matrices have square-free initial ideals. We give a description of the facets of their Stanley-Reisner complex. With this description, we characterize the configuration of ladders that yield Cohen-Macaulay ideals. In the special case where both ladders are triangles, we show that the corresponding complex is vertex decomposable. Also in this case, we compute the height and multiplicity of the ideals.


title: Adaptive Group Testing with Mismatched Models

abstract: Accurate detection of infected individuals is one of the critical steps in stopping any pandemic. When the underlying infection rate of the disease is low, testing people in groups, instead of testing each individual in the population, can be more efficient. In this work, we consider noisy adaptive group testing design with specific test sensitivity and specificity that select the optimal group given previous test results based on pre-selected utility function. As in prior studies on group testing, we model this problem as a sequential Bayesian Optimal Experimental Design (BOED) to adaptively design the groups for each test. We analyze the required number of group tests when using the updated posterior on the infection status and the corresponding Mutual Information (MI) as our utility function for selecting new groups. More importantly, we study how the potential bias on the ground-truth noise of group tests may affect the group testing sample complexity.


title: Time-averaged velocity and scalar fields of the flow surrounding a group of cylinders

abstract: We examine the characteristics of mean flow and scalar concentration characteristics of a turbulent boundary layer flow impinging on a cluster of tall obstacles (which can also be interpreted as a porous obstruction). The cluster is created with a group of cylinders of diameter $d$ and height $h$ arranged in a circular patch of diameter $D$.The solidity of the patch/obstruction is defined by $\phi$ (the total planar area covered by cylinders), which is systematically varied ($0.098 \leq \phi \leq 1$) by increasing the number of cylinders in a patch ($N_c$). A point source is placed at ground level upstream of the patch and its transport around the patch is examined. Time-averaged velocity and scalar fields, obtained from simultaneous planar PIV-PLIF (particle image velocimetry-planar laser-induced fluorescence) measurements, reveal that the characteristics of wake and flow above porous patches are heavily influenced by $\phi$. In particular, we observe that the horizontal and vertical extent of the wake and scalar concentration downstream of the patches decreases and increases with $\phi$, respectively. Here, the recirculation bubble is shifted closer to the trailing edge of the patches as $\phi$ increases, limiting the flow from convecting downstream, decreasing the scalar concentration and virtually `extending' the patch in streamwise direction. As the bubble forms in the trailing edge, vertical bleeding increases and hence the concentration increases above patch where the cylinders appear to be `extend' vertically towards the freestream.


title: BERT Attends the Conversation: Improving Low-Resource Conversational ASR

abstract: We propose new, data-efficient training tasks for BERT models that improve performance of automatic speech recognition (ASR) systems on conversational speech. We include past conversational context and fine-tune BERT on transcript disambiguation without external data to rescore ASR candidates. Our results show word error rate recoveries up to 37.2%. We test our methods in low-resource data domains, both in language (Norwegian), tone (spontaneous, conversational), and topics (parliament proceedings and customer service phone calls). These techniques are applicable to any ASR system and do not require any additional data, provided a pre-trained BERT model. We also show how the performance of our context-augmented rescoring methods strongly depends on the degree of spontaneity and nature of the conversation.


title: Investigation of unsteady secondary flows and large-scale turbulence in heterogeneous turbulent boundary layers

abstract: Following the findings in \cite{wangsawijaya2020}, we re-examine the turbulent boundary layers developing over surfaces with spanwise heterogeneous roughness of various roughness wavelengths $0.32 \leq S/\overline{\delta} \leq 3.63$, where $S$ is the width of the roughness strips and $\overline{\delta}$ is the spanwise-averaged boundary-layer thickness. The heterogeneous cases induce counter-rotating secondary flows, and these are compared to the large-scale turbulent structures that occur naturally over the smooth wall. Both appear as meandering elongated high- and low-momentum streaks in the instantaneous flow field. Results suggest that the secondary flows might be spanwise-locked turbulent structures, with $S/\overline{\delta}$ governing the strength of the turbulent structures and possibly the efficacy of the surface in locking the structures in place (most effective when $S/\overline{\delta} \approx 1$). Conditional averages of the fluctuating velocity fields of both spanwise heterogeneous and smooth wall cases result in structures that are strongly reminiscent of the streak-vortex instability model. Secondary flows and large-scale structures coexist in the limits where either $S/\overline{\delta} \gg 1$ or $S/\overline{\delta} \ll 1$, where the secondary flows scale on $\delta$ or $S$, respectively. When $S/\overline{\delta} \gg 1$, the secondary flows are locked about the roughness transition, while relatively unaltered large-scale structures occur further from the transition. In the case where $S/\overline{\delta} \ll 1$, $S$-scaled secondary flows are confined close to the surface, coexisting with unaltered larger scale turbulent structures that penetrate much deeper into the layer.


title: Raining Rocks: An analytical formulation for collision timescales in planetary systems

abstract: The dynamical interaction of minor bodies (such as comets or asteroids) with planets plays an essential role in the planetary system's architecture and evolution. As a result of these interactions, structures like the Kuiper belt and the Oort cloud can be created. In particular, the collision of minor bodies with planets can drastically change the planet's internal and orbital evolution. We present an analytic formulation to determine the collision timescale for a minor body to impact a planet for arbitrary geometry. By comparing with a suite of detailed N-body simulations and an analytical method for collision timescales in the solar system, we confirmed the accuracy of our analytic formulation. As a proof of concept, we focused on the collision rate of minor bodies randomly distributed around a Jupiter-like planet, emulating a Kuiper belt-like disk. We show that our analytical method yields in good agreement with the numerical simulations. The formalism presented here thus provides a succinct and accurate alternative to numerical calculations.


title: Transformer Assisted Convolutional Network for Cell Instance Segmentation

abstract: Region proposal based methods like R-CNN and Faster R-CNN models have proven to be extremely successful in object detection and segmentation tasks. Recently, Transformers have also gained popularity in the domain of Computer Vision, and are being utilised to improve the performance of conventional models. In this paper, we present a relatively new transformer based approach to enhance the performance of the conventional convolutional feature extractor in the existing region proposal based methods. Our approach merges the convolutional feature maps with transformer-based token embeddings by applying a projection operation similar to self-attention in transformers. The results of our experiments show that transformer assisted feature extractor achieves a significant improvement in mIoU (mean Intersection over Union) scores compared to vanilla convolutional backbone.


title: Networked Time Series Prediction with Incomplete Data via Generative Adversarial Network

abstract: A networked time series (NETS) is a family of time series on a given graph, one for each node. It has a wide range of applications from intelligent transportation, environment monitoring to smart grid management. An important task in such applications is to predict the future values of a NETS based on its historical values and the underlying graph. Most existing methods require complete data for training. However, in real-world scenarios, it is not uncommon to have missing data due to sensor malfunction, incomplete sensing coverage, etc. In this paper, we study the problem of NETS prediction with incomplete data. We propose NETS-ImpGAN, a novel deep learning framework that can be trained on incomplete data with missing values in both history and future. Furthermore, we propose Graph Temporal Attention Networks, which incorporate the attention mechanism to capture both inter-time series and temporal correlations. We conduct extensive experiments on four real-world datasets under different missing patterns and missing rates. The experimental results show that NETS-ImpGAN outperforms existing methods, reducing the MAE by up to 25%.


title: An Analytic Treatment of Underdamped Axionic Blue Isocurvature Perturbations

abstract: Previous computations of strongly blue tilted axionic isocurvature spectra were computed in the parametric region in which the lightest time-dependent mass is smaller than the Hubble expansion rate during inflation, leading to an overdamped time evolution. Here we present the strongly blue tilted axionic isocurvature spectrum in an underdamped time evolution parametric regime. Somewhat surprisingly, there exist parametric regions with a strong resonant spectral behavior that leads to a rich isocurvature spectral shape. We focus on computing this resonant spectrum analytically in a large parametric region amenable to such computations. Because the spectrum is sensitive to nonperturbative classical field dynamics, a wide variety of analytic techniques are used including a time-space effective potential obtained by integrating out high frequency fluctuations.


title: Bilevel Imaging Learning Problems as Mathematical Programs with Complementarity Constraints: Reformulation and Theory

abstract: We investigate a family of bilevel imaging learning problems where the lower-level instance corresponds to a convex variational model involving first- and second-order nonsmooth sparsity-based regularizers. By using geometric properties of the primal-dual reformulation of the lower-level problem and introducing suitable auxiliar variables, we are able to reformulate the original bilevel problems as Mathematical Programs with Complementarity Constraints (MPCC). For the latter, we prove tight constraint qualification conditions (MPCC-RCPLD and partial MPCC-LICQ) and derive Mordukhovich (M-) and Strong (S-) stationarity conditions. The stationarity systems for the MPCC turn also into stationarity conditions for the original formulation. Second-order sufficient optimality conditions are derived as well, together with a local uniqueness result for stationary points. The proposed reformulation may be extended to problems in function spaces, leading to MPCC's with constraints on the gradient of the state. The MPCC reformulation also leads to the efficient use of available large-scale nonlinear programming solvers, as shown in a companion paper, where different imaging applications are studied.


title: Design of Spiral-Cable Forearm Exoskeleton to Assist Supination for Hemiparetic Stroke Subjects

abstract: We present the development of a cable-based passive forearm exoskeleton that is designed to assist supination for hemiparetic stroke survivors. Our device uniquely provides torque sufficient for counteracting spasticity within a below-elbow apparatus. The mechanism consists of a spiral single-tendon routing embedded in a rigid forearm brace and terminated at the hand and upper-forearm. A spool with an internal releasable-ratchet mechanism allows the user to manually retract the tendon and rotate the hand to counteract involuntary pronation synergies due to stroke. We characterize the mechanism with benchtop testing and five healthy subjects, and perform a preliminary assessment of the exoskeleton with a single chronic stroke subject having minimal supination ability. The mechanism can be integrated into an existing active hand-opening orthosis to enable supination support during grasping tasks, and also allows for a future actuated supination strategy.


title: Uniformly frustrated XY model: strengthening of the vortex lattice by intrinsic disorder

abstract: In superconducting films, the role of intrinsic disorder is typically to compete with superconductivity by fragmenting the global phase coherence and lowering the superfluid density. Nonetheless, when a transverse magnetic field is applied to the system and an Abrikosov vortex lattice forms, the presence of disorder can actually strengthen the superconducting state against thermal fluctuations. By means of Monte Carlo simulations on the uniformly frustrated XY model in two dimensions, we show that while for weak pinning the superconducting critical temperature $T_c$ increases with the applied field $H$, for strong enough pinning the experimental decreasing dependence between $T_c$ and $H$ is recovered with a resulting more robust vortex lattice.


title: SeanNet: Semantic Understanding Network for Localization Under Object Dynamics

abstract: We aim for domestic robots to perform long-term indoor service. Under the object-level scene dynamics induced by daily human activities, a robot needs to robustly localize itself in the environment subject to scene uncertainties. Previous works have addressed visual-based localization in static environments, yet the object-level scene dynamics challenge existing methods for the long-term deployment of the robot. This paper proposes a SEmantic understANding Network (SeanNet) architecture that enables an effective learning process with coupled visual and semantic inputs. With a dataset that contains object dynamics, we propose a cascaded contrastive learning scheme to train the SeanNet for learning a vector scene embedding. Subsequently, we can measure the similarity between the current observed scene and the target scene, whereby enables robust localization under object-level dynamics. In our experiments, we benchmark SeanNet against state-of-the-art image-encoding networks (baselines) on scene similarity measures. The SeanNet architecture with the proposed training method can achieve an 85.02\% accuracy which is higher than baselines. We further integrate the SeanNet and the other networks as the localizers into a visual navigation application. We demonstrate that SeanNet achieves higher success rates compared to the baselines.


title: Scaling up instance annotation via label propagation

abstract: Manually annotating object segmentation masks is very time-consuming. While interactive segmentation methods offer a more efficient alternative, they become unaffordable at a large scale because the cost grows linearly with the number of annotated masks. In this paper, we propose a highly efficient annotation scheme for building large datasets with object segmentation masks. At a large scale, images contain many object instances with similar appearance. We exploit these similarities by using hierarchical clustering on mask predictions made by a segmentation model. We propose a scheme that efficiently searches through the hierarchy of clusters and selects which clusters to annotate. Humans manually verify only a few masks per cluster, and the labels are propagated to the whole cluster. Through a large-scale experiment to populate 1M unlabeled images with object segmentation masks for 80 object classes, we show that (1) we obtain 1M object segmentation masks with an total annotation time of only 290 hours; (2) we reduce annotation time by 76x compared to manual annotation; (3) the segmentation quality of our masks is on par with those from manually annotated datasets. Code, data, and models are available online.


title: A note on inflation in dRGT massive gravity

abstract: Although the dRGT massive gravity successfully explains the late-time cosmic acceleration, it cannot justify inflation. On the other hand, and in the frameworks of General Relativity and modified gravity, the interests and attempts to describe dark energy and inflation by using Lagranginas, which may have pole, have recently been enhanced. Subsequently, we are going to show that this kind of Lagrangian may justify inflation in the framework of dRGT massive gravity. The study is done focusing on the power and exponential potentials, and the results show a plausible consistency with the Planck 2018 data and its combination with BK18 and BAO.


title: Turing approximations, toric isometric embeddings & manifold convolutions

abstract: Convolutions are fundamental elements in deep learning architectures. Here, we present a theoretical framework for combining extrinsic and intrinsic approaches to manifold convolution through isometric embeddings into tori. In this way, we define a convolution operator for a manifold of arbitrary topology and dimension. We also explain geometric and topological conditions that make some local definitions of convolutions which rely on translating filters along geodesic paths on a manifold, computationally intractable. A result of Alan Turing from 1938 underscores the need for such a toric isometric embedding approach to achieve a global definition of convolution on computable, finite metric space approximations to a smooth manifold.


title: Distributed Privacy-Preserving Electric Vehicle Charging Control Based on Secret Sharing

abstract: Cooperative electric vehicle (EV) charging control has emerged as a key component in grid-edge resource (GER) management. However, customers' privacy remains a major barrier to large-scale implementation of EV charging control. In this paper, we develop a distributed privacy-preserving EV charging control protocol based on secret sharing (SS) that 1) achieves scalability over EV population size; 2) enjoys higher computation efficiency compared to homomorphic encryption (HE) based methods; and 3) secures the privacy of the participating EVs against honest-but-curious adversaries and external eavesdroppers. The cooperative EV charging control is formulated to achieve overnight valley-filling and framed into the projected gradient algorithm (PGA) structure as a distributed optimization problem. SS is integrated into PGA to achieve secure updates of both primal and dual variables. Theoretical security analyses and simulations in a residential area are conducted to prove the privacy preservation guarantee as well as the efficacy and efficiency of the proposed privacy preservation method. Broadly, the proposed method can be readily extended to various GER control applications.


title: A Rate Splitting Strategy for Uplink CR-NOMA Systems

abstract: In uplink non-orthogonal multiple access (NOMA) channels, the existing cooperative successive interference cancellation (SIC) and power control (PC) schemes lack the capability of achieving the full capacity region, which restricts the outage performance of uplink NOMA users. For the uplink cognitive radio (CR) inspired NOMA system, we propose a new rate splitting (RS) strategy to maximize the achievable rate of the secondary user without deteriorating the primary user's outage performance. Based on the interference threshold and its own channel gain, the secondary user adaptively conducts RS, transmit power allocation and SIC, which utilizes the transmit power efficiently. Closed-form expression of the outage probability is derived for the secondary user. Numerical results show that the proposed RS scheme achieves the best outage performance for the secondary user among the existing cooperative SIC and PC schemes.


title: First beams at Neutrons For Science

abstract: The neutrons for science facility (NFS), the first operational experimental area of the new GANIL/SPIRAL-2 facility, received its first beams in December 2019. Proton-induced reaction cross-sections as well as neutron beam characteristics were measured during the first commissioning phases. The first results, showing the features of the facility, are presented here and compared with previously published data. The physics cases and the first accepted experiments are presented as well.


title: Co-training an Unsupervised Constituency Parser with Weak Supervision

abstract: We introduce a method for unsupervised parsing that relies on bootstrapping classifiers to identify if a node dominates a specific span in a sentence. There are two types of classifiers, an inside classifier that acts on a span, and an outside classifier that acts on everything outside of a given span. Through self-training and co-training with the two classifiers, we show that the interplay between them helps improve the accuracy of both, and as a result, effectively parse. A seed bootstrapping technique prepares the data to train these classifiers. Our analyses further validate that such an approach in conjunction with weak supervision using prior branching knowledge of a known language (left/right-branching) and minimal heuristics injects strong inductive bias into the parser, achieving 63.1 F$_1$ on the English (PTB) test set. In addition, we show the effectiveness of our architecture by evaluating on treebanks for Chinese (CTB) and Japanese (KTB) and achieve new state-of-the-art results. Our code and pre-trained models are available at https://github.com/Nickil21/weakly-supervised-parsing.


title: An explorative study on how human-robot interaction is taken into account by robot developers in praxis

abstract: How is human-robot interaction considered within the development of new robotic systems by practitioners? This study sets out to inquire, whether the development teams of robotic products have been considering human factor methods in their design and implementation process. We were specifically interested in the non-verbal communication methods they were aiming to implement, and how they have approached the design process for these. Although valuable insights on tasks and communication needs during the different phases of robot operation could be gathered, the results of this study indicate, that the perspective of the human user or bystander is very often neglected and that knowledge on methods for engineering human-robot interaction is missing. The study was conducted with eleven development teams consisting of robot manufacturers and students within a robot building course representing overall 68 individual participants.


title: Modelling of the Fender Bassman 5F6-A Tone Stack

abstract: This paper outlines the procedure for the effective modelling of a complex analogue filter circuit. The Fender Bassman 5F6-A is a circuit commonly employed in guitar amplifiers to shape the tonal characteristics of the amplifier output. On first inspection this circuit may look rather simple, however the controls are not orthogonal, resulting in complicated filter coefficients as the controls are varied. This in turn can make the circuit difficult to analyse without the use of mathematical emulation tools such as PSPICE or MATLAB. First the circuit is described, a method of analysis is proposed and general expressions for continuous-time coefficients are given. A MATLAB model is then produced and the frequency responses of which are shown.


title: Flavor isospin waves in one-dimensional axisymmetric neutrino gases

abstract: Flavor oscillations can occur on very short spatial and temporal scales in the dense neutrino media in a core-collapse supernova (CCSN) or binary neutron star merger (BNSM) event. Although the dispersion relations (DRs) of the fast neutrino oscillations can be obtained by linearizing the equations of motion (EoM) before the emergence of any significant flavor conversion, one largely depends on numerical calculations to understand this interesting phenomenon in the nonlinear regime. In this work we demonstrate that there exist nontrivial solutions to the flavor EoM that govern the fast oscillations in one-dimensional axisymmetric neutrino gases. These solutions represent a coherent flavor isospin wave similar to the magnetic spin wave in a lattice of magnetic dipoles. We also compute the DRs of such waves in some example cases which are closely related to the DRs of the fast neutrino oscillations obtained in the linear regime. This result sheds new light on the long-term behavior of fast neutrino oscillations which can have various implications for the CCSN and BNSM events.


title: $BC_2$ type multivariable matrix functions and matrix spherical functions

abstract: Matrix spherical functions associated to the compact symmetric pair $(\mathrm{SU}(m+2), \mathrm{S}(\mathrm{U}(2)\times \mathrm{U}(m))$, having reduced root system of type $\mathrm{BC}_2$, are studied. We consider an irreducible $K$-representation $(\pi,V)$ arising from the $\mathrm{U}(2)$-part of $K$, and the induced representation $\mathrm{Ind}_K^G \pi$ splits multiplicity free. The corresponding spherical functions, i.e. $\Phi \colon G \to \mathrm{End}(V)$ satisfying $\Phi(k_1gk_2)=\pi(k_1)\Phi(g)\pi(k_2)$ for all $g\in G$, $k_1,k_2\in K$, are studied by studying certain leading coefficients which involve hypergeometric functions. This is done explicitly using the action of the radial part of the Casimir operator on these functions and their leading coefficients. To suitably grouped matrix spherical functions we associate two-variable matrix orthogonal polynomials giving a matrix analogue of Koornwinder's 1970s two-variable orthogonal polynomials, which are Heckman-Opdam polynomials for $\mathrm{BC}_2$. In particular, we find explicit orthogonality relations and the matrix polynomials being eigenfunctions to an explicit second order matrix partial differential operator. The scalar part of the matrix weight is less general than Koornwinder's weight.


title: Evolutionary Algorithms for Solving Unconstrained, Constrained and Multi-objective Noisy Combinatorial Optimisation Problems

abstract: We present an empirical study of a range of evolutionary algorithms applied to various noisy combinatorial optimisation problems. There are three sets of experiments. The first looks at several toy problems, such as OneMax and other linear problems. We find that UMDA and the Paired-Crossover Evolutionary Algorithm (PCEA) are the only ones able to cope robustly with noise, within a reasonable fixed time budget. In the second stage, UMDA and PCEA are then tested on more complex noisy problems: SubsetSum, Knapsack and SetCover. Both perform well under increasing levels of noise, with UMDA being the better of the two. In the third stage, we consider two noisy multi-objective problems (CountingOnesCountingZeros and a multi-objective formulation of SetCover). We compare several adaptations of UMDA for multi-objective problems with the Simple Evolutionary Multi-objective Optimiser (SEMO) and NSGA-II. We conclude that UMDA, and its variants, can be highly effective on a variety of noisy combinatorial optimisation, outperforming many other evolutionary algorithms.


title: An approximate expectation-maximization for two-dimensional multi-target detection

abstract: We consider the two-dimensional multi-target detection (MTD) problem of estimating a target image from a noisy measurement that contains multiple copies of the image, each randomly rotated and translated. The MTD model serves as a mathematical abstraction of the structure reconstruction problem in single-particle cryo-electron microscopy, the chief motivation of this study. We focus on high noise regimes, where accurate detection of image occurrences within a measurement is impossible. To estimate the image, we develop an expectation-maximization framework that aims to maximize an approximation of the likelihood function. We demonstrate image recovery in highly noisy environments, and show that our framework outperforms the previously studied autocorrelation analysis in a wide range of parameters.


title: Error metric for non-trace-preserving quantum operations

abstract: We study the problem of measuring errors in non-trace-preserving quantum operations, with a focus on their impact on quantum computing. We propose an error metric that efficiently provides an upper bound on the trace distance between the normalized output states from imperfect and ideal operations, while remaining compatible with the diamond distance. As a demonstration of its application, we apply our metric in the analysis of a lossy beam splitter and a nondeterministic conditional sign-flip gate, two primary non-trace-preserving operations in the Knill-Laflamme-Milburn protocol. We then turn to the leakage errors of neutral-atom quantum computers, finding that these errors scale worse than previously anticipated, implying a more stringent fault-tolerant threshold. We also assess the quantum Zeno gate's error using our metric. In a broader context, we discuss the potential of our metric to analyze general postselected protocols, where it can be employed to study error propagation and estimate thresholds in fault-tolerant quantum computing. The results highlight the critical role of our proposed error metric in understanding and addressing challenges in practical quantum information processing.


title: FedDQ: Communication-Efficient Federated Learning with Descending Quantization

abstract: Federated learning (FL) is an emerging learning paradigm without violating users' privacy. However, large model size and frequent model aggregation cause serious communication bottleneck for FL. To reduce the communication volume, techniques such as model compression and quantization have been proposed. Besides the fixed-bit quantization, existing adaptive quantization schemes use ascending-trend quantization, where the quantization level increases with the training stages. In this paper, we first investigate the impact of quantization on model convergence, and show that the optimal quantization level is directly related to the range of the model updates. Given the model is supposed to converge with the progress of the training, the range of the model updates will gradually shrink, indicating that the quantization level should decrease with the training stages. Based on the theoretical analysis, a descending quantization scheme named FedDQ is proposed. Experimental results show that the proposed descending quantization scheme can save up to 65.2% of the communicated bit volume and up to 68% of the communication rounds, when compared with existing schemes.


title: On relation between statistical ideal and ideal generated by a modulus function

abstract: We characterize those modulus functions $f$ for which the ideal generated by $f$ is equal to the statistical ideal.


title: Sources in the Weyl double copy

abstract: The Weyl double copy relates exact solutions in general relativity to exact solutions in gauge theory, formulated in the spinorial language. To date, the Weyl double copy is understood and employed only for vacuum spacetimes, and hence only to vacuum gauge theories. In this letter, we propose an extension to the Weyl double copy that provides a systematic procedure for treating gravitational sources. We show this extended Weyl double copy gives a new perspective on the Kerr-Newman black hole and the general class of Petrov type D electro-vac spacetimes.


title: Dynamic clustering of passive colloids in dense suspensions of motile bacteria

abstract: Mixtures of active and passive particles are predicted to exhibit a variety of nonequilibrium phases. Here we report a dynamic clustering phase in mixtures of colloids and motile bacteria. We show that colloidal clustering results from a balance between bond breaking due to persistent active motion and bond stabilization due to torques that align active particle velocity tangentially to the passive particle surface. Furthermore, dynamic clustering spans a broad regime between diffusivity-based and motility-induced phase separation that subsumes typical bacterial motility parameters.


title: Kundt wave geometries in Eddington-inspired Born-Infeld gravity: New solutions and memory effects

abstract: We explore memory effects for novel Kundt wave spacetimes in the recently proposed Eddington-inspired Born-Infeld (EiBI) theory of gravity. First, we construct new, exact Kundt wave geometries in this theory for two different matter sources--(i) generic matter designed to satisfy the field equations as well as energy conditions, and (ii) electromagnetic field. For both sources we find that the EiBI theory parameter $\kappa$ couples only with the nonradiative part of the physical metric solution. Thereafter, we solve the geodesic and the geodesic deviation equations, in the above spacetimes with the aim of arriving at memory effects. This analysis is carried out numerically and reveals unique memory features depending on the type of matter source present and the signature of the spacetime scalar curvature. The role of $\kappa$ in influencing the memory effect, for a given background spacetime, is also noted. Thus, apart from providing novel radiative solutions in EiBI gravity, we also show how different matter configurations are responsible for distinct memory characteristics.


title: On the Correspondence between Gaussian Processes and Geometric Harmonics

abstract: We discuss the correspondence between Gaussian process regression and Geometric Harmonics, two similar kernel-based methods that are typically used in different contexts. Research communities surrounding the two concepts often pursue different goals. Results from both camps can be successfully combined, providing alternative interpretations of uncertainty in terms of error estimation, or leading towards accelerated Bayesian Optimization due to dimensionality reduction.


title: Robustness modularity in complex networks

abstract: A basic question in network community detection is how modular a given network is. This is usually addressed by evaluating the quality of partitions detected in the network. The Girvan-Newman (GN) modularity function is the standard way to make this assessment, but it has a number of drawbacks. Most importantly, it is not clearly interpretable, given that the measure can take relatively large values on partitions of random networks without communities. Here we propose a new measure based on the concept of robustness: modularity is the probability to find trivial partitions when the structure of the network is randomly perturbed. This concept can be implemented for any clustering algorithm capable of telling when a group structure is absent. Tests on artificial and real graphs reveal that robustness modularity can be used to assess and compare the strength of the community structure of different networks. We also introduce two other quality functions: modularity difference, a suitably normalized version of the GN modularity; information modularity, a measure of distance based on information compression. Both measures are strongly correlated with robustness modularity, and are promising options as well.

categories: physics.soc-ph cs.LG cs.SI physics.data-an


title: Tradeoffs in Hierarchical Voting Systems

abstract: Condorcet's jury theorem states that the correct outcome is reached in direct majority voting systems with sufficiently large electorates as long as each voter's independent probability of voting for that outcome is greater than 0.5. Yet, in situations where direct voting systems are infeasible, such as due to high implementation and infrastructure costs, hierarchical voting systems provide a reasonable alternative. We study differences in outcome precision between hierarchical and direct voting systems for varying group sizes, abstention rates, and voter competencies. Using asymptotic expansions of the derivative of the reliability function (or Banzhaf number), we first prove that indirect systems differ most from their direct counterparts when group size and number are equal to each other, and therefore to $\sqrt{N_{\rm d}}$, where $N_{\rm d}$ is the total number of voters in the direct system. In multitier systems, we prove that this difference is maximized when group size equals $\sqrt[n]{N_{\rm d}}$, where $n$ is the number of hierarchical levels. Second, we show that while direct majority rule always outperforms hierarchical voting for homogeneous electorates that vote with certainty, as group numbers and size increase, hierarchical majority voting gains in its ability to represent all eligible voters. Furthermore, when voter abstention and competency are correlated within groups, hierarchical systems often outperform direct voting, which we show by using a generating function approach that is able to analytically characterize heterogeneous voting systems.


title: Twisted two-dimensional material stacks for polarization optics

abstract: The ability to control light polarization state is critically important for diverse applications in information processing, telecommunications, and spectroscopy. Here, we propose that a stack of anisotropic van der Waals materials can facilitate the building of optical elements with Jones matrices of unitary, Hermitian, non-normal, singular, degenerate, and defective classes. We show that the twisted stack with electrostatic control can function as arbitrary-birefringent wave-plate or arbitrary polarizer with tunable degree of non-normality, which in turn give access to plethora of polarization transformers including rotators, pseudorotators, symmetric and ambidextrous polarizers. Moreover, we discuss an electrostatic-reconfigurable stack which can be tuned to operate as four different polarizers and be used for Stokes polarimetry.


title: Complexity of Traveling Tournament Problem with Trip Length More Than Three

abstract: The Traveling Tournament Problem is a sports-scheduling problem where the goal is to minimize the total travel distance of teams playing a double round-robin tournament. The constraint 'k' is an imposed upper bound on the number of consecutive home or away matches. It is known that TTP is NP-Hard for k=3 as well as k=infinity. In this work, the general case has been settled by proving that TTP-k is NP-Complete for any fixed k>3.


title: Wronskians, total positivity, and real Schubert calculus

abstract: A complete flag in $\mathbb{R}^n$ is a sequence of nested subspaces $V_1 \subset \cdots \subset V_{n-1}$ such that each $V_k$ has dimension $k$. It is called totally nonnegative if all its Pl\"ucker coordinates are nonnegative. We may view each $V_k$ as a subspace of polynomials in $\mathbb{R}[x]$ of degree at most $n-1$, by associating a vector $(a_1, \dots, a_n)$ in $\mathbb{R}^n$ to the polynomial $a_1 + a_2x + \cdots + a_nx^{n-1}$. We show that a complete flag is totally nonnegative if and only if each of its Wronskian polynomials $\mathsf{Wr}(V_k)$ is nonzero on the interval $(0, \infty)$. In the language of Chebyshev systems, this means that the flag forms a Markov system or $ECT$-system on $(0, \infty)$. This gives a new characterization and membership test for the totally nonnegative flag variety. Similarly, we show that a complete flag is totally positive if and only if each $\mathsf{Wr}(V_k)$ is nonzero on $[0, \infty]$. We use these results to show that a conjecture of Eremenko (2015) in real Schubert calculus is equivalent to the following conjecture: if $V$ is a finite-dimensional subspace of polynomials such that all complex zeros of $\mathsf{Wr}(V)$ lie in the interval $(-\infty, 0)$, then all Pl\"ucker coordinates of $V$ are real and positive. This conjecture is a totally positive strengthening of a result of Mukhin, Tarasov, and Varchenko (2009), and can be reformulated as saying that all complex solutions to a certain family of Schubert problems in the Grassmannian are real and totally positive. We also show that our conjecture is equivalent to a totally positive version of the secant conjecture of Sottile (2003).


title: Processing, electronic properties and applications of Carbon-based nanocomposites

abstract: The last two decades, in particular, have witnessed a large volume of research revolving around structure-property correlation in Carbon based nanocomposites, synthesized by several methods.In the simplest of terms, the electronic properties of these nanomaterials, which form the present context of discussion, vary mainly as a function of three parameters, out of which 2 are process parameters (viz. the kind of reinforcement and method of synthesis) and one is a structure-dependent parameter which is highly influenced by the two aforementioned process parameters and plays a vital role in determining the ionic and electronic transport phenomenon in these materials. To be more specific in the context of structure dependent parameter, the interaction between electrons and the equilibrium 0D (point) defects, along with different types of interfaces, plays a crucial role in the understanding of electronic properties, apart from physical and chemical properties in these materials. These materials find applications in a large number of sectors, such as nanoelectronics, supercapacitors, absorbent, photocatalytic, photovoltaic, solar cells, sensors, and energy storage due to their exceptional electronic, mechanical and chemical properties. The present chapter will intend to focus primarily on providing a brief overview of the state-of-the-art on research along with detailed discussions on some recent developments in determining and understanding electronic properties of some conventional Carbon-based nanocomposites (synthesized by several different techniques) based on the structure property correlation of these materials and finally, address some of the significant challenges in this field.


title: On pseudobounded and premeager paratopological groups

abstract: Let $G$ be a paratopological group. Following F. Lin and S. Lin, we say that the group $G$ is pseudobounded, if for any neighborhood $U$ of the identity of $G$, there exists a natural number $n$ such that $U^n=G$. The group $G$ is $\omega$-pseudobounded, if for any neighborhood $U$ of the identity of $G$, the group $G$ is a union of sets $U^n$, where $n$ is a natural number. The group $G$ is premeager, if $G\ne N^n$ for any nowhere dense subset $N$ of $G$ and any positive integer $n$. In this paper we investigate relations between the above classes of groups and answer some questions posed by F. Lin, S. Lin, and S\'anchez.


title: You Only Evaluate Once: a Simple Baseline Algorithm for Offline RL

abstract: The goal of offline reinforcement learning (RL) is to find an optimal policy given prerecorded trajectories. Many current approaches customize existing off-policy RL algorithms, especially actor-critic algorithms in which policy evaluation and improvement are iterated. However, the convergence of such approaches is not guaranteed due to the use of complex non-linear function approximation and an intertwined optimization process. By contrast, we propose a simple baseline algorithm for offline RL that only performs the policy evaluation step once so that the algorithm does not require complex stabilization schemes. Since the proposed algorithm is not likely to converge to an optimal policy, it is an appropriate baseline for actor-critic algorithms that ought to be outperformed if there is indeed value in iterative optimization in the offline setting. Surprisingly, we empirically find that the proposed algorithm exhibits competitive and sometimes even state-of-the-art performance in a subset of the D4RL offline RL benchmark. This result suggests that future work is needed to fully exploit the potential advantages of iterative optimization in order to justify the reduced stability of such methods.


title: Oxide two-dimensional electron gas with high mobility at room-temperature

abstract: The prospect of 2-dimensional electron gases (2DEGs) possessing high mobility at room temperature in wide-bandgap perovskite stannates is enticing for oxide electronics, particularly to realize transparent and high-electron mobility transistors. Nonetheless only a small number of studies to date report 2DEGs in BaSnO3-based heterostructures. Here, we report 2DEG formation at the LaScO3/BaSnO3 (LSO/BSO) interface with a room-temperature mobility of 60 cm2/V s at a carrier concentration of 1.7x1013 cm-2. This is an order of magnitude higher mobility at room temperature than achieved in SrTiO3-based 2DEGs. We achieved this by combining a thick BSO buffer layer with an ex-situ high-temperature treatment, which not only reduces the dislocation density but also produces a SnO2-terminated atomically flat surface, followed by the growth of an overlying BSO/LSO interface. Using weak-beam dark field imaging and in-line electron holography technique, we reveal a reduction of the threading dislocation density, and provide direct evidence for the spatial confinement of a 2DEG at the BSO/LSO interface. Our work opens a new pathway to explore the exciting physics of stannate-based 2DEGs at application-relevant temperatures for oxide nanoelectronics.


title: Stability of a flexible missile described by asymptotics of the eigenvalues of fourth order boundary value problems

abstract: Fourth order problems, with the differential equation $y^{(4)}-(gy')'=\lambda^2y$, where $g\in C^1[0,a]$ and $a>0$, occur in engineering on stability of elastic rods. They occur as well in aeronautics to describe the stability of a flexible missile. Fourth order Birkhoff regular problems with the differential equation $y^{(4)}-(gy')'=\lambda^2y$ and eigenvalue dependent boundary conditions are considered. These problems have quadratic operator representations with non self-adjoint operators. The first four terms of the asymptotics of the eigenvalues of the problems as well as those of the eigenvalues of the problem describing the stability of a flexible missile are evaluated explicitly.


title: Coarsening Optimization for Differentiable Programming

abstract: This paper presents a novel optimization for differentiable programming named coarsening optimization. It offers a systematic way to synergize symbolic differentiation and algorithmic differentiation (AD). Through it, the granularity of the computations differentiated by each step in AD can become much larger than a single operation, and hence lead to much reduced runtime computations and data allocations in AD. To circumvent the difficulties that control flow creates to symbolic differentiation in coarsening, this work introduces phi-calculus, a novel method to allow symbolic reasoning and differentiation of computations that involve branches and loops. It further avoids "expression swell" in symbolic differentiation and balance reuse and coarsening through the design of reuse-centric segment of interest identification. Experiments on a collection of real-world applications show that coarsening optimization is effective in speeding up AD, producing several times to two orders of magnitude speedups.


title: The search for gas in debris discs: ALMA detection of CO gas in HD 36546

abstract: Debris discs represent the last stages of planet formation and as such are expected to be depleted of primordial gas. Nonetheless, in the last few years the presence of cold gas has been reported in $\sim$ 20 debris discs from far-IR to (sub-)mm observations and hot gas has been observed in the optical spectra of debris discs for decades. While the origin of this gas is still uncertain, most evidences point towards a secondary origin, as a result of collisions and evaporation of small bodies in the disc. In this paper, we present ALMA observations aimed at the detection of CO gas in a sample of 8 debris discs with optical gas detections. We report the detection of CO ($^{12}$CO and $^{13}$CO) gas in HD 36546, the brightest and youngest disc in our sample, and provide upper limits to the presence of gas in the remaining seven discs.


title: On the toroidal-velocity anti-dynamo theorem under the presence of nonuniform electric conductivity

abstract: Laminar electrically conducting Couette flows with the hydrodynamically stable quasi-Keplerian rotation profile and non-uniform conductivity are probed for dynamo instability. In spherical geometry the equations for the poloidal and the toroidal field components completely decouple, resulting in free decay, regardless of the spatial distribution of the electric conductivity. In cylindrical geometry the poloidal and toroidal components do not decouple, but here also we do not find dynamo excitations for the cases that the electric conductivity only depends on the radius or -- much more complex -- that it only depends on the azimuthal or the axial coordinate. The transformation of the plane-flow dynamo model of Busse \& Wicht (1992) to cylindrical or spherical geometry therefore fails. It is also shown that even the inclusion of axial flows of both directions does {\em not} support the dynamo mechanism. The Elsasser toroidal-velocity antidynamo theorem, according to which dynamos without any radial velocity component cannot work, is thus not softened by non-uniform conductivity distributions.


title: Stability of uniformly eventually positive $C_0$-semigroups on $L_p$-spaces

abstract: We give a short and elementary proof of the theorem of Lutz Weis that the growth bound of a positive $C_0$-semigroup on $L_p(\mu)$ equals the spectral bound of its generator. In addition, we generalise the result to the case of uniformly eventually positive semigroups.


title: COVID-19 India Dataset: Parsing COVID-19 Data in Daily Health Bulletins from States in India

abstract: While India has been one of the hotspots of COVID-19, data about the pandemic from the country has proved to be largely inaccessible at scale. Much of the data exists in unstructured form on the web, and limited aspects of such data are available through public APIs maintained manually through volunteer effort. This has proved to be difficult both in terms of ease of access to detailed data and with regards to the maintenance of manual data-keeping over time. This paper reports on our effort at automating the extraction of such data from public health bulletins with the help of a combination of classical PDF parsers and state-of-the-art machine learning techniques. In this paper, we will describe the automated data-extraction technique, the nature of the generated data, and exciting avenues of ongoing work.


title: Symplectic embeddings into disk cotangent bundles

abstract: In this paper, we compute the embedded contact homology (ECH) capacities of the disk cotangent bundles $D^*S^2$ and $D^*\mathbb{R} P^2$. We also find sharp symplectic embeddings into these domains. In particular, we compute their Gromov widths. In order to do that, we explicitly calculate the ECH chain complexes of $S^*S^2$ and $S^* \mathbb{R} P^2$ using a direct limit argument on the action inspired by Bourgeois's Morse-Bott approach and ideas from Nelson-Weiler's work on the ECH of prequantization bundles. Moreover, we use integrable systems techniques to find explicit symplectic embeddings. In particular, we prove that the disk cotangent bundles of a hemisphere and of a punctured sphere are symplectomorphic to an open ball and a symplectic bidisk, respectively.


title: Phoebe: A Learning-based Checkpoint Optimizer

abstract: Easy-to-use programming interfaces paired with cloud-scale processing engines have enabled big data system users to author arbitrarily complex analytical jobs over massive volumes of data. However, as the complexity and scale of analytical jobs increase, they encounter a number of unforeseen problems, hotspots with large intermediate data on temporary storage, longer job recovery time after failures, and worse query optimizer estimates being examples of issues that we are facing at Microsoft. To address these issues, we propose Phoebe, an efficient learning-based checkpoint optimizer. Given a set of constraints and an objective function at compile-time, Phoebe is able to determine the decomposition of job plans, and the optimal set of checkpoints to preserve their outputs to durable global storage. Phoebe consists of three machine learning predictors and one optimization module. For each stage of a job, Phoebe makes accurate predictions for: (1) the execution time, (2) the output size, and (3) the start/end time taking into account the inter-stage dependencies. Using these predictions, we formulate checkpoint optimization as an integer programming problem and propose a scalable heuristic algorithm that meets the latency requirement of the production environment. We demonstrate the effectiveness of Phoebe in production workloads, and show that we can free the temporary storage on hotspots by more than 70% and restart failed jobs 68% faster on average with minimum performance impact. Phoebe also illustrates that adding multiple sets of checkpoints is not cost-efficient, which dramatically reduces the complexity of the optimization.


title: Give Me a Few Hours: Exploring Short Timescales in Rubin Observatory Cadence Simulations

abstract: The limiting temporal resolution of a time-domain survey in detecting transient behavior is set by the time between observations of the same sky area. We analyze the distribution of visit separations for a range of Vera C. Rubin Observatory cadence simulations. Current simulations are strongly peaked at the 22 minute visit pair separation and provide effectively no constraint on temporal evolution within the night. This choice will necessarily prevent Rubin from discovering a wide range of astrophysical phenomena in time to trigger rapid followup. We present a science-agnostic metric to supplement detailed simulations of fast-evolving transients and variables and suggest potential approaches for improving the range of timescales explored.


title: Cryogenic characterization and modeling of a CMOS floating-gate device for quantum control hardware

abstract: We perform the characterization and modeling of a floating-gate device realized with a commercial 350-nm CMOS technology at cryogenic temperature. The programmability of the device offers a solution in the realization of a precise and flexible cryogenic system for qubits control in large-scale quantum computers. The device stores onto a floating-gate node a non-volatile charge, which can be bidirectionally modified by Fowler-Nordheim tunneling and impact-ionized hotelectron injection. These two injection mechanisms are characterized and modeled in compact equations both at 300 K and 15 K. At cryogenic temperature, we show a fine-tuning of the stored charge compatible with the operation of a precise analog memory. Moreover, we developed accurate simulation models of the proposed floating-gate device that set the stage for designing a programmable analog circuit with better performances and accuracy at a few Kelvin. This work offers a solution in the design of configurable analog electronics to be employed for accurately read out the qubit state at deep-cryogenic temperature.


title: Prediction of the Facial Growth Direction is Challenging

abstract: Facial dysmorphology or malocclusion is frequently associated with abnormal growth of the face. The ability to predict facial growth (FG) direction would allow clinicians to prepare individualized therapy to increase the chance for successful treatment. Prediction of FG direction is a novel problem in the machine learning (ML) domain. In this paper, we perform feature selection and point the attribute that plays a central role in the abovementioned problem. Then we successfully apply data augmentation (DA) methods and improve the previously reported classification accuracy by 2.81%. Finally, we present the results of two experienced clinicians that were asked to solve a similar task to ours and show how tough is solving this problem for human experts.


title: Cartesian dictionary-based native T1 and T2 mapping of the myocardium

abstract: Purpose: To implement and evaluate a new dictionary-based technique for native myocardial T1 and T2 mapping using Cartesian sampling. Methods: The proposed technique (Multimapping) consisted of single-shot Cartesian image acquisitions in 10 consecutive cardiac cycles, with inversion pulses in cycle 1 and 5, and T2 preparation (TE: 30ms, 50ms and 70ms) in cycles 8-10. Multimapping was simulated for different T1 and T2, where entries corresponding to the k-space centers were matched to acquired data. Experiments were performed in a phantom, 16 healthy subjects and three patients with cardiovascular disease. Results: Multimapping phantom measurements showed good agreement with reference values for both T1 and T2, with no discernable heart-rate dependency for T1 and T2 within the range of myocardium. In vivo mean T1 in healthy subjects was significantly higher using Multimapping (T1=1114+/-14ms) compared to the reference (T1=991+/-26ms) (p<0.01). Mean Multimapping T2 (47.1+/-1.3ms) and T2 spatial variability (5.8+/-1.0ms) was significantly lower compared to the reference (T2=54.7+/-2.2ms, p<0.001; spatial variability=8.4+/-2.0ms, p<0.01). Increased T1 and T2 was detected in all patients using Multimapping. Conclusions: Multimapping allows for simultaneous native myocardial T1 and T2 mapping with a conventional Cartesian trajectory, demonstrating promising in vivo image quality and parameter quantification results.


title: Approximate Message Passing for orthogonally invariant ensembles: Multivariate non-linearities and spectral initialization

abstract: We study a class of Approximate Message Passing (AMP) algorithms for symmetric and rectangular spiked random matrix models with orthogonally invariant noise. The AMP iterates have fixed dimension $K \geq 1$, a multivariate non-linearity is applied in each AMP iteration, and the algorithm is spectrally initialized with $K$ super-critical sample eigenvectors. We derive the forms of the Onsager debiasing coefficients and corresponding AMP state evolution, which depend on the free cumulants of the noise spectral distribution. This extends previous results for such models with $K=1$ and an independent initialization. Applying this approach to Bayesian principal components analysis, we introduce a Bayes-OAMP algorithm that uses as its non-linearity the posterior mean conditional on all preceding AMP iterates. We describe a practical implementation of this algorithm, where all debiasing and state evolution parameters are estimated from the observed data, and we illustrate the accuracy and stability of this approach in simulations.


title: Carroll symmetry, dark energy and inflation

abstract: Carroll symmetry arises from Poincar\'e symmetry upon taking the limit of vanishing speed of light. We determine the constraints on the energy-momentum tensor implied by Carroll symmetry and show that for energy-momentum tensors of perfect fluid form, these imply an equation of state ${\cal E}+P=0$ for energy density plus pressure. Therefore Carroll symmetry might be relevant for dark energy and inflation. In the Carroll limit, the Hubble radius goes to zero and outside it recessional velocities are naturally large compared to the speed of light. The de Sitter group of isometries, after the limit, becomes the conformal group in Euclidean flat space. We also study the Carroll limit of chaotic inflation, and show that the scalar field is naturally driven to have an equation of state with $w=-1$. Finally we show that the freeze-out of scalar perturbations in the two point function at horizon crossing is a consequence of Carroll symmetry. To make the paper self-contained, we include a brief pedagogical review of Carroll symmetry, Carroll particles and Carroll field theories that contains some new material as well. In particular we show, using an expansion around speed of light going to zero, that for scalar and Maxwell type theories one can take two different Carroll limits at the level of the action. In the Maxwell case these correspond to the electric and magnetic limit. For point particles we show that there are two types of Carroll particles: those that cannot move in space and particles that cannot stand still.


title: Characterizations of weakly $\mathcal{K}$-analytic and Va\v{s}\'ak spaces using projectional skeletons and separable PRI

abstract: We find characterizations of Va\v{s}\'ak spaces and weakly $\mathcal{K}$-analytic spaces using the notions of separable projectional resolution of the identity (SPRI) and of projectional skeleton. This in particular addresses a recent challenge suggested by M. Fabian and V. Montesinos in \cite{FM18}. Our method of proof also gives similar characterizations of WCG spaces and their subspaces (some aspects of which were known, some are new). Moreover we show that for countably many projectional skeletons $\{\mathfrak{s}_n: n \in \omega\}$ on a Banach space inducing the same set, there exists a projectional skeleton on the space (indexed by ranges of the corresponding projections) which is isomorphic to a subskeleton of each $\mathfrak{s}_n$, $n \in \omega$.


title: Enhancement of Anime Imaging Enlargement using Modified Super-Resolution CNN

abstract: Anime is a storytelling medium similar to movies and books. Anime images are a kind of artworks, which are almost entirely drawn by hand. Hence, reproducing existing Anime with larger sizes and higher quality images is expensive. Therefore, we proposed a model based on convolutional neural networks to extract outstanding features of images, enlarge those images, and enhance the quality of Anime images. We trained the model with a training set of 160 images and a validation set of 20 images. We tested the trained model with a testing set of 20 images. The experimental results indicated that our model successfully enhanced the image quality with a larger image-size when compared with the common existing image enlargement and the original SRCNN method.


title: Experimental study of transport properties of Weyl semimetal LaAlGe thin films grown by molecular beam epitaxy

abstract: Rare earth compounds display diverse electronic, magnetic, and magneto-transport properties. Recently these compounds of the type RAlGe (R = La, Ce, Pr) have been shown to exhibit Weyl semimetallic behavior. In this work, we have investigated the crystal structure, electronic, and magneto-transport properties of the Weyl semimetal LaAlGe thin films grown by molecular beam epitaxy. The temperature dependence of longitudinal resistivity at different magnetic fields is discussed. Observations of magnetoresistances and Hall effect at different temperatures and their evolution with magnetic field up to 6 T are also discussed with relevant mechanisms. We have observed positive unsaturated magnetoresistances, with a small quadratic contribution at low temperatures, which tends to saturate at higher fields. The Hall measurements confirm the electron-dominated semimetallic conduction with an average charge carrier density of ~ 9.68*10^21 cm^(-3) at room temperature.


title: A two-vertex theorem for normal tilings

abstract: We regard a smooth, $d=2$-dimensional manifold $\mathcal{M}$ and its normal tiling $M$, the cells of which may have non-smooth or smooth vertices (at the latter, two edges meet at 180 degrees.) We denote the average number (per cell) of non-smooth vertices by $\bar v^{\star}$ and we prove that if $M$ is periodic then $v^{\star} \geq 2$ and we show the same result for the monohedral case by an entirely different argument. Our theory also makes a closely related prediction for non-periodic tilings. In 3 dimensions we show a monohedral construction with $\bar v^{\star}=0$.


title: On the dimension of Bergman spaces on $\mathbb{P}^1$

abstract: Inspired by a result by Sz\H{o}ke, we give potential-theoretic characterizations of the dimension of the Bergman space of holomorphic sections of a restriction of a holomorphic line bundle of $\mathbb{P}^1$ to some open set $D\subset\mathbb{P}^1$.


title: Unifying AI Algorithms with Probabilistic Programming using Implicitly Defined Representations

abstract: We introduce Scruff, a new framework for developing AI systems using probabilistic programming. Scruff enables a variety of representations to be included, such as code with stochastic choices, neural networks, differential equations, and constraint systems. These representations are defined implicitly using a set of standardized operations that can be performed on them. General-purpose algorithms are then implemented using these operations, enabling generalization across different representations. Zero, one, or more operation implementations can be provided for any given representation, giving algorithms the flexibility to use the most appropriate available implementations for their purposes and enabling representations to be used in ways that suit their capabilities. In this paper, we explain the general approach of implicitly defined representations and provide a variety of examples of representations at varying degrees of abstraction. We also show how a relatively small set of operations can serve to unify a variety of AI algorithms. Finally, we discuss how algorithms can use policies to choose which operation implementations to use during execution.


title: Higher order first integrals of autonomous dynamical systems

abstract: A theorem is derived which determines higher order first integrals of autonomous holonomic dynamical systems in a general space, provided the collineations and the Killing tensors -- up to the order of the first integral -- of the kinetic metric, defined by the kinetic energy of the system, can be computed. The theorem is applied in the case of Newtonian autonomous conservative dynamical systems of two degrees of freedom, where known and new integrable and superintegrable potentials that admit cubic first integrals are determined.


title: Distcomp: Comparing distributions

abstract: The distcomp command is introduced and illustrated. The command assesses whether or not two distributions differ at each possible value while controlling the probability of any false positive, even in finite samples. Syntax and the underlying methodology (from Goldman and Kaplan, 2018) are discussed. Multiple examples illustrate the distcomp command, including revisiting the experimental data of Gneezy and List (2006) and the regression discontinuity design of Cattaneo, Frandsen, and Titiunik (2015).


title: Time-discretization of a plasma-neutral MHD model with a semi-implicit leapfrog algorithm

abstract: The semi-implicit leapfrog time-discretization is a workhorse algorithm for initial-value MHD codes to bridge between vastly separated time scales. Inclusion of atomic interactions with neutrals breaks the functional structure of the MHD equations that exploited by the leapfrog. We address how to best integrate atomic physics into the semi-implicit leapfrog. Following the Crank-Nicolson method, one approach is to time-center the atomic interactions in the linear solver and use a Newton method to include the nonlinear contributions. Alternatively, another family of methods are based on operator-splitting the terms associated with the atomic interactions using a Strang-splitting technique. These methods naturally break equations into constituent ODE and PDE parts and preserve the structure exploited by the semi-implicit leapfrog. We study the accuracy and efficiency of these methods through a battery of 0D and 1D cases and show that a second-order-in-time Douglas-Rachford inspired coupling between the ODE and PDE advances is effective in reducing the time-discretization error to be comparable to that of Crank-Nicolson with Newton iteration of the nonlinear terms. Splitting ODE and PDE parts results in independent matrix solves for each field which reduces the computational cost considerably and provides parallelization over species relative to Crank-Nicolson.


title: Task-aware Privacy Preservation for Multi-dimensional Data

abstract: Local differential privacy (LDP) can be adopted to anonymize richer user data attributes that will be input to sophisticated machine learning (ML) tasks. However, today's LDP approaches are largely task-agnostic and often lead to severe performance loss -- they simply inject noise to all data attributes according to a given privacy budget, regardless of what features are most relevant for the ultimate task. In this paper, we address how to significantly improve the ultimate task performance with multi-dimensional user data by considering a task-aware privacy preservation problem. The key idea is to use an encoder-decoder framework to learn (and anonymize) a task-relevant latent representation of user data. We obtain an analytical near-optimal solution for the linear setting with mean-squared error (MSE) task loss. We also provide an approximate solution through a gradient-based learning algorithm for general nonlinear cases. Extensive experiments demonstrate that our task-aware approach significantly improves ultimate task accuracy compared to standard benchmark LDP approaches with the same level of privacy guarantee.


title: Shape-aware Multi-Person Pose Estimation from Multi-View Images

abstract: In this paper we contribute a simple yet effective approach for estimating 3D poses of multiple people from multi-view images. Our proposed coarse-to-fine pipeline first aggregates noisy 2D observations from multiple camera views into 3D space and then associates them into individual instances based on a confidence-aware majority voting technique. The final pose estimates are attained from a novel optimization scheme which links high-confidence multi-view 2D observations and 3D joint candidates. Moreover, a statistical parametric body model such as SMPL is leveraged as a regularizing prior for these 3D joint candidates. Specifically, both 3D poses and SMPL parameters are optimized jointly in an alternating fashion. Here the parametric models help in correcting implausible 3D pose estimates and filling in missing joint detections while updated 3D poses in turn guide obtaining better SMPL estimations. By linking 2D and 3D observations, our method is both accurate and generalizes to different data sources because it better decouples the final 3D pose from the inter-person constellation and is more robust to noisy 2D detections. We systematically evaluate our method on public datasets and achieve state-of-the-art performance. The code and video will be available on the project page: https://ait.ethz.ch/projects/2021/multi-human-pose/.


title: A Formal Characterization of Black-Box System Safety Performance with Scenario Sampling

abstract: A typical scenario-based evaluation framework seeks to characterize a black-box system's safety performance (e.g., failure rate) through repeatedly sampling initialization configurations (scenario sampling) and executing a certain test policy for scenario propagation (scenario testing) with the black-box system involved as the test subject. In this letter, we first present a novel safety evaluation criterion that seeks to characterize the actual operational domain within which the test subject would remain safe indefinitely with high probability. By formulating the black-box testing scenario as a dynamic system, we show that the presented problem is equivalent to finding a certain "almost" robustly forward invariant set for the given system. Second, for an arbitrary scenario testing strategy, we propose a scenario sampling algorithm that is provably asymptotically optimal in obtaining the safe invariant set with arbitrarily high accuracy. Moreover, as one considers different testing strategies (e.g., biased sampling of safety-critical cases), we show that the proposed algorithm still converges to the unbiased approximation of the safety characterization outcome if the scenario testing satisfies a certain condition. Finally, the effectiveness of the presented scenario sampling algorithms and various theoretical properties are demonstrated in a case study of the safety evaluation of a control barrier function-based mobile robot collision avoidance system.


title: OTTR: Off-Road Trajectory Tracking using Reinforcement Learning

abstract: In this work, we present a novel Reinforcement Learning (RL) algorithm for the off-road trajectory tracking problem. Off-road environments involve varying terrain types and elevations, and it is difficult to model the interaction dynamics of specific off-road vehicles with such a diverse and complex environment. Standard RL policies trained on a simulator will fail to operate in such challenging real-world settings. Instead of using a naive domain randomization approach, we propose an innovative supervised-learning based approach for overcoming the sim-to-real gap problem. Our approach efficiently exploits the limited real-world data available to adapt the baseline RL policy obtained using a simple kinematics simulator. This avoids the need for modeling the diverse and complex interaction of the vehicle with off-road environments. We evaluate the performance of the proposed algorithm using two different off-road vehicles, Warthog and Moose. Compared to the standard ILQR approach, our proposed approach achieves a 30% and 50% reduction in cross track error in Warthog and Moose, respectively, by utilizing only 30 minutes of real-world driving data.