title: Distributed Contrastive Learning for Medical Image Segmentation

abstract: Supervised deep learning needs a large amount of labeled data to achieve high performance. However, in medical imaging analysis, each site may only have a limited amount of data and labels, which makes learning ineffective. Federated learning (FL) can learn a shared model from decentralized data. But traditional FL requires fully-labeled data for training, which is very expensive to obtain. Self-supervised contrastive learning (CL) can learn from unlabeled data for pre-training, followed by fine-tuning with limited annotations. However, when adopting CL in FL, the limited data diversity on each site makes federated contrastive learning (FCL) ineffective. In this work, we propose two federated self-supervised learning frameworks for volumetric medical image segmentation with limited annotations. The first one features high accuracy and fits high-performance servers with high-speed connections. The second one features lower communication costs, suitable for mobile devices. In the first framework, features are exchanged during FCL to provide diverse contrastive data to each site for effective local CL while keeping raw data private. Global structural matching aligns local and remote features for a unified feature space among different sites. In the second framework, to reduce the communication cost for feature exchanging, we propose an optimized method FCLOpt that does not rely on negative samples. To reduce the communications of model download, we propose the predictive target network update (PTNU) that predicts the parameters of the target network. Based on PTNU, we propose the distance prediction (DP) to remove most of the uploads of the target network. Experiments on a cardiac MRI dataset show the proposed two frameworks substantially improve the segmentation and generalization performance compared with state-of-the-art techniques.


title: Effects of network topology and trait distribution on collective decision making

abstract: Social networks play an important role in analyzing the impact of individual-level interactions on societal or economic outcomes. We model interactive decision making for a community of individuals with different traits, represented by a social network with trait-attributed nodes. We develop a deterministic process generating a sequence of choices for each individual based on a trait-attributed social network, initial choices of individuals and a set of predetermined trait-dependent rules for making decisions. The object of interest is the sequence of cumulative sum of choices over all individuals, which we call the cumulative sequence and consider as an index of collective decisions. We observe that, in a time period, a cumulative sequence can be unpredictable or predictable showing a repeated pattern either escalating to an extreme or constantly oscillating. We consider that predictable cumulative sequences represent unstable collective decisions of communities either extremizing or internally conflicting, while unpredictable cumulative sequences show stable changes. We analyze the effects of network topology and trait distribution on the probability of cumulative sequences being predictable, escalating and oscillating by simulations. Our findings include that unstable collective decisions are more probable as network density increases, that centralized networks are more likely to have unstable collective decisions and that networks with excessively clustered or scattered conformists and rebels tend to produce unstable cumulative sequences. We discuss the potential of the model as a framework for studying individuals with different traits on a social network directly and indirectly interacting in decision making.


title: Adaptivity Gaps for the Stochastic Boolean Function Evaluation Problem

abstract: We consider the Stochastic Boolean Function Evaluation (SBFE) problem where the task is to efficiently evaluate a known Boolean function $f$ on an unknown bit string $x$ of length $n$. We determine $f(x)$ by sequentially testing the variables of $x$, each of which is associated with a cost of testing and an independent probability of being true. If a strategy for solving the problem is adaptive in the sense that its next test can depend on the outcomes of previous tests, it has lower expected cost but may take up to exponential space to store. In contrast, a non-adaptive strategy may have higher expected cost but can be stored in linear space and benefit from parallel resources. The adaptivity gap, the ratio between the expected cost of the optimal non-adaptive and adaptive strategies, is a measure of the benefit of adaptivity. We present lower bounds on the adaptivity gap for the SBFE problem for popular classes of Boolean functions, including read-once DNF formulas, read-once formulas, and general DNFs. Our bounds range from $\Omega(\log n)$ to $\Omega(n/\log n)$, contrasting with recent $O(1)$ gaps shown for symmetric functions and linear threshold functions.


title: Decomposable Non-Smooth Convex Optimization with Nearly-Linear Gradient Oracle Complexity

abstract: Many fundamental problems in machine learning can be formulated by the convex program \[ \min_{\theta\in R^d}\ \sum_{i=1}^{n}f_{i}(\theta), \] where each $f_i$ is a convex, Lipschitz function supported on a subset of $d_i$ coordinates of $\theta$. One common approach to this problem, exemplified by stochastic gradient descent, involves sampling one $f_i$ term at every iteration to make progress. This approach crucially relies on a notion of uniformity across the $f_i$'s, formally captured by their condition number. In this work, we give an algorithm that minimizes the above convex formulation to $\epsilon$-accuracy in $\widetilde{O}(\sum_{i=1}^n d_i \log (1 /\epsilon))$ gradient computations, with no assumptions on the condition number. The previous best algorithm independent of the condition number is the standard cutting plane method, which requires $O(nd \log (1/\epsilon))$ gradient computations. As a corollary, we improve upon the evaluation oracle complexity for decomposable submodular minimization by Axiotis et al. (ICML 2021). Our main technical contribution is an adaptive procedure to select an $f_i$ term at every iteration via a novel combination of cutting-plane and interior-point methods.


title: When can I Speak? Predicting initiation points for spoken dialogue agents

abstract: Current spoken dialogue systems initiate their turns after a long period of silence (700-1000ms), which leads to little real-time feedback, sluggish responses, and an overall stilted conversational flow. Humans typically respond within 200ms and successfully predicting initiation points in advance would allow spoken dialogue agents to do the same. In this work, we predict the lead-time to initiation using prosodic features from a pre-trained speech representation model (wav2vec 1.0) operating on user audio and word features from a pre-trained language model (GPT-2) operating on incremental transcriptions. To evaluate errors, we propose two metrics w.r.t. predicted and true lead times. We train and evaluate the models on the Switchboard Corpus and find that our method outperforms features from prior work on both metrics and vastly outperforms the common approach of waiting for 700ms of silence.


title: Correlation of Anomaly Rates in the Slovak Electric Transmission Grid with Geomagnetic Activity

abstract: The induction of electric currents in electric power distribution networks is a well-known effect of Earth-directed eruptive events. Inspired by recent studies showing that the rate of power-grid anomalies may increase after exposure to strong geomagnetically induced currents in mid-latitude countries in the middle of Europe, we decided to investigate such effects for Slovak electric-power distribution grid managed by the Slovak Electricity Transmission System (SEPS). We obtained a list of disturbances recorded in the maintenance logs by the SEPS company technicians with their dates and other details. Working with unadjusted data we found unexpectedly strong seasonal variations of failure rates. We show that this is a consequence of the strong dependence on meteorological conditions. In particular, mean temperature showed the strongest correlation with anomaly rate, while precipitation amount was just a little weaker in terms of the correlation. After the suppression of the meteorological dependence in the failure-rate series, we obtained a new temperature-independent dataset. Using this clean dataset we found an increase in failure rates in the periods of maxima of geomagnetic activity compared to the adjacent minima of activity.


title: Topological States in Two-Dimensional Su-Schrieffer-Heeger Models

abstract: We study the topological properties of the generalized two-dimensional (2D) Su-Schrieffer-Heeger (SSH) models. We show that a pair of Dirac points appear in the Brillouin zone (BZ), consisting a semimetallic phase. Interestingly, the locations of these Dirac points are not pinned to any high-symmetry points of the BZ but tunable by model parameters. Moreover, the merging of two Dirac points undergoes a novel topological phase transition, which leads to either a weak topological insulator or a nodal-line metallic phase. We demonstrate these properties by constructing two specific models, which we referred as type-I and type-II 2D SSH models. The feasible experimental platforms to realize our models are also discussed.


title: Selection on moral hazard in the Swiss market for mandatory health insurance: Empirical evidence from Swiss Household Panel data

abstract: Selection on moral hazard represents the tendency to select a specific health insurance coverage depending on the heterogeneity in utilisation ''slopes''. I use data from the Swiss Household Panel and from publicly available regulatory data to explore the extent of selection on slopes in the Swiss managed competition system. I estimate responses in terms of (log) doctor visits to lowest and highest deductible levels using Roy-type models, identifying marginal treatment effects with local instrumental variables. The response to high coverage plans (i.e. plans with the lowest deductible level) among high moral hazard types is 25-35 percent higher than average.


title: Can accreting primordial black holes explain the excess radio background?

abstract: The excess radio background seen at $\simeq 0.1-10\,{\rm GHz}$ has stimulated much scientific debate in the past years. Recently, it was pointed out that the soft photon emission from accreting primordial black holes may be able to explain this signal. We show that the expected ultraviolet photon emission from these accreting black holes would ionize the universe completely at $z>6$ and thus wash out the 21 cm absorption signature at $z\simeq$ 20 as well as be in tension with existing cosmic microwave background anisotropy and average spectral distortion limits. We discuss possible augmentations of the model; however, it seems that an explanation of radio excess by accreting primordial black holes is not well-justified.


title: Programs as Diagrams: From Categorical Computability to Computable Categories

abstract: This is a draft of the textbook/monograph that presents computability theory using string diagrams. The introductory chapters have been taught as graduate and undergraduate courses and evolved through 8 years of lecture notes. The later chapters contain new ideas and results about categorical computability and some first steps into computable category theory. The underlying categorical view of computation is based on monoidal categories with program evaluators, called *monoidal computers*. This categorical structure can be viewed as a single-instruction diagrammatic programming language called Run, whose only instruction is called RUN. This version: improved text, moved the final chapter to the next volume. (The final version will continue lots of exercises and workouts, but already this version has severely degraded graphics to meet the size bounds.)


title: Lipschitz means and mixers on metric spaces

abstract: The standard arithmetic measures of center, the mean and median, have natural topological counterparts which have been widely used in continuum theory. In the context of metric spaces it is natural to consider the Lipschitz continuous versions of the mean and median. We show that they are related to familiar concepts of the geometry of metric spaces: the bounded turning property, the existence of quasisymmetric parameterization, and others.


title: Cross-Skeleton Interaction Graph Aggregation Network for Representation Learning of Mouse Social Behaviour

abstract: Automated social behaviour analysis of mice has become an increasingly popular research area in behavioural neuroscience. Recently, pose information (i.e., locations of keypoints or skeleton) has been used to interpret social behaviours of mice. Nevertheless, effective encoding and decoding of social interaction information underlying the keypoints of mice has been rarely investigated in the existing methods. In particular, it is challenging to model complex social interactions between mice due to highly deformable body shapes and ambiguous movement patterns. To deal with the interaction modelling problem, we here propose a Cross-Skeleton Interaction Graph Aggregation Network (CS-IGANet) to learn abundant dynamics of freely interacting mice, where a Cross-Skeleton Node-level Interaction module (CS-NLI) is used to model multi-level interactions (i.e., intra-, inter- and cross-skeleton interactions). Furthermore, we design a novel Interaction-Aware Transformer (IAT) to dynamically learn the graph-level representation of social behaviours and update the node-level representation, guided by our proposed interaction-aware self-attention mechanism. Finally, to enhance the representation ability of our model, an auxiliary self-supervised learning task is proposed for measuring the similarity between cross-skeleton nodes. Experimental results on the standard CRMI13-Skeleton and our PDMB-Skeleton datasets show that our proposed model outperforms several other state-of-the-art approaches.


title: Gamma-ray flares from pulsar wind nebulae in the Large Magellanic Cloud

abstract: High-energy radiation of young pulsar wind nebulae (PWNe) is known to be variable. This is most prominently exemplified by the Crab nebula which can undergo both rapid brightenings and dimmings. Two pulsars in the Large Magellanic Cloud, PSR J0540-6919 and PSR J0537-6910 are evolutionally very close to Crab, so one may expect the same kind of variability from the PWNe around them as from the Crab nebula. In this work we search for variability in these PWNe in gamma rays using the data from the Fermi Large Area Telescope in the range 100 MeV-10 GeV collected from August 2008 to December 2021. We construct light curves of these sources in the three bands, 100-300 MeV, 300-1000 MeV and 1-10 GeV with the one week time resolution. We find evidence of flaring activity in all the bands, in contrast with Crab, where no flares at $E>1$ GeV were observed. Analysis of the flaring episode in the 100-300 and 300-1000 MeV bands indicates that the flux of one of the PWNe could grew by a factor of $\approx 5-10$. We are not confident about which of the two PWNe flared because of their proximity in the sky. However, in the 1-10 GeV band where the angular resolution of the LAT is better, we find several episodes of enhanced brightness in both the PWNe. We check possible contaminants which could be responsible for the observed variability, but find their contribution not to be relevant.


title: Norm inequalities for Dunkl-type fractional integral and fractional maximal operators in the Dunkl-Fofana spaces

abstract: We establish some new properties of the Dunkl-Wiener amalgam spaces defined on the real line. These results allow us to obtain the boundedness of Dunkl-type fractional integral and fractional maximal operators in the Dunkl-Fofana spaces.


title: Garbled EDA: Privacy Preserving Electronic Design Automation

abstract: The complexity of modern integrated circuits (ICs) necessitates collaboration between multiple distrusting parties, including thirdparty intellectual property (3PIP) vendors, design houses, CAD/EDA tool vendors, and foundries, which jeopardizes confidentiality and integrity of each party's IP. IP protection standards and the existing techniques proposed by researchers are ad hoc and vulnerable to numerous structural, functional, and/or side-channel attacks. Our framework, Garbled EDA, proposes an alternative direction through formulating the problem in a secure multi-party computation setting, where the privacy of IPs, CAD tools, and process design kits (PDKs) is maintained. As a proof-of-concept, Garbled EDA is evaluated in the context of simulation, where multiple IP description formats (Verilog, C, S) are supported. Our results demonstrate a reasonable logical-resource cost and negligible memory overhead. To further reduce the overhead, we present another efficient implementation methodology, feasible when the resource utilization is a bottleneck, but the communication between two parties is not restricted. Interestingly, this implementation is private and secure even in the presence of malicious adversaries attempting to, e.g., gain access to PDKs or in-house IPs of the CAD tool providers.


title: Automatically Finding Optimal Index Structure

abstract: Existing learned indexes (e.g., RMI, ALEX, PGM) optimize the internal regressor of each node, not the overall structure such as index height, the size of each layer, etc. In this paper, we share our recent findings that we can achieve significantly faster lookup speed by optimizing the structure as well as internal regressors. Specifically, our approach (called AirIndex) expresses the end-to-end lookup time as a novel objective function, and searches for optimal design decisions using a purpose-built optimizer. In our experiments with state-of-the-art methods, AirIndex achieves 3.3x-7.7x faster lookup for the data stored on local SSD, and 1.4x-3.0x faster lookup for the data on Azure Cloud Storage.


title: Towards Graph Representation Learning Based Surgical Workflow Anticipation

abstract: Surgical workflow anticipation can give predictions on what steps to conduct or what instruments to use next, which is an essential part of the computer-assisted intervention system for surgery, e.g. workflow reasoning in robotic surgery. However, current approaches are limited to their insufficient expressive power for relationships between instruments. Hence, we propose a graph representation learning framework to comprehensively represent instrument motions in the surgical workflow anticipation problem. In our proposed graph representation, we maps the bounding box information of instruments to the graph nodes in the consecutive frames and build inter-frame/inter-instrument graph edges to represent the trajectory and interaction of the instruments over time. This design enhances the ability of our network on modeling both the spatial and temporal patterns of surgical instruments and their interactions. In addition, we design a multi-horizon learning strategy to balance the understanding of various horizons indifferent anticipation tasks, which significantly improves the model performance in anticipation with various horizons. Experiments on the Cholec80 dataset demonstrate the performance of our proposed method can exceed the state-of-the-art method based on richer backbones, especially in instrument anticipation (1.27 v.s. 1.48 for inMAE; 1.48 v.s. 2.68 for eMAE). To the best of our knowledge, we are the first to introduce a spatial-temporal graph representation into surgical workflow anticipation.


title: Integrated Metasurface-based Wavelengths Division Demultiplexers

abstract: We present a design approach for realizing on-chip wavelength division demultiplexing (WDD) schemes by integrating all-dielectric metasurfaces of TiO2 nanorod arrays into a SiN waveguide. The designed metasurface locally modifies the effective refractive index of the SiN waveguide, creating an effective WDD that selectively passes a certain band of wavelengths into a specific output port. A set of representative 2-channel and 3-channel WDDs schemes were demonstrated for input TE00/TM00 modes and operating in different bands, showing the flexibility of our design approach. The proposed WDD schemes are compatible with visible to infrared wavelengths, photolithography-based fabrication, high efficiency with maximum transmission of 91%, and a small footprint at a few microns. Our design method paves the way for realizing several on-chip integrated devices for applications in optical data processing and biological sensing.


title: Fine-Grained Egocentric Hand-Object Segmentation: Dataset, Model, and Applications

abstract: Egocentric videos offer fine-grained information for high-fidelity modeling of human behaviors. Hands and interacting objects are one crucial aspect of understanding a viewer's behaviors and intentions. We provide a labeled dataset consisting of 11,243 egocentric images with per-pixel segmentation labels of hands and objects being interacted with during a diverse array of daily activities. Our dataset is the first to label detailed hand-object contact boundaries. We introduce a context-aware compositional data augmentation technique to adapt to out-of-distribution YouTube egocentric video. We show that our robust hand-object segmentation model and dataset can serve as a foundational tool to boost or enable several downstream vision applications, including hand state classification, video activity recognition, 3D mesh reconstruction of hand-object interactions, and video inpainting of hand-object foregrounds in egocentric videos. Dataset and code are available at: https://github.com/owenzlz/EgoHOS


title: On the Cosmic Web Elongation in Fuzzy Dark Matter Cosmologies: Effects on Density Profiles, Shapes and Alignments of Halos

abstract: The fuzzy dark matter (FDM) scenario has received increased attention in recent years due to the small-scale challenges of the vanilla Lambda cold dark matter ($\Lambda$CDM) cosmological model and the lack of any experimental evidence for any candidate particle. In this study, we use cosmological $N$-body simulations to investigate high-redshift dark matter halos and their responsiveness to an FDM-like power spectrum cutoff on small scales in the primordial density perturbations. We study halo density profiles, shapes and alignments in FDM-like cosmologies (the latter two for the first time) by providing fits and quantifying departures from $\Lambda$CDM as a function of the particle mass $m$. Compared to $\Lambda$CDM, the concentrations of FDM-like halos are lower, peaking at an $m$-dependent halo mass and thus breaking the approximate universality of density profiles in $\Lambda$CDM. The intermediate-to-major and minor-to-major shape parameter profiles are monotonically increasing with ellipsoidal radius in $N$-body simulations of $\Lambda$CDM. In FDM-like cosmologies, the monotonicity is broken, halos are more elongated around the virial radius than their $\Lambda$CDM counterparts and less elongated closer to the center. Finally, intrinsic alignment correlations, stemming from the deformation of initially spherically collapsing halos in an ambient gravitational tidal field, become stronger with decreasing $m$. At $z\sim 4$, we find a $6.4 \sigma$-significance in the fractional differences between the isotropised linear alignment magnitudes $D_{\text{iso}}$ in the $m=10^{-22}$ eV model and $\Lambda$CDM. Such FDM-like imprints on the internal properties of virialised halos are expected to be strikingly visible in the high-$z$ Universe.


title: How Do Viewers Synthesize Conflicting Information from Data Visualizations?

abstract: Scientific knowledge develops through cumulative discoveries that build on, contradict, contextualize, or correct prior findings. Scientists and journalists often communicate these incremental findings to lay people through visualizations and text (e.g., the positive and negative effects of caffeine intake). Consequently, readers need to integrate diverse and contrasting evidence from multiple sources to form opinions or make decisions. However, the underlying mechanism for synthesizing information from multiple visualizations remains underexplored. To address this knowledge gap, we conducted a series of four experiments (N = 1166) in which participants synthesized empirical evidence from a pair of line charts presented sequentially. In Experiment 1, we administered a baseline condition with charts depicting no specific context where participants held no strong belief. To test for the generalizability, we introduced real-world scenarios to our visualizations in Experiment 2, and added accompanying text descriptions similar to on-line news articles or blog posts in Experiment 3. In all three experiments, we varied the relative direction and magnitude of line slopes within the chart pairs. We found that participants tended to weigh the positive slope more when the two charts depicted relationships in the opposite direction (e.g., one positive slope and one negative slope). Participants tended to weigh the less steep slope when the two charts depicted relationships in the same direction (e.g., both positive). Through these experiments, we characterize participants' synthesis behaviors depending on the relationship between the information they viewed, contribute to theories describing underlying cognitive mechanisms in information synthesis, and describe design implications for data storytelling.


title: Revisiting Random Points: Combinatorial Complexity and Algorithms

abstract: Consider a set $P$ of $n$ points picked uniformly and independently from $[0,1]^d$ for a constant dimension $d$ -- such a point set is extremely well behaved in many aspects. For example, for a fixed $r \in [0,1]$, we prove a new concentration result on the number of pairs of points of $P$ at a distance at most $r$ -- we show that this number lies in an interval that contains only $O(n \log n)$ numbers. We also present simple linear time algorithms to construct the Delaunay triangulation, Euclidean MST, and the convex hull of the points of $P$. The MST algorithm is an interesting divide-and-conquer algorithm which might be of independent interest. We also provide a new proof that the expected complexity of the Delaunay triangulation of $P$ is linear -- the new proof is simpler and more direct, and might be of independent interest. Finally, we present a simple $\tilde{O}(n^{4/3})$ time algorithm for the distance selection problem for $d=2$.


title: A family of quintic Thue equations via Skolem's $p$-adic method

abstract: In this semi-expository article we solve the diophantine equation $m^5+(4 \cdot 5^4 b^4)mn^4 - n^5=1$ for all integers $b \neq 0$. This gives an example of a family of quintic Thue equations that can be solved completely by using nothing more than Skolem's $p$-adic method. We also give a general introduction to Skolem's method from a modern perspective.


title: Limits on Wave Optics Simulations of Plane Wave Propagation in Non-Kolmogorov Turbulence

abstract: We derive limits for wave optics simulations of plane wave propagation in non-Kolmogorov turbulence using the split-step method and thin phase screens. These limits are used to inform two simulation campaigns where the relationship between volume turbulence strength and normalized intensity variance for various non-Kolmogorov power-laws. We find that simulations of smaller power laws are limited turbulence strengths with Rytov numbers of 7 when the simulation side-length sampling rate is 8192. Under these same conditions it is possible to simulate volume turbulence strength out to Rytov numbers of 12 for Kolmogorov power-laws and higher. We show that the peak scintillation and Rytov number where peak scintillation occur increases monotonically with power-law. Also, that if turbulence strength is fixed, the relationship between scintillation index and power-law depends on the operating regime. In weak turbulence the relationship is negative, and it is positive in stronger turbulence. This work also emphasizes the importance of properly scaling turbulence strength with comparing results with different power-laws and the influence and importance of defining inner and outer scales in these simulations.


title: Sampling QCD field configurations with gauge-equivariant flow models

abstract: Machine learning methods based on normalizing flows have been shown to address important challenges, such as critical slowing-down and topological freezing, in the sampling of gauge field configurations in simple lattice field theories. A critical question is whether this success will translate to studies of QCD. This Proceedings presents a status update on advances in this area. In particular, it is illustrated how recently developed algorithmic components may be combined to construct flow-based sampling algorithms for QCD in four dimensions. The prospects and challenges for future use of this approach in at-scale applications are summarized.


title: Bubble nuclei: single-particle versus Coulomb interaction effects

abstract: The detailed investigation of microscopic mechanisms leading to the formation of bubble structures in the nuclei has been performed in the framework of covariant density functional theory. The main emphasis of this study is on the role of single-particle degrees of freedom and Coulomb interaction. In general, the formation of bubbles lowers the Coulomb energy. However, in nuclei this trend is counteracted by the quantum nature of the single-particle states: only specific single-particle states with specific density profiles can be occupied with increasing proton and neutron numbers. A significant role of central classically forbidden region at the bottom of the wine bottle potentials in the formation of nuclear bubbles (via primarily the reduction of the densities of the $s$ states at $r=0$) has been revealed for the first time. Their formation also depends on the availability of low-$l$ single-particle states for occupation since single-particle densities represent the basic building blocks of total densities. Nucleonic potentials disfavor the occupation of such states in hyperheavy nuclei and this contributes to the formation of bubbles in such nuclei. Additivity rule for densities has been proposed for the first time. It was shown that the differences in the densities of bubble and flat density nuclei follow this rule in the $A\approx 40$ mass region and in superheavy nuclei with comparable accuracy. This strongly suggests the same mechanism of the formation of central depression in bubble nuclei of these two mass regions. Nuclear saturation mechanisms and self-consistency effects also affect the formation of bubble structures. The detailed analysis of different aspects of bubble physics strongly suggests that the formation of bubble structures in superheavy nuclei is dominated by single-particle effects.


title: Lattice QCD calculation of $\pi^0\rightarrow e^+ e^-$ decay

abstract: We extend the application of lattice QCD to the two-photon-mediated, order $\alpha^2$ rare decay $\pi^0\rightarrow e^+ e^-$. By combining Minkowski- and Euclidean-space methods we are able to calculate the complex amplitude describing this decay directly from the underlying theories (QCD and QED) which predict this decay. The leading connected and disconnected diagrams are considered; a continuum limit is evaluated and the systematic errors are estimated. We find $\mathrm{Re} \mathcal{A} = 18.60(1.19)(1.04)\,$eV, $\mathrm{Im} \mathcal{A} = 32.59(1.50)(1.65)\,$eV, a more accurate value for the ratio $\frac{\mathrm{Re} \mathcal{A}}{\mathrm{Im} \mathcal{A}}=0.571(10)(4)$ and a result for the partial width $\Gamma(\pi^0\to\gamma\gamma) = 6.60(0.61)(0.67)\,$eV. Here the first errors are statistical and the second systematic. This calculation is the first step in determining the more challenging, two-photon-mediated decay amplitude that contributes to the rare decay $K\to\mu^+\mu^-$.


title: On Transfer of Adversarial Robustness from Pretraining to Downstream Tasks

abstract: As large-scale training regimes have gained popularity, the use of pretrained models for downstream tasks has become common practice in machine learning. While pretraining has been shown to enhance the performance of models in practice, the transfer of robustness properties from pretraining to downstream tasks remains poorly understood. In this study, we demonstrate that the robustness of a linear predictor on downstream tasks can be constrained by the robustness of its underlying representation, regardless of the protocol used for pretraining. We prove (i) a bound on the loss that holds independent of any downstream task, as well as (ii) a criterion for robust classification in particular. We validate our theoretical results in practical applications, show how our results can be used for calibrating expectations of downstream robustness, and when our results are useful for optimal transfer learning. Taken together, our results offer an initial step towards characterizing the requirements of the representation function for reliable post-adaptation performance.


title: Artificial Intelligence and Machine Learning for Quantum Technologies

abstract: In recent years, the dramatic progress in machine learning has begun to impact many areas of science and technology significantly. In the present perspective article, we explore how quantum technologies are benefiting from this revolution. We showcase in illustrative examples how scientists in the past few years have started to use machine learning and more broadly methods of artificial intelligence to analyze quantum measurements, estimate the parameters of quantum devices, discover new quantum experimental setups, protocols, and feedback strategies, and generally improve aspects of quantum computing, quantum communication, and quantum simulation. We highlight open challenges and future possibilities and conclude with some speculative visions for the next decade.


title: Automatic Security Assessment of GitHub Actions Workflows

abstract: The demand for quick and reliable DevOps operations pushed distributors of repository platforms to implement workflows. Workflows allow automating code management operations directly on the repository hosting the software. However, this feature also introduces security issues that directly affect the repository, its content, and all the software supply chains in which the hosted code is involved in. Hence, an attack exploiting vulnerable workflows can affect disruptively large software ecosystems. To empirically assess the importance of this problem, in this paper, we focus on the de-facto main distributor (i.e., GitHub), and we developed a security assessment methodology for GitHub Actions workflows, which are widely adopted in software supply chains. We implemented the methodology in a tool (GHAST) and applied it on 50 open-source projects. The experimental results are worrisome as they allowed identifying a total of 24,905 security issues (all reported to the corresponding stakeholders), thereby indicating that the problem is open and demands further research and investigation.


title: Remarks on totally positive flag manifolds

abstract: We study the intersection of the totally positive part of a split semisimple group over the real numbers with a totally positive parabolic subgroup.


title: Does the Catalog of California Earthquakes, with Aftershocks Included, Contain Information about Future Large Earthquakes?

abstract: Yes. Interval statistics have been used to conclude that major earthquakes are random events in time and cannot be anticipated or predicted. Machine learning is a powerful new technique that enhances our ability to understand the information content of earthquake catalogs. We show that catalogs contain significant information on current hazard and future predictability for large earthquakes.


title: Justification of the use of Stokes-V Cryo-NIRSP/DKIST observations for the 3D Reconstruction of the Coronal Magnetic Field

abstract: This study presents a justification of the use of Stokes-V Cryo-NIRSP/DKIST observations for the 3D Reconstruction of the Coronal Magnetic Field. A magnetohydrodynamic (MHD) model of the solar corona during solar minimum generated by Predictive Science Inc. was used to synthesize spectropolarimetric measurements of the Fe XIII 1075 nm coronal emission line as observed from the Earth. Stokes-Q,U data is "taken" for half a solar rotation period (about two weeks) by Upgraded Coronal Multichannel Polarimeter (UCoMP) with field of view (FOV) up to $2\ R_\odot$. Stokes-V data is "taken" once a day for two days by Cryo-NIRSP/DKIST with a total FOV equal to the two Cryo-NIRSP FOVs. We demonstrated that even this amount Stokes-V observations with Cryo-NIRSP FOV coverage can remarkably improve the 3D coronal magnetic field reconstruction over an active region compared with tomography reconstruction based only on UCoMP linear polarization data.


title: Partial actions on quotient spaces and globalization

abstract: Given a partial action of a topological group $G$ on a space $X$, we determine properties $\mathcal P$ which can be extended from $X$ to its globalization. We treat the cases when $\mathcal P$ is any of the following: Hausdorff, regular, metrizable, second countable and having invariant metric. Further, for a normal subgroup $H$ we introduce and study a partial action of $G/H$ on the orbit space $X/\!\sim,$ applications to invariant metrics and inverse limits are presented.


title: The cold gas and dust properties of red star-forming galaxies

abstract: We study the cold gas and dust properties for a sample of red star forming galaxies called "red misfits." We collect single-dish CO observations and HI observations from representative samples of low-redshift galaxies, as well as our own JCMT CO observations of red misfits. We also obtain SCUBA-2 850 um observations for a subset of these galaxies. With these data we compare the molecular gas, total cold gas, and dust properties of red misfits against those of their blue counterparts ("blue actives") taking non-detections into account using a survival analysis technique. We compare these properties at fixed position in the log SFR-log M* plane, as well as versus offset from the star-forming main sequence. Compared to blue actives, red misfits have slightly longer molecular gas depletion times, similar total gas depletion times, significantly lower molecular- and total-gas mass fractions, lower dust-to-stellar mass ratios, similar dust-to-gas ratios, and a significantly flatter slope in the $\log M_\mathrm{mol}$-$\log M_\star$ plane. Our results suggest that red misfits as a population are likely quenching due to a shortage in gas supply.


title: Counterfactual Fairness Is Basically Demographic Parity

abstract: Making fair decisions is crucial to ethically implementing machine learning algorithms in social settings. In this work, we consider the celebrated definition of counterfactual fairness [Kusner et al., NeurIPS, 2017]. We begin by showing that an algorithm which satisfies counterfactual fairness also satisfies demographic parity, a far simpler fairness constraint. Similarly, we show that all algorithms satisfying demographic parity can be trivially modified to satisfy counterfactual fairness. Together, our results indicate that counterfactual fairness is basically equivalent to demographic parity, which has important implications for the growing body of work on counterfactual fairness. We then validate our theoretical findings empirically, analyzing three existing algorithms for counterfactual fairness against three simple benchmarks. We find that two simple benchmark algorithms outperform all three existing algorithms -- in terms of fairness, accuracy, and efficiency -- on several data sets. Our analysis leads us to formalize a concrete fairness goal: to preserve the order of individuals within protected groups. We believe transparency around the ordering of individuals within protected groups makes fair algorithms more trustworthy. By design, the two simple benchmark algorithms satisfy this goal while the existing algorithms for counterfactual fairness do not.


title: Type-Theoretic Approaches to Ordinals

abstract: In a constructive setting, no concrete formulation of ordinal numbers can simultaneously have all the properties one might be interested in; for example, being able to calculate limits of sequences is constructively incompatible with deciding extensional equality. Using homotopy type theory as the foundational setting, we develop an abstract framework for ordinal theory and establish a collection of desirable properties and constructions. We then study and compare three concrete implementations of ordinals in homotopy type theory: first, a notation system based on Cantor normal forms (binary trees); second, a refined version of Brouwer trees (infinitely-branching trees); and third, extensional well-founded orders. Each of our three formulations has the central properties expected of ordinals, such as being equipped with an extensional and well-founded ordering as well as allowing basic arithmetic operations, but they differ with respect to what they make possible in addition. For example, for finite collections of ordinals, Cantor normal forms have decidable properties, but suprema of infinite collections cannot be computed. In contrast, extensional well-founded orders work well with infinite collections, but almost all properties are undecidable. Brouwer trees take the sweet spot in the middle by combining a restricted form of decidability with the ability to work with infinite increasing sequences. Our three approaches are connected by canonical order-preserving functions from the "more decidable" to the "less decidable" notions. We have formalised the results on Cantor normal forms and Brouwer trees in cubical Agda, while extensional well-founded orders have been studied and formalised thoroughly by Escardo and his collaborators. Finally, we compare the computational efficiency of our implementations with the results reported by Berger.


title: Pion inspired by QCD: Nakanishi and Light-Front Integral Representations

abstract: The pion structure in Minkowski space is explored using the Nakanishi integral representation. A general framework is developed for the pion Bethe-Salpeter amplitude based on the Kallen-Lehmann representation of the dressed quarks with an ansatz for the pseudo-scalar $\pi-q\bar q$ vertex fulfilling the axial Ward-Takahashi identity. The Nakanishi weight functions are derived for the scalar amplitudes associated with the decomposition of the pion Bethe-Salpeter amplitude in different operator bases in the Dirac spinor space in terms of the involved spectral densities. The approach is applied to an analytical model of the pion Bethe-Salpeter amplitude, which is combined with Landau gauge lattice QCD results for the quark running mass at space-like momentum. From the Nakanishi integral representation several pion observables were calculated, such as the decay constant, the spin decomposition of the valence probabilities, longitudinal and transverse momentum distributions from the valence component of the light-front wave function.


title: Error Profile for Discontinuous Galerkin Time Stepping of Parabolic PDEs

abstract: We consider the time discretization of a linear parabolic problem by the discontinuous Galerkin (DG) method using piecewise polynomials of degree at most $r-1$ in $t$, for $r\ge1$ and with maximum step size~$k$. It is well known that the spatial $L_2$-norm of the DG error is of optimal order $k^r$ globally in time, and is, for $r\ge2$, superconvergent of order $k^{2r-1}$ at the nodes. We show that on the $n$th subinterval $(t_{n-1},t_n)$, the dominant term in the DG error is proportional to the local right Radau polynomial of degree $r$. This error profile implies that the DG error is of order $k^{r+1}$ at the right-hand Gauss--Radau quadrature points in each interval. We show that the norm of the jump in the DG solution at the left end point $t_{n-1}$ provides an accurate \emph{a posteriori} estimate for the maximum error over the subinterval $(t_{n-1},t_n)$. Furthermore, a simple post-processing step yields a \emph{continuous} piecewise polynomial of degree $r$ with the optimal global convergence rate of order $k^{r+1}$. We illustrate these results with some numerical experiments.


title: Effect of Geometric Confinement on the Stabilization of Skyrmions

abstract: In this paper, we study the stability of skyrmions (SKYs) caused by the geometric confinement (GC) effect observed in nano-domains in recent experiments, where SKYs appear only inside the boundary and is stable at the low magnetic field region. However, the mechanism of the GC effect is unclear for skyrmions. We numerically find that this effect is not observed in the standard model, in which the Dzyaloshinskii-Moriya interaction (DMI) energy is discretized on lattice edges, while the effect is observed in a discrete model where the DMI energy is defined on lattice volume elements. In the latter model, the DMI energy on the surfaces effectively becomes very small compared with that of the inside. From this observation, we study a model with zero DMI energy on the surfaces parallel to the external magnetic field and find that SKY is significantly stabilized in the sample in the low magnetic field region.


title: Information bottleneck theory of high-dimensional regression: relevancy, efficiency and optimality

abstract: Avoiding overfitting is a central challenge in machine learning, yet many large neural networks readily achieve zero training loss. This puzzling contradiction necessitates new approaches to the study of overfitting. Here we quantify overfitting via residual information, defined as the bits in fitted models that encode noise in training data. Information efficient learning algorithms minimize residual information while maximizing the relevant bits, which are predictive of the unknown generative models. We solve this optimization to obtain the information content of optimal algorithms for a linear regression problem and compare it to that of randomized ridge regression. Our results demonstrate the fundamental trade-off between residual and relevant information and characterize the relative information efficiency of randomized regression with respect to optimal algorithms. Finally, using results from random matrix theory, we reveal the information complexity of learning a linear map in high dimensions and unveil information-theoretic analogs of double and multiple descent phenomena.

categories: cs.IT cond-mat.stat-mech cs.LG math.IT physics.data-an stat.ML


title: RadSegNet: A Reliable Approach to Radar Camera Fusion

abstract: Perception systems for autonomous driving have seen significant advancements in their performance over last few years. However, these systems struggle to show robustness in extreme weather conditions because sensors like lidars and cameras, which are the primary sensors in a sensor suite, see a decline in performance under these conditions. In order to solve this problem, camera-radar fusion systems provide a unique opportunity for all weather reliable high quality perception. Cameras provides rich semantic information while radars can work through occlusions and in all weather conditions. In this work, we show that the state-of-the-art fusion methods perform poorly when camera input is degraded, which essentially results in losing the all-weather reliability they set out to achieve. Contrary to these approaches, we propose a new method, RadSegNet, that uses a new design philosophy of independent information extraction and truly achieves reliability in all conditions, including occlusions and adverse weather. We develop and validate our proposed system on the benchmark Astyx dataset and further verify these results on the RADIATE dataset. When compared to state-of-the-art methods, RadSegNet achieves a 27% improvement on Astyx and 41.46% increase on RADIATE, in average precision score and maintains a significantly better performance in adverse weather conditions


title: Asymptotic behavior and Liouville-type theorems for axisymmetric stationary Navier-Stokes equations outside of an infinite cylinder with a periodic boundary condition

abstract: We study the asymptotic behavior of solutions to the steady Navier-Stokes equations outside of an infinite cylinder in $\mathbb{R}^3$. We assume that the flow is periodic in $x_3$-direction and has no swirl. This problem is closely related with two-dimensional exterior problem. Under a condition on the generalized finite Dirichlet integral, we give a pointwise decay estimate of the vorticity at the spatial infinity. Moreover, we prove a Liouville-type theorem only from the condition of the generalized finite Dirichlet integral.


title: Numerical Evaluation of Mittag-Leffler Functions

abstract: The Mittag-Leffler function is computed via a quadrature approximation of a contour integral representation. We compare results for parabolic and hyperbolic contours, and give special attention to evaluation on the real line. The main point of difference with respect to similar approaches from the literature is the way that poles in the integrand are handled. Rational approximation of the Mittag-Leffler function on the negative real axis is also discussed.


title: Reformulation of the twist-3 gluon Sivers effect toward the application to the heavy quarkonium production

abstract: In this paper, we propose a new calculation method for the twist-3 gluon Sivers effect within the collinear factorization approach. The method called pole calculation has been used to derive the cross section formula for the single transverse-spin asymmetry(SSA) as a standard method. We point out that we encounter a problem when we try to apply the pole calculation to the SSA in the heavy quarkonium production whose hadronization mechanism is described by non-relativistic QCD(NRQCD) framework. We show that the new calculation method solves this problem and successfully reproduces known results derived by the pole calculation. Our new method extends the applicability of the twist-3 calculation technique to heavy quarkonium productions which are ideal observables for the investigation of the gluon Sivers effect.


title: Superlinear stochastic heat equation on $\mathbb{R}^d$

abstract: In this paper, we study the stochastic heat equation (SHE) on $\mathbb{R}^d$ subject to a centered Gaussian noise that is white in time and colored in space. We establish the existence and uniqueness of the random field solution in the presence of locally Lipschitz drift and diffusion coefficients, which can have certain superlinear growth. This is a nontrivial extension of the recent work by Dalang, Khoshnevisan and Zhang (2019), where the one-dimensional SHE on $[0,1]$ subject to space-time white noise has been studied.


title: Towards Fair Conversational Recommender Systems

abstract: Conversational recommender systems have demonstrated great success. They can accurately capture a user's current detailed preference -- through a multi-round interaction cycle -- to effectively guide users to a more personalized recommendation. Alas, conversational recommender systems can be plagued by the adverse effects of bias, much like traditional recommenders. In this work, we argue for increased attention on the presence of and methods for counteracting bias in these emerging systems. As a starting point, we propose three fundamental questions that should be deeply examined to enable fairness in conversational recommender systems.


title: Galactic bar resonances with diffusion: an analytic model with implications for bar-dark matter halo dynamical friction

abstract: The secular evolution of disk galaxies is largely driven by resonances between the orbits of 'particles' (stars or dark matter) and the rotation of non-axisymmetric features (spiral arms or a bar). Such resonances may also explain kinematic and photometric features observed in the Milky Way and external galaxies. In simplified cases, these resonant interactions are well understood: for instance, the dynamics of a test particle trapped near a resonance of a steadily rotating bar is easily analyzed using the angle-action tools pioneered by Binney, Monari and others. However, such treatments do not address the stochasticity and messiness inherent to real galaxies - effects which have, with few exceptions, been previously explored only with complex N-body simulations. In this paper, we propose a simple kinetic equation describing the distribution function of particles near an orbital resonance with a rigidly rotating bar, allowing for diffusion of the particles' slow actions. We solve this equation for various values of the dimensionless diffusion strength $\Delta$, and then apply our theory to the calculation of bar-halo dynamical friction. For $\Delta = 0$ we recover the classic result of Tremaine & Weinberg that friction ultimately vanishes, owing to the phase-mixing of resonant orbits. However, for $\Delta > 0$ we find that diffusion suppresses phase-mixing, leading to a finite torque. Our results suggest that stochasticity - be it physical or numerical - tends to increase bar-halo friction, and that bars in cosmological simulations might experience significant artificial slowdown, even if the numerical two-body relaxation time is much longer than a Hubble time.


title: Irreducible polynomials in quadratic semigroups

abstract: We construct many irreducible polynomials within semigroups generated by sets of the form $S=\{x^2+c_1,\dots,x^2+c_s\}$ under composition.


title: Mechanism of paramagnetic spin Seebeck effect

abstract: We have theoretically investigated the spin Seebeck effect (SSE) in a normal metal (NM)/paramagnetic insulator (PI) bilayer system. Through a linear response approach, we calculated the thermal spin pumping from PI to NM and backflow spin current from NM to PI, where the spin-flip scattering via the interfacial exchange coupling between conduction-electron spin in NM and localized spin in PI is taken into account. We found a finite spin current appears at the interface under the difference in the effective temperatures between spins in NM and PI, and its intensity increases by increasing the density of the localized spin $S$. Our model well reproduces the magnetic-field-induced reduction of the paramagnetic SSE in Pt/Gd$_3$Ga$_5$O$_{12}$ experimentally observed when the Zeeman energy is comparable to the thermal energy, which can be interpreted as the suppression of the interfacial spin-flip scattering. The present finding provides an insight into the mechanism of paramagnetic SSEs and the thermally induced spin-current generation in magnetic materials.


title: A Riemannian Optimization Approach to Clustering Problems

abstract: This paper considers the optimization problem in the form of $\min_{X \in \mathcal{F}_v} f(x) + \lambda \|X\|_1,$ where $f$ is smooth, $\mathcal{F}_v = \{X \in \mathbb{R}^{n \times q} : X^T X = I_q, v \in \mathrm{span}(X)\}$, and $v$ is a given positive vector. The clustering models including but not limited to the models used by $k$-means, community detection, and normalized cut can be reformulated as such optimization problems. It is proven that the domain $\mathcal{F}_v$ forms a compact embedded submanifold of $\mathbb{R}^{n \times q}$ and optimization-related tools including a family of computationally efficient retractions and an orthonormal basis of any normal space of $\mathcal{F}_v$ are derived. An inexact accelerated Riemannian proximal gradient method that allows adaptive step size is proposed and its global convergence is established. Numerical experiments on community detection in networks and normalized cut for image segmentation are used to demonstrate the performance of the proposed method.


title: Dihedral rigidity in hyperbolic 3-space

abstract: We prove a comparison theorem for certain types of polyhedra in a 3-manifold with its scalar curvature bounded below by $-6$. The result confirms in some cases the Gromov dihedral rigidity conjecture in hyperbolic $3$-space.


title: Rankability and Linear Ordering Problem: New Probabilistic Insight and Algorithms

abstract: The linear ordering problem (LOP), which consists in ordering M objects from their pairwise comparisons, is commonly applied in many areas of research. While efforts have been made to devise efficient LOP algorithms, verification of whether the data are rankable, that is, if the linear ordering problem (LOP) solutions have a meaningful interpretation, received much less attention. To address this problem, we adopt a probabilistic perspective where the results of pairwise comparisons are modeled as Bernoulli variables with a common parameter and we estimate the latter from the observed data. The brute-force approach to the required enumeration has a prohibitive complexity of O(M !), so we reformulate the problem and introduce a concept of the Slater spectrum that generalizes the Slater index, and then devise an algorithm to find the spectrum with complexity O(M^3 2^M) that is manageable for moderate values of M. Furthermore, with a minor modification of the algorithm, we are able to find all solutions of the LOP with the complexity O(M 2^M). Numerical examples are shown on synthetic and real-world data, and the algorithms are publicly available.


title: Learning and predicting photonic responses of plasmonic nanoparticle assemblies via dual variational autoencoders

abstract: We demonstrate the application of machine learning for rapid and accurate extraction of plasmonic particles cluster geometries from hyperspectral image data via a dual variational autoencoder (dual-VAE). In this approach, the information is shared between the latent spaces of two VAEs acting on the particle shape data and spectral data, respectively, but enforcing a common encoding on the shape-spectra pairs. We show that this approach can establish the relationship between the geometric characteristics of nanoparticles and their far-field photonic responses, demonstrating that we can use hyperspectral darkfield microscopy to accurately predict the geometry (number of particles, arrangement) of a multiparticle assemblies below the diffraction limit in an automated fashion with high fidelity (for monomers (0.96), dimers (0.86), and trimers (0.58). This approach of building structure-property relationships via shared encoding is universal and should have applications to a broader range of materials science and physics problems in imaging of both molecular and nanomaterial systems.


title: DeepTLS: comprehensive and high-performance feature extraction for encrypted traffic

abstract: Feature extraction is critical for TLS traffic analysis using machine learning techniques, which it is also very difficult and time-consuming requiring huge engineering efforts. We designed and implemented DeepTLS, a system which extracts full spectrum of features from pcaps across meta, statistical, SPLT, byte distribution, TLS header and certificates. The backend is written in C++ to achieve high performance, which can analyze a GB-size pcap in a few minutes. DeepTLS was thoroughly evaluated against two state-of-the-art tools Joy and Zeek with four well-known malicious traffic datasets consisted of 160 pcaps. Evaluation results show DeepTLS has advantage of analyzing large pcaps with half analysis time, and identified more certificates with acceptable performance loss compared with Joy. DeepTLS can significantly accelerate machine learning pipeline by reducing feature extraction time from hours even days to minutes. The system is online at https://deeptls.com, where test artifacts can be viewed and validated. In addition, two open source tools Pysharkfeat and Tlsfeatmark are also released.


title: Creating Reverse Bilingual Dictionaries

abstract: Bilingual dictionaries are expensive resources and not many are available when one of the languages is resource-poor. In this paper, we propose algorithms for creation of new reverse bilingual dictionaries from existing bilingual dictionaries in which English is one of the two languages. Our algorithms exploit the similarity between word-concept pairs using the English Wordnet to produce reverse dictionary entries. Since our algorithms rely on available bilingual dictionaries, they are applicable to any bilingual dictionary as long as one of the two languages has Wordnet type lexical ontology.


title: Minimal Binary Linear Codes from Vectorial Boolean Functions

abstract: Recently, much progress has been made to construct minimal linear codes due to their preference in secret sharing schemes and secure two-party computation. In this paper, we put forward a new method to construct minimal linear codes by using vectorial Boolean functions. Firstly, we give a necessary and sufficient condition for a generic class of linear codes from vectorial Boolean functions to be minimal. Based on that, we derive some new three-weight minimal linear codes and determine their weight distributions. Secondly, we obtain a necessary and sufficient condition for another generic class of linear codes from vectorial Boolean functions to be minimal and to be violated the AB condition. As a result, we get three infinite families of minimal linear codes violating the AB condition. To the best of our knowledge, this is the first time that minimal liner codes are constructed from vectorial Boolean functions. Compared with other known ones, in general the minimal liner codes obtained in this paper have higher dimensions.


title: A Python-based tool for constructing observables from the DSN's closed-loop archival tracking data files

abstract: Radio science data collected from NASA's Deep Space Networks (DSNs) are made available in various formats through NASA's Planetary Data System (PDS). The majority of these data are packed in complex formats, making them inaccessible to users without specialized knowledge. In this paper, we present a Python-based tool that can preprocess the closed-loop archival tracking data files (ATDFs), produce Doppler and range observables, and write them in an ASCII table along with ancillary information. ATDFs are primitive closed-loop radio science products with limited available documentation. Early in the 2000s, DSN deprecated ATDF and replaced it with the Tracking and Navigation Service Data Files (TNF) to keep up with the evolution of the radio science system. Most data processing software (e.g., orbit determination software) cannot use them directly, thus limiting the utilization of these data. As such, the vast majority of historical closed-loop radio science data have not yet been processed with modern software and with our improved understanding of the solar system. The preprocessing tool presented in this paper makes it possible to revisit such historical data using modern techniques and software to conduct crucial radio science experiments.


title: Valid Widgets Contain Legal Subwidgets

abstract: This paper proves a linear algebra result that has to do with the geometry of "widgets". For us a widget is a collection of n pairs of points in a vector space. (The pairs represent the different possible spin states of a particle.) We investigate linear relations among such collections. A corollary of our theorem was conjectured in arXiv:2208.02478v1 where it arose in an attempt to understand some issues in super string theory. In that paper an investigation of perturbative superstring theory with Ramond punctures required the special case when the ambient dimension is n. Here we prove the general case.


title: Elastic nucleon-pion scattering at $m_{\pi} = 200~{\rm MeV}$ from lattice QCD

abstract: Elastic nucleon-pion scattering amplitudes are computed using lattice QCD on a single ensemble of gauge field configurations with $N_{\rm f} = 2+1$ dynamical quark flavors and $m_{\pi} = 200~{\rm MeV}$. The $s$-wave scattering lengths with both total isospins $I=1/2$ and $I=3/2$ are inferred from the finite-volume spectrum below the inelastic threshold together with the $I=3/2$ $p$-wave containing the $\Delta(1232)$ resonance. The amplitudes are well-described by the effective range expansion with parameters constrained by fits to the finite-volume energy levels enabling a determination of the $I=3/2$ scattering length with statistical errors below $5\%$, while the $I = 1/2$ is somewhat less precise. Systematic errors due to excited states and the influence of higher partial waves are controlled, providing a pathway for future computations down to the physical light quark masses with multiple lattice spacings and physical volumes.


title: Data-centric AI approach to improve optic nerve head segmentation and localization in OCT en face images

abstract: The automatic detection and localization of anatomical features in retinal imaging data are relevant for many aspects. In this work, we follow a data-centric approach to optimize classifier training for optic nerve head detection and localization in optical coherence tomography en face images of the retina. We examine the effect of domain knowledge driven spatial complexity reduction on the resulting optic nerve head segmentation and localization performance. We present a machine learning approach for segmenting optic nerve head in 2D en face projections of 3D widefield swept source optical coherence tomography scans that enables the automated assessment of large amounts of data. Evaluation on manually annotated 2D en face images of the retina demonstrates that training of a standard U-Net can yield improved optic nerve head segmentation and localization performance when the underlying pixel-level binary classification task is spatially relaxed through domain knowledge.


title: Animated Vega-Lite: Unifying Animation with a Grammar of Interactive Graphics

abstract: We present Animated Vega-Lite, a set of extensions to Vega-Lite that model animated visualizations as time-varying data queries. In contrast to alternate approaches for specifying animated visualizations, which prize a highly expressive design space, Animated Vega-Lite prioritizes unifying animation with the language's existing abstractions for static and interactive visualizations to enable authors to smoothly move between or combine these modalities. Thus, to compose animation with static visualizations, we represent time as an encoding channel. Time encodings map a data field to animation keyframes, providing a lightweight specification for animations without interaction. To compose animation and interaction, we also represent time as an event stream; Vega-Lite selections, which provide dynamic data queries, are now driven not only by input events but by timer ticks as well. We evaluate the expressiveness of our approach through a gallery of diverse examples that demonstrate coverage over taxonomies of both interaction and animation. We also critically reflect on the conceptual affordances and limitations of our contribution by interviewing five expert developers of existing animation grammars. These reflections highlight the key motivating role of in-the-wild examples, and identify three central tradeoffs: the language design process, the types of animated transitions supported, and how the systems model keyframes.


title: Automatically constructing Wordnet synsets

abstract: Manually constructing a Wordnet is a difficult task, needing years of experts' time. As a first step to automatically construct full Wordnets, we propose approaches to generate Wordnet synsets for languages both resource-rich and resource-poor, using publicly available Wordnets, a machine translator and/or a single bilingual dictionary. Our algorithms translate synsets of existing Wordnets to a target language T, then apply a ranking method on the translation candidates to find best translations in T. Our approaches are applicable to any language which has at least one existing bilingual dictionary translating from English to it.


title: Optimal Dynamic Orchestration in NDN-based Computing Networks

abstract: Named Data Networking (NDN) offers promising advantages in deploying next-generation service applications over distributed computing networks. We consider the problem of dynamic orchestration over a NDN-based computing network, in which nodes can be equipped with communication, computation, and data producing resources. Given a set of services with function-chaining structures, we address the design of distributed online algorithm that controls each node to make adaptive decisions on flowing service requests, committing function implementations, and/or producing data. We design a Service Discovery Assisted Dynamic Orchestration (SDADO) algorithm that reduces the end-to-end (E2E) delay of delivering the services, while providing optimal throughput performance. The proposed algorithm hybrids queuing-based flexibility and topology-based discipline, where the topological information is not pre-available but obtained through our proposed service discovery mechanism. We provide throughput-optimality analysis for SDADO, and then provide numerical results that confirm our analysis and demonstrates reduced round-trip E2E delay.


title: SHAARP: An Open-Source Package for Analytical and Numerical Modeling of Optical Second Harmonic Generation in Anisotropic Crystals

abstract: Optical second harmonic generation is a second-order nonlinear process that combines two photons of a given frequency into a third photon at twice the frequency. Due to the symmetry constraints, it is widely used as a sensitive probe to detect broken inversion symmetry and local polar order. Analytical modeling of the electric-dipole SHG response is essential to extract fundamental properties of materials from experiments. However, complexity builds up dramatically in the analytical model when the probed crystal is of a low bulk crystal symmetry, with a low-symmetry surface orientation, exhibits absorption and dispersion, and consists of multiple interfaces. As a result, there is a largely uneven landscape in the literature on the SHG modeling of new materials, involving numerous approximations and a wide range of (in)accuracies, leading to a rather scattered dataset of reported SHG nonlinear susceptibility. Towards streamlining the reliability and accuracy of this process, we have developed an open-source package called the Second Harmonic Analysis of Anisotropic Rotational Polarimetry (SHAARP) which derives analytical solutions and performs numerical simulations of reflection SHG from a single interface for homogeneous crystals. Five key generalizations in SHG modeling are implemented, including all crystal symmetries down to triclinic, any crystal orientation, complex dielectric tensor (refractive indices) with frequency dispersion, and general polarization states of the light. SHAARP enables accurate anisotropic modeling of SHG response for a broad range of materials systems. The method is extendible to multiple interfaces. The code is free to download from https://github.com/Rui-Zu/SHAARP


title: CheXRelNet: An Anatomy-Aware Model for Tracking Longitudinal Relationships between Chest X-Rays

abstract: Despite the progress in utilizing deep learning to automate chest radiograph interpretation and disease diagnosis tasks, change between sequential Chest X-rays (CXRs) has received limited attention. Monitoring the progression of pathologies that are visualized through chest imaging poses several challenges in anatomical motion estimation and image registration, i.e., spatially aligning the two images and modeling temporal dynamics in change detection. In this work, we propose CheXRelNet, a neural model that can track longitudinal pathology change relations between two CXRs. CheXRelNet incorporates local and global visual features, utilizes inter-image and intra-image anatomical information, and learns dependencies between anatomical region attributes, to accurately predict disease change for a pair of CXRs. Experimental results on the Chest ImaGenome dataset show increased downstream performance compared to baselines. Code is available at https://github.com/PLAN-Lab/ChexRelNet


title: Scaling limits and fluctuations of a family of $N$-urn branching processes

abstract: In this paper we are concerned with a family of $N$-urn branching processes, where some particles are put into $N$ urns initially and then each particle gives birth to several new particles in some urn when dies. This model includes the $N$-urn Ehrenfest model and the $N$-urn branching random walk as special cases. We show that the scaling limit of the process is driven by a $C(\mathbb{T})$-valued linear ordinary differential equation and the fluctuation of the process is driven by a generalized Ornstein-Uhlenbeck process in the dual of $C^\infty(\mathbb{T})$, where $\mathbb{T}=(0, 1]$ is the one-dimensional torus. A crucial step for proofs of above main results is to show that numbers of particles in different urns are approximately independent. As applications of our main results, limit theorems of hitting times of the process are also discussed.


title: Explicit K-symplectic methods for nonseparable non-canonical Hamiltonian systems

abstract: We propose efficient numerical methods for nonseparable non-canonical Hamiltonian systems which are explicit, K-symplectic in the extended phase space with long time energy conservation properties. They are based on extending the original phase space to several copies of the phase space and imposing a mechanical restraint on the copies of the phase space. Explicit K-symplectic methods are constructed for three non-canonical Hamiltonian systems. Numerical results show that they outperform the higher order Runge-Kutta methods in preserving the phase orbit and the energy of the system over long time.


title: Creating Lexical Resources for Endangered Languages

abstract: This paper examines approaches to generate lexical resources for endangered languages. Our algorithms construct bilingual dictionaries and multilingual thesauruses using public Wordnets and a machine translator (MT). Since our work relies on only one bilingual dictionary between an endangered language and an "intermediate helper" language, it is applicable to languages that lack many existing resources.


title: Learning Entity Linking Features for Emerging Entities

abstract: Entity linking (EL) is the process of linking entity mentions appearing in text with their corresponding entities in a knowledge base. EL features of entities (e.g., prior probability, relatedness score, and entity embedding) are usually estimated based on Wikipedia. However, for newly emerging entities (EEs) which have just been discovered in news, they may still not be included in Wikipedia yet. As a consequence, it is unable to obtain required EL features for those EEs from Wikipedia and EL models will always fail to link ambiguous mentions with those EEs correctly as the absence of their EL features. To deal with this problem, in this paper we focus on a new task of learning EL features for emerging entities in a general way. We propose a novel approach called STAMO to learn high-quality EL features for EEs automatically, which needs just a small number of labeled documents for each EE collected from the Web, as it could further leverage the knowledge hidden in the unlabeled data. STAMO is mainly based on self-training, which makes it flexibly integrated with any EL feature or EL model, but also makes it easily suffer from the error reinforcement problem caused by the mislabeled data. Instead of some common self-training strategies that try to throw the mislabeled data away explicitly, we regard self-training as a multiple optimization process with respect to the EL features of EEs, and propose both intra-slot and inter-slot optimizations to alleviate the error reinforcement problem implicitly. We construct two EL datasets involving selected EEs to evaluate the quality of obtained EL features for EEs, and the experimental results show that our approach significantly outperforms other baseline methods of learning EL features.


title: On the long-time asymptotic of the modified Camassa-Holm equation with nonzero boundary conditions in space-time solitonic regions

abstract: We investigate the long-time asymptotic behavior for the Cauchy problem of the modified Camassa-Holm (mCH) equation with nonzero boundary conditions in different regions \begin{align*} &m_{t}+\left((u^2-u_x^2)m\right)_{x}=0,~~ m=u-u_{xx}, ~~ (x,t)\in\mathbb{R}\times\mathbb{R}^{+},\\ &u(x,0)=u_{0}(x),~~\lim_{x\to\pm\infty} u_{0}(x)=1,~~u_{0}(x)-1\in H^{4,1}(\mathbb{R}), \end{align*} where $m(x,t=0):=m_{0}(x)$ and $m_{0}(x)-1\in H^{2,1}(\mathbb{R})$. Through spectral analysis, the initial value problem of the mCH equation is transformed into a matrix RH problem on a new plane $(y,t)$, and then using the $\overline{\partial}$-nonlinear steepest descent method, we analyze the different asymptotic behaviors of the four regions divided by the interval of $\xi=y/t$ on plane $\{(y,t)|y\in(-\infty,+\infty), t>0\}$. There is no steady-state phase point corresponding to the regions $\xi\in(-\infty,-1/4)\cup(2,\infty)$. We prove that the solution of mCH equation is characterized by $N$-soliton solution and error on these two regions. In $\xi\in(-1/4,0)$ and $\xi\in(0,2)$, the phase function $\theta(z)$ has eight and four steady-state phase points, respectively. We prove that the soliton resolution conjecture holds, that is, the solution of the mCH equation can be expressed as the soliton solution on the discrete spectrum, the leading term on the continuous spectrum, and the residual error. Our results also show that soliton solutions of the mCH equation with nonzero condition boundary are asymptotically stable.


title: Clear Memory-Augmented Auto-Encoder for Surface Defect Detection

abstract: In surface defect detection, due to the extreme imbalance in the number of positive and negative samples, positive-samples-based anomaly detection methods have received more and more attention. Specifically, reconstruction-based methods are the most popular. However, existing methods are either difficult to repair abnormal foregrounds or reconstruct clear backgrounds. Therefore, we propose a clear memory-augmented auto-encoder (CMA-AE). At first, we propose a novel clear memory-augmented module (CMAM), which combines the encoding and memoryencoding in a way of forgetting and inputting, thereby repairing abnormal foregrounds and preserving clear backgrounds. Secondly, a general artificial anomaly generation algorithm (GAAGA) is proposed to simulate anomalies that are as realistic and feature-rich as possible. At last, we propose a novel multi scale feature residual detection method (MSFR) for defect segmentation, which makes the defect location more accurate. Extensive comparison experiments demonstrate that CMA-AE achieves state-of-the-art detection accuracy and shows great potential in industrial applications.


title: BASS XXXII: Studying the Nuclear Mm-wave Continuum Emission of AGNs with ALMA at Scales $\lesssim$ 100-200 pc

abstract: To understand the origin of nuclear ($\lesssim$ 100 pc) millimeter-wave (mm-wave) continuum emission in active galactic nuclei (AGNs), we systematically analyzed sub-arcsec resolution Band-6 (211-275 GHz) ALMA data of 98 nearby AGNs ($z <$ 0.05) from the 70-month Swift/BAT catalog. The sample, almost unbiased for obscured systems, provides the largest number of AGNs to date with high mm-wave spatial resolution sampling ($\sim$ 1-200 pc), and spans broad ranges of 14-150 keV luminosity {$40 < \log[L_{\rm 14-150}/({\rm erg\,s^{-1}})] < 45$}, black hole mass [$5 < \log(M_{\rm BH}/M_\odot) < 10$], and Eddington ratio ($-4 < \log \lambda_{\rm Edd} < 2$). We find a significant correlation between 1.3 mm (230 GHz) and 14-150 keV luminosities. Its scatter is $\approx$ 0.36 dex, and the mm-wave emission may serve as a good proxy of the AGN luminosity, free of dust extinction up to $N_{\rm H} \sim 10^{26}$ cm$^{-2}$. While the mm-wave emission could be self-absorbed synchrotron radiation around the X-ray corona according to past works, we also discuss different possible origins of the mm-wave emission; AGN-related dust emission, outflow-driven shocks, and a small-scale ($<$ 200 pc) jet. The dust emission is unlikely to be dominant, as the mm-wave slope is generally flatter than expected. Also, due to no increase in the mm-wave luminosity with the Eddington ratio, a radiation-driven outflow model is possibly not the common mechanism. Furthermore, we find independence of the mm-wave luminosity on indicators of the inclination angle from the polar axis of the nuclear structure, which is inconsistent with a jet model whose luminosity depends only on the angle.


title: Network Critical Slowing Down: Data-Driven Detection of Critical Transitions in Nonlinear Networks

abstract: In a Nature article, Scheffer et al. presented a novel data-driven framework to predict critical transitions in complex systems. These transitions, which may stem from failures, degradation, or adversarial actions, have been attributed to bifurcations in the nonlinear dynamics. Their approach was built upon the phenomenon of critical slowing down, i.e., slow recovery in response to small perturbations near bifurcations. We extend their approach to detect and localize critical transitions in nonlinear networks. By introducing the notion of network critical slowing down, the objective of this paper is to detect that the network is undergoing a bifurcation only by analyzing its signatures from measurement data. We focus on two classes of widely-used nonlinear networks: (1) Kuramoto model for the synchronization of coupled oscillators and (2) attraction-repulsion dynamics in swarms, each of which presents a specific type of bifurcation. Based on the phenomenon of critical slowing down, we study the asymptotic behavior of the perturbed system away and close to the bifurcation and leverage this fact to develop a deterministic method to detect and identify critical transitions in nonlinear networks. Furthermore, we study the state covariance matrix subject to a stochastic noise process away and close to the bifurcation and use it to develop a stochastic framework for detecting critical transitions. Our simulation results show the strengths and limitations of the methods.


title: On a generalization of Busemann's intersection inequality

abstract: Busemann's intersection inequality gives an upper bound for the volume of the intersection body of a star body in terms of the volume of the body itself. Koldobsky, Paouris, and Zymonopoulou asked if there is a similar result for $k$-intersection bodies. We solve this problem for star bodies that are close to the Euclidean ball in the Banach-Mazur distance. We also improve a bound obtained by Koldobsky, Paouris, and Zymonopoulou for general star bodies in the case when $k$ is proportional to the dimension.


title: The double thresholds distort the lineshapes of the $P_{\psi s}^\Lambda(4338)^0$ resonance

abstract: Very recently, the LHCb Collaboration reported the first observation of the hidden charm pentaquark with strangeness, $P_{\psi s}^\Lambda(4338)^0$. Considering this state is very close to the $\Xi_c^0\bar{D}^0$ and $\Xi_c^+D^-$ thresholds, we explore the possible bias of the Breit-Wigner parameterization, with emphasis on the effect of its coupling to the double thresholds $\Xi_c^0\bar{D}^0$ and $\Xi_c^+D^-$. We first use a qualitative picture based on the ``uniformization" of the Riemann surface of the two-channel system to understand the positions of the enhancement. Then we use the Lippmann-Schwinger equation (LSE) formalism (equivalent to the $K$-matrix parameterization) with two models, the zero-range model and Flatt\'e model to investigate the $J/\psi \Lambda$ lineshapes. Our results show that the nominal peak of the $P_{\psi s}^\Lambda(4338)^0$ could arise either from the pole well above the $\Xi_c^+ D^-$ threshold on the $(-,+)$ sheet or from the pole well below the $\Xi_c^0 \bar{D}^0$ threshold on the $(-,-)$ sheet in the two-channel system. Using the Breit-Wigner distribution to depict the above two lineshapes could be misleading. We also find a novel type of lineshapes with the enhancement constrained by the threshold difference. We urge the LHCb collaboration to perform the refined experimental analysis considering the unitarity and analyticity, e.g. using the K-matrix parameterization. As a by-product, we obtain that the ratio of the isospin violating decay ${\Gamma_{P_{\psi s}^{\Lambda}\to J/\psi\Sigma}}/{\Gamma_{P_{\psi s}^{\Lambda}\to J/\psi\Lambda}}$ could be up to $10\%$.


title: Differential biases, $c$-differential uniformity, and their relation to differential attacks

abstract: Differential cryptanalysis famously uses statistical biases in the propagation of differences in a block cipher to attack the cipher. In this paper, we investigate the existence of more general statistical biases in the differences. To this end, we discuss the $c$-differential uniformity of S-boxes, which is a concept that was recently introduced in Ellingsen et. al. to measure certain statistical biases that could potentially be used in attacks similar to differential attacks. Firstly, we prove that a large class of potential candidates for S-boxes necessarily has large $c$-differential uniformity for all but at most $B$ choices of $c$, where $B$ is a constant independent of the size of the finite field $q$. This result implies that for a large class of functions, certain statistical differential biases are inevitable. In a second part, we discuss the practical possibility of designing a differential attack based on weaknesses of S-boxes related to their $c$-differential uniformity.


title: Statistical Properties of the Probabilistic Numeric Linear Solver BayesCG

abstract: We analyse the calibration of BayesCG under the Krylov prior, a probabilistic numeric extension of the Conjugate Gradient (CG) method for solving systems of linear equations with symmetric positive definite coefficient matrix. Calibration refers to the statistical quality of the posterior covariances produced by a solver. Since BayesCG is not calibrated in the strict existing notion, we propose instead two test statistics that are necessary but not sufficient for calibration: the Z-statistic and the new S-statistic. We show analytically and experimentally that under low-rank approximate Krylov posteriors, BayesCG exhibits desirable properties of a calibrated solver, is only slightly optimistic, and is computationally competitive with CG.


title: What can we know about that which we cannot even imagine?

abstract: In this essay I will consider a sequence of questions. The first questions concern the biological function of intelligence in general, and cognitive prostheses of human intelligence in particular. These will lead into questions concerning human language, perhaps the most important cognitive prosthesis humanity has ever developed. While it is traditional to rhapsodize about the cognitive power encapsulated in human language, I will emphasize how horribly limited human language is - and therefore how limited our cognitive abilities are, despite their being augmented with language. This will lead to questions of whether human mathematics, being ultimately formulated in terms of human language, is also deeply limited. I will then combine these questions to pose a partial, sort-of, sideways answer to the guiding concern of this essay: what we can ever discern about that we cannot even conceive?


title: Sensitivity Analyses of Clinical Trial Designs: Selecting Scenarios and Summarizing Operating Characteristics

abstract: The use of simulation-based sensitivity analyses is fundamental to evaluate and compare candidate designs for future clinical trials. In this context, sensitivity analyses are especially useful to assess the dependence of important design operating characteristics (OCs) with respect to various unknown parameters (UPs). Typical examples of OCs include the likelihood of detecting treatment effects and the average study duration, which depend on UPs that are not known until after the onset of the clinical study, such as the distributions of the primary outcomes and patient profiles. Two crucial components of sensitivity analyses are (i) the choice of a set of plausible simulation scenarios $\{\boldsymbol{\theta}_1,...,\boldsymbol{\theta}_K\}$ and (ii) the list of OCs of interest. We propose a new approach to choose the set of scenarios for inclusion in design sensitivity analyses. Our approach balances the need for simplicity and interpretability of OCs computed across several scenarios with the need to faithfully summarize -- through simulations -- how the OCs vary across all plausible values of the UPs. Our proposal also supports the selection of the number of simulation scenarios to be included in the final sensitivity analysis report. To achieve these goals, we minimize a loss function $\mathcal{L}(\boldsymbol{\theta}_1,...,\boldsymbol{\theta}_K)$ that formalizes whether a specific set of $K$ sensitivity scenarios $\{\boldsymbol{\theta}_1,...,\boldsymbol{\theta}_K\}$ is adequate to summarize how the OCs of the trial design vary across all plausible values of the UPs. Then, we use optimization techniques to select the best set of simulation scenarios to exemplify the OCs of the trial design.


title: Spin Hall conductivity of interacting two-dimensional electron systems

abstract: We consider a two-dimensional electron system subjected to a short-ranged nonmagnetic disorder potential, Coulomb interactions, and Rashba spin-orbit coupling. The path-integral approach incorporated within the Keldysh formalism is used to derive the kinetic equation for the semiclassical Green's function and applied to compute the spin current within the linear response theory. We discuss the frequency dependence of the spin Hall conductivity and further elucidate the role of electron interactions at finite temperatures for both the ballistic and diffusive regimes of transport. We argue that interaction corrections to the spin Hall effect stem from the quantum interference processes whose magnitude is estimated in terms of parameters of the considered model.


title: Robust Training and Verification of Implicit Neural Networks: A Non-Euclidean Contractive Approach

abstract: This paper proposes a theoretical and computational framework for training and robustness verification of implicit neural networks based upon non-Euclidean contraction theory. The basic idea is to cast the robustness analysis of a neural network as a reachability problem and use (i) the $\ell_{\infty}$-norm input-output Lipschitz constant and (ii) the tight inclusion function of the network to over-approximate its reachable sets. First, for a given implicit neural network, we use $\ell_{\infty}$-matrix measures to propose sufficient conditions for its well-posedness, design an iterative algorithm to compute its fixed points, and provide upper bounds for its $\ell_\infty$-norm input-output Lipschitz constant. Second, we introduce a related embedded network and show that the embedded network can be used to provide an $\ell_\infty$-norm box over-approximation of the reachable sets of the original network. Moreover, we use the embedded network to design an iterative algorithm for computing the upper bounds of the original system's tight inclusion function. Third, we use the upper bounds of the Lipschitz constants and the upper bounds of the tight inclusion functions to design two algorithms for the training and robustness verification of implicit neural networks. Finally, we apply our algorithms to train implicit neural networks on the MNIST dataset and compare the robustness of our models with the models trained via existing approaches in the literature.


title: Fourth-order conservative non-splitting semi-Lagrangian Hermite WENO schemes for kinetic and fluid simulations

abstract: We present fourth-order conservative non-splitting semi-Lagrangian (SL) Hermite essentially non-oscillatory (HWENO) schemes for linear transport equations with applications for nonlinear problems including the Vlasov-Poisson system, the guiding center Vlasov model, and the incompressible Euler equations in the vorticity-stream function formulation. The proposed SL HWENO schemes combine a weak formulation of the characteristic Galerkin method with two newly constructed HWENO reconstruction methods. Fourth-order accuracy is accomplished in both space and time under a non-splitting setting. Mass conservation naturally holds due to the weak formulation of the characteristic Galerkin method and the design of the HWENO reconstructions. We apply a positive-preserving limiter to maintain the positivity of numerical solutions when needed. Although the proposed SL framework allows us to take large time steps for improving computational efficiency, it also brings challenges to the spatial reconstruction technique; we construct two kind of novel HWENO reconstructions to fit the need for the proposed SL framework. Abundant benchmark tests are performed to verify the effectiveness of the proposed SL HWENO schemes.


title: The chemical abundance pattern of the extremely metal-poor thin disk star 2MASS J1808-5104 and its origins

abstract: We present a high-resolution ($R\sim35,000$), high signal-to-noise ($S/N=350$) Magellan/MIKE spectrum of the bright extremely metal-poor star 2MASS~J1808$-$5104. We find [Fe/H] = $-$4.01 (spectroscopic LTE stellar parameters), [Fe/H] = $-$3.8 (photometric stellar parameters), [Fe/H] = $-$3.7 (spectroscopic NLTE stellar parameters). We measured a carbon-to-iron ratio of $\mbox{[C/Fe]}= 0.38$ from the CH G-band. J1808$-$5104 is thus not carbon-enhanced, contrary to many other stars with similarly low iron abundances. We also determine, for the first time, a barium abundance ($\mbox{[Ba/Fe]} =-0.78$), and obtain a significantly reduced upper limit for the nitrogen abundance ([N/Fe]$ < - 0.2$). J1808$-$5104 has low ratio of $\mbox{[Sr/Ba]}=-0.17$, which is consistent with that of stars in ultra-faint dwarf galaxies. We also fit the abundance pattern of J1808$-$5104 with nucleosynthesis yields from a grid of Population\,III supernova models. There is a good fit to the abundance pattern which suggests J1808$-$5104 originated from gas enriched by a single massive supernova with a high explosion energy of E $=10\times10^{51}$\,erg and a progenitor stellar mass of M$=29.5$\,M$_{\odot}$. Interestingly, J1808$-$5104 is a member of the Galactic thin disk, as confirmed by our detailed kinematic analysis and calculated stellar actions and velocities. Finally, we also established the orbital history of J1808$-$5104 using our time-dependent Galactic potential the \texttt{ORIENT}. J1808$-$5104 appears to have a stable quasi-circular orbit and been largely confined to the thin disk. This unique orbital history, the star's very old age ($\sim13.5$\,Gyr), and the low [C/Fe] and [Sr/Ba] ratios suggest that J1808$-$5104 may have formed at the earliest epoch of the hierarchical assembly of the Milky Way, and it is most likely associated with the primordial thin disk.


title: Composition-differentiation operators on $S^2(\mathbb{D})$

abstract: We investigate composition-differentiation operators acting on the space $S^2$, the space of analytic functions on the open unit disk whose first derivative is in $H^2$. Specifically, we determine characterizations for bounded and compact composition-differentiation operators acting on $S^p$. In addition, for particular classes of inducing maps, we compute the norm, and identify the spectrum. Finally, for particular linear fractional inducing maps, we determine the adjoint of the composition-differentiation operator acting on weighted Bergman spaces which include $S^2, H^2$, and the Dirichlet space.


title: Category $\mathcal{O}$ for the Lie algebra of vector fields on the line

abstract: Let $\mathfrak{W}$ be the Lie algebra of vector fields on the line. Via computing extensions between all simple modules in the category $\mathcal{O}$, we give the block decomposition of $\mathcal{O}$, and show that the representation type of each block of $\mathcal{O}$ is wild using the Ext-quiver. Each block of $\mathcal{O}$ has infinite simple objects. This result is very different from that of $\mathcal{O}$ for complex semisimple Lie algebras. To find a connection between $\mathcal{O}$ and the module category over some associative algebra, we define a subalgebra $H_1$ of $U(\mathfrak{b})$. We give an exact functor from $\mathcal{O}$ to the category $\Omega$ of finite dimensional modules over $H_1$. We also construct new simple $\mathfrak{W}$-modules from Weyl modules and modules over the Borel subalgebra $\mathfrak{b}$ of $\mathfrak{W}$.


title: Experimental study of secure quantum key distribution with source and detection imperfections

abstract: The quantum key distribution (QKD), guaranteed by the principle of quantum physics, is a promising solution for future secure information and communication technology. However, device imperfections compromise the security of real-life QKD systems, restricting the wide deployment of QKD. This study reports a decoy-state BB84 QKD experiment that considers both source and detection imperfections. In particular, we achieved a rigorous finite-key security bound over fiber links of up to 75 km by applying a systematic performance analysis. Furthermore, our study considers more device imperfections than most previous experiments, and the proposed theory can be extended to other discrete-variable QKD systems. These features constitute a crucial step toward securing QKD with imperfect practical devices.