title: Secure Video Streaming Using Dedicated Hardware

abstract: Purpose: The purpose of this article is to present a system that enhances the security, efficiency, and reconfigurability of an Internet-of-Things (IoT) system used for surveillance and monitoring. Methods: A Multi-Processor System-On-Chip (MPSoC) composed of Central Processor Unit (CPU) and Field-Programmable Gate Array (FPGA) is proposed for increasing the security and the frame rate of a smart IoT edge device. The private encryption key is safely embedded in the FPGA unit to avoid being exposed in the Random Access Memory (RAM). This allows the edge device to securely store and authenticate the key, protecting the data transmitted from the same Integrated Circuit (IC). Additionally, the edge device can simultaneously publish and route a camera stream using a lightweight communication protocol, achieving a frame rate of 14 frames per Second (fps). The performance of the MPSoC is compared to a NVIDIA Jetson Nano (NJN) and a Raspberry Pi 4 (RPI4) and it is found that the RPI4 is the most cost-effective solution but with lower frame rate, the NJN is the fastest because it can achieve higher frame-rate but it is not secure, and the MPSoC is the optimal solution because it offers a balanced frame rate and it is secure because it never exposes the secure key into the memory. Results: The proposed system successfully addresses the challenges of security, scalability, and efficiency in an IoT system used for surveillance and monitoring. The encryption key is securely stored and authenticated, and the edge device is able to simultaneously publish and route a camera stream feed high-definition images at 14 fps.


title: Early Detection for Multiversion Concurrency Control Conflicts in Hyperledger Fabric

abstract: Hyperledger Fabric is a popular permissioned blockchain system that features a highly modular and extensible system for deploying permissioned blockchains which are expected to have a major effect on a wide range of sectors. Unlike traditional blockchain systems such as Bitcoin and Ethereum, Hyperledger Fabric uses the EOV model for transaction processing: the submitted transactions are executed by the endorsing peer, ordered and batched by the ordering services, and validated by the validating peers. Due to this EOV workflow, a well-documented issue that arises is the multi-version concurrency control conflict. This happens when two transactions try to writes and read the same key in the ledger at the same time. Existing solutions to address this problem includes eliminating blocks in favor of streaming transactions, repairing conflicts during the ordering phase, and automatically merging the conflicting transactions using CRDT (Conflict Free Replicated Data) techniques. In this paper, we propose a novel solution called Early Detection for MVCC Conflicts. Our solution detects the conflicting transactions at an early stage of the transaction execution instead of processing them until the validation phase to be aborted. The advantage of our solution is that it detects conflict as soon as possible to minimize the overhead of conflicting transaction on the network resulting in the reduction of the end-to-end transaction latency and the increase of the system's effective throughput. We have successfully implemented our solution in Hyperledger Fabric. We propose three different implementations which realize early detection. Our results show that our solutions all perform better than the baseline Fabric, with our best solution SyncMap which improves the goodput by up to 23% and reduces the latency by up to 80%.


title: Bayesian Models of Functional Connectomics and Behavior

abstract: The problem of jointly analysing functional connectomics and behavioral data is extremely challenging owing to the complex interactions between the two domains. In addition, clinical rs-fMRI studies often have to contend with limited samples, especially in the case of rare disorders. This data-starved regimen can severely restrict the reliability of classical machine learning or deep learning designed to predict behavior from connectivity data. In this work, we approach this problem from the lens of representation learning and bayesian modeling. To model the distributional characteristics of the domains, we first examine the ability of approaches such as Bayesian Linear Regression, Stochastic Search Variable Selection after performing a classical covariance decomposition. Finally, we present a fully bayesian formulation for joint representation learning and prediction. We present preliminary results on a subset of a publicly available clinical rs-fMRI study on patients with Autism Spectrum Disorder.


title: On abstract results of operator representation of frames in Hilbert spaces

abstract: In this paper, we give a multiplication operator representation of bounded self-adjoint operators T on a Hilbert space H such that -- is a frame for H, for some -- . We state a necessary condition in order for a frame -- to have a representation of the form -- . Frame sequence -- with the synthesis operator U, have also been characterized in terms of the behavior spectrum of -- . Next, using the operator response of an element with respect to a unit vector in H, frames -- of the form -- are characterized. We also consider stability frames as -- . Finally, we conclude this note by raising a conjecture connecting frame theory and operator theory


title: LitAR: Visually Coherent Lighting for Mobile Augmented Reality

abstract: An accurate understanding of omnidirectional environment lighting is crucial for high-quality virtual object rendering in mobile augmented reality (AR). In particular, to support reflective rendering, existing methods have leveraged deep learning models to estimate or have used physical light probes to capture physical lighting, typically represented in the form of an environment map. However, these methods often fail to provide visually coherent details or require additional setups. For example, the commercial framework ARKit uses a convolutional neural network that can generate realistic environment maps; however the corresponding reflective rendering might not match the physical environments. In this work, we present the design and implementation of a lighting reconstruction framework called LitAR that enables realistic and visually-coherent rendering. LitAR addresses several challenges of supporting lighting information for mobile AR. First, to address the spatial variance problem, LitAR uses two-field lighting reconstruction to divide the lighting reconstruction task into the spatial variance-aware near-field reconstruction and the directional-aware far-field reconstruction. The corresponding environment map allows reflective rendering with correct color tones. Second, LitAR uses two noise-tolerant data capturing policies to ensure data quality, namely guided bootstrapped movement and motion-based automatic capturing. Third, to handle the mismatch between the mobile computation capability and the high computation requirement of lighting reconstruction, LitAR employs two novel real-time environment map rendering techniques called multi-resolution projection and anchor extrapolation. These two techniques effectively remove the need of time-consuming mesh reconstruction while maintaining visual quality.


title: Inhomogeneous superconductivity onset in FeSe studied by transport properties

abstract: Heterogeneous superconductivity onset is a common phenomenon in high-$T_c$ superconductors of both the cuprate and iron-based families. It is manifested by a fairly wide transition from the metallic to zero-resistance state. Usually, in these strongly anisotropic materials, superconductivity (SC) first appears as isolated domains. This leads to anisotropic excess conductivity above $T_c$, and the transport measurements provide valuable information about the SC domain structure deep within the sample. In bulk samples, this anisotropic SC onset gives an approximate average shape of SC grains, while in thin samples it also indicates the average size of SC grains. In this work, both interlayer and intralayer resistivity are measured as a function of temperature in FeSe samples of various thickness. To measure the interlayer resistivity, FeSe mesa structures oriented across the layers were fabricated using FIB. As the sample thickness decreases, a significant increase in superconducting transition temperature $T_c$ is observed : $T_c$ raises from 8 K in bulk material to 12 K in microbridges of thickness $\sim 40$ nm. We apply analytical and numerical calculations to analyze these and earlier data and find the aspect ratio and size of the SC domains in FeSe consistent with our resistivity and diamagnetic response measurements. We propose a simple and fairly accurate method for estimating the aspect ratio of SC domains from $T_c$ anisotropy in samples of various small thickness. The relationship between nematic and superconducting domains in FeSe is discussed. We also generalize the analytical formulas for conductivity in heterogeneous anisotropic superconductors to the case of elongated SC domains of two perpendicular orientations with equal volume fractions, corresponding to the nematic domain structure in various Fe-based superconductors.


title: A Multiwavelength Study of Active Galactic Nuclei in Post-Merger Remnants

abstract: We investigate the role of galaxy mergers in triggering AGN in the nearby Universe. Our analysis is based on a sample of 79 post-merger remnant galaxies with deep X-ray observations from Chandra/XMM-Newton capable of detecting a low-luminosity AGN of > 10^40.5 erg s^-1. This sample is derived from a visually classified, volume-limited sample of 807 post-mergers identified in the Sloan Digital Sky Survey Data Release 14 with log M*/M_sun > 10.5 and 0.02 < z < 0.06. We find the X-ray AGN fraction in this sample is 55.7% +\- 5.6% compared to 23.6% +\- 2.8% for a mass and redshift matched non-interacting control sample. The multi-wavelength AGN fraction (identified as an AGN in one of X-ray, IR, radio or optical diagnostics) for post-mergers is 76.6% +\- 4.8% compared to 39.1% +\- 3.2% for controls. Thus post-mergers exhibit a high overall AGN fraction with an excess between 2 - 4 depending on the AGN diagnostics used. In addition, we find most optical, IR, and radio AGN are also identified as X-ray AGN while a large fraction of X-ray AGN are not identified in any other diagnostic. This highlights the importance of deep X-ray imaging to identify AGN. We find the X-ray AGN fraction of post-mergers is independent of stellar mass above log M*/M_sun > 10.5 unlike the trend seen in control galaxies. Overall, our results show that post-merger galaxies are a good tracer of the merger-AGN connection and strongly support the theoretical expectations that mergers trigger AGN.


title: CNN-Based Action Recognition and Pose Estimation for Classifying Animal Behavior from Videos: A Survey

abstract: Classifying the behavior of humans or animals from videos is important in biomedical fields for understanding brain function and response to stimuli. Action recognition, classifying activities performed by one or more subjects in a trimmed video, forms the basis of many of these techniques. Deep learning models for human action recognition have progressed significantly over the last decade. Recently, there is an increased interest in research that incorporates deep learning-based action recognition for animal behavior classification. However, human action recognition methods are more developed. This survey presents an overview of human action recognition and pose estimation methods that are based on convolutional neural network (CNN) architectures and have been adapted for animal behavior classification in neuroscience. Pose estimation, estimating joint positions from an image frame, is included because it is often applied before classifying animal behavior. First, we provide foundational information on algorithms that learn spatiotemporal features through 2D, two-stream, and 3D CNNs. We explore motivating factors that determine optimizers, loss functions and training procedures, and compare their performance on benchmark datasets. Next, we review animal behavior frameworks that use or build upon these methods, organized by the level of supervision they require. Our discussion is uniquely focused on the technical evolution of the underlying CNN models and their architectural adaptations (which we illustrate), rather than their usability in a neuroscience lab. We conclude by discussing open research problems, and possible research directions. Our survey is designed to be a resource for researchers developing fully unsupervised animal behavior classification systems of which there are only a few examples in the literature.


title: Kinked tracks from $\Sigma^+$ baryons as a probe of light quark polarizations

abstract: Jets initiated by energetic quarks carry in an obvious way approximate information about the original quark energy and momentum. In a less obvious way, they also carry information about the quark polarization. Previous works have shown how the polarization information can be extracted by ATLAS and CMS using events in which the quark hadronizes to an energetic baryon - $\Lambda_b$ for $b$ quarks, $\Lambda_c$ for $c$ quarks and $\Lambda$ for $s$ quarks. In this paper we extend these proposals to $\Sigma^+$ baryons, which can provide sensitivity to polarizations of $s$ and $u$ quarks. We analyze the various aspects of how the kinked track signature of the $\Sigma^+ \to p \pi^0$ decay can be used for this purpose. We evaluate the feasibility of such measurements in $t\bar t$ samples.


title: Cosmological simulations of the same spiral galaxy: connecting the dark matter distribution of the host halo with the subgrid baryonic physics

abstract: The role of baryonic physics, star formation, and stellar feedback, in shaping the galaxies and their host halos is an evolving topic. The dark matter aspects are illustrated in this work by showing distribution features in a Milky-Way-sized halo. We focus on the halo morphology, geometry, and profile as well as the phase space distribution using one dark matter only and five hydrodynamical cosmological high-resolution simulations of the same halo with different subgrid prescriptions for the baryonic physics (Kennicut versus multi-freefall star formation and delayed cooling versus mechanical supernovae feedback). If some general properties like the relative halo-galaxy orientation are similar, the modifications of the gravitational potential due to the presence of baryons are found to induce different dark matter distributions (rounder and more concentrated halo). The mass density profile as well as the velocity distribution are modified distinctively according to the specific resulting baryonic distribution highlighting the variability of those properties (e.g inner power index from 1.3 to 1.8, broader speed distribution). The uncertainties on those features are of paramount importance for dark matter phenomenology, particularly when dealing with dark matter dynamics or direct and indirect detection searches. As a consequence, dark matter properties and prospects using cosmological simulations require improvement on baryonic physics description. Modeling such processes is a key issue not only for galaxy formation but also for dark matter investigations.


title: BuildSeg: A General Framework for the Segmentation of Buildings

abstract: Building segmentation from aerial images and 3D laser scanning (LiDAR) is a challenging task due to the diversity of backgrounds, building textures, and image quality. While current research using different types of convolutional and transformer networks has considerably improved the performance on this task, even more accurate segmentation methods for buildings are desirable for applications such as automatic mapping. In this study, we propose a general framework termed \emph{BuildSeg} employing a generic approach that can be quickly applied to segment buildings. Different data sources were combined to increase generalization performance. The approach yields good results for different data sources as shown by experiments on high-resolution multi-spectral and LiDAR imagery of cities in Norway, Denmark and France. We applied ConvNeXt and SegFormer based models on the high resolution aerial image dataset from the MapAI-competition. The methods achieved an IOU of 0.7902 and a boundary IOU of 0.6185. We used post-processing to account for the rectangular shape of the objects. This increased the boundary IOU from 0.6185 to 0.6189.


title: Elizabethan vortices

abstract: Radial solutions to the elliptic sinh-Gordon and Tzitzeica equations can be interpreted as Abelian vortices on certain surfaces of revolution. These surfaces have a conical excess angle at infinity (in a way which makes them similar to Elizabethan ruff collars). While they can not be embedded in the Euclidean 3-space, we will show that they can be globally embedded in the hyperbolic space. The existence of these hyperbolic embeddings follows from the asymptotic analysis of a Painleve III ODE.


title: Do We Have Any Viable Solution to the Measurement Problem?

abstract: Wallace (2022) has recently argued that a number of popular approaches to the measurement problem can't be fully extended to relativistic quantum mechanics and quantum field theory; Wallace thus contends that as things currently stand, only the unitary-only approaches to the measurement problem are viable. However, the unitary-only approaches face serious epistemic problems which may threaten their viability as solutions, and thus we consider that it remains an urgent outstanding problem to find a viable solution to the measurement problem which can be extended to relativistic quantum mechanics. In this article we seek to understand in general terms what such a thing might look like. We argue that in order to avoid serious epistemic problems, the solution must be a single-world realist approach, and we further argue that any single-world realist approach which is able to reproduce the predictions of relativistic quantum mechanics will most likely have the property that our observable reality does not supervene on dynamical, precisely-defined microscopic beables. Thus we suggest three possible routes for further exploration: observable reality could be approximate and emergent, as in relational quantum mechanics with the addition of cross-perspective links, or observable reality could supervene on beables which are not microscopically defined, as in the consistent histories approach, or observable reality could supervene on beables which are not dynamical, as in Kent's solution to the Lorentzian classical reality problem. We conclude that once all of these issues are taken into account, the options for a viable solution to the measurement problem are significantly narrowed down.


title: RedBit: An End-to-End Flexible Framework for Evaluating the Accuracy of Quantized CNNs

abstract: In recent years, Convolutional Neural Networks (CNNs) have become the standard class of deep neural network for image processing, classification and segmentation tasks. However, the large strides in accuracy obtained by CNNs have been derived from increasing the complexity of network topologies, which incurs sizeable performance and energy penalties in the training and inference of CNNs. Many recent works have validated the effectiveness of parameter quantization, which consists in reducing the bit width of the network's parameters, to enable the attainment of considerable performance and energy efficiency gains without significantly compromising accuracy. However, it is difficult to compare the relative effectiveness of different quantization methods. To address this problem, we introduce RedBit, an open-source framework that provides a transparent, extensible and easy-to-use interface to evaluate the effectiveness of different algorithms and parameter configurations on network accuracy. We use RedBit to perform a comprehensive survey of five state-of-the-art quantization methods applied to the MNIST, CIFAR-10 and ImageNet datasets. We evaluate a total of 2300 individual bit width combinations, independently tuning the width of the network's weight and input activation parameters, from 32 bits down to 1 bit (e.g., 8/8, 2/2, 1/32, 1/1, for weights/activations). Upwards of 20000 hours of computing time in a pool of state-of-the-art GPUs were used to generate all the results in this paper. For 1-bit quantization, the accuracy losses for the MNIST, CIFAR-10 and ImageNet datasets range between [0.26%, 0.79%], [9.74%, 32.96%] and [10.86%, 47.36%] top-1, respectively. We actively encourage the reader to download the source code and experiment with RedBit, and to submit their own observed results to our public repository, available at https://github.com/IT-Coimbra/RedBit.


title: Geometric Graph Learning with Extended Atom-Types Features for Protein-Ligand Binding Affinity Prediction

abstract: Understanding and accurately predicting protein-ligand binding affinity are essential in the drug design and discovery process. At present, machine learning-based methodologies are gaining popularity as a means of predicting binding affinity due to their efficiency and accuracy, as well as the increasing availability of structural and binding affinity data for protein-ligand complexes. In biomolecular studies, graph theory has been widely applied since graphs can be used to model molecules or molecular complexes in a natural manner. In the present work, we upgrade the graph-based learners for the study of protein-ligand interactions by integrating extensive atom types such as SYBYL and extended connectivity interactive features (ECIF) into multiscale weighted colored graphs (MWCG). By pairing with the gradient boosting decision tree (GBDT) machine learning algorithm, our approach results in two different methods, namely $^\text{sybyl}\text{GGL}$-Score and $^\text{ecif}\text{GGL}$-Score. Both of our models are extensively validated in their scoring power using three commonly used benchmark datasets in the drug design area, namely CASF-2007, CASF-2013, and CASF-2016. The performance of our best model $^\text{sybyl}\text{GGL}$-Score is compared with other state-of-the-art models in the binding affinity prediction for each benchmark. While both of our models achieve state-of-the-art results, the SYBYL atom-type model $^\text{sybyl}\text{GGL}$-Score outperforms other methods by a wide margin in all benchmarks.


title: Calibrated Data-Dependent Constraints with Exact Satisfaction Guarantees

abstract: We consider the task of training machine learning models with data-dependent constraints. Such constraints often arise as empirical versions of expected value constraints that enforce fairness or stability goals. We reformulate data-dependent constraints so that they are calibrated: enforcing the reformulated constraints guarantees that their expected value counterparts are satisfied with a user-prescribed probability. The resulting optimization problem is amendable to standard stochastic optimization algorithms, and we demonstrate the efficacy of our method on a fairness-sensitive classification task where we wish to guarantee the classifier's fairness (at test time).


title: Young Male and Female Scientists: A Quantitative Exploratory Study of the Changing Demographics of the Global Scientific Workforce

abstract: In this study, the global scientific workforce is explored through large-scale, generational, cross-sectional, and longitudinal approaches. We examine 4.3 million nonoccasional scientists from 38 OECD countries publishing in 1990-2021. Our interest is in the changing distribution of young male and female scientists over time across 16 STEMM (science, technology, engineering, mathematics, medicine) disciplines. We unpack the details of the changing scientific workforce using age groups. Some disciplines are already numerically dominated by women, and the change is fast in some and slow in other disciplines. In one-third of disciplines, there are already more youngest female than male scientists. Across all disciplines combined, the majority of women are young women. And more than half of women scientists (55.02%) are located in medicine. The usefulness of global bibliometric data sources in analyzing the scientific workforce along gender, age, discipline, and time is tested. Traditional aggregated data about scientists in general hide a nuanced picture of the changing gender dynamics within and across disciplines and age groups. The limitations of bibliometric datasets are explored, and global studies are compared with national-level studies. The methodological choices and their implications are shown, and new opportunities for how to study scientists globally are discussed.


title: Who Should Predict? Exact Algorithms For Learning to Defer to Humans

abstract: Automated AI classifiers should be able to defer the prediction to a human decision maker to ensure more accurate predictions. In this work, we jointly train a classifier with a rejector, which decides on each data point whether the classifier or the human should predict. We show that prior approaches can fail to find a human-AI system with low misclassification error even when there exists a linear classifier and rejector that have zero error (the realizable setting). We prove that obtaining a linear pair with low error is NP-hard even when the problem is realizable. To complement this negative result, we give a mixed-integer-linear-programming (MILP) formulation that can optimally solve the problem in the linear setting. However, the MILP only scales to moderately-sized problems. Therefore, we provide a novel surrogate loss function that is realizable-consistent and performs well empirically. We test our approaches on a comprehensive set of datasets and compare to a wide range of baselines.


title: Generalized Neural Closure Models with Interpretability

abstract: Improving the predictive capability and computational cost of dynamical models is often at the heart of augmenting computational physics with machine learning (ML). However, most learning results are limited in interpretability and generalization over different computational grid resolutions, initial and boundary conditions, domain geometries, and physical or problem-specific parameters. In the present study, we simultaneously address all these challenges by developing the novel and versatile methodology of unified neural partial delay differential equations. We augment existing/low-fidelity dynamical models directly in their partial differential equation (PDE) forms with both Markovian and non-Markovian neural network (NN) closure parameterizations. The melding of the existing models with NNs in the continuous spatiotemporal space followed by numerical discretization automatically allows for the desired generalizability. The Markovian term is designed to enable extraction of its analytical form and thus provides interpretability. The non-Markovian terms allow accounting for inherently missing time delays needed to represent the real world. We obtain adjoint PDEs in the continuous form, thus enabling direct implementation across differentiable and non-differentiable computational physics codes, different ML frameworks, and treatment of nonuniformly-spaced spatiotemporal training data. We demonstrate the new generalized neural closure models (gnCMs) framework using four sets of experiments based on advecting nonlinear waves, shocks, and ocean acidification models. Our learned gnCMs discover missing physics, find leading numerical error terms, discriminate among candidate functional forms in an interpretable fashion, achieve generalization, and compensate for the lack of complexity in simpler models. Finally, we analyze the computational advantages of our new framework.

categories: cs.LG math.DS physics.flu-dyn physics.geo-ph


title: Doubly Robust Counterfactual Classification

abstract: We study counterfactual classification as a new tool for decision-making under hypothetical (contrary to fact) scenarios. We propose a doubly-robust nonparametric estimator for a general counterfactual classifier, where we can incorporate flexible constraints by casting the classification problem as a nonlinear mathematical program involving counterfactuals. We go on to analyze the rates of convergence of the estimator and provide a closed-form expression for its asymptotic distribution. Our analysis shows that the proposed estimator is robust against nuisance model misspecification, and can attain fast $\sqrt{n}$ rates with tractable inference even when using nonparametric machine learning approaches. We study the empirical performance of our methods by simulation and apply them for recidivism risk prediction.


title: Efficiently Computing Sparse Fourier Transforms of $q$-ary Functions

abstract: Fourier transformations of pseudo-Boolean functions are popular tools for analyzing functions of binary sequences. Real-world functions often have structures that manifest in a sparse Fourier transform, and previous works have shown that under the assumption of sparsity the transform can be computed efficiently. But what if we want to compute the Fourier transform of functions defined over a $q$-ary alphabet? These types of functions arise naturally in many areas including biology. A typical workaround is to encode the $q$-ary sequence in binary, however, this approach is computationally inefficient and fundamentally incompatible with the existing sparse Fourier transform techniques. Herein, we develop a sparse Fourier transform algorithm specifically for $q$-ary functions of length $n$ sequences, dubbed $q$-SFT, which provably computes an $S$-sparse transform with vanishing error as $q^n \rightarrow \infty$ in $O(Sn)$ function evaluations and $O(S n^2 \log q)$ computations, where $S = q^{n\delta}$ for some $\delta < 1$. Under certain assumptions, we show that for fixed $q$, a robust version of $q$-SFT has a sample complexity of $O(Sn^2)$ and a computational complexity of $O(Sn^3)$ with the same asymptotic guarantees. We present numerical simulations on synthetic and real-world RNA data, demonstrating the scalability of $q$-SFT to massively high dimensional $q$-ary functions.


title: Risk-aware Vehicle Motion Planning Using Bayesian LSTM-Based Model Predictive Control

abstract: Understanding the probabilistic traffic environment is a vital challenge for the motion planning of autonomous vehicles. To make feasible control decisions, forecasting future trajectories of adjacent cars is essential for intelligent vehicles to assess potential conflicts and react to reduce the risk. This paper first introduces a Bayesian Long Short-term Memory (BLSTM) model to learn human drivers' behaviors and habits from their historical trajectory data. The model predicts the probability distribution of surrounding vehicles' positions, which are used to estimate dynamic conflict risks. Next, a hybrid automaton is built to model the basic motions of a car, and the conflict risks are assessed for real-time state-space transitions based on environmental information. Finally, a BLSTM-based Model Predictive Control (MPC) is built to navigate vehicles through safe paths with the least predicted conflict risk. By merging BLSTM with MPC, the designed neural-based MPC overcomes the defect that traditional MPC is hard to model uncertain conflict risks. The simulation results show that our proposed BLSTM-based MPC performs better than human drivers because it can foresee potential conflicts and take action to avoid them.


title: Shades of Iteration: from Elgot to Kleene

abstract: Notions of iteration range from the arguably most general Elgot iteration to a very specific Kleene iteration. The fundamental nature of Elgot iteration has been extensively explored by Bloom and Esik in the form of iteration theories, while Kleene iteration became extremely popular as an integral part of (untyped) formalisms, such as automata theory, regular expressions and Kleene algebra. Here, we establish a formal connection between Elgot iteration and Kleene iteration in the form of Elgot monads and Kleene monads, respectively. We also introduce a novel class of while-monads, which like Kleene monads admit a relatively simple description in algebraic terms. Like Elgot monads, while-monads cover a large variety of models that meaningfully support while-loops, but may fail the Kleene algebra laws, or even fail to support a Kleen iteration operator altogether.


title: Symplectic Perturbation Theory in Massive Ambitwistor Space: A Zig-Zag Theory of Massive Spinning Particles

abstract: We develop a theory of massive spinning particles interacting with background fields in four spacetime dimensions in which holomorphy and chirality play a central role. Applying a perturbation theory of symplectic forms to the massive twistor space as a K\"ahler manifold, we find that the spin precession behavior of a massive spinning particle is directly determined from the manner in which self-dual and anti-self-dual field strengths permeate into "complex spacetime." Especially, the particle shows the minimally coupled precession behavior if self-dual field strength continues holomorphically into the complex: the Newman-Janis shift. In general, computing the momentum impulse shows that the parameters that control generic non-holomorphic continuations are directly related to the coupling constants in the massive-massive-massless spinning on-shell amplitude of Arkani-Hamed, Huang, and Huang, and thus they are interpreted as the single-curvature Wilson coefficients given by Levi and Steinhoff, redefined on complex worldlines. Finally, exact expressions for Kerr and $\sqrt{\text{Kerr}}$ actions are bootstrapped in monochromatic self-dual plane-wave backgrounds from symplectivity and a matching between classical scattering and the on-shell amplitude, from which we obtain all-order exact impulses of classical observables.


title: First-principles study of hydrogen dynamics in monoclinic TiO

abstract: The existence of intrinsic vacancies in cubic (monoclinic) TiO suggests opportunity for hydrogen absorption, which was addressed in recent experiments. In the present work, based on first principle calculations, the preferences are studied for the hydrogen absorption sites and diffusion paths between them. The oxygen vacancies are found to be primary hydrogen traps with absorption energy of -2.87 eV. The plausible channels for hydrogen diffusion between adjacent vacancy sites (ordered in the monoclinic TiO structure) are compared with the help of calculations done with the nudge elastic band method. Several competitive channels are identified, with barrier heights varying from 2.87 to 3.71 eV, that is high enough to ensure relative stability of trapped hydrogen atoms at oxygen vacancy sites. Moreover, the possibility of adsorption of molecular hydrogen was tested and found improbable, in the sense that the H2 molecules penetrating the TiO crystal are easily dissociated (and released atoms tend to proceed towards oxygen vacancy sites). These results suggest that hydrogen may persist in oxygen vacancy sites up to high enough temperatures.


title: Remarks on the oscillation energy of Legendrian isotopies

abstract: We construct non-compact contact manifolds containing compact Legendrians which can be displaced from their Reeb flow with arbitrarily small oscillation energy. We use this to show the Shelukhin-Chekanov-Hofer pseudo-metric considered by Rosen and Zhang is degenerate on the isotopy class of the constructed Legendrians. Other aspects related to oscillation energy of Legendrian isotopies are explored.


title: Efficiency in Collective Decision-Making via Quadratic Transfers

abstract: Consider the following collective choice problem: a group of budget constrained agents must choose one of several alternatives. Is there a budget balanced mechanism that: i) does not depend on the specific characteristics of the group, ii) does not require unaffordable transfers, and iii) implements utilitarianism if the agents' preferences are quasilinear and their private information? We study the following procedure: every agent can express any intensity of support or opposition to each alternative, by transferring to the rest of the agents wealth equal to the square of the intensity expressed; and the outcome is determined by the sums of the expressed intensities. We prove that as the group grows large, in every equilibrium of this quadratic-transfers mechanism, each agent's transfer converges to zero, and the probability that the efficient outcome is chosen converges to one.


title: The Binary Linearization Complexity of Pseudo-Boolean Functions

abstract: We consider the problem of linearizing a pseudo-Boolean function $f : \{0,1\}^n \to \mathbb{R}$ by means of $k$ Boolean functions. Such a linearization yields an integer linear programming formulation with only $k$ auxiliary variables. This motivates the definition of the linarization complexity of $f$ as the minimum such $k$. Our theoretical contributions are the proof that random polynomials almost surely have a high linearization complexity and characterizations of its value in case we do or do not restrict the set of admissible Boolean functions. The practical relevance is shown by devising and evaluating integer linear programming models of two such linearizations for the low auto-correlation binary sequences problem. Still, many problems around this new concept remain open.


title: On thermodynamically consistent quasiparticle model at finite chemical potential

abstract: We explore the quasiparticle model at finite chemical potential related to Ru-Keng Su's distinguished contributions to the topic. Besides, we discuss recent developments in the model, and in particular, one argues that the effective mass of the quasiparticle might attain a specific form as a function of momentum, in addition to its dependence on temperature and chemical potential. Unlike the approaches based on the properties of underlying symmetry or renormalization group, the momentum dependence emerges as a special solution to an integro-differential equation resulting from the underlying thermodynamic consistency. Moreover, this special solution to the problem is shown to be more general than previously explored in the literature. Instead of fitting to the lattice QCD data at vanishing chemical potential, in this work, we adopt a ``bottom-up'' approach by assuming some analytic ansatzes that are manifestly thermodynamically consistent. The remaining physical quantities are subsequently derived, and possible implications are also addressed.


title: Efficient Loop Conditions for Bounded Model Checking Hyperproperties

abstract: Bounded model checking (BMC) is an effective technique for hunting bugs by incrementally exploring the state space of a system. To reason about infinite traces through a finite structure and to ultimately obtain completeness, BMC incorporates loop conditions that revisit previously observed states. This paper focuses on developing loop conditions for BMC of HyperLTL- a temporal logic for hyperproperties that allows expressing important policies for security and consistency in concurrent systems, etc. Loop conditions for HyperLTL are more complicated than for LTL, as different traces may loop inconsistently in unrelated moments. Existing BMC approaches for HyperLTL only considered linear unrollings without any looping capability, which precludes both finding small infinite traces and obtaining a complete technique. We investigate loop conditions for HyperLTL BMC, where the HyperLTL formula can contain up to one quantifier alternation. We first present a general complete automata-based technique which is based on bounds of maximum unrollings. Then, we introduce alternative simulation-based algorithms that allow exploiting short loops effectively, generating SAT queries whose satisfiability guarantees the outcome of the original model checking problem. We also report empirical evaluation of the prototype implementation of our BMC techniques using Z3py.


title: V-Guard: An Efficient Permissioned Blockchain for Achieving Consensus under Dynamic Memberships in V2X

abstract: This paper presents V-Guard, a new permissioned blockchain that achieves consensus for vehicular data under changing memberships, targeting the problem in V2X networks where vehicles are often intermittently connected on the roads. To achieve this goal, V-Guard integrates membership management into the consensus process for agreeing on data entries. It binds a data entry with a membership configuration profile that describes responsible vehicles for achieving consensus for the data entry. As such, V-Guard produces chained consensus results of both data entries and their residing membership profiles, which enables consensus to be achieved seamlessly under changing memberships. In addition, V-Guard separates the ordering of transactions from consensus, allowing concurrent ordering instances and periodic consensus instances to order and commit data entries. These features make V-Guard efficient for achieving consensus under dynamic memberships with high throughput and latency performance.


title: What artificial intelligence might teach us about the origin of human language

abstract: This study explores an interesting pattern emerging from research that combines artificial intelligence with sound symbolism. In these studies, supervised machine learning algorithms are trained to classify samples based on the sounds of referent names. Machine learning algorithms are efficient learners of sound symbolism, but they tend to bias one category over the other. The pattern is this: when a category arguably represents greater threat, the algorithms tend to overpredict to that category. A hypothesis, framed by error management theory, is presented that proposes that this may be evidence of an adaptation to preference cautious behaviour. This hypothesis is tested by constructing extreme gradient boosted (XGBoost) models using the sounds that make up the names of Chinese, Japanese and Korean Pokemon and observing classification error distribution.


title: Lower bounds for multicolor van der Waerden numbers

abstract: We give an exponential improvement to the diagonal van der Waerden numbers for $r\ge 5$ colors.


title: Robust and Sparse M-Estimation of DOA

abstract: A robust and sparse Direction of Arrival (DOA) estimator is derived for array data that follows a Complex Elliptically Symmetric (CES) distribution with zero-mean and finite second-order moments. The derivation allows to choose the loss function and four loss functions are discussed in detail: the Gauss loss which is the Maximum-Likelihood (ML) loss for the circularly symmetric complex Gaussian distribution, the ML-loss for the complex multivariate $t$-distribution (MVT) with $\nu$ degrees of freedom, as well as Huber and Tyler loss functions. For Gauss loss, the method reduces to Sparse Bayesian Learning (SBL). The root mean square DOA error of the derived estimators is discussed for Gaussian, MVT, and $\epsilon$-contaminated data. The robust SBL estimators perform well for all cases and nearly identical with classical SBL for Gaussian noise.


title: Quantum-jump analysis of frequency up-conversion amplification without inversion in a four-level scheme

abstract: In this work, we use the quantum-jump approach to study light-matter interactions in a four-level scheme giving rise to frequency up-conversion amplification without inversion (AWI). The results obtained apply to the case of neutral Hg vapour where, recently, it has been reported AWI of a probe field in the UV regime, opening the way for lasing without inversion in this range of frequencies. We show that, in this scheme, the key element to obtain maximum amplification is the fulfillment of the three-photon resonance condition, which favors coherent evolution periods associated with the probe field gain. We also investigate the parameter values in order to optimize amplification. The present study extends the theoretical understanding of the underlying mechanisms involved in AWI.


title: Automated Software Testing Starting from Static Analysis: Current State of the Art

abstract: The mass production of complex software has made it impossible to manually test it for security vulnerabilities. Automated security testing tools come in a variety of flavors, function at various stages of software development, and target different categories of software vulnerabilities. It is great that we have a plethora of automated tools to choose from, but it is a problem that their adoption and recognition are not prominent. The purpose of this study is to explore the possibilities of existing techniques while also broadening the horizon by exploring the future of security testing tools and related techniques.


title: Modelling human logical reasoning process in dynamic environmental stress with cognitive agents

abstract: Modelling human cognition can provide key insights into behavioral dynamics under changing conditions. This enables synthetic data generation and guides adaptive interventions for cognitive regulation. Challenges arise when environments are highly dynamic, obscuring stimulus-behavior relationships. We propose a cognitive agent integrating drift-diffusion with deep reinforcement learning to simulate granular stress effects on logical reasoning process. Leveraging a large dataset of 21,157 logical responses, we investigate performance impacts of dynamic stress. This prior knowledge informed model design and evaluation. Quantitatively, the framework improves cognition modelling by capturing both subject-specific and stimuli-specific behavioural differences. Qualitatively, it captures general trends in human logical reasoning under stress. Our approach is extensible to examining diverse environmental influences on cognition and behavior. Overall, this work demonstrates a powerful, data-driven methodology to simulate and understand the vagaries of human logical reasoning process in dynamic contexts.


title: Analogy between Boltzmann machines and Feynman path integrals

abstract: We provide a detailed exposition of the connections between Boltzmann machines commonly utilized in machine learning problems and the ideas already well known in quantum statistical mechanics through Feynman's description of the same. We find that this equivalence allows the interpretation that the hidden layers in Boltzmann machines and other neural network formalisms are in fact discrete versions of path elements that are present within the Feynman path-integral formalism. Since Feynman paths are the natural and elegant depiction of interference phenomena germane to quantum mechanics, it appears that in machine learning, the goal is to find an appropriate combination of ``paths'', along with accumulated path-weights, through a network that cumulatively capture the correct $x \rightarrow y$ map for a given mathematical problem. As a direct consequence of this analysis, we are able to provide general quantum circuit models that are applicable to both Boltzmann machines and to Feynman path integral descriptions. Connections are also made to inverse quantum scattering problems which allow a robust way to define ``interpretable'' hidden layers.


title: Fixed points of the sum of divisors function on $F_2[x]$

abstract: We work an analogue of a classical arithmetic problem over polynomials. More precisely, we study the fixed points $F$ of the sum of divisors function $\sigma : F_2[x] \mapsto F_2[x]$ (defined \emph{mutatis mutandi} like the usual sum of divisors over the integers) of the form $F := A^2 \cdot S$, $S$ square-free, with $\omega(S) \leq 3$, coprime with $A$, for $A$ even, of whatever degree, under some conditions. This gives a characterization of $5$ of the $11$ known fixed points of $\sigma$ in $F_2[x]$


title: Applying causal inference to inform early-childhood policy from administrative data

abstract: Improving public policy is one of the key roles of governments, and they can do this in an evidence-based way using administrative data. Causal inference for observational data improves on current practice of using descriptive or predictive analyses to inform policy decisions. Causal inference allows analysts to estimate the impact a policy change would have on the population if the encoded assumptions about the data generation process are valid. In this paper, we discuss the importance of causal analysis methods when analysing data to inform policy decisions. We take the education sector as a case study and provide examples of when to use a causal analysis. We use simulation to demonstrate the vital role causal diagrams play in variable selection and how bias can be introduced if extraneous variables are included in the model. Our exploration provides clear evidence for the utility of causal methods and practical examples of how to conduct such analyses. The paper promotes the incorporation of these methods in policy both for improved educational outcomes and scientific understanding.


title: New bounds on the cardinality of Hausdorff spaces and regular spaces

abstract: Using weaker versions of the cardinal function $\psi_c(X)$, we derive a series of new bounds for the cardinality of Hausdorff spaces and regular spaces that do not involve $\psi_c(X)$ nor its variants at all. For example, we show if $X$ is regular then $|X|\leq 2^{c(X)^{\pi\chi(X)}}$ and $|X|\leq 2^{c(X)\pi\chi(X)^{ot(X)}}$, where the cardinal function $ot(X)$, introduced by Tkachenko, has the property $ot(X)\leq\min\{t(X),c(X)\}$. It follows from the latter that a regular space with cellularity at most $\mathfrak{c}$ and countable $\pi$-character has cardinality at most $2^\mathfrak{c}$. For a Hausdorff space $X$ we show $|X|\leq 2^{d(X)^{\pi\chi(X)}}$, $|X|\leq d(X)^{\pi\chi(X)^{ot(X)}}$, and $|X|\leq 2^{\pi w(X)^{dot(X)}}$, where $dot(X)\leq\min\{ot(X),\pi\chi(X)\}$. None of these bounds involve $\psi_c(X)$ or $\psi(X)$. By introducing the cardinal functions $w\psi_c(X)$ and $d\psi_c(X)$ with the property $w\psi_c(X)d\psi_c(X)\leq\psi_c(X)$ for a Hausdorff space $X$, we show $|X|\leq\pi\chi(X)^{c(X)w\psi_c(X)}$ if $X$ is regular and $|X|\leq\pi\chi(X)^{c(X)d\psi_c(X)w\psi_c(X)}$ if $X$ is Hausdorff. This improves results of Sapirovskii and Sun. It is also shown that if $X$ is Hausdorff then $|X|\leq 2^{d(X)w\psi_c(X)}$, which appears to be new even in the case where $w\psi_c(X)$ is replaced with $\psi_c(X)$. Compact examples show that $\psi(X)$ cannot be replaced with $d\psi_c(X)w\psi_c(X)$ in the bound $2^{\psi(X)}$ for the cardinality of a compact Hausdorff space $X$. Likewise, $\psi(X)$ cannot be replaced with $d\psi_c(X)w\psi_c(X)$ in the Arhangel'skii-Sapirovskii bound $2^{L(X)t(X)\psi(X)}$ for the cardinality of a Hausdorff space $X$. Finally, we make several observations concerning homogeneous spaces in this connection.


title: First demonstration of neural sensing and control in a kilometer-scale gravitational wave observatory

abstract: Suspended optics in gravitational wave (GW) observatories are susceptible to alignment perturbations, particularly slow drifts over time, due to variations in temperature and seismic levels. Such misalignments affect the coupling of the incident laser beam into the optical cavities, degrade both circulating power and optomechanical photon squeezing and thus decrease the astrophysical sensitivity to merging binaries. Traditional alignment techniques involve differential wavefront sensing using multiple quadrant photodiodes but are often restricted in bandwidth and are limited by the sensing noise. We present the first-ever successful implementation of neural network-based sensing and control at a gravitational wave observatory and demonstrate low-frequency control of the signal recycling mirror at the GEO 600 detector. Alignment information for three critical optics is simultaneously extracted from the interferometric dark port camera images via a CNN-LSTM network architecture and is then used for MIMO control using soft actor-critic-based deep reinforcement learning. Overall sensitivity improvement achieved using our scheme demonstrates deep learning's capabilities as a viable tool for real-time sensing and control for current and next-generation GW interferometers.


title: Kinetic Insights into Bridge Cleavage Pathways in Periodic Mesoporous Organosilicas

abstract: Bridging functionalities in periodic mesoporous organosilicas (PMOs) enable new functionalities for a wide range of applications. Bridge cleavage is frequently observed during anneals required to form porous structures, yet the mechanism of these bridge cleavages has not been completely resolved. Here, we reveal these chemical transformations and their kinetic pathways on sub-millisecond timescales induced by laser heating. By varying anneal times and temperatures, the transformation dynamics of bridge cleavage and structural transformations, and their activation energies, are determined. The structural relaxation time for individual reactions and their effective local heating time are determined and compared, and results directly demonstrate the manipulation of different molecules through kinetic control of the sequence of reactions. By isolating and understanding the earliest stage of structural transformations, this study identifies the kinetic principles for new synthesis and post-processing routes to control individual molecules and reactions in PMOs and other material systems with multi-functionalities.


title: Learning-based Intelligent Surface Configuration, User Selection, Channel Allocation, and Modulation Adaptation for Jamming-resisting Multiuser OFDMA Systems

abstract: Reconfigurable intelligent surfaces (RISs) can potentially combat jamming attacks by diffusing jamming signals. This paper jointly optimizes user selection, channel allocation, modulation-coding, and RIS configuration in a multiuser OFDMA system under a jamming attack. This problem is non-trivial and has never been addressed, because of its mixed-integer programming nature and difficulties in acquiring channel state information (CSI) involving the RIS and jammer. We propose a new deep reinforcement learning (DRL)-based approach, which learns only through changes in the received data rates of the users to reject the jamming signals and maximize the sum rate of the system. The key idea is that we decouple the discrete selection of users, channels, and modulation-coding from the continuous RIS configuration, hence facilitating the RIS configuration with the latest twin delayed deep deterministic policy gradient (TD3) model. Another important aspect is that we show a winner-takes-all strategy is almost surely optimal for selecting the users, channels, and modulation-coding, given a learned RIS configuration. Simulations show that the new approach converges fast to fulfill the benefit of the RIS, due to its substantially small state and action spaces. Without the need of the CSI, the approach is promising and offers practical value.


title: Finding meaningful and workable applied mathematics problems in science

abstract: In this short review, I will summarize my research experience in three fields in applied mathematics: mathematical biology, applied probability, and applied discrete mathematics. Specifically, I will show how each project was initiated, and what wrong approaches were applied. Such details are important in learning how to do research, but they cannot be read out from research papers. I wish that students and junior researchers in applied mathematics could learn a lesson from this summary.


title: Non-phononic density of states of two-dimensional glasses revealed by random pinning

abstract: The vibrational density of states of glasses is considerably different from that of crystals. In particular, there exist spatially localized vibrational modes in glasses. The density of states of these non-phononic modes has been observed to follow $g(\omega) \propto \omega^4$, where $\omega$ is the frequency. However, in two-dimensional systems, the abundance of phonons makes it difficult to accurately determine this non-phononic density of states because they are strongly coupled to non-phononic modes and yield strong system-size and preparation-protocol dependencies. In this article, we utilize the random pinning method to suppress phonons and disentangle their coupling with non-phononic modes and successfully calculate their density of states as $g(\omega) \propto \omega^4$. We also study their localization properties and confirm that low-frequency non-phononic modes in pinned systems are truly localized without far-field contributions. We finally discuss the excess density of states over the Debye value that results from the hybridization of phonons and non-phononic modes.


title: Deep Learning based Novel Cascaded Approach for Skin Lesion Analysis

abstract: Automatic lesion analysis is critical in skin cancer diagnosis and ensures effective treatment. The computer aided diagnosis of such skin cancer in dermoscopic images can significantly reduce the clinicians workload and help improve diagnostic accuracy. Although researchers are working extensively to address this problem, early detection and accurate identification of skin lesions remain challenging. This research focuses on a two step framework for skin lesion segmentation followed by classification for lesion analysis. We explored the effectiveness of deep convolutional neural network based architectures by designing an encoder-decoder architecture for skin lesion segmentation and CNN based classification network. The proposed approaches are evaluated quantitatively in terms of the Accuracy, mean Intersection over Union and Dice Similarity Coefficient. Our cascaded end to end deep learning based approach is the first of its kind, where the classification accuracy of the lesion is significantly improved because of prior segmentation.


title: General Distribution Steering: A Sub-Optimal Solution by Convex Optimization

abstract: General distribution steering is intrinsically an infinite-dimensional problem, where the distributions to steer are arbitrary. In the literature, the distribution steering problem governed by system dynamics is usually treated by assuming the distribution of the system states to be Gaussian. In our previous paper, we considered the distribution steering problem where the initial and terminal distributions are arbitrary (only required to have first several orders of power moments), and proposed to use the moments to turn this problem into a finite-dimensional one. We put forward a moment representation of the primal system for control. However, the control law in that paper was an empirical one without optimization towards a design criterion, which doesn't always ensure a most satisfactory solution. In this paper, which is a very preliminary version, we propose a convex optimization approach to the general distribution steering problem of the first-order discrete-time linear system, i.e., an optimal control law for the corresponding moment system. The optimal control inputs of the moment system are obtained by convex optimization, of which the convexity of the domain is proved. An algorithm of distribution steering is then put forward by extending a realization scheme of control inputs using the Kullback-Leibler distance to one realization using the squared Hellinger distance, of which the performance has been shown to be better than the former one. Experiments on different types of cost functions are given to validate the performance of our proposed algorithm. Since the moment system is a dimension-reduced counterpart of the primal system, and we are not optimizing the cost function over all feasible control inputs, we call this solution a sub-optimal one to the primal general distribution steering problem.


title: An information-theoretic branch-and-prune algorithm for discrete phase optimization of RIS in massive MIMO

abstract: In this paper, we consider passive RIS-assisted multi-user communication between wireless nodes to improve the blocked line-of-sight (LOS) link performance. The wireless nodes are assumed to be equipped with Massive Multiple-Input Multiple-Output antennas, hybrid precoder, combiner, and low-resolution analog-to-digital converters (ADCs). We first derive the expression for the Cramer-Rao lower bound (CRLB) of the Mean Squared Error (MSE) of the received and combined signal at the intended receiver under interference. By appropriate design of the hybrid precoder, combiner, and RIS phase settings, it can be shown that the MSE achieves the CRLB. We further show that minimizing the MSE w.r.t. the phase settings of the RIS is equivalent to maximizing the throughput and energy efficiency of the system. We then propose a novel Information-Directed Branch-and-Prune (IDBP) algorithm to derive the phase settings of the RIS. We, for the first time in the literature, use an information-theoretic measure to decide on the pruning rules in a tree-search algorithm to arrive at the RIS phase-setting solution, which is vastly different compared to the traditional branch-and-bound algorithm that uses bounds of the cost function to define the pruning rules. In addition, we provide the theoretical guarantees of the near-optimality of the RIS phase-setting solution thus obtained using the Asymptotic Equipartition property. This also ensures near-optimal throughput and MSE performance.


title: Novelty Detection in Network Traffic: Using Survival Analysis for Feature Identification

abstract: Intrusion Detection Systems are an important component of many organizations' cyber defense and resiliency strategies. However, one downside of these systems is their reliance on known attack signatures for detection of malicious network events. When it comes to unknown attack types and zero-day exploits, modern Intrusion Detection Systems often fall short. In this paper, we introduce an unconventional approach to identifying network traffic features that influence novelty detection based on survival analysis techniques. Specifically, we combine several Cox proportional hazards models and implement Kaplan-Meier estimates to predict the probability that a classifier identifies novelty after the injection of an unknown network attack at any given time. The proposed model is successful at pinpointing PSH Flag Count, ACK Flag Count, URG Flag Count, and Down/Up Ratio as the main features to impact novelty detection via Random Forest, Bayesian Ridge, and Linear Support Vector Regression classifiers.


title: Swarm-SLAM : Sparse Decentralized Collaborative Simultaneous Localization and Mapping Framework for Multi-Robot Systems

abstract: Collaborative Simultaneous Localization And Mapping (C-SLAM) is a vital component for successful multi-robot operations in environments without an external positioning system, such as indoors, underground or underwater. In this paper, we introduce Swarm-SLAM, an open-source C-SLAM system that is designed to be scalable, flexible, decentralized, and sparse, which are all key properties in swarm robotics. Our system supports inertial, lidar, stereo, and RGB-D sensing, and it includes a novel inter-robot loop closure prioritization technique that reduces communication and accelerates convergence. We evaluated our ROS-2 implementation on five different datasets, and in a real-world experiment with three robots communicating through an ad-hoc network. Our code is publicly available: https://github.com/MISTLab/Swarm-SLAM


title: Different spin relaxation property observed in linearly and circularly polarized laser induced terahertz emission from Bi/Co bilayer

abstract: Recently, helicity-dependent photocurrent was reported in Bi single thin fi lms. It is proposed that the origin of this photocurrent is the combination of photo-spin conversion and spin-charge conversion effects in Bi and efficient spin conversion in Bi is expected. In this study, we measured two types of terahertz (THz) emissions from Bi/Co bilayer films induced by spin current generation using laser-induced demagnetization of the Co layer and photo-spin conversion effect in the Bi layer to investigate the spin current induced by the two mechanisms simultaneously. We clearly observed diff erent Bi thickness dependence of peak intensity and that of bandwidth for THz spin current in two experiments, i.e., spin current induced by demagnetization of Co and that by photo-spin conversion in Bi. The different Bi thickness dependence of spin current intensity and bandwidth in two experiments is caused by different spin relaxation properties of optically excited spin currents in Bi layers.


title: Optimal/Nearly-optimal simulation of multi-periodic time-dependent Hamiltonians

abstract: Simulating Hamiltonian dynamics is one of the most fundamental and significant tasks for characterising quantum materials. Recently, a series of quantum algorithms employing block-encoding of Hamiltonians have succeeded in providing efficient simulation of time-evolution operators on quantum computers. While time-independent Hamiltonians can be simulated by the quantum eigenvalue transformation (QET) or quantum singularvalue transformation with the optimal query complexity in time $t$ and desirable accuracy $\varepsilon$, generic time-dependent Hamiltonians face at larger query complexity and more complicated oracles due to the difficulty of handling time-dependency. In this paper, we establish a QET-based approach for simulating time-dependent Hamiltonians with multiple time-periodicity. Such time-dependent Hamiltonians involve a variety of nonequilibrium systems such as time-periodic systems (Floquet systems) and time-quasiperiodic systems. Overcoming the difficulty of time-dependency, our protocol can simulate the dynamics under multi-periodic time-dependent Hamiltonians with optimal/nearly-optimal query complexity both in time $t$ and desirable accuracy $\varepsilon$, and simple oracles as well as the optimal algorithm for time-independent cases.


title: Dimension approximation in smooth dynamical systems

abstract: For a non-conformal repeller $\Lambda$ of a $C^{1+\alpha}$ map $f$ preserving an ergodic measure $\mu$ of positive entropy, this paper shows that the Lyapunov dimension of $\mu$ can be approximated gradually by the Carath\'{e}odory singular dimension of a sequence of horseshoes. For a $C^{1+\alpha}$ diffeomorphism $f$ preserving a hyperbolic ergodic measure $\mu$ of positive entropy, if $(f, \mu)$ has only two Lyapunov exponents $\lambda_u(\mu)>0>\lambda_s(\mu)$, then the Hausdorff or lower box or upper box dimension of $\mu$ can be approximated by the corresponding dimension of the horseshoes $\{\Lambda_n\}$. The same statement holds true if $f$ is a $C^1$ diffeomorphism with a dominated Oseledet's splitting with respect to $\mu$.


title: Regularity results for mixed local and nonlocal double phase functionals

abstract: We investigate the De Giorgi-Nash-Moser theory for minimizers of mixed local and nonlocal functionals modeled after \[ v \mapsto \int_{\mathbb{R}^{n}}\int_{\mathbb{R}^{n}}\dfrac{|v(x)-v(y)|^{p}}{|x-y|^{n+sp}}\,dxdy+\int_{\Omega}a(x)|Dv|^{q}\,dx, \] where $0<s<1<p \le q$ and $a(\cdot) \ge 0$. In particular, we prove H\"older regularity and Harnack's inequality under possibly sharp assumptions on $s,p,q$ and $a(\cdot)$.


title: Almost global well-posedness of 2-D Ericksen-Leslie's hyperbolic liquid crystal model for small data

abstract: This article is devoted to the two dimensional simplified Ericksen-Leslie's hyperbolic system for incompressible liquid crystal model, where the direction $d$ of liquid crystal molecules satisfies a wave map equation with an acoustical metric. We established the almost global well-posedness for small and smooth initial data near the constant equilibrium. Our proof relies on the idea of vector-field method and ghost weight method. There are two key ingredients in our proof: (i) Inspired by the gauge theory in Tataru \cite{Tataru,Tataru05}, we reformulate the wave map equation into a free wave equation with acoustical metric, where the nonlinearity is annihilated due to the geometry of $\mathbb S^1$; (ii) Motivated by the ghost weight method in Alinhac \cite{A01}, we introduce a new and important ``good unknown", the velocity $u$, which provides the additional dissipation $u/\langle t-r\rangle\in L^2_tL^2_x$. These new observations turn out to be extremely crucial in resolving the system in low dimensions.


title: Search for Extremely Metal-Poor Stars with GEMINI-N/GRACES I. Chemical-Abundance Analysis

abstract: We present stellar parameters and abundances of 13 elements for 18 very metal-poor (VMP; [Fe/H] $<$ -2.0) stars, selected as extremely metal-poor (EMP; [Fe/H] $<$ -3.0) candidates from SDSS and LAMOST survey. High-resolution spectroscopic observations were performed using GEMINI-N/GRACES. We find ten EMP stars among our candidates, and we newly identify three carbon-enhanced metal poor (CEMP) stars with [Ba/Fe] $<$ 0. Although chemical abundances of our VMP/EMP stars generally follow the overall trend of other Galactic halo stars, there are a few exceptions. One Na-rich star ([Na/Fe] = +1.14) with low [Mg/Fe] suggests a possible chemical connection with second-generation stars in a globular cluster. The progenitor of an extremely Na-poor star ([Na/Fe] = -1.02) with an enhancement of K- and Ni-abundance ratios may have undergone a distinct nucleosynthesis episode, associated with core-collapse supernovae (CCSNe) having a high explosion energy. We have also found a Mg-rich star ([Mg/Fe] = +0.73) with slightly enhanced Na and extremely low [Ba/Fe], indicating that its origin is not associated with neutron-capture events. On the other hand, the origin of the lowest Mg abundance ([Mg/Fe] = -0.61) star could be explained by accretion from a dwarf galaxy, or formation in a gas cloud largely polluted by SNe Ia. We have also explored the progenitor masses of our EMP stars by comparing their chemical-abundance patterns with those predicted by Population III SNe models, and find a mass range of 10 - 26 $M_\odot$, suggesting that such stars were primarily responsible for the chemical enrichment of the early Milky Way.


title: A separation logic for sequences in pointer programs and its decidability

abstract: Separation logic and its variants can describe various properties on pointer programs. However, when it comes to properties on sequences, one may find it hard to formalize. To deal with properties on variable-length sequences and multilevel data structures, we propose sequence-heap separation logic which integrates sequences into logical reasoning on heap-manipulated programs. Quantifiers over sequence variables and singleton heap storing sequence (sequence singleton heap) are new members in our logic. Further, we study the satisfiability problem of two fragments. The propositional fragment of sequence-heap separation logic is decidable, and the fragment with 2 alternations on program variables and 1 alternation on sequence variables is undecidable. In addition, we explore boundaries between decidable and undecidable fragments of the logic with prenex normal form.


title: Supporting Secure Dynamic Alert Zones Using Searchable Encryption and Graph Embedding

abstract: Location-based alerts have gained increasing popularity in recent years, whether in the context of healthcare (e.g., COVID-19 contact tracing), marketing (e.g., location-based advertising), or public safety. However, serious privacy concerns arise when location data are used in clear in the process. Several solutions employ Searchable Encryption (SE) to achieve secure alerts directly on encrypted locations. While doing so preserves privacy, the performance overhead incurred is high. We focus on a prominent SE technique in the public-key setting -- Hidden Vector Encryption (HVE), and propose a graph embedding technique to encode location data in a way that significantly boosts the performance of processing on ciphertexts. We show that finding the optimal encoding is NP-hard, and provide several heuristics that are fast and obtain significant performance gains. Furthermore, we investigate the more challenging case of dynamic alert zones, where the area of interest changes over time. Our extensive experimental evaluation shows that our solutions can significantly improve computational overhead compared to existing baselines.


title: Controlling solar radiation forces with graphene in plasmonic metasurface

abstract: Controlling and harvesting solar radiation pressure is a significant challenge, however, there are few potential solutions, which are suitable for several key applications. In this study, an electrically tunable plasmonic metasurface is designed for the visible spectrum. Moreover, the normal and the tangential optical forces acting on the metasurface are calculated. Whilst presenting high efficiency in the anomalous reflection, the designed active metasurface provides tunability of optical forces acting on the metasurface. The metasurface is composed of tapered silver cells embedded on top of the graphene layer with 20 layers of graphene sheets. Hence, by tuning the Fermi level of graphene sheets, the transferred momentum to the metasurface can be controlled. Our results can provide a suitable platform for optical force control desired in tunable radiation pressure harvesting, micro vehicles, solar sailing and optical tweezers.


title: Kernel-based off-policy estimation without overlap: Instance optimality beyond semiparametric efficiency

abstract: We study optimal procedures for estimating a linear functional based on observational data. In many problems of this kind, a widely used assumption is strict overlap, i.e., uniform boundedness of the importance ratio, which measures how well the observational data covers the directions of interest. When it is violated, the classical semi-parametric efficiency bound can easily become infinite, so that the instance-optimal risk depends on the function class used to model the regression function. For any convex and symmetric function class $\mathcal{F}$, we derive a non-asymptotic local minimax bound on the mean-squared error in estimating a broad class of linear functionals. This lower bound refines the classical semi-parametric one, and makes connections to moduli of continuity in functional estimation. When $\mathcal{F}$ is a reproducing kernel Hilbert space, we prove that this lower bound can be achieved up to a constant factor by analyzing a computationally simple regression estimator. We apply our general results to various families of examples, thereby uncovering a spectrum of rates that interpolate between the classical theories of semi-parametric efficiency (with $\sqrt{n}$-consistency) and the slower minimax rates associated with non-parametric function estimation.


title: BEAGLE: Forensics of Deep Learning Backdoor Attack for Better Defense

abstract: Deep Learning backdoor attacks have a threat model similar to traditional cyber attacks. Attack forensics, a critical counter-measure for traditional cyber attacks, is hence of importance for defending model backdoor attacks. In this paper, we propose a novel model backdoor forensics technique. Given a few attack samples such as inputs with backdoor triggers, which may represent different types of backdoors, our technique automatically decomposes them to clean inputs and the corresponding triggers. It then clusters the triggers based on their properties to allow automatic attack categorization and summarization. Backdoor scanners can then be automatically synthesized to find other instances of the same type of backdoor in other models. Our evaluation on 2,532 pre-trained models, 10 popular attacks, and comparison with 9 baselines show that our technique is highly effective. The decomposed clean inputs and triggers closely resemble the ground truth. The synthesized scanners substantially outperform the vanilla versions of existing scanners that can hardly generalize to different kinds of attacks.


title: Periodic dimensions and some homological properties of eventually periodic algebras

abstract: For an eventually periodic module, we have the degree and the period of its first periodic syzygy. This paper studies the former under the name \lq\lq periodic dimension\rq\rq. We give a bound for the periodic dimension of an eventually periodic module with finite Gorenstein projective dimension. We also provide a method of computing the Gorenstein projective dimension of an eventually periodic module under certain conditions. Besides, motivated by recent results of Dotsenko, G\'elinas and Tamaroff and of the author, we determine the bimodule periodic dimension of an eventually periodic Gorenstein algebra. Another aim of this paper is to obtain some of the basic homological properties of eventually periodic algebras. We show that a lot of homological conjectures hold for this class of algebras. As an application, we characterize eventually periodic Gorenstein algebras in terms of bimodules Gorenstein projective dimensions.


title: Black Hole Ultra-compact X-ray Binaries: Galactic Low-Frequency Gravitational Wave Sources

abstract: In the Galaxy, close binaries with compact objects are important low-frequency gravitational wave (GW) sources. As potential low-frequency GW sources, neutron star/white dwarf (WD) ultra-compact X-ray binaries (UCXBs) have been investigated extensively. Using the MESA code, we systematically explored the evolution of black hole (BH)-main sequence star (MS) binaries to diagnose whether their descendants can be detected by space-borne GW detectors. Our simulations show that BH-MS binaries with an initial orbital period less than the bifurcation period can evolve into BH UCXBs that can be detected by LISA. Such an evolutionary channel would form compact mass-transferring BH-WD systems rather than detached BH-WD systems. The calculated X-ray luminosities of BH UCXBs detected by LISA at a distance $d=1$ kpc are $\sim10^{33}-10^{35}~\rm erg\,s^{-1}$ ($\sim10^{34}-10^{35}~\rm erg\,s^{-1}$ for $d=10$ kpc), hence it is possible to detect their electromagnetic counterparts. It is worth emphasizing only some BH-MS systems with an initial orbital period very close to the bifurcation period can evolve toward low-frequency GW sources whose chirp masses can be measured. The maximum GW frequency of BH UCXBs forming by BH-MS pathway is about 3 mHz, which is smaller than the minimum GW frequency (6.4 mHz) of mass-transferring BH-WD originating from a dynamic process. Furthermore, we obtain an initial parameter space (donor-star masses and orbital periods) of progenitors of BH UCXB-GW sources, which can be applied to future population synthesis simulations. By a rough estimation, we predict that LISA could detect only a few BH UCXB-GW sources forming by the BH-MS channel.


title: Haptic Transparency and Interaction Force Control for a Lower-Limb Exoskeleton

abstract: Controlling the interaction forces between a human and an exoskeleton is crucial for providing transparency or adjusting assistance or resistance levels. However, it is an open problem to control the interaction forces of lower-limb exoskeletons designed for unrestricted overground walking. For these types of exoskeletons, it is challenging to implement force/torque sensors at every contact between the user and the exoskeleton for direct force measurement. Moreover, it is important to compensate for the exoskeleton's whole-body gravitational and dynamical forces, especially for heavy lower-limb exoskeletons. Previous works either simplified the dynamic model by treating the legs as independent double pendulums, or they did not close the loop with interaction force feedback. The proposed whole-exoskeleton closed-loop compensation (WECC) method calculates the interaction torques during the complete gait cycle by using whole-body dynamics and joint torque measurements on a hip-knee exoskeleton. Furthermore, it uses a constrained optimization scheme to track desired interaction torques in a closed loop while considering physical and safety constraints. We evaluated the haptic transparency and dynamic interaction torque tracking of WECC control on three subjects. We also compared the performance of WECC with a controller based on a simplified dynamic model and a passive version of the exoskeleton. The WECC controller results in a consistently low absolute interaction torque error during the whole gait cycle for both zero and nonzero desired interaction torques. In contrast, the simplified controller yields poor performance in tracking desired interaction torques during the stance phase.


title: Deformations of $\mathbb Z_2$-Harmonic Spinors on 3-Manifolds

abstract: A $\mathbb Z_2$-harmonic spinor on a 3-manifold $Y$ is a solution of the Dirac equation on a bundle that is twisted around a submanifold $\mathcal Z$ of codimension 2 called the singular set. This article investigates the local structure of the universal moduli space of $\mathbb Z_2$-harmonic spinors over the space of parameters $(g,B)$ consisting of a metric and perturbation to the spin connection. The main result states that near a $\mathbb Z_2$-harmonic spinor with $\mathcal Z$ smooth, the universal moduli space projects to a codimension 1 submanifold in the space of parameters. The analysis is complicated by the presence of an infinite-dimensional obstruction bundle and a loss of regularity in the first variation of the Dirac operator with respect to deformations of the singular set $\mathcal Z$, necessitating the use of the Nash-Moser Implicit Function Theorem.


title: Mobility Data in Operations: The Facility Location Problem

abstract: The recent large scale availability of mobility data, which captures individual mobility patterns, poses novel operational problems that are exciting and challenging. Motivated by this, we introduce and study a variant of the (cost-minimization) facility location problem where each individual is endowed with two locations (hereafter, her home and work locations), and the connection cost is the minimum distance between any of her locations and its closest facility. We design a polynomial-time algorithm whose approximation ratio is at most 2.497. We complement this positive result by showing that the proposed algorithm is at least a 2.428-approximation, and there exists no polynomial-time algorithm with approximation ratio $2-\epsilon$ under UG-hardness. We further extend our results and analysis to the model where each individual is endowed with K locations. Finally, we conduct numerical experiments over both synthetic data and US census data (for NYC, greater LA, greater DC, Research Triangle) and evaluate the performance of our algorithms.


title: Circle action of the punctured mapping class group and cross homomorphism

abstract: In the following short note, we give a new geometric interpretation of the generator of the infinite cyclic group $H^1(\text{Mod}(S_{g,1});H^1(S_g;\mathbb{Z}))$ (this computation is proved by Morita). There are several construction of this class given by Earle, Morita, Trapp and Furuta. The construction we give here uses the action of $\text{Mod}(S_{g,1})$ on the circle and its rotation numbers. We suspect that our construction is the same as the construction by Furuta and Trapp using winding numbers and provide half of the proof.


title: Unconventional polarization switching mechanism in (Hf, Zr)O2 ferroelectrics

abstract: HfO$_{2}$-based ferroelectric thin films are promising for their application in ferroelectric devices. Predicting the ultimate magnitude of polarization and understanding its switching mechanism are critical to realize the optimal performance of these devices. Here, a generalized solid-state variable cell nudged elastic band (VCNEB) method is employed to predict the switching pathway associated with domain-wall motion in (Hf, Zr)O$_{2}$ ferroelectrics. It is found that the polarization reversal pathway, where three-fold coordinated O atoms pass across the nominal unit-cell boundaries defined by the Hf/Zr atomic planes, is energetically more favorable than the conventional pathway where the O atoms do not pass through these planes. This finding implies that the polarization orientation in the orthorhombic Pca2$_{1}$ phase of HfO$_{2}$ nd its derivatives is opposite to that normally assumed, predicts the spontaneous polarization magnitude of about 70 ${\mu}$C/cm$^{2}$ that is nearly 50% larger than the commonly accepted value, signifies a positive intrinsic longitudinal piezoelectric coefficient, and suggests growth of ferroelectric domains, in response to an applied electric field, structurally reversed to those usually anticipated. These results provide important insights into the understanding of ferroelectricity in HfO$_{2}$-based ferroelectrics.


title: DisPad: Flexible On-Body Displacement of Fabric Sensors for Robust Joint-Motion Tracking

abstract: The last few decades have witnessed an emerging trend of wearable soft sensors; however, there are important signal-processing challenges for soft sensors that still limit their practical deployment. They are error-prone when displaced, resulting in significant deviations from their ideal sensor output. In this work, we propose a novel prototype that integrates an elbow pad with a sparse network of soft sensors. Our prototype is fully bio-compatible, stretchable, and wearable. We develop a learning-based method to predict the elbow orientation angle and achieve an average tracking error of 9.82 degrees for single-user multi-motion experiments. With transfer learning, our method achieves the average tracking errors of 10.98 degrees and 11.81 degrees across different motion types and users, respectively. Our core contributions lie in a solution that realizes robust and stable human joint motion tracking across different device displacements.


title: High sensitivity silicon carbide divacancy-based thermometer

abstract: Color centers in silicon carbide have become potentially versatile quantum sensors. Particularly, wide temperature range temperature sensing has been realized in recent years. However, the sensitivity is limited due to the short dephasing time of the color centers. In this work, we realize a high sensitivity silicon carbide divacancy-based thermometer using the thermal Carr-Purcell-Meiboom-Gill (TCPMG) method. First, the zero field splitting D of PL6 divacancy as a function of temperature is measured with a linear slope of -99.7 kHz/K. The coherence times of TCPMG pulses linearly increase with the pulse number and the longest coherence time is about 21 us, which is ten times larger than dephasing time. The corresponding temperature sensing sensitivity is 13.4 mK/Hz1/2, which is about 15 times higher than previous results. Finally, we monitor the laboratory temperature variations for 24 hours using the TCMPG pulse. The experiments pave the way for the applications of silicon carbide-based high sensitivity thermometer in the semiconductor industry, biology, and materials sciences.


title: Machine Learning-Aided Efficient Decoding of Reed-Muller Subcodes

abstract: Reed-Muller (RM) codes achieve the capacity of general binary-input memoryless symmetric channels and are conjectured to have a comparable performance to that of random codes in terms of scaling laws. However, such results are established assuming maximum-likelihood decoders for general code parameters. Also, RM codes only admit limited sets of rates. Efficient decoders such as successive cancellation list (SCL) decoder and recently-introduced recursive projection-aggregation (RPA) decoders are available for RM codes at finite lengths. In this paper, we focus on subcodes of RM codes with flexible rates. We first extend the RPA decoding algorithm to RM subcodes. To lower the complexity of our decoding algorithm, referred to as subRPA, we investigate different approaches to prune the projections. Next, we derive the soft-decision based version of our algorithm, called soft-subRPA, that not only improves upon the performance of subRPA but also enables a differentiable decoding algorithm. Building upon the soft-subRPA algorithm, we then provide a framework for training a machine learning (ML) model to search for \textit{good} sets of projections that minimize the decoding error rate. Training our ML model enables achieving very close to the performance of full-projection decoding with a significantly smaller number of projections. We also show that the choice of the projections in decoding RM subcodes matters significantly, and our ML-aided projection pruning scheme is able to find a \textit{good} selection, i.e., with negligible performance degradation compared to the full-projection case, given a reasonable number of projections.


title: Estimating the feasibility of `standard speed-gun' distances

abstract: In a previous paper, we demonstrated a single-rung method for measuring cosmological distances in active galactic nuclei (AGN) that can be used from low redshift (z < 0.1) to high redshift (z > 3). This method relies on the assumption that the variability seen in AGN is constrained by the speed of light during a flare event and can therefore be used to estimate the size of an emitting region. A limitation of this method is that previously, the Doppler factor was required to be known. In this paper, we derive an extension of the `standard speed-gun' method for measuring cosmological distances that depends on the maximum intrinsic brightness temperature that a source can reach, rather than the Doppler factor. If the precise value of the intrinsic brightness temperature does not evolve with redshift and flares are statistically independent, we can in principle improve the errors in measurements of the matter content of the universe (in a flat LambdaCDM model) statistically. We then explored how well a future observing program would constrain cosmological parameters. We found that recovering the input cosmology depends critically on the uncertainty of the intrinsic brightness temperature and the number of flares observed.


title: Switching integrators reversibly in the astrophysical $N$-body problem

abstract: We present a simple algorithm to switch between $N$-body time integrators in a reversible way. We apply it to planetary systems undergoing arbitrarily close encounters and highly eccentric orbits, but the potential applications are broader. Upgrading an ordinary non-reversible switching integrator to a reversible one is straightforward and introduces no appreciable computational burden in our tests. Our method checks if the integrator during the time step violates a time-symmetric selection condition and redoes the step if necessary. In our experiments a few percent of steps would have violated the condition without our corrections. By eliminating them the algorithm avoids long-term error accumulation, of several orders magnitude in some cases.

categories: astro-ph.EP astro-ph.GA astro-ph.IM physics.comp-ph


title: Optical conductivity of the threefold Hopf semimetal

abstract: A multifold Hopf semimetal is a topological point node semimetal possessing an anisotropy in the internal electronic structure, e.g., the dipole structure of the Berry curvature. In this paper, the unique features of threefold Hopf semimetals in terms of the optical conductivity are theoretically investigated with a minimal theoretical model by using linear response theory for a linearly polarized photon. The frequency spectrum of the optical conductivity shows an anisotropic dependence on the polarization angle of the incident photon even if the electronic band structure is completely isotropic. The longitudinal optical conductivity linearly depends on the photon frequency and possesses step-like changes in the frequency spectrum. The anisotropic electronic structure has a varying the number of steps with the orientation of the photon polarization axis. We reveal that the anisotropy is attributed to symmetries preserving the point node in threefold Hopf semimetals. The linearly polarized photon also induces a Hall current but it vanishes with the photon polarization axis parallel to the Berry dipole axis. The numerical calculations show that these characteristic features can be observed even with a non-zero temperature and disorder.


title: Stability of time-periodic $\mathcal{PT}$ and anti-$\mathcal{PT}$-symmetric Hamiltonians with different periodicities

abstract: Hermitian Hamiltonians with time-periodic coefficients can be analyzed via Floquet theory, and have been extensively used for engineering Floquet Hamiltonians in standard quantum simulators. Generalized to non-Hermitian Hamiltonians, time-periodicity offers avenues to engineer the landscape of Floquet quasi-energies across the complex plane. We investigate two-level non-Hermitian Hamiltonians with coefficients that have different periodicities using Floquet theory. By analytical and numerical calculations, we obtain their regions of stability, defined by real Floquet quasi-energies, and contours of exceptional point (EP) degeneracies. We extend our analysis to study the phases that accompany the cyclic changes. Our results demonstrate that time-periodic, non-Hermitian Hamiltonians generate a rich landscape of stable and unstable regions.


title: A low Mach enthalpy method to model non-isothermal gas-liquid-solid flows with melting and solidification

abstract: Modeling phase change problems numerically is vital for understanding many natural (e.g., ice formation, steam generation) and engineering processes (e.g., casting, welding, additive manufacturing). Almost all phase change materials (PCMs) exhibit density/volume changes during melting, solidification, boiling, or condensation, causing additional fluid flow during this transition. Most numerical works consider only two phase flows (either solid-liquid or liquid-gas) for modeling phase change phenomena and some also neglect volume/density change of PCMs in the models. This paper presents a novel low Mach enthalpy method for simulating solidification and melting problems with variable thermophysical properties, including density. Additionally, this formulation allows coupling a solid-liquid PCM with a gas phase in order to simulate the free surface dynamics of PCMs undergoing melting and solidification. We revisit the two-phase Stefan problem involving a density jump between two material phases. We propose a possible means to include the kinetic energy jump in the Stefan condition while still allowing for an analytical solution. The new low Mach enthalpy method is validated against analytical solutions for a PCM undergoing a large density change during its phase transition. Additionally, a few simple sanity checks are proposed to benchmark computational fluid dynamics (CFD) algorithms that aim to capture the volume change effects of PCMs.


title: On the complexity of analyticity in semi-definite optimization

abstract: It is well-known that the central path of semi-definite optimization, unlike linear optimization, has no analytic extension to $\mu = 0$ in the absence of the strict complementarity condition. In this paper, we show the existence of a positive integer $\rho$ by which the reparametrization $\mu \mapsto \mu^{\rho}$ recovers the analyticity of the central path at $\mu = 0$. We investigate the complexity of computing $\rho$ using algorithmic real algebraic geometry and the theory of complex algebraic curves. We prove that the optimal $\rho$ is bounded by $2^{O(m^2+n^2m+n^4)}$, where $n$ is the matrix size and $m$ is the number of affine constraints. Our approach leads to a symbolic algorithm, based on the Newton-Puiseux algorithm, which computes a feasible $\rho$ using $2^{O(m+n^2)}$ arithmetic operations.


title: Attractors for the Navier-Stokes-Cahn-Hilliard System with Chemotaxis and Singular Potential in 2D

abstract: We analyze the long-time behavior of solutions to a Navier-Stokes-Cahn--Hilliard-Oono system with chemotaxis effects and physically relevant singular potential. This model also includes some significant mechanisms such as active transport and chemotaxis effects. The system couples the Navier-Stokes equations for the fluid velocity, a convective Cahn-Hilliard equation for the phase-field variable and a diffusion equation for the nutrient density. For the initial boundary value problem in a smooth bounded domain $\Omega\subset \mathbb{R}^2$, we prove the the existence of the global attractor in a suitable phase space. Furthermore, we obtain the existence of an exponential attractor, and it can be implied that the global attractor is of finite fractal dimension.


title: Understanding Best Subset Selection: A Tale of Two C(omplex)ities

abstract: For decades, best subset selection (BSS) has eluded statisticians mainly due to its computational bottleneck. However, until recently, modern computational breakthroughs have rekindled theoretical interest in BSS and have led to new findings. Recently, \cite{guo2020best} showed that the model selection performance of BSS is governed by a margin quantity that is robust to the design dependence, unlike modern methods such as LASSO, SCAD, MCP, etc. Motivated by their theoretical results, in this paper, we also study the variable selection properties of best subset selection for high-dimensional sparse linear regression setup. We show that apart from the identifiability margin, the following two complexity measures play a fundamental role in characterizing the margin condition for model consistency: (a) complexity of \emph{residualized features}, (b) complexity of \emph{spurious projections}. In particular, we establish a simple margin condition that depends only on the identifiability margin and the dominating one of the two complexity measures. Furthermore, we show that a margin condition depending on similar margin quantity and complexity measures is also necessary for model consistency of BSS. For a broader understanding, we also consider some simple illustrative examples to demonstrate the variation in the complexity measures that refines our theoretical understanding of the model selection performance of BSS under different correlation structures.