title: Shafarevich-Tate groups of holomorphic Lagrangian fibrations

abstract: Consider a Lagrangian fibration $\pi\colon X\to \mathbb P^n$ on a hyperk\"ahler manifold $X$. There are two ways to construct a holomorphic family of deformations of $\pi$ over $\mathbb C$. The first one is known under the name Shafarevich-Tate family while the second one is the degenerate twistor family constructed by Verbitsky. We show that both families coincide. We prove that for a very general $X$ all members of the Shafarevich-Tate family are K\"ahler. There is a related notion of the Shafarevich-Tate group associated to a Lagrangian fibration. Its connected component of unity can be shown to be isomorphic to $\mathbb C/\Lambda$ where $\Lambda$ is a finitely generated subgroup of $\mathbb C$ and $\mathbb C$ is thought of as the base of the Shafarevich-Tate family. We show that for a very general $X$, projective deformations in the Shafarevich-Tate family correspond to the torsion points in the connected component of unity of the Shafarevich-Tate group. A sufficient condition for a Lagrangian fibration $X$ to be projective is existence of a holomorphic section. We find sufficient cohomological conditions for existence of a deformation in the Shafarevich-Tate family that admits a section.


title: Running primordial perturbations: Inflationary Dynamics and Observational Constraints

abstract: Inflationary cosmology proposes that the early Universe undergoes accelerated expansion, driven, in simple scenarios, by a single scalar field, or inflaton. The form of the inflaton potential determines the initial spectra of density perturbations and gravitational waves. We show that constraints on the duration of inflation together with the BICEP3/Keck bounds on the gravitational wave background imply that higher derivatives of the potential are nontrivial with a confidence of 99%. Such terms contribute to the scale-dependence, or running, of the density perturbation spectrum. We clarify the ``universality classes'' of inflation in this limit showing that a very small gravitational wave background can be correlated with a larger running. If pending experiments do not observe a gravitational wave background the running will be at the threshold of detectability if inflation is well-described at third-order in the slow roll expansion.


title: One-dimensional discrete Hardy and Rellich inequalities on integers

abstract: In this paper, we consider a weighted version of one-dimensional discrete Hardy inequalities with power weights of the form $n^\alpha$. We prove the inequality when $\alpha$ is an even natural number with the sharp constant and remainder terms. We also find explicit constants in standard and weighted Rellich inequalities and its higher order versions. As a by-product of this work we derive a combinatorial identity using purely analytic methods. This suggests a correlation between combinatorial identities and functional identities.


title: A family tree of two-dimensional magnetic materials with tunable topological properties

abstract: Two-dimensional (2D) magnetic materials empowered with nontrivial band topology may lead to the emergence of exotic quantum states with significant application potentials. Here we predict a family tree of 2D magnetic materials with tunable topological properties, starting from the parental materials of CrI3 and CrBr3. The underlying design principle is that, by substituting the alternating sites of the Cr honeycomb lattice sandwiched between the halogen layers with V or Mn, the parental materials of trivial ferromagnetic insulators are ripe to be converted into topological systems. Specifically, our first-principles calculations show that, due to the elegant interplay between bandgap narrowing and spin-orbital coupling, CrI3 branches into high-temperature quantum anomalous Hall insulators of CrVI6 and CrMnI6 with different topological invariants, while CrBr3 branches into topological half-metals of CrVBr6 and CrMnBr6. Those novel 2D magnets are also shown to be easily exfoliated from their bulk counterparts. The present study is thus geared to advance the field of 2D magnetic materials into the topologically nontrivial realm.


title: DB-BERT: a Database Tuning Tool that "Reads the Manual"

abstract: DB-BERT is a database tuning tool that exploits information gained via natural language analysis of manuals and other relevant text documents. It uses text to identify database system parameters to tune as well as recommended parameter values. DB-BERT applies large, pre-trained language models (specifically, the BERT model) for text analysis. During an initial training phase, it fine-tunes model weights in order to translate natural language hints into recommended settings. At run time, DB-BERT learns to aggregate, adapt, and prioritize hints to achieve optimal performance for a specific database system and benchmark. Both phases are iterative and use reinforcement learning to guide the selection of tuning settings to evaluate (penalizing settings that the database system rejects while rewarding settings that improve performance). In our experiments, we leverage hundreds of text documents about database tuning as input for DB-BERT. We compare DB-BERT against various baselines, considering different benchmarks (TPC-C and TPC-H), metrics (throughput and run time), as well as database systems (Postgres and MySQL). In all cases, DB-BERT finds the best parameter settings among all compared methods. The code of DB-BERT is available online at https://itrummer.github.io/dbbert/.


title: 100 GHz Micrometer compact broadband Monolithic ITO Mach Zehnder Interferometer Modulator enabling 3500 times higher Packing Density

abstract: Electro-optic modulators provide a key function in optical transceivers and increasingly in photonic programmable Application Specific Integrated Circuits (ASICs) for machine learning and signal processing. However, both foundry ready silicon based modulators and conventional material based devices utilizing Lithium niobate fall short in simultaneously providing high chip packaging density and fast speed. Current driven ITO based modulators have the potential to achieve both enabled by efficient light matter interactions. Here, we introduce micrometer compact Mach Zehnder Interferometer (MZI) based modulators capable of exceeding 100 GHz switching rates. Integrating ITO thin films atop a photonic waveguide, spectrally broadband, and compact MZI phase shifter. Remarkably, this allows integrating more than 3500 of these modulators within the same chip area as only one single silicon MZI modulator. The modulator design introduced here features a holistic photonic, electronic, and RF-based optimization and includes an asymmetric MZI tuning step to optimize the Extinction Ratio (ER) to Insertion Loss (IL) and dielectric thickness sweep to balance the tradeoffs between ER and speed. Driven by CMOS compatible bias voltage levels, this device is the first to address next generation modulator demands for processors of the machine intelligence revolution, in addition to the edge and cloud computing demands as well as optical transceivers alike.


title: CNN-tuned spatial filters for P- and S-wave decomposition and applications in elastic imaging

abstract: P- and S-wave decomposition is essential for imaging multi-component seismic data in elastic media. A data-driven workflow is proposed to obtain a set of spatial filters that are highly accurate and artifact-free in decomposing the P- and S-waves in 2D isotropic elastic wavefields. The filters are formulated initially by inverse Fourier transforms of the wavenumber-domain operators, and then are tuned in a convolutional neural network to improve accuracy using synthetic snapshots. Spatial filters are flexible for decomposing P-and S-waves at any time step without performing Fourier transforms over the entire wavefield snapshots, and thus are suitable for target-oriented imaging. Snapshots from synthetic data show that the network-tuned spatial filters can decompose P- and S-waves with improved decomposition accuracy compared with other space-domain PS decomposition methods. Elastic reverse-time migration using P- and S-waves decomposed from the proposed algorithm shows reduced artifacts in the presence of a high velocity contrast.


title: Bialgebras, Frobenius algebras and associative Yang-Baxter equations for Rota-Baxter algebras

abstract: Rota-Baxter operators and bialgebras go hand in hand in their applications, such as in the Connes-Kreimer approach to renormalization and the operator approach to the classical Yang-Baxter equation. We establish a bialgebra structure that is compatible with the Rota-Baxter operator, called the Rota-Baxter antisymmetric infinitesimal (ASI) bialgebra. This bialgebra is characterized by generalizations of matched pairs of algebras and double constructions of Frobenius algebras to the context of Rota-Baxter algebras. The study of the coboundary case leads to an enrichment of the associative Yang-Baxter equation (AYBE) to Rota-Baxter algebras. Antisymmetric solutions of the equation are used to construct Rota-Baxter ASI bialgebras. The notions of an $\mathcal{O}$-operator on a Rota-Baxter algebra and a Rota-Baxter dendriform algebra are also introduced to produce solutions of the AYBE in Rota-Baxter algebras and thus to provide Rota-Baxter ASI bialgebras. An unexpected byproduct is that a Rota-Baxter ASI bialgebra of weight zero gives rise to a quadri-bialgebra instead of bialgebra constructions for the dendriform algebra.


title: Quantum probability from temporal structure

abstract: The Born probability measure describes the statistics of measurements in which observers self-locate themselves in some region of reality. In $\psi$-ontic quantum theories, reality is directly represented by the wavefunction. We show that quantum probabilities may be identified with fractions of a universal multiple-time wavefunction containing both causal and retrocausal temporal parts. This wavefunction is defined in an appropriately generalized history space on the Keldysh time contour. Our deterministic formulation of quantum mechanics replaces the initial condition of standard Schr\"odinger dynamics with a network of `fixed points' defining quantum histories on the contour. The Born measure is derived by summing up the wavefunction along these histories. We then apply the same technique to the derivation of the statistics of measurements with pre- and post-selection.

categories: quant-ph cond-mat.stat-mech math-ph math.MP physics.hist-ph


title: Compact Multi-level Sparse Neural Networks with Input Independent Dynamic Rerouting

abstract: Deep neural networks (DNNs) have shown to provide superb performance in many real life applications, but their large computation cost and storage requirement have prevented them from being deployed to many edge and internet-of-things (IoT) devices. Sparse deep neural networks, whose majority weight parameters are zeros, can substantially reduce the computation complexity and memory consumption of the models. In real-use scenarios, devices may suffer from large fluctuations of the available computation and memory resources under different environment, and the quality of service (QoS) is difficult to maintain due to the long tail inferences with large latency. Facing the real-life challenges, we propose to train a sparse model that supports multiple sparse levels. That is, a hierarchical structure of weights are satisfied such that the locations and the values of the non-zero parameters of the more-sparse sub-model area subset of the less-sparse sub-model. In this way, one can dynamically select the appropriate sparsity level during inference, while the storage cost is capped by the least sparse sub-model. We have verified our methodologies on a variety of DNN models and tasks, including the ResNet-50, PointNet++, GNMT, and graph attention networks. We obtain sparse sub-models with an average of 13.38% weights and 14.97% FLOPs, while the accuracies are as good as their dense counterparts. More-sparse sub-models with 5.38% weights and 4.47% of FLOPs, which are subsets of the less-sparse ones, can be obtained with only 3.25% relative accuracy loss.


title: Hiding Signal Strength Interference from Outside Adversaries

abstract: The presence of people can be detected by passively observing the signal strength of Wifi and related forms of communication. This paper tackles the question of how and when can this be prevented by adjustments to the transmitted signal strength, and other similar measures. The main contribution of this paper is a formal framework to analyze this problem, and the identification of several scenarios and corresponding protocols which can prevent or limit the inference from passive signal strength snooping.


title: $p$-energies on p.c.f. self-similar sets

abstract: We study $p$-energies on post critically finite (p.c.f.) self-similar sets for $1<p<\infty$, as limits of discrete $p$-energies on approximation graphs, extending the construction of Dirichlet forms, the $p=2$ setting. By suitably enlarging the choices of discrete $p$-energies, and employing the energy averaging method developed by Kusuoka-Zhou, we prove the existence of symmetric $p$-energies on affine nested fractals, and extend Sabot's celebrated criterion for existence and non-existence of Dirichlet forms on p.c.f. self-similar sets to the $1<p<\infty$ setting.


title: On the Size and Width of the Decoder of a Boolean Threshold Autoencoder

abstract: In this paper, we study the size and width of autoencoders consisting of Boolean threshold functions, where an autoencoder is a layered neural network whose structure can be viewed as consisting of an encoder, which compresses an input vector to a lower dimensional vector, and a decoder which transforms the low-dimensional vector back to the original input vector exactly (or approximately). We focus on the decoder part, and show that $\Omega(\sqrt{Dn/d})$ and $O(\sqrt{Dn})$ nodes are required to transform $n$ vectors in $d$-dimensional binary space to $D$-dimensional binary space. We also show that the width can be reduced if we allow small errors, where the error is defined as the average of the Hamming distance between each vector input to the encoder part and the resulting vector output by the decoder.


title: Understanding Software Architecture Erosion: A Systematic Mapping Study

abstract: Architecture erosion (AEr) can adversely affect software development and has received significant attention in the last decade. However, there is an absence of a comprehensive understanding of the state of research about the reasons and consequences of AEr, and the countermeasures to address AEr. This work aims at systematically investigating, identifying, and analyzing the reasons, consequences, and ways of detecting and handling AEr. With 73 studies included, the main results are as follows: (1) AEr manifests not only through architectural violations and structural issues but also causing problems in software quality and during software evolution; (2) non-technical reasons that cause AEr should receive the same attention as technical reasons, and practitioners should raise awareness of the grave consequences of AEr, thereby taking actions to tackle AEr-related issues; (3) a spectrum of approaches, tools, and measures has been proposed and employed to detect and tackle AEr; and (4) three categories of difficulties and five categories of lessons learned on tackling AEr were identified. The results can provide researchers a comprehensive understanding of AEr and help practitioners handle AEr and improve the sustainability of their architecture. More empirical studies are required to investigate the practices of detecting and addressing AEr in industrial settings.


title: Nearly Optimal Policy Optimization with Stable at Any Time Guarantee

abstract: Policy optimization methods are one of the most widely used classes of Reinforcement Learning (RL) algorithms. However, theoretical understanding of these methods remains insufficient. Even in the episodic (time-inhomogeneous) tabular setting, the state-of-the-art theoretical result of policy-based method in \citet{shani2020optimistic} is only $\tilde{O}(\sqrt{S^2AH^4K})$ where $S$ is the number of states, $A$ is the number of actions, $H$ is the horizon, and $K$ is the number of episodes, and there is a $\sqrt{SH}$ gap compared with the information theoretic lower bound $\tilde{\Omega}(\sqrt{SAH^3K})$. To bridge such a gap, we propose a novel algorithm Reference-based Policy Optimization with Stable at Any Time guarantee (\algnameacro), which features the property "Stable at Any Time". We prove that our algorithm achieves $\tilde{O}(\sqrt{SAH^3K} + \sqrt{AH^4K})$ regret. When $S > H$, our algorithm is minimax optimal when ignoring logarithmic factors. To our best knowledge, RPO-SAT is the first computationally efficient, nearly minimax optimal policy-based algorithm for tabular RL.


title: Watch Those Words: Video Falsification Detection Using Word-Conditioned Facial Motion

abstract: In today's era of digital misinformation, we are increasingly faced with new threats posed by video falsification techniques. Such falsifications range from cheapfakes (e.g., lookalikes or audio dubbing) to deepfakes (e.g., sophisticated AI media synthesis methods), which are becoming perceptually indistinguishable from real videos. To tackle this challenge, we propose a multi-modal semantic forensic approach to discover clues that go beyond detecting discrepancies in visual quality, thereby handling both simpler cheapfakes and visually persuasive deepfakes. In this work, our goal is to verify that the purported person seen in the video is indeed themselves by detecting anomalous facial movements corresponding to the spoken words. We leverage the idea of attribution to learn person-specific biometric patterns that distinguish a given speaker from others. We use interpretable Action Units (AUs) to capture a person's face and head movement as opposed to deep CNN features, and we are the first to use word-conditioned facial motion analysis. We further demonstrate our method's effectiveness on a range of fakes not seen in training including those without video manipulation, that were not addressed in prior work.


title: CrSbSe3: a pseudo one-dimensional ferromagnetic semiconductor

abstract: Low-dimensional magnetic materials have attracted much attention due to their novel properties and high potential for spintronic applications. In this work, we study the electronic structure and magnetic properties of the pseudo one-dimensional compound CrSbSe$_3$, using density functional calculations, superexchange model analyses, and Monte Carlo simulations. We find that CrSbSe$_3$ is a ferromagnetic (FM) semiconductor with a band gap of about 0.65 eV. The FM couplings within each zig-zag spin chain are due to the Cr-Se-Cr superexchange with the near-90$^\circ$ bonds, and the inter-chain FM couplings are one order of magnitude weaker. By inclusion of the spin-orbit coupling (SOC) effects, our calculations reproduce the experimental observation of the easy magnetization $a$ axis and the hard $b$ axis (the spin-chain direction), with the calculated moderate magnetic anisotropy of 0.19 meV/Cr. Moreover, we identify the nearly equal contributions from the single ion anisotropy of Cr, and from the exchange anisotropy due to the strong SOC of the heavy elements Sb and Se and their couplings with Cr. Using the parameters of the magnetic coupling and anisotropy from the above calculations, our Monte Carlo simulations yield the Curie temperature $T_{\rm{C}}$ of 108 K.


title: CADV: A software visualization approach for code annotations distribution

abstract: Code annotations is a widely used feature in Java systems to configure custom metadata on programming elements. Their increasing presence creates the need for approaches to assess and comprehend their usage and distribution. In this context, software visualization has been studied and researched to improve program comprehension in different aspects. This study aimed at designing a software visualization approach that graphically displays how code annotations are distributed and organized in a software system and developing a tool, as a reference implementation of the approach, to generate views and interact with users. We conducted an empirical evaluation through questionnaires and interviews to evaluate our visualization approach considering four aspects: effectiveness for program comprehension, perceived usefulness, perceived ease of use, and suitability for the intended audience. The resulting data was used to perform a qualitative and quantitative analysis. The tool identifies package responsibilities providing visual information about their annotations at different levels. Using the developed tool, the participants achieved a high correctness rate in the program comprehension tasks and performed very well in questions about the overview of the system under analysis. Finally, participants perceived that the tool outperforms existing approaches for code inspection when searching for information related to code annotations. The results show that the visualization approach using the developed tool is effective in program comprehension tasks related to code annotations, which can also be used to identify responsibilities in the application packages. Moreover, it was evaluated as suitable for newcomers to overview the usage of annotations in the system and for architects to perform a deep analysis that can potentially detect misplaced annotations and abnormal growths on their usage.


title: $CP$-violating observables in four-body $B\rightarrow \phi(\rightarrow K\bar K)K^*(\rightarrow K\pi)$ decays

abstract: We analyse the four-body $B\rightarrow \phi(\rightarrow K\bar K)K^*(\rightarrow K\pi)$ decays in the perturbative QCD approach,where the invariant mass of $K\bar K$($K\pi$) system is limited in a window of $\pm 15$ MeV ($\pm150$ MeV) around the $\phi(K^*(892))$ mass.In addition to the P wave resonances,two important S wave backgrounds in the selected invariant mass region are also accounted for. Angular momentum conservation allows six helicity amplitudes to contribute,including three P waves, two single S waves,and one double S wave. We calculated the branching ratio for each component and found sizable S wave contributions,coincide with the experimental observation.The obtained branching ratios of $B^{0(+)}\rightarrow \phi K^{*0(+)}$ are comparable with the previous predictions and support the measurements, whereas the predicted $\mathcal{B}(B^0_s\rightarrow \phi \bar K^{*0})$ is smaller than the world average. The longitudinal polarizations are predicted to be around 0.7,consistent with previous PQCD results but larger than the data. Aside from the direct CP asymmetries,the true and fake triple product asymmetries(TPAs) are calculated in this work. In the case of neutral modes, both direct CP asymmetries and true TPAs are expected to be zero due to the vanishing weak phase difference. The direct CP asymmetries for the $B^+$ mode are predicted to be tiny,since the tree contributions are suppressed with respect to the penguin ones. The true asymmetries have shown no significant deviations from zero.In contrast,large fake asymmetries are observed in these decays,indicating the presence of significant final state interactions.We give the predictions of the S wave induced TPAs for the first time,which is consistent with LHCb data and would be checked with future measurements from Belle and BABAR experiments if the S wave components can be properly taken into account in angular analysis.


title: Origin of cool cores, cold fronts and spiral structures in cool core clusters of galaxies

abstract: We consider a situation in which a brightest cluster galaxy (BCG) moves in ambient hot gas in the central region of a cool core cluster of galaxies, following the study by Inoue (2014, PASJ, 66, 60). In the rest frame of the BCG, the hot gas is supposed to flow toward the BCG in parallel from a sufficiently large distance. Then, it is expected that only the gas flowing with the impact parameter less than a critical value is trapped by the gravitation field of the BCG because of the efficient radiative cooling, getting a cooling flow, and that the remaining outer gas can get over the potential well. In such a circumstance, we can draw the following picture: A boundary layer between the out-flowing gas and the trapped gas arises around the stagnation point at the back side of the BCG. Since the boundary temperature is so low as to be X-ray dim, the boundary could be observed as the cold front in X-rays. The trapped gas once stagnates on the inner side of the boundary and starts in-falling toward the BCG. Since the wandering motion of the BCG is likely to have a rotational component, the Coriolis force induces a rotational motion in the in-falling flow from the stagnation place to the BCG, forming a spiral structure around the BCG. The spiraling flow turns the BCG on the up stream side of the main flow from the far outside, and arises another boundary layer having contact discontinuity with the main hot gas flow. These pictures well reproduce the observed features such as cool cores, cold fronts, and spiral structures. It can also be resolved how the cooling flow is suppressed from what the cooling flow hypothesis predicts, without any heating mechanism.


title: Structured Semantic Transfer for Multi-Label Recognition with Partial Labels

abstract: Multi-label image recognition is a fundamental yet practical task because real-world images inherently possess multiple semantic labels. However, it is difficult to collect large-scale multi-label annotations due to the complexity of both the input images and output label spaces. To reduce the annotation cost, we propose a structured semantic transfer (SST) framework that enables training multi-label recognition models with partial labels, i.e., merely some labels are known while other labels are missing (also called unknown labels) per image. The framework consists of two complementary transfer modules that explore within-image and cross-image semantic correlations to transfer knowledge of known labels to generate pseudo labels for unknown labels. Specifically, an intra-image semantic transfer module learns image-specific label co-occurrence matrix and maps the known labels to complement unknown labels based on this matrix. Meanwhile, a cross-image transfer module learns category-specific feature similarities and helps complement unknown labels with high similarities. Finally, both known and generated labels are used to train the multi-label recognition models. Extensive experiments on the Microsoft COCO, Visual Genome and Pascal VOC datasets show that the proposed SST framework obtains superior performance over current state-of-the-art algorithms. Codes are available at https://github.com/HCPLab-SYSU/HCP-MLR-PL.


title: Towards a discovery of BSM physics from the Cabibbo angle anomaly

abstract: New developments in both the theories and experiments related to the extraction of the top-row Cabibbo-Kobayashi-Maskawa matrix elements $V_{ud}$ and $V_{us}$ led to a series of new anomalies, for instance the apparent violation of the top-row unitarity relation. It is important to further reduce all the associated Standard Model theory uncertainties in order to better understand whether such observations point towards the possibility of physics beyond the Standard Model, or rather some unexpectedly large Standard Model effects. This requires improved studies of tree-level and higher-order Standard Model corrections that enter the beta decays of pions, neutron, nuclei and kaons. We will briefly review the recent progress along this direction and discuss possible improvements in the future.


title: Shortcuts to Quantum Approximate Optimization Algorithm

abstract: The Quantum Approximate Optimization Algorithm (QAOA) is a quantum-classical hybrid algorithm intending to find the ground state of a target Hamiltonian. Theoretically, QAOA can obtain the approximate solution if the quantum circuit is deep enough. Actually, the performance of QAOA decreases practically if the quantum circuit is deep since near-term devices are not noise-free and the errors caused by noise accumulate as the quantum circuit increases. In order to reduce the depth of quantum circuits, we propose a new ansatz dubbed as "Shortcuts to QAOA" (S-QAOA), S-QAOA provides shortcuts to the ground state of target Hamiltonian by including more two-body interactions and releasing the parameter freedoms. To be specific, besides the existing ZZ interaction in the QAOA ansatz, other two-body interactions are introduced in the S-QAOA ansatz such that the approximate solutions could be obtained with smaller circuit depth. Considering the MaxCut problem and Sherrington-Kirkpatrick (SK) model, numerically computation shows the YY interaction has the best performance. The reason for this might arise from the counterdiabatic effect generated by YY interaction. On top of this, we release the freedom of parameters of two-body interactions, which a priori do not necessarily have to be fully identical, and numerical results show that it is worth paying the extra cost of having more parameter freedom since one has a greater improvement on success rate.


title: Reinforcement Learning based Sequential Batch-sampling for Bayesian Optimal Experimental Design

abstract: Engineering problems that are modeled using sophisticated mathematical methods or are characterized by expensive-to-conduct tests or experiments, are encumbered with limited budget or finite computational resources. Moreover, practical scenarios in the industry, impose restrictions, based on logistics and preference, on the manner in which the experiments can be conducted. For example, material supply may enable only a handful of experiments in a single-shot or in the case of computational models one may face significant wait-time based on shared computational resources. In such scenarios, one usually resorts to performing experiments in a manner that allows for maximizing one's state-of-knowledge while satisfying the above mentioned practical constraints. Sequential design of experiments (SDOE) is a popular suite of methods, that has yielded promising results in recent years across different engineering and practical problems. A common strategy, that leverages Bayesian formalism is the Bayesian SDOE, which usually works best in the one-step-ahead or myopic scenario of selecting a single experiment at each step of a sequence of experiments. In this work, we aim to extend the SDOE strategy, to query the experiment or computer code at a batch of inputs. To this end, we leverage deep reinforcement learning (RL) based policy gradient methods, to propose batches of queries that are selected taking into account entire budget in hand. The algorithm retains the sequential nature, inherent in the SDOE, while incorporating elements of reward based on task from the domain of deep RL. A unique capability of the proposed methodology is its ability to be applied to multiple tasks, for example optimization of a function, once its trained. We demonstrate the performance of the proposed algorithm on a synthetic problem, and a challenging high-dimensional engineering problem.


title: Pixel-Stega: Generative Image Steganography Based on Autoregressive Models

abstract: In this letter, we explored generative image steganography based on autoregressive models. We proposed Pixel-Stega, which implements pixel-level information hiding with autoregressive models and arithmetic coding algorithm. Firstly, one of the autoregressive models, PixelCNN++, is utilized to produce explicit conditional probability distribution of each pixel. Secondly, secret messages are encoded to the selection of pixels through steganographic sampling (stegosampling) based on arithmetic coding. We carried out qualitative and quantitative assessment on gray-scale and colour image datasets. Experimental results show that Pixel-Stega is able to embed secret messages adaptively according to the entropy of the pixels to achieve both high embedding capacity (up to 4.3 bpp) and nearly perfect imperceptibility (about 50% detection accuracy).


title: Cram\'er-type Moderate deviations under local dependence

abstract: We establish Cram\'er-type moderate deviation theorems for sums of locally dependent random variables and combinatorial central limit theorems. Under some mild exponential moment conditions, optimal error bounds and convergence ranges are obtained. Our main results are more general or shaper than the existing results in the literature. The main results follows from a more general Cram\'er-type moderate deviation theorem for dependent random variables without any boundedness assumptions, which is of independent interest. The proofs couple Stein's method with a recursive argument.


title: The entropic barrier is $n$-self-concordant

abstract: For any convex body $K \subseteq \mathbb R^n$, S. Bubeck and R. Eldan introduced the entropic barrier on $K$ and showed that it is a $(1+o(1)) \, n$-self-concordant barrier. In this note, we observe that the optimal bound of $n$ on the self-concordance parameter holds as a consequence of the dimensional Brascamp-Lieb inequality.


title: Task-Oriented Image Transmission for Scene Classification in Unmanned Aerial Systems

abstract: The vigorous developments of Internet of Things make it possible to extend its computing and storage capabilities to computing tasks in the aerial system with collaboration of cloud and edge, especially for artificial intelligence (AI) tasks based on deep learning (DL). Collecting a large amount of image/video data, Unmanned aerial vehicles (UAVs) can only handover intelligent analysis tasks to the back-end mobile edge computing (MEC) server due to their limited storage and computing capabilities. How to efficiently transmit the most correlated information for the AI model is a challenging topic. Inspired by the task-oriented communication in recent years, we propose a new aerial image transmission paradigm for the scene classification task. A lightweight model is developed on the front-end UAV for semantic blocks transmission with perception of images and channel conditions. In order to achieve the tradeoff between transmission latency and classification accuracy, deep reinforcement learning (DRL) is used to explore the semantic blocks which have the best contribution to the back-end classifier under various channel conditions. Experimental results show that the proposed method can significantly improve classification accuracy compared to the fixed transmission strategy and traditional content perception methods.


title: District Cooling System Control for Providing Operating Reserve based on Safe Deep Reinforcement Learning

abstract: Heating, ventilation, and air conditioning (HVAC) systems are well proved to be capable to provide operating reserve for power systems. As a type of large-capacity and energy-efficient HVAC system (up to 100 MW), district cooling system (DCS) is emerging in modern cities and has huge potential to be regulated as a flexible load. However, strategically controlling a DCS to provide flexibility is challenging, because one DCS services multiple buildings with complex thermal dynamics and uncertain cooling demands. Improper control may lead to significant thermal discomfort and even deteriorate the power system's operation security. To address the above issues, we propose a model-free control strategy based on the deep reinforcement learning (DRL) without the requirement of accurate system model and uncertainty distribution. To avoid damaging "trial & error" actions that may violate the system's operation security during the training process, we further propose a safe layer combined to the DRL to guarantee the satisfaction of critical constraints, forming a safe-DRL scheme. Moreover, after providing operating reserve, DCS increases power and tries to recover all the buildings' temperature back to set values, which may probably cause an instantaneous peak-power rebound and bring a secondary impact on power systems. Therefore, we design a self-adaption reward function within the proposed safe-DRL scheme to constrain the peak-power effectively. Numerical studies based on a realistic DCS demonstrate the effectiveness of the proposed methods.


title: Augmented Contrastive Self-Supervised Learning for Audio Invariant Representations

abstract: Improving generalization is a major challenge in audio classification due to labeled data scarcity. Self-supervised learning (SSL) methods tackle this by leveraging unlabeled data to learn useful features for downstream classification tasks. In this work, we propose an augmented contrastive SSL framework to learn invariant representations from unlabeled data. Our method applies various perturbations to the unlabeled input data and utilizes contrastive learning to learn representations robust to such perturbations. Experimental results on the Audioset and DESED datasets show that our framework significantly outperforms state-of-the-art SSL and supervised learning methods on sound/event classification tasks.


title: General-purpose cooperativeness and altruism in humans: elements of the mathematical framework for the Interdependence Hypothesis

abstract: We propose a decision-making model for joint intentionality by interpreting it as group-mindedness at the microlevel. We apply this model to give a formal justification of the first part of the Interdependence Hypothesis due to Tomasello et al. [Current Anthropology, 2012] which asserts that the emergence of joint intentionality evolved due to the challenges of difficult collaborative foraging practices among early humans, and that its evolution led to robust collaboration and some form of altruism. In another application of the microlevel group-mindedness we consider the problem of establishing cooperation in high-risk-of-defection strategic conflicts and we show that the emergence of cooperation in such situations can be explained in the context of cultural group selection as the result of adaptive learning.


title: Mitigating Barren Plateaus with Transfer-learning-inspired Parameter Initializations

abstract: Variational quantum algorithms (VQAs) are widely applied in the noisy intermediate-scale quantum era and are expected to demonstrate quantum advantage. However, training VQAs faces difficulties, one of which is the so-called barren plateaus (BP) phenomenon, where gradients of cost functions vanish exponentially with the number of qubits. In this paper, inspired by transfer learning, where knowledge of pre-solved tasks could be further used in a different but related work with training efficiency improved, we report a parameter initialization method to mitigate BP. In the method, a small-sized task is solved with a VQA. Then the ansatz and its optimum parameters are transferred to tasks with larger sizes. Numerical simulations show that this method could mitigate BP and improve training efficiency. A brief discussion on how this method can work well is also provided. This work provides a reference for mitigating BP, and therefore, VQAs could be applied to more practical problems.


title: An adaptation of InfoMap to absorbing random walks using absorption-scaled graphs

abstract: InfoMap is a popular approach to detect densely connected "communities" of nodes in networks. To detect such communities, InfoMap uses random walks and ideas from information theory. Motivated by the dynamics of disease spread on networks, whose nodes can have heterogeneous disease-removal rates, we adapt InfoMap to absorbing random walks. To do this, we use absorption-scaled graphs (in which edge weights are scaled according to absorption rates) and Markov time sweeping. One of our adaptations of InfoMap converges to the standard version of InfoMap in the limit in which the node-absorption rates approach $0$. We demonstrate that the community structure that one obtains using our adaptations of InfoMap can differ markedly from the community structure that one detects using methods that do not account for node-absorption rates. We also illustrate that the community structure that is induced by heterogeneous absorption rates can have important implications for susceptible-infected-recovered (SIR) dynamics on ring-lattice networks. For example, in some situations, the outbreak duration is maximized when a moderate number of nodes have large node-absorption rates.

categories: cs.SI cs.LG math.PR nlin.AO physics.soc-ph


title: A Kapitza Pendulum for Ultracold Atoms

abstract: We report on the experimental realization of a Kapitza pendulum for ultracold atoms. Using time-periodic attractive and repulsive Gaussian potentials, we create an effective trap for ultracold neutral atoms in a regime where the time average of the potential is equal to zero. We analyze the role of experimental imperfections, the stability of the trapped atomic cloud, and the magnitude of the effective potential. We find good agreement with the high-frequency expansion of the underlying system dynamics. Our experimental approach opens up new possibilities to study Floquet systems of neutral atoms.


title: Joint Learning of Linear Time-Invariant Dynamical Systems

abstract: Linear time-invariant systems are very popular models in system theory and applications. A fundamental problem in system identification that remains rather unaddressed in extant literature is to leverage commonalities amongst related linear systems to estimate their transition matrices more accurately. To address this problem, the current paper investigates methods for jointly estimating the transition matrices of multiple systems. It is assumed that the transition matrices are unknown linear functions of some unknown shared basis matrices. We establish finite-time estimation error rates that fully reflect the roles of trajectory lengths, dimension, and number of systems under consideration. The presented results are fairly general and show the significant gains that can be achieved by pooling data across systems in comparison to learning each system individually. Further, they are shown to be robust against model misspecifications. To obtain the results, we develop novel techniques that are of interest for addressing similar joint-learning problems. They include tightly bounding estimation errors in terms of the eigen-structures of transition matrices, establishing sharp high probability bounds for singular values of dependent random matrices, and capturing effects of misspecified transition matrices as the systems evolve over time.


title: The Fourier Transform of Anisotropic Hardy Spaces with Variable Exponents and Their Applications

abstract: Let $A$ be an expansive dilation on $\mathbb{R}^n$, and $p(\cdot):\mathbb{R}^n\rightarrow(0,\,\infty)$ be a variable exponent function satisfying the globally log-H\"{o}lder continuous condition. Let $\mathcal{H}^{p(\cdot)}_A({\mathbb {R}}^n)$ be the variable anisotropic Hardy space defined via the non-tangential grand maximal function. In this paper, the authors obtain that the Fourier transform of $f\in \mathcal{H}^{p(\cdot)}_A({\mathbb {R}}^n)$ coincides with a continuous function $F$ on $\mathbb{R}^n$ in the sense of tempered distributions. As applications, the authors further conclude a higher order convergence of the continuous function $F$ at the origin and then give a variant of the Hardy-Littlewood inequality in the setting of anisotropic Hardy spaces with variable exponents.


title: RSSI prediction using Machine Learning models

abstract: In this study, we present a method to predict the Received signal strength indication (RSSI) in an area of the base station. Traditional attenuated wave propagation models are often time consuming as well as computationally complex, depending on the unique factors of the medium. This study focuses on providing a solution to predict signal quality using coordinate values of many points in the considering area. We apply machine learning models such as linear regression, Support Vector Machine (SVM) or Decision tree model, to directly predict the RSSI of many points in the range of a base station without computing the complex parameters of the attenuated propagation model. The effectiveness of RSSI prediction was evaluated by mean square error (MSE) and mean absolute error (MAE). The stage of training and testing machine learning models in the research uses data that are the actual measurement results during the research process.


title: Short range vs long range dependence. An hyppothesis test based on Fractional Iterated Ornstein--Uhlenbeck processes

abstract: In this work, which is based on the family of Fractional Iterated Ornstein Uhlenbeck processes, we propose a new hypothesis test to contrast short memory versus long memory in time series. This family includes short memory and long memory processes, and has the ability to approximate a long memory processes by a short memory processes. Based on the asymptotic results of the estimators of its parameters, we will present the test and show how it can be implemented. Also, we will show a comparison with other tests widely used under both short memory and long memory scenarios. The main conclusion is that this new test is the one with best performance under the null hypothesis, and has the maximum power in some alternatives.


title: Bias of reconstructing the dark energy equation of state from the Pad$\acute{e}$ cosmography

abstract: Pad$\acute{e}$ cosmography has been widely used to probe the cosmic evolution and to investigate the nature of dark energy. In this paper, we find that the Pad$\acute{e}$ approximant can describe the cosmic evolution better than the standard cosmography, and if the luminosity distance $d_L(z)$ described by the Pad$\acute{e}$ approximant is used to reconstruct the dark energy equation of state $w(z)$, then the reconstructed $w(z)$ will approach a constant, i.e. $1/3$ or $0$, when the redshift is very high. This result is general since it is independent of the coefficients in the Pad$\acute{e}$ approximant and the value of the present dimensionless matter density parameter. This intrinsic character will bias the $w(z)$ reconstruction and lead to misconception of the property of dark energy. Therefore, one must exercise caution in investigating the property of dark energy from Pad$\acute{e}$ cosmography when the high redshift data, i.e. $z>2$, are included.


title: Continuous-Time Video Generation via Learning Motion Dynamics with Neural ODE

abstract: In order to perform unconditional video generation, we must learn the distribution of the real-world videos. In an effort to synthesize high-quality videos, various studies attempted to learn a mapping function between noise and videos, including recent efforts to separate motion distribution and appearance distribution. Previous methods, however, learn motion dynamics in discretized, fixed-interval timesteps, which is contrary to the continuous nature of motion of a physical body. In this paper, we propose a novel video generation approach that learns separate distributions for motion and appearance, the former modeled by neural ODE to learn natural motion dynamics. Specifically, we employ a two-stage approach where the first stage converts a noise vector to a sequence of keypoints in arbitrary frame rates, and the second stage synthesizes videos based on the given keypoints sequence and the appearance noise vector. Our model not only quantitatively outperforms recent baselines for video generation, but also demonstrates versatile functionality such as dynamic frame rate manipulation and motion transfer between two datasets, thus opening new doors to diverse video generation applications.


title: Nonlinear Transform Source-Channel Coding for Semantic Communications

abstract: In this paper, we propose a class of high-efficiency deep joint source-channel coding methods that can closely adapt to the source distribution under the nonlinear transform, it can be collected under the name nonlinear transform source-channel coding (NTSCC). In the considered model, the transmitter first learns a nonlinear analysis transform to map the source data into latent space, then transmits the latent representation to the receiver via deep joint source-channel coding. Our model incorporates the nonlinear transform as a strong prior to effectively extract the source semantic features and provide side information for source-channel coding. Unlike existing conventional deep joint source-channel coding methods, the proposed NTSCC essentially learns both the source latent representation and an entropy model as the prior on the latent representation. Accordingly, novel adaptive rate transmission and hyperprior-aided codec refinement mechanisms are developed to upgrade deep joint source-channel coding. The whole system design is formulated as an optimization problem whose goal is to minimize the end-to-end transmission rate-distortion performance under established perceptual quality metrics. Across test image sources with various resolutions, we find that the proposed NTSCC transmission method generally outperforms both the analog transmission using the standard deep joint source-channel coding and the classical separation-based digital transmission. Notably, the proposed NTSCC method can potentially support future semantic communications due to its content-aware ability and perceptual optimization goal.


title: Magnetically tunable Shubnikov-de Hass oscillations in MnBi2Te4

abstract: Shubnikov-de Hass oscillations are directly observed in undoped antiferromagnetic topological insulator MnBi2Te4. With increasing magnetic fields, the oscillation period decreases gradually in the magnetic transition from canted antiferromagnetism to ferromagnetism and then saturates in high magnetic fields, indicating the field-induced evolution of the band structure. From the analysis of the high-field oscillations, a nontrivial Berry phase and a small effective mass are extracted, in agreement with the predicted Weyl semimetal phase in ferromagnetic MnBi2Te4. Furthermore, rotating the magnetization of MnBi2Te4 can lead to a splitting of the high-field oscillations, which suggests the enhanced asymmetry of the Weyl cones in tilted fields. Therefore, the observation of these magnetically tunable quantum oscillations clearly demonstrates the indispensable role of field in tuning the band structure or physical properties of MnBi2Te4.


title: DRPN: Making CNN Dynamically Handle Scale Variation

abstract: Based on our observations of infrared targets, serious scale variation along within sequence frames has high-frequently occurred. In this paper, we propose a dynamic re-parameterization network (DRPN) to deal with the scale variation and balance the detection precision between small targets and large targets in infrared datasets. DRPN adopts the multiple branches with different sizes of convolution kernels and the dynamic convolution strategy. Multiple branches with different sizes of convolution kernels have different sizes of receptive fields. Dynamic convolution strategy makes DRPN adaptively weight multiple branches. DRPN can dynamically adjust the receptive field according to the scale variation of the target. Besides, in order to maintain effective inference in the test phase, the multi-branch structure is further converted to a single-branch structure via the re-parameterization technique after training. Extensive experiments on FLIR, KAIST, and InfraPlane datasets demonstrate the effectiveness of our proposed DRPN. The experimental results show that detectors using the proposed DRPN as the basic structure rather than SKNet or TridentNet obtained the best performances.


title: Upper Field Strength Limit of Fast Radio Bursts

abstract: Fast radio bursts (FRBs) are cosmological radio transients with unclear generation mechanism. Known characteristics such as their luminosity, duration, spectrum and repetition rate, etc. suggest that FRBs are powerful coherent radio signals at GHz frequencies, but the status of FRBs near source remain unknown. As an extreme astronomical event, FRBs should be accompanied by energy -- comparable or even more powerful x/{\gamma}-ray counterparts. Here, particle-in-cell simulations of ultra-strong GHz radio pulse interaction with GeV photons show that at 3*10^12V/cm field-strengths, quantum cascade can generate dense pair plasmas, which greatly dampen the radio pulse. Thus, in the presence of GeV photons in the source region, GHz radio pulses stronger than 3*10^12V/cm cannot escape. This result indicates an upper field-strength limit of FRB at the source.


title: Fermi arcs and pseudogap phase in a minimal microscopic model of $d$-wave superconductivity

abstract: We show conclusively that a pseudogap state can arise at $T > T_c$, for reasonable pairing interaction strength, from order parameter fluctuations in a two dimensional minimal model of $d$-wave superconductivity. The occurrence of the pseudogap requires neither strong correlation nor the presence of competing order. We study a model with attractive nearest neighbor interaction and establish our result using a combination of cluster based Monte Carlo for the order parameter field and a twisted-boundary scheme to compute the momentum-resolved spectral function. Apart from a dip in the density of states that characterizes the pseudogap, the momentum and frequency resolution on our effective lattice size $\sim 160 \times 160$ allows two major conclusions: (i)~at $T < T_c$, despite the presence of thermal phase fluctuations the superconductor has only nodal Fermi points while all non nodal points on the normal state Fermi surface show a two peak spectral function with a dip at $\omega =0$, and (ii)~for $T > T_c$ the Fermi points develops into arcs, characterized by a single quasiparticle peak, and the arcs connect up to recover the normal state Fermi surface at a temperature $T^* > T_c$. We show the variation of $T_c$ and $T^*$ with coupling strength and provide detailed spectral results at a coupling where $T^* \sim 1.5T_c$.


title: Mass spectra, wave functions and mixing effects of the $(bcq)$ baryon

abstract: Mass spectra and wave functions of the $J^P=\frac{1}{2}^+$ $(bcq)$ baryons are calculated by the relativistic Bethe-Salpeter equation\,(BSE) with considering the mixing effects between the $1^+$ and $0^+$ $(bc)$-diquarks inside. Based on the diquark picture, the three-body problem of baryons is transformed into two two-body problems. The BSE and wave functions of the $0^+$ diquark are given, and then solved numerically to obtain the effective mass spectra and form factors. Also we present the wave functions at zero point for the $(bc)$-diquark. Considering the obtained diquark form factors, the $(bcq)$ baryons are then described by the BSE as the bound state of a diquark and a light quark, where the interaction kernel includes the inner transitions between the $0^+$ and $1^+$ diquarks. The general wave function of the $\frac{1}{2}^+$ $(bcq)$ baryons is constructed and solved to obtain the corresponding mass spectra. Especially, by using the obtained wave functions, the mixing effects between $\Xi_{bc}(\Omega_{bc})$ and $\Xi_{bc}'(\Omega'_{bc})$ in ground states are computed and determined to be small ($\sim \!1\%$). The numerical results indicate that it is a good choice to take $\Xi_{bc}$ and $\Xi'_{bc}$ as the baryon states with the inside $(bc)$-diquarks occupying the definite spin.


title: Improved Efficiency for Cross-Arm Comparisons via Platform Designs

abstract: Though platform trials have been touted for their flexibility and streamlined use of trial resources, their statistical efficiency is not well understood. We fill this gap by establishing their greater efficiency for comparing the relative efficacy of multiple interventions over using several separate, two-arm trials, where the relative efficacy of an arbitrary pair of interventions is evaluated by contrasting their relative risks as compared to control. In theoretical and numerical studies, we demonstrate that the inference of such a contrast using data from a platform trial enjoys identical or better precision than using data from separate trials, even when the former enrolls substantially fewer participants. This benefit is attributed to the sharing of controls among interventions under contemporaneous randomization, which is a key feature of platform trials. We further provide a novel procedure for establishing the non-inferiority of a given intervention relative to the most efficacious of the other interventions under evaluation, where this procedure is adaptive in the sense that it need not be \textit{a priori} known which of these other interventions is most efficacious. Our numerical studies show that this testing procedure can attain substantially better power when the data arise from a platform trial rather than multiple separate trials. Our results are illustrated using data from two monoclonal antibody trials for the prevention of HIV.


title: Efficient method for accelerating line searches using a combined Schur complement domain decomposition and Born series expansions in photonic-based adjoint optimization

abstract: A line search in gradient-based optimization algorithm solves the problem of determining the optimal learning rate for a given gradient or search direction in a single iteration. For most problems, this is determined by evaluating different candidate learning rates to find the optimum, which can be expensive. Recent work has provided an efficient way to perform a line search with the use of the Shanks transformation of a Born series derived from the Lippman-Schwinger formalism. In this paper we show that the cost for performing such a line search can be further reduced with the use of the method of the Schur complement domain decomposition, which can lead to a 10-fold total speed-up resulting from the reduced number of iterations to convergence and reduced wall-clock time per iteration.


title: Generalizing Interactive Backpropagating Refinement for Dense Prediction

abstract: As deep neural networks become the state-of-the-art approach in the field of computer vision for dense prediction tasks, many methods have been developed for automatic estimation of the target outputs given the visual inputs. Although the estimation accuracy of the proposed automatic methods continues to improve, interactive refinement is oftentimes necessary for further correction. Recently, feature backpropagating refinement scheme (f-BRS) has been proposed for the task of interactive segmentation, which enables efficient optimization of a small set of auxiliary variables inserted into the pretrained network to produce object segmentation that better aligns with user inputs. However, the proposed auxiliary variables only contain channel-wise scale and bias, limiting the optimization to global refinement only. In this work, in order to generalize backpropagating refinement for a wide range of dense prediction tasks, we introduce a set of G-BRS (Generalized Backpropagating Refinement Scheme) layers that enable both global and localized refinement for the following tasks: interactive segmentation, semantic segmentation, image matting and monocular depth estimation. Experiments on SBD, Cityscapes, Mapillary Vista, Composition-1k and NYU-Depth-V2 show that our method can successfully generalize and significantly improve performance of existing pretrained state-of-the-art models with only a few clicks.


title: A deterministic-particle-based scheme for micro-macro viscoelastic flows

abstract: In this article, we introduce a new method for discretizing micro-macro models of dilute polymeric fluids based on deterministic particles. Our approach integrates a finite element discretization for the macroscopic fluid dynamic equation with a deterministic variational particle scheme for the microscopic Fokker-Planck equation. To address challenges arising from micro-macro coupling, we employ a discrete energetic variational approach to derive a coarse-grained micro-macro model with a particle approximation first and then develop a particle-FEM discretization for the coarse-grained model. The accuracy of our method is evaluated for a Hookean dumbbell model in a Couette flow by comparing the computed velocity field with existing analytical solutions. We also use our method to study nonlinear FENE dumbbell models in different scenarios, such as extension flow, pure shear flow, and lid-driven cavity flow. Numerical examples demonstrate that the proposed deterministic particle approach can accurately capture the various key rheological phenomena in the original FENE model, including hysteresis and $\delta$-function-like spike behavior in extension flows, velocity overshoot phenomenon in pure shear flow, symmetries breaking, vortex center shifting and vortices weakening in the lid-driven cavity flow, with a small number of particles.

categories: math.NA cond-mat.soft cs.NA physics.flu-dyn


title: Differentiated uniformization: A new method for inferring Markov chains on combinatorial state spaces including stochastic epidemic models

abstract: Motivation: We consider continuous-time Markov chains that describe the stochastic evolution of a dynamical system by a transition-rate matrix $Q$ which depends on a parameter $\theta$. Computing the probability distribution over states at time $t$ requires the matrix exponential $\exp(tQ)$, and inferring $\theta$ from data requires its derivative $\partial\exp\!(tQ)/\partial\theta$. Both are challenging to compute when the state space and hence the size of $Q$ is huge. This can happen when the state space consists of all combinations of the values of several interacting discrete variables. Often it is even impossible to store $Q$. However, when $Q$ can be written as a sum of tensor products, computing $\exp(tQ)$ becomes feasible by the uniformization method, which does not require explicit storage of $Q$. Results: Here we provide an analogous algorithm for computing $\partial\exp\!(tQ)/\partial\theta$, the differentiated uniformization method. We demonstrate our algorithm for the stochastic SIR model of epidemic spread, for which we show that $Q$ can be written as a sum of tensor products. We estimate monthly infection and recovery rates during the first wave of the COVID-19 pandemic in Austria and quantify their uncertainty in a full Bayesian analysis. Availability: Implementation and data are available at https://github.com/spang-lab/TenSIR.


title: Search for baryon and lepton number violating decays $D^{0}\to \bar{p}e^{+}$ and $D^{0}\to pe^{-}$

abstract: Using an electron-positron collision data sample corresponding to an integrated luminosity of 2.93~fb$^{-1}$ collected with the BESIII detector at a center-of-mass energy of 3.773 GeV, we search for the baryon and lepton number violating decays $D^{0}\to \bar{p}e^{+}$ and $D^{0}\to pe^{-}$. No obvious signals are found with the current statistics. The upper limits on the branching fractions for $D^{0}\to \bar{p}e^{+}$ and $D^{0}\to pe^{-}$ are set to be $1.2\times 10^{-6}$ and $2.2\times 10^{-6}$ at 90\% confidence level, respectively.


title: Minimizing Congestion for Balanced Dominators

abstract: A primary challenge in metagenomics is reconstructing individual microbial genomes from the mixture of short fragments created by sequencing. Recent work leverages the sparsity of the assembly graph to find $r$-dominating sets which enable rapid approximate queries through a dominator-centric graph partition. In this paper, we consider two problems related to reducing uncertainty and improving scalability in this setting. First, we observe that nodes with multiple closest dominators necessitate arbitrary tie-breaking in the existing pipeline. As such, we propose finding $\textit{sparse}$ dominating sets which minimize this effect via a new $\textit{congestion}$ parameter. We prove minimizing congestion is NP-hard, and give an $\mathcal{O}(\sqrt{\Delta^r})$ approximation algorithm, where $\Delta$ is the max degree. To improve scalability, the graph should be partitioned into uniformly sized pieces, subject to placing vertices with a closest dominator. This leads to $\textit{balanced neighborhood partitioning}$: given an $r$-dominating set, find a partition into connected subgraphs with optimal uniformity so that each vertex is co-assigned with some closest dominator. Using variance of piece sizes to measure uniformity, we show this problem is NP-hard iff $r$ is greater than $1$. We design and analyze several algorithms, including a polynomial-time approach which is exact when $r=1$ (and heuristic otherwise). We complement our theoretical results with computational experiments on a corpus of real-world networks showing sparse dominating sets lead to more balanced neighborhood partitionings. Further, on the metagenome $\textsf{HuSB1}$, our approach maintains high query containment and similarity while reducing piece size variance.


title: What are Attackers after on IoT Devices? An approach based on a multi-phased multi-faceted IoT honeypot ecosystem and data clustering

abstract: The growing number of Internet of Things (IoT) devices makes it imperative to be aware of the real-world threats they face in terms of cybersecurity. While honeypots have been historically used as decoy devices to help researchers/organizations gain a better understanding of the dynamic of threats on a network and their impact, IoT devices pose a unique challenge for this purpose due to the variety of devices and their physical connections. In this work, by observing real-world attackers' behavior in a low-interaction honeypot ecosystem, we (1) presented a new approach to creating a multi-phased, multi-faceted honeypot ecosystem, which gradually increases the sophistication of honeypots' interactions with adversaries, (2) designed and developed a low-interaction honeypot for cameras that allowed researchers to gain a deeper understanding of what attackers are targeting, and (3) devised an innovative data analytics method to identify the goals of adversaries. Our honeypots have been active for over three years. We were able to collect increasingly sophisticated attack data in each phase. Furthermore, our data analytics points to the fact that the vast majority of attack activities captured in the honeypots share significant similarity, and can be clustered and grouped to better understand the goals, patterns, and trends of IoT attacks in the wild.


title: A Linear Outer Approximation of Line Losses for DC-based Optimal Power Flow Problems

abstract: This paper proposes a novel and simple linear model to capture line losses for use in linearized DC models, such as optimal power flow (DC-OPF) and security-constrained economic dispatch (SCED). The \textit{Line Loss Outer Approximation} (LLOA) model implements an outer approximation of the line losses lazily and typically terminates in a small number of iterations. Experiments on large-scale power systems demonstrate the accuracy and computational efficiency of LLOA and contrast it with classical line loss approaches. The results seem to indicate that LLOA is a practical and useful model for real-world applications, providing a good tradeoff between accuracy, computational efficiency, and implementation simplicity. In particular, the LLOA method may have significant advantages compared to the traditional loss factor formulation for multi-period, stochastic optimization problems where good reference points may not be available. The paper also provides a comprehensive overview and evaluation of approximation methods for line losses.


title: A Linearized Boundary Control Method for the Acoustic Inverse Boundary Value Problem

abstract: We develop a linearized boundary control method for the inverse boundary value problem of determining a potential in the acoustic wave equation from the Neumann-to-Dirichlet map. When the linearization is at the zero potential, we derive a reconstruction formula based on the boundary control method and prove that it is of Lipschitz-type stability. When the linearization is at a nonzero potential, we prove that the problem is of H\"{o}lder-type stability. The proposed reconstruction formula is implemented and evaluated using several numerical experiments to validate its feasibility.


title: ACGNet: Action Complement Graph Network for Weakly-supervised Temporal Action Localization

abstract: Weakly-supervised temporal action localization (WTAL) in untrimmed videos has emerged as a practical but challenging task since only video-level labels are available. Existing approaches typically leverage off-the-shelf segment-level features, which suffer from spatial incompleteness and temporal incoherence, thus limiting their performance. In this paper, we tackle this problem from a new perspective by enhancing segment-level representations with a simple yet effective graph convolutional network, namely action complement graph network (ACGNet). It facilitates the current video segment to perceive spatial-temporal dependencies from others that potentially convey complementary clues, implicitly mitigating the negative effects caused by the two issues above. By this means, the segment-level features are more discriminative and robust to spatial-temporal variations, contributing to higher localization accuracies. More importantly, the proposed ACGNet works as a universal module that can be flexibly plugged into different WTAL frameworks, while maintaining the end-to-end training fashion. Extensive experiments are conducted on the THUMOS'14 and ActivityNet1.2 benchmarks, where the state-of-the-art results clearly demonstrate the superiority of the proposed approach.


title: Tree-informed Bayesian multi-source domain adaptation: cross-population probabilistic cause-of-death assignment using verbal autopsy

abstract: Determining causes of deaths (COD) occurred outside of civil registration and vital statistics systems is challenging. A technique called verbal autopsy (VA) is widely adopted to gather information on deaths in practice. A VA consists of interviewing relatives of a deceased person about symptoms of the deceased in the period leading to the death, often resulting in multivariate binary responses. While statistical methods have been devised for estimating the cause-specific mortality fractions (CSMFs) for a study population, continued expansion of VA to new populations (or "domains") necessitates approaches that recognize between-domain differences while capitalizing on potential similarities. In this paper, we propose such a domain-adaptive method that integrates external between-domain similarity information encoded by a pre-specified rooted weighted tree. Given a cause, we use latent class models to characterize the conditional distributions of the responses that may vary by domain. We specify a logistic stick-breaking Gaussian diffusion process prior along the tree for class mixing weights with node-specific spike-and-slab priors to pool information between the domains in a data-driven way. Posterior inference is conducted via a scalable variational Bayes algorithm. Simulation studies show that the domain adaptation enabled by the proposed method improves CSMF estimation and individual COD assignment. We also illustrate and evaluate the method using a validation data set. The paper concludes with a discussion on limitations and future directions.


title: Mott Insulator-Density Ordered Superfluid Transition and "Shamrock Transition" in a Frustrated Triangle Lattice

abstract: Density order is usually a consequence of the competition between long-range and short-range interactions. Here we report a density ordered superfluid emergent from a homogeneous Mott insulator due to the competition between frustrations and local interactions. This transition is found in a Bose-Hubbard model on a frustrated triangle lattice with an extra pairing term. Further, we find a quantum phase transition between two different density ordered superfluids, which is beyond the Landau-Ginzburg paradigm. Across this transition, a U(1) symmetry is emergent, while the symmetry in each density ordered superfluid is Z2*Z3. Because there emerges a shamrock-like degenerate ground state in parameter space, we call the transition "shamrock transition". Effective low energy theories are established for the two transitions mentioned above and we find their resemblance and differences with clock models.


title: Alternating knots do not admit cosmetic crossings

abstract: By examining the homology groups of a 4-manifold associated to an integral surgery on a knot $K$ in a rational homology 3-sphere $Y$ yielding a rational homology 3-sphere $Y^*$ with surgery dual knot $K^*$, we show that the subgroups generated by $[K]$ and $[K^*]$ in $H_1(Y)$ and $H_1(Y^*)$, respectively, have co-prime orders. We obtain an immediate corollary that, in conjunction with an argument of Lidman and Moore, proves the cosmetic crossing conjecture for knots whose branched double covers are Heegaard Floer L-spaces.


title: Ballooning in Spiders using Multiple Silk Threads

abstract: In this paper, three-dimensional numerical simulations of ballooning in spiders using multiple silk threads are performed using the discrete elastic rods method. The ballooning of spiders is hypothesised to be caused by the presence of the negative electric charge of the spider silk threads and the positive electric potential field in the earth's atmosphere. The numerical model presented here is first validated against experimental data from the open literature. After which, two cases are examined, in the first it is assumed that the electric charge is uniformly distributed along the threads while in the second, the electric charge is located at the thread tip. It is shown that the normalized terminal ballooning velocity, i.e. the velocity at which the spiders balloon after they reach steady-state, decrease linearly with the normalized lift force, especially for the tip located charge case. For the uniform electric charge case, this velocity shows a slightly weaker dependence on the normalized lift force. Moreover, it is shown in both cases that the normalized terminal ballooning velocity has no dependence on the normalized elastic bending stiffness of the threads and on the normalized viscous forces. Finally, the multi-thread bending process shows a three-dimensional conical sheet. Here we show that this behavior is caused by the Coulomb repelling forces owing to the threads electric charge which leads to dispersing the threads apart and thus avoid entanglement.


title: Generalized Few-Shot Semantic Segmentation: All You Need is Fine-Tuning

abstract: Generalized few-shot semantic segmentation was introduced to move beyond only evaluating few-shot segmentation models on novel classes to include testing their ability to remember base classes. While the current state-of-the-art approach is based on meta-learning, it performs poorly and saturates in learning after observing only a few shots. We propose the first fine-tuning solution, and demonstrate that it addresses the saturation problem while achieving state-of-the-art results on two datasets, PASCAL-5i and COCO-20i. We also show that it outperforms existing methods, whether fine-tuning multiple final layers or only the final layer. Finally, we present a triplet loss regularization that shows how to redistribute the balance of performance between novel and base categories so that there is a smaller gap between them.


title: Hybrid Modeling of Regional COVID-19 Transmission Dynamics in the U.S

abstract: The fast transmission rate of COVID-19 worldwide has made this virus the most important challenge of year 2020. Many mitigation policies have been imposed by the governments at different regional levels (country, state, county, and city) to stop the spread of this virus. Quantifying the effect of such mitigation strategies on the transmission and recovery rates, and predicting the rate of new daily cases are two crucial tasks. In this paper, we propose a hybrid modeling framework which not only accounts for such policies but also utilizes the spatial and temporal information to characterize the pattern of COVID-19 progression. Specifically, a piecewise susceptible-infected-recovered (SIR) model is developed while the dates at which the transmission/recover rates change significantly are defined as "break points" in this model. A novel and data-driven algorithm is designed to locate the break points using ideas from fused lasso and thresholding. In order to enhance the forecasting power and to describe additional temporal dependence among the daily number of cases, this model is further coupled with spatial smoothing covariates and vector auto-regressive (VAR) model. The proposed model is applied to several U.S. states and counties, and the results confirm the effect of "stay-at-home orders" and some states' early "re-openings" by detecting break points close to such events. Further, the model provided satisfactory short-term forecasts of the number of new daily cases at regional levels by utilizing the estimated spatio-temporal covariance structures. They were also better or on par with other proposed models in the literature, including flexible deep learning ones. Finally, selected theoretical results and empirical performance of the proposed methodology on synthetic data are reported which justify the good performance of the proposed method.


title: Analytic Bootstrap in 2D Boundary Conformal Field Theory: Towards Braneworld Holography

abstract: Recently, boundary conformal field theories (BCFTs) have attracted much attention in the context of quantum gravity. This is because a BCFT can be dual to gravity coupled to a heat bath CFT, known as the island model. On this background, it would be interesting to explore the duality between the boundary and the braneworld. However, this seems to be a challenging problem. The reason is because although there has been much study of rational BCFTs, there has been comparatively little study of irrational BCFTs, and irrational BCFTs are expected to be the boundary duals of the braneworlds. For this reason, we explore properties of boundary ingredients: the boundary primary spectrum, the boundary-boundary-boundary OPE coefficients and the bulk-boundary OPE coefficients. For this purpose, the conformal bootstrap is extremely useful. This is the first step in providing an understanding of BCFTs in the context of braneworld holography by using the conformal bootstrap. The techniques developed in this paper may be useful for further investigation of irrational BCFTs.


title: Learned ISTA with Error-based Thresholding for Adaptive Sparse Coding

abstract: Drawing on theoretical insights, we advocate an error-based thresholding (EBT) mechanism for learned ISTA (LISTA), which utilizes a function of the layer-wise reconstruction error to suggest a specific threshold for each observation in the shrinkage function of each layer. We show that the proposed EBT mechanism well disentangles the learnable parameters in the shrinkage functions from the reconstruction errors, endowing the obtained models with improved adaptivity to possible data variations. With rigorous analyses, we further show that the proposed EBT also leads to a faster convergence on the basis of LISTA or its variants, in addition to its higher adaptivity. Extensive experimental results confirm our theoretical analyses and verify the effectiveness of our methods.


title: Shared Frailty Models Based on Cancer Data

abstract: Traditional survival analysis techniques focus on the occurrence of failures over the time. During analysis of such events, ignoring the related unobserved covariates or heterogeneity involved in data sample may leads us to adverse consequences. In this context, frailty models are the viable choice to investigate the effect of the unobserved covariates. In this article, we assume that frailty acts multiplicatively to hazard rate. We propose inverse Gaussian (IG) and generalized Lindley (GL) shared frailty models with generalized Weibull (GW) as baseline distribution in order to analyze the unobserved heterogeneity. To estimate the parameters in models, Bayesian paradigm of Markov Chain Monte Carlo technique has been proposed. Model selection criteria have been used for the comparison of models. Three different cancer data sets have been analyzed using the shared frailty models. Better models have been suggested for the data sets.


title: Lower Bounds for Sparse Oblivious Subspace Embeddings

abstract: An oblivious subspace embedding (OSE), characterized by parameters $m,n,d,\epsilon,\delta$, is a random matrix $\Pi\in \mathbb{R}^{m\times n}$ such that for any $d$-dimensional subspace $T\subseteq \mathbb{R}^n$, $\Pr_\Pi[\forall x\in T, (1-\epsilon)\|x\|_2 \leq \|\Pi x\|_2\leq (1+\epsilon)\|x\|_2] \geq 1-\delta$. For $\epsilon$ and $\delta$ at most a small constant, we show that any OSE with one nonzero entry in each column must satisfy that $m = \Omega(d^2/(\epsilon^2\delta))$, establishing the optimality of the classical Count-Sketch matrix. When an OSE has $1/(9\epsilon)$ nonzero entries in each column, we show it must hold that $m = \Omega(\epsilon^{O(\delta)} d^2)$, improving on the previous $\Omega(\epsilon^2 d^2)$ lower bound due to Nelson and Nguyen (ICALP 2014).


title: Mapping industrial poultry operations at scale with deep learning and aerial imagery

abstract: Concentrated Animal Feeding Operations (CAFOs) pose serious risks to air, water, and public health, but have proven to be challenging to regulate. The U.S. Government Accountability Office notes that a basic challenge is the lack of comprehensive location information on CAFOs. We use the USDA's National Agricultural Imagery Program (NAIP) 1m/pixel aerial imagery to detect poultry CAFOs across the continental United States. We train convolutional neural network (CNN) models to identify individual poultry barns and apply the best performing model to over 42 TB of imagery to create the first national, open-source dataset of poultry CAFOs. We validate the model predictions against held-out validation set on poultry CAFO facility locations from 10 hand-labeled counties in California and demonstrate that this approach has significant potential to fill gaps in environmental monitoring.


title: Interpreting how nonlinear diffusion affects the fate of bistable populations using a discrete modelling framework

abstract: Understanding whether a population will survive and flourish or become extinct is a central question in population biology. One way of exploring this question is to study population dynamics using reaction-diffusion equations, where migration is usually represented as a linear diffusion term, and birth-death is represented with a bistable source term. While linear diffusion is most commonly employed to study migration, there are several limitations of this approach, such as the inability of linear diffusion-based models to predict a well-defined population front. One way to overcome this is to generalise the constant diffusivity, $D$, to a nonlinear diffusivity function $D(C)$, where $C>0$ is the density. While it has been formally established that the choice of $D(C)$ affects long-term survival or extinction of a bistable population, working solely in a classical continuum framework makes it difficult to understand precisely how the choice of $D(C)$ affects survival or extinction. Here, we address this question by working with a simple discrete simulation model that is easy to interpret. The continuum limit of the discrete model is a nonlinear reaction-diffusion equation, where the flux involves a nonlinear diffusion term and the source term is given by the strong Allee effect bistable model. We study population extinction/survival using this very intuitive discrete framework together with numerical solutions of the reaction-diffusion continuum limit equation. This approach provides clear insight into how the choice of $D(C)$ either encourages or suppresses population extinction relative to the classical linear diffusion model.


title: Narrowband searches for continuous and long-duration transient gravitational waves from known pulsars in the LIGO-Virgo third observing run

abstract: Isolated neutron stars that are asymmetric with respect to their spin axis are possible sources of detectable continuous gravitational waves. This paper presents a fully-coherent search for such signals from eighteen pulsars in data from LIGO and Virgo's third observing run (O3). For known pulsars, efficient and sensitive matched-filter searches can be carried out if one assumes the gravitational radiation is phase-locked to the electromagnetic emission. In the search presented here, we relax this assumption and allow the frequency and frequency time-derivative of the gravitational waves to vary in a small range around those inferred from electromagnetic observations. We find no evidence for continuous gravitational waves, and set upper limits on the strain amplitude for each target. These limits are more constraining for seven of the targets than the spin-down limit defined by ascribing all rotational energy loss to gravitational radiation. In an additional search we look in O3 data for long-duration (hours-months) transient gravitational waves in the aftermath of pulsar glitches for six targets with a total of nine glitches. We report two marginal outliers from this search, but find no clear evidence for such emission either. The resulting duration-dependent strain upper limits do not surpass indirect energy constraints for any of these targets.


title: Regularizing End-to-End Speech Translation with Triangular Decomposition Agreement

abstract: End-to-end speech-to-text translation (E2E-ST) is becoming increasingly popular due to the potential of its less error propagation, lower latency, and fewer parameters. Given the triplet training corpus $\langle speech, transcription, translation\rangle$, the conventional high-quality E2E-ST system leverages the $\langle speech, transcription\rangle$ pair to pre-train the model and then utilizes the $\langle speech, translation\rangle$ pair to optimize it further. However, this process only involves two-tuple data at each stage, and this loose coupling fails to fully exploit the association between triplet data. In this paper, we attempt to model the joint probability of transcription and translation based on the speech input to directly leverage such triplet data. Based on that, we propose a novel regularization method for model training to improve the agreement of dual-path decomposition within triplet data, which should be equal in theory. To achieve this goal, we introduce two Kullback-Leibler divergence regularization terms into the model training objective to reduce the mismatch between output probabilities of dual-path. Then the well-trained model can be naturally transformed as the E2E-ST models by the pre-defined early stop tag. Experiments on the MuST-C benchmark demonstrate that our proposed approach significantly outperforms state-of-the-art E2E-ST baselines on all 8 language pairs, while achieving better performance in the automatic speech recognition task. Our code is open-sourced at https://github.com/duyichao/E2E-ST-TDA.


title: Expansion-Squeeze-Excitation Fusion Network for Elderly Activity Recognition

abstract: This work focuses on the task of elderly activity recognition, which is a challenging task due to the existence of individual actions and human-object interactions in elderly activities. Thus, we attempt to effectively aggregate the discriminative information of actions and interactions from both RGB videos and skeleton sequences by attentively fusing multi-modal features. Recently, some nonlinear multi-modal fusion approaches are proposed by utilizing nonlinear attention mechanism that is extended from Squeeze-and-Excitation Networks (SENet). Inspired by this, we propose a novel Expansion-Squeeze-Excitation Fusion Network (ESE-FN) to effectively address the problem of elderly activity recognition, which learns modal and channel-wise Expansion-Squeeze-Excitation (ESE) attentions for attentively fusing the multi-modal features in the modal and channel-wise ways. Furthermore, we design a new Multi-modal Loss (ML) to keep the consistency between the single-modal features and the fused multi-modal features by adding the penalty of difference between the minimum prediction losses on single modalities and the prediction loss on the fused modality. Finally, we conduct experiments on a largest-scale elderly activity dataset, i.e., ETRI-Activity3D (including 110,000+ videos, and 50+ categories), to demonstrate that the proposed ESE-FN achieves the best accuracy compared with the state-of-the-art methods. In addition, more extensive experimental results show that the proposed ESE-FN is also comparable to the other methods in terms of normal action recognition task.


title: Learning in Random Utility Models Via Online Decision Problems

abstract: This paper studies the Random Utility Model (RUM) in a repeated stochastic choice situation, in which the decision maker is imperfectly informed about the payoffs of each available alternative. We develop a gradient-based learning algorithm by embedding the RUM into an online decision problem. We show that a large class of RUMs are Hannan consistent (\citet{Hahn1957}); that is, the average difference between the expected payoffs generated by a RUM and that of the best-fixed policy in hindsight goes to zero as the number of periods increase. In addition, we show that our gradient-based algorithm is equivalent to the Follow the Regularized Leader (FTRL) algorithm, which is widely used in the machine learning literature to model learning in repeated stochastic choice problems. Thus, we provide an economically grounded optimization framework to the FTRL algorithm. Finally, we apply our framework to study recency bias, no-regret learning in normal form games, and prediction markets.


title: Detection of Anticipatory Dynamics Between a Pair of Zebrafish

abstract: Trajectories from a pair of interacting zebrafish are used to test for the existence of anticipatory dynamics in natural systems. Anticipatory dynamics (AD) is unusual in that causal events are not necessarily ordered by their temporal order. However, their causal order can still be established if the direction of information flow (DIF) is known. In order to obtain DIF between trajectories of the two fish, we have made use of the difference of the transfer entropy between the trajectories with a history length established by experiments with known DIF. Our experimental results indicate that AD can be observed much more often in fish pairs of different genders. The use of DIF to determine causal order is further verified by the simulation of two chaotic Lorenz oscillators with anticipatory coupling; mimicking the interaction between the fish. Our simulation results further suggest that the two fish are interacting with their own internal dynamics, not by adaptation.


title: Asymptotic behavior of bifurcation curves of one-dimensional nonlocal elliptic equations

abstract: We study the one-dimensional nonlocal elliptic equation \begin{eqnarray*} -\left(\int_0^1 \vert u(x)\vert^p dx + b\right)^q u''(x) &=& \lambda u(x)^p, \quad x \in I:= (0,1), \ u(x) > 0, \ x\in I, \\ u(0) &=& u(1) = 0, \end{eqnarray*} where $b \ge 0, p \ge 1, q > 1 - \frac{1}{p}$ are given constants and $\lambda > 0$ is a bifurcation parameter. We establish the global behavior of bifurcation diagrams and precise asymptotic formulas for $u_\lambda(x)$ as $\lambda \to \infty$.


title: Efficient Estimation of the Maximal Association between Multiple Predictors and a Survival Outcome

abstract: This paper develops a new approach to post-selection inference for screening high-dimensional predictors of survival outcomes. Post-selection inference for right-censored outcome data has been investigated in the literature, but much remains to be done to make the methods both reliable and computationally-scalable in high-dimensions. Machine learning tools are commonly used to provide {\it predictions} of survival outcomes, but the estimated effect of a selected predictor suffers from confirmation bias unless the selection is taken into account. The new approach involves construction of semi-parametrically efficient estimators of the linear association between the predictors and the survival outcome, which are used to build a test statistic for detecting the presence of an association between any of the predictors and the outcome. Further, a stabilization technique reminiscent of bagging allows a normal calibration for the resulting test statistic, which enables the construction of confidence intervals for the maximal association between predictors and the outcome and also greatly reduces computational cost. Theoretical results show that this testing procedure is valid even when the number of predictors grows superpolynomially with sample size, and our simulations support that this asymptotic guarantee is indicative the performance of the test at moderate sample sizes. The new approach is applied to the problem of identifying patterns in viral gene expression associated with the potency of an antiviral drug.


title: Hyperscaling analysis of a disorder-induced ferromagnetic quantum critical point in Ni$_{1-x}$Rh$_{x}$ with $x = 0.375$

abstract: Here we report on hyperscaling analysis on thermodynamic measurements as a function of temperature and magnetic field for Ni$_{1-x}$Rh$_{x}$ with $x = 0.375$ where a ferromagnetic quantum critical point has been recently identified [Phys. Rev. Lett. $\textbf{124}$, 117203 (2020)]. The obtained critical exponents agree well with the theory proposed by Belitz, Kirkpatrick, and Vojta for disorder tuned quantum critical point in the preasymptotic region.


title: Thermal Emission and Scattering by Aligned Grains: Plane-Parallel Model and Application to Multiwavelength Polarization of the HL Tau Disk

abstract: Telescopes are now able to resolve dust polarization across circumstellar disks at multiple wavelengths, allowing the study of the polarization spectrum. Most disks show clear evidence of dust scattering through their unidirectional polarization pattern typically at the shorter wavelength of $\sim 870 \mu$m. However, certain disks show an elliptical pattern at $\sim 3$mm, which is likely due to aligned grains. With HL Tau, its polarization pattern at $\sim 1.3$mm shows a transition between the two patterns making it the first example to reveal such transition. We use the T-matrix method to model elongated dust grains and properly treat scattering of aligned non-spherical grains with a plane-parallel slab model. We demonstrate that a change in optical depth can naturally explain the polarization transition of HL Tau. At low optical depths, the thermal polarization dominates, while at high optical depths, dichroic extinction effectively takes out the thermal polarization and scattering polarization dominates. Motivated by results from the plane-parallel slab, we develop a simple technique to disentangle thermal polarization of the aligned grains $T_{0}$ and polarization due to scattering $S$ using the azimuthal variation of the polarization fraction. We find that, with increasing wavelength, the fractional polarization spectrum of the scattering component $S$ decreases, while the thermal component $T_{0}$ increases, which is expected since the optical depth decreases. We find several other sources similar to HL Tau that can be explained by azimuthally aligned scattering prolate grains when including optical depth effects. In addition, we explore how spirally aligned grains with scattering can appear in polarization images.


title: X-ray dark-field and phase retrieval without optics, via the Fokker-Planck equation

abstract: Emerging methods of x-ray imaging that capture phase and dark-field effects are equipping medicine with complementary sensitivity to conventional radiography. These methods are being applied over a wide range of scales, from virtual histology to clinical chest imaging, and typically require the introduction of optics such as gratings. Here, we consider extracting x-ray phase and dark-field signals from bright-field images collected using nothing more than a coherent x-ray source and detector. Our approach is based on the Fokker--Planck equation for paraxial imaging, which is the diffusive generalization of the transport-of-intensity equation. Specifically, we utilize the Fokker--Planck equation in the context of propagation-based phase-contrast imaging, where we show that two intensity images are sufficient for successful retrieval of the projected thickness and dark-field signals associated with the sample. We show the results of our algorithm using both a simulated dataset and an experimental dataset. These demonstrate that the x-ray dark-field signal can be extracted from propagation-based images, and that x-ray phase can be retrieved with better spatial resolution when dark-field effects are taken into account. We anticipate the proposed algorithm will be of benefit in biomedical imaging, industrial settings, and other non-invasive imaging applications.


title: Continuity of the Spectrum of Dirac Operators of Spectral Triples for the Spectral Propinquity

abstract: The spectral propinquity is a distance, up to unitary equivalence, on the class of metric spectral triples. We prove in this paper that if a sequence of metric spectral triples converges for the propinquity, then the spectra of the Dirac operators for these triples do converge to the spectrum of the Dirac operator at the limit. We obtain this result by first proving that, in an appropriate sense induced by some natural metric, the bounded, continuous functional calculi defined by the Dirac operators also converge. As an application of our work, we see, in particular, that action functionals of a wide class of metric spectral triples are continuous for the spectral propinquity, which clearly connects convergence for the spectral propinquity with the applications of noncommutative geometry to mathematical physics. This fact is a consequence of results on the continuity of the multiplicity of eigenvalues of Dirac operators. In particular, we formalize convergence of adjoinable operators of different C*-correspondences, endowed with appropriate quantum metric data.