title: Active Disturbance Rejection Control (ADRC) Toolbox for MATLAB/Simulink

abstract: In this study, an active disturbance rejection control (ADRC) toolbox for MATLAB/Simulink is introduced. Although ADRC has already been established as a powerful robust control framework with successful industrial implementations and strong theoretical foundations, a comprehensive tool for computer-aided design of ADRC has not been developed until now. The proposed open-source ADRC Toolbox is a response to the growing need in the scientific community and the control industry for a straightforward software application of the ADRC methodology. Its main purpose is to fill the gap between the current theories and applications of ADRC and to provide an easy-to-use solution for users in various control fields who want to employ the ADRC scheme in their applications. The ADRC Toolbox contains a single, general-purpose, drag-and-drop function block that allows the synthesis of a predefined ADRC-based strategy with minimal design effort. Additionally, its open structure allows creation of custom control solutions. The efficacy of the ADRC Toolbox is validated through both simulations and hardware experiments, which were conducted using a variety of problems known in the motion, process, and power control areas.


title: The average number of integral points on the congruent number curves

abstract: We show that the total number of non-torsion integral points on the elliptic curves $\mathcal{E}_D:y^2=x^3-D^2x$, where $D$ ranges over positive squarefree integers less than $N$, is $O( N(\log N)^{-1/4+\epsilon})$. The proof involves a discriminant-lowering procedure on integral binary quartic forms and an application of Heath-Brown's method on estimating the average size of the $2$-Selmer group of the curves in this family.


title: Evaluator for Emotionally Consistent Chatbots

abstract: One challenge for evaluating current sequence- or dialogue-level chatbots, such as Empathetic Open-domain Conversation Models, is to determine whether the chatbot performs in an emotionally consistent way. The most recent work only evaluates on the aspects of context coherence, language fluency, response diversity, or logical self-consistency between dialogues. This work proposes training an evaluator to determine the emotional consistency of chatbots.


title: Label noise detection under the Noise at Random model with ensemble filters

abstract: Label noise detection has been widely studied in Machine Learning because of its importance in improving training data quality. Satisfactory noise detection has been achieved by adopting ensembles of classifiers. In this approach, an instance is assigned as mislabeled if a high proportion of members in the pool misclassifies it. Previous authors have empirically evaluated this approach; nevertheless, they mostly assumed that label noise is generated completely at random in a dataset. This is a strong assumption since other types of label noise are feasible in practice and can influence noise detection results. This work investigates the performance of ensemble noise detection under two different noise models: the Noisy at Random (NAR), in which the probability of label noise depends on the instance class, in comparison to the Noisy Completely at Random model, in which the probability of label noise is entirely independent. In this setting, we investigate the effect of class distribution on noise detection performance since it changes the total noise level observed in a dataset under the NAR assumption. Further, an evaluation of the ensemble vote threshold is conducted to contrast with the most common approaches in the literature. In many performed experiments, choosing a noise generation model over another can lead to different results when considering aspects such as class imbalance and noise level ratio among different classes.


title: Bayesian supervised predictive classification and hypothesis testing toolkit for partition exchangeability

abstract: Bayesian supervised predictive classifiers, hypothesis testing, and parametric estimation under Partition Exchangeability are implemented. The two classifiers presented are the marginal classifier (that assumes test data is i.i.d.) next to a more computationally costly but accurate simultaneous classifier (that finds a labelling for the entire test dataset at once based on simultanous use of all the test data to predict each label). We also provide the Maximum Likelihood Estimation (MLE) of the only underlying parameter of the partition exchangeability generative model as well as hypothesis testing statistics for equality of this parameter with a single value, alternative, or multiple samples. We present functions to simulate the sequences from Ewens Sampling Formula as the realisation of the Poisson-Dirichlet distribution and their respective probabilities.


title: A Differential Representation of Cosmological Wavefunctions

abstract: Our understanding of quantum field theory rests largely on explicit and controlled calculations in perturbation theory. Because of this, much recent effort has been devoted to improve our grasp of perturbative techniques on cosmological spacetimes. While scattering amplitudes in flat space at tree level are obtained from simple algebraic operations, things are harder for cosmological observables. Indeed, computing cosmological correlation functions or the associated wavefunction coefficients requires evaluating a growing number of nested time integrals already at tree level, which is computationally challenging. Here, we present a new "differential" representation of the cosmological wavefunction in de Sitter spacetime that obviates this problem for a large class of phenomenologically relevant theories. Given any tree-level Feynman-Witten diagram, we give simple algebraic rules to write down a seed function and a differential operator that transforms it into the desired wavefunction coefficient for any scale-invariant, parity-invariant theory of massless scalars and gravitons with general boost-breaking interactions. In particular, this applies to large classes of phenomenologically relevant theories such as those described by the effective field theory of inflation or solid inflation. Trading nested bulk time integrals for derivatives on boundary kinematical data provides a great computational advantage, especially for processes involving many vertices.


title: Observing, calibrating and reducing near-infrared imaging mosaics

abstract: In near-infrared bands, co-adding and tiling of astronomical imaging datasets require a sufficiently high calibration quality (flat fielding, background subtraction). Here we present a complete workflow for obtaining imaging mosaics with the MMT and Magellan Infrared Spectrograph (MMIRS) operated at the 6.5-m MMT in Arizona and open-source add-on tools developed for the MMIRS pipeline for preparation and data reduction of mosaic observations. We describe pre-observing actions, such as design of dithering patterns and mosaic layouts and post-processing steps to perform absolute astrometric and photometric calibration, and also generate HiPS maps to display the final data product in Aladin / Aladin Lite.


title: Compensating the delayed electro-optical deflector response

abstract: Due to its immediate response, the electro-optic effect is exploited for high-frequency phase and power modulation and for beam deflection. Besides the immediate response, electro-optic materials exhibit creep due to relaxation processes in the material, which results in a gradual settling of the response when a voltage is applied for tens of milliseconds. A bi-exponential filter can compensate this behavior, reported here for electro-optic deflectors made of AD*P, such that these deflectors can be used over the entire frequency range of the driver with a residual settling error <0.001 of the deflection.


title: Explicit Quantum Green Function for Scattering Problems in 2-D Potential

abstract: In this work, we present a new result which concerns the derivation of the Green function relative to the time-independent Schrodinger equation in two dimensional space. The system considered in this work is a quantum particle that have an energy E and moves in an axi-symmetrical potential. Precisely, we have assumed that the potential V(r), in which the quantum particle moves, to be equal to zero inside a disk (radius b) and to be equal a positive constant V0 in a crown of internal radius b and external radius a (b < a) and equal zero outside the crown (r > a). We have explored the diffusion states regime for which E > V0. We have used, to obtain the Green function, the continuity of the solution and of its first derivative at r = b and r = a. We have obtained the associate Green function showing the resonance energies (absence of the reflected waves) for the case E > V0.


title: Time integrator based on rescaled Rodrigues parameters

abstract: We develop an explicit, second-order, variational time integrator for full body dynamics that preserves the momenta of the continuous dynamics, such as linear and angular momenta, and exhibits near-conservation of total energy over exponentially long times. In order to achieve these properties, we parametrize the space of rotations using exponential local coordinates represented by a rescaled form of the Rodrigues rotation vector and we systematically derive the time integrator from a discrete Lagrangian function that yields discrete Euler-Lagrange equations amenable to explicit, closed-form solutions. By restricting attention to spherical bodies and Lagrangian functions with a quadratic kinetic energy and potential energies that solely depend on positions and attitudes, we show that the discrete Lagrangian map exhibits the same mathematical structure, up to terms of second order, of explicit Newmark or velocity Verlet algorithms, both known to be variational time integrators. These preserving properties, together with linear convergence of trajectories and quadratic convergence of total energy, are born out by two examples, namely the dynamics on $\mathrm{SO}(3)$ of a three-dimensional pendulum, and the nonlinear dynamics on $\mathrm{SE}(3)^n$ that results from the impact of a particle-binder torus against a rigid wall.


title: A Survey on Awesome Korean NLP Datasets

abstract: English based datasets are commonly available from Kaggle, GitHub, or recently published papers. Although benchmark tests with English datasets are sufficient to show off the performances of new models and methods, still a researcher need to train and validate the models on Korean based datasets to produce a technology or product, suitable for Korean processing. This paper introduces 15 popular Korean based NLP datasets with summarized details such as volume, license, repositories, and other research results inspired by the datasets. Also, I provide high-resolution instructions with sample or statistics of datasets. The main characteristics of datasets are presented on a single table to provide a rapid summarization of datasets for researchers.


title: Sample-Efficient Generation of Novel Photo-acid Generator Molecules using a Deep Generative Model

abstract: Photo-acid generators (PAGs) are compounds that release acids ($H^+$ ions) when exposed to light. These compounds are critical components of the photolithography processes that are used in the manufacture of semiconductor logic and memory chips. The exponential increase in the demand for semiconductors has highlighted the need for discovering novel photo-acid generators. While de novo molecule design using deep generative models has been widely employed for drug discovery and material design, its application to the creation of novel photo-acid generators poses several unique challenges, such as lack of property labels. In this paper, we highlight these challenges and propose a generative modeling approach that utilizes conditional generation from a pre-trained deep autoencoder and expert-in-the-loop techniques. The validity of the proposed approach was evaluated with the help of subject matter experts, indicating the promise of such an approach for applications beyond the creation of novel photo-acid generators.


title: Finite-temperature electron-capture rates for neutron-rich nuclei around N=50 and effects on core-collapse supernovae simulations

abstract: The temperature dependence of stellar electron-capture (EC) rates is investigated, with a focus on nuclei around $N=50$, just above $Z=28$, which play an important role during the collapse phase of core-collapse supernovae (CCSN). Two new microscopic calculations of stellar EC rates are obtained from a relativistic and a non-relativistic finite-temperature quasiparticle random-phase approximation approaches, for a conventional grid of temperatures and densities. In both approaches, EC rates due to Gamow-Teller transitions are included. In the relativistic calculation contributions from first-forbidden transitions are also included, and add strongly to the EC rates. The new EC rates are compared with large-scale shell model calculations for the specific case of $^{86}$Kr, providing insight into the finite-temperature effects on the EC rates. At relevant thermodynamic conditions for core-collapse, the discrepancies between the different calculations of this work are within about one order of magnitude. Numerical simulations of CCSN are performed with the spherically-symmetric GR1D simulation code to quantify the impact of such differences on the dynamics of the collapse. These simulations also include EC rates based on two parametrized approximations. A comparison of the neutrino luminosities and enclosed mass at core bounce shows that differences between simulations with different sets of EC rates are relatively small ($\approx 5\%$), suggesting that the EC rates used as inputs for these simulations have become well constrained.


title: High-Precision Inversion of Dynamic Radiography Using Hydrodynamic Features

abstract: Radiography is often used to probe complex, evolving density fields in dynamic systems and in so doing gain insight into the underlying physics. This technique has been used in numerous fields including materials science, shock physics, inertial confinement fusion, and other national security applications. In many of these applications, however, complications resulting from noise, scatter, complex beam dynamics, etc. prevent the reconstruction of density from being accurate enough to identify the underlying physics with sufficient confidence. As such, density reconstruction from static/dynamic radiography has typically been limited to identifying discontinuous features such as cracks and voids in a number of these applications. In this work, we propose a fundamentally new approach to reconstructing density from a temporal sequence of radiographic images. Using only the robust features identifiable in radiographs, we combine them with the underlying hydrodynamic equations of motion using a machine learning approach, namely, conditional generative adversarial networks (cGAN), to determine the density fields from a dynamic sequence of radiographs. Next, we seek to further enhance the hydrodynamic consistency of the ML-based density reconstruction through a process of parameter estimation and projection onto a hydrodynamic manifold. In this context, we note that the distance from the hydrodynamic manifold given by the training data to the test data in the parameter space considered both serves as a diagnostic of the robustness of the predictions and serves to augment the training database, with the expectation that the latter will further reduce future density reconstruction errors. Finally, we demonstrate the ability of this method to outperform a traditional radiographic reconstruction in capturing allowable hydrodynamic paths even when relatively small amounts of scatter are present.


title: Spacetime foam, midisuperspace, and the cosmological constant

abstract: Perhaps the cosmological constant really is huge at the Planck scale, but is "hidden" by Planck scale quantum fluctuations of spacetime. I briefly review this proposal and provide some evidence, coming from a simplified midisuperspace model, that an appropriate "foamy" structure can do the job of hiding a large cosmological constant, and can persist under time evolution.


title: Engineering AI Tools for Systematic and Scalable Quality Assessment in Magnetic Resonance Imaging

abstract: A desire to achieve large medical imaging datasets keeps increasing as machine learning algorithms, parallel computing, and hardware technology evolve. Accordingly, there is a growing demand in pooling data from multiple clinical and academic institutes to enable large-scale clinical or translational research studies. Magnetic resonance imaging (MRI) is a frequently used, non-invasive imaging modality. However, constructing a big MRI data repository has multiple challenges related to privacy, data size, DICOM format, logistics, and non-standardized images. Not only building the data repository is difficult, but using data pooled from the repository is also challenging, due to heterogeneity in image acquisition, reconstruction, and processing pipelines across MRI vendors and imaging sites. This position paper describes challenges in constructing a large MRI data repository and using data downloaded from such data repositories in various aspects. To help address the challenges, the paper proposes introducing a quality assessment pipeline, with considerations and general design principles.


title: Achieving the Capacity of a DNA Storage Channel with Linear Coding Schemes

abstract: Due to the redundant nature of DNA synthesis and sequencing technologies, a basic model for a DNA storage system is a multi-draw "shuffling-sampling" channel. In this model, a random number of noisy copies of each sequence is observed at the channel output. Recent works have characterized the capacity of such a DNA storage channel under different noise and sequencing models, relying on sophisticated typicality-based approaches for the achievability. Here, we consider a multi-draw DNA storage channel in the setting of noise corruption by a binary erasure channel. We show that, in this setting, the capacity is achieved by linear coding schemes. This leads to a considerably simpler derivation of the capacity expression of a multi-draw DNA storage channel than existing results in the literature.


title: The numerical solution of semidiscrete linear evolution problems on the finite interval using the Unified Transform Method

abstract: We study a semidiscrete analogue of the Unified Transform Method introduced by A. S. Fokas, to solve initial-boundary-value problems for linear evolution partial differential equations with constant coefficients on the finite interval $x \in (0,L)$. The semidiscrete method is applied to various spatial discretizations of several first and second-order linear equations, producing the exact solution for the semidiscrete problem, given appropriate initial and boundary data. From these solutions, we derive alternative series representations that are better suited for numerical computations. In addition, we show how the Unified Transform Method treats derivative boundary conditions and ghost points introduced by the choice of discretization stencil and we propose the notion of "natural" discretizations. We consider the continuum limit of the semidiscrete solutions and compare with standard finite-difference schemes.


title: An analytical reconstruction formula with efficient implementation for a modality of Compton Scattering Tomography with translational geometry

abstract: In this paper, we address an alternative formulation for the exact inverse formula of the Radon transform on circle arcs arising in a modality of Compton Scattering Tomography in translational geometry proposed by Webber and Miller (Inverse Problems (36)2, 025007, 2020). The original study proposes a first method of reconstruction, using the theory of Volterra integral equations. The numerical realization of such a type of inverse formula may exhibit some difficulties, mainly due to stability issues. Here, we provide a suitable formulation for exact inversion that can be straightforwardly implemented in the Fourier domain. Simulations are carried out to illustrate the efficiency of the proposed reconstruction algorithm.


title: Graph Self-Supervised Learning for Optoelectronic Properties of Organic Semiconductors

abstract: The search for new high-performance organic semiconducting molecules is challenging due to the vastness of the chemical space, machine learning methods, particularly deep learning models like graph neural networks (GNNs), have shown promising potential to address such challenge. However, practical applications of GNNs for chemistry are often limited by the availability of labelled data. Meanwhile, unlabelled molecular data is abundant and could potentially be utilized to alleviate the scarcity issue of labelled data. Here, we advocate the use of self-supervised learning to improve the performance of GNNs by pre-training them with unlabeled molecular data. We investigate regression problems involving ground and excited state properties, both relevant for optoelectronic properties of organic semiconductors. Additionally, we extend the self-supervised learning strategy to molecules in non-equilibrium configurations which are important for studying the effects of disorder. In all cases, we obtain considerable performance improvement over results without pre-training, in particular when labelled training data is limited, and such improvement is attributed to the capability of self-supervised learning in identifying structural similarity among unlabeled molecules.


title: Optimizing frequency allocation for fixed-frequency superconducting quantum processors

abstract: Fixed-frequency superconducting quantum processors are one of the most mature quantum computing architectures with high-coherence qubits and simple controls. However, high-fidelity multi-qubit gates pose tight requirements on individual qubit frequencies in these processors , and these constraints are difficult to satisfy when constructing larger processors due to the large dispersion in the fabrication of Josephson junctions. In this article, we propose a mixed-integer-programming-based optimization approach that determines qubit frequencies to maximize the fabrication yield of quantum processors. We study traditional qubit and qutrit (three-level) architectures with cross-resonance interaction processors. We compare these architectures to a differential AC-Stark shift based on entanglement gates and show that our approach greatly improves the fabrication yield and also increases the scalability of these devices. Our approach is general and can be adapted to problems where one must avoid specific frequency collisions.


title: A Grounded Theory Based Approach to Characterize Software Attack Surfaces

abstract: The notion of Attack Surface refers to the critical points on the boundary of a software system which are accessible from outside or contain valuable content for attackers. The ability to identify attack surface components of software system has a significant role in effectiveness of vulnerability analysis approaches. Most prior works focus on vulnerability techniques that use an approximation of attack surfaces and there has not been many attempt to create a comprehensive list of attack surface components. Although limited number of studies have focused on attack surface analysis, they defined attack surface components based on project specific hypotheses to evaluate security risk of specific types of software applications. In this study, we leverage a qualitative analysis approach to empirically identify an extensive list of attack surface components. To this end, we conduct a Grounded Theory (GT) analysis on 1444 previously published vulnerability reports and weaknesses with a team of three software developers and security experts. We extract vulnerability information from two publicly available repositories: 1) Common Vulnerabilities and Exposures, and 2) Common Weakness Enumeration. We ask three key questions: where the attacks come from, what they target, and how they emerge, and to help answer these questions we define three core categories for attack surface components: Entry points, Targets, and Mechanisms. We extract attack surface concepts related to each category from collected vulnerability information using the GT analysis and provide a comprehensive categorization that represents attack surface components of software systems from various perspectives. The comparison of the proposed attack surface model with the literature shows in the best case previous works cover only 50% of the attack surface components at network level and only 6.7% of the components at code level.


title: Empirical phi-divergence test statistics in the logistic regression model

abstract: In this paper we apply divergence measures to empirical likelihood applied to logistic regression models. We define a family of empirical test statistics based on divergence measures, called empirical phi-divergence test statistics, extending the empirical likelihood ratio test. We study the asymptotic distribution of these empirical test statistics, showing that it is the same for all the test statistics in this family, and the same as the classical empirical likelihood ratio test. Next, we study the power function for the members in this family, showing that the empirical phi-divergence tests introduced in the paper are consistent in the Fraser sense. In order to compare the differences in behavior among the empirical phi-divergence test statistics in this new family, considered for the first time in this paper, we carry out a simulation study.


title: AdaSplit: Adaptive Trade-offs for Resource-constrained Distributed Deep Learning

abstract: Distributed deep learning frameworks like federated learning (FL) and its variants are enabling personalized experiences across a wide range of web clients and mobile/IoT devices. However, FL-based frameworks are constrained by computational resources at clients due to the exploding growth of model parameters (eg. billion parameter model). Split learning (SL), a recent framework, reduces client compute load by splitting the model training between client and server. This flexibility is extremely useful for low-compute setups but is often achieved at cost of increase in bandwidth consumption and may result in sub-optimal convergence, especially when client data is heterogeneous. In this work, we introduce AdaSplit which enables efficiently scaling SL to low resource scenarios by reducing bandwidth consumption and improving performance across heterogeneous clients. To capture and benchmark this multi-dimensional nature of distributed deep learning, we also introduce C3-Score, a metric to evaluate performance under resource budgets. We validate the effectiveness of AdaSplit under limited resources through extensive experimental comparison with strong federated and split learning baselines. We also present a sensitivity analysis of key design choices in AdaSplit which validates the ability of AdaSplit to provide adaptive trade-offs across variable resource budgets.


title: Mixed Stochastic-Deterministic Time-Dependent Density Functional Theory: Application to Stopping Power of Warm Dense Carbon

abstract: Warm dense matter (WMD) describes an intermediate phase, between condensed matter and classical plasmas, found in natural and man-made systems. In a laboratory setting, WDM needs to be created dynamically. It is typically laser or pulse-power generated and can be difficult to characterize experimentally. Measuring the energy loss of high energy ions, caused by a WDM target, is both a promising diagnostic and of fundamental importance to inertial confinement fusion research. However, electron coupling, degeneracy, and quantum effects limit the accuracy of easily calculable kinetic models for stopping power, while high temperatures make the traditional tools of condensed matter, e.g. Time-Dependent Density Functional Theory (TD-DFT), often intractable. We have developed a mixed stochastic-deterministic approach to TD-DFT which provides more efficient computation while maintaining the required precision for model discrimination. Recently, this approach showed significant improvement compared to models when compared to experimental energy loss measurements in WDM carbon. Here, we describe this approach and demonstrate its application to warm dense carbon stopping acr

categories: physics.plasm-ph cond-mat.other hep-th physics.comp-ph


title: Patient-Centered Appraisal of Race-Free Clinical Risk Assessment

abstract: Until recently, there has been a consensus that clinicians should condition patient risk assessments on all observed patient covariates with predictive power. The broad idea is that knowing more about patients enables more accurate predictions of their health risks and, hence, better clinical decisions. This consensus has recently unraveled with respect to a specific covariate, namely race. There have been increasing calls for race-free risk assessment, arguing that using race to predict patient outcomes contributes to racial disparities and inequities in health care. Writers calling for race-free risk assessment have not studied how it would affect the quality of clinical decisions. Considering the matter from the patient-centered perspective of medical economics yields a disturbing conclusion: Race-free risk assessment would harm patients of all races.


title: MultiVerS: Improving scientific claim verification with weak supervision and full-document context

abstract: The scientific claim verification task requires an NLP system to label scientific documents which Support or Refute an input claim, and to select evidentiary sentences (or rationales) justifying each predicted label. In this work, we present MultiVerS, which predicts a fact-checking label and identifies rationales in a multitask fashion based on a shared encoding of the claim and full document context. This approach accomplishes two key modeling goals. First, it ensures that all relevant contextual information is incorporated into each labeling decision. Second, it enables the model to learn from instances annotated with a document-level fact-checking label, but lacking sentence-level rationales. This allows MultiVerS to perform weakly-supervised domain adaptation by training on scientific documents labeled using high-precision heuristics. Our approach outperforms two competitive baselines on three scientific claim verification datasets, with particularly strong performance in zero / few-shot domain adaptation experiments. Our code and data are available at https://github.com/dwadden/multivers.


title: Hamiltonian latent operators for content and motion disentanglement in image sequences

abstract: We introduce \textit{HALO} -- a deep generative model utilising HAmiltonian Latent Operators to reliably disentangle content and motion information in image sequences. The \textit{content} represents summary statistics of a sequence, and \textit{motion} is a dynamic process that determines how information is expressed in any part of the sequence. By modelling the dynamics as a Hamiltonian motion, important desiderata are ensured: (1) the motion is reversible, (2) the symplectic, volume-preserving structure in phase space means paths are continuous and are not divergent in the latent space. Consequently, the nearness of sequence frames is realised by the nearness of their coordinates in the phase space, which proves valuable for disentanglement and long-term sequence generation. The sequence space is generally comprised of different types of dynamical motions. To ensure long-term separability and allow controlled generation, we associate every motion with a unique Hamiltonian that acts in its respective subspace. We demonstrate the utility of \textit{HALO} by swapping the motion of a pair of sequences, controlled generation, and image rotations.


title: Probabilistic Contrastive Loss for Self-Supervised Learning

abstract: This paper proposes a probabilistic contrastive loss function for self-supervised learning. The well-known contrastive loss is deterministic and involves a temperature hyperparameter that scales the inner product between two normed feature embeddings. By reinterpreting the temperature hyperparameter as a quantity related to the radius of the hypersphere, we derive a new loss function that involves a confidence measure which quantifies uncertainty in a mathematically grounding manner. Some intriguing properties of the proposed loss function are empirically demonstrated, which agree with human-like predictions. We believe the present work brings up a new prospective to the area of contrastive learning.


title: The critical layer in quadratic flow boundary layers over acoustic linings

abstract: A straight cylindrical duct is considered containing an axial mean flow that is uniform everywhere except within a boundary layer near the wall, which need not be thin. Within this boundary layer the mean flow varies parabolically. The linearized Euler equations are Fourier transformed to give the Pridmore-Brown equation, for which the Greens function is constructed using Frobenius series. Inverting the spatial Fourier transform, the critical layer contribution is given as the non-modal contribution from integrating around the continuous spectrum branch cut. This contribution is found to be the dominant downstream contribution to the pressure perturbation in certain cases, particularly for thicker boundary layers. The continuous spectrum branch cut is also found to stabilize what are otherwise convectively unstable modes by hiding them behind the branch cut. Overall, the contribution from the critical layer is found to give a neutrally stable non-modal wave with a phase velocity equal to the mean flow velocity at the source when the source is located within the sheared-flow region, and to decay algebraically along the duct as $O(x^{-5/2})$ for a source located with the uniform flow region. The Frobenius expansion, in addition to being numerically accurate close to the critical layer where other numerical methods loose accuracy, is also able to locate modal poles hidden behind the branch cut, which other methods are unable to find; this includes the stabilized hydrodynamic instability. Matlab code is provided to compute the Greens function.


title: Systematically reviewing the layered architectural pattern principles and their use to reconstruct software architectures

abstract: Architectural reconstruction is a reverse engineering activity aiming at recovering the missing decisions on a system. It can help identify the components, within a legacy software application, according to the application's architectural pattern. It is useful to identify architectural technical debt. We are interested in identifying layers within a layered application since the layered pattern is one of the most used patterns to structure large systems. Earlier component reconstruction work focusing on that pattern relied on generic component identification criteria, such as cohesion and coupling. Recent work has identified architectural-pattern specific criteria to identify components within that pattern. However, the architectural-pattern specific criteria that the layered pattern embodies are loosely defined. In this paper, we present a first systematic literature review (SLR) of the literature aiming at inventorying such criteria for layers within legacy applications and grouping them under four principles that embody the fundamental design principles under-lying the architectural pattern. We identify six such criteria in the form of design rules. We also perform a second systematic literature review to synthesize the literature on software architecture reconstruction in the light of these criteria. We report those principles, the rules they encompass, their representation, and their usage in software architecture reconstruction.


title: Winding and intersection of Brownian motions

abstract: We study the set of points $\mathcal{D}_{n,m}$ around which two independent Brownian motions wind at least $n$ (resp. $m$) times. We prove that its area is asymptotically equivalent, in $L^p$ and almost surely, to $\frac{\ell(\mathbb{R}^2)}{4\pi^2 n m}$, where $\ell$ is the intersection measure of the two trajectories. We also prove that the properly scaled Lebesgue measure carried by $\mathcal{D}_{n,m}$ converges almost surely weakly toward $\ell$.


title: Investigating the usefulness of Quantum Blur

abstract: Though some years remain before quantum computation can fully outperform conventional computation, it already provides resources that can be used for exploratory purposes in various fields. This includes certain tasks for procedural generation in computer games, music and art. The so-called `Quantum Blur' method represents the first step on this journey, providing a simple proof-of-principle example of how quantum software can be useful in these areas today. Here we analyse the `Quantum Blur' method and compare it to conventional blur effects. This investigation was guided by discussions with the most prominent user of the method, to determine which features were found most useful. In particular we determine how these features depend on the quantum phenomena of superposition and entanglement.


title: Explicit Abelian Lifts and Quantum LDPC Codes

abstract: For an abelian group $H$ acting on the set $[\ell]$, an $(H,\ell)$-lift of a graph $G_0$ is a graph obtained by replacing each vertex by $\ell$ copies, and each edge by a matching corresponding to the action of an element of $H$. In this work, we show the following explicit constructions of expanders obtained via abelian lifts. For every (transitive) abelian group $H \leqslant \text{Sym}(\ell)$, constant degree $d \ge 3$ and $\epsilon > 0$, we construct explicit $d$-regular expander graphs $G$ obtained from an $(H,\ell)$-lift of a (suitable) base $n$-vertex expander $G_0$ with the following parameters: (i) $\lambda(G) \le 2\sqrt{d-1} + \epsilon$, for any lift size $\ell \le 2^{n^{\delta}}$ where $\delta=\delta(d,\epsilon)$, (ii) $\lambda(G) \le \epsilon \cdot d$, for any lift size $\ell \le 2^{n^{\delta_0}}$ for a fixed $\delta_0 > 0$, when $d \ge d_0(\epsilon)$, or (iii) $\lambda(G) \le \widetilde{O}(\sqrt{d})$, for lift size ``exactly'' $\ell = 2^{\Theta(n)}$. As corollaries, we obtain explicit quantum lifted product codes of Panteleev and Kalachev of almost linear distance (and also in a wide range of parameters) and explicit classical quasi-cyclic LDPC codes with wide range of circulant sizes. Items $(i)$ and $(ii)$ above are obtained by extending the techniques of Mohanty, O'Donnell and Paredes [STOC 2020] for $2$-lifts to much larger abelian lift sizes (as a byproduct simplifying their construction). This is done by providing a new encoding of special walks arising in the trace power method, carefully "compressing'" depth-first search traversals. Result $(iii)$ is via a simpler proof of Agarwal et al. [SIAM J. Discrete Math 2019] at the expense of polylog factors in the expansion.


title: Heralded single-pixel imaging with high loss-resistance and noise-robustness

abstract: Imaging with non-classically correlated photon-pairs takes advantages over classical limits in terms of sensitivity and signal-to-noise ratio. However, it is still a challenge to achieve a strong resilience to background noise and losses for practical applications. In this work, we present heralded single-pixel imaging that is remarkably robust against bright background noise and severe signal losses. Using a strong temporal correlation between a photon-pair and joint measurement-based imaging method, we achieve the suppression of noise up to 1000 times larger than the signal and also demonstrate the correlation-induced SNR enhancement factor of over 200 against 70 times larger noise and a 90% signal loss compared to non-time-gated classical imaging. Our work enables correlated imaging with a highly scalable photon capacity.


title: Floquet Engineering of Two Dimensional Photonic Waveguide Arrays with $\pi$ or $\pm2\pi/3$ Corner states

abstract: In this paper, we theoretically study the Floquet engineering of two dimensional photonic waveguide arrays in three types of lattices: honeycomb lattice with Kekule distortion, breathing square lattice and breathing Kagome lattice. The Kekule distortion factor or the breathing factor in the corresponding lattice is periodically changed along the axial direction of the photonic waveguide with frequency $\omega$. Within certain ranges of $\omega$, the Floquet corner states in the Floquet band gap of quasi-energy spectrum are found, which are localized at the corner of the finite two-dimensional arrays. Due to particle-hole symmetric in the model of honeycomb and square lattice, the quasi-energy level of the Floquet $\pi$ corner states is $\pm\omega/2$. On the other hand, Kagome lattice does not have particle-hole symmetric, so that the quasi-energy level of the Floquet $\pm2\pi/3$ corner states is near to $\pm1\omega/3$. The corner states are either protected by crystalline symmetry or reflection symmetry. The finding of Floquet fractional-$\pi$ corner states could provide more options for engineering of on-chip photonic devices.


title: The impact of varying electrical stimulation parameters on neuromuscular response

abstract: High density neurostimulation systems are coming to market to help spinal cord injury patients by stimulating and recording neuromuscular function. However, the parameter space that these systems have to explore is exceedingly large, and would need an artificial intelligence (AI) system to optimize. We need a platform that will allow us to determine the optimal parameter space for these systems. Our project aims to build a platform for mapping and controlling neuromuscular activity, as a high-throughput testbed for implementing and testing closed-loop neuromuscular activity. This abstract presents the first phase (the mapping phase) of building that testbed by combining multi-electrode stimulation/recording with visual motion-tracking. A 3D-printed rectangular raceway was used with 4 pairs of differential recording electrodes, and two stimulation electrodes embedded in the raceway bed. Non-anesthetized earthworms were placed on the raceway with their head section on the stimulating electrodes. Bipolar sinusoidal stimulation pulses of a range of voltages (2 to 6Vp-p), pulse durations (2 ms to 6.7 ms), and a burst rate of 1 pulse per second were applied, and action potentials and physical motion were recorded and analyzed. Action potentials were found to correlate with expansion/contraction displacements of worm segments, and voltage increases were shown to increase action potential propagation amplitude. Using the multiple electrode recording allowed us to capture the wave propagation of action potential pulse over the length of the worm. Feasibility of a platform to simultaneously monitor action potentials and motion of earthworms with real-time mapping was demonstrated.


title: Multi-modal application: Image Memes Generation

abstract: Meme is an interesting word. Internet memes offer unique insights into the changes in our perception of the world, the media and our own lives. If you surf the Internet for long enough, you will see it somewhere on the Internet. With the rise of social media platforms and convenient image dissemination, Image Meme has gained fame. Image memes have become a kind of pop culture and they play an important role in communication over social media, blogs, and open messages. With the development of artificial intelligence and the widespread use of deep learning, Natural Language Processing (NLP) and Computer Vision (CV) can also be used to solve more problems in life, including meme generation. An Internet meme commonly takes the form of an image and is created by combining a meme template (image) and a caption (natural language sentence). In our project, we propose an end-to-end encoder-decoder architecture meme generator. For a given input sentence, we use the Meme template selection model to determine the emotion it expresses and select the image template. Then generate captions and memes through to the meme caption generator. Code and models are available at github


title: Data-Enabled Gradient Flow as Feedback Controller: Regulation of Linear Dynamical Systems to Minimizers of Unknown Functions

abstract: This paper considers the problem of regulating a linear dynamical system to the solution of a convex optimization problem with an unknown or partially-known cost. We design a data-driven feedback controller - based on gradient flow dynamics - that (i) is augmented with learning methods to estimate the cost function based on infrequent (and possibly noisy) functional evaluations; and, concurrently, (ii) is designed to drive the inputs and outputs of the dynamical system to the optimizer of the problem. We derive sufficient conditions on the learning error and the controller gain to ensure that the error between the optimizer of the problem and the state of the closed-loop system is ultimately bounded; the error bound accounts for the functional estimation errors and the temporal variability of the unknown disturbance affecting the linear dynamical system. Our results directly lead to exponential input-to-state stability of the closed-loop system. The proposed method and the theoretical bounds are validated numerically.


title: Learning Curves for Continual Learning in Neural Networks: Self-Knowledge Transfer and Forgetting

abstract: Sequential training from task to task is becoming one of the major objects in deep learning applications such as continual learning and transfer learning. Nevertheless, it remains unclear under what conditions the trained model's performance improves or deteriorates. To deepen our understanding of sequential training, this study provides a theoretical analysis of generalization performance in a solvable case of continual learning. We consider neural networks in the neural tangent kernel (NTK) regime that continually learn target functions from task to task, and investigate the generalization by using an established statistical mechanical analysis of kernel ridge-less regression. We first show characteristic transitions from positive to negative transfer. More similar targets above a specific critical value can achieve positive knowledge transfer for the subsequent task while catastrophic forgetting occurs even with very similar targets. Next, we investigate a variant of continual learning which supposes the same target function in multiple tasks. Even for the same target, the trained model shows some transfer and forgetting depending on the sample size of each task. We can guarantee that the generalization error monotonically decreases from task to task for equal sample sizes while unbalanced sample sizes deteriorate the generalization. We respectively refer to these improvement and deterioration as self-knowledge transfer and forgetting, and empirically confirm them in realistic training of deep neural networks as well.


title: A new family of minimal ideal triangulations of cusped hyperbolic 3-manifolds

abstract: Previous work of the authors with Bus Jaco determined a lower bound on the complexity of cusped hyperbolic 3-manifolds and showed that it is attained by the monodromy ideal triangulations of once-punctured torus bundles. This paper exhibits an infinite family of minimal ideal triangulations of Dehn fillings on the link $8^3_9$ that also attain this lower bound on complexity.


title: Quantum Algorithm for Solving a Quadratic Nonlinear System of Equations

abstract: Solving a quadratic nonlinear system of equations (QNSE) is a fundamental, but important, task in nonlinear science. We propose an efficient quantum algorithm for solving $n$-dimensional QNSE. Our algorithm embeds QNSE into a finite-dimensional system of linear equations using the homotopy perturbation method and a linearization technique; then we solve the linear equations with a quantum linear system solver and obtain a state which is $\epsilon$-close to the normalized exact solution of the QNSE with success probability $\Omega(1)$. The complexity of our algorithm is $O({\rm polylog}(n/\epsilon))$, which provides an exponential improvement over the optimal classical algorithm in dimension $n$, and the dependence on $\epsilon$ is almost optimal. Therefore, our algorithm exponentially accelerates the solution of QNSE and has wide applications in all kinds of nonlinear problems, contributing to the research progress of nonlinear science.


title: Multistability of globally coupled Duffing oscillators

abstract: We analyze the collective dynamics of an ensemble of globally coupled, externally forced, identical mechanical oscillators with cubic nonlinearity. Focus is put on solutions where the ensemble splits into two internally synchronized clusters, as a consequence of the bistability of individual oscillators. The multiplicity of these solutions, induced by the many possible ways of distributing the oscillators between the two clusters, implies that the ensemble can exhibit multistability. As the strength of coupling grows, however, the two-cluster solutions are replaced by a state of full synchronization. By a combination of analytical and numerical techniques, we study the existence and stability of two-cluster solutions. The role of the distribution of oscillators between the clusters and the relative prevalence of the two stable solutions are disclosed.


title: Limitations of Linear Cross-Entropy as a Measure for Quantum Advantage

abstract: Demonstrating quantum advantage requires experimental implementation of a computational task that is hard to achieve using state-of-the-art classical systems. One approach is to perform sampling from a probability distribution associated with a class of highly entangled many-body wavefunctions. It has been suggested that this approach can be certified with the Linear Cross-Entropy Benchmark (XEB). We critically examine this notion. First, in a "benign" setting where an honest implementation of noisy quantum circuits is assumed, we characterize the conditions under which the XEB approximates the fidelity. Second, in an "adversarial" setting where all possible classical algorithms are considered for comparison, we show that achieving relatively high XEB values does not imply faithful simulation of quantum dynamics. We present an efficient classical algorithm that, with 1 GPU within 2s, yields high XEB values, namely 2-12% of those obtained in experiments. By identifying and exploiting several vulnerabilities of the XEB, we achieve high XEB values without full simulation of quantum circuits. Remarkably, our algorithm features better scaling with the system size than noisy quantum devices for commonly studied random circuit ensembles. To quantitatively explain the success of our algorithm and the limitations of the XEB, we use a theoretical framework in which the average XEB and fidelity are mapped to statistical models. We illustrate the relation between the XEB and the fidelity for quantum circuits in various architectures, with different gate choices, and in the presence of noise. Our results show that XEB's utility as a proxy for fidelity hinges on several conditions, which must be checked in the benign setting but cannot be assumed in the adversarial setting. Thus, the XEB alone has limited utility as a benchmark for quantum advantage. We discuss ways to overcome these limitations.


title: Virtual Coset Coding for Encrypted Non-Volatile Memories with Multi-Level Cells

abstract: PCM is a popular backing memory for DRAM main memory in tiered memory systems. PCM has asymmetric access energy; writes dominate reads. MLC asymmetry can vary by an order of magnitude. Many schemes have been developed to take advantage of the asymmetric patterns of 0s and 1s in the data to reduce write energy. Because the memory is non-volatile, data can be recovered via physical attack or across system reboot cycles. To protect information stored in PCM against these attacks requires encryption. Unfortunately, most encryption algorithms scramble 0s and 1s in the data, effectively removing any patterns and negatively impacting schemes that leverage data bias and similarity to reduce write energy. In this paper, we introduce Virtual Coset Coding (VCC) as a workload-independent approach that reduces costly symbol transitions for storing encrypted data. VCC is based on two ideas. First, using coset encoding with random coset candidates, it is possible to effectively reduce the frequency of costly bit/symbol transitions when writing encrypted data. Second, a small set of random substrings can be used to achieve the same encoding efficiency as a large number of random coset candidates, but at a much lower encoding/decoding cost. Additionally, we demonstrate how VCC can be leveraged for energy reduction in combination with fault-mitigation and fault-tolerance to dramatically increase the lifetimes of endurance-limited NVMs, such as PCM. We evaluate the design of VCC and demonstrate that it can be implemented on-chip with only a nominal area overhead. VCC reduces dynamic energy by 22-28% while maintaining the same performance. Using our multi-objective optimization approach achieves at least a 36% improvement in lifetime over the state-of-the-art and at least a 50% improvement in lifetime vs. an unencoded memory, while maintaining its energy savings and system performance.


title: On real-analytic Levi-flat hypersurfaces and holomorphic webs

abstract: We investigate holomorphic webs tangent to real-analytic Levi-flat hypersurfaces on compact complex surfaces. Under certain conditions, we prove that a holomorphic web tangent to a real-analytic Levi-flat hypersurface admits a multiple-valued meromorphic first integral. We also prove that the Levi foliation of a Levi-flat hypersurface induced by an irreducible real-analytic curve in the Grassmannian $G(n+1,n)$ extends to an algebraic web on the complex projective space.


title: The Influence of Data Pre-processing and Post-processing on Long Document Summarization

abstract: Long document summarization is an important and hard task in the field of natural language processing. A good performance of the long document summarization reveals the model has a decent understanding of the human language. Currently, most researches focus on how to modify the attention mechanism of the transformer to achieve a higher ROUGE score. The study of data pre-processing and post-processing are relatively few. In this paper, we use two pre-processing methods and a post-processing method and analyze the effect of these methods on various long document summarization models.


title: Models for characterizing Darrieus turbines and synthesizing wind speed scenarios

abstract: With the evolving global energy dynamics, new and innovative ways to integrate renewable energy technologies into residential and commercial complexes are being adopted. Techno-economic models such as NREL-SAM and HOMER are handy for evaluating such integration of renewables. However, the unavailability of libraries for vertical axis wind turbines (VAWT) and wind speed data limitations have been identified in the two platforms (NREL-SAM and HOMER). Accordingly, this project looks into developing models for evaluating power curve tables for Darrieus turbines and synthesizing wind speed scenarios. Analytical and numerical models for Darrieus turbines have been presented, and the ability to evaluate power curve tables from both models has been demonstrated. Double-Multiple Streamtube model was implemented for the analytical modelling, and 2-D CFD simulations on OpenFOAM were developed for the numerical modelling. The analytical model was found to fail for small-sized turbines, likely due to poor lift and drag coefficient values at high local angles of attack, but matched sufficiently with experimental results for large-sized turbines. Numerical simulations for Darrieus turbines were of two types - torque evaluation at constant wind speed and angular velocity, and starting characteristics for constant wind speed. Auto-regressive moving average models were implemented to characterize historical wind speed measurements, thus synthesizing wind speed scenarios. For this, the wind speeds or wind speed residuals (following detrending) were first normalized using the cumulative distribution functions, following which an ARMA model was fit.


title: FP-Radar: Longitudinal Measurement and Early Detection of Browser Fingerprinting

abstract: Browser fingerprinting is a stateless tracking technique that attempts to combine information exposed by multiple different web APIs to create a unique identifier for tracking users across the web. Over the last decade, trackers have abused several existing and newly proposed web APIs to further enhance the browser fingerprint. Existing approaches are limited to detecting a specific fingerprinting technique(s) at a particular point in time. Thus, they are unable to systematically detect novel fingerprinting techniques that abuse different web APIs. In this paper, we propose FP-Radar, a machine learning approach that leverages longitudinal measurements of web API usage on top-100K websites over the last decade, for early detection of new and evolving browser fingerprinting techniques. The results show that FP-Radar is able to early detect the abuse of newly introduced properties of already known (e.g., WebGL, Sensor) and as well as previously unknown (e.g., Gamepad, Clipboard) APIs for browser fingerprinting. To the best of our knowledge, FP-Radar is also the first to detect the abuse of the Visibility API for ephemeral fingerprinting in the wild.


abstract: We study a Lorentzian version of the well-known Calder\'{o}n problem that is concerned with determination of lower order coefficients in a wave equation on a smooth Lorentzian manifold, given the associated Dirichlet-to-Neumann map. In the earlier work of the authors it was shown that zeroth order coefficients can be uniquely determined under a two-sided spacetime curvature bound and the additional assumption that there are no conjugate points along null or spacelike geodesics. In this paper we show that uniqueness for the zeroth order coefficient holds for manifolds satisfying a weaker curvature bound as well as spacetime perturbations of such manifolds. This relies on a new optimal unique continuation principle for the wave equation in the exterior regions of double null cones. In particular, we solve the Lorentzian Calder\'{o}n problem near the Minkowski geometry.


title: Linking Solar Minimum, Space Weather, and Night Sky Brightness

abstract: This paper presents time series observations and analysis of broadband night sky airglow intensity 4 September 2018 through 30 April 2020. Data were obtained at 5 sites spanning more than 8500 km during the historically deep minimum of Solar Cycle 24 into the beginning of Solar Cycle 25. New time series observations indicate previously unrecognized significant sources of broadband night sky brightness variations, not involving corresponding changes in the Sun's 10.7cm solar flux, occur during deep solar minimum. Even during a deep solar minimum the natural night sky is rarely, if ever, constant in brightness. Changes with time scales of minutes, hours, days, and months are observed. Semiannual night sky brightness variations are coincident with changes in the orientation of Earth's magnetic field relative to the interplanetary magnetic field. Solar wind plasma streams from solar coronal holes arriving at Earth's bow shock nose are coincident with major night sky brightness increase events. Sites more than 8500 km along the Earth's surface experience nights in common with either very bright or very faint night sky airglow emissions. The reason for this observational fact remains an open question. It is plausible, terrestrial night airglow and geomagnetic indices have similar responses to the solar energy input into Earth's magnetosphere. Our empirical results contribute to a quantitative basis for understanding and predicting broadband night sky brightness variations. They are applicable in astronomical, planetary science, space weather, light pollution, biological, and recreational studies.


title: On the stability of isothermal shocks in black hole accretion disks

abstract: Most black holes possess accretion disks. Models of such disks inform observations and constrain the properties of the black holes and their surrounding medium. Here, we study isothermal shocks in a thin black hole accretion flow. Modelling infinitesimal molecular viscosity allows the use of multiple-scales matched asymptotic methods. We thus derive the first explicit calculations of isothermal shock stability. We find that the inner shock is always unstable, and the outer shock is always stable. The growth/decay rates of perturbations depend only on an effective potential and the incoming--outgoing flow difference at the shock location. We give a prescription of accretion regimes in terms of angular momentum and black hole radius. Accounting for angular momentum dissipation implies unstable outer shocks in much of parameter space, even for realistic viscous Reynolds numbers of the order $\approx 10^{20}$.

categories: astro-ph.HE astro-ph.SR math.DS physics.flu-dyn


title: Computing the quantum guesswork: a quadratic assignment problem

abstract: The quantum guesswork quantifies the minimum number of queries needed to guess the state of a quantum ensemble if one is allowed to query only one state at a time. Previous approaches to the computation of the guesswork were based on standard semi-definite programming techniques and therefore lead to approximated results. In contrast, we show that computing the quantum guesswork of qubit ensembles with uniform probability distribution corresponds to solving a quadratic assignment problem and we provide an algorithm that, upon the input of any qubit ensemble over a discrete ring, after finitely many steps outputs the exact closed-form expression of its guesswork. While in general the complexity of our guesswork-computing algorithm is factorial in the number of states, our main result consists of showing a more-than-quadratic speedup for symmetric ensembles, a scenario corresponding to the three-dimensional analog of the maximization version of the turbine-balancing problem. To find such symmetries, we provide an algorithm that, upon the input of any point set over a discrete ring, after finitely many steps outputs its exact symmetries. The complexity of our symmetries-finding algorithm is polynomial in the number of points. As examples, we compute the guesswork of regular and quasi-regular sets of qubit states.


title: The covering numbers of rings

abstract: A cover of an associative (not necessarily commutative nor unital) ring $R$ is a collection of proper subrings of $R$ whose set-theoretic union equals $R$. If such a cover exists, then the covering number $\sigma(R)$ of $R$ is the cardinality of a minimal cover, and a ring $R$ is called $\sigma$-elementary if $\sigma(R) < \sigma(R/I)$ for every nonzero two-sided ideal $I$ of $R$. If $R$ is a ring with unity, then we define the unital covering number $\sigma_u(R)$ to be the size of a minimal cover of $R$ by subrings that contain $1_R$ (if such a cover exists), and $R$ is $\sigma_u$-elementary if $\sigma_u(R) < \sigma_u(R/I)$ for every nonzero two-sided ideal of $R$. In this paper, we classify all $\sigma$-elementary unital rings and determine their covering numbers. Building on this classification, we are further able to classify all $\sigma_u$-elementary rings and prove $\sigma_u(R) = \sigma(R)$ for every $\sigma_u$-elementary ring $R$. We also prove that, if $R$ is a ring without unity with a finite cover, then there exists a unital ring $R'$ such that $\sigma(R) = \sigma_u(R')$, which in turn provides a complete list of all integers that are the covering number of a ring. Moreover, if \[\mathscr{E}(N) := \{m : m \le N, \sigma(R) = m \text{ for some ring } R\},\] then we show that $|\mathscr{E}(N)| = \Theta(N/\log(N))$, which proves that almost all integers are not covering numbers of a ring.


title: Fundamental component enhancement via adaptive nonlinear activation functions

abstract: In many real world oscillatory signals, the fundamental component of a signal $f(t)$ might be weak or does not exist. This makes it difficult to estimate the instantaneous frequency of the signal. Traditionally, researchers apply the rectification trick, working with $|f(t)|$ or $\mbox{ReLu}(f(t))$ instead, to enhance the fundamental component. This raises an interesting question: what type of nonlinear function $g:\mathbb{R} \rightarrow \mathbb{R}$ has the property that $g(f(t))$ has a more pronounced fundamental frequency? $g(t) = |t|$ and $g(t) = \mbox{ReLu}(t)$ seem to work well in practice; we propose a variant of $g(t) = 1/(1-|t|)$ and provide a theoretical guarantee. Several simulated signals and real signals are analyzed to demonstrate the performance of the proposed solution.


title: SDSS-IV MaStar: Data-driven Parameter Derivation for the MaStar Stellar Library

abstract: The MaNGA Stellar Library (MaStar) is a large collection of high-quality empirical stellar spectra designed to cover all spectral types and ideal for use in the stellar population analysis of galaxies observed in the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey. The library contains 59,266 spectra of 24,130 unique stars with spectral resolution $R\sim1800$ and covering a wavelength range of $3,622-10,354$ \r{A}. In this work, we derive five physical parameters for each spectrum in the library: effective temperature ($T_{eff}$), surface gravity ($\log g$), metalicity ($[Fe/H]$), micro-turbulent velocity ($\log(v_{micro})$), and alpha-element abundance ($[\alpha/Fe]$). These parameters are derived with a flexible data-driven algorithm that uses a neural network model. We train a neural network using the subset of 1,675 MaStar targets that have also been observed in the Apache Point Observatory Galactic Evolution Experiment (APOGEE), adopting the independently-derived APOGEE Stellar Parameter and Chemical Abundance Pipeline (ASPCAP) parameters for this reference set. For the regions of parameter space not well represented by the APOGEE training set ($7,000 \leq T \leq 30,000$ K), we supplement with theoretical model spectra. We present our derived parameters along with an analysis of the uncertainties and comparisons to other analyses from the literature.


title: Wave-optical effects in the microlensing of continuous gravitational waves by star clusters

abstract: Rapidly rotating neutron stars are promising sources for existing and upcoming gravitational-wave interferometers. While relatively dim, these systems are expected to emit continuously, allowing for signal to be accumulated through persistent monitoring over year-long timescales. If, at some point during the observational window, the source comes to lie behind a dense collection of stars, transient gravitational lensing may occur. Such events, though rare, would modulate the waveform, induce phase drifts, and ultimately affect parameter inferences concerning the nuclear equation of state and/or magnetic field structure of the neutron star. Importantly, the radiation wavelength will typically exceed the Schwarzschild radius of the individual perturbers in this scenario, implying that (micro-)lensing occurs in the diffractive regime where geometric optics does not apply. In this paper, we make use of numerical tools that borrow from Picard-Lefschetz theory to efficiently evaluate the relevant Fresnel-Kirchhoff integrals for $n \gtrsim 10^{2}$ microlenses. Modulated strain profiles are constructed both in general and for particular neutron star trajectories relative to some simulated macrolenses.


title: An Automatic Approach for Generating Rich, Linked Geo-Metadata from Historical Map Images

abstract: Historical maps contain detailed geographic information difficult to find elsewhere covering long-periods of time (e.g., 125 years for the historical topographic maps in the US). However, these maps typically exist as scanned images without searchable metadata. Existing approaches making historical maps searchable rely on tedious manual work (including crowd-sourcing) to generate the metadata (e.g., geolocations and keywords). Optical character recognition (OCR) software could alleviate the required manual work, but the recognition results are individual words instead of location phrases (e.g., "Black" and "Mountain" vs. "Black Mountain"). This paper presents an end-to-end approach to address the real-world problem of finding and indexing historical map images. This approach automatically processes historical map images to extract their text content and generates a set of metadata that is linked to large external geospatial knowledge bases. The linked metadata in the RDF (Resource Description Framework) format support complex queries for finding and indexing historical maps, such as retrieving all historical maps covering mountain peaks higher than 1,000 meters in California. We have implemented the approach in a system called mapKurator. We have evaluated mapKurator using historical maps from several sources with various map styles, scales, and coverage. Our results show significant improvement over the state-of-the-art methods. The code has been made publicly available as modules of the Kartta Labs project at https://github.com/kartta-labs/Project.


title: An Empirical Estimation of the Galaxy Radio Number Counts model at 3~GHz and the Effect of Gravitational Lensing

abstract: We analyzed the evolution of 3GHz luminosity function of radio selected galaxies through cosmic time, up to about 1 Gyr after Big Bang, with making use of the VLA-COSMOS3GHz large project. We estimated the luminosity function for SFGs and AGNs in the VLA-COSMOS field respectively with two independent methods, the non-parametric C-method for accuracy to avoid the effect of cosmic variance and parametric MCMC method assuming pure luminosity evolution. We obtained pure luminosity evolution parameter as L*=(1+z)(3.36-0.31)z for SFG and L* = (1+z)(2.91-0.99)z for AGN. The resultant model agreed well with observed number counts down to sub-muJy level. Additionally, we compared the model with cosmic star formation rate density history with multiple observations in various wavelength and confirmed that our result is consistent with literature up to z=3. We further analyzed the statistical effect of gravitational lensing due to the dark matter halos and found that the effect is maximized at S_3GHz = 100muJy with about 1% contribution for SFG and less than 0.5% for AGN for the flux range considered in this study. In order to evaluate the effect of lensing magnification on parameter estimation in the SKA observations, we performed Fisher analysis and find that while the bias is limited, there is a distinguishable effect on SFG LF evolution parameter with SKA I MID wide survey, and all-sky and wide survey for AGN LF evolution parameters.


abstract: Coherent states for equally spaced, homogeneous waveguide arrays are defined, in the infinite, semiinfinite and finite cases, and resolutions of the identity are constructed, using different methods. In the infinite case, which corresponds to Euclidean coherent states, a resolution of the identity with coherent states on the circle and involving a nonlocal inner product is reviewed. In the semiinfinite case, which corresponds to London coherent states, various construction are given (restricting to the circle with a non-local scalar product, rescaling the coherent states, modifying them, or using a non-tight frame). In the finite case, a construction in terms of coherent states on the circle is given, and this construction is shown to be a regularization of the infinite and semiinfinite cases.


title: Further study on Horozov-Iliev's method of estimating the number of limit cycles

abstract: In the study of the number of limit cycles of near-Hamiltonian systems, the first order Melnikov function plays an important role. This paper aims to establish a development of a known method to estimate the upper bound of the number of zeros of the function.


title: Challenges and Opportunities in Approximate Bayesian Deep Learning for Intelligent IoT Systems

abstract: Approximate Bayesian deep learning methods hold significant promise for addressing several issues that occur when deploying deep learning components in intelligent systems, including mitigating the occurrence of over-confident errors and providing enhanced robustness to out of distribution examples. However, the computational requirements of existing approximate Bayesian inference methods can make them ill-suited for deployment in intelligent IoT systems that include lower-powered edge devices. In this paper, we present a range of approximate Bayesian inference methods for supervised deep learning and highlight the challenges and opportunities when applying these methods on current edge hardware. We highlight several potential solutions to decreasing model storage requirements and improving computational scalability, including model pruning and distillation methods.


title: A spectroscopic follow-up for Gaia19bld

abstract: Due to their scarcity, microlensing events in the Galactic disk are of great interest and high-cadence photometric observations, supplemented by spectroscopic follow-up, are necessary for constraining the physical parameters of the lensing system. In particular, a precise estimate of the source characteristics is required to accurately measure the lens distance and mass. We conducted a spectroscopic follow-up of microlensing event Gaia19bld to derive the properties of the microlensing source and, ultimately, to estimate the mass and distance of the lens. We obtained low- and high-resolution spectroscopy from multiple sites around the world during the course of the event. The spectral lines and template matching analysis has led to two independent, consistent characterizations of the source. We found that the source is a red giant located at about 8.5 kpc from the Earth. Combining our results with the photometric analysis has led to a lens mass of Ml=1.1 M at a distance of D l = 5.5 kpc. We did not find any significant blend light in the spectra (with an upper detection limit of V < 17 mag), which is in agreement with photometric observations. Therefore, we cannot exclude the possibility that the lens is a main-sequence star. Indeed, we predict in this scenario a lens brightness of V about 20 mag, a value that would make it much fainter than the detection limit.


title: Microscopic Theory of Superconducting Phase Diagram in Infinite-Layer Nickelates

abstract: Since the discovery of superconductivity in infinite-layer nickelates RNiO$_2$ (R=La, Pr, Nd), great research efforts have been paid to unveil its underlying superconducting mechanism. However, the physical origin of the intriguing hole-doped superconductivity phase diagram, characterized by a superconductivity dome sandwiched between two weak insulators, is still unclear. Here, we present a microscopic theory for electronic structure of nickelates from a fundamental model-based perspective. We found that the appearance of weak insulator phase in lightly and heavily hole-doped regime is dominated by Mottness and Hundness, respectively, exhibiting a unique orbital-selective doping originated from the competition of Hund interaction and crystal field splitting. Moreover, the superconducting phase can also be created in the "mixed" transition regime between Mott-insulator and Hund-induced insulator, exactly reproducing the experimentally observed superconducting phase diagram. Our findings not only demonstrate the orbital-dependent strong-correlation physics in Ni 3$d$ states, but also provide a unified understanding of superconducting phase diagram in hole-doped infinite-layer nickelates, which are distinct from the well-established paradigms in cuprates and iron pnictides.


title: Universal Thermal Entanglement of Multichannel Kondo Effects

abstract: Quantum entanglement between an impurity and its environment is expected to be central in quantum impurity problems. We develop a method to compute the entanglement in spin-1/2 impurity problems, based on the entanglement negativity and the boundary conformal field theory (BCFT). Using the method, we study the thermal decay of the entanglement in the multichannel Kondo effects. At zero temperature, the entanglement has the maximal value independent of the number of the screening channels. At low temperature, the entanglement exhibits a power-law thermal decay. The power-law exponent equals two times of the scaling dimension of the BCFT boundary operator describing the impurity spin, and it is attributed to the energy-dependent scaling behavior of the entanglement in energy eigenstates. These agree with numerical renormalization group results, unveiling quantum coherence inside the Kondo screening length.


title: Parametric Resonance of a charged pendulum with suspension point oscillating between two vertical charged lines

abstract: In this work, we study a mathematical planar pendulum whose support point is positioned equidistant between two vertical and uniformly electrically charged wires. Its bob carries an electric charge and, its support point oscillates vertically, following a harmonic law of motion. We study the dynamics of such phenomenon and the parametric resonances of the equilibria. Moreover, we obtain the surface in the parameter space (since such system presents three parameters) which separates the region of stability from the region of instability. On the particular case of zero charge, we obtain the boundary curves of the stability/instability of Matheiu equation.


title: Yukawa textures in modular symmetric vacuum of magnetized orbifold models

abstract: We study quark mass matrices derived from magnetized $T^2/\mathbb{Z}_2$ orbifold models. Yukawa matrices at three modular fixed points, $\tau=i, e^{2\pi i/3}$ and $i\infty$ are invariant under $S$, $ST$ and $T$-transformations. We study these invariances on $T^2/\mathbb{Z}_2$ twisted orbifold. We find that Yukawa matrices have a kind of texture structures although ones at $\tau=i\infty$ are not realistic. We classify Yukawa textures at $\tau=i$ and $e^{2\pi i/3}$. Moreover we investigate the conditions such that quark mass matrix constructed by Yukawa textures becomes approximately rank one matrix, which is favorable to lead to hierarchical masses between the third generation and the others. It is found that realistic quark mass matrices can be obtained around the $S$-invariant vacuum and $ST$-invariant vacuum. As an illustrating example, we show the realization of the quark mass ratios and mixing based on Fritzch and Fritzch-Xing mass matrices.


title: Recent formation and likely cometary activity of near-Earth asteroid pair 2019 PR2 -- 2019 QR6

abstract: Asteroid pairs are genetically related asteroids that recently separated ($<$few million years), but still reside on similar heliocentric orbits. A few hundred of these systems have been identified, primarily in the asteroid main-belt. Here we studied a newly discovered pair of near-Earth objects (NEOs): 2019 PR2 and 2019 QR6. Based on broad-band photometry, we found these asteroids to be spectrally similar to D-types, a type rare amongst NEOs. We recovered astrometric observations for both asteroids from the Catalina Sky Survey from 2005, which significantly improved their fitted orbits. With these refinements we ran backwards orbital integrations to study formation and evolutionary history. We found that neither a pure gravitational model nor a model with the Yarkovsky effect could explain their current orbits. We thus implemented two models of comet-like non-gravitational forces based on water or CO sublimation. The first model assumed quasi-continuous, comet-like activity after separation, which suggested a formation time of the asteroid pair $300^{+120}_{-70}$ years ago. The second model assumed short-term activity for up to one heliocentric orbit ($\sim$13.9 years) after separation, which suggested that the pair formed 272$\pm$7 years ago. Image stacks showed no activity for 2019~PR2 during its last perihelion passage. These results strongly argue for a common origin that makes these objects the youngest asteroid pair known to date. Questions remain regarding whether these objects derived from a parent comet or asteroid, and how activity may have evolved since their separation.


title: Most charming dibaryon near unitarity

abstract: We present a first study on a pair of triply charmed baryons, $\Omega_{ccc}\Omega_{ccc}$ in the $^1S_0$ channel, on the basis of the HAL QCD method. The measurements are perfomed on the $(2+1)$-flavor lattice QCD configurations with nearly physical light-quark masses and physical charm-quark mass. We show that the system with the Coulomb repulsion taking into account the charge form factor of $\Omega_{ccc}$ leads to the scattering length $a^\mathrm{C}_0\simeq-19$ fm and the effective range $r^\mathrm{C}_\mathrm{eff}\simeq0.45$ fm, which indicates $\Omega_{ccc}\Omega_{ccc}$ is located in the unitary regime.


title: TransZero: Attribute-guided Transformer for Zero-Shot Learning

abstract: Zero-shot learning (ZSL) aims to recognize novel classes by transferring semantic knowledge from seen classes to unseen ones. Semantic knowledge is learned from attribute descriptions shared between different classes, which act as strong priors for localizing object attributes that represent discriminative region features, enabling significant visual-semantic interaction. Although some attention-based models have attempted to learn such region features in a single image, the transferability and discriminative attribute localization of visual features are typically neglected. In this paper, we propose an attribute-guided Transformer network, termed TransZero, to refine visual features and learn attribute localization for discriminative visual embedding representations in ZSL. Specifically, TransZero takes a feature augmentation encoder to alleviate the cross-dataset bias between ImageNet and ZSL benchmarks, and improves the transferability of visual features by reducing the entangled relative geometry relationships among region features. To learn locality-augmented visual features, TransZero employs a visual-semantic decoder to localize the image regions most relevant to each attribute in a given image, under the guidance of semantic attribute information. Then, the locality-augmented visual features and semantic vectors are used to conduct effective visual-semantic interaction in a visual-semantic embedding network. Extensive experiments show that TransZero achieves the new state of the art on three ZSL benchmarks. The codes are available at: \url{https://github.com/shiming-chen/TransZero}.


title: A Wide and Deep Exploration of Radio Galaxies with Subaru HSC (WERGS). VI. Distant Filamentary Structures Pointed by High-z Radio Galaxies at z~4

abstract: We present the environmental properties around high-$z$ radio galaxies (HzRGs) at $z\sim4$, which have been poorly investigated because of their rarity. We use the largest samples of HzRGs and $g$-dropout galaxy overdense regions at $z\sim4$, which were constructed from Hyper Suprime-Cam Subaru Strategic Program, to characterize the HzRG environments statistically. We measure the $g$-dropout galaxy overdensities around 21 HzRGs whose rest-frame 1.4 GHz radio luminosities ($L_{1.4\mathrm{GHz}}$) are $10^{26-27}$ W Hz$^{-1}$. We find that the overdensities around the faint HzRGs with $L_{1.4\mathrm{GHz}}\sim10^{26.0-26.5}$ W Hz$^{-1}$ tend to be higher than that of the $g$-dropout galaxies. On the other hand, no significant difference of density environments is found between the luminous HzRGs with $ L_{1.4\mathrm{GHz}}\sim10^{26.5-27.0} $ W Hz$^{-1}$ and the $g$-dropout galaxies. The HzRGs are found to occupy more massive halos than $g$-dropout galaxies through a cross-correlation between the HzRGs and $g$-dropout galaxies. This trend is more pronounced in the faint HzRGs. These results are consistent with a scenario where HzRGs get older and more massive as the radio-luminosity decreases. The HzRGs are expected to trace the progenitors of local cluster halos from their calculated halo mass. In addition, we find that surrounding galaxies tend to distribute along the radio-jet major axis of the HzRGs at angular distances less than $\lesssim500$ physical kpc. Our findings imply the onset of the filamentary structures around the HzRGs at $z\sim4$.


title: Fault-tolerant Identifying Codes in Special Classes of Graphs

abstract: A detection system, modeled in a graph, is composed of "detectors" positioned at a subset of vertices in order to uniquely locate an ``intruder" at any vertex. \emph{Identifying codes} use detectors that can sense the presence or absence of an intruder within distance one. We introduce a fault-tolerant identifying code called a \emph{redundant identifying code}, which allows at most one detector to go offline or be removed without disrupting the detection system. In real-world applications, this would be a necessary feature, as it would allow for maintenance on individual components without disabling the entire system. Specifically, we prove that the problem of determining the lowest cardinality of an identifying code for an arbitrary graph is NP-complete, we determine the bounds on the lowest cardinality for special classes of graphs, including trees, cylinders, and cubic graphs.


title: Make A Long Image Short: Adaptive Token Length for Vision Transformers

abstract: The vision transformer splits each image into a sequence of tokens with fixed length and processes the tokens in the same way as words in natural language processing. More tokens normally lead to better performance but considerably increased computational cost. Motivated by the proverb "A picture is worth a thousand words" we aim to accelerate the ViT model by making a long image short. To this end, we propose a novel approach to assign token length adaptively during inference. Specifically, we first train a ViT model, called Resizable-ViT (ReViT), that can process any given input with diverse token lengths. Then, we retrieve the "token-length label" from ReViT and use it to train a lightweight Token-Length Assigner (TLA). The token-length labels are the smallest number of tokens to split an image that the ReViT can make the correct prediction, and TLA is learned to allocate the optimal token length based on these labels. The TLA enables the ReViT to process the image with the minimum sufficient number of tokens during inference. Thus, the inference speed is boosted by reducing the token numbers in the ViT model. Our approach is general and compatible with modern vision transformer architectures and can significantly reduce computational expanse. We verified the effectiveness of our methods on multiple representative ViT models (DeiT, LV-ViT, and TimesFormer) across two tasks (image classification and action recognition).


title: Differential Property Prediction: A Machine Learning Approach to Experimental Design in Advanced Manufacturing

abstract: Advanced manufacturing techniques have enabled the production of materials with state-of-the-art properties. In many cases however, the development of physics-based models of these techniques lags behind their use in the lab. This means that designing and running experiments proceeds largely via trial and error. This is sub-optimal since experiments are cost-, time-, and labor-intensive. In this work we propose a machine learning framework, differential property classification (DPC), which enables an experimenter to leverage machine learning's unparalleled pattern matching capability to pursue data-driven experimental design. DPC takes two possible experiment parameter sets and outputs a prediction of which will produce a material with a more desirable property specified by the operator. We demonstrate the success of DPC on AA7075 tube manufacturing process and mechanical property data using shear assisted processing and extrusion (ShAPE), a solid phase processing technology. We show that by focusing on the experimenter's need to choose between multiple candidate experimental parameters, we can reframe the challenging regression task of predicting material properties from processing parameters, into a classification task on which machine learning models can achieve good performance.


title: Machine Learning Subsystem for Autonomous Collision Avoidance on a small UAS with Embedded GPU

abstract: Interest in unmanned aerial system (UAS) powered solutions for 6G communication networks has grown immensely with the widespread availability of machine learning based autonomy modules and embedded graphical processing units (GPUs). While these technologies have revolutionized the possibilities of UAS solutions, designing an operable, robust autonomy framework for UAS remains a multi-faceted and difficult problem. In this work, we present our novel, modular framework for UAS autonomy, entitled MR-iFLY, and discuss how it may be extended to enable 6G swarm solutions. We begin by detailing the challenges associated with machine learning based UAS autonomy on resource constrained devices. Next, we describe in depth, how MR-iFLY's novel depth estimation and collision avoidance technology meets these challenges. Lastly, we describe the various evaluation criteria we have used to measure performance, show how our optimized machine vision components provide up to 15X speedup over baseline models and present a flight demonstration video of MR-iFLY's vision-based collision avoidance technology. We argue that these empirical results substantiate MR-iFLY as a candidate for use in reducing communication overhead between nodes in 6G communication swarms by providing standalone collision avoidance and navigation capabilities.


title: Dataset of gold nanoparticle sizes and morphologies extracted from literature-mined microscopy images

abstract: The factors controlling the size and morphology of nanoparticles have so far been poorly understood. Data-driven techniques are an exciting avenue to explore this field through the identification of trends and correlations in data. However, for these techniques to be utilized, large datasets annotated with the structural attributes of nanoparticles are required. While experimental SEM/TEM images collected from controlled experiments are reliable sources of this information, large-scale collection of these images across a variety of experimental conditions is expensive and infeasible. Published scientific literature, which provides a vast source of high-quality figures including SEM/TEM images, can provide a large amount of data at a lower cost if effectively mined. In this work, we develop an automated pipeline to retrieve and analyse microscopy images from gold nanoparticle literature and provide a dataset of 4361 SEM/TEM images of gold nanoparticles along with automatically extracted size and morphology information. The dataset can be queried to obtain information about the physical attributes of gold nanoparticles and their statistical distributions.


title: Quantum time dynamics of 1D-Heisenberg models employing the Yang-Baxter equation for circuit compression

abstract: Quantum time dynamics (QTD) is considered a promising problem for quantum supremacy on near-term quantum computers. However, QTD quantum circuits grow with increasing time simulations. This study focuses on simulating the time dynamics of 1-D integrable spin chains with nearest neighbor interactions. We show how the quantum Yang-Baxter equation can be exploited to compress and produce a shallow quantum circuit. With this compression scheme, the depth of the quantum circuit becomes independent of step size and only depends on the number of spins. We show that the compressed circuit scales quadratically with system size, which allows for the simulations of time dynamics of very large 1-D spin chains. We derive the compressed circuit representations for different special cases of the Heisenberg Hamiltonian. We compare and demonstrate the effectiveness of this approach by performing simulations on quantum computers.


title: Characterization of Low Light Performance of a CMOS sensor for Ultraviolet Astronomical Applications

abstract: CMOS detectors offer many advantages over CCDs for optical and UV astronomical applications, especially in space where high radiation tolerance is required. However, astronomical instruments are most often designed for low light-level observations demanding low dark current and read noise, good linearity and high dynamic range, characteristics that have not been widely demonstrated for CMOS imagers. We report the performance, over temperatures from 140 - 240 K, of a radiation hardened SRI 4Kx2K back-side illuminated CMOS image sensor with surface treatments that make it highly sensitive in blue and UV bands. After suppressing emission from glow sites resulting from defects in the engineering grade device examined in this work, a 0.077 me-/s dark current floor is reached at 160 K, rising to 1 me$^-$/s at 184 K, rivaling that of the best CCDs. We examine the trade-off between readout speed and read noise, finding that 1.43 e$^-$ median read noise is achieved using line-wise digital correlated double sampling at 700 kpix/s/ch corresponding to a 1.5 s readout time. The 15 ke$^-$ well capacity in high gain mode extends to 120 ke$^-$ in dual gain mode. Continued collection of photo-generated charge during readout enables a further dynamic range extension beyond $10^6$ e$^-$ effective well capacity with only 1% loss of exposure efficiency by combining short and long exposures. A quadratic fit to correct for non-linearity reduces gain correction residuals from 1.5% to 0.2% in low gain mode and to 0.4% in high gain mode. Cross-talk to adjacent pixels is only 0.4% vertically, 0.6% horizontally and 0.1% diagonally. These characteristics plus the relatively large 10 $\mu$m pixel size, quasi 4-side buttability, electronic shutter and sub-array readout make this sensor an excellent choice for wide field astronomical imaging in space, even at FUV wavelengths where sky background is very low.


title: An inverse averaging finite element method for solving the size-modified Poisson-Nernst-Planck equations in ion channel simulations

abstract: In this work, we introduce an inverse averaging finite element method (IAFEM) for solving the size-modified Poisson-Nernst-Planck (SMPNP) equations. Comparing with the classical Poisson-Nernst-Planck (PNP) equations, the SMPNP equations add a nonlinear term to each of the Nernst-Planck (NP) fluxes to describe the steric repulsion which can treat multiple nonuniform particle sizes in simulations. Since the new terms include sums and gradients of ion concentrations, the nonlinear coupling of SMPNP equations is much stronger than that of PNP equations. By introducing a generalized Slotboom transform, each of the size-modified NP equation is transformed into a self-adjoint equation with exponentially behaved coefficient, which has similar simple form to the standard NP equation with the Slotboom transformation. This treatment enables employing our recently developed inverse averaging technique to deal with the exponential coefficients of the reformulated formulations, featured with advantages of numerical stability and flux conservation especially in strong nonlinear and convection-dominated cases. Comparing with previous stabilization methods, the IAFEM proposed in this paper can still possess the numerical stability when dealing with convection-dominated problems. And it is more concise and easier to be numerically implemented. Numerical experiments about a model problem with analytic solutions are presented to verify the accuracy and order of IAFEM for SMPNP equations. Studies about the size-effects of a sphere model and an ion channel system are presented to show that our IAFEM is more effective and robust than the traditional finite element method (FEM) when solving SMPNP equations in simulations of biological systems.


title: Optimism brings accurate perception in Iterated Prisoner's Dilemma

abstract: We analyze an extended model of the Iterated Prisoner's Dilemma where agents decide to play based on the data from their limited memory or recommendations. The cooperators can decide whether to play with the matched opponent or not. The agents' decisions are directly linked to their level of optimism since they decide to play if they believe the opponent has a high probability of cooperating. Optimism is precisely tuned by parameters optimism threshold and tolerance. Our experiment showed that being optimistic is better for cooperators as it leads to more accurate exploration in the multi-agent system, which tolerates the vulnerability against defectors.


title: On the Existence of the Adversarial Bayes Classifier (Extended Version)

abstract: Adversarial robustness is a critical property in a variety of modern machine learning applications. While it has been the subject of several recent theoretical studies, many important questions related to adversarial robustness are still open. In this work, we study a fundamental question regarding Bayes optimality for adversarial robustness. We provide general sufficient conditions under which the existence of a Bayes optimal classifier can be guaranteed for adversarial robustness. Our results can provide a useful tool for a subsequent study of surrogate losses in adversarial robustness and their consistency properties. This manuscript is the extended and corrected version of the paper \emph{On the Existence of the Adversarial Bayes Classifier} published in NeurIPS 2021. There were two errors in theorem statements in the original paper -- one in the definition of pseudo-certifiable robustness and the other in the measurability of $A^\e$ for arbitrary metric spaces. In this version we correct the errors. Furthermore, the results of the original paper did not apply to some non-strictly convex norms and here we extend our results to all possible norms.