title: Fast Spline Trajectory Planning: Minimum Snap and Beyond

abstract: In this paper, we study spline trajectory generation via the solution of two optimisation problems: (i) a quadratic program (QP) with linear equality constraints and (ii) a nonlinear and nonconvex optimisation program. We propose an efficient algorithm to solve (i), which we then leverage to use in an iterative algorithm to solve (ii). Both the first algorithm and each iteration of the second algorithm have linear computational complexity in the number of spline segments. The scaling of each algorithm is such that we are able to solve the two problems faster than state-of-the-art methods and in times amenable to real-time trajectory generation requirements. The trajectories we generate are applicable to differentially flat systems, a broad class of mechanical systems, which we demonstrate by planning trajectories for a quadrotor.


title: Characterization of Thermal Infrared Dust Emission and Refinements to the Nucleus Properties of Centaur 29P/Schwassmann-Wachmann 1

abstract: We present analyses of Spitzer observations of 29P/Schwassmann-Wachmann 1 using 16 $\mu$m IRS "blue" peak-up (PU) and 24 $\mu$m and 70 $\mu$m MIPS images obtained on UT 2003 November 23 and 24 that characterize the Centaur's large-grain (10-100 $\mu$m) dust coma during a time of non-outbursting "quiescent" activity. Estimates of $\epsilon f \rho$ for each band (16 $\mu$m (2600 $\pm$ 43 cm), 24 $\mu$m (5800 $\pm$ 63 cm), and 70 $\mu$m (1800 $\pm$ 900 cm)) follow the trend between nucleus size vs. $\epsilon f \rho$ that was observed for the WISE/NEOWISE comet ensemble. A coma model was used to derive a dust production rate in the range of 50-100 kg/s. For the first time, a color temperature map of SW1's coma was constructed using the 16 $\mu$m and 24 $\mu$m imaging data. With peaks at $\sim$ 140K, this map implies that coma water ice grains should be slowly sublimating and producing water gas in the coma. We analyzed the persistent 24 $\mu$m "wing" (a curved southwestern coma) feature at 352,000 km (90$''$) from the nucleus attributed by Stansberry et al. (2004) to nucleus rotation and instead propose that it is largely created by solar radiation pressure and gravity acting on micron sized grains. We performed coma removal to the 16 $\mu$m PU image in order to refine the nucleus' emitted thermal flux. A new application of the Near Earth Asteroid Thermal Model (NEATM; Harris 1998) at five wavelengths (5.730 $\mu$m, 7.873 $\mu$m, 15.80 $\mu$m, 23.68 $\mu$m, and 71.42 $\mu$m) was then used to refine SW1's effective radius measurement to $R = 32.3 \pm 3.1$ km and infrared beaming parameter to $\eta = 1.1 \pm 0.2$, respectively.


title: Phase Diagram of YbZnGaO4 in Applied Magnetic Field

abstract: Recently, Yb-based triangular lattice antiferromagnets have garnered significant interest as possible quantum spin liquid candidates. One example is YbMgGaO4, which showed many promising spin liquid features, but also possesses a high degree of disorder owing to site-mixing between the non-magnetic cations. To further elucidate the role of chemical disorder and to explore the phase diagram of these materials in applied field, we present neutron scattering and sensitive magnetometry measurements of the closely related compound, YbZnGaO4. Our results suggest a difference in magnetic anisotropy between the two compounds, and we use key observations of the magnetic phase crossover to motivate an exploration of the field- and exchange parameter-dependent phase diagram, providing an expanded view of the available magnetic states in applied field. This enriched map of the phase space serves as a basis to restrict the values of parameters describing the magnetic Hamiltonian with broad application to recently discovered related materials.


title: Optomechanical Synchronization across Multi-Octaves Frequency Spans

abstract: Experimental exploration of synchronization in scalable oscillator micro systems has unfolded a deeper understanding of networks, collective phenomena, and signal processing. Cavity optomechanical devices have played an important role in this scenario, with the perspective of bridging optical and radio frequencies through nonlinear classical and quantum synchronization concepts. In its simplest form, synchronization occurs when an oscillator is entrained by a signal with frequency nearby the oscillator's tone, and becomes increasingly challenging as their frequency detuning increases. Here, we experimentally demonstrate entrainment of a silicon-nitride optomechanical oscillator driven up to the fourth harmonic of its 32 MHz fundamental frequency. Exploring this effect, we also experimentally demonstrate for the first time a purely optomechanical RF frequency divider, where we performed frequency division up to a 4:1 ratio, i.e., from 128 MHz to 32 MHz. Further developments could harness these effects towards frequency synthesizers, phase-sensitive amplification and nonlinear sensing.


title: Aggregate Cyber-Risk Management in the IoT Age: Cautionary Statistics for (Re)Insurers and Likes

abstract: In this paper, we provide (i) a rigorous general theory to elicit conditions on (tail-dependent) heavy-tailed cyber-risk distributions under which a risk management firm might find it (non)sustainable to provide aggregate cyber-risk coverage services for smart societies, and (ii)a real-data driven numerical study to validate claims made in theory assuming boundedly rational cyber-risk managers, alongside providing ideas to boost markets that aggregate dependent cyber-risks with heavy-tails.To the best of our knowledge, this is the only complete general theory till date on the feasibility of aggregate cyber-risk management.


title: Intensity Harmonization for Airborne LiDAR

abstract: Constructing a point cloud for a large geographic region, such as a state or country, can require multiple years of effort. Often several vendors will be used to acquire LiDAR data, and a single region may be captured by multiple LiDAR scans. A key challenge is maintaining consistency between these scans, which includes point density, number of returns, and intensity. Intensity in particular can be very different between scans, even in areas that are overlapping. Harmonizing the intensity between scans to remove these discrepancies is expensive and time consuming. In this paper, we propose a novel method for point cloud harmonization based on deep neural networks. We evaluate our method quantitatively and qualitatively using a high quality real world LiDAR dataset. We compare our method to several baselines, including standard interpolation methods as well as histogram matching. We show that our method performs as well as the best baseline in areas with similar intensity distributions, and outperforms all baselines in areas with different intensity distributions. Source code is available at https://github.com/mvrl/lidar-harmonization .


title: Real-time Deep Dynamic Characters

abstract: We propose a deep videorealistic 3D human character model displaying highly realistic shape, motion, and dynamic appearance learned in a new weakly supervised way from multi-view imagery. In contrast to previous work, our controllable 3D character displays dynamics, e.g., the swing of the skirt, dependent on skeletal body motion in an efficient data-driven way, without requiring complex physics simulation. Our character model also features a learned dynamic texture model that accounts for photo-realistic motion-dependent appearance details, as well as view-dependent lighting effects. During training, we do not need to resort to difficult dynamic 3D capture of the human; instead we can train our model entirely from multi-view video in a weakly supervised manner. To this end, we propose a parametric and differentiable character representation which allows us to model coarse and fine dynamic deformations, e.g., garment wrinkles, as explicit space-time coherent mesh geometry that is augmented with high-quality dynamic textures dependent on motion and view point. As input to the model, only an arbitrary 3D skeleton motion is required, making it directly compatible with the established 3D animation pipeline. We use a novel graph convolutional network architecture to enable motion-dependent deformation learning of body and clothing, including dynamics, and a neural generative dynamic texture model creates corresponding dynamic texture maps. We show that by merely providing new skeletal motions, our model creates motion-dependent surface deformations, physically plausible dynamic clothing deformations, as well as video-realistic surface textures at a much higher level of detail than previous state of the art approaches, and even in real-time.


title: NeuroXplorer 1.0: An Extensible Framework for Architectural Exploration with Spiking Neural Networks

abstract: Recently, both industry and academia have proposed many different neuromorphic architectures to execute applications that are designed with Spiking Neural Network (SNN). Consequently, there is a growing need for an extensible simulation framework that can perform architectural explorations with SNNs, including both platform-based design of today's hardware, and hardware-software co-design and design-technology co-optimization of the future. We present NeuroXplorer, a fast and extensible framework that is based on a generalized template for modeling a neuromorphic architecture that can be infused with the specific details of a given hardware and/or technology. NeuroXplorer can perform both low-level cycle-accurate architectural simulations and high-level analysis with data-flow abstractions. NeuroXplorer's optimization engine can incorporate hardware-oriented metrics such as energy, throughput, and latency, as well as SNN-oriented metrics such as inter-spike interval distortion and spike disorder, which directly impact SNN performance. We demonstrate the architectural exploration capabilities of NeuroXplorer through case studies with many state-of-the-art machine learning models.


title: A GPU-based Kalman Filter for Track Fitting

abstract: Computing centres, including those used to process High-Energy Physics data and simulations, are increasingly providing significant fractions of their computing resources through hardware architectures other than x86 CPUs, with GPUs being a common alternative. GPUs can provide excellent computational performance at a good price point for tasks that can be suitably parallelized. Charged particle (track) reconstruction is a computationally expensive component of HEP data reconstruction, and thus needs to use available resources in an efficient way. In this paper, an implementation of Kalman filter-based track fitting using CUDA and running on GPUs is presented. This utilizes the ACTS (A Common Tracking Software) toolkit; an open source and experiment-independent toolkit for track reconstruction. The implementation details and parallelization approach are described, along with the specific challenges for such an implementation. Detailed performance benchmarking results are discussed, which show encouraging performance gains over a CPU-based implementation for representative configurations. Finally, a perspective on the challenges and future directions for these studies is outlined. These include more complex and realistic scenarios which can be studied, and anticipated developments to software frameworks and standards which may open up possibilities for greater flexibility and improved performance.


title: Star cluster survival in dark matter halos: An old cluster in Eridanus II?

abstract: The star formation history and the internal dynamics of Milky Way satellite galaxies are often complicated. In the last years, a substantial fraction of the known faint dwarf satellites have been studied. Some of them show embedded stellar substructures, such as star clusters and even globular star clusters. In this work we study Eridanus II, a dwarf spheroidal satellite which hosts a star cluster, using published and archival data from the Hubble Space Telescope Advanced Camera for Surveys. We employ a Bayesian hierarchical method to infer the star formation history of Eridanus II. We find that the bulk of the stars in Eridanus II are very old ($13.5_{-1}^{+0.5}$ Gyr) and quite metal poor ($Z$\,=\,$0.00001$). We do not find any evidence of the presence of an intermediate age or young population in Eri II. We cannot date the embedded star cluster as a separate entity, but we find it likely that the cluster has a similar age and metallicity as the bulk of the stars in Eri II. The existence of an old star cluster in a dark matter dominated old metal poor dwarf galaxy is of major importance to cast light on the dark matter distribution within dwarf galaxies. The existence of intermediate age stars is required by the recent detection of carbon stars in Eri II. Since no recent star formation is detected, {\it blue-straggler} fusions of lower mass stars are the most likely origin of the carbon star progenitors.


title: Pervasive AI for IoT applications: A Survey on Resource-efficient Distributed Artificial Intelligence

abstract: Artificial intelligence (AI) has witnessed a substantial breakthrough in a variety of Internet of Things (IoT) applications and services, spanning from recommendation systems to robotics control and military surveillance. This is driven by the easier access to sensory data and the enormous scale of pervasive/ubiquitous devices that generate zettabytes (ZB) of real-time data streams. Designing accurate models using such data streams, to predict future insights and revolutionize the decision-taking process, inaugurates pervasive systems as a worthy paradigm for a better quality-of-life. The confluence of pervasive computing and artificial intelligence, Pervasive AI, expanded the role of ubiquitous IoT systems from mainly data collection to executing distributed computations with a promising alternative to centralized learning, presenting various challenges. In this context, a wise cooperation and resource scheduling should be envisaged among IoT devices (e.g., smartphones, smart vehicles) and infrastructure (e.g. edge nodes, and base stations) to avoid communication and computation overheads and ensure maximum performance. In this paper, we conduct a comprehensive survey of the recent techniques developed to overcome these resource challenges in pervasive AI systems. Specifically, we first present an overview of the pervasive computing, its architecture, and its intersection with artificial intelligence. We then review the background, applications and performance metrics of AI, particularly Deep Learning (DL) and online learning, running in a ubiquitous system. Next, we provide a deep literature review of communication-efficient techniques, from both algorithmic and system perspectives, of distributed inference, training and online learning tasks across the combination of IoT devices, edge devices and cloud servers. Finally, we discuss our future vision and research challenges.


title: Towards End-to-End Deep Learning for Autonomous Racing: On Data Collection and a Unified Architecture for Steering and Throttle Prediction

abstract: Deep Neural Networks (DNNs) which are trained end-to-end have been successfully applied to solve complex problems that we have not been able to solve in past decades. Autonomous driving is one of the most complex problems which is yet to be completely solved and autonomous racing adds more complexity and exciting challenges to this problem. Towards the challenge of applying end-to-end learning to autonomous racing, this paper shows results on two aspects: (1) Analyzing the relationship between the driving data used for training and the maximum speed at which the DNN can be successfully applied for predicting steering angle, (2) Neural network architecture and training methodology for learning steering and throttle without any feedback or recurrent connections.


title: Generative Adversarial Networks (GAN) Powered Fast Magnetic Resonance Imaging -- Mini Review, Comparison and Perspectives

abstract: Magnetic Resonance Imaging (MRI) is a vital component of medical imaging. When compared to other image modalities, it has advantages such as the absence of radiation, superior soft tissue contrast, and complementary multiple sequence information. However, one drawback of MRI is its comparatively slow scanning and reconstruction compared to other image modalities, limiting its usage in some clinical applications when imaging time is critical. Traditional compressive sensing based MRI (CS-MRI) reconstruction can speed up MRI acquisition, but suffers from a long iterative process and noise-induced artefacts. Recently, Deep Neural Networks (DNNs) have been used in sparse MRI reconstruction models to recreate relatively high-quality images from heavily undersampled k-space data, allowing for much faster MRI scanning. However, there are still some hurdles to tackle. For example, directly training DNNs based on L1/L2 distance to the target fully sampled images could result in blurry reconstruction because L1/L2 loss can only enforce overall image or patch similarity and does not take into account local information such as anatomical sharpness. It is also hard to preserve fine image details while maintaining a natural appearance. More recently, Generative Adversarial Networks (GAN) based methods are proposed to solve fast MRI with enhanced image perceptual quality. The encoder obtains a latent space for the undersampling image, and the image is reconstructed by the decoder using the GAN loss. In this chapter, we review the GAN powered fast MRI methods with a comparative study on various anatomical datasets to demonstrate the generalisability and robustness of this kind of fast MRI while providing future perspectives.


title: Fair and Truthful Mechanism with Limited Subsidy

abstract: The notion of \emph{envy-freeness} is a natural and intuitive fairness requirement in resource allocation. With indivisible goods, such fair allocations are unfortunately not guaranteed to exist. Classical works have avoided this issue by introducing an additional divisible resource, i.e., money, to subsidize envious agents. In this paper, we aim to design a truthful allocation mechanism of indivisible goods to achieve both fairness and efficiency criteria with a limited amount of subsidy. Following the work of Halpern and Shah, our central question is as follows: to what extent do we need to rely on the power of money to accomplish these objectives? For general valuations, the impossibility theorem of combinatorial auction translates to our setting: even if an arbitrarily large amount of money is available for use, no mechanism can achieve truthfulness, envy-freeness, and utilitarian optimality simultaneously when agents have general monotone submodular valuations. By contrast, we show that, when agents have matroidal valuations, there is a truthful allocation mechanism that achieves envy-freeness and utilitarian optimality by subsidizing each agent with at most $1$, the maximum marginal contribution of each item for each agent. The design of the mechanism rests crucially on the underlying matroidal M-convexity of the Lorenz dominating allocations.


title: Independent Spaces of q-Polymatroids

abstract: This paper is devoted to the study of independent spaces of q-polymatroids. With the aid of an auxiliary q-matroid it is shown that the collection of independent spaces satisfies the same properties as for q-matroids. However, in contrast to q-matroids, the rank value of an independent space does not agree with its dimension. Nonetheless, the rank values of the independent spaces fully determine the q-polymatroid, and this fact can be exploited to derive a cryptomorphism of q-polymatroids. Finally, the notions of minimal spanning spaces, maximally strongly independent spaces, and bases will be elaborated on.


title: DeepRT: A Soft Real Time Scheduler for Computer Vision Applications on the Edge

abstract: The ubiquity of smartphone cameras and IoT cameras, together with the recent boom of deep learning and deep neural networks, proliferate various computer vision driven mobile and IoT applications deployed on the edge. This paper focuses on applications which make soft real time requests to perform inference on their data - they desire prompt responses within designated deadlines, but occasional deadline misses are acceptable. Supporting soft real time applications on a multi-tenant edge server is not easy, since the requests sharing the limited GPU computing resources of an edge server interfere with each other. In order to tackle this problem, we comprehensively evaluate how latency and throughput respond to different GPU execution plans. Based on this analysis, we propose a GPU scheduler, DeepRT, which provides latency guarantee to the requests while maintaining high overall system throughput. The key component of DeepRT, DisBatcher, batches data from different requests as much as possible while it is proven to provide latency guarantee for requests admitted by an Admission Control Module. DeepRT also includes an Adaptation Module which tackles overruns. Our evaluation results show that DeepRT outperforms state-of-the-art works in terms of the number of deadline misses and throughput.


title: The Polarized Image of a Synchrotron Emitting Ring of Gas Orbiting a Black Hole

abstract: Synchrotron radiation from hot gas near a black hole results in a polarized image. The image polarization is determined by effects including the orientation of the magnetic field in the emitting region, relativistic motion of the gas, strong gravitational lensing by the black hole, and parallel transport in the curved spacetime. We explore these effects using a simple model of an axisymmetric, equatorial accretion disk around a Schwarzschild black hole. By using an approximate expression for the null geodesics derived by Beloborodov (2002) and conservation of the Walker-Penrose constant, we provide analytic estimates for the image polarization. We test this model using currently favored general relativistic magnetohydrodynamic simulations of M87*, using ring parameters given by the simulations. For a subset of these with modest Faraday effects, we show that the ring model broadly reproduces the polarimetric image morphology. Our model also predicts the polarization evolution for compact flaring regions, such as those observed from Sgr A* with GRAVITY. With suitably chosen parameters, our simple model can reproduce the EVPA pattern and relative polarized intensity in Event Horizon Telescope images of M87*. Under the physically motivated assumption that the magnetic field trails the fluid velocity, this comparison is consistent with the clockwise rotation inferred from total intensity images.


title: Newton Polyhedra and Whitney Equisingularity for Isolated Determinantal Singularities

abstract: Using Newton polyhedra and non-degeneracy of matrices we present conditions which guarantee the Whitney equisingularity of families of isolated determinantal singularities.


title: Thom's counterexamples for the Steenrod problem

abstract: The present paper deals with integral classes $\xi_p\in H_{2p+1}(L^{2p+1}\times L^{2p+1})$ which are counterexamples for the Steenrod realization problem, where $L^{2p+1}$ is the $(2p+1)$-dimensional lens space and $p\geq 3$ is a prime number. For $p=3$, this is Thom's famous counterexample. We give a geometric description of this class using the theory of stratifolds. As a consequence, we obtain a geometric interpretation of the obstruction to realizability in terms of the Atiyah--Hirzebruch spectral sequence.


title: Bayesian inference of an uncertain generalized diffusion operator

abstract: This paper defines a novel Bayesian inverse problem to infer an infinite-dimensional uncertain operator appearing in a differential equation, whose action on an observable state variable affects its dynamics. Inference is made tractable by parametrizing the operator using its eigendecomposition. The plausibility of operator inference in the sparse data regime is explored in terms of an uncertain, generalized diffusion operator appearing in an evolution equation for a contaminant's transport through a heterogeneous porous medium. Sparse data are augmented with prior information through the imposition of deterministic constraints on the eigendecomposition and the use of qualitative information about the system in the definition of the prior distribution. Limited observations of the state variable's evolution are used as data for inference, and the dependence on the solution of the inverse problem is studied as a function of the frequency of observations, as well as on whether or not the data is collected as a spatial or time series.


title: Herd immunity for Covid-19 in homogeneous populations

abstract: In this note we estimate herd immunity levels for the Covid-19 epidemic using a standard SEIR system that models the disease dynamics in homogeneous populations. The results obtained are indicative of values between 80% and 90% for unprotected, fully susceptible populations. Basic protective measures such as hand hygiene and mask wearing may be effective to lower herd immunity levels down to values between 50% and 60% of the total population.


title: The diverse nature of small-scale turbulence

abstract: The self-similar Richardson cascade admits two logically possible scenarios of small-scale turbulence at high Reynolds numbers. In the first scenario, eddies' population densities vary as a function of eddies' scales. As a result, one or a few eddy types dominate at small scales, and small-scale turbulence lacks diversity. In the second scenario, eddies' population densities are scale-invariant across the inertial range, resulting in small-scale diversity. That is, there are as many types of eddies at the small scales as at the large scales. In this letter, we measure eddies' population densities in three-dimensional isotropic turbulence and determine the nature of small-scale turbulence. The result shows that eddies' population densities are scale-invariant.


title: The $o\tau$-continuous, Lebesgue, KB, and Levi operators between vector lattices and topological vector spaces

abstract: We investigate the $o\tau$-continuous/bounded/compact and Lebesgue operators from vector lattices to topological vector spaces; the KB operators between locally solid lattices and topological vector spaces; and the Levi operators from locally solid lattices to vector lattices. The main idea of operator versions of notions related to vector lattices lies in redistributing topological and order properties of a topological vector lattice between the domain and range of an operator under investigation. Domination properties for these classes of operators are studied.


title: Training Structured Mechanical Models by Minimizing Discrete Euler-Lagrange Residual

abstract: Model-based paradigms for decision-making and control are becoming ubiquitous in robotics. They rely on the ability to efficiently learn a model of the system from data. Structured Mechanical Models (SMMs) are a data-efficient black-box parameterization of mechanical systems, typically fit to data by minimizing the error between predicted and observed accelerations or next states. In this work, we propose a methodology for fitting SMMs to data by minimizing the discrete Euler-Lagrange residual. To study our methodology, we fit models to joint-angle time-series from undamped and damped double-pendulums, studying the quality of learned models fit to data with and without observation noise. Experiments show that our methodology learns models that are better in accuracy to those of the conventional schemes for fitting SMMs. We identify use cases in which our method is a more appropriate methodology. Source code for reproducing the experiments is available at https://github.com/sisl/delsmm.


title: Constructions of BiHom-X algebras and bimodules of some BiHom-dialgebras

abstract: The aim of this paper is to introduce and to give some constructions results of BiHom-X algebras by using Yau's twisting, Rota Baxter and Some elements of centroids. Next, we define the bimodules of BiHom-left symmetric dialgebras, BiHom-associative dialgebras and BiHom-tridendriform algebras. A sequence of this kind of bimodules can be constructed. Their matched pairs are also introduced and related relevant properties are given.


title: Phenomenological model for the third-harmonic magnetic response of superconductors: application to Sr$_{2}$RuO$_{4}$

abstract: We employ the phenomenological Lawrence-Doniach model to compute the contributions of the superconducting fluctuations to the third-harmonic magnetic response, denoted here by $\overline{M_{3}}$, which can be measured in a precise way using ac magnetic fields and lock-in techniques. We show that, in an intermediate temperature regime, this quantity behaves as the third-order nonlinear susceptibility, which shows a power-law dependence with the reduced temperature $\epsilon=\frac{T-T_{c}}{T_{c}}$ as $\epsilon^{-5/2}$. Very close to $T_{c}$, however, $\overline{M_{3}}$ saturates due to the nonzero amplitude of the ac field. We compare our theoretical results with experimental data for three conventional superconductors -- lead, niobium, and vanadium -- and for the unconventional superconductor Sr$_{2}$RuO$_{4}$ (SRO). We find good agreement between theory and experiment for the elemental superconductors, although the theoretical values for the critical field systematically deviate from the experimental ones. In the case of SRO, however, the phenomenological model completely fails to describe the data, as the third-harmonic response remains sizable over a much wider reduced temperature range compared to Pb, Nb, and V. We show that an inhomogeneous distribution of $T_{c}$ can partially account for this discrepancy, since regions with a locally higher $T_{c}$ contribute to the fluctuation $\overline{M_{3}}$ significantly more than regions with the "nominal" $T_{c}$ of the clean system. However, the exponential temperature dependence of $\overline{M_{3}}$ first reported in Ref. [3] is not captured by the model with inhomogeneity. We conclude that, while inhomogeneity is an important ingredient to understand the superconducting fluctuations of SRO and other perovskite superconductors, additional effects may be at play, such as non-Gaussian fluctuations or rare-region effects.


title: On the minimal extension and structure of weakly group-theoretical braided fusion categories

abstract: We show that any slightly degenerate weakly group-theoretical fusion category admits a minimal non-degenerate extension. Let $d$ be a positive square-free integer, given a weakly group-theoretical non-degenerate fusion category $\mathcal{C}$, assume that $\text{FPdim}(\mathcal{C})=nd$ and $(n,d)=1$. If $(\text{FPdim}(X)^2,d)=1$ for all simple objects $X$ of $\mathcal{C}$, then we show that $\mathcal{C}$ contains a non-degenerate fusion subcategory $\mathcal{C}(\mathbb{Z}_d,q)$. In particular, we obtain that integral fusion categories of FP-dimensions $p^md$ such that $\mathcal{C}'\subseteq \text{sVec}$ are nilpotent and group-theoretical, where $p$ is a prime and $(p,d)=1$.


title: Vulnerability of Blockchain Technologies to Quantum Attacks

abstract: Quantum computation represents a threat to many cryptographic protocols in operation today. It has been estimated that by 2035, there will exist a quantum computer capable of breaking the vital cryptographic scheme RSA2048. Blockchain technologies rely on cryptographic protocols for many of their essential sub-routines. Some of these protocols, but not all, are open to quantum attacks. Here we analyze the major blockchain-based cryptocurrencies deployed today -- including Bitcoin, Ethereum, Litecoin and ZCash, and determine their risk exposure to quantum attacks. We finish with a comparative analysis of the studied cryptocurrencies and their underlying blockchain technologies and their relative levels of vulnerability to quantum attacks.


title: Real-time Face Mask Detection in Video Data

abstract: In response to the ongoing COVID-19 pandemic, we present a robust deep learning pipeline that is capable of identifying correct and incorrect mask-wearing from real-time video streams. To accomplish this goal, we devised two separate approaches and evaluated their performance and run-time efficiency. The first approach leverages a pre-trained face detector in combination with a mask-wearing image classifier trained on a large-scale synthetic dataset. The second approach utilizes a state-of-the-art object detection network to perform localization and classification of faces in one shot, fine-tuned on a small set of labeled real-world images. The first pipeline achieved a test accuracy of 99.97% on the synthetic dataset and maintained 6 FPS running on video data. The second pipeline achieved a mAP(0.5) of 89% on real-world images while sustaining 52 FPS on video data. We have concluded that if a larger dataset with bounding-box labels can be curated, this task is best suited using object detection architectures such as YOLO and SSD due to their superior inference speed and satisfactory performance on key evaluation metrics.


title: Discretization of 4d Poincar\'e BF theory: from groups to 2-groups

abstract: We study the discretization of a Poincar\'e/Euclidean BF theory. Upon the addition of a boundary term, this theory is equivalent to the BFCG theory defined in terms of the Poincar\'e/Euclidean 2-group. At an intermediate step in the discretization, we note that there are multiple options for how to proceed. One option brings us back to recovering the discrete variables and phase space of BF theory. Another option allows us to rediscover the phase space related to the G-networks given in [2]. Indeed, our main result is that we are now able to relate the continuum fields with the discrete variables in [2]. This relation is important to determine how to implement the simplicity constraints to recover gravity using the BFCG action. In fact, we show that such relation is not as simple as in the BF discretization: the discretized variables on the triangles actually depend on several of the continuum fields instead of solely the continuum B-field. We also compare and contrast the discretized BF and BFCG models as pairs of dual 2-groups. This work highlights (again) how the choice of boundary term influences the resulting symmetry structure of the \textit{discretized theory} -- and hence ultimately the choice of quantum states.


title: Dynamic Enumeration of Similarity Joins

abstract: This paper considers enumerating answers to similarity-join queries under dynamic updates: Given two sets of $n$ points $A,B$ in $\mathbb{R}^d$, a metric $\phi(\cdot)$, and a distance threshold $r > 0$, report all pairs of points $(a, b) \in A \times B$ with $\phi(a,b) \le r$. Our goal is to store $A,B$ into a dynamic data structure that, whenever asked, can enumerate all result pairs with worst-case delay guarantee, i.e., the time between enumerating two consecutive pairs is bounded. Furthermore, the data structure can be efficiently updated when a point is inserted into or deleted from $A$ or $B$. We propose several efficient data structures for answering similarity-join queries in low dimension. For exact enumeration of similarity join, we present near-linear-size data structures for $\ell_1, \ell_\infty$ metrics with $\log^{O(1)} n$ update time and delay. We show that such a data structure is not feasible for the $\ell_2$ metric for $d \ge 4$. For approximate enumeration of similarity join, where the distance threshold is a soft constraint, we obtain a unified linear-size data structure for $\ell_p$ metric, with $\log^{O(1)} n$ delay and update time. In high dimensions, we present an efficient data structure with worst-case delay-guarantee using locality sensitive hashing (LSH).


title: ExcavatorCovid: Extracting Events and Relations from Text Corpora for Temporal and Causal Analysis for COVID-19

abstract: Timely responses from policy makers to mitigate the impact of the COVID-19 pandemic rely on a comprehensive grasp of events, their causes, and their impacts. These events are reported at such a speed and scale as to be overwhelming. In this paper, we present ExcavatorCovid, a machine reading system that ingests open-source text documents (e.g., news and scientific publications), extracts COVID19 related events and relations between them, and builds a Temporal and Causal Analysis Graph (TCAG). Excavator will help government agencies alleviate the information overload, understand likely downstream effects of political and economic decisions and events related to the pandemic, and respond in a timely manner to mitigate the impact of COVID-19. We expect the utility of Excavator to outlive the COVID-19 pandemic: analysts and decision makers will be empowered by Excavator to better understand and solve complex problems in the future. An interactive TCAG visualization is available at http://afrl402.bbn.com:5050/index.html. We also released a demonstration video at https://vimeo.com/528619007.


title: Calibration of Human Driving Behavior and Preference Using Naturalistic Traffic Data

abstract: Understanding human driving behaviors quantitatively is critical even in the era when connected and autonomous vehicles and smart infrastructure are becoming ever more prevalent. This is particularly so as that mixed traffic settings, where autonomous vehicles and human driven vehicles co-exist, are expected to persist for quite some time. Towards this end it is necessary that we have a comprehensive modeling framework for decision-making within which human driving preferences can be inferred statistically from observed driving behaviors in realistic and naturalistic traffic settings. Leveraging a recently proposed computational framework for smart vehicles in a smart world using multi-agent based simulation and optimization, we first recapitulate how the forward problem of driving decision-making is modeled as a state space model. We then show how the model can be inverted to estimate driver preferences from naturalistic traffic data using the standard Kalman filter technique. We explicitly illustrate our approach using the vehicle trajectory data from Sugiyama experiment that was originally meant to demonstrate how stop-and-go shockwave can arise spontaneously without bottlenecks. Not only the estimated state filter can fit the observed data well for each individual vehicle, the inferred utility functions can also re-produce quantitatively similar pattern of the observed collective behaviors. One distinct advantage of our approach is the drastically reduced computational burden. This is possible because our forward model treats driving decision process, which is intrinsically dynamic with multi-agent interactions, as a sequence of independent static optimization problems contingent on the state with a finite look ahead anticipation. Consequently we can practically sidestep solving an interacting dynamic inversion problem that would have been much more computationally demanding.


title: Quantum Advantage on Proof of Work

abstract: Proof-of-Work (PoW) is a fundamental underlying technology behind most major blockchain cryptocurrencies. It has been previously pointed out that quantum devices provide a computational advantage in performing PoW in the context of Bitcoin. Here we make the case that this quantum advantage extends not only to all existing PoW mechanisms, but to any possible PoW as well. This has strong consequences regarding both quantum-based attacks on the integrity of the entirety of the blockchain, as well as more legitimate uses of quantum computation for the purpose of mining Bitcoin and other cryptocurrencies. For the first case, we estimate when these quantum attacks will become feasible, for various cryptocurrencies, and discuss the impact of such attacks. For the latter, we derive a precise formula to calculate the economic incentive for switching to quantum-based cryptocurrency miners. Using this formula, we analyze several test scenarios, and conclude that investing in quantum hardware for cryptocurrency mining has the potential to pay off immensely.


title: On the Spatial and Temporal Order of Convergence of Hyperbolic PDEs

abstract: In this work, we determine the full expression for the global truncation error of hyperbolic partial differential equations (PDEs). In particular, we use theoretical analysis and symbolic algebra to find exact expressions for the coefficients of the generic global truncation error. Our analysis is valid for any hyperbolic PDE, be it linear or non-linear, and employing finite difference, finite volume, or finite element discretization in space, and advanced in time with a predictor-corrector, multistep, or a deferred correction method, belonging to the Method of Lines. Furthermore, we discuss the practical implications of this analysis. If we employ a stable numerical scheme and the orders of accuracy of the global solution error and the global truncation error agree, we make the following asymptotic observations: (a) the order of convergence at constant ratio of $\Delta t$ to $\Delta x$ is governed by the minimum of the orders of the spatial and temporal discretizations, and (b) convergence cannot even be guaranteed under only spatial or temporal refinement. An implication of (a) is that it is impractical to invest in a time-stepping method of order higher than the spatial discretization. In addition to (b), we demonstrate that under certain circumstances, the error can even monotonically increase with refinement only in space or only in time, and explain why this phenomenon occurs. To verify our theoretical findings, we conduct convergence studies of linear and non-linear advection equations using finite difference and finite volume spatial discretizations, and predictor-corrector and multistep time-stepping methods. Finally, we study the effect of slope limiters and monotonicity-preserving strategies on the order of accuracy.


title: TransHash: Transformer-based Hamming Hashing for Efficient Image Retrieval

abstract: Deep hamming hashing has gained growing popularity in approximate nearest neighbour search for large-scale image retrieval. Until now, the deep hashing for the image retrieval community has been dominated by convolutional neural network architectures, e.g. \texttt{Resnet}\cite{he2016deep}. In this paper, inspired by the recent advancements of vision transformers, we present \textbf{Transhash}, a pure transformer-based framework for deep hashing learning. Concretely, our framework is composed of two major modules: (1) Based on \textit{Vision Transformer} (ViT), we design a siamese vision transformer backbone for image feature extraction. To learn fine-grained features, we innovate a dual-stream feature learning on top of the transformer to learn discriminative global and local features. (2) Besides, we adopt a Bayesian learning scheme with a dynamically constructed similarity matrix to learn compact binary hash codes. The entire framework is jointly trained in an end-to-end manner.~To the best of our knowledge, this is the first work to tackle deep hashing learning problems without convolutional neural networks (\textit{CNNs}). We perform comprehensive experiments on three widely-studied datasets: \textbf{CIFAR-10}, \textbf{NUSWIDE} and \textbf{IMAGENET}. The experiments have evidenced our superiority against the existing state-of-the-art deep hashing methods. Specifically, we achieve 8.2\%, 2.6\%, 12.7\% performance gains in terms of average \textit{mAP} for different hash bit lengths on three public datasets, respectively.


title: Deep learning-based tumor segmentation on digital images of histopathology slides for microdosimetry applications

abstract: $\bf{Purpose:}$ The goal of this study was (i) to use artificial intelligence to automate the traditionally labor-intensive process of manual segmentation of tumor regions in pathology slides performed by a pathologist and (ii) to validate the use of a well-known and readily available deep learning architecture. Automation will reduce the human error involved in manual delineation, increase efficiency, and result in accurate and reproducible segmentation. This advancement will alleviate the bottleneck in the workflow in clinical and research applications due to a lack of pathologist time. Our application is patient-specific microdosimetry and radiobiological modeling, which builds on the contoured pathology slides. $\bf{Methods:}$ A U-Net architecture was used to segment tumor regions in pathology core biopsies of lung tissue with adenocarcinoma stained using hematoxylin and eosin. A pathologist manually contoured the tumor regions in 56 images with binary masks for training. Overlapping patch extraction with various patch sizes and image downsampling were investigated individually. Data augmentation and 8-fold cross-validation were used. $\bf{Results:}$ The U-Net achieved accuracy of 0.91$\pm$0.06, specificity of 0.90$\pm$0.08, sensitivity of 0.92$\pm$0.07, and precision of 0.8$\pm$0.1. The F1/DICE score was 0.85$\pm$0.07, with a segmentation time of 3.24$\pm$0.03 seconds per image, achieving a 370$\pm$3 times increased efficiency over manual segmentation. In some cases, the U-Net correctly delineated the tumor's stroma from its epithelial component in regions that were classified as tumor by the pathologist. $\bf{Conclusion:}$ The U-Net architecture can segment images with a level of efficiency and accuracy that makes it suitable for tumor segmentation of histopathological images in fields such as radiotherapy dosimetry, specifically in the subfields of microdosimetry.


title: Inconsequential results on the Merino-Welsh conjecture for Tutte polynomials

abstract: The Merino-Welsh conjectures say that subject to conditions, there is an inequality among the Tutte-polynomial evaluations $T(M;2,0)$, $T(M;0,2)$, and $T(M;1,1)$. We present three results on a Merino-Welsh conjecture. These results are "inconsequential" in the sense that although they imply a version of the conjecture for many matroids, they seem to be dead ends.


title: First implications of Tibet AS$_\gamma$ data for heavy dark matter

abstract: Extensive air shower detectors of gamma rays in the sub-PeV energy region provide a new and relatively unexplored window for dark matter searches. Here we derive some implications of the recently published Tibet AS$_\gamma$ data for decaying dark matter candidates. The available spectral information is already useful in obtaining competitive constraints, surpassing existing limits above 10 PeV mass for hadronic or massive boson final states. This is particularly true if accounting for a benchmark astrophysical background of Galactic cosmic rays in the (0.1-1) PeV range. By relying on the arrival distribution of the photons, we show that significantly better sensitivity can be attained, comparable or better than IceCube also for most leptonic final states. Full data exploitation requires however further information disclosure.


title: GALA: Greedy ComputAtion for Linear Algebra in Privacy-Preserved Neural Networks

abstract: Machine Learning as a Service (MLaaS) is enabling a wide range of smart applications on end devices. However, privacy-preserved computation is still expensive. Our investigation has found that the most time-consuming component of the HE-based linear computation is a series of Permutation (Perm) operations that are imperative for dot product and convolution in privacy-preserved MLaaS. To this end, we propose GALA: Greedy computAtion for Linear Algebra in privacy-preserved neural networks, which views the HE-based linear computation as a series of Homomorphic Add, Mult and Perm operations and chooses the least expensive operation in each linear computation step to reduce the overall cost. GALA makes the following contributions: (1) It introduces a row-wise weight matrix encoding and combines the share generation that is needed for the GC-based nonlinear computation, to reduce the Perm operations for the dot product; (2) It designs a first-Add-second-Perm approach (named kernel grouping) to reduce Perm operations for convolution. As such, GALA efficiently reduces the cost for the HE-based linear computation, which is a critical building block in almost all of the recent frameworks for privacy-preserved neural networks, including GAZELLE (Usenix Security'18), DELPHI (Usenix Security'20), and CrypTFlow2 (CCS'20). With its deep optimization of the HE-based linear computation, GALA can be a plug-and-play module integrated into these systems to further boost their efficiency. Our experiments show that it achieves a significant speedup up to 700x for the dot product and 14x for the convolution computation under different data dimensions. Meanwhile, GALA demonstrates an encouraging runtime boost by 2.5x, 2.7x, 3.2x, 8.3x, 7.7x, and 7.5x over GAZELLE and 6.5x, 6x, 5.7x, 4.5x, 4.2x, and 4.1x over CrypTFlow2, on AlexNet, VGG, ResNet-18, ResNet-50, ResNet-101, and ResNet-152, respectively.


title: Lesion Segmentation and RECIST Diameter Prediction via Click-driven Attention and Dual-path Connection

abstract: Measuring lesion size is an important step to assess tumor growth and monitor disease progression and therapy response in oncology image analysis. Although it is tedious and highly time-consuming, radiologists have to work on this task by using RECIST criteria (Response Evaluation Criteria In Solid Tumors) routinely and manually. Even though lesion segmentation may be the more accurate and clinically more valuable means, physicians can not manually segment lesions as now since much more heavy laboring will be required. In this paper, we present a prior-guided dual-path network (PDNet) to segment common types of lesions throughout the whole body and predict their RECIST diameters accurately and automatically. Similar to [1], a click guidance from radiologists is the only requirement. There are two key characteristics in PDNet: 1) Learning lesion-specific attention matrices in parallel from the click prior information by the proposed prior encoder, named click-driven attention; 2) Aggregating the extracted multi-scale features comprehensively by introducing top-down and bottom-up connections in the proposed decoder, named dual-path connection. Experiments show the superiority of our proposed PDNet in lesion segmentation and RECIST diameter prediction using the DeepLesion dataset and an external test set. PDNet learns comprehensive and representative deep image features for our tasks and produces more accurate results on both lesion segmentation and RECIST diameter prediction.


title: Who Are I: Time Inconsistency and Intrapersonal Conflict and Reconciliation

abstract: Time inconsistency is prevalent in dynamic choice problems: a plan of actions to be taken in the future that is optimal for an agent today may not be optimal for the same agent in the future. If the agent is aware of this intra-personal conflict but unable to commit herself in the future to following the optimal plan today, the rational strategy for her today is to reconcile with her future selves, namely to correctly anticipate her actions in the future and then act today accordingly. Such a strategy is named intra-personal equilibrium and has been studied since as early as in the 1950s. A rigorous treatment in continuous-time settings, however, had not been available until a decade ago. Since then, the study on intra-personal equilibrium for time-inconsistent problems in continuous time has grown rapidly. In this chapter, we review the classical results and some recent development in this literature.


title: Contemporary COBOL: Developers' Perspectives on Defects and Defect Location

abstract: Mainframe systems are facing a critical shortage of developer workforce as the current generation of COBOL developers retires. Furthermore, due to the limited availability of public COBOL resources, entry-level developers, who assume the mantle of legacy COBOL systems maintainers, face significant difficulties during routine maintenance tasks, such as code comprehension and defect location. While we made substantial advances in the field of software maintenance for modern programming languages yearly, mainframe maintenance has received limited attention. With this study, we aim to direct the attention of researchers and practitioners towards investigating and addressing challenges associated with mainframe development. Specifically, we explore the scope of defects affecting COBOL systems and defect location strategies commonly followed by COBOL developers and compare them with the modern programming language counterparts. To this end, we surveyed 30 COBOL and 74 modern Programming Language (PL) developers to understand the differences in defects and defect location strategies employed by the two groups. Our preliminary results show that: (1) major defect categories affecting the COBOL ecosystem are different than defects encountered in modern PL software projects; (2) the most challenging defect types in COBOL are also the ones that occur most frequently; and (3) COBOL and modern PL developers follow similar strategies to locate defective code.


title: What happens once an accelerating observer has detected a Rindler particle?

abstract: In a seminal paper, Unruh and Wald found that the detection of a right Rindler particle by a linearly uniformly accelerated detector coupled to a Klein-Gordon field in the Minkowski vacuum leads to the creation of a Minkowski particle from the inertial viewpoint. In this paper, we revisit the framework studied by Unruh and Wald, but now consider in addition what happens once the particle has been measured somewhere in the right Rindler wedge. From an orthodox point of view, the change in the field state induced by the measurement is non-local and occurs both in the left and right Rindler wedges. If one takes semiclassical gravity seriously in this context, this seems to open the possiblity for designing superluminal communication protocols between two spacelike separated observers confined to the right and left Rindler wedges respectively. We discuss the possible ways in which physics could prevent such measurement-induced, faster-than-light signaling protocols.


title: Thermal conductivity minimum of graded superlattices due to phonon localization

abstract: The Anderson localization of thermal phonons has been shown only in few nano-structures with strong random disorder by the exponential decay of transmission to zero and a thermal conductivity maximum when increasing system length. In this work, we present a path to demonstrate the phonon localization with distinctive features in graded superlattices with short-range order and long-range disorder. A thermal conductivity minimum with system length appears due to the exponential decay of transmission to a non-zero constant, which is a feature of partial phonon localization caused by the moderate disorder. We provide clear evidence of localization through the combined analysis of the participation ratio, transmission, and real-space phonon number density distribution based on our quantum transport simulation. The present work would promote heat conduction engineering by localization via the wave nature of phonons.


title: The Complexity of Symmetry Breaking in Massive Graphs

abstract: The goal of this paper is to understand the complexity of symmetry breaking problems, specifically maximal independent set (MIS) and the closely related $\beta$-ruling set problem, in two computational models suited for large-scale graph processing, namely the $k$-machine model and the graph streaming model. We present a number of results. For MIS in the $k$-machine model, we improve the $\tilde{O}(m/k^2 + \Delta/k)$-round upper bound of Klauck et al. (SODA 2015) by presenting an $\tilde{O}(m/k^2)$-round algorithm. We also present an $\tilde{\Omega}(n/k^2)$ round lower bound for MIS, the first lower bound for a symmetry breaking problem in the $k$-machine model. For $\beta$-ruling sets, we use hierarchical sampling to obtain more efficient algorithms in the $k$-machine model and also in the graph streaming model. More specifically, we obtain a $k$-machine algorithm that runs in $\tilde{O}(\beta n\Delta^{1/\beta}/k^2)$ rounds and, by using a similar hierarchical sampling technique, we obtain one-pass algorithms for both insertion-only and insertion-deletion streams that use $O(\beta \cdot n^{1+1/2^{\beta-1}})$ space. The latter result establishes a clear separation between MIS, which is known to require $\Omega(n^2)$ space (Cormode et al., ICALP 2019), and $\beta$-ruling sets, even for $\beta = 2$. Finally, we present an even faster 2-ruling set algorithm in the $k$-machine model, one that runs in $\tilde{O}(n/k^{2-\epsilon} + k^{1-\epsilon})$ rounds for any $\epsilon$, $0 \le \epsilon \le 1$.


title: Curvatures of Stiefel manifolds with deformation metrics

abstract: We compute curvatures of a family of tractable metrics on Stiefel manifolds, introduced recently by H{\"u}per, Markina and Silva Leite, which includes the well-known embedded and canonical metrics on Stiefel manifolds as special cases. The metrics could be identified with the Cheeger deformation metrics. We identify parameter values in the family to make a Stiefel manifold an Einstein manifold and show Stiefel manifolds always carry an Einstein metric. We analyze the sectional curvature range and identify the parameter range where the manifold has non-negative sectional curvature. We provide the exact sectional curvature range when the number of columns in a Stiefel matrix is $2$, and a conjectural range for other cases. We derive the formulas from two approaches, one from a global curvature formula derived in our recent work, another using curvature formulas for left-invariant metrics. The second approach leads to curvature formulas for Cheeger deformation metrics on normal homogeneous spaces.

categories: math.DG cs.CV cs.LG cs.SY eess.SY math.OC


title: Twin-domain formation in epitaxial triangular lattice delafossites

abstract: Twin domains are often found as structural defects in symmetry mismatched epitaxial thin films. The delafossite ABO2, which has a rhombohedral structure, is a good example that often forms twin domains. Although bulk metallic delafossites are known to be the most conducting oxides, the high conductivity is yet to be realized in thin film forms. Suppressed conductivity found in thin films is mainly caused by the formation of twin domains, and their boundaries can be a source of scattering centers for charge carriers. To overcome this challenge, the underlying mechanism for their formation must be understood, so that such defects can be controlled and eliminated. Here, we report the origin of structural twins formed in a CuCrO2 delafossite thin film on a substrate with hexagonal or triangular symmetries. A robust heteroepitaxial relationship is found for the delafossite film with the substrate, and the surface termination turns out to be critical to determine and control the domain structure of epitaxial delafossites. Based on such discoveries, we also demonstrate a twin-free epitaxial thin films grown on high-miscut substrates. This finding provides an important synthesis strategy for growing single domain delafossite thin films and can be applied to other delafossites for epitaxial synthesis of high-quality thin films.


title: Acoustic Scene Classification Using Multichannel Observation with Partially Missing Channels

abstract: Sounds recorded with smartphones or IoT devices often have partially unreliable observations caused by clipping, wind noise, and completely missing parts due to microphone failure and packet loss in data transmission over the network. In this paper, we investigate the impact of the partially missing channels on the performance of acoustic scene classification using multichannel audio recordings, especially for a distributed microphone array. Missing observations cause not only losses of time-frequency and spatial information on sound sources but also a mismatch between a trained model and evaluation data. We thus investigate how a missing channel affects the performance of acoustic scene classification in detail. We also propose simple data augmentation methods for scene classification using multichannel observations with partially missing channels and evaluate the scene classification performance using the data augmentation methods.


title: Strain-induced atomic-scale building blocks for ferromagnetism in epitaxial LaCoO3

abstract: The origin of strain-induced ferromagnetism, which is robust regardless of the type and degree of strain in LaCoO3 (LCO) thin films, is enigmatic despite intensive research efforts over the past decade. Here, by combining scanning transmission electron microscopy with ab initio density functional theory plus U calculations, we report that the ferromagnetism does not emerge directly from the strain itself, but rather from the creation of compressed structural units within ferroelastically formed twin-wall domains. The compressed structural units are magnetically active with the rocksalt-type high-spin/low-spin order. Our study highlights that the ferroelastic nature of ferromagnetic structural units is important for understanding the intriguing ferromagnetic properties in LCO thin films.


title: Improved Surrogate Modeling of Fluid Dynamics with Physics-Informed Neural Networks

abstract: Physics-Informed Neural Networks (PINNs) have recently shown great promise as a way of incorporating physics-based domain knowledge, including fundamental governing equations, into neural network models for many complex engineering systems. They have been particularly effective in the area of inverse problems, where boundary conditions may be ill-defined, and data-absent scenarios, where typical supervised learning approaches will fail. Here, we further explore the use of this modeling methodology to surrogate modeling of a fluid dynamical system, and demonstrate additional undiscussed and interesting advantages of such a modeling methodology over conventional data-driven approaches: 1) improving the model's predictive performance even with incomplete description of the underlying physics; 2) improving the robustness of the model to noise in the dataset; 3) reduced effort to convergence during optimization for a new, previously unseen scenario by transfer optimization of a pre-existing model. Hence, we noticed the inclusion of a physics-based regularization term can substantially improve the equivalent data-driven surrogate model in many substantive ways, including an order of magnitude improvement in test error when the dataset is very noisy, and a 2-3x improvement when only partial physics is included. In addition, we propose a novel transfer optimization scheme for use in such surrogate modeling scenarios and demonstrate an approximately 3x improvement in speed to convergence and an order of magnitude improvement in predictive performance over conventional Xavier initialization for training of new scenarios.


title: Encoder Fusion Network with Co-Attention Embedding for Referring Image Segmentation

abstract: Recently, referring image segmentation has aroused widespread interest. Previous methods perform the multi-modal fusion between language and vision at the decoding side of the network. And, linguistic feature interacts with visual feature of each scale separately, which ignores the continuous guidance of language to multi-scale visual features. In this work, we propose an encoder fusion network (EFN), which transforms the visual encoder into a multi-modal feature learning network, and uses language to refine the multi-modal features progressively. Moreover, a co-attention mechanism is embedded in the EFN to realize the parallel update of multi-modal features, which can promote the consistent of the cross-modal information representation in the semantic space. Finally, we propose a boundary enhancement module (BEM) to make the network pay more attention to the fine structure. The experiment results on four benchmark datasets demonstrate that the proposed approach achieves the state-of-the-art performance under different evaluation metrics without any post-processing.


title: CUAB: Convolutional Uncertainty Attention Block Enhanced the Chest X-ray Image Analysis

abstract: In recent years, convolutional neural networks (CNNs) have been successfully implemented to various image recognition applications, such as medical image analysis, object detection, and image segmentation. Many studies and applications have been working on improving the performance of CNN algorithms and models. The strategies that aim to improve the performance of CNNs can be grouped into three major approaches: (1) deeper and wider network architecture, (2) automatic architecture search, and (3) convolutional attention block. Unlike approaches (1) and (2), the convolutional attention block approach is more flexible with lower cost. It enhances the CNN performance by extracting more efficient features. However, the existing attention blocks focus on enhancing the significant features, which lose some potential features in the uncertainty information. Inspired by the test time augmentation and test-time dropout approaches, we developed a novel convolutional uncertainty attention block (CUAB) that can leverage the uncertainty information to improve CNN-based models. The proposed module discovers potential information from the uncertain regions on feature maps in computer vision tasks. It is a flexible functional attention block that can be applied to any position in the convolutional block in CNN models. We evaluated the CUAB with notable backbone models, ResNet and ResNeXt, on a medical image segmentation task. The CUAB achieved a dice score of 73% and 84% in pneumonia and pneumothorax segmentation, respectively, thereby outperforming the original model and other notable attention approaches. The results demonstrated that the CUAB can efficiently utilize the uncertainty information to improve the model performance.


title: Counting closed geodesics on rank one manifolds without focal points

abstract: In this article, we consider a closed rank one Riemannian manifold $M$ without focal points. Let $P(t)$ be the set of free-homotopy classes containing a closed geodesic on $M$ with length at most $t$, and $\# P(t)$ its cardinality. We obtain the following Margulis-type asymptotic estimates: \[\lim_{t\to \infty}\#P(t)/\frac{e^{ht}}{ht}=1\] where $h$ is the topological entropy of the geodesic flow. In the appendix, we also show that the unique measure of maximal entropy of the geodesic flow has the Bernoulli property.


title: Electromagnetic properties of neutrinos from scattering on bound electrons in atom

abstract: We consider the effects of bound atomic electrons scattered by solar neutrinos due to the electromagnetic properties of neutrinos. This necessiate considering the recoil of atomic nucleus, which should be considered in the momentum conservation, but that effect to the energy conservation is negligible. This effect changes the kinematic behavior of the scattered electron compared to that scattered on free electrons. We apply this effect to the recent XENON1T data, but the bounds obtained from this is not very restrictive. We obtained the bounds: the (transition) magnetic moment $|f_{\alpha\beta}|\le 0.86\times 10^{-7}$ (times the electron Bohr magneton) and the charge radius $|\tilde{r}|< 4.30\times 10^{-17\,}{\rm cm}$. For a non-vanishing millicharge ($\varepsilon$), the allowed bound is shown in the $\tilde{r}^2-\varepsilon$ plane.


title: Promoting global stability in data-driven models of quadratic nonlinear dynamics

abstract: Modeling realistic fluid and plasma flows is computationally intensive, motivating the use of reduced-order models for a variety of scientific and engineering tasks. However, it is challenging to characterize, much less guarantee, the global stability (i.e., long-time boundedness) of these models. The seminal work of Schlegel and Noack (JFM, 2015) provided a theorem outlining necessary and sufficient conditions to ensure global stability in systems with energy-preserving, quadratic nonlinearities, with the goal of evaluating the stability of projection-based models. In this work, we incorporate this theorem into modern data-driven models obtained via machine learning. First, we propose that this theorem should be a standard diagnostic for the stability of projection-based and data-driven models, examining the conditions under which it holds. Second, we illustrate how to modify the objective function in machine learning algorithms to promote globally stable models, with implications for the modeling of fluid and plasma flows. Specifically, we introduce a modified "trapping SINDy" algorithm based on the sparse identification of nonlinear dynamics (SINDy) method. This method enables the identification of models that, by construction, only produce bounded trajectories. The effectiveness and accuracy of this approach are demonstrated on a broad set of examples of varying model complexity and physical origin, including the vortex shedding in the wake of a circular cylinder.


title: Joint Registration and Segmentation via Multi-Task Learning for Adaptive Radiotherapy of Prostate Cancer

abstract: Medical image registration and segmentation are two of the most frequent tasks in medical image analysis. As these tasks are complementary and correlated, it would be beneficial to apply them simultaneously in a joint manner. In this paper, we formulate registration and segmentation as a joint problem via a Multi-Task Learning (MTL) setting, allowing these tasks to leverage their strengths and mitigate their weaknesses through the sharing of beneficial information. We propose to merge these tasks not only on the loss level, but on the architectural level as well. We studied this approach in the context of adaptive image-guided radiotherapy for prostate cancer, where planning and follow-up CT images as well as their corresponding contours are available for training. The study involves two datasets from different manufacturers and institutes. The first dataset was divided into training (12 patients) and validation (6 patients), and was used to optimize and validate the methodology, while the second dataset (14 patients) was used as an independent test set. We carried out an extensive quantitative comparison between the quality of the automatically generated contours from different network architectures as well as loss weighting methods. Moreover, we evaluated the quality of the generated deformation vector field (DVF). We show that MTL algorithms outperform their Single-Task Learning (STL) counterparts and achieve better generalization on the independent test set. The best algorithm achieved a mean surface distance of $1.06 \pm 0.3$ mm, $1.27 \pm 0.4$ mm, $0.91 \pm 0.4$ mm, and $1.76 \pm 0.8$ mm on the validation set for the prostate, seminal vesicles, bladder, and rectum, respectively. The high accuracy of the proposed method combined with the fast inference speed, makes it a promising method for automatic re-contouring of follow-up scans for adaptive radiotherapy.


title: Photoinduced anomalous Hall effect in the interacting Haldane model: targeting topological states with pump pulses

abstract: We investigate the nonequilibrium dynamics of the spinless Haldane model with nearest-neighbor interactions on the honeycomb lattice by employing an unbiased numerical method. In this system, a first-order transition from the Chern insulator (CI) at weak coupling to the charge-density-wave (CDW) phase at strong coupling can be characterized by a level crossing of the lowest energy levels. Here we show that adiabatically following the eigenstates across this level crossing, their Chern numbers are preserved, leading to the identification of a topologically-nontrivial low-energy excited state in the CDW regime. By promoting a resonant energy excitation via an ultrafast circularly polarized pump pulse, we find that the system acquires a non-vanishing Hall response as a result of the large overlap enhancement between the time-dependent wave-function and the topologically non-trivial excited state. This is suggestive of a photoinduced topological phase transition via unitary dynamics, despite a proper definition of the Chern number remaining elusive for an out-of-equilibrium interacting system. We contrast these results with more common quench protocols, where such features are largely absent in the dynamics even if the post-quench Hamiltonian displays a topologically nontrivial ground state.


title: PingAn-VCGroup's Solution for ICDAR 2021 Competition on Scientific Table Image Recognition to Latex

abstract: This paper presents our solution for the ICDAR 2021 Competition on Scientific Table Image Recognition to LaTeX. This competition has two sub-tasks: Table Structure Reconstruction (TSR) and Table Content Reconstruction (TCR). We treat both sub-tasks as two individual image-to-sequence recognition problems. We leverage our previously proposed algorithm MASTER \cite{lu2019master}, which is originally proposed for scene text recognition. We optimize the MASTER model from several perspectives: network structure, optimizer, normalization method, pre-trained model, resolution of input image, data augmentation, and model ensemble. Our method achieves 0.7444 Exact Match and 0.8765 Exact Match @95\% on the TSR task, and obtains 0.5586 Exact Match and 0.7386 Exact Match 95\% on the TCR task.


title: Localized modes and dark solitons sustained by nonlinear defects

abstract: Dark solitons and localized defect modes against periodic backgrounds are considered in arrays of waveguides with defocusing Kerr nonlinearity constituting a nonlinear lattice. Bright defect modes are supported by local increase of the nonlinearity, while dark defect modes are supported by a local decrease of the nonlinearity. Dark solitons exist for both types of the defect, although in the case of weak nonlinearity they feature side bright humps making the total energy propagating through the system larger than the energy transferred by the constant background. All considered defect modes are found stable. Dark solitons are characterized by relatively narrow windows of stability. Interactions of unstable dark solitons with bright and dark modes are described.


title: PingAn-VCGroup's Solution for ICDAR 2021 Competition on Scientific Literature Parsing Task B: Table Recognition to HTML

abstract: This paper presents our solution for ICDAR 2021 competition on scientific literature parsing taskB: table recognition to HTML. In our method, we divide the table content recognition task into foursub-tasks: table structure recognition, text line detection, text line recognition, and box assignment.Our table structure recognition algorithm is customized based on MASTER [1], a robust image textrecognition algorithm. PSENet [2] is used to detect each text line in the table image. For text linerecognition, our model is also built on MASTER. Finally, in the box assignment phase, we associatedthe text boxes detected by PSENet with the structure item reconstructed by table structure prediction,and fill the recognized content of the text line into the corresponding item. Our proposed methodachieves a 96.84% TEDS score on 9,115 validation samples in the development phase, and a 96.32%TEDS score on 9,064 samples in the final evaluation phase.


title: Analytical results for a minimalist thermal diode

abstract: We consider a system consisting of two interacting classical particles, each one subject to an on-site potential and to a Langevin thermal bath. We analytically calculate the heat current that can be established through the system when the bath temperatures are different, for weak nonlinear forces. We explore the conditions under which the diode effect emerges when inverting the temperature difference. Despite the simplicity of this two-particle diode, an intricate dependence on the system parameters is put in evidence. Moreover, behaviors reported for long chains of particles can be extracted, for instance, the dependence of the flux with the interfacial stiffness and type of forces present, as well as the dependencies on the temperature required for rectification. These analytical results can be a tool to foresee the distinct role that diverse types of nonlinearity and asymmetry play in thermal conduction and rectification.


title: Preference learning along multiple criteria: A game-theoretic perspective

abstract: The literature on ranking from ordinal data is vast, and there are several ways to aggregate overall preferences from pairwise comparisons between objects. In particular, it is well known that any Nash equilibrium of the zero sum game induced by the preference matrix defines a natural solution concept (winning distribution over objects) known as a von Neumann winner. Many real-world problems, however, are inevitably multi-criteria, with different pairwise preferences governing the different criteria. In this work, we generalize the notion of a von Neumann winner to the multi-criteria setting by taking inspiration from Blackwell's approachability. Our framework allows for non-linear aggregation of preferences across criteria, and generalizes the linearization-based approach from multi-objective optimization. From a theoretical standpoint, we show that the Blackwell winner of a multi-criteria problem instance can be computed as the solution to a convex optimization problem. Furthermore, given random samples of pairwise comparisons, we show that a simple plug-in estimator achieves near-optimal minimax sample complexity. Finally, we showcase the practical utility of our framework in a user study on autonomous driving, where we find that the Blackwell winner outperforms the von Neumann winner for the overall preferences.


title: Associativity of fusion products of $C_1$-cofinite ${\mathbb N}$-gradable modules of vertex operator algebra

abstract: We prove an associative law of the fusion products $\boxtimes$ of $C_1$-cofinite ${\mathbb N}$-gradable modules for a vertex operator algebra $V$. To be more precise, for $C_1$-cofinite ${\mathbb N}$-gradable $V$-modules $A,B,C$ and their fusion products $(A\!\boxtimes\! B, {\cal Y}^{AB})$, $((A\!\boxtimes\! B)\!\boxtimes\! C, {\cal Y}^{(AB)C})$, $(B\!\boxtimes\! C, {\cal Y}^{BC})$, $(A\!\boxtimes\! (B\!\boxtimes\! C),{\cal Y}^{A(BC)})$ with logarithmic intertwining operators ${\cal Y}^{AB},\ldots,{\cal Y}^{A(BC)}$ satisfying the universal properties for ${\mathbb N}$-gradable modules, we prove that four-point correlation functions $\langle \theta, {\cal Y}^{A(BC)}(v,x){\cal Y}^{BC}(u,y)w\rangle$ and $\langle \theta', {\cal Y}^{(AB)C}({\cal Y}^{AB}(v,x-y)u,y)w\rangle$ are locally normally convergent over $\{(x,y)\in {\mathbb C}^2 \mid 0\!<\!|x\!-\!y|\!<\!|y|\!<\!|x|\}$. We then take their respective principal branches $\tilde{F}(\langle \theta,{\cal Y}^{A(BC)}(v,x){\cal Y}^{BC}(u,y)w\rangle)$ and $\tilde{F}(\langle \theta,{\cal Y}^{(AB)C}({\cal Y}^{AB)}(v,x-y)u,y)w\rangle)$ on ${\cal D}^2\!=\!\{(x,y)\in {\mathbb C}^2 \mid 0\!<\!|x\!-\!y|\!<\!|y|\!<\!|x|, \mbox{ and } x,y,x\!-\!y\not\in {\mathbb R}^{\leq 0}\}$ and then show that there is an isomorphism $\phi_{[AB]C}:(A\boxtimes B)\boxtimes C \to A\boxtimes (B\boxtimes C)$ such that $$ \widetilde{F}(\langle \theta, {\cal Y}^{A(BC)}(v,x){\cal Y}^{BC}(u,y)w\rangle) =\tilde{F}(\langle \phi_{[AB]C}^{\ast}(\theta), {\cal Y}^{(AB)C}({\cal Y}^{AB}(v,x-y)u,y)w)\rangle $$ on ${\cal D}^2$ for $\theta\in (A\boxtimes (B\boxtimes C))^{\vee}$, $v\in A$, $u\in B$, and $w\in C$, where $W^{\vee}$ denotes the contragredient module of $W$ and $\phi_{[AB]C}^{\ast}$ denotes the dual of $\phi_{[AB]C}$. We also prove the pentagon identity.


title: Deep Learning for Needle Detection in a Cannulation Simulator

abstract: Cannulation for hemodialysis is the act of inserting a needle into a surgically created vascular access (e.g., an arteriovenous fistula) for the purpose of dialysis. The main risk associated with cannulation is infiltration, the puncture of the wall of the vascular access after entry, which can cause medical complications. Simulator-based training allows clinicians to gain cannulation experience without putting patients at risk. In this paper, we propose to use deep-learning-based techniques for detecting, based on video, whether the needle tip is in or has infiltrated the simulated fistula. Three categories of deep neural networks are investigated in this work: modified pre-trained models based on VGG-16 and ResNet-50, light convolutional neural networks (light CNNs), and convolutional recurrent neural networks (CRNNs). CRNNs consist of convolutional layers and a long short-term memory (LSTM) layer. A data set of cannulation experiments was collected and analyzed. The results show that both the light CNN and the CRNN achieve better performance than the pre-trained baseline models. The CRNN was implemented in real time on commodity hardware for use in the cannulation simulator, and the performance was verified. Deep-learning video analysis is a viable method for detecting needle state in a low cost cannulation simulator. Our data sets and code are released at https://github.com/axin233/DL_for_Needle_Detection_Cannulation


title: Two-Stage Stochastic Optimization via Primal-Dual Decomposition and Deep Unrolling

abstract: We consider a two-stage stochastic optimization problem, in which a long-term optimization variable is coupled with a set of short-term optimization variables in both objective and constraint functions. Despite that two-stage stochastic optimization plays a critical role in various engineering and scientific applications, there still lack efficient algorithms, especially when the long-term and short-term variables are coupled in the constraints. To overcome the challenge caused by tightly coupled stochastic constraints, we first establish a two-stage primal-dual decomposition (PDD) method to decompose the two-stage problem into a long-term problem and a family of short-term subproblems. Then we propose a PDD-based stochastic successive convex approximation (PDD-SSCA) algorithmic framework to find KKT solutions for two-stage stochastic optimization problems. At each iteration, PDD-SSCA first runs a short-term sub-algorithm to find stationary points of the short-term subproblems associated with a mini-batch of the state samples. Then it constructs a convex surrogate for the long-term problem based on the deep unrolling of the short-term sub-algorithm and the back propagation method. Finally, the optimal solution of the convex surrogate problem is solved to generate the next iterate. We establish the almost sure convergence of PDD-SSCA and customize the algorithmic framework to solve two important application problems. Simulations show that PDD-SSCA can achieve superior performance over existing solutions.


title: Parameterizing the Angular Distribution of Emission: A Model for TOF-PET Low-Count Reconstruction

abstract: Low-count reconstruction remains a challenge for Positron Emission Tomography (PET) even with the recent progress in time-of-flight (TOF) resolution. In this context, the bias between the acquired histogram, consisting of low values or zeros, and the expected histogram, obtained from the forward projector, is propagated to the image, resulting in a biased reconstruction. This situation could be exacerbated with finer resolution of the TOF information, which further sparsifies the acquired histogram. We propose a new approach to circumvent this limitation of the classical reconstruction model. It consists of extending the description of the reconstruction scheme to also explicitly include the projection domain by Parameterizing the Angular Distribution of Emission (PADE). This parametrization has greater degrees of freedom than the log-likelihood model, which can not be harnessed in classical circumstances. We hypothesize that with ultra-fast TOF this new approach would not only be viable for low-count reconstruction, but also more adequate than the classical reconstruction model. As a proof of concept, an implementation of this approach is compared to the log-likelihood model by investigating two-dimensional simulations of a hot spots phantom. The proposed model achieves similar contrast recovery coefficients as MLEM except for the smallest structures where the low-count nature of the simulations makes it difficult to draw conclusions. However, this new model seems to converge toward less noisy solutions than MLEM. These results suggest that the PADE approach has potential for low-count reconstruction with ultra-fast TOF.


title: Hennessy-Milner Properties via Topological Compactness

abstract: We give Hennessy-Milner classes for intuitionistic, dual-intuitionistic and bi-intuitionistic logic interpreted in intuitionistic Kripke models, and generalise these results to modal (dual- and bi-)intuitionistic logics. Our main technical tools are image-compact and pre-image-compact relations that provide a semantical description of modal saturation properties.


title: Identity testing under label mismatch

abstract: Testing whether the observed data conforms to a purported model (probability distribution) is a basic and fundamental statistical task, and one that is by now well understood. However, the standard formulation, identity testing, fails to capture many settings of interest; in this work, we focus on one such natural setting, identity testing under promise of permutation. In this setting, the unknown distribution is assumed to be equal to the purported one, up to a relabeling (permutation) of the model: however, due to a systematic error in the reporting of the data, this relabeling may not be the identity. The goal is then to test identity under this assumption: equivalently, whether this systematic labeling error led to a data distribution statistically far from the reference model.


title: Magnetic Properties of NbSi2N4, VSi2N4, and VSi2P4 Monolayers

abstract: The recent demonstration of MoSi2N4 and its exceptional stability to air, water, acid, and heat has generated intense interest in this family of two-dimensional (2D) materials. Among these materials, NbSi2N4, VSi2N4, and VSi2P4 are semiconducting, easy-plane ferromagnets with negligible in-plane magnetic anisotropy. They thus satisfy a necessary condition for exhibiting a dissipationless spin superfluid mode. The Curie temperatures of monolayer VSi2P4 and VSi2N4 are determined to be above room temperature based on Monte Carlo and density functional theory calculations. The magnetic moments of VSi2N4 can be switched from in-plane to out-of-plane by applying tensile biaxial strain or electron doping.


title: Performance Analysis of Free-space Quantum Key Distribution Using Multiple Spatial Modes

abstract: In the diffraction-limited near-field propagation regime, free-space optical quantum key distribution (QKD) systems can employ multiple spatial modes to improve their key rate. Here, we analyze QKD using the non-orthogonal flat-top focused beams. Although they suffer from a rate penalty, their ease of implementation makes them an attractive alternative to the well-studied orthonormal Laguerre-Gauss (LG) modes. Indeed, in the presence of turbulence, the non-orthogonal modes may achieve higher QKD rate than the LG modes.


title: Function4D: Real-time Human Volumetric Capture from Very Sparse Consumer RGBD Sensors

abstract: Human volumetric capture is a long-standing topic in computer vision and computer graphics. Although high-quality results can be achieved using sophisticated off-line systems, real-time human volumetric capture of complex scenarios, especially using light-weight setups, remains challenging. In this paper, we propose a human volumetric capture method that combines temporal volumetric fusion and deep implicit functions. To achieve high-quality and temporal-continuous reconstruction, we propose dynamic sliding fusion to fuse neighboring depth observations together with topology consistency. Moreover, for detailed and complete surface generation, we propose detail-preserving deep implicit functions for RGBD input which can not only preserve the geometric details on the depth inputs but also generate more plausible texturing results. Results and experiments show that our method outperforms existing methods in terms of view sparsity, generalization capacity, reconstruction quality, and run-time efficiency.


title: SemperFi: A Spoofer Eliminating GPS Receiver for UAVs

abstract: It is well-known that GPS is vulnerable to signal spoofing attacks. Although several spoofing detection techniques exist, they are incapable of mitigation and recovery from stealthy attackers. In this work, we present SemperFi, a single antenna GPS receiver capable of tracking legitimate GPS satellite signals and estimating the true location even during a spoofing attack. The main challenge in building SemperFi is, unlike most wireless systems where \emph{the data} contained in the wireless signals is important, GPS relies on the time of arrival (ToA) of satellite signals. SemperFi is capable of distinguishing spoofing signals and recovering legitimate GPS signals that are even completely overshadowed by a strong adversary. We exploit the short-term stability of inertial sensors to identify the spoofing signal and extend the successive interference cancellation algorithm to preserve the legitimate signal's ToA. We implement SemperFi in GNSS-SDR, an open-source software-defined GNSS receiver, and evaluate its performance using UAV simulators, real drones, a variety of real-world GPS datasets, and various embedded platforms. Our evaluation results indicate that in many scenarios, SemperFi can identify adversarial peaks by executing flight patterns that are less than 50 m long and recover the true location within 10 seconds (Jetson Xavier). We show that our receiver is secure against stealthy attackers who exploit inertial sensor errors and execute seamless takeover attacks. We design SemperFi as a pluggable module capable of generating a spoofer-free GPS signal for processing on any commercial-off-the-shelf GPS receiver available today. Finally, we release our implementation to the community for usage and further research.


title: Higher-Order Model Checking Step by Step

abstract: We show a new simple algorithm that solves the model-checking problem for recursion schemes: check whether the tree generated by a given higher-order recursion scheme is accepted by a given alternating parity automaton. The algorithm amounts to a procedure that transforms a recursion scheme of order $n$ to a recursion scheme of order $n-1$, preserving acceptance, and increasing the size only exponentially. After repeating the procedure $n$ times, we obtain a recursion scheme of order $0$, for which the problem boils down to solving a finite parity game. Since the size grows exponentially at each step, the overall complexity is $n$-EXPTIME, which is known to be optimal. More precisely, the transformation is linear in the size of the recursion scheme, assuming that the arity of employed nonterminals and the size of the automaton are bounded by a constant; this results in an FPT algorithm for the model-checking problem. Our transformation is a generalization of a previous transformation of the author (2020), working for reachability automata in place of parity automata. The step-by-step approach can be opposed to previous algorithms solving the considered problem "in one step", being compulsorily more complicated.


title: Observations of near-perfect nonclassical correlation using coherent light

abstract: Complementarity theory is the essence of the Copenhagen interpretation. Since the Hanbury Brown and Twiss experiments, the particle nature of photons has been intensively studied for various quantum phenomena such as anticorrelation and Bell inequality violation in terms of two-photon correlation. Regarding the fundamental question on these quantum features, however, no clear answer exists for how to generate such an entanglement photon pair and what causes the maximum correlation between them. Here, we experimentally demonstrate the physics of anticorrelation on a beam splitter using sub-Poisson distributed coherent photons, where a particular photon number is post-selected using a multiphoton resolving coincidence measurement technique. According to Born rule regarding self-interference in an interferometric scheme, a photon does not interact with others, but can interfere by itself. This is the heart of anticorrelation, where a particular phase relation between paired photons is unveiled for anticorrelation, satisfying the complementarity theory of quantum mechanics.


title: Quantum speed-up in collisional battery charging

abstract: We present a collision model for the charging of a quantum battery by identical nonequilibrium qubit units. When the units are prepared in a mixture of energy eigenstates, the energy gain in the battery can be described by a classical random walk, where both average energy and variance grow linearly with time. Conversely, when the qubits contain quantum coherence, interference effects buildup in the battery and lead to a faster spreading of the energy distribution, reminiscent of a quantum random walk. This can be exploited for faster and more efficient charging of a battery initialized in the ground state. Specifically, we show that coherent protocols can yield higher charging power than any possible incoherent strategy, demonstrating a quantum speed-up at the level of a single battery. Finally, we characterize the amount of extractable work from the battery through the notion of ergotropy.


title: Tree Path Minimum Query Oracle via Boruvka Trees

abstract: Tree path minimum query problem is a fundamental problem while processing trees, and is used widely in minimum spanning tree verification and randomized minimum spanning tree algorithms. In this paper, we study the possibility of building an oracle in advance, which is able to answer the queries efficiently. We present an algorithm based on Boruvka trees. Our algorithm is the first to achieve a near-optimal bound on query time, while matching the currently optimal trade-off between construction time and the number of comparisons required at query. Particularly, in order to answer each query within $2k$ comparisons, our algorithm requires $O(n \log \lambda_k(n))$ time and space to construct the oracle, and the oracle can answer queries in $O(k + \log \lambda_k(n))$ time. Here $\lambda_k(n)$ is the inverse of the Ackermann function along the $k$-th column. This algorithm not only is simpler than the previous ones, but also gives a completely different method of solving this problem.


title: An efficient phase-field method for turbulent multiphase flows

abstract: With the aim of efficiently simulating three-dimensional multiphase turbulent flows with a phase-field method, we propose a new discretization scheme for the biharmonic term (the 4th-order derivative term) of the Cahn-Hilliard equation. This novel scheme can significantly reduce the computational cost while retaining the same accuracy as the original procedure. Our phase-field method is built on top of a direct numerical simulation solver, named AFiD (www.afid.eu) and open-sourced by our research group. It relies on a pencil distributed parallel strategy and a FFT-based Poisson solver. To deal with large density ratios between the two phases, a pressure split method [1] has been applied to the Poisson solver. To further reduce computational costs, we implement a multiple-resolution algorithm which decouples the discretizations for the Navier-Stokes equations and the scalar equation: while a stretched wall-resolving grid is used for the Navier-Stokes equations, for the Cahn-Hilliard equation we use a fine uniform mesh. The present method shows excellent computational performance for large-scale computation: on meshes up to 8 billion nodes and 3072 CPU cores, a multiphase flow needs only slightly less than 1.5 times the CPU time of the single-phase flow solver on the same grid. The present method is validated by comparing the results to previous studies for the cases of drop deformation in shear flow, including the convergence test with mesh refinement, and breakup of a rising buoyant bubble with density ratio up to 1000. Finally, we simulate the breakup of a big drop and the coalescence of O(10^3) drops in turbulent Rayleigh-B\'enard convection at a Rayleigh number of $10^8$, observing good agreement with theoretical results.


title: Necessity of selfconsistent calculations for the electromagnetic field in probing the nuclear symmetry energy using pion observables in heavy-ion collisions

abstract: Within an isospin- and momentum-dependent transport model, we investigate the necessity of selfconsistent calculations for the electromagnetic field in probing the nuclear symmetry energy using pion observables in heavy-ion collisions at intermediate energies. To this end, we perform the $^{96}$Ru + $^{96}$Ru collisions at 400 MeV/nucleon with two calculations scenarios for the electromagnetic field including the selfconsistent calculation and the most used Li\'{e}nard-Wiechert formula, while the latter is a simplified one of the complete Li\'{e}nard-Wiechert formula by neglecting the radiation field for practical calculations in heavy-ion collisions at intermediate and/or relativistic energies. As a comparison, we also consider the static Coulomb field formula for calculations of the electromagnetic field in heavy-ion collisions. It is shown that the most used simplified Li\'{e}nard-Wiechert formula is not enough for the electromagnetic field calculation because the absent radiation field in this formula also affects significantly the charged pions as well as their $\pi^{-}/\pi^{+}$ ratio. Moreover, we also examine effects of the electromagnetic field in these scenarios on the double $\pi^{-}/\pi^{+}$ ratio of two isobar reaction systems of $^{96}$Ru + $^{96}$Ru and $^{96}$Zr + $^{96}$Zr at 400 MeV/nucleon. It is shown that the double $\pi^{-}/\pi^{+}$ ratio of two reactions tends to be less affected by the electromagnetic field calculation scenario and thus can still be an effective probe of the nuclear symmetry energy in heavy-ion collisions. Therefore, according to these findings, it is suggested that the selfconsistent calculation for the electromagnetic field should be carefully taken into account when using the pion observables to probe the nuclear symmetry energy in heavy-ion collisions.