title: CREATIVESUMM: Shared Task on Automatic Summarization for Creative Writing

abstract: This paper introduces the shared task of summarizing documents in several creative domains, namely literary texts, movie scripts, and television scripts. Summarizing these creative documents requires making complex literary interpretations, as well as understanding non-trivial temporal dependencies in texts containing varied styles of plot development and narrative structure. This poses unique challenges and is yet underexplored for text summarization systems. In this shared task, we introduce four sub-tasks and their corresponding datasets, focusing on summarizing books, movie scripts, primetime television scripts, and daytime soap opera scripts. We detail the process of curating these datasets for the task, as well as the metrics used for the evaluation of the submissions. As part of the CREATIVESUMM workshop at COLING 2022, the shared task attracted 18 submissions in total. We discuss the submissions and the baselines for each sub-task in this paper, along with directions for facilitating future work in the field.


title: Compact pseudoconcave sets

abstract: Replying to three questions posed by N. Shcherbina, we show that a compact psudoconcave set can have the core smaller than itself, that the core of a compact set must be pseudoconcave, and that it can be decomposed into compact pseudoconcave sets on which all smooth plurisubharmonic functions are constant.


title: Finite $3$-connected-set-homogeneous locally $2\K_n$ graphs and $s$-arc-transitive graphs

abstract: In this paper, all graphs are assumed to be finite. For $s\geq 1$ and a graph $\G$, if for every pair of isomorphic connected induced subgraphs on at most $s$ vertices there exists an automorphism of $\G$ mapping the first to the second, then we say that $\G$ is $s$-connected-set-homogeneous, and if every isomorphism between two isomorphic connected induced subgraphs on at most $s$ vertices can be extended to an automorphism of $\G$, then we say that $\G$ is $s$-connected-homogeneous. For $n\geq 1$, a graph $\G$ is said to be locally $2\K_n$ if the subgraph $[\G(u)]$ induced on the set of vertices of $\G$ adjacent to a given vertex $u$ is isomorphic to $2\K_n$. Note that $2$-connected-set-homogeneous but not $2$-connected-homogeneous graphs are just the half-arc-transitive graphs which are a quite active topic in algebraic graph theory. Motivated by this, we posed the problem of characterizing or classifying $3$-connected-set-homogeneous graphs of girth $3$ which are not $3$-connected-homogeneous in (Eur. J. Combin. 93 (2021) 103275). Until now, there have been only two known families of $3$-connected-set-homogeneous graphs of girth $3$ which are not $3$-connected-homogeneous, and these graphs are locally $2\K_n$ with $n=2$ or $4$. In this paper, we complete the classification of finite $3$-connected-set-homogeneous graphs which are locally $2\K_n$ with $n\geq 2$, and all such graphs are line graphs of some specific $2$-arc-transitive graphs. Furthermore, we give a good description of finite $3$-connected-set-homogeneous but not $3$-connected-homogeneous graphs which are locally $2\K_n$ and have solvable automorphism groups. This is then used to construct some new $3$-connected-set-homogeneous but not $3$-connected-homogeneous graphs as well as some new $2$-arc-transitive graphs.


title: Some remarks on simplified double porosity model of immiscible incompressible two-phase flow

abstract: The paper is devoted to the derivation, by linearization, of simplified (fully homogenized) homogenized models of an immiscible incompressible two-phase flow in double porosity media in the case of thin fissures. In a simplified dual porosity model derived previously by the authors the matrix-fracture source term is approximated by a convolution type source term. This approach enables to exclude the cell problem, in form of the imbibition equation, from the global double porosity model. In this paper we propose a new linear version of the imbibition equation which leads to a new simplified dual porosity model. We also present numerical simulations which show that the matrix-fracture exchange term based on this new linearization procedure gives a better approximation of the exact one than the corresponding exchange term obtained earlier by the authors.


title: Sparse Equidistribution of Geometric Invariants of Real Quadratic Fields

abstract: Duke, Imamo\={g}lu, and T\'oth have recently constructed a new geometric invariant, a hyperbolic orbifold, associated to each narrow ideal class of a real quadratic field. Furthermore, they have shown that the projection of these hyperbolic orbifolds onto the modular surface $\Gamma \backslash \mathbb{H}$ equidistributes on average over a genus of the narrow class group as the fundamental discriminant $D$ of the real quadratic field tends to infinity. We extend this construction of hyperbolic orbifolds to allow for a level structure, akin to Heegner points and closed geodesics of level $q$. Additionally, we refine this equidistribution result in several directions. First, we investigate sparse equidistribution in the level aspect, where we prove the equidistribution of level $q$ hyperbolic orbifolds when restricted to a translate of $\Gamma \backslash \mathbb{H}$ in $\Gamma_0(q) \backslash \mathbb{H}$, which presents some new interesting features. Second, we explore sparse equidistribution in the subgroup aspect, namely equidistribution on average over small subgroups of the narrow class group. Third, we prove small scale equidistribution and give upper bounds for the discrepancy. Behind these refinements is a new interpretation of the Weyl sums arising in these equidistribution problems in terms of ad\`elic period integrals, which in turn are related to Rankin-Selberg $L$-functions via Waldspurger's formula. The key remaining inputs are hybrid subconvex bounds for these $L$-functions and a certain homological version of the sup-norm problem.


title: A Closer Cut: Computing Near-Optimal Lawn Mowing Tours

abstract: For a given polygonal region $P$, the Lawn Mowing Problem (LMP) asks for a shortest tour $T$ that gets within Euclidean distance 1 of every point in $P$; this is equivalent to computing a shortest tour for a unit-disk cutter $C$ that covers all of $P$. As a geometric optimization problem of natural practical and theoretical importance, the LMP generalizes and combines several notoriously difficult problems, including minimum covering by disks, the Traveling Salesman Problem with neighborhoods (TSPN), and the Art Gallery Problem (AGP). In this paper, we conduct the first study of the Lawn Mowing Problem with a focus on practical computation of near-optimal solutions. We provide new theoretical insights: Optimal solutions are polygonal paths with a bounded number of vertices, allowing a restriction to straight-line solutions; on the other hand, there can be relatively simple instances for which optimal solutions require a large class of irrational coordinates. On the practical side, we present a primal-dual approach with provable convergence properties based on solving a special case of the TSPN restricted to witness sets. In each iteration, this establishes both a valid solution and a valid lower bound, and thereby a bound on the remaining optimality gap. As we demonstrate in an extensive computational study, this allows us to achieve provably optimal and near-optimal solutions for a large spectrum of benchmark instances with up to 2000 vertices.


title: The electron cyclotron drift instability: a comparison of particle-in-cell and continuum Vlasov simulations

abstract: The linear and nonlinear characteristics of the electron cyclotron drift instability (ECDI) have been studied through the particle-in-cell (PIC) and continuum Vlasov simulation methods in connection with the effects of the azimuthal length (in the $E \times B$ direction) on the simulations. Simulation results for a long azimuthal length (17.82 cm $= 627\;v_d/\omega_{ce}$, where $\omega_{ce}$ is the electron cyclotron frequency and $v_d$ is the $E\times B$ drift of the electrons) are reported, for which a high resolution is achieved in Fourier space. For simulations with a long azimuthal length, the linear growth rates of the PIC simulations show a considerable discrepancy with the theory, whereas the linear growth rate of the Vlasov simulations remains close to the theory. In the nonlinear regime, the inverse cascade is shown in both PIC and Vlasov simulations with a sufficiently large azimuthal length. In simulations with a short azimuthal length, however, the inverse cascade is barely observed. Instead, the PIC simulations with a short azimuthal length (0.5625 cm $=19.8\;v_d/\omega_{ce}$) show an essentially continuous nonlinear dispersion, similar to what is predicted by the ion-sound turbulence theory. It is shown that, in the PIC and Vlasov simulations, the inverse cascade coincides with the formation and merging of electron structures in phase space. This process, however, terminates differently in the PIC simulations compared with the Vlasov simulations. Larger amplitudes of ECDI fluctuations are observed in the PIC simulations compared with the Vlasov simulations, leading to an intensified electron heating and anomalous current. This suggests that the statistical noise of PIC simulations might contribute to the extreme electron heating that has been observed in previous studies.


title: A new conservative discontinuous Galerkin method via implicit penalization for the generalized KdV equation

abstract: We design, analyze, and implement a new conservative Discontinuous Galerkin (DG) method for the simulation of solitary wave solutions to the generalized Korteweg-de Vries (KdV) Equation. The key feature of our method is the conservation, at the numerical level, of the mass, energy and Hamiltonian that are conserved by exact solutions of all KdV equations. To our knowledge, this is the first DG method that conserves all these three quantities, a property critical for the accurate long-time evolution of solitary waves. To achieve the desired conservation properties, our novel idea is to introduce two stabilization parameters in the numerical fluxes as new unknowns which then allow us to enforce the conservation of energy and Hamiltonian in the formulation of the numerical scheme. We prove the conservation properties of the scheme which are corroborated by numerical tests. This idea of achieving conservation properties by implicitly defining penalization parameters, that are traditionally specified a priori, can serve as a framework for designing physics-preserving numerical methods for other types of problems.


title: Spectral bounds for exit times on metric measure Dirichlet spaces and applications

abstract: Assuming the heat kernel on a doubling Dirichlet metric measure space has a generalized sub-Gaussian bound, we prove an asymptotically sharp spectral upper bound on the survival probability of the associated diffusion process. As a consequence, we can show that the supremum of the mean exit time over all starting points is finite if and only if the bottom of the spectrum is positive. Among several applications, we show that the spectral upper bound on the survival probability implies a bound for the Hot Spots constant for Riemannian manifolds. Our results apply to many interesting examples including Carnot groups, sub-Riemannian manifolds with transverse symmetries, and fractals.


title: Understanding ME? Multimodal Evaluation for Fine-grained Visual Commonsense

abstract: Visual commonsense understanding requires Vision Language (VL) models to not only understand image and text but also cross-reference in-between to fully integrate and achieve comprehension of the visual scene described. Recently, various approaches have been developed and have achieved high performance on visual commonsense benchmarks. However, it is unclear whether the models really understand the visual scene and underlying commonsense knowledge due to limited evaluation data resources. To provide an in-depth analysis, we present a Multimodal Evaluation (ME) pipeline to automatically generate question-answer pairs to test models' understanding of the visual scene, text, and related knowledge. We then take a step further to show that training with the ME data boosts the model's performance in standard VCR evaluation. Lastly, our in-depth analysis and comparison reveal interesting findings: (1) semantically low-level information can assist the learning of high-level information but not the opposite; (2) visual information is generally under utilization compared with text.


title: An improved method of delta summation for faster current value selection across filtered subsets of interval and temporal relational data

abstract: Aggregation in relational databases is accomplished through hashing and sorting interval data, which is computationally expensive and scales poorly as the data volumes grow. In this paper, we show how quantitative interval and time-series data in relational attributes can be represented using delta summary values rather than absolute values. The need for sorting to determine the row corresponding to some maximum timestamp is negated, reducing the time complexity from at least O(n log(n)) towards O(n) and improving query execution times. We illustrate this new method in the relational algebra, present the implementation algorithmically, and test an implementation in two leading RDBMS products against the use of normal equivalents. We found this delta summation technique to be most effective for use cases with additive, numerical data upon which it is necessary to frequently obtain the latest values, and where the row cardinalities are in the order of 10^5. Our experiments found the proposed new delta summation technique could execute faster than the equivalent standard selection method by up to 22.4%, reducing the overall query cost in some circumstances by up to 24.0%, reducing I/O load by up to 60.6%, but with increased query costs for write operations, an increase in CPU time and memory allocation, uncertain performance with very low or very high cardinalities and inconsistent results across different RDBMS platforms.


title: The Star-Planet Activity Research CubeSat (SPARCS): Determining Inputs to Planetary Habitability

abstract: Seventy-five billion low-mass stars in our galaxy host at least one small planet in their habitable zone (HZ). The stellar ultraviolet (UV) radiation received by the planets is strong and highly variable, and has consequences for atmospheric loss, composition, and habitability. SPARCS is a NASA-funded mission to characterize the quiescent and flare UV emission from low-mass stars, by observing 10 to 20 low-mass stars, over timescales of days, simultaneously in two UV bands: 153-171 nm and 260-300 nm. SPARCS Sun-synchronous terminator orbit allows for long periods of uninterrupted observations, reaching 10s of days for some targets. The payload consists of a 10 cm-class telescope, a dichroic element, UV detectors and associated electronics, a thermal control system, and an on-board processor. The payload is hosted on a Blue Canyon Technologies 6U CubeSat. SPARCS hosts several technology innovations that have broad applicability to other missions. The payload demonstrates the use of "2D-doped" (i.e., delta- and superlattice-doped) detectors and detector-integrated metal dielectric filters in space. This detector technology provides ~5x larger quantum efficiency than NASA's GALEX detectors. In addition, SPARCS' payload processor provides dynamic exposure control, automatically adjusting the exposure time to avoid flare saturation and to time-resolve the strongest stellar flares. A simple passive cooling system maintains the detector temperature under 238K to minimize dark current. The spacecraft bus provides pointing jitter smaller than 6", minimizing the impact of flat-field errors, dark current, and read-noise. All these elements enable competitive astrophysics science within a CubeSat platform. SPARCS is currently in the final design and fabrication phase (Phase C in the NASA context). It will be launched in 2024, for a primary science mission of one year.


title: Influence of extragalactic magnetic fields on extragalactic cascade gamma-ray emission

abstract: We discuss the influence of extragalactic magnetic fields on the intensity of gamma-ray emission produced in electromagnetic cascades from ultra-high energy cosmic rays propagating in extragalactic space. Both cosmic rays and cascade particles propagate mostly out of galaxies, galactic clusters, and large scale structures as their relative volume is small. Therefore their magnetic fields weakly affect emission produced in cascades. Yet estimates of this influence can be useful searching for dark matter particles, when components of extragalactic gamma-ray background should be known, including cascade gamma-ray emission. To study magnetic field influence on cascade emission we analyze cosmic particle propagation in fields of ~10^(-6) and 10^(-12) G (the former is typical inside galaxies and clusters and the latter is common in voids and outside galaxies and clusters). The calculated spectra of cascade gamma-ray emission are similar in the range of ~10^7-10^9 eV. So analyzing cascade emission in this range it is not necessary to specify models of extragalactic magnetic field.


title: Successive interacting coronal mass ejections: How to create a perfect storm?

abstract: Coronal mass ejections (CMEs) are the largest type of eruptions on the Sun and the main driver of severe space weather at the Earth. In this study, we implement a force-free spheromak CME description within 3-D magnetohydrodynamic simulations to parametrically evaluate successive interacting CMEs within a representative heliosphere. We explore CME-CME interactions for a range of orientations, launch time variations and CME handedness and quantify their geo-effectiveness via the primary solar wind variables and empirical measures of the disturbance storm time index and subsolar magnetopause standoff distance. We show how the interaction of two moderate CMEs between the Sun and the Earth can translate into extreme conditions at the Earth and how CME-CME interactions at different radial distances can maximise different solar wind variables that induce different geophysical impacts. In particular, we demonstrate how the orientation and handedness of a given CME can have a significant impact on the conservation and loss of magnetic flux, and consequently B$_z$, due to magnetic reconnection with the interplanetary magnetic field. This study thus implicates identification of CME chirality in the solar corona as an early diagnostic for forecasting geomagnetic storms involving multiple CMEs.


title: How do tidal waves interact with convective vortices in rapidly-rotating planets and stars?

abstract: The dissipation of tidal inertial waves in planetary and stellar convective regions is one of the key mechanisms that drive the evolution of star-planet/planet-moon systems. In this context, the interaction between tidal inertial waves and turbulent convective flows must be modelled in a realistic and robust way. In the state-of-the-art simulations, the friction applied by convection on tidal waves is modelled most of the time by an effective eddy-viscosity. This approach may be valid when the characteristic length scales of convective eddies are smaller than those of tidal waves. However, it becomes highly questionable in the case where tidal waves interact with potentially stable large-scale vortices, as those observed at the pole of Jupiter and Saturn. They are potentially triggered by convection in rapidly-rotating bodies in which the Coriolis acceleration forms the flow in columnar vortical structures along the direction of the rotation axis. In this paper, we investigate the complex interactions between a tidal inertial wave and a columnar convective vortex. We use a quasi-geostrophic semi-analytical model of a convective columnar vortex. We perform linear stability analysis to identify the unstable regime and conduct linear numerical simulations for the interactions between the convective vortex and an incoming tidal inertial wave. We verify that in the unstable regime, an incoming tidal inertial wave triggers the most unstable mode of the vortex leading to turbulent dissipation. For stable vortices, the wave-vortex interaction leads to the momentum mixing while it creates a low-velocity region around the vortex core and a new wave-like perturbation in the form of a progressive wave radiating in the far field. The emission of this secondary wave is the strongest when the wavelength of the incoming wave is close to the characteristic size of the vortex.


title: Archaeoastronomical study on the north-central coast of Peru

abstract: The Caral civilization developed on the north-central coast of Peru and had an occupation period between 2870 and 1970 years BC. The first studies carried out in the field of archaeoastronomy showed evidence of possible astronomical orientations in some buildings of its capital city, the Ciudad Sagrada de Caral. However, methodological issues cast doubt on these conclusions. A recent study carried out a more general statistical analysis, which covered a total of 55 architectural structures distributed in ten urban settlements that were part of this civilization, thus managing to identify topographic and astronomical orientation patterns. Based on this evidence, we propose to carry out a new study focused on the capital city, with the objective of analyzing the orientation pattern of the city, placing emphasis on the analysis of the most important religious and administrative structures in order to determine their functionality and their possible links with relevant astronomical objects. The study will include field work to measure the various structures and the subsequent statistical analysis of the data, using declination histograms, density functions and probability tests.


title: Blockchain-Assisted Intelligent Symbiotic Radio in Space-Air-Ground Integrated Networks

abstract: In a space-air-ground integrated network (SAGIN), managing resources for the growing number of highly-dynamic and heterogeneous radios is a challenging task. Symbiotic communication (SC) is a novel paradigm, which leverages the analogy of the natural ecosystem in biology to create a radio ecosystem in wireless networks that achieves cooperative service exchange and resource sharing, i.e., service/resource trading, among numerous radios. As a result, the potential of symbiotic communication can be exploited to enhance resource management in SAGIN. Despite the fact that different radio resource bottlenecks can complement each other via symbiotic relationships, unreliable information sharing among heterogeneous radios and multi-dimensional resources managing under diverse service requests impose critical challenges on trusted trading and intelligent decision-making. In this article, we propose a secure and smart symbiotic SAGIN (S^4) framework by using blockchain for ensuring trusted trading among heterogeneous radios and machine learning (ML) for guiding complex service/resource trading. A case study demonstrates that our proposed S^4 framework provides better service with rational resource management when compared with existing schemes. Finally, we discuss several potential research directions for future symbiotic SAGIN.


title: Two-Stage Distributionally Robust Conic Linear Programming over 1-Wasserstein Balls

abstract: This paper studies two-stage distributionally robust conic linear programming under constraint uncertainty over type-1 Wasserstein balls. We present optimality conditions for the dual of the worst-case expectation problem, which characterizes worst-case uncertain parameters for its inner maximization problem. This condition offers an alternative proof, a counter-example, and an extension to previous works. Additionally, the condition highlights the potential advantage of a specific distance metric for out-of-sample performance, as exemplified in a numerical study on a facility location problem with demand uncertainty. A cutting-plane-based algorithm and a variety of algorithmic enhancements are proposed with a finite convergence proof under less stringent assumptions.


title: 4DVarNet-SSH: end-to-end learning of variational interpolation schemes for nadir and wide-swath satellite altimetry

abstract: The reconstruction of sea surface currents from satellite altimeter data is a key challenge in spatial oceanography, especially with the upcoming wide-swath SWOT (Surface Ocean and Water Topography) altimeter mission. Operational systems however generally fail to retrieve mesoscale dynamics for horizontal scales below 100km and time-scale below 10 days. Here, we address this challenge through the 4DVarnet framework, an end-to-end neural scheme backed on a variational data assimilation formulation. We introduce a parametrization of the 4DVarNet scheme dedicated to the space-time interpolation of satellite altimeter data. Within an observing system simulation experiment (NATL60), we demonstrate the relevance of the proposed approach both for nadir and nadir+swot altimeter configurations for two contrasted case-study regions in terms of upper ocean dynamics. We report relative improvement with respect to the operational optimal interpolation between 30% and 60% in terms of reconstruction error. Interestingly, for the nadir+swot altimeter configuration, we reach resolved space-time scales below 70km and 7days. The code is open-source to enable reproductibility and future collaborative developments. Beyond its applicability to large-scale domains, we also address uncertainty quantification issues and generalization properties of the proposed learning setting. We discuss further future research avenues and extensions to other ocean data assimilation and space oceanography challenges.


title: Hunting ewinos and a light scalar of $Z_3$-NMSSM with a bino-like dark matter in top squark decays at the LHC

abstract: We study the prospects of a simultaneous hunt at the Large Hadron Collider (LHC) of relatively light electroweakinos and a singlet-like scalar of the $Z_3$-symmetric Next-to-Minimal Supersymmetric Standard Model ($Z_3$-NMSSM) in the cascade decays of not so heavy ($\lesssim 1.5$ TeV) top squarks that are produced in pairs at the LHC which characteristically involve the singlet-like states. We work in a scenario where the lightest (next-to-lightest) SUSY particle is bino (singlino)-like with a mass below 100 GeV ($\gtrsim 100$ GeV), whereas a pair of immediately heavier neutralinos and the lighter chargino are higgsino-like with masses in the range $\sim 500$ GeV -- 1 TeV. Further, the singlet-like scalar present in the spectrum provides a funnel for a rapid enough mutual annihilation of the LSP thus making the latter meet the experimental upper bound on its relic abundance. The scenario is motivated by its ability to offer a bino-like dark matter with such a mass unlike what is now disfavored in the MSSM while avoiding the stringent lower bounds from the LHC experiments on the masses of these involved particles and still remaining reasonably `natural'. We find that while a usual cut-based analysis (CBA) on LHC data worth 300 $\text{fb}^{-1}$ would be unable to discover such excitations, a multivariate analysis (MVA) can be reasonably sensitive to higgsino-like electroweakinos having masses $\gtrsim 650$ GeV when the lighter top squark has a mass $\gtrsim 1$ TeV. On the other hand, with 3000 $\text{fb}^{-1}$ of data these masses become accessible in a CBA while even an MVA on such a data set is unlikely to find these electroweakinos with masses around 1 TeV when the mass of the lighter top squark hits $\sim 1.5$ TeV.


title: A New Conjecture on Hardness of Low-Degree 2-CSP's with Implications to Hardness of Densest $k$-Subgraph and Other Problems

abstract: We propose a new conjecture on hardness of low-degree $2$-CSP's, and show that new hardness of approximation results for Densest $k$-Subgraph and several other problems, including a graph partitioning problem, and a variation of the Graph Crossing Number problem, follow from this conjecture. The conjecture can be viewed as occupying a middle ground between the $d$-to-$1$ conjecture, and hardness results for $2$-CSP's that can be obtained via standard techniques, such as Parallel Repetition combined with standard $2$-prover protocols for the 3SAT problem. We hope that this work will motivate further exploration of hardness of $2$-CSP's in the regimes arising from the conjecture. We believe that a positive resolution of the conjecture will provide a good starting point for further hardness of approximation proofs. Another contribution of our work is proving that the problems that we consider are roughly equivalent from the approximation perspective. Some of these problems arose in previous work, from which it appeared that they may be related to each other. We formalize this relationship in this work.


title: Microscopics of de Sitter Entropy from Precision Holography

abstract: We calculate quantum corrections to the entropy of four-dimensional de Sitter space induced by higher-derivative terms in the gravitational action and by one-loop effects. Employing the intertwinement in semiclassical gravity of Euclidean de Sitter and anti-de Sitter saddles, we embed effective de Sitter gravity theories in M-theory and express the entropy in terms of the regularized Euclidean anti-de Sitter action on an auxiliary $\mathrm{EAdS}_4 \times S^7/\mathbb{Z}_k$ background. We conjecture that the partition function of the holographically dual 3d ABJM CFT determines the explicit form of the corrections to the de Sitter entropy. This includes a logarithmic term, the coefficient of which, we show, agrees with an independent one-loop calculation around the $-S^4 \times S^7/\mathbb{Z}_k$ Euclidean de Sitter saddle. This provides evidence that the microscopic degrees of freedom behind the entropy of four-dimensional de Sitter space in gravitational theories with a holographic dual description are encapsulated by the path integral of the Euclidean CFT on the three-sphere.


title: Accelerating Irregular Applications via Efficient Synchronization and Data Access Techniques

abstract: Irregular applications comprise an increasingly important workload domain for many fields, including bioinformatics, chemistry, physics, social sciences and machine learning. Therefore, achieving high performance and energy efficiency in the execution of emerging irregular applications is of vital importance. This dissertation studies the root causes of inefficiency of irregular applications in modern computing systems, and fundamentally addresses such inefficiencies, by proposing low-overhead synchronization techniques among parallel threads in cooperation with well-crafted data access policies. We make four major contributions to accelerating irregular applications in different contexts including CPU and Near-Data-Processing (NDP) (or Processing-In-Memory (PIM)) systems. First, we design ColorTM, a novel parallel graph coloring algorithm for CPU systems that trades off using synchronization with lower data access costs. Second, we propose SmartPQ, an adaptive priority queue that achieves high performance under all various contention scenarios in Non-Uniform Memory Access CPU systems. Third, we introduce SynCron, the first practical hardware synchronization mechanism tailored for NDP systems. Fourth, we design SparseP, the first library for high-performance Sparse Matrix Vector Multiplication on real PIM systems. We demonstrate that the execution of irregular applications in CPU and NDP/PIM architectures can be significantly accelerated by co-designing lightweight synchronization approaches along with well-crafted data access policies. This dissertation bridges the gap between processor-centric CPU systems and memory-centric PIM systems in the critically-important area of irregular applications. We hope that this dissertation inspires future work in co-designing software algorithms with cutting-edge computing platforms to significantly accelerate emerging irregular applications.


title: Four-body singlet potential energy surface for reactions of calcium monofluoride

abstract: A full six-dimensional Born-Oppenheimer singlet potential energy surface is constructed for the reaction CaF + CaF $\rightarrow$ CaF$_2$ + Ca using a multireference configuration interaction (MRCI) electronic structure calculation. The {\it ab initio} data thus calculated are interpolated by Gaussian process (GP) regression. The four-body potential energy surface features one $D_{2h}$ global minimum and one $C_s$ local minimum, connected by a barrierless transition state that lends insight to the reaction mechanism. This surface is intended to be of use in understanding ultracold chemistry of CaF molecules.


title: Efficient and Accurate Quantized Image Super-Resolution on Mobile NPUs, Mobile AI & AIM 2022 challenge: Report

abstract: Image super-resolution is a common task on mobile and IoT devices, where one often needs to upscale and enhance low-resolution images and video frames. While numerous solutions have been proposed for this problem in the past, they are usually not compatible with low-power mobile NPUs having many computational and memory constraints. In this Mobile AI challenge, we address this problem and propose the participants to design an efficient quantized image super-resolution solution that can demonstrate a real-time performance on mobile NPUs. The participants were provided with the DIV2K dataset and trained INT8 models to do a high-quality 3X image upscaling. The runtime of all models was evaluated on the Synaptics VS680 Smart Home board with a dedicated edge NPU capable of accelerating quantized neural networks. All proposed solutions are fully compatible with the above NPU, demonstrating an up to 60 FPS rate when reconstructing Full HD resolution images. A detailed description of all models developed in the challenge is provided in this paper.


title: Weak approximation on the norm one torus

abstract: For any abelian group $A$, we prove an asymptotic formula for the number of $A$-extensions $K/\mathbb{Q}$ of bounded discriminant such that the associated norm one torus $R_{K/\mathbb{Q}}^1 \mathbb{G}_m$ satisfies weak approximation. We are also able to produce new results on the Hasse norm principle and to provide new explicit values for the leading constant in some instances of Malle's conjecture.


title: Set-based state estimation for discrete-time constrained nonlinear systems: an approach based on constrained zonotopes and DC programming

abstract: This paper proposes a new state estimator for discrete-time nonlinear dynamical systems with unknown-but-bounded uncertainties and state linear inequality and nonlinear equality constraints. Our algorithm is based on constrained zonotopes (CZs) and on a DC programming approach (DC stands for difference of convex functions). Recently, mean value extension and first-order Taylor extension have been adapted from zonotopes to propagate CZs over nonlinear mappings. Although the resulting algorithms (called CZMV and CZFO) reach better precision than the original zonotopic versions, they carry the sensitivity to the wrapping and dependency effects inherited from interval arithmetic. These interval issues can be mitigated with DC programming since the approximation error bounds are obtained solving optimization problems. A direct benefit of this technique is the elimination of the dependency effect. Our set-membership filter (called CZDC) offers an alternative solution to CZMV and CZFO. In order to demonstrate the effectiveness of the proposed approach, CZDC is experimented over two numerical examples.


title: Twitter Spam and False Accounts Prevalence, Detection and Characterization: A Survey

abstract: The issue of quantifying and characterizing various forms of social media manipulation and abuse has been at the forefront of the computational social science research community for over a decade. In this paper, I provide a (non-comprehensive) survey of research efforts aimed at estimating the prevalence of spam and false accounts on Twitter, as well as characterizing their use, activity, and behavior. I propose a taxonomy of spam and false accounts, enumerating known techniques used to create and detect them. Then, I summarize studies estimating the prevalence of spam and false accounts on Twitter. Finally, I report on research that illustrates how spam and false accounts are used for scams and frauds, stock market manipulation, political disinformation and deception, conspiracy amplification, coordinated influence, public health misinformation campaigns, radical propaganda and recruitment, and more. I will conclude with a set of recommendations aimed at charting the path forward to combat these problems.


title: Integrability and BRST invariance from BF topological theory

abstract: We consider the BRST invariant effective action of the non-abelian BF topological theory in $1+1$ dimensions with gauge group $Sl(2,\mathbb{R})$. By considering different gauge fixing conditions, the zero-curvature field equation give rise to several well known integrable equations. We prove that each integrable equation together with the associated ghost field evolution equation, obtained from the BF theory, is a BRST invariant system with an infinite sequence of BRST invariant conserved quantities. We construct explicitly the systems and the BRST transformation laws for the KdV sequence (including the KdV, mKdV and CKdV equations) and Harry Dym integrable equation.


title: Badly approximable infinite products of quadratic polynomials

abstract: We provide a number of conditions on the rational numbers $u$ and $v$ which ensure that the Laurent series $g_{u,v}(x):=\prod_{t=0}^\infty (1+ux^{-3^t} + vx^{-2\cdot 3^t})$ is badly approximable.


title: Joint localization of directed polymers

abstract: We consider $(1+1)$-dimensional directed polymers in a random potential and provide sufficient conditions guaranteeing joint localization. Joint localization means that for typical realizations of the environment, and for polymers started at different starting points, all the associated endpoint distributions localize in a common random region that does not grow with the length of the polymer. In particular, we prove that joint localization holds when the reference random walk of the polymer model is either a simple symmetric lattice walk or a Gaussian random walk. We also prove that the very strong disorder property holds for a large class of space-continuous polymer models, implying the usual single polymer localization.


title: Unstable dimension variability, heterodimensional cycles, and blenders in the border-collision normal form

abstract: Chaotic attractors commonly contain periodic solutions with unstable manifolds of different dimensions. This allows for a zoo of dynamical phenomena not possible for hyperbolic attractors. The purpose of this Letter is to demonstrate these phenomena in the border-collision normal form. This is a continuous, piecewise-linear family of maps that is physically relevant as it captures the dynamics created in border-collision bifurcations in diverse applications. Since the maps are piecewise-linear they are relatively amenable to an exact analysis and we are able to explicitly identify parameter values for heterodimensional cycles and blenders. For a one-parameter subfamily we identify bifurcations involved in a transition through unstable dimension variability. This is facilitated by being able to compute periodic solutions quickly and accurately, and the piecewise-linear form should provide a useful test-bed for further study.


title: Derivative-based SINDy (DSINDy): Addressing the challenge of discovering governing equations from noisy data

abstract: Recent advances in the field of data-driven dynamics allow for the discovery of ODE systems using state measurements. One approach, known as Sparse Identification of Nonlinear Dynamics (SINDy), assumes the dynamics are sparse within a predetermined basis in the states and finds the expansion coefficients through linear regression with sparsity constraints. This approach requires an accurate estimation of the state time derivatives, which is not necessarily possible in the high-noise regime without additional constraints. We present an approach called Derivative-based SINDy (DSINDy) that combines two novel methods to improve ODE recovery at high noise levels. First, we denoise the state variables by applying a projection operator that leverages the assumed basis for the system dynamics. Second, we use a second order cone program (SOCP) to find the derivative and governing equations simultaneously. We derive theoretical results for the projection-based denoising step, which allow us to estimate the values of hyperparameters used in the SOCP formulation. This underlying theory helps limit the number of required user-specified parameters. We present results demonstrating that our approach leads to improved system recovery for the Van der Pol oscillator, the Duffing oscillator, and the R\"{o}ssler attractor.


title: Correlation Analysis of Decaying Sterile Neutrino Dark Matter in the Context of the SRG Mission

abstract: We provide a correlation analysis of signatures associated with traces of the dark matter decay and the galaxy spatial distribution according to the 2MRS catalog of galaxies. Signature data analysis plays an important role in the context of current and future observations and cosmological constraints. Attention is paid to the constraints that can be obtained for decaying sterile neutrinos when analyzing observations in the context of the Spectr-Roentegn-Gamma (SRG) mission. We study the correlation spectra of dark matter and galaxies, which can be obtained both for the eROSITA telescope and for the first time for the ART-XC telescope. The analysis is carried out both within the framework of the Limber approximation and within the framework of the extended Limber approximation, which makes it possible to more accurately study the power spectra in the region of small multipoles. We calculate the power spectra in both approaches and examine the contribution of different ranges of multipoles to the resulting constraints on sterile neutrino parameters.


title: When Less is More: On the Value of "Co-training" for Semi-Supervised Software Defect Predictors

abstract: Labeling a module defective or non-defective is an expensive task. Hence, there are often limits on how much-labeled data is available for training. Semi-supervised classifiers use far fewer labels for training models. However, there are numerous semi-supervised methods, including self-labeling, co-training, maximal-margin, and graph-based methods, to name a few. Only a handful of these methods have been tested in SE for (e.g.) predicting defects and even there, those methods have been tested on just a handful of projects. This paper applies a wide range of 55 semi-supervised learners to over 714 projects. We find that semi-supervised "co-training methods" work significantly better than other approaches. Specifically, after labeling, just 2.5% of data, then make predictions that are competitive to those using 100% of the data. That said, co-training needs to be used cautiously since the specific choice of co-training methods needs to be carefully selected based on a user's specific goals. Also, we warn that a commonly-used co-training method ("multi-view"-- where different learners get different sets of columns) does not improve predictions (while adding too much to the run time costs 11 hours vs. 1.8 hours). It is an open question, worthy of future work, to test if these reductions can be seen in other areas of software analytics. To assist with exploring other areas, all the codes used are available at https://github.com/ai-se/Semi-Supervised.


title: Static and dynamic magnetic properties of the spin-5/2 triangle lattice antiferromagnet Na3Fe(PO4)2 studied by 31P NMR

abstract: $^{31}$P nuclear magnetic resonance (NMR) measurements have been carried out to investigate the magnetic properties and spin dynamics of Fe$^{3+}$ ($S$ = 5/2) spins in the two-dimensional triangular lattice (TL) compound Na$_3$Fe(PO$_4$)$_2$. The temperature ($T$) dependence of nuclear spin-lattice relaxation rates ($1/T_1$) shows a clear peak around N\'eel temperature, $T_{\rm N} = 10.9$~K, corresponding to an antiferromagnetic (AFM) transition. From the temparature dependence of NMR shift ($K$) above $T_{\rm N}$, an exchange coupling between Fe$^{3+}$ spins was estimated to be $J/k_{\rm B}\simeq 1.9$~K using the spin-5/2 Heisenberg isotropic-TL model. The temperature dependence of $1/T_1T$ divided by the magnetic susceptibility ($\chi$), $1/T_1T\chi$, above $T_{\rm N}$ proves the AFM nature of spin fluctuations below $\sim$ 50 K in the paramagnetic state. In the magnetically ordered state below $T_{\rm N}$, the characteristic rectangular shape of the NMR spectra is observed, indicative of a commensurate AFM state in its ground state. The strong temperature dependence of 1/$T_1$ in the AFM state is well explained by the two-magnon (Raman) process of the spin waves in a 3D antiferromagnet with a spin-anisotropy energy gap of 5.7 K. The temperature dependence of sublattice magnetization is also well reproduced by the spin waves. Those results indicate that the magnetically ordered state of Na$_3$Fe(PO$_4$)$_2$ is a conventional 3D AFM state, and no obvious spin frustration effects were detected in its ground state.


title: Internal feedback in the cortical perception-action loop enables fast and accurate behavior

abstract: Animals move smoothly and reliably in unpredictable environments. Models of sensorimotor control have assumed that sensory information from the environment leads to actions, which then act back on the environment, creating a single, unidirectional perception-action loop. This loop contains internal delays in sensory and motor pathways, which can lead to unstable control. We show here that these delays can be compensated by internal feedback signals that flow backwards, from motor towards sensory areas. Internal feedback is ubiquitous in neural sensorimotor systems and recent advances in control theory show how internal feedback compensates internal delays. This is accomplished by filtering out self-generated and other predictable changes in early sensory areas so that unpredicted, actionable information can be rapidly transmitted toward action by the fastest components. For example, fast, giant neurons are necessarily less accurate than smaller neurons, but they are crucial for fast and accurate behavior. We use a mathematically tractable control model to show that internal feedback has an indispensable role in achieving state estimation, localization of function -- how different parts of cortex control different parts of the body -- and attention, all of which are crucial for effective sensorimotor control. This control model can explain anatomical, physiological and behavioral observations, including motor signals in visual cortex, heterogeneous kinetics of sensory receptors and the presence of giant Betz cells in motor cortex, Meynert cells in visual cortex and giant von Economo cells in the prefrontal cortex of humans as well as internal feedback patterns and unexplained heterogeneity in other neural systems.


title: Calogero model revisited, commuting Hamiltonians, Hurwitz numbers

abstract: The generalized Mironov-Morozov-Natanson (MMN) equation includes a set of commuting operators, which can be considered as Hamiltonians for the quantum Calogero-Sutherland problem with a special value of the coupling constant (free fermion point). These Hamiltonians can be considered as the center of the enveloping algebra of the group $GL_N(C)$. Another commuting series of Hamiltonians is presented, parametrized by an arbitrary matrix $A\in GL_N$. These Hamiltonians are related to the Hurwitz numbers in the same way as in the case of the MMN equation and generate a generalized variant of the Calogero-Surtheland model.


title: Slightly Altruistic Nash Equilibrium for Multi-agent Pursuit-Evasion Games With Input Constraints

abstract: This is an initial manuscript that presents the basic idea of "slightly altruistic Nash equilibrium", "bi-layer game topology", "rolling horizon target selection". This manuscript is just used for peer discussion and joint Ph.D. application affairs rather than submission to any journal. Thus some references are not all provided. The complete paper for submission will be provided in the next version.


title: Resonances for rational Anosov maps on the torus

abstract: A complete description of resonances for rational toral Anosov diffeomorphisms preserving certain Reinhardt domains is presented. As a consequence it is shown that every homotopy class of two-dimensional Anosov diffeomorphisms contains maps with the sequence of resonances decaying stretched-exponentially. This is achieved by introducing a certain group of rational toral diffeomorphisms and computing the resonances of the respective composition operator considered on suitable anisotropic spaces of hyperfunctions. The class of examples is sufficiently rich to also include non-linear Anosov maps with trivial resonances, or resonances decaying exponentially, as well as with or without area-preservation or reversing symmetries.


title: Non-singular and singular flat bands in tunable acoustic metamaterials

abstract: Dispersionless flat bands can be classified into two types: (1) non-singular flat bands whose eigenmodes are completely characterized by compact localized states; and (2) singular flat bands that have a discontinuity in their Bloch eigenfunctions at a band touching point with an adjacent dispersive band, thereby requiring additional extended states to span their eigenmode space. In this study, we design and numerically demonstrate two-dimensional thin-plate acoustic metamaterials in which tunable flat bands of both kinds can be achieved. Non-singular flat bands are achieved by fine-tuning the ratio of the global tension and the bending stiffness in triangular and honeycomb lattices of plate resonators. A singular flat band arises in a kagome lattice due to the underlying lattice geometry, which can be made degenerate with two additional flat bands by tuning the plate tension. A discrete model of the continuum thin-plate system reveals the interplay of geometric and mechanical factors in determining the existence of flat bands of both types. The singular nature of the kagome lattice flat band is established via a metric called the Hilbert-Schmidt distance calculated between a pair of eigenstates infinitesimally close to the quadratic band touching point. We also simulate an acoustic manifestation of a robust boundary mode arising from the singular flat band and protected by real-space topology in a finite system. Our theoretical and computational study establishes a framework for exploring flat-band physics in a tunable classical system, and for designing acoustic metamaterials with potentially useful sound manipulation capabilities.


title: Optimal Condition Training for Target Source Separation

abstract: Recent research has shown remarkable performance in leveraging multiple extraneous conditional and non-mutually exclusive semantic concepts for sound source separation, allowing the flexibility to extract a given target source based on multiple different queries. In this work, we propose a new optimal condition training (OCT) method for single-channel target source separation, based on greedy parameter updates using the highest performing condition among equivalent conditions associated with a given target source. Our experiments show that the complementary information carried by the diverse semantic concepts significantly helps to disentangle and isolate sources of interest much more efficiently compared to single-conditioned models. Moreover, we propose a variation of OCT with condition refinement, in which an initial conditional vector is adapted to the given mixture and transformed to a more amenable representation for target source extraction. We showcase the effectiveness of OCT on diverse source separation experiments where it improves upon permutation invariant models with oracle assignment and obtains state-of-the-art performance in the more challenging task of text-based source separation, outperforming even dedicated text-only conditioned models.


title: Assessing the Lognormal Distribution Assumption For the Crude Odds Ratio: Implications For Point and Interval Estimation

abstract: The assumption that the sampling distribution of the crude odds ratio (ORcrude) is a log-normal distribution with parameters mu and sigma leads to the incorrect conclusion that the expectation of the log of ORcrude is equal to the parameter mu. In fact, exp(mu) is the median of the lognormal distribution, not the mean. If a different parameter is obtained as the expected value of the lognormal distribution, then this quantity can be used to obtain a new estimate of the true odds ratio (ORtrue). Here, simulations are conducted based on a simple randomized clinical trial study design. The simulations demonstrate that the new estimate of ORtrue (based on the expectation of the lognormal distribution function) yields interval estimates that are more statistically valid than the standard method. These interval estimates are obtained by both a parametric bootstrap method and a calculated percentile method. The statistical conclusion validity of the estimated confidence intervals are based on the intended coverage probability (ie the probability the confidence interval contains ORtrue). Additionally, an interval based hypothesis test based on the improved confidence interval estimate has higher power to reject the null hypothesis that ORtrue is equal to one when the alternative hypothesis is true (ie ORtrue is not equal to one) than the standard hypothesis test when the intervention is protective.


title: Structured Singular Value of a Repeated Complex Full-Block Uncertainty

abstract: The structured singular value (SSV), or mu, is used to assess the robust stability and performance of an uncertain linear time-invariant system. Existing algorithms compute upper and lower bounds on the SSV for structured uncertainties that contain repeated (real or complex) scalars and/or non-repeated complex full blocks. This paper presents algorithms to compute bounds on the SSV for the case of repeated complex full blocks. This specific class of uncertainty is relevant for the input output analysis of many convective systems, such as fluid flows. Specifically, we present a power iteration to compute a lower bound on SSV for the case of repeated complex full blocks. This generalizes existing power iterations for repeated complex scalar and non-repeated complex full blocks. The upper bound can be formulated as a semi-definite program (SDP), which we solve using a standard interior-point method to compute optimal scaling matrices associated with the repeated full blocks. Our implementation of the method only requires gradient information, which improves the computational efficiency of the method. Finally, we test our proposed algorithms on an example model of incompressible fluid flow. The proposed methods provide less conservative bounds as compared to prior results, which ignore the repeated full block structure.


title: Decomposition of class II graphs into two class I graphs

abstract: Mkrtchyan and Steffen [J. Graph Theory, 70 (4), 473--482, 2012] showed that every class II simple graph can be decomposed into a maximum $\Delta$-edge-colorable subgraph and a matching. They further conjectured that every graph $G$ with chromatic index $\Delta(G)+k$ ($k\geq 1$) can be decomposed into a maximum $\Delta(G)$-edge-colorable subgraph (not necessarily class I) and a $k$-edge-colorable subgraph. In this paper, we first generalize their result to multigraphs and show that every multigraph $G$ with multiplicity $\mu$ can be decomposed into a maximum $\Delta(G)$-edge-colorable subgraph and a subgraph with maximum degree at most $\mu$. Then we prove that every graph $G$ with chromatic index $\Delta(G)+k$ can be decomposed into two class I subgraphs $H_1$ and $H_2$ such that $\Delta(H_1) = \Delta(G)$ and $\Delta(H_2) = k$, which is a variation of their conjecture.


title: EEG-based performance-driven adaptive automated hazard alerting system in security surveillance support

abstract: Computer-vision technologies have emerged to assist security surveillance. However, automation alert/alarm systems often apply a low-beta threshold to avoid misses and generates excessive false alarms. This study proposed an adaptive hazard diagnosis and alarm system with adjustable alert threshold levels based on environmental scenarios and operator's hazard recognition performance. We recorded electroencephalogram (EEG) data during hazard recognition tasks. The linear ballistic accumulator model was used to decompose the response time into several psychological subcomponents, which were further estimated by a Markov chain Monte Carlo algorithm and compared among different types of hazardous scenarios. Participants were most cautious about falling hazards, followed by electricity hazards, and had the least conservative attitude toward structural hazards. Participants were classified into three performance-level subgroups using a latent profile analysis based on task accuracy. We applied the transfer learning paradigm to classify subgroups based on their time-frequency representations of EEG data. Additionally, two continual learning strategies were investigated to ensure a robust adaptation of the model to predict participants' performance levels in different hazardous scenarios. These findings can be leveraged in real-world brain-computer interface applications, which will provide human trust in automation and promote the successful implementation of alarm technologies.


title: Deep Appearance Prefiltering

abstract: Physically based rendering of complex scenes can be prohibitively costly with a potentially unbounded and uneven distribution of complexity across the rendered image. The goal of an ideal level of detail (LoD) method is to make rendering costs independent of the 3D scene complexity, while preserving the appearance of the scene. However, current prefiltering LoD methods are limited in the appearances they can support due to their reliance of approximate models and other heuristics. We propose the first comprehensive multi-scale LoD framework for prefiltering 3D environments with complex geometry and materials (e.g., the Disney BRDF), while maintaining the appearance with respect to the ray-traced reference. Using a multi-scale hierarchy of the scene, we perform a data-driven prefiltering step to obtain an appearance phase function and directional coverage mask at each scale. At the heart of our approach is a novel neural representation that encodes this information into a compact latent form that is easy to decode inside a physically based renderer. Once a scene is baked out, our method requires no original geometry, materials, or textures at render time. We demonstrate that our approach compares favorably to state-of-the-art prefiltering methods and achieves considerable savings in memory for complex scenes.


title: Teaching Blockchain in K9-12: Instruction materials and their assessment

abstract: This paper analyses the feasibility of including emerging IT-topics into secondary education (K9-12). We developed a class on Blockchain for K9-12 students which delivers relevant content but is also indented to motivate them to further engage with the topic. Our course consists of theory material as well as a hands-on application. Gamification techniques are utilised to encourage students and improve the learning experience. The assessment, consisting of 166 students, shows a significant knowledge gain of students exposed to our materials especially when the theory and the hands-on materials are combined. Both, i.e., teaching material and chat application, are freely available.


title: Decarbonizing Indian Electricity Grid

abstract: India, being one of the fastest growing economies of the world, must take a sustainable path for development. India is responsible for 7 percent of global CO2 emissions. The electricity sector accounts for nearly 35 percent of emissions from the country. The switch from fossil fuels to renewable sources is the key in decarbonizing this sector and is considered as the crucial step for climate mitigation. This research investigates the potential of renewable energy sources; wind, solar and hydro. The optimization model developed in this study analyzes various scenarios for the transition to a sustainable future. The results show that India aims to achieve 450 GW of installed capacity from RES is far from a Net Zero future. Results confirm that India has the potential to meet 100 percent of electricity demand in 2030 from RES including wind, solar and hydro. Introducing Social Cost of Carbon is a viable option to reduce emissions in India. However, due to the low cost of coal, high coal taxes do not lead to reduced emissions.


title: Inferring probabilistic Boolean networks from steady-state gene data samples

abstract: Probabilistic Boolean Networks have been proposed for estimating the behaviour of dynamical systems as they combine rule-based modelling with uncertainty principles. Inferring PBNs directly from gene data is challenging however, especially when data is costly to collect and/or noisy, e.g., in the case of gene expression profile data. In this paper, we present a reproducible method for inferring PBNs directly from real gene expression data measurements taken when the system was at a steady state. The steady-state dynamics of PBNs is of special interest in the analysis of biological machinery. The proposed approach does not rely on reconstructing the state evolution of the network, which is computationally intractable for larger networks. We demonstrate the method on samples of real gene expression profiling data from a well-known study on metastatic melanoma. The pipeline is implemented using Python and we make it publicly available.


title: Frequency and Amplitude Optimizations for Magnetic Particle Spectroscopy Applications

abstract: Nowadays, there is a growing interest in the field of magnetic particle spectroscopy (MPS)-based bioassays. MPS monitors the dynamic magnetic response of surface-functionalized magnetic nanoparticles (MNPs) upon excitation by an alternating magnetic field (AMF) to detect various target analytes. This technology has flourished in the past decade due to its low cost, low background magnetic noise interference from biomatrix, and fast response time. A large number of MPS variants have been reported by different groups around the world, with applications ranging from disease diagnosis to foodborne pathogen detection, and virus detection. However, there is an urgent need for guidance on how to optimize the sensitivity of MPS detection by choosing different types of MNPs, AMF modalities, and MPS assay strategies (i.e., volume- and surface-based assays). In this work, we systematically study the effect of AMF frequencies and amplitudes on the responses of single- and multi-core MNPs under two extreme conditions, namely, the bound and unbound states. Our results show that some modalities such as dual-frequency MPS utilizing multicore MNPs are more suitable for surface-based bioassay applications, whereas, single-frequency MPS systems using single- or multi-core MNPs are better suited for volumetric bioassay applications. Furthermore, the bioassay sensitivities for these modalities can be further improved by careful selection of AMF frequencies and amplitudes.


title: Optimal Designs of Two-Phase Case-Control Studies for General Predictor Effects

abstract: Under two-phase designs, the outcome and several covariates and confounders are measured in the first phase, and a new predictor of interest, which may be costly to collect, can be measured on a subsample in the second phase, without incurring the costs of recruiting subjects. By using the information gathered in the first phase, the second-phase subsample can be selected to enhance the efficiency of testing and estimating the effect of the new predictor on the outcome. Past studies have focused on optimal two-phase sampling schemes for statistical inference on local ($\beta = o(1)$) effects of the predictor of interest. In this study, we propose an extension of the two-phase designs that employs an optimal sampling scheme for estimating predictor effects with pseudo conditional likelihood estimators in case-control studies. This approach is applicable to both local and non-local effects. We demonstrate the effectiveness of the proposed sampling scheme through simulation studies and analysis of data from 170 patients hospitalized for treatment of COVID-19. The results show a significant improvement in the estimation of the parameter of interest.


title: The Merger Fraction of Ultramassive White Dwarfs

abstract: We search for merger products among the 25 most massive white dwarfs in the Montreal White Dwarf Database 100 pc sample through follow-up spectroscopy and high-cadence photometry. We find an unusually high fraction, 40%, of magnetic white dwarfs among this population. In addition, we identify four outliers in transverse velocity and detect rapid rotation in five objects. Our results show that $56^{+9}_{-10}$\% of the $M\approx1.3~M_{\odot}$ ultramassive white dwarfs form through mergers. This fraction is significantly higher than expected from the default binary population synthesis calculations using the $\alpha$-prescription (with $\alpha \lambda = 2$), and provides further support for efficient orbital shrinkage, such as with low values of the common envelope efficiency.


title: pyRDDLGym: From RDDL to Gym Environments

abstract: We present pyRDDLGym, a Python framework for auto-generation of OpenAI Gym environments from RDDL declerative description. The discrete time step evolution of variables in RDDL is described by conditional probability functions, which fits naturally into the Gym step scheme. Furthermore, since RDDL is a lifted description, the modification and scaling up of environments to support multiple entities and different configurations becomes trivial rather than a tedious process prone to errors. We hope that pyRDDLGym will serve as a new wind in the reinforcement learning community by enabling easy and rapid development of benchmarks due to the unique expressive power of RDDL. By providing explicit access to the model in the RDDL description, pyRDDLGym can also facilitate research on hybrid approaches for learning from interaction while leveraging model knowledge. We present the design and built-in examples of pyRDDLGym, and the additions made to the RDDL language that were incorporated into the framework.


title: M Subdwarf Research III. Spectroscopic Diagnostics for Breaking Parameter Degeneracy

abstract: To understand the parameter degeneracy of M subdwarf spectra at low resolution, we assemble a large number of spectral features in the wavelength range of 0.6-2.5 {\mu}m with band strength quantified by narrowband indices. Based on the index trends of BT-Settl model sequences, we illustrate how the main atmospheric parameters (Teff, log g, [M/H], and [alpha/Fe]) affect each spectral feature differently. Furthermore, we propose a four-step process to determine the four parameters sequentially, which extends the basic idea proposed by Jao et al. Each step contains several spectral features that break the degeneracy effect when determining a specific stellar parameter. Finally, the feasibility of each spectroscopic diagnostic with different spectral qualities is investigated. The result is resolution-independent down to R~200.


title: Negative GMR Effect in current perpendicular-to-plane (Zn,Cr)Te/Cu/Co spin salves

abstract: Magnetic and transport properties are explored in the current perpendicular-to-plane (CPP) spin salves with Cr doped wide band gap semiconductor ZnTe as one of the ferromagnetic electrodes. A negative magnetoresistance is observed in these CPP spin valves at low temperature, with a strong temperature dependence. This effect can be explained by the large difference of spin scattering asymmetry coefficients in (Zn,Cr)Te and Cobalt, due to the very different spin polarizations of the two materials as revealed by the DFT calculation.


title: Knowledge Distillation from Cross Teaching Teachers for Efficient Semi-Supervised Abdominal Organ Segmentation in CT

abstract: For more clinical applications of deep learning models for medical image segmentation, high demands on labeled data and computational resources must be addressed. This study proposes a coarse-to-fine framework with two teacher models and a student model that combines knowledge distillation and cross teaching, a consistency regularization based on pseudo-labels, for efficient semi-supervised learning. The proposed method is demonstrated on the abdominal multi-organ segmentation task in CT images under the MICCAI FLARE 2022 challenge, with mean Dice scores of 0.8429 and 0.8520 in the validation and test sets, respectively.


title: Deep equilibrium models as estimators for continuous latent variables

abstract: Principal Component Analysis (PCA) and its exponential family extensions have three components: observations, latents and parameters of a linear transformation. We consider a generalised setting where the canonical parameters of the exponential family are a nonlinear transformation of the latents. We show explicit relationships between particular neural network architectures and the corresponding statistical models. We find that deep equilibrium models -- a recently introduced class of implicit neural networks -- solve maximum a-posteriori (MAP) estimates for the latents and parameters of the transformation. Our analysis provides a systematic way to relate activation functions, dropout, and layer structure, to statistical assumptions about the observations, thus providing foundational principles for unsupervised DEQs. For hierarchical latents, individual neurons can be interpreted as nodes in a deep graphical model. Our DEQ feature maps are end-to-end differentiable, enabling fine-tuning for downstream tasks.


title: A Gait Triaging Toolkit for Overlapping Acoustic Events in Indoor Environments

abstract: Gait has been used in clinical and healthcare applications to assess the physical and cognitive health of older adults. Acoustic based gait detection is a promising approach to collect gait data of older adults passively and non-intrusively. However, there has been limited work in developing acoustic based gait detectors that can operate in noisy polyphonic acoustic scenes of homes and care homes. We attribute this to the lack of good quality gait datasets from the real-world to train a gait detector on. In this paper, we put forward a novel machine learning based filter which can triage gait audio samples suitable for training machine learning models for gait detection. The filter achieves this by eliminating noisy samples at an f(1) score of 0.85 and prioritising gait samples with distinct spectral features and minimal noise. To demonstrate the effectiveness of the filter, we train and evaluate a deep learning model on gait datasets collected from older adults with and without applying the filter. The model registers an increase of 25 points in its f(1) score on unseen real-word gait data when trained with the filtered gait samples. The proposed filter will help automate the task of manual annotation of gait samples for training acoustic based gait detection models for older adults in indoor environments.


title: Searching for candidates of coalescing binary black holes formed through chemically homogeneous evolution in GWTC-3

abstract: The LIGO, Virgo, and KAGRA (LVK) collaboration has announced 90 coalescing binary black holes (BBHs) with $p_{\rm astro} > 50\%$ to date, however, the origin of their formation channels is still an open scientific question. Given various properties of BBHs (BH component masses and individual spins) inferred using the default priors by the LVK, independent groups have been trying to explain the formation of the BBHs with different formation channels. Of all formation scenarios, the chemically homogeneous evolution (CHE) channel has stood out with distinguishing features, namely, nearly-equal component masses and preferentially high individual spins aligned with the orbital angular momentum. We perform Bayesian inference on the BBH events officially reported in GWTC-3 with astrophysically-predicted priors representing different formation channels of the isolated binary evolution (CEE: common-envelope evolution channel; CHE; SMT: stable mass transfer). Given assumed models, we report strong evidence for GW190517\_055101 being most likely to have formed through the CHE channel. Assuming the BBH events in the subsample are all formed through one of the isolated binary evolution channels, we obtain the lower limits on the local merger rate density of these channels at $11.45 ~\mathrm{Gpc^{-3}~yr^{-1}}$ (CEE), $0.18 ~\mathrm{Gpc^{-3}~yr^{-1}}$ (CHE), and $0.63 ~\mathrm{Gpc^{-3}~yr^{-1}}$ (SMT) at $90\%$ credible level.


title: Deep Reinforcement Learning Microgrid Optimization Strategy Considering Priority Flexible Demand Side

abstract: As an efficient way to integrate multiple distributed energy resources and the user side, a microgrid is mainly faced with the problems of small-scale volatility, uncertainty, intermittency and demand-side uncertainty of DERs. The traditional microgrid has a single form and cannot meet the flexible energy dispatch between the complex demand side and the microgrid. In response to this problem, the overall environment of wind power, thermostatically controlled loads, energy storage systems, price-responsive loads and the main grid is proposed. Secondly, the centralized control of the microgrid operation is convenient for the control of the reactive power and voltage of the distributed power supply and the adjustment of the grid frequency. However, there is a problem in that the flexible loads aggregate and generate peaks during the electricity price valley. The existing research takes into account the power constraints of the microgrid and fails to ensure a sufficient supply of electric energy for a single flexible load. This paper considers the response priority of each unit component of TCLs and ESSs on the basis of the overall environment operation of the microgrid so as to ensure the power supply of the flexible load of the microgrid and save the power input cost to the greatest extent. Finally, the simulation optimization of the environment can be expressed as a Markov decision process process. It combines two stages of offline and online operations in the training process. The addition of multiple threads with the lack of historical data learning leads to low learning efficiency. The asynchronous advantage actor-critic with the experience replay pool memory library is added to solve the data correlation and nonstatic distribution problems during training.


title: Tractable Evaluation of Stein's Unbiased Risk Estimate with Convex Regularizers

abstract: Stein's unbiased risk estimate (SURE) gives an unbiased estimate of the $\ell_2$ risk of any estimator of the mean of a Gaussian random vector. We focus here on the case when the estimator minimizes a quadratic loss term plus a convex regularizer. For these estimators SURE can be evaluated analytically for a few special cases, and generically using recently developed general purpose methods for differentiating through convex optimization problems; these generic methods however do not scale to large problems. In this paper we describe methods for evaluating SURE that handle a wide class of estimators, and also scale to large problem sizes.


title: Terminal orders on arithmetic surfaces

abstract: The local structure of terminal Brauer classes on arithmetic surfaces were classified in [CI21] generalising the classification on geometric surfaces carried out in [CI05]. Part of the interest in these classifications is that it enables the minimal model program to be applied to the noncommutative setting of orders on surfaces. In this paper, we give etale local structure theorems for terminal orders on arithemtic surfaces, at least when the degree is a prime p >5. This generalises the structure theorem given in the geometric case. They can all be explicitly constructed as algebras of matrices over symbols. From this description one sees that such terminal orders all have global dimension two, thus generalising the fact that terminal (commutative) surfaces are smooth and hence homologically regular.


title: MetaBayesDTA: Codeless Bayesian meta-analysis of test accuracy, with or without a gold standard

abstract: Introduction: Despite their applicability, statistical models used for the meta-analysis of test accuracy require specialised knowledge to implement, with the necessary level of expertise having recently increased. This is due to the development and recommendation to use more sophisticated methods; such as those in Version 2 of the Cochrane Handbook for Systematic Reviews of Diagnostic Test Accuracy. This paper describes a web-based application that extends the functionality of previous applications, making many advanced analysis methods more accessible. Methods: We sought to create an extended, stand-alone, Bayesian version of MetaDTA, which (i) has the benefits of previously proposed applications and addresses key limitations of them, (ii) is accessible to researchers who do not have the specific expertise required to fit such models, and (iii) is suitable for experienced analysts. We created the application using Shiny and Stan. Results: We created MetaBayesDTA (https://crsu.shinyapps.io/MetaBayesDTA/), which allows users to conduct meta-analysis of test accuracy, with or without a gold standard. The application addresses several key limitations of other applications. For instance, for the bivariate model, one can conduct subgroup analysis, univariate meta-regression, and comparative test accuracy evaluation. Meanwhile, for the model which does not assume a perfect gold standard, the application can account for the fact that studies use different reference tests. Conclusions: Due to its user-friendliness and broad array of features, MetaBayesDTA should appeal to a wide variety of researchers. We anticipate that the application will encourage wider use of more advanced methods, which ultimately should improve the quality of test accuracy reviews.


title: CR-LSO: Convex Neural Architecture Optimization in the Latent Space of Graph Variational Autoencoder with Input Convex Neural Networks

abstract: In neural architecture search (NAS) methods based on latent space optimization (LSO), a deep generative model is trained to embed discrete neural architectures into a continuous latent space. In this case, different optimization algorithms that operate in the continuous space can be implemented to search neural architectures. However, the optimization of latent variables is challenging for gradient-based LSO since the mapping from the latent space to the architecture performance is generally non-convex. To tackle this problem, this paper develops a convexity regularized latent space optimization (CR-LSO) method, which aims to regularize the learning process of latent space in order to obtain a convex architecture performance mapping. Specifically, CR-LSO trains a graph variational autoencoder (G-VAE) to learn the continuous representations of discrete architectures. Simultaneously, the learning process of latent space is regularized by the guaranteed convexity of input convex neural networks (ICNNs). In this way, the G-VAE is forced to learn a convex mapping from the architecture representation to the architecture performance. Hereafter, the CR-LSO approximates the performance mapping using the ICNN and leverages the estimated gradient to optimize neural architecture representations. Experimental results on three popular NAS benchmarks show that CR-LSO achieves competitive evaluation results in terms of both computational complexity and architecture performance.


title: First $\ell^2$-Betti numbers and proper proximality

abstract: We show that for a countable exact group, having positive first $\ell^2$-Betti number implies proper proximality in this sense of \cite{BoIoPe21}. This is achieved by showing a cocycle superrigidty result for Bernoulli shifts of non-properly proximal groups. We also obtain that Bernoulli shifts of countable, nonamenable, i.c.c., exact, non-properly proximal groups are OE-superrigid.


title: Efficient Domain Coverage for Vehicles with Second-Order Dynamics via Multi-Agent Reinforcement Learning

abstract: Collaborative autonomous multi-agent systems covering a specified area have many potential applications, such as UAV search and rescue, forest fire fighting, and real-time high-resolution monitoring. Traditional approaches for such coverage problems involve designing a model-based control policy based on sensor data. However, designing model-based controllers is challenging, and the state-of-the-art classical control policy still exhibits a large degree of sub-optimality. In this paper, we present a reinforcement learning (RL) approach for the multi-agent efficient domain coverage problem involving agents with second-order dynamics. Our approach is based on the Multi-Agent Proximal Policy Optimization Algorithm (MAPPO). Our proposed network architecture includes the incorporation of LSTM and self-attention, which allows the trained policy to adapt to a variable number of agents. Our trained policy significantly outperforms the state-of-the-art classical control policy. We demonstrate our proposed method in a variety of simulated experiments.


title: Breadth-First Pipeline Parallelism

abstract: We introduce Breadth-First Pipeline Parallelism, a novel training schedule which optimizes the combination of pipeline and data parallelism. Breadth-First Pipeline Parallelism lowers training time, cost and memory usage by combining a high GPU utilization with a small batch size per GPU, and by making use of fully sharded data parallelism. Experimentally, we observed an increase of up to 43% in training throughput for a 52 billion-parameter model using a small batch size per GPU compared to Megatron-LM, which would reduce the training time and cost by the same amount on a large GPU cluster.


title: Signal-to-noise ratio aware minimaxity and higher-order asymptotics

abstract: Since its development, the minimax framework has been one of the corner stones of theoretical statistics, and has contributed to the popularity of many well-known estimators, such as the regularized M-estimators for high-dimensional problems. In this paper, we will first show through the example of sparse Gaussian sequence model, that the theoretical results under the classical minimax framework are insufficient for explaining empirical observations. In particular, both hard and soft thresholding estimators are (asymptotically) minimax, however, in practice they often exhibit sub-optimal performances at various signal-to-noise ratio (SNR) levels. The first contribution of this paper is to demonstrate that this issue can be resolved if the signal-to-noise ratio is taken into account in the construction of the parameter space. We call the resulting minimax framework the signal-to-noise ratio aware minimaxity. The second contribution of this paper is to showcase how one can use higher-order asymptotics to obtain accurate approximations of the SNR-aware minimax risk and discover minimax estimators. The theoretical findings obtained from this refined minimax framework provide new insights and practical guidance for the estimation of sparse signals.


title: MEE: A Novel Multilingual Event Extraction Dataset

abstract: Event Extraction (EE) is one of the fundamental tasks in Information Extraction (IE) that aims to recognize event mentions and their arguments (i.e., participants) from text. Due to its importance, extensive methods and resources have been developed for Event Extraction. However, one limitation of current research for EE involves the under-exploration for non-English languages in which the lack of high-quality multilingual EE datasets for model training and evaluation has been the main hindrance. To address this limitation, we propose a novel Multilingual Event Extraction dataset (MEE) that provides annotation for more than 50K event mentions in 8 typologically different languages. MEE comprehensively annotates data for entity mentions, event triggers and event arguments. We conduct extensive experiments on the proposed dataset to reveal challenges and opportunities for multilingual EE.


title: Constraining neutrino mass in dynamical dark energy cosmologies with the logarithm parametrization and the oscillating parametrization

abstract: We constrain two dynamical dark energy models that are parametrized by the logarithm form of $w(z)=w_{0}+w_{1}\left(\frac{\ln (2+z)}{1+z}-\ln 2\right)$ and the oscillating form of $w(z)=w_{0}+w_{1}\left(\frac{\sin(1+z)}{1+z}-\sin(1)\right)$. Comparing with the Chevallier-Polarski-Linder (CPL) model, the two parametrizations for dark energy can explore the whole evolution history of the universe properly. Using the current mainstream observational data including the cosmic microwave background data and the baryon acoustic oscillation data as well as the type Ia supernovae data, we perform the $\chi^2$ statistic analysis to global fit these models, finding that the logarithm parametrization and the oscillating parametrization are almost as well as the CPL scenario in fitting these data. We make a comparison for the impacts of the dynamical dark energy on the cosmological constraints on the total mass of active neutrinos. We find that the dark energy properties could significantly change the fitting results of neutrino mass. Looser constraints on $\sum m_{\nu}$ are obtained in the logarithm and oscillating models than those derived in the CPL model. Consideration of the possible mass ordering of neutrinos reveals that the most stringent constraint on $\sum m_{\nu}$ appears in the degenerate hierarchy case.


title: Linking numbers of modular knots

abstract: The modular group PSL(2;Z) acts on the hyperbolic plane HP with quotient the modular surface M, whose unit tangent bundle U is a 3-manifold homeomorphic to the complement of the trefoil knot in the 3-sphere. The hyperbolic conjugacy classes of PSL(2;Z) correspond to the closed oriented geodesics in M. Those lift to the periodic orbits for the geodesic flow in U, which define the modular knots. The linking numbers between modular knots and the trefoil is well understood. Indeed, Etienne Ghys showed in 2006 that they are given by the Rademacher invariant of the corresponding conjugacy classes. The Rademacher function is a homogeneous quasi-morphism of PSL(2;Z) which he had recognised with Jean Barge in 1992 as half the primitive of the bounded euler class. This shed light on the 1987 work of Michael Atiyah concerning the logarithm of the Dedekind eta function which identified it with no less than that six other important functions appearing in diverse areas of mathematics. We are concerned with the linking numbers between modular knots and derive several formulae with arithmetical, combinatorial, topological and group theoretical flavours. In particular we associate to a pair of modular knots a function defined on the character variety of PSL(2;Z), whose limit at the boundary point recovers their linking number. Moreover, we show that the linking number with a modular knot minus that with its inverse yields a homogeneous quasi-morphism on the modular group, and how to extract a free basis out of these. For this we prove that the linking pairing is non degenerate.


title: MINION: a Large-Scale and Diverse Dataset for Multilingual Event Detection

abstract: Event Detection (ED) is the task of identifying and classifying trigger words of event mentions in text. Despite considerable research efforts in recent years for English text, the task of ED in other languages has been significantly less explored. Switching to non-English languages, important research questions for ED include how well existing ED models perform on different languages, how challenging ED is in other languages, and how well ED knowledge and annotation can be transferred across languages. To answer those questions, it is crucial to obtain multilingual ED datasets that provide consistent event annotation for multiple languages. There exist some multilingual ED datasets; however, they tend to cover a handful of languages and mainly focus on popular ones. Many languages are not covered in existing multilingual ED datasets. In addition, the current datasets are often small and not accessible to the public. To overcome those shortcomings, we introduce a new large-scale multilingual dataset for ED (called MINION) that consistently annotates events for 8 different languages; 5 of them have not been supported by existing multilingual datasets. We also perform extensive experiments and analysis to demonstrate the challenges and transferability of ED across languages in MINION that in all call for more research effort in this area.


title: Fast model averaging via buffered states and first-order accelerated optimization algorithms

abstract: In this letter, we study the problem of accelerating reaching average consensus over connected graphs in a discrete-time communication setting. Literature has shown that consensus algorithms can be accelerated by increasing the graph connectivity or optimizing the weights agents place on the information received from their neighbors. Here, instead of altering the communication graph, we investigate two methods that use buffered states to accelerate reaching average consensus over a given graph. In the first method, we study how convergence rate of the well-known first-order Laplacian average consensus algorithm changes when agreement feedback is generated from buffered states. For this study, we obtain a sufficient condition on the ranges of buffered state that leads to faster convergence. In the second proposed method, we show how the average consensus problem can be cast as a convex optimization problem and solved by first-order accelerated optimization algorithms for strongly-convex cost functions. We construct an accelerated average consensus algorithm using the so-called Triple Momentum optimization algorithm. The first approach requires less global knowledge for choosing the step size, whereas the second one converges faster in our numerical results by using extra information from the graph topology. We demonstrate our results by implementing the proposed algorithms in a Gaussian Mixture Model (GMM) estimation problem used in sensor networks.


title: A $\mathrm{GL}(\mathbb{F}_q)$-compatible Hopf algebra of unitriangular class functions

abstract: This paper constructs a novel Hopf algebra $\mathsf{cf}(\mathrm{UT}_{\bullet})$ on the class functions of the unipotent upper triangular groups $\mathrm{UT}_{n}(\mathbb{F}_{q})$ over a finite field. This construction is representation theoretic in nature and uses the machinery of Hopf monoids in the category of vector species. In contrast with a similar known construction, this Hopf algebra has the property that induction to the finite general linear group induces a homomorphism to Zelevinsky's Hopf algebra of $\mathrm{GL}_{n}(\mathbb{F}_{q})$ class functions. Furthermore, $\mathsf{cf}(\mathrm{UT}_{\bullet})$ contains a Hopf subalgebra which is isomorphic to a known combiantorial Hopf algebra, previously used to prove a conjecture about chromatic quasisymmetric functions. Some additional Hopf algebraic properties are also established.


title: Inverse Kernel Decomposition

abstract: The state-of-the-art dimensionality reduction approaches largely rely on complicated optimization procedures. On the other hand, closed-form approaches requiring merely eigen-decomposition do not have enough sophistication and nonlinearity. In this paper, we propose a novel nonlinear dimensionality reduction method -- Inverse Kernel Decomposition (IKD) -- based on an eigen-decomposition of the sample covariance matrix of data. The method is inspired by Gaussian process latent variable models (GPLVMs) and has comparable performance with GPLVMs. To deal with very noisy data with weak correlations, we propose two solutions -- blockwise and geodesic -- to make use of locally correlated data points and provide better and numerically more stable latent estimations. We use synthetic datasets and four real-world datasets to show that IKD is a better dimensionality reduction method than other eigen-decomposition-based methods, and achieves comparable performance against optimization-based methods with faster running speeds. Open-source IKD implementation in Python can be accessed at this \url{https://github.com/JerrySoybean/ikd}.


title: Feature-aggregated spatiotemporal spine surface estimation for wearable patch ultrasound volumetric imaging

abstract: Clear identification of bone structures is crucial for ultrasound-guided lumbar interventions, but it can be challenging due to the complex shapes of the self-shadowing vertebra anatomy and the extensive background speckle noise from the surrounding soft tissue structures. Therefore, we propose to use a patch-like wearable ultrasound solution to capture the reflective bone surfaces from multiple imaging angles and create 3D bone representations for interventional guidance. In this work, we will present our method for estimating the vertebra bone surfaces by using a spatiotemporal U-Net architecture learning from the B-Mode image and aggregated feature maps of hand-crafted filters. The methods are evaluated on spine phantom image data collected by our proposed miniaturized wearable "patch" ultrasound device, and the results show that a significant improvement on baseline method can be achieved with promising accuracy. Equipped with this surface estimation framework, our wearable ultrasound system can potentially provide intuitive and accurate interventional guidance for clinicians in augmented reality setting.


title: JSRNN: Joint Sampling and Reconstruction Neural Networks for High Quality Image Compressed Sensing

abstract: Most Deep Learning (DL) based Compressed Sensing (DCS) algorithms adopt a single neural network for signal reconstruction, and fail to jointly consider the influences of the sampling operation for reconstruction. In this paper, we propose unified framework, which jointly considers the sampling and reconstruction process for image compressive sensing based on well-designed cascade neural networks. Two sub-networks, which are the sampling sub-network and the reconstruction sub-network, are included in the proposed framework. In the sampling sub-network, an adaptive full connected layer instead of the traditional random matrix is used to mimic the sampling operator. In the reconstruction sub-network, a cascade network combining stacked denoising autoencoder (SDA) and convolutional neural network (CNN) is designed to reconstruct signals. The SDA is used to solve the signal mapping problem and the signals are initially reconstructed. Furthermore, CNN is used to fully recover the structure and texture features of the image to obtain better reconstruction performance. Extensive experiments show that this framework outperforms many other state-of-the-art methods, especially at low sampling rates.


title: Thompson Sampling for High-Dimensional Sparse Linear Contextual Bandits

abstract: We consider the stochastic linear contextual bandit problem with high-dimensional features. We analyze the Thompson sampling algorithm using special classes of sparsity-inducing priors (e.g., spike-and-slab) to model the unknown parameter and provide a nearly optimal upper bound on the expected cumulative regret. To the best of our knowledge, this is the first work that provides theoretical guarantees of Thompson sampling in high-dimensional and sparse contextual bandits. For faster computation, we use variational inference instead of Markov Chain Monte Carlo (MCMC) to approximate the posterior distribution. Extensive simulations demonstrate the improved performance of our proposed algorithm over existing ones.


title: Using dynamic circles and squares to visualize spatio-temporal variation

abstract: Visualizations such as bar charts, scatter plots, and objects on geographical maps often convey critical information, including exact and relative numeric values, using shapes. The choice of shape and method of encoding information is often arbitrarily, or based on convention. However, past studies have shown that the human eye can be fooled by visual representations. The Ebbinghaus illusion demonstrates that the perceived relative sizes of shapes depends on their configuration, which in turn can affect judgements, especially in visualizations like proportional symbol maps. In this study we evaluate the effects of varying the type of shapes and metrics for encoding data in visual representations on a spatio-temporal map interface. We find that some combinations of shape and metric are more conducive to accurate human judgements than others, and provide recommendations for applying these findings in future visualization designs.


title: Dynamical downfolding for localized quantum states

abstract: We introduce an approach to treat localized correlated electronic states in the otherwise weakly correlated host medium. Here, the environment is dynamically downfolded on the correlated subspace. It is captured via renormalization of one and two quasiparticle interaction terms which are evaluated using many-body perturbation theory. We outline the strategy on how to take the dynamical effects into account by going beyond the static limit approximation. Further, we introduce an efficient stochastic implementation that enables treating the host environment with a large number of electrons at a minimal computational cost. For small explicitly correlated subspace, the dynamical effects are critical. We demonstrate the methodology by reproducing optical excitations in the negatively charged NV center defect in diamond, that are in excellent agreement with experimental results.