title: Ising formulation of integer optimization problems for utilizing quantum annealing in iterative improvement strategy

abstract: Quantum annealing is a heuristic algorithm for searching the ground state of an Ising model. Heuristic algorithms aim to obtain near-optimal solutions with a reasonable computation time. Accordingly, many algorithms have so far been proposed. In general, the performance of heuristic algorithms strongly depends on the instance of the combinatorial optimization problem to be solved because they escape the local minima in different ways. Therefore, combining several algorithms to exploit their complementary strength is effective for obtaining highly accurate solutions for a wide range of combinatorial optimization problems. However, quantum annealing cannot be used to improve a candidate solution obtained by other algorithms because it starts from an initial state where all spin configurations are found with a uniform probability. In this study, we propose an Ising formulation of integer optimization problems to utilize quantum annealing in the iterative improvement strategy. Our formulation exploits the biased sampling of degenerated ground states in transverse magnetic field quantum annealing. We also analytically show that a first-order phase transition is successfully avoided for a fully connected ferromagnetic Potts model if the overlap between a ground state and a candidate solution exceeds a threshold. The proposed formulation is applicable to a wide range of integer optimization problems and enables us to hybridize quantum annealing with other optimization algorithms.


title: Effective fluid mixture of tensor-multi-scalar gravity

abstract: We apply to tensor-multi-scalar gravity the effective fluid analysis based on the representation of the gravitational scalar field as a dissipative effective fluid. This generalization poses new challenges as the effective fluid is now a complicated mixture of individual fluids mutually coupled to each other and many reference frames are possible for its description. They are all legitimate, although not all convenient for specific problems, and they give rise to different physical interpretations. Two of these frames are highlighted.


title: Pretraining in Deep Reinforcement Learning: A Survey

abstract: The past few years have seen rapid progress in combining reinforcement learning (RL) with deep learning. Various breakthroughs ranging from games to robotics have spurred the interest in designing sophisticated RL algorithms and systems. However, the prevailing workflow in RL is to learn tabula rasa, which may incur computational inefficiency. This precludes continuous deployment of RL algorithms and potentially excludes researchers without large-scale computing resources. In many other areas of machine learning, the pretraining paradigm has shown to be effective in acquiring transferable knowledge, which can be utilized for a variety of downstream tasks. Recently, we saw a surge of interest in Pretraining for Deep RL with promising results. However, much of the research has been based on different experimental settings. Due to the nature of RL, pretraining in this field is faced with unique challenges and hence requires new design principles. In this survey, we seek to systematically review existing works in pretraining for deep reinforcement learning, provide a taxonomy of these methods, discuss each sub-field, and bring attention to open problems and future directions.


title: A size-consistent Gr\"uneisen-quasiharmonic approach for lattice thermal conductivity

abstract: We propose a size-consistent Gr\"uneisen-quasiharmonic approach (GQA) to calculate the lattice thermal conductivity $\kappa_l$ where the Gr\"uneisen parameters that measure the degree of phonon anharmonicity are calculated directly using first-principles calculations. This is achieved by identifying and modifying two existing equations related to the Slack formulae for $\kappa_l$ that suffer from the size-inconsistency problem when dealing with non-monoatomic primitive cells (where the number of atoms in the primitive cell $n$ is greater than one). In conjunction with other thermal parameters such as the acoustic Debye temperature $\theta_a$ that can also be obtained within the GQA, we predict $\kappa_l$ for a range of materials taken from the diamond, zincblende, rocksalt, and wurtzite compounds. The results are compared with that from the experiment and the quasiharmonic Debye model (QDM). We find that in general the prediction of $\theta_a$ is rather consistent among the GQA, experiment, and QDM. However, while the QDM somewhat overestimates the Gr\"uneisen parameters and hence underestimates $\kappa_l$ for most materials, the GQA predicts the experimental trends of Gr\"uneisen parameters and $\kappa_l$ more closely. We expect the GQA with the modified Slack formulae could be used as an effective and practical predictor for $\kappa_l$, especially for crystals with large $n$.


title: Equilibrium thermodynamic properties of binary hard-sphere mixtures from integral equation theory

abstract: The binary additive hard-sphere mixtures have been studied by the Ornstein-Zernike integral equation coupled with the Martynov-Sarkisov (MS) closure approximation. Virial equation of state is computed in the MS approximation. The excess chemical potential for the mixture is evaluated with a closed-form expression based on correlation functions. The excess Helmholtz free energy is obtained using the Euler relation of thermodynamics. Moreover, these thermodynamic quantities are obtained by the Boubl\'{i}k-Mansoori-Carnahan-Starling-Leland (BMCSL) formulas. Our findings for pressure and excess chemical potential for a number of binary sets of the mixtures from the MS approximation show good agreements with those obtained by the BMCSL formulas and available data in literature, having a maximum deviation of $5\%$ for a packing fraction up to 0.5. The maximum deviation of the excess free energy obtained for the mixtures is shown to be $\sim 16\%$ for a packing fraction of 0.5. To our knowledge, this work presents an initial calculation of an excess chemical potential of the system in the MS approximation.


title: Overlapping time of a virtual customer in time-varying many-server queues

abstract: Motivated by the ongoing COVID-19 pandemic, this paper investigates customers' infection risk by evaluating the overlapping time of a virtual customer with others in queueing systems. Most of the current methodologies focus on characterizing the risk in stationary systems, which may not apply to the more practical time-varying systems. As such, we propose an approximation framework that relies on the fluid limit to compute the expected overlapping time in time-varying queueing systems. Simulation experiments verify the accuracy of our approach.


title: Fast Algorithms for $\ell_p$-Regression

abstract: The $\ell_p$-norm regression problem is a classic problem in optimization with wide ranging applications in machine learning and theoretical computer science. The goal is to compute $x^{\star} =\arg\min_{Ax=b}\|x\|_p^p$, where $x^{\star}\in \mathbb{R}^n, A\in \mathbb{R}^{d\times n},b \in \mathbb{R}^d$ and $d\leq n$. Efficient high-accuracy algorithms for the problem have been challenging both in theory and practice and the state of the art algorithms require $poly(p)\cdot n^{\frac{1}{2}-\frac{1}{p}}$ linear system solves for $p\geq 2$. In this paper, we provide new algorithms for $\ell_p$-regression (and a more general formulation of the problem) that obtain a high-accuracy solution in $O(p n^{\frac{(p-2)}{(3p-2)}})$ linear system solves. We further propose a new inverse maintenance procedure that speeds-up our algorithm to $\widetilde{O}(n^{\omega})$ total runtime, where $O(n^{\omega})$ denotes the running time for multiplying $n \times n$ matrices. Additionally, we give the first Iteratively Reweighted Least Squares (IRLS) algorithm that is guaranteed to converge to an optimum in a few iterations. Our IRLS algorithm has shown exceptional practical performance, beating the currently available implementations in MATLAB/CVX by 10-50x.


title: Numerical solution to the time-dependent Gross-Pitaevskii equation

abstract: In this work we employ the split-step technique combined with a Legendre pseudospectral representation to solve various time-dependent Gross-Pitaevskii equations (GPE). Our findings based on the numerical accuracy of this approach applied for one-dimensional (1D) and two-dimensional (2D) problems show that it can provide accurate and stable solutions. Moreover, this approach has been applied to study the dynamics of the Bose-Einstein condensate which is modeled with the GPE. The breathing of condensate with an repulsive and attractive interactions trapped in 1D and 2D harmonic potentials has been simulated as well.


title: CP violation in the interference between $\rho(770)^0$ and $S$-wave in $B^\pm\to \pi^+ \pi^-\pi^\pm$ decays and its implication for CP asymmetry in $B^\pm\to \rho^0\pi^\pm$

abstract: The decay amplitude analyses of $B^\pm\to \pi^+ \pi^-\pi^\pm$ decays in the Dalitz plot performed by LHCb indicate that CP asymmetry for the dominant quasi-two-body decay $B^\pm\to\rho(770)^0\pi^\pm$ was found to be consistent with zero in all three approaches for the $S$-wave component and that CP-violation effects related to the interference between the $\rho(770)^0$ resonance and the $S$-wave were clearly observed. We show that the nearly vanishing CP violation in $B^\pm\to\rho^0\pi^\pm$ is understandable within the framework of QCD factorization. There are two $1/m_b$ power corrections, one from penguin annihilation and the other from hard spectator interactions. They contribute destructively to $A_{CP}(B^\pm\to\rho^0\pi^\pm)$ to render it compatible with zero, in contrast to the sizable negative CP asymmetry predicted in most of the existing models. We next show that the measured interference pattern between the $\rho$ and $S$-wave contributions in the low $m_{\pi^+\pi^-}$ region separated by the sign of value of $\cos\theta$ with $\theta$ being the angle between the two same charged pions measured in the rest frame of the $\rho$ resonance can be explained in terms of the smallness of $A_{CP}(B^+\to\rho^0\pi^+)$ and the interference between $\rho(770)$ and $\sigma/f_0(500)$. If CP asymmetry in $B^\pm\to\rho^0\pi^\pm$ is not negligible as predicted in many existing models, the observed interference pattern will be destroyed. We conclude that the experimental observation of the interference pattern between $P$- and $S$-waves in the low-$m_{\pi^+\pi^-}$ region between 0.5 and 1.0 GeV is consistent with a nearly vanishing CP violation in $B^\pm\to\rho^0\pi^\pm$.


title: Parameter and Data Efficient Continual Pre-training for Robustness to Dialectal Variance in Arabic

abstract: The use of multilingual language models for tasks in low and high-resource languages has been a success story in deep learning. In recent times, Arabic has been receiving widespread attention on account of its dialectal variance. While prior research studies have tried to adapt these multilingual models for dialectal variants of Arabic, it still remains a challenging problem owing to the lack of sufficient monolingual dialectal data and parallel translation data of such dialectal variants. It remains an open problem on whether the limited dialectical data can be used to improve the models trained in Arabic on its dialectal variants. First, we show that multilingual-BERT (mBERT) incrementally pretrained on Arabic monolingual data takes less training time and yields comparable accuracy when compared to our custom monolingual Arabic model and beat existing models (by an avg metric of +$6.41$). We then explore two continual pre-training methods -- (1) using small amounts of dialectical data for continual finetuning and (2) parallel Arabic to English data and a Translation Language Modeling loss function. We show that both approaches help improve performance on dialectal classification tasks ($+4.64$ avg. gain) when used on monolingual models.


title: Noncommutative Schur functions for posets

abstract: The machinery of noncommutative Schur functions is a general approach to Schur positivity of symmetric functions initiated by Fomin-Greene. Hwang recently adapted this theory to posets to give a new approach to the Stanley-Stembridge conjecture. We further develop this theory to prove that the symmetric function associated to any $P$-Knuth equivalence graph is Schur positive. This settles a conjecture of Kim and the third author, and refines results of Gasharov, Shareshian-Wachs, and Hwang on the Schur positivity of chromatic symmetric functions.


title: Electromagnetic fields in ultra-peripheral relativistic heavy-ion collisions

abstract: Ultraperipheral heavy-ion collisions (UPCs) offer unique opportunities to study processes under strong electromagnetic fields. In these collisions, highly charged fast-moving ions carry strong electromagnetic fields that can be effectively treated as photon fluxes. The exchange of photons can induce photonuclear and two-photon interactions, and excite ions. This excitation of the ions results in Coulomb dissociation with the emission of photons, neutrons, and other particles. Additionally, the electromagnetic fields generated by the ions can be sufficiently strong to enforce mutual interactions between the two colliding ions. Consequently, the two colliding ions experienced an electromagnetic force that pushed them in opposite directions, causing a back-to-back correlation in the emitted neutrons. Using a Monte Carlo simulation, we qualitatively demonstrated that the above electromagnetic effect is large enough to be observed in UPCs, which would provide a clear means to study strong electromagnetic fields and their effects.


title: Pitfalls of Zero Voltage Values in Optimal Power Flow Problems

abstract: The existence of strictly positive lower bounds on voltage magnitude is taken for granted in optimal power flow problems. Nevertheless, it is not possible to rely on such bounds for a variety of real-world network optimization problems. This paper discusses a few issues related to 0 V assumptions made during the process of deriving optimization formulations in the current-voltage, power-voltage and power-lifted-voltage variable spaces. The differences between the assumptions are illustrated for a 2-bus 2-wire test case, where the feasible sets are visualized. A nonzero relaxation gap is observed for the canonical multiconductor nonlinear power-voltage formulation. A zero gap can be obtained for the branch flow model semi-definite relaxation, using newly proposed valid equalities.


title: On the Algorithmic Stability and Generalization of Adaptive Optimization Methods

abstract: Despite their popularity in deep learning and machine learning in general, the theoretical properties of adaptive optimizers such as Adagrad, RMSProp, Adam or AdamW are not yet fully understood. In this paper, we develop a novel framework to study the stability and generalization of these optimization methods. Based on this framework, we show provable guarantees about such properties that depend heavily on a single parameter $\beta_2$. Our empirical experiments support our claims and provide practical insights into the stability and generalization properties of adaptive optimization methods.


title: Charge Crowding in Graphene-Silicon Diodes

abstract: The performance of nanoscale electronic devices based on a two-three dimensional (2D-3D) interface is significantly affected by the electrical contacts that interconnect these materials with external circuitry. This work investigates charge transport effects at the 2D-3D ohmic contact coupled with the thermionic injection model for graphene/Si Schottky junction. Here, w e focus on the intrinsic properties of graphene-metal contacts, paying particular attention to the nature of the contact failure mechanism under high electrical stress. According to our findings, severe current crowding (CC) effects in highly conductive electrical contact significantly affect device failure that can be reduced by spatially varying the contact properties and geometry. The impact of electrical breakdown on material degradation is systematically analyzed by atomic force, Raman, scanning electron, and energy dispersive X-ray spectroscopies. Our devices withstand high electrostatic discharge spikes over a longer period, manifesting high robustness and operational stability. This research paves the way towards a highly robust and reliable graphene/Si heterostructure in futuristic on-chip integration in dynamic switching. The methods we employed here can be extended for other nanoscale electronic devices based on 2D-3D interfaces


title: Quantization-Based Optimization: Alternative Stochastic Approximation of Global Optimization

abstract: In this study, we propose a global optimization algorithm based on quantizing the energy level of an objective function in an NP-hard problem. According to the white noise hypothesis for a quantization error with a dense and uniform distribution, we can regard the quantization error as i.i.d. white noise. From stochastic analysis, the proposed algorithm converges weakly only under conditions satisfying Lipschitz continuity, instead of local convergence properties such as the Hessian constraint of the objective function. This shows that the proposed algorithm ensures global optimization by Laplace's condition. Numerical experiments show that the proposed algorithm outperforms conventional learning methods in solving NP-hard optimization problems such as the traveling salesman problem.


title: The Wentaculus: Density Matrix Realism Meets the Arrow of Time

abstract: Two of the most difficult problems in the foundations of physics are (1) what gives rise to the arrow of time and (2) what the ontology of quantum mechanics is. They are difficult because the fundamental dynamical laws of physics do not privilege an arrow of time, and the quantum-mechanical wave function describes a high-dimensional reality that is radically different from our ordinary experiences. In this paper, I characterize and elaborate on the ''Wentaculus'' theory, a new approach to time's arrow in a quantum universe that offers a unified solution to both problems. Central to the Wentaculus are (i) Density Matrix Realism, the idea that the quantum state of the universe is objective but can be impure, and (ii) the Initial Projection Hypothesis, a new law of nature that selects a unique initial quantum state. On the Wentaculus, the quantum state of the universe is sufficiently simple to be a law, and the arrow of time can be traced back to an exact boundary condition. It removes the intrinsic vagueness of the Past Hypothesis, eliminates the Statistical Postulate, provides a higher degree of theoretical unity, and contains a natural realization of ''strong determinism.'' I end by responding to four recent objections. In a companion paper, I elaborate on Density Matrix Realism.

categories: quant-ph cond-mat.stat-mech gr-qc physics.hist-ph


title: Limits on Simultaneous and Delayed Optical Emission from Well-localized Fast Radio Bursts

abstract: We present the largest compilation to date of optical observations during and following fast radio bursts (FRBs). The data set includes our dedicated simultaneous and follow-up observations, as well as serendipitous archival survey observations, for a sample of 15 well-localized FRBs: eight repeating and seven one-off sources. Our simultaneous (and nearly simultaneous with a $0.4$ s delay) optical observations of 13 (1) bursts from the repeating FRB 20220912A provide the deepest such limits to date for any extragalactic FRB, reaching a luminosity limit of $\nu L_\nu\lesssim 10^{42}$ erg s$^{-1}$ ($\lesssim 2\times10^{41}$ erg s$^{-1}$) with $15-400$ s exposures; an optical-flux-to-radio-fluence ratio of $f_{\rm opt}/F_{\rm radio}\lesssim 10^{-7}$ ms$^{-1}$ ($\lesssim 10^{-8}$ ms$^{-1}$); and flux ratio of $f_{\rm opt}/f_{\rm radio}\lesssim 0.02-\lesssim 2\times 10^{-5}$ ($\lesssim 10^{-6}$) on millisecond to second timescales. These simultaneous limits provide useful constraints in the context of FRB emission models, such as the pulsar magnetosphere and pulsar nebula models. Interpreting all available optical limits in the context of the synchrotron maser model, we find that they constrain the flare energies to $\lesssim 10^{43}-10^{49}$ erg (depending on the distances of the various repeating FRBs, with $\lesssim 10^{39}$ erg for the Galactic SGR 1935+2154). These limits are generally at least an order of magnitude larger than those inferred from the FRBs themselves, although in the case of FRB 20220912A our simultaneous and rapid follow-up observations severely restrict the model parameter space. We conclude by exploring the potential of future simultaneous and rapid-response observations with large optical telescopes.


title: Optimal Smoothed Analysis and Quantitative Universality for the Smallest Singular Value of Random Matrices

abstract: The smallest singular value and condition number play important roles in numerical linear algebra and the analysis of algorithms. In numerical analysis with randomness, many previous works make Gaussian assumptions, which are not general enough to reflect the arbitrariness of the input. To overcome this drawback, we prove the first quantitative universality for the smallest singular value and condition number of random matrices. Moreover, motivated by the study of smoothed analysis that random perturbation makes deterministic matrices well-conditioned, we consider an analog for random matrices. For a random matrix perturbed by independent Gaussian noise, we show that this matrix quickly becomes approximately Gaussian. In particular, we derive an optimal smoothed analysis for random matrices in terms of a sharp Gaussian approximation.

categories: math.PR cs.DS cs.NA math.NA math.ST stat.TH


title: The Logic of Cardinality Comparison Without the Axiom of Choice

abstract: We work in the setting of Zermelo-Fraenkel set theory without assuming the Axiom of Choice. We consider sets with the Boolean operations together with the additional structure of comparing cardinality (in the Cantorian sense of injections). What principles does one need to add to the laws of Boolean algebra to reason not only about intersection, union, and complementation of sets, but also about the relative size of sets? We give a complete axiomatization. A particularly interesting case is when one restricts to the Dedekind-finite sets. In this case, one needs exactly the same principles as for reasoning about imprecise probability comparisons, the central principle being Generalized Finite Cancellation (which includes, as a special case, division-by-$m$). In the general case, the central principle is a restricted version of Generalized Finite Cancellation within Archimedean classes which we call Covered Generalized Finite Cancellation.


title: Assemble Them All: Physics-Based Planning for Generalizable Assembly by Disassembly

abstract: Assembly planning is the core of automating product assembly, maintenance, and recycling for modern industrial manufacturing. Despite its importance and long history of research, planning for mechanical assemblies when given the final assembled state remains a challenging problem. This is due to the complexity of dealing with arbitrary 3D shapes and the highly constrained motion required for real-world assemblies. In this work, we propose a novel method to efficiently plan physically plausible assembly motion and sequences for real-world assemblies. Our method leverages the assembly-by-disassembly principle and physics-based simulation to efficiently explore a reduced search space. To evaluate the generality of our method, we define a large-scale dataset consisting of thousands of physically valid industrial assemblies with a variety of assembly motions required. Our experiments on this new benchmark demonstrate we achieve a state-of-the-art success rate and the highest computational efficiency compared to other baseline algorithms. Our method also generalizes to rotational assemblies (e.g., screws and puzzles) and solves 80-part assemblies within several minutes.


title: Robust Functional Magnetoencephalographic Brain Measures with 1.0 Millimeter Spatial Separation

abstract: Neuroelectric currents were extracted from free-running magnetoencephalographic (MEG) rest and task recordings from 617 normative subjects (ages: 18-87). State-dependent neuroelectric differential activation (DA) with spatial resolution comparable to that of local field potentials was detected in the majority of this cohort. Rest-high (rest greater than task) or task-high DA was found in the majority of individual subjects in more than 13,000 1 mm^3 voxels per subject. On average, 6% of the DA voxels bordered a second voxel whose DA was opposite, i.e., one was rest-high and the other was task-high. 516 subjects showed more than 100 such opposite voxel pairs 1 mm apart; 226 subjects showed more than 1000. The number of bordering voxel pairs with the same DA was consistently higher for almost all subjects and averaged 20%, ruling out the possibility that opposite bordering voxels occur simply by chance. For 65 brain regions, more than 10% of the cohort showed significantly more same than opposite pairs. These findings taken together support the conclusion that neuroelectric DA is consistently distinguishable at single 1 mm^3 brain voxels with 1-mm spatial separation. When restricted to voxels with near-zero rest or task counts, significantly more rest-high than task-high voxels were found in 35 regions for at least 10 percent of the subjects. This inequality was not found when all DA-voxels were included. This supports the conclusion that the DA found in many rest-high voxels with near-zero task counts is due in part to task-dependent inhibition.


title: AI Testing Framework for Next-G O-RAN Networks: Requirements, Design, and Research Opportunities

abstract: Openness and intelligence are two enabling features to be introduced in next generation wireless networks, e.g. Beyond 5G and 6G, to support service heterogeneity, open hardware, optimal resource utilization, and on-demand service deployment. The open radio access network (O-RAN) is a promising RAN architecture to achieve both openness and intelligence through virtualized network elements and well-defined interfaces. While deploying artificial intelligence (AI) models is becoming easier in O-RAN, one significant challenge that has been long neglected is the comprehensive testing of their performance in realistic environments. This article presents a general automated, distributed and AI-enabled testing framework to test AI models deployed in O-RAN in terms of their decision-making performance, vulnerability and security. This framework adopts a master-actor architecture to manage a number of end devices for distributed testing. More importantly, it leverages AI to automatically and intelligently explore the decision space of AI models in O-RAN. Both software simulation testing and software-defined radio hardware testing are supported, enabling rapid proof of concept research and experimental research on wireless research platforms.


title: Vortex solitons in quasi-phase-matched photonic crystals

abstract: We report solutions for stable compound solitons in a three-dimensional quasi-phase-matched photonic crystal with the quadratic ($\chi ^{(2)}$) nonlinearity. The photonic crystal is introduced with a checkerboard structure, which can be realized by means of the available technology. The solitons are built as four-peak vortex modes of two types, rhombuses and squares (intersite- and onsite-centered self-trapped states, respectively). Their stability areas are identified in the system's parametric space (rhombuses occupy an essentially broader stability domain), while all bright vortex solitons are subject to strong azimuthal instability in uniform $\chi^{(2)}$ media. Possibilities for experimental realization of the solitons are outlined.


title: Ab initio study of proton-exchanged LiNbO3(I): Structural, thermodynamic, dielectric, and optical properties

abstract: Using first principles calculations, we study the ground-state structure of bulk proton-exchanged lithium niobate, which is also called hydrogen niobate and is widely used in waveguides. Thermodynamics helps to establish the most favorable nonpolar surface as well as the water-deficient and water-rich phases under different ambient conditions, which we refer to as "dehydrated" and "rehydrated" phases, respectively. We compute the low-frequency dielectric response and the optical refractive indices of hydrogen niobate in different phases. The dielectric constant is greatly enhanced compared to lithium niobate. At shorter wavelengths, the refractive indices vary between each phase and have a sharp contrast to lithium niobate. Our study characterizes the structures and thermal instabilities of this compound and reveals its excellent dielectric and optical properties, which can be important in the future application in waveguides.


title: Low regularity integrators for semilinear parabolic equations with maximum bound principles

abstract: This paper is concerned with conditionally structure-preserving, low regularity time integration methods for a class of semilinear parabolic equations of Allen-Cahn type. Important properties of such equations include maximum bound principle (MBP) and energy dissipation law; for the former, that means the absolute value of the solution is pointwisely bounded for all the time by some constant imposed by appropriate initial and boundary conditions. The model equation is first discretized in space by the central finite difference, then by iteratively using Duhamel's formula, first- and second-order low regularity integrators (LRIs) are constructed for time discretization of the semi-discrete system. The proposed LRI schemes are proved to preserve the MBP and the energy stability in the discrete sense. Furthermore, their temporal error estimates are also successfully derived under a low regularity requirement that the exact solution of the semi-discrete problem is only assumed to be continuous in time. Numerical results show that the proposed LRI schemes are more accurate and have better convergence rates than classic exponential time differencing schemes, especially when the interfacial parameter approaches zero.


title: Doubly Inhomogeneous Reinforcement Learning

abstract: This paper studies reinforcement learning (RL) in doubly inhomogeneous environments under temporal non-stationarity and subject heterogeneity. In a number of applications, it is commonplace to encounter datasets generated by system dynamics that may change over time and population, challenging high-quality sequential decision making. Nonetheless, most existing RL solutions require either temporal stationarity or subject homogeneity, which would result in sub-optimal policies if both assumptions were violated. To address both challenges simultaneously, we propose an original algorithm to determine the ``best data chunks" that display similar dynamics over time and across individuals for policy learning, which alternates between most recent change point detection and cluster identification. Our method is general, and works with a wide range of clustering and change point detection algorithms. It is multiply robust in the sense that it takes multiple initial estimators as input and only requires one of them to be consistent. Moreover, by borrowing information over time and population, it allows us to detect weaker signals and has better convergence properties when compared to applying the clustering algorithm per time or the change point detection algorithm per subject. Empirically, we demonstrate the usefulness of our method through extensive simulations and a real data application.


title: Causal Discovery in Linear Latent Variable Models Subject to Measurement Error

abstract: We focus on causal discovery in the presence of measurement error in linear systems where the mixing matrix, i.e., the matrix indicating the independent exogenous noise terms pertaining to the observed variables, is identified up to permutation and scaling of the columns. We demonstrate a somewhat surprising connection between this problem and causal discovery in the presence of unobserved parentless causes, in the sense that there is a mapping, given by the mixing matrix, between the underlying models to be inferred in these problems. Consequently, any identifiability result based on the mixing matrix for one model translates to an identifiability result for the other model. We characterize to what extent the causal models can be identified under a two-part faithfulness assumption. Under only the first part of the assumption (corresponding to the conventional definition of faithfulness), the structure can be learned up to the causal ordering among an ordered grouping of the variables but not all the edges across the groups can be identified. We further show that if both parts of the faithfulness assumption are imposed, the structure can be learned up to a more refined ordered grouping. As a result of this refinement, for the latent variable model with unobserved parentless causes, the structure can be identified. Based on our theoretical results, we propose causal structure learning methods for both models, and evaluate their performance on synthetic data.

categories: cs.LG cs.AI cs.IT math.IT stat.ME stat.ML


title: Adaptive Data Depth via Multi-Armed Bandits

abstract: Data depth, introduced by Tukey (1975), is an important tool in data science, robust statistics, and computational geometry. One chief barrier to its broader practical utility is that many common measures of depth are computationally intensive, requiring on the order of $n^d$ operations to exactly compute the depth of a single point within a data set of $n$ points in $d$-dimensional space. Often however, we are not directly interested in the absolute depths of the points, but rather in their relative ordering. For example, we may want to find the most central point in a data set (a generalized median), or to identify and remove all outliers (points on the fringe of the data set with low depth). With this observation, we develop a novel and instance-adaptive algorithm for adaptive data depth computation by reducing the problem of exactly computing $n$ depths to an $n$-armed stochastic multi-armed bandit problem which we can efficiently solve. We focus our exposition on simplicial depth, developed by Liu (1990), which has emerged as a promising notion of depth due to its interpretability and asymptotic properties. We provide general instance-dependent theoretical guarantees for our proposed algorithms, which readily extend to many other common measures of data depth including majority depth, Oja depth, and likelihood depth. When specialized to the case where the gaps in the data follow a power law distribution with parameter $\alpha<2$, we show that we can reduce the complexity of identifying the deepest point in the data set (the simplicial median) from $O(n^d)$ to $\tilde{O}(n^{d-(d-1)\alpha/2})$, where $\tilde{O}$ suppresses logarithmic factors. We corroborate our theoretical results with numerical experiments on synthetic data, showing the practical utility of our proposed methods.


title: PANDORA project: photo-nuclear reactions below $A=60$

abstract: Photo-nuclear reactions of light nuclei below a mass of $A=60$ are studied experimentally and theoretically by the PANDORA (Photo-Absorption of Nuclei and Decay Observation for Reactions in Astrophysics) project. Two experimental methods, virtual-photon excitation by proton scattering and real-photo absorption by a high-brilliance gamma-ray beam produced by laser Compton scattering, will be applied to measure the photo-absorption cross sections and the decay branching ratio of each decay channel as a function of the photon energy. Several nuclear models, e.g. anti-symmetrized molecular dynamics, mean-field type models, a large-scale shell model, and ab initio models, will be employed to predict the photo-nuclear reactions. The uncertainty in the model predictions will be evaluated from the discrepancies between the model predictions and the experimental data. The data and the predictions will be implemented in a general reaction calculation code TALYS . The results will be applied to the simulation of the photo-disintegration process of ultra-high-energy cosmic rays in inter-galactic propagation.


title: Theta series of ternary quadratic lattice cosets

abstract: In this paper, we consider the decomposition of theta series for lattice cosets of ternary lattices. We show that the natural decomposition into an Eisenstein series, a unary theta function, and a cuspidal form which is orthogonal to unary theta functions correspond to the theta series for the genus, the deficiency of the theta series for the spinor genus from that of the genus, and the deficiency of the theta series for the class from that of the spinor genus, respectively. These three pieces are hence invariants of the genus, spinor genus, and class, respectively, extending known results for lattices and verifying a conjecture of the first author and Haensch. We furthermore extend the definition of $p$-neighbors to include lattice cosets and construct an algorithm to compute respresentatives for the classes in the genus or spinor genus via the $p$-neighborhoods.


title: Unsupervised Domain Adaptation for Sparse Retrieval by Filling Vocabulary and Word Frequency Gaps

abstract: IR models using a pretrained language model significantly outperform lexical approaches like BM25. In particular, SPLADE, which encodes texts to sparse vectors, is an effective model for practical use because it shows robustness to out-of-domain datasets. However, SPLADE still struggles with exact matching of low-frequency words in training data. In addition, domain shifts in vocabulary and word frequencies deteriorate the IR performance of SPLADE. Because supervision data are scarce in the target domain, addressing the domain shifts without supervision data is necessary. This paper proposes an unsupervised domain adaptation method by filling vocabulary and word-frequency gaps. First, we expand a vocabulary and execute continual pretraining with a masked language model on a corpus of the target domain. Then, we multiply SPLADE-encoded sparse vectors by inverse document frequency weights to consider the importance of documents with lowfrequency words. We conducted experiments using our method on datasets with a large vocabulary gap from a source domain. We show that our method outperforms the present stateof-the-art domain adaptation method. In addition, our method achieves state-of-the-art results, combined with BM25.


title: $BT^2$: Backward-compatible Training with Basis Transformation

abstract: Modern retrieval system often requires recomputing the representation of every piece of data in the gallery when updating to a better representation model. This process is known as backfilling and can be especially costly in the real world where the gallery often contains billions of samples. Recently, researchers have proposed the idea of Backward Compatible Training (BCT) where the new representation model can be trained with an auxiliary loss to make it backward compatible with the old representation. In this way, the new representation can be directly compared with the old representation, in principle avoiding the need for any backfilling. However, followup work shows that there is an inherent tradeoff where a backward compatible representation model cannot simultaneously maintain the performance of the new model itself. This paper reports our ``not-so-surprising'' finding that adding extra dimensions to the representation can help here. However, we also found that naively increasing the dimension of the representation did not work. To deal with this, we propose Backward-compatible Training with a novel Basis Transformation ($BT^2$). A basis transformation (BT) is basically a learnable set of parameters that applies an orthonormal transformation. Such a transformation possesses an important property whereby the original information contained in its input is retained in its output. We show in this paper how a BT can be utilized to add only the necessary amount of additional dimensions. We empirically verify the advantage of $BT^2$ over other state-of-the-art methods in a wide range of settings. We then further extend $BT^2$ to other challenging yet more practical settings, including significant change in model architecture (CNN to Transformers), modality change, and even a series of updates in the model architecture mimicking the evolution of deep learning models.


title: Robust Unstructured Knowledge Access in Conversational Dialogue with ASR Errors

abstract: Performance of spoken language understanding (SLU) can be degraded with automatic speech recognition (ASR) errors. We propose a novel approach to improve SLU robustness by randomly corrupting clean training text with an ASR error simulator, followed by self-correcting the errors and minimizing the target classification loss in a joint manner. In the proposed error simulator, we leverage confusion networks generated from an ASR decoder without human transcriptions to generate a variety of error patterns for model training. We evaluate our approach on the DSTC10 challenge targeted for knowledge-grounded task-oriented conversational dialogues with ASR errors. Experimental results show the effectiveness of our proposed approach, boosting the knowledge-seeking turn detection (KTD) F1 significantly from 0.9433 to 0.9904. Knowledge cluster classification is boosted from 0.7924 to 0.9333 in Recall@1. After knowledge document re-ranking, our approach shows significant improvement in all knowledge selection metrics, from 0.7358 to 0.7806 in Recall@1, from 0.8301 to 0.9333 in Recall@5, and from 0.7798 to 0.8460 in MRR@5 on the test set. In the recent DSTC10 evaluation, our approach demonstrates significant improvement in knowledge selection, boosting Recall@1 from 0.495 to 0.7144 compared to the official baseline. Our source code is released in GitHub https://github.com/yctam/dstc10_track2_task2.git.


title: Dynamic Interpretable Change Point Detection

abstract: Identifying change points (CPs) in a time series is crucial to guide better decision making across various fields like finance and healthcare and facilitating timely responses to potential risks or opportunities. Existing Change Point Detection (CPD) methods have a limitation in tracking changes in the joint distribution of multidimensional features. In addition, they fail to generalize effectively within the same time series as different types of CPs may require different detection methods. As the volume of multidimensional time series continues to grow, capturing various types of complex CPs such as changes in the correlation structure of the time-series features has become essential. To overcome the limitations of existing methods, we propose TiVaCPD, an approach that uses a Time-Varying Graphical Lasso (TVGL) to identify changes in correlation patterns between multidimensional features over time, and combines that with an aggregate Kernel Maximum Mean Discrepancy (MMD) test to identify changes in the underlying statistical distributions of dynamic time windows with varying length. The MMD and TVGL scores are combined using a novel ensemble method based on similarity measures leveraging the power of both statistical tests. We evaluate the performance of TiVaCPD in identifying and characterizing various types of CPs and show that our method outperforms current state-of-the-art methods in real-world CPD datasets. We further demonstrate that TiVaCPD scores characterize the type of CPs and facilitate interpretation of change dynamics, offering insights into real-life applications.


title: Splitting expands the application range of Vision Transformer -- variable Vision Transformer (vViT)

abstract: Vision Transformer (ViT) has achieved outstanding results in computer vision. Although there are many Transformer-based architectures derived from the original ViT, the dimension of patches are often the same with each other. This disadvantage leads to a limited application range in the medical field because in the medical field, datasets whose dimension is different from each other; e.g. medical image, patients' personal information, laboratory test and so on. To overcome this limitation, we develop a new derived type of ViT termed variable Vision Transformer (vViT). The aim of this study is to introduce vViT and to apply vViT to radiomics using T1 weighted magnetic resonance image (MRI) of glioma. In the prediction of 365 days of survival among glioma patients using radiomics,vViT achieved 0.83, 0.82, 0.81, and 0.76 in sensitivity, specificity, accuracy, and AUC-ROC, respectively. vViT has the potential to handle different types of medical information at once.


title: The Radul cocycle, the Chern-Connes character, and manifolds with conical singularities

abstract: This short note establishes a relationship between a generalized version of the Radul residue cocycle introduced in former works of the author and the Connes-Moscovici residue cocycle, and discusses the applicability of such a formula to manifolds with conical singularities, where zeta functions of Fuchs-type pseudodifferential operators may exhibit double or triple poles.


title: Reinforcement Learning with Stepwise Fairness Constraints

abstract: AI methods are used in societally important settings, ranging from credit to employment to housing, and it is crucial to provide fairness in regard to algorithmic decision making. Moreover, many settings are dynamic, with populations responding to sequential decision policies. We introduce the study of reinforcement learning (RL) with stepwise fairness constraints, requiring group fairness at each time step. Our focus is on tabular episodic RL, and we provide learning algorithms with strong theoretical guarantees in regard to policy optimality and fairness violation. Our framework provides useful tools to study the impact of fairness constraints in sequential settings and brings up new challenges in RL.


title: Computing palindromes on a trie in linear time

abstract: A trie $\mathcal{T}$ is a rooted tree such that each edge is labeled by a single character from the alphabet, and the labels of out-going edges from the same node are mutually distinct. Given a trie $\mathcal{T}$ with $n$ edges, we show how to compute all distinct palindromes and all maximal palindromes on $\mathcal{T}$ in $O(n)$ time, in the case of integer alphabets of size polynomial in $n$. This improves the state-of-the-art $O(n \log h)$-time algorithms by Funakoshi et al. [PCS 2019], where $h$ is the height of $\mathcal{T}$. Using our new algorithms, the eertree with suffix links for a given trie $\mathcal{T}$ can readily be obtained in $O(n)$ time. Further, our trie-based $O(n)$-space data structure allows us to report all distinct palindromes and maximal palindromes in a query string represented in the trie $\mathcal{T}$, in output optimal time. This is an improvement over an existing (na\"ive) solution that precomputes and stores all distinct palindromes and maximal palindromes for each and every string in the trie $\mathcal{T}$ separately, using a total $O(n^2)$ preprocessing time and space, and reports them in output optimal time upon query.


title: Algebra cochains, the bivariant JLO cocycle and the Mathai-Quillen form

abstract: This is a first investigation by the author of the similarity between Quillen's superconnection formalism, his constructions of (periodic) cyclic cocycles via algebra cochains on a bar construction, and Kasparov bimodules for KK-theory. In this article, we do so by deriving a slight extension of the Mathai-Quillen Thom form via a bivariant JLO cocycle. The main idea (which is in fact not really new) is that KK-cycles should be thought of as superconnection forms; these methods will be applied to other contexts elsewhere.


title: A Simple Algorithm for Online Decision Making

abstract: Motivated by recent progress on online linear programming (OLP), we study the online decision making problem (ODMP) as a natural generalization of OLP. In ODMP, there exists a single decision maker who makes a series of decisions spread out over a total of $T$ time stages. At each time stage, the decision maker makes a decision based on information obtained up to that point without seeing into the future. The task of the decision maker is to maximize the accumulated reward while overall meeting some predetermined $m$-dimensional long-term goal (linking) constraints. ODMP significantly broadens the modeling framework of OLP by allowing more general feasible regions (for local and goal constraints) potentially involving both discreteness and nonlinearity in each local decision making problem. We propose a Fenchel dual-based online algorithm for ODMP. At each time stage, the proposed algorithm requires solving a potentially nonconvex optimization problem over the local feasible set and a convex optimization problem over the goal set. Under the uniform random permutation model, we show that our algorithm achieves $O(\sqrt{mT})$ constraint violation deterministically in meeting the long-term goals, and $O(\sqrt{m\log m}\sqrt{T})$ competitive difference in expected reward with respect to the optimal offline decisions. We also extend our results to the grouped random permutation model.


title: A transverse Chern character for compact Lie group actions

abstract: This article aims to explore new perspectives offered by Kasparov's recent work on transverse index theory in the context of actions of compact Lie groups on a manifold, and hints at potential connections to the work of Berline-Vergne and Paradan-Vergne.


title: A polynomial-time classical algorithm for noisy random circuit sampling

abstract: We give a polynomial time classical algorithm for sampling from the output distribution of a noisy random quantum circuit in the regime of anti-concentration to within inverse polynomial total variation distance. This gives strong evidence that, in the presence of a constant rate of noise per gate, random circuit sampling (RCS) cannot be the basis of a scalable experimental violation of the extended Church-Turing thesis. Our algorithm is not practical in its current form, and does not address finite-size RCS based quantum supremacy experiments.


title: When Less Is More: Consequence-Finding in a Weak Theory of Arithmetic

abstract: This paper presents a theory of non-linear integer/real arithmetic and algorithms for reasoning about this theory. The theory can be conceived as an extension of linear integer/real arithmetic with a weakly-axiomatized multiplication symbol, which retains many of the desirable algorithmic properties of linear arithmetic. In particular, we show that the conjunctive fragment of the theory can be effectively manipulated (analogously to the usual operations on convex polyhedra, the conjunctive fragment of linear arithmetic). As a result, we can solve the following consequence-finding problem: given a ground formula F, find the strongest conjunctive formula that is entailed by F. As an application of consequence-finding, we give a loop invariant generation algorithm that is monotone with respect to the theory and (in a sense) complete. Experiments show that the invariants generated from the consequences are effective for proving safety properties of programs that require non-linear reasoning.


title: Progress and summary of reinforcement learning on energy management of MPS-EV

abstract: The high emission and low energy efficiency caused by internal combustion engines (ICE) have become unacceptable under environmental regulations and the energy crisis. As a promising alternative solution, multi-power source electric vehicles (MPS-EVs) introduce different clean energy systems to improve powertrain efficiency. The energy management strategy (EMS) is a critical technology for MPS-EVs to maximize efficiency, fuel economy, and range. Reinforcement learning (RL) has become an effective methodology for the development of EMS. RL has received continuous attention and research, but there is still a lack of systematic analysis of the design elements of RL-based EMS. To this end, this paper presents an in-depth analysis of the current research on RL-based EMS (RL-EMS) and summarizes the design elements of RL-based EMS. This paper first summarizes the previous applications of RL in EMS from five aspects: algorithm, perception scheme, decision scheme, reward function, and innovative training method. The contribution of advanced algorithms to the training effect is shown, the perception and control schemes in the literature are analyzed in detail, different reward function settings are classified, and innovative training methods with their roles are elaborated. Finally, by comparing the development routes of RL and RL-EMS, this paper identifies the gap between advanced RL solutions and existing RL-EMS. Finally, this paper suggests potential development directions for implementing advanced artificial intelligence (AI) solutions in EMS.


title: The free algebra in R

abstract: The free algebra is an interesting and useful algebraic object. Here I introduce "freealg", an R package which furnishes computational support for free algebras. The package uses the standard template library's "map" class for efficiency, which uses the fact that the order of the terms is algebraically immaterial. The package follows "disordR" discipline. I demonstrate some properties of free algebra using the package, and showcase package idiom. The package is available on CRAN at https://CRAN.R-project.org/package=freealg.


title: Calcul du cocycle JLO pour l'op\'erateur de Dirac classique

abstract: We provide a short account of a classical folklore result in connection to the local index theorem, which identifies the short-time limit of the JLO cocycle of the Dirac operator to de Rham current obtained by the cap product of the fundamental class with the \^A-genus.


title: State-Insensitive Trapping of Alkaline-Earth Atoms in a Nanofiber-Based Optical Dipole Trap

abstract: Neutral atoms trapped in the evanescent optical potentials of nanotapered optical fibers are a promising platform for developing quantum technologies and exploring fundamental science, such as quantum networks and quantum electrodynamics. Building on the successful advancements with trapped alkali atoms, here we demonstrate a state-insensitive optical dipole trap for strontium-88, an alkaline-earth atom, using the evanescent fields of a nanotapered optical fiber. Leveraging the low laser-cooling temperatures of $\sim\!\!1~\mu$K readily achievable with strontium, we demonstrate trapping in record low trap depths corresponding to $\sim\!\!3~\mu$K. Further, employing a double magic wavelength trapping scheme, we realize state-insensitive trapping on the kilohertz-wide $5s^{2}\;^{1}\!S_{0}-5s5p\;^{3}\!P_{1,|m|=1}$ cooling transition, which we verify by performing near-surface high-resolution spectroscopy of the atomic transition. This allows us to experimentally find and verify the state insensitivity of the trap nearby a theoretically predicted magic wavelength of 435.827(25) nm. Given the non-magnetic ground state and low collisional scattering length of strontium-88, this work also lays the foundation for developing versatile and robust matter-wave atomtronic circuits over nanophotonic waveguides.


title: ABC: Adversarial Behavioral Cloning for Offline Mode-Seeking Imitation Learning

abstract: Given a dataset of expert agent interactions with an environment of interest, a viable method to extract an effective agent policy is to estimate the maximum likelihood policy indicated by this data. This approach is commonly referred to as behavioral cloning (BC). In this work, we describe a key disadvantage of BC that arises due to the maximum likelihood objective function; namely that BC is mean-seeking with respect to the state-conditional expert action distribution when the learner's policy is represented with a Gaussian. To address this issue, we introduce a modified version of BC, Adversarial Behavioral Cloning (ABC), that exhibits mode-seeking behavior by incorporating elements of GAN (generative adversarial network) training. We evaluate ABC on toy domains and a domain based on Hopper from the DeepMind Control suite, and show that it outperforms standard BC by being mode-seeking in nature.


title: Improving the science process skills of physics education students by using guided inquiry practicum

abstract: This research investigate that science process skills significantly improve after doing some practicum activities. The research population are fifth semester physics education students and the research sample are fifth semester physics education students who was doing electricity and magnetism experiment C class course. We used two questionnaires, the first one is given to the students after doing three experiment activities and the second one is given after doing six experiment activities. This research is quantitative research using paired sample t test analysis that compared the first questionnaires score and the second questionnaires by using SPSS software. The result indicates that the number of practicum activities is able to improve the science process skills significantly.


title: Lattice local integrable regularization of the Sine-Gordon model

abstract: We study the local lattice integrable regularization of the Sine-Gordon model written down in terms of the lattice Bose-operators. We show that the local spin Hamiltonian obtained from the six-vertex model with alternating inhomogeneities in fact leads to the Sine-Gordon model in the low-energy limit. We show that the Bethe Ansatz results for this model lead to the correct general relations for different critical exponents of the coupling constant.


title: A Toy Model for Low Energy Nuclear Fusion

abstract: We study the fusion of a proton with a nucleus with the emission of two photons at low incident energy of the order of eV or smaller. We use a step model for the repulsive potential between proton and the nuclei. We consider the reaction both in free space and inside a medium. We make a simple model for the medium by assuming a hard wall potential beyond a certain length scale. This essentially leads to discretization of the energy spectrum which is expected inside a medium and is seen both for a crystalline lattice structure and for amorphous materials. We use second order perturbation theory to compute the transition rate. We find that the rate in free space is very small. However in medium, the rate may be substantial. Hence, we conclude that nuclear fusion reactions may take place at low energies at observable rates.


title: SOTIF Entropy: Online SOTIF Risk Quantification and Mitigation for Autonomous Driving

abstract: Autonomous driving confronts great challenges in complex traffic scenarios, where the risk of Safety of the Intended Functionality (SOTIF) can be triggered by the dynamic operational environment and system insufficiencies. The SOTIF risk is reflected not only intuitively in the collision risk with objects outside the autonomous vehicles (AVs), but also inherently in the performance limitation risk of the implemented algorithms themselves. How to minimize the SOTIF risk for autonomous driving is currently a critical, difficult, and unresolved issue. Therefore, this paper proposes the "Self-Surveillance and Self-Adaption System" as a systematic approach to online minimize the SOTIF risk, which aims to provide a systematic solution for monitoring, quantification, and mitigation of inherent and external risks. The core of this system is the risk monitoring of the implemented artificial intelligence algorithms within the AV. As a demonstration of the Self-Surveillance and Self-Adaption System, the risk monitoring of the perception algorithm, i.e., YOLOv5 is highlighted. Moreover, the inherent perception algorithm risk and external collision risk are jointly quantified via SOTIF entropy, which is then propagated downstream to the decision-making module and mitigated. Finally, several challenging scenarios are demonstrated, and the Hardware-in-the-Loop experiments are conducted to verify the efficiency and effectiveness of the system. The results demonstrate that the Self-Surveillance and Self-Adaption System enables dependable online monitoring, quantification, and mitigation of SOTIF risk in real-time critical traffic environments.


title: Numerical analysis of Givens rotation

abstract: Generating 2-by-2 unitary matrices in floating-precision arithmetic is a delicate task. One way to reduce the accumulation error is to use less floating-point operations to compute each of the entries in the 2-by-2 unitary matrix. This paper shows an algorithm that reduces the number of operations to compute the entries of a Givens rotation. Overall, the new algorithm has more operations in total when compared to algorithms in different releases of LAPACK, but less operations per entry. Numerical tests show that the new algorithm is more accurate on average.


title: An Incremental Phase Mapping Approach for X-ray Diffraction Patterns using Binary Peak Representations

abstract: Despite the huge advancement in knowledge discovery and data mining techniques, the X-ray diffraction (XRD) analysis process has mostly remained untouched and still involves manual investigation, comparison, and verification. Due to the large volume of XRD samples from high-throughput XRD experiments, it has become impossible for domain scientists to process them manually. Recently, they have started leveraging standard clustering techniques, to reduce the XRD pattern representations requiring manual efforts for labeling and verification. Nevertheless, these standard clustering techniques do not handle problem-specific aspects such as peak shifting, adjacent peaks, background noise, and mixed phases; hence, resulting in incorrect composition-phase diagrams that complicate further steps. Here, we leverage data mining techniques along with domain expertise to handle these issues. In this paper, we introduce an incremental phase mapping approach based on binary peak representations using a new threshold based fuzzy dissimilarity measure. The proposed approach first applies an incremental phase computation algorithm on discrete binary peak representation of XRD samples, followed by hierarchical clustering or manual merging of similar pure phases to obtain the final composition-phase diagram. We evaluate our method on the composition space of two ternary alloy systems- Co-Ni-Ta and Co-Ti-Ta. Our results are verified by domain scientists and closely resembles the manually computed ground-truth composition-phase diagrams. The proposed approach takes us closer towards achieving the goal of complete end-to-end automated XRD analysis.


title: A multivariate functional-data mixture model for spatio-temporal data: inference and cokriging

abstract: In this paper, we introduce a model for multivariate, spatio-temporal functional data. Specifically, this work proposes a mixture model that is used to perform spatio-temporal prediction (cokriging) when both the response and the additional covariates are functional data. The estimation of such models in the context of expansive data poses many methodological and computational challenges. We propose a new Monte Carlo Expectation Maximization algorithm based on importance sampling to estimate model parameters. We validate our methodology using simulation studies and provide a comparison to previously proposed functional clustering methodologies. To tackle computational challenges, we describe a variety of advances that enable application on large spatio-temporal datasets. The methodology is applied on Argo oceanographic data in the Southern Ocean to predict oxygen concentration, which is critical to ocean biodiversity and reflects fundamental aspects of the carbon cycle. Our model and implementation comprehensively provide oxygen predictions and their estimated uncertainty as well as recover established oceanographic fronts.


title: COV19IR : COVID-19 Domain Literature Information Retrieval

abstract: Increasing number of COVID-19 research literatures cause new challenges in effective literature screening and COVID-19 domain knowledge aware Information Retrieval. To tackle the challenges, we demonstrate two tasks along withsolutions, COVID-19 literature retrieval, and question answering. COVID-19 literature retrieval task screens matching COVID-19 literature documents for textual user query, and COVID-19 question answering task predicts proper text fragments from text corpus as the answer of specific COVID-19 related questions. Based on transformer neural network, we provided solutions to implement the tasks on CORD-19 dataset, we display some examples to show the effectiveness of our proposed solutions.


title: Visualization of Strain-Induced Landau Levels in a Graphene - Black Phosphorus Heterostructure

abstract: Strain-induced pseudo magnetic fields offer the possibility of realizing zero magnetic field Quantum Hall effect in graphene, possibly up to room temperature, representing a promising avenue for lossless charge transport applications. Strain engineering on graphene has been achieved via random nanobubbles or artificial nanostructures on the substrate, but the highly localized and non-uniform pseudomagnetic fields can make spectroscopic probes of electronic structure difficult. Heterostructure engineering offers an alternative approach: By stacking graphene on top of another van der Waals material with large lattice mismatch at a desired twist angle, it is possible to generate large strain-induced pseudo magnetic fields uniformly over the entire heterostructure. Here, we report using nano-angle resolved photoemission spectroscopy (nano-ARPES) to probe the electronic bandstructure of a graphene/black phosphorus heterostructure (G/BP). By directly measuring the iso-energy contours of graphene and black phosphorus we determine a twist angle of 20-degrees in our heterostructure. High-resolution nano-ARPES of the graphene bands near the Fermi level reveals the emergence of flat bands located within the Dirac cone. The spacing of the flat bands is consistent with Landau level formation in graphene, and corresponds to a pseudo-field of 11.36 T. Our work provides a new way to study quantum Hall phases induced by strain in 2D materials and heterostructures.


title: A Classification of Two-dimensional Endo-commutative Algebras over F_2

abstract: We introduce a new class of algebras called endo-commutative algebras in which the square mapping preserves multiplication, and provide a complete classification of endo-commutative algebras of dimension 2 over the field F_2 of two elements. We list all multiplication tables of the algebras up to isomorphism. This clarifies the difference between commutativity and endo-commutativity of algebras.


title: Single-Photon Signal Sideband Detection for High-Power Michelson Interferometers

abstract: The Michelson interferometer is a cornerstone of experimental physics. Its applications range from providing first impressions of wave interference in educational settings to probing spacetime at minuscule precision scales. Interferometer precision provides a unique view of the fundamental medium of matter and energy, enabling tests for new physics as well as searches for the gravitational wave signatures of distant astrophysical events. Optical interferometers are typically operated by continuously measuring the power at their output port. Signal perturbations then create sideband fields, forming a beat-note with the fringe light that modulates that power. When operated at a nearly-dark destructive interference fringe, this readout is a form of homodyne detection, with an imprecision set by a ``standard quantum limit'' attributed to shot noise from quantum vacuum fluctuations. The sideband signal fields carry energy which can, alternatively, be directly observed as photons distinct from the source laser. Without signal energy, vacuum does not form sidebands and cannot spuriously create photon counts or shot noise. Thus, counting can offer improved statistics when searching for weak signals when classical backgrounds are below the standard quantum limit. Here, photon counting statistics are described for optical interferometry, relating the two forms of measurement and showing cases where counting greatly outperforms homodyne readout, even with squeezed state quantum enhancement. The most immediate application for photon counting is improving searches of stochastic signals, such as from quantum gravity or from new particle fields. The advantages of counting may extend to wider applications, such as gravitational wave detectors, and the concept of Fisher-information representative spectral density is introduced to motivate further study.


title: On the topological surface states of the intrinsic magnetic topological insulator Mn-Bi-Te family

abstract: We review recent progress in the electronic structure study of intrinsic magnetic topological insulators (MnBi$_2$Te$_4$)(Bi$_2$Te$_3$)$_n$ ($n=0,1,2,3$) family. Specifically, we focus on the ubiquitously (nearly) gapless behavior of the topological surface state Dirac cone observed by photoemission spectroscopy, even though a large Dirac gap is expected because of surface ferromagnetic order. The dichotomy between experiment and theory concerning this gap behavior is perhaps the most critical and puzzling question in this frontier. We discuss various proposals accounting for the lack of magnetic effect on the topological surface state Dirac cone, which are mainly categorized into two pictures, magnetic reconfiguration, and topological surface state redistribution. Band engineering towards opening a magnetic gap of topological surface states provides great opportunities to realize quantized topological transport and axion electrodynamics at higher temperatures.


title: Thermally Induced Structural Evolution and Nanoscale Interfacial Dynamics in Layered Metal Chalcogenides

abstract: Layered chalcogenides including Bi2Te3, Sb2Te3, Bi-Sb-Te ternary alloys and heterostructures are known as great thermoelectric, topological insulators and recently highlighted as plasmonic building blocks beyond noble metals. Here, we conduct a joint in situ transmission electron microscopy (in situ TEM) and density functional theory (DFT) calculations to investigate the temperature dependent nanoscale dynamics, interfacial properties and further identifying the role of native defects and edge configurations in anisotropic sublimations of Bi2Te3-Sb2Te3 in-plane heterostructure and Sb2-xBixTe3 alloy. Structural dynamics including edge evolution, formation, expansion, and coalescence of thermally induced polygonal nanopores are reported. The nanopores appear to be initiated by preferential dissociation of chalcogenide species (Te) from the center, heterointerface and edges in the heterostructure and only from the outer edges in the alloy counterpart. This results in a reduced thermal stability and significantly different sublimation pathways of the heterostructure. Furthermore, triangular and quasi hexagonal configurations are observed to be the dominant nanopores configurations in the heterostructure. Additionally, our DFT calculations provide a mechanistic understanding on the role of native defects and edge formation energies, revealing the antisite defects TeBi to be the dominant native defect in a Te-rich condition and playing a key role on the defects assisted sublimation. These findings significantly impact our understanding of controlling the nanoscale sublimation dynamics and can ultimately assist us in designing tunable low-dimensional chalcogenides.


title: Dynamic Sensor Placement Based on Graph Sampling Theory

abstract: In this paper, we consider a sensor placement problem where sensors can move within a network over time. Sensor placement problem aims to select K sensor positions from N candidates where K < N. Most existing methods assume that sensor positions are static, i.e., they do not move, however, many mobile sensors like drones, robots, and vehicles can change their positions over time. Moreover, underlying measurement conditions could also be changed, which are difficult to cover with statically placed sensors. We tackle the problem by allowing the sensors to change their positions in their neighbors on the network. We sequentially learn the dictionary from a pool of observed signals on the network based on graph sampling theory. Using the learned dictionary, we dynamically determine the sensor positions such that the non-observed signals on the network can be best recovered from the observations. In experiments, we validate the effectiveness of the proposed method via the mean squared error (MSE) of the reconstructed signals. The proposed dynamic sensor placement outperforms the existing static ones for both synthetic and real data.


title: Generating counterfactual explanations of tumor spatial proteomes to discover effective strategies for enhancing immune infiltration

abstract: The tumor microenvironment (TME) significantly impacts cancer prognosis due to its immune composition. While therapies for altering the immune composition, including immunotherapies, have shown exciting results for treating hematological cancers, they are less effective for immunologically-cold, solid tumors. Spatial omics technologies capture the spatial organization of the TME with unprecedented molecular detail, revealing the relationship between immune cell localization and molecular signals. Here, we formulate T-cell infiltration prediction as a self-supervised machine learning problem and develop a counterfactual optimization strategy that leverages large scale spatial omics profiles of patient tumors to design tumor perturbations predicted to boost T-cell infiltration. A convolutional neural network predicts T-cell distribution based on signaling molecules in the TME provided by imaging mass cytometry. Gradient-based counterfactual generation, then, computes perturbations predicted to boost T-cell abundance. We apply our framework to melanoma, colorectal cancer liver metastases, and breast tumor data, discovering combinatorial perturbations predicted to support T-cell infiltration across tens to hundreds of patients. This work presents a paradigm for counterfactual-based prediction and design of cancer therapeutics using spatial omics data.


title: Space Weather Observations, Modeling, and Alerts in Support of Human Exploration of Mars

abstract: Space weather observations and modeling at Mars have begun but they must be significantly increased to support the future of Human Exploration on the Red Planet. A comprehensive space weather understanding of a planet without a global magnetosphere and a thin atmosphere is very different from our situation at Earth so there is substantial fundamental research remaining. It is expected that the development of suitable models will lead to a comprehensive operational Mars space weather alert (MSWA) system that would provide rapid dissemination of information to Earth controllers, astronauts in transit, and those in the exploration zone (EZ) on the surface by producing alerts that are delivered rapidly and are actionable. To illustrate the importance of such a system, we use a magnetohydrodynamic code to model an extreme Carrington-type coronal mass ejection (CME) event at Mars. The results show a significant induced surface field of nearly 3000 nT on the dayside that could radically affect unprotected electrical systems that would dramatically impact human survival on Mars. Other associated problems include coronal mass ejection (CME) shock-driven acceleration of solar energetic particles producing large doses of ionizing radiation at the Martian surface. In summary, along with working more closely with international partners, the next Heliophysics Decadal Survey must include a new initiative to meet expected demands for space weather forecasting in support of humans living and working on the surface of Mars. It will require significant effort to coordinate NASA and the international community contributions.


title: Integrated Sensing, Computation, and Communication: System Framework and Performance Optimization

abstract: Integrated sensing, computation, and communication (ISCC) has been recently considered as a promising technique for beyond 5G systems. In ISCC systems, the competition for communication and computation resources between sensing tasks for ambient intelligence and computation tasks from mobile devices becomes an increasingly challenging issue. To address it, we first propose an efficient sensing framework with a novel action detection module. In this module, a threshold is used for detecting whether the sensing target is static and thus the overhead can be reduced. Subsequently, we mathematically analyze the sensing performance of the proposed framework and theoretically prove its effectiveness with the help of the sampling theorem. Based on sensing performance models, we formulate a sensing performance maximization problem while guaranteeing the quality-of-service (QoS) requirements of tasks. To solve it, we propose an optimal resource allocation strategy, in which the minimum resource is allocated to computation tasks, and the rest is devoted to the sensing task. Besides, a threshold selection policy is derived and the results further demonstrate the necessity of the proposed sensing framework. Finally, a real-world test of action recognition tasks based on USRP B210 is conducted to verify the sensing performance analysis. Extensive experiments demonstrate the performance improvement of our proposal by comparing it with some benchmark schemes.


title: A Dynamic Graph Interactive Framework with Label-Semantic Injection for Spoken Language Understanding

abstract: Multi-intent detection and slot filling joint models are gaining increasing traction since they are closer to complicated real-world scenarios. However, existing approaches (1) focus on identifying implicit correlations between utterances and one-hot encoded labels in both tasks while ignoring explicit label characteristics; (2) directly incorporate multi-intent information for each token, which could lead to incorrect slot prediction due to the introduction of irrelevant intent. In this paper, we propose a framework termed DGIF, which first leverages the semantic information of labels to give the model additional signals and enriched priors. Then, a multi-grain interactive graph is constructed to model correlations between intents and slots. Specifically, we propose a novel approach to construct the interactive graph based on the injection of label semantics, which can automatically update the graph to better alleviate error propagation. Experimental results show that our framework significantly outperforms existing approaches, obtaining a relative improvement of 13.7% over the previous best model on the MixATIS dataset in overall accuracy.


title: Comparing Two Counting Methods for Estimating the Probabilities of Strings

abstract: There are two methods for counting the number of occurrences of a string in another large string. One is to count the number of places where the string is found. The other is to determine how many pieces of string can be extracted without overlapping. The difference between the two becomes apparent when the string is part of a periodic pattern. This research reports that the difference is significant in estimating the occurrence probability of a pattern. In this study, the strings used in the experiments are approximated from time-series data. The task involves classifying strings by estimating the probability or computing the information quantity. First, the frequencies of all substrings of a string are computed. Each counting method may sometimes produce different frequencies for an identical string. Second, the probability of the most probable segmentation is selected. The probability of the string is the product of all probabilities of substrings in the selected segmentation. The classification results demonstrate that the difference in counting methods is statistically significant, and that the method without overlapping is better.


title: Directed Steiner path packing and directed path connectivity

abstract: For a digraph $D=(V(D), A(D))$, and a set $S\subseteq V(D)$ with $r\in S$ and $|S|\geq 2$, a directed $(S, r)$-Steiner path or, simply, an $(S, r)$-path is a directed path $P$ started at $r$ with $S\subseteq V(P)$. Two $(S, r)$-paths are said to be arc-disjoint if they have no common arc. Two arc-disjoint $(S, r)$-paths are said to be internally disjoint if the set of common vertices of them is exactly $S$. Let $\kappa^p_{S,r}(D)$ (resp. $\lambda^p_{S,r}(D)$) be the maximum number of internally disjoint (resp. arc-disjoint) $(S, r)$-paths in $D$. The directed path $k$-connectivity of $D$ is defined as $$\kappa^p_k(D)= \min \{\kappa^p_{S,r}(D)\mid S\subseteq V(D), |S|=k, r\in S\}.$$ Similarly, the directed path $k$-arc-connectivity of $D$ is defined as $$\lambda^p_k(D)= \min \{\lambda^p_{S,r}(D)\mid S\subseteq V(D), |S|=k, r\in S\}.$$ The directed path $k$-connectivity and directed path $k$-arc-connectivity are also called directed path connectivity which extends the path connectivity on undirected graphs to directed graphs and could be seen as a generalization of classical connectivity of digraphs. In this paper, we obtain complexity results for $\kappa^p_{S,r}(D)$ on Eulerian digraphs and symmetric digraphs, and $\lambda^p_{S,r}(D)$ on general digraphs. We also give bounds for the parameters $\kappa^p_k(D)$ and $\lambda^p_k(D)$.


title: Domain-decomposed Bayesian inversion based on local Karhunen-Lo\`{e}ve expansions

abstract: In many Bayesian inverse problems the goal is to recover a spatially varying random field. Such problems are often computationally challenging especially when the forward model is governed by complex partial differential equations (PDEs). The challenge is particularly severe when the spatial domain is large and the unknown random field needs to be represented by a high-dimensional parameter. In this paper, we present a domain-decomposed method to attack the dimensionality issue and the method decomposes the spatial domain and the parameter domain simultaneously. On each subdomain, a local Karhunen-Lo`eve (KL) expansion is constructed, and a local inversion problem is solved independently in a parallel manner, and more importantly, in a lower-dimensional space. After local posterior samples are generated through conducting Markov chain Monte Carlo (MCMC) simulations on subdomains, a novel projection procedure is developed to effectively reconstruct the global field. In addition, the domain decomposition interface conditions are dealt with an adaptive Gaussian process-based fitting strategy. Numerical examples are provided to demonstrate the performance of the proposed method.


title: Bootstraps for Dynamic Panel Threshold Models

abstract: This paper develops valid bootstrap inference methods for the dynamic panel threshold regression. For the first-differenced generalized method of moments (GMM) estimation for the dynamic short panel, we show that the standard nonparametric bootstrap is inconsistent. The inconsistency is due to an $n^{1/4}$-consistent non-normal asymptotic distribution for the threshold estimate when the parameter resides within the continuity region of the parameter space, which stems from the rank deficiency of the approximate Jacobian of the sample moment conditions on the continuity region. We propose a grid bootstrap to construct confidence sets for the threshold, a residual bootstrap to construct confidence intervals for the coefficients, and a bootstrap for testing continuity. They are shown to be valid under uncertain continuity. A set of Monte Carlo experiments demonstrate that the proposed bootstraps perform well in the finite samples and improve upon the asymptotic normal approximation even under a large jump at the threshold. An empirical application to firms' investment model illustrates our methods.


title: Impact of Radiation and Slip Conditions on MHD Flow of Nanofluid Past an Exponentially Stretched Surface

abstract: The current research establishes magnetohydrodynamics (MHD) boundary layer flow with heat and mass transfer of a nanofluid over an exponentially extending sheet embedded in a porous medium. During this exploration, nanoparticles, single-wall carbon nanotubes (SWCNTs) and multi-wall carbon nanotubes (MWCNTs) are recruited, while lamp fuel oil is being utilised as a base fluid for the diffusion of nano materials. The effects of warm radiation and an inclined magnetic field are included. In addition, rather than no-slip assumptions at the surface, velocity slides as well as thermal upsurge are incorporated in this study. Similarity transformations are implemented to adapt a set of partial differential equations into a system of non-linear ordinary differential equations. The bvp4c solver and Keller-box approach are employed to tackle nonlinear ordinary differential equations numerically. The significance of prominent parameters such as the Darcy Forchheimer model, magnetic field, radiation, suction, velocity slip, and temperature jump is visually probed and addressed in depth. In fact, the evolution of the coefficient of skin friction and percentage of heat shipping (Nusselt number) for both SWCNTs and MWCNTs is presented in tabular form. The temperature goes up as the magnetic parameter rises. Temperature has been seen to be decreased as the thermal slip parameter is improved. The results indicate that SWCNTs yield a higher coefficient of skin friction and speed of heat transformation than MWCNTs.


title: Recent studies on the super edge-magic deficiency of graphs

abstract: A graph $G$ is called edge-magic if there exists a bijective function $f:V\left(G\right) \cup E\left(G\right)\rightarrow \left\{1, 2, \ldots , \left\vert V\left( G\right) \right\vert +\left\vert E\left( G\right) \right\vert \right\}$ such that $f\left(u\right) + f\left(v\right) + f\left(uv\right)$ is a constant for each $uv\in E\left( G\right) $. Also, $G$ is said to be super edge-magic if $f\left(V \left(G\right)\right) =\left\{1, 2, \ldots , \left\vert V\left( G\right) \right\vert \right\}$. Furthermore, the super edge-magic deficiency $ \mu_{s}\left(G\right)$ of a graph $G$ is defined to be either the smallest nonnegative integer $n$ with the property that $G \cup nK_{1}$ is super edge-magic or $+ \infty$ if there exists no such integer $n$. In this paper, we introduce the parameter $l\left(n\right)$ as the minimum size of a graph $G$ of order $n$ for which all graphs of order $n$ and size at least $l\left(n\right)$ have $\mu_{s} \left( G \right)=+\infty $, and provide lower and upper bounds for $l\left(G\right)$. Imran, Baig, and Fe\u{n}ov\u{c}\'{i}kov\'{a} established that for integers $n$ with $n\equiv 0\pmod{4}$, $ \mu_{s}\left(D_{n}\right) \leq 3n/2-1$, where $D_{n}$ is the cartesian product of the cycle $C_{n}$ of order $n$ and the complete graph $K_{2}$ of order $2$. We improve this bound by showing that $ \mu_{s}\left(D_{n}\right) \leq n+1$ when $n \geq 4$ is even. Enomoto, Llad\'{o}, Nakamigawa, and Ringel posed the conjecture that every nontrivial tree is super edge-magic. We propose a new approach to attak this conjecture. This approach may also help to resolve another labeling conjecture on trees by Graham and Sloane.


title: Reduction theorem for support $\tau$-tilting modules over group algebras

abstract: In studying the structure of derived categories of module categories of group algebras or their blocks, it is fundamental to classify support $\tau$-tilting modules. Koshio and Kozakai showed that the structure of support $\tau$-tilting modules over blocks of finite groups can be reduced to that of their subgroups under suitable conditions. We show that the 'quotient reduction' is also valid under suitable conditions.


title: Hilbert Distillation for Cross-Dimensionality Networks

abstract: 3D convolutional neural networks have revealed superior performance in processing volumetric data such as video and medical imaging. However, the competitive performance by leveraging 3D networks results in huge computational costs, which are far beyond that of 2D networks. In this paper, we propose a novel Hilbert curve-based cross-dimensionality distillation approach that facilitates the knowledge of 3D networks to improve the performance of 2D networks. The proposed Hilbert Distillation (HD) method preserves the structural information via the Hilbert curve, which maps high-dimensional (>=2) representations to one-dimensional continuous space-filling curves. Since the distilled 2D networks are supervised by the curves converted from dimensionally heterogeneous 3D features, the 2D networks are given an informative view in terms of learning structural information embedded in well-trained high-dimensional representations. We further propose a Variable-length Hilbert Distillation (VHD) method to dynamically shorten the walking stride of the Hilbert curve in activation feature areas and lengthen the stride in context feature areas, forcing the 2D networks to pay more attention to learning from activation features. The proposed algorithm outperforms the current state-of-the-art distillation techniques adapted to cross-dimensionality distillation on two classification tasks. Moreover, the distilled 2D networks by the proposed method achieve competitive performance with the original 3D networks, indicating the lightweight distilled 2D networks could potentially be the substitution of cumbersome 3D networks in the real-world scenario.


title: Non-existence of a short algorithm for multiplication of $3\times3$ matrices with group $S_4\times S_3$, II

abstract: It is proved that there is no an algorithm for multiplication of $3\times3$ matrices of multiplicative length $\leq23$ that is invariant under a certain group isomorphic to $S_4\times S_3$. The proof makes use of description of the orbits of this group on decomposable tensors in the tensor cube $(M_3({\mathbb C}))^{\otimes3}$ which was obtained earlier.


title: AEDNet: Adaptive Edge-Deleting Network For Subgraph Matching

abstract: Subgraph matching is to find all subgraphs in a data graph that are isomorphic to an existing query graph. Subgraph matching is an NP-hard problem, yet has found its applications in many areas. Many learning-based methods have been proposed for graph matching, whereas few have been designed for subgraph matching. The subgraph matching problem is generally more challenging, mainly due to the different sizes between the two graphs, resulting in considerable large space of solutions. Also the extra edges existing in the data graph connecting to the matched nodes may lead to two matched nodes of two graphs having different adjacency structures and often being identified as distinct objects. Due to the extra edges, the existing learning based methods often fail to generate sufficiently similar node-level embeddings for matched nodes. This study proposes a novel Adaptive Edge-Deleting Network (AEDNet) for subgraph matching. The proposed method is trained in an end-to-end fashion. In AEDNet, a novel sample-wise adaptive edge-deleting mechanism removes extra edges to ensure consistency of adjacency structure of matched nodes, while a unidirectional cross-propagation mechanism ensures consistency of features of matched nodes. We applied the proposed method on six datasets with graph sizes varying from 20 to 2300. Our evaluations on six open datasets demonstrate that the proposed AEDNet outperforms six state-of-the-arts and is much faster than the exact methods on large graphs.


title: Structured Mixture of Continuation-ratio Logits Models for Ordinal Regression

abstract: We develop a nonparametric Bayesian modeling approach to ordinal regression based on priors placed directly on the discrete distribution of the ordinal responses. The prior probability models are built from a structured mixture of multinomial distributions. We leverage a continuation-ratio logits representation to formulate the mixture kernel, with mixture weights defined through the logit stick-breaking process that incorporates the covariates through a linear function. The implied regression functions for the response probabilities can be expressed as weighted sums of parametric regression functions, with covariate-dependent weights. Thus, the modeling approach achieves flexible ordinal regression relationships, avoiding linearity or additivity assumptions in the covariate effects. Model flexibility is formally explored through the Kullback-Leibler support of the prior probability model. A key model feature is that the parameters for both the mixture kernel and the mixture weights can be associated with a continuation-ratio logits regression structure. Hence, an efficient and relatively easy to implement posterior simulation method can be designed, using P\'olya-Gamma data augmentation. Moreover, the model is built from a conditional independence structure for category-specific parameters, which results in additional computational efficiency gains through partial parallel sampling. In addition to the general mixture structure, we study simplified model versions that incorporate covariate dependence only in the mixture kernel parameters or only in the mixture weights. For all proposed models, we discuss approaches to prior specification and develop Markov chain Monte Carlo methods for posterior simulation. The methodology is illustrated with several synthetic and real data examples.


title: Quantum metallic state at the titanium sesquioxide heterointerface

abstract: The emergence of the quantum metallic state marked by a saturating finite electrical resistance in the zero-temperature limit in a variety of two-dimensional superconductors injects an exciting momentum to the realm of unconventional superconductivity. Despite much research efforts over last few decades, there is not yet a general consensus on the nature of this unexpected quantum metal. Here, we report the observation of a unique quantum metallic state within the hallmark of Bose-metal characterized by a saturating resistance and a simultaneously vanishing Hall resistance in the titanium sesquioxide heterointerface superconductor Ti$_2$O$_3$/GaN. Remarkably, the quantum bosonic metallic state tuned by a magnetic field in the vicinity of the two-dimensional superconductivity-metal transition persists in the normal phase, suggesting the existence of composite bosons, formed by electron Cooper pairs, even in the normal phase. Our findings indicate the existence of preformed electron Cooper pairs in the normal phase of heterointerface superconductors, and shed new light on the pairing mechanism underlying unconventional superconductivity.


title: Performance Analysis of LEO Satellite-Based IoT Networks in the Presence of Interference

abstract: This paper presents a star-of-star topology for internet-of-things (IoT) networks using mega low-Earth-orbit constellations. The proposed topology enables IoT users to broadcast their sensed data to multiple satellites simultaneously over a shared channel, which is then relayed to the ground station (GS) using amplify-and-forward relaying. The GS coherently combines the signals from multiple satellites using maximal ratio combining. To analyze the performance of the proposed topology in the presence of interference, a comprehensive outage probability (OP) analysis is performed, assuming imperfect channel state information at the GS. The paper employs stochastic geometry to model the random locations of satellites, making the analysis general and independent of any specific constellation. Furthermore, the paper examines successive interference cancellation (SIC) and capture model (CM)-based decoding schemes at the GS to mitigate interference. The average OP for the CM-based scheme and the OP of the best user for the SIC scheme are derived analytically. The paper also presents simplified expressions for the OP under a high signal-to-noise ratio (SNR) assumption, which are utilized to optimize the system parameters for achieving a target OP. The simulation results are consistent with the analytical expressions and provide insights into the impact of various system parameters, such as mask angle, altitude, number of satellites, and decoding order. The findings of this study demonstrate that the proposed topology can effectively leverage the benefits of multiple satellites to achieve the desired OP and enable burst transmissions without coordination among IoT users, making it an attractive choice for satellite-based IoT networks.


title: A $C^0$ Linear Finite Element Method for a Second Order Elliptic Equation in Non-Divergence Form with Cordes Coefficients

abstract: In this paper, we develop a gradient recovery based linear (GRBL) finite element method (FEM) and a Hessian recovery based linear (HRBL) FEM for second order elliptic equations in non-divergence form. The elliptic equation is casted into a symmetric non-divergence weak formulation, in which second order derivatives of the unknown function are involved. We use gradient and Hessian recovery operators to calculate the second order derivatives of linear finite element approximations. Although, thanks to low degrees of freedom (DOF) of linear elements, the implementation of the proposed schemes is easy and straightforward, the performances of the methods are competitive. The unique solvability and the $H^2$ seminorm error estimate of the GRBL scheme are rigorously proved. Optimal error estimates in both the $L^2$ norm and the $H^1$ seminorm have been proved when the coefficient is diagonal, which have been confirmed by numerical experiments. Superconvergence in errors has also been observed. Moreover, our methods can handle computational domains with curved boundaries without loss of accuracy from approximation of boundaries. Finally, the proposed numerical methods have been successfully applied to solve fully nonlinear Monge-Amp\`{e}re equations.