title: Class-Incremental Learning with Repetition

abstract: Real-world data streams naturally include the repetition of previous concepts. From a Continual Learning (CL) perspective, repetition is a property of the environment and, unlike replay, cannot be controlled by the agent. Nowadays, the Class-Incremental (CI) scenario represents the leading test-bed for assessing and comparing CL strategies. This scenario type is very easy to use, but it never allows revisiting previously seen classes, thus completely neglecting the role of repetition. We focus on the family of Class-Incremental with Repetition (CIR) scenario, where repetition is embedded in the definition of the stream. We propose two stochastic stream generators that produce a wide range of CIR streams starting from a single dataset and a few interpretable control parameters. We conduct the first comprehensive evaluation of repetition in CL by studying the behavior of existing CL strategies under different CIR streams. We then present a novel replay strategy that exploits repetition and counteracts the natural imbalance present in the stream. On both CIFAR100 and TinyImageNet, our strategy outperforms other replay approaches, which are not designed for environments with repetition.


title: Automating Knowledge-Driven Model Recommendation: Methodology, Evaluation, and Key Challenges

abstract: There is significant interest in using existing repositories of biological entities, relationships, and models to automate biological model assembly and extension. Current methods aggregate human-curated biological information into executable, simulatable models, but these models do not resemble human curated models and do not recapitulate experimental results. Here, we outline the process of automated model assembly and extension, while demonstrating it on both synthetic models and human-curated models of biological signaling networks. We begin with an iterative, greedy, and combinatoric approach to automated assembly and demonstrate the key difficulties inherent to contextless assembly. We publicly release the software used in this paper to enable further exploration of this problem.


title: Smigoc's glue for universal realizability in the left half-plane

abstract: A list of A list {\Lambda} of complex numbers is said to be realizable if it is the spectrum of a nonnegative matrix. {\Lambda} is said to be universally realizable (UR) if it is realizable for each possible Jordan canonical form allowed by {\Lambda}. In this paper, using companion matrices and applying a procedure by \v{S}migoc, is provides a sufficient condition for the universal realizability of left half-plane spectra. It is also shown how the effect of adding a negative real number to a not UR left half-plane list of complex numbers, makes the new list UR, and a family of left half-plane lists that are UR is characterized.


title: Distributional outcome regression via quantile functions and its application to modelling continuously monitored heart rate and physical activity

abstract: Modern clinical and epidemiological studies widely employ wearables to record parallel streams of real-time data on human physiology and behavior. With recent advances in distributional data analysis, these high-frequency data are now often treated as distributional observations resulting in novel regression settings. Motivated by these modelling setups, we develop a distributional outcome regression via quantile functions (DORQF) that expands existing literature with three key contributions: i) handling both scalar and distributional predictors, ii) ensuring jointly monotone regression structure without enforcing monotonicity on individual functional regression coefficients, iii) providing statistical inference via asymptotic projection-based joint confidence bands and a statistical test of global significance to quantify uncertainty of the estimated functional regression coefficients. The method is motivated by and applied to Actiheart component of Baltimore Longitudinal Study of Aging that collected one week of minute-level heart rate (HR) and physical activity (PA) data on 781 older adults to gain deeper understanding of age-related changes in daily life heart rate reserve, defined as a distribution of daily HR, while accounting for daily distribution of physical activity, age, gender, and body composition. Intriguingly, the results provide novel insights in epidemiology of daily life heart rate reserve.


title: Effect of different selected food processing techniques on control release kinetics of encapsulated lycopene in simulated human gastrointestinal tract

abstract: Nanoencapsulation has become a widespread technique to improve the bioavailability of bioactive compounds like lycopene. Consumption of lycopene rich foods is effective in preventing cancer, diabetes, and cardiovascular diseases due to its strong, oxygen-quenching ability. But the functional activity of lycopene is compromised by light, oxygen, and heat. A biodegradable polymer such as PLA or PLGA is the most effective carrier for encapsulating lycopene due to its excellent biodegradability, biocompatibility, and nontoxigenic effect on the human metabolic system. Hence, the primary objective of this study was to evaluate the effect of pasteurization on bioaccessibility and control release kinetics of encapsulated lycopene nanoparticles in invitro human GIT. In the first objective, sonication time, surfactant, and polymer concentration were considered as the three factors to synthesize polymeric lycopene nanoparticles (LNP) whereas in the second objective, type of pasteurization, encapsulation, and juice concentration were three factors. In Objectives 4 and 5, the type of encapsulation and pasteurization, and digestion time were evaluated as the three factors used to evaluate the in-vitro bioaccessibility of encapsulated lycopene NP. The study evidenced that encapsulation improved the lycopene's bioaccessibility by 70% and more than 60% for conventional pasteurized (CP) and microwave pasteurized (MP) nanoemulsions, respectively without compromising the physicochemical properties of both PLA and PLGA-lycopene nanoparticles. The invitro bioaccessibility study also showed that CP reduced the functional activity of PLALNP by 20% whereas MP had no significant effect on the bioaccessibility of PLA LNP. It is unlikely that the PLGA LNP was more sensitive against MP nanoemulsion than the CP. In conclusion, it was shown that the PLA LNP treated with MP provided the highest bioaccessibility.


title: Causal Bandits without Graph Learning

abstract: We study the causal bandit problem when the causal graph is unknown and develop an efficient algorithm for finding the parent node of the reward node using atomic interventions. We derive the exact equation for the expected number of interventions performed by the algorithm and show that under certain graphical conditions it could perform either logarithmically fast or, under more general assumptions, slower but still sublinearly in the number of variables. We formally show that our algorithm is optimal as it meets the universal lower bound we establish for any algorithm that performs atomic interventions. Finally, we extend our algorithm to the case when the reward node has multiple parents. Using this algorithm together with a standard algorithm from bandit literature leads to improved regret bounds.


title: A Hybrid Deep Neural Operator/Finite Element Method for Ice-Sheet Modeling

abstract: One of the most challenging and consequential problems in climate modeling is to provide probabilistic projections of sea level rise. A large part of the uncertainty of sea level projections is due to uncertainty in ice sheet dynamics. At the moment, accurate quantification of the uncertainty is hindered by the cost of ice sheet computational models. In this work, we develop a hybrid approach to approximate existing ice sheet computational models at a fraction of their cost. Our approach consists of replacing the finite element model for the momentum equations for the ice velocity, the most expensive part of an ice sheet model, with a Deep Operator Network, while retaining a classic finite element discretization for the evolution of the ice thickness. We show that the resulting hybrid model is very accurate and it is an order of magnitude faster than the traditional finite element model. Further, a distinctive feature of the proposed model compared to other neural network approaches, is that it can handle high-dimensional parameter spaces (parameter fields) such as the basal friction at the bed of the glacier, and can therefore be used for generating samples for uncertainty quantification. We study the impact of hyper-parameters, number of unknowns and correlation length of the parameter distribution on the training and accuracy of the Deep Operator Network on a synthetic ice sheet model. We then target the evolution of the Humboldt glacier in Greenland and show that our hybrid model can provide accurate statistics of the glacier mass loss and can be effectively used to accelerate the quantification of uncertainty.

categories: physics.comp-ph cs.LG physics.ao-ph physics.geo-ph


title: Detecting Pump&Dump Stock Market Manipulation from Online Forums

abstract: The intersection of social media, low-cost trading platforms, and naive investors has created an ideal situation for information-based market manipulations, especially pump&dumps. Manipulators accumulate small-cap stocks, disseminate false information on social media to inflate their price, and sell at the peak. We collect a dataset of stocks whose price and volume profiles have the characteristic shape of a pump&dump, and social media posts for those same stocks that match the timing of the initial price rises. From these we build predictive models for pump&dump events based on the language used in the social media posts. There are multiple difficulties: not every post will cause the intended market reaction, some pump&dump events may be triggered by posts in other forums, and there may be accidental confluences of post timing and market movements. Nevertheless, our best model achieves a prediction accuracy of 85% and an F1-score of 62%. Such a tool can provide early warning to investors and regulators that a pump&dump may be underway.


title: A Quantum Monte Carlo study of the structural, energetic, and magnetic properties of two-dimensional (2D) H and T phase VSe$_2$

abstract: Previous works have controversially claimed near-room temperature ferromagnetism in two-dimensional (2D) VSe$_2$, with conflicting results throughout the literature. These discrepancies in magnetic properties between both phases (T and H phase) of 2D VSe$_2$ are most likely due to the structural parameters being coupled to the magnetic properties. Specifically, both phases have a close lattice match and similar total energies, which makes it difficult to determine which phase is being observed experimentally. In this study, we used a combination of density functional theory (DFT), highly accurate diffusion Monte Carlo (DMC) and a surrogate Hessian line-search optimization technique to resolve the previously reported discrepancy in structural parameters and relative phase stability. With DMC accuracy, we determined the freestanding geometry of both phases and constructed a phase diagram. Our findings demonstrate the successes of the DMC method coupled with the surrogate Hessian structural optimization technique when applied to a 2D magnetic system.


title: Discriminative Entropy Clustering and its Relation to K-means and SVM

abstract: Maximization of mutual information between the model's input and output is formally related to "decisiveness" and "fairness" of the softmax predictions, motivating these unsupervised entropy-based criteria for clustering. First, in the context of linear softmax models, we discuss some general properties of entropy-based clustering. Disproving some earlier claims, we point out fundamental differences with K-means. On the other hand, we prove the margin maximizing property for decisiveness establishing a relation to SVM-based clustering. Second, we propose a new self-labeling formulation of entropy clustering for general softmax models. The pseudo-labels are introduced as auxiliary variables "splitting" the fairness and decisiveness. The derived self-labeling loss includes the reverse cross-entropy robust to pseudo-label errors and allows an efficient EM solver for pseudo-labels. Our algorithm improves the state of the art on several standard benchmarks for deep clustering.


title: Beyond Arabic: Software for Perso-Arabic Script Manipulation

abstract: This paper presents an open-source software library that provides a set of finite-state transducer (FST) components and corresponding utilities for manipulating the writing systems of languages that use the Perso-Arabic script. The operations include various levels of script normalization, including visual invariance-preserving operations that subsume and go beyond the standard Unicode normalization forms, as well as transformations that modify the visual appearance of characters in accordance with the regional orthographies for eleven contemporary languages from diverse language families. The library also provides simple FST-based romanization and transliteration. We additionally attempt to formalize the typology of Perso-Arabic characters by providing one-to-many mappings from Unicode code points to the languages that use them. While our work focuses on the Arabic script diaspora rather than Arabic itself, this approach could be adopted for any language that uses the Arabic script, thus providing a unified framework for treating a script family used by close to a billion people.


title: SN1987A neutrino burst: limits on flavor conversion

abstract: In this paper, we revisit the SN1987A neutrino data to see its constraints on flavor conversion. We are motivated by the fact that most works that analyze this data consider a specific conversion mechanism, such as the MSW (Mikheyev-Smirnov-Wolfenstein) effect, although flavor conversion is still an open question in supernovae due to the presence of neutrino-neutrino interactions. In our analysis, instead of considering a specific conversion mechanism, we let the electron antineutrino survival probability $P_{\overline{e}\overline{e}}$ be a free parameter. We fit the data from Kamiokande-II, Baksan, and IMB detected spectrum with two classes of models: time-integrated and time-dependent. For the time-integrated model, it is not possible to put limits above $1\sigma$ (68% confidence level) on the survival probability. The same happens for the time-dependent model when cooling is the only mechanism of antineutrino emission. However, for models considering an accretion phase, $P_{\overline{e}\overline{e}}\sim0$ is strongly rejected, showing a preference for the existence of an accretion component in the detected antineutrino flux, and a preference for normal mass ordering when only the MSW is present.


title: DBGDGM: Dynamic Brain Graph Deep Generative Model

abstract: Graphs are a natural representation of brain activity derived from functional magnetic imaging (fMRI) data. It is well known that clusters of anatomical brain regions, known as functional connectivity networks (FCNs), encode temporal relationships which can serve as useful biomarkers for understanding brain function and dysfunction. Previous works, however, ignore the temporal dynamics of the brain and focus on static graphs. In this paper, we propose a dynamic brain graph deep generative model (DBGDGM) which simultaneously clusters brain regions into temporally evolving communities and learns dynamic unsupervised node embeddings. Specifically, DBGDGM represents brain graph nodes as embeddings sampled from a distribution over communities that evolve over time. We parameterise this community distribution using neural networks that learn from subject and node embeddings as well as past community assignments. Experiments demonstrate DBGDGM outperforms baselines in graph generation, dynamic link prediction, and is comparable for graph classification. Finally, an analysis of the learnt community distributions reveals overlap with known FCNs reported in neuroscience literature.


title: Self-testing composite measurements and bound entangled state in a unified framework

abstract: Within the quantum networks scenario, we introduce a single scheme allowing to certify three different types of composite projective measurements acting on a three-qubit Hilbert space: one constructed from genuinely entangled GHZ-like states, one constructed from fully product vectors that exhibit the phenomenon of nonlocality without entanglement (NLWE), and a hybrid measurement obtained from an unextendible product basis (UPB). Noticeably, we certify a basis exhibiting NLWE in the smallest dimension capable of supporting this phenomenon. On the other hand, the possibility of certification of a measurement obtained from a UPB has an interesting implication that one can also self-test a bound entangled state in the considered quantum network. Such a possibility does not seem to exist in the standard Bell scenario.


title: Automatic differentiation as an effective tool in Electrical Impedance Tomography

abstract: Determining physical properties inside an object without access to direct measurements of target regions can be formulated as a specific type of \textit{inverse problem}. One of such problems is applied in \textit{Electrical Impedance Tomography} (EIT). In general, EIT can be posed as a minimization problem and solved by iterative methods, which require knowledge of derivatives of the objective function. In practice, this can be challenging because analytical closed-form solutions for them are hard to derive and implement efficiently. In this paper, we study the effectiveness of \textit{automatic differentiation (AD)} to solve EIT in a minimization framework. We devise a case study where we compare solutions of the inverse problem obtained with AD methods and with the manually-derived formulation of the derivative against the true solution. Furthermore, we study the viability of AD for large scale inverse problems by checking the memory and load requirements of AD as the resolution of the model increases. With powerful infrastructure, AD can pave the way for faster and simpler inverse solvers and provide better results than classical methods.


title: Anisotropic Electron Heating in an Electron Cyclotron Resonance Thruster with Magnetic Nozzle

abstract: In a grid-less Electron Cyclotron Resonance (ECR) plasma thruster with a diverging magnetic nozzle, the magnitude of the ambipolar field accelerating the positive ions depends of the perpendicular energy gained by the electrons. This work investigates the heating of the electrons by electromagnetic waves, taking their bouncing motion into account in a confining well formed by the magnetic mirror force and the electrostatic potential of the thruster. An electromagnetic Particle-In-Cell (PIC) code is used to simulate the plasma in a magnetic field tube. The code's Maxwell solver is based on a semi-Lagrangian scheme known as the Constrained Interpolation Profile (CIP) which enables larger time steps. The results show that anisotropic plasma heating takes place exclusively inside the coaxial chamber, along a Doppler-broadened zone. It is also shown that a trapped population of electrons with a larger perpendicular energy exists in the plume.


title: Accretion of sub-stellar companions as the origin of chemical abundance inhomogeneities in globular clusters

abstract: Globular clusters exhibit abundance variations, defining `multiple populations', which have prompted a protracted search for their origin. Properties requiring explanation include: the high fraction of polluted stars ($\sim 40{-}90$~percent, correlated with cluster mass), the absence of pollution in young clusters and the lower pollution rate with binarity and distance from the cluster centre. We present a novel mechanism for late delivery of pollutants into stars via accretion of sub-stellar companions. In this scenario, stars move through a medium polluted with AGB and massive star ejecta, accreting material to produce companions with typical mass ratio $q\sim 0.1$. These companions undergo eccentricity excitation due to dynamical perturbations by passing stars, culminating in a merger with their host star. The accretion of the companion alters surface abundances via injected pollutant. Alongside other self-enrichment models, the companion accretion model can explain the dilution of pollutant and correlation with intra-cluster location. The model also explains the ubiquity and discreteness of the populations and correlations of enrichment rates with cluster mass, cluster age and stellar binarity. Abundance variations in some clusters can be broadly reproduced using AGB and massive binary ejecta abundances from the literature. In other clusters, some high companion mass ratios ($q\gtrsim 1$) are required. In these cases, the available mass budget necessitates a variable degree of mixing of the polluted material with the primary star, deviations from model ejecta abundances or mixing of internal burning products. We highlight the avenues of further investigation which are required to explore some of the key processes invoked in this model.


title: A massive quiescent galaxy at redshift 4.658

abstract: The extremely rapid assembly of the earliest galaxies during the first billion years of cosmic history is a major challenge for our understanding of galaxy formation physics. The advent of JWST has exacerbated this issue by confirming the existence of galaxies in significant numbers as early as the first few hundred million years. Perhaps even more surprisingly, in some galaxies, this initial highly efficient star formation rapidly shuts down, or quenches, giving rise to massive quiescent galaxies as little as 1.5 billion years after the Big Bang. However, due to their faintness and red colour, it has proven extremely challenging to learn about these extreme quiescent galaxies, or to confirm whether any exist at earlier times. Here we report the spectroscopic confirmation of a massive quiescent galaxy, GS-9209, at redshift $z=4.658$, just 1.25 billion years after the Big Bang, using JWST NIRSpec. From these data we infer a stellar mass of $M_* = 3.8\pm0.2\times10^{10}\ M_\odot$, which formed over a $\simeq200$ Myr period before this galaxy quenched its star formation activity at $z=6.5^{+0.2}_{-0.5}$, when the Universe was $\simeq800$ million years old. Based on the presence of broad H$\alpha$ in the spectrum and a high narrow-line [NII]/H$\alpha$ ratio, we infer the presence of an accreting supermassive black hole, with a mass of $M_\bullet = 5\pm1\times10^{8}\ M_\odot$. This large black hole mass relative to the stellar mass suggests that active galactic nucleus (AGN) feedback may have been responsible for quenching this galaxy. GS-9209 is also extremely compact, with an effective radius, $r_e=215\pm20$ parsecs. This galaxy is both a likely descendent of the highest-redshift submillimetre galaxies and quasars, and a likely progenitor for the dense, ancient cores of the most massive local galaxies.


title: A Simple Algorithm For Scaling Up Kernel Methods

abstract: The recent discovery of the equivalence between infinitely wide neural networks (NNs) in the lazy training regime and Neural Tangent Kernels (NTKs) (Jacot et al., 2018) has revived interest in kernel methods. However, conventional wisdom suggests kernel methods are unsuitable for large samples due to their computational complexity and memory requirements. We introduce a novel random feature regression algorithm that allows us (when necessary) to scale to virtually infinite numbers of random features. We illustrate the performance of our method on the CIFAR-10 dataset.


title: Approximate Bilevel Difference Convex Programming for Bayesian Risk Markov Decision Processes

abstract: We consider infinite-horizon Markov Decision Processes where parameters, such as transition probabilities, are unknown and estimated from data. The popular distributionally robust approach to addressing the parameter uncertainty can sometimes be overly conservative. In this paper, we utilize the recently proposed formulation, Bayesian risk Markov Decision Process (BR-MDP), to address parameter (or epistemic) uncertainty in MDPs (Lin et al., 2022). To solve the infinite-horizon BR-MDP with a class of convex risk measures, we propose a computationally efficient approach of approximate bilevel difference convex programming (ABDCP). The optimization is performed offline and produces the optimal policy that is represented as a finite state controller with desirable performance guarantees. We also demonstrate the empirical performance of the infinite-horizon BR-MDP formulation and proposed algorithms.


title: Feature space exploration as an alternative for design space exploration beyond the parametric space

abstract: This paper compares the parametric design space with a feature space generated by the extraction of design features using deep learning (DL) as an alternative way for design space exploration. In this comparison, the parametric design space is constructed by creating a synthetic dataset of 15.000 elements using a parametric algorithm and reducing its dimensions for visualization. The feature space - reduced-dimensionality vector space of embedded data features - is constructed by training a DL model on the same dataset. We analyze and compare the extracted design features by reducing their dimension and visualizing the results. We demonstrate that parametric design space is narrow in how it describes the design solutions because it is based on the combination of individual parameters. In comparison, we observed that the feature design space can intuitively represent design solutions according to complex parameter relationships. Based on our results, we discuss the potential of translating the features learned by DL models to provide a mechanism for intuitive design exploration space and visualization of possible design solutions.


title: Are Labels Needed for Incremental Instance Learning?

abstract: In this paper, we learn to classify visual object instances, incrementally and via self-supervision (self-incremental). Our learner observes a single instance at a time, which is then discarded from the dataset. Incremental instance learning is challenging, since longer learning sessions exacerbate forgetfulness, and labeling instances is cumbersome. We overcome these challenges via three contributions: i. We propose VINIL, a self-incremental learner that can learn object instances sequentially, ii. We equip VINIL with self-supervision to by-pass the need for instance labelling, iii. We compare VINIL to label-supervised variants on two large-scale benchmarks, and show that VINIL significantly improves accuracy while reducing forgetfulness.


title: Parkinson gait modelling from an anomaly deep representation

abstract: Parkinson's Disease (PD) is associated with gait movement disorders, such as bradykinesia, stiffness, tremors and postural instability, caused by progressive dopamine deficiency. Today, some approaches have implemented learning representations to quantify kinematic patterns during locomotion, supporting clinical procedures such as diagnosis and treatment planning. These approaches assumes a large amount of stratified and labeled data to optimize discriminative representations. Nonetheless these considerations may restrict the approaches to be operable in real scenarios during clinical practice. This work introduces a self-supervised generative representation to learn gait-motion-related patterns, under the pretext of video reconstruction and an anomaly detection framework. This architecture is trained following a one-class weakly supervised learning to avoid inter-class variance and approach the multiple relationships that represent locomotion. The proposed approach was validated with 14 PD patients and 23 control subjects, and trained with the control population only, achieving an AUC of 95%, homocedasticity level of 70% and shapeness level of 70% in the classification task considering its generalization.


title: Efficiently predicting high resolution mass spectra with graph neural networks

abstract: Identifying a small molecule from its mass spectrum is the primary open problem in computational metabolomics. This is typically cast as information retrieval: an unknown spectrum is matched against spectra predicted computationally from a large database of chemical structures. However, current approaches to spectrum prediction model the output space in ways that force a tradeoff between capturing high resolution mass information and tractable learning. We resolve this tradeoff by casting spectrum prediction as a mapping from an input molecular graph to a probability distribution over molecular formulas. We discover that a large corpus of mass spectra can be closely approximated using a fixed vocabulary constituting only 2% of all observed formulas. This enables efficient spectrum prediction using an architecture similar to graph classification - GrAFF-MS - achieving significantly lower prediction error and orders-of-magnitude faster runtime than state-of-the-art methods.


title: Classical Algorithm for the Mean Value problem over Short-Time Hamiltonian Evolutions

abstract: Simulating physical systems has been an important application of classical and quantum computers. In this article we present an efficient classical algorithm for simulating time-dependent quantum mechanical Hamiltonians over constant periods of time. The algorithm presented here computes the mean value of an observable over the output state of such short-time Hamiltonian evolutions. In proving the performance of this algorithm we use Lieb-Robinson type bounds to limit the evolution of local operators within a lightcone. This allows us to divide the task of simulating a large quantum system into smaller systems that can be handled on normal classical computers.


title: Non-invasive and noise-robust light focusing using confocal wavefront shaping

abstract: Wavefront-shaping is a promising approach for imaging fluorescent targets deep inside scattering tissue despite strong aberrations. It enables focusing an incoming illumination into a single spot inside tissue, as well as correcting the outgoing light scattered from the tissue, by modulating the incoming and/or outgoing wavefronts. Previously, wavefront shaping modulations have been successively estimated using feedback from strong fluorescent beads, which have been manually added to a sample. However, ideally, such feedback should be provided by the fluorescent components of the tissue itself, whose emission is orders of magnitude weaker than the one provided by beads. When a low number of photons is spread over multiple sensor pixels, the image is highly susceptible to noise, and the feedback signal required for previous algorithms cannot be detected. In this work, we suggest a wavefront shaping approach that works with a confocal modulation of both the illumination and imaging arms. Since the aberrations are corrected in the optics before the detector, the low photon budget can be directed into a single sensor spot and detected with high SNR. We derive a score function for modulation evaluation from mathematical principles, and successfully use it to image EGFP labeled neurons, despite scattering through thick tissue.


title: RMSim: Controlled Respiratory Motion Simulation on Static Patient Scans

abstract: This work aims to generate realistic anatomical deformations from static patient scans. Specifically, we present a method to generate these deformations/augmentations via deep learning driven respiratory motion simulation that provides the ground truth for validating deformable image registration (DIR) algorithms and driving more accurate deep learning based DIR. We present a novel 3D Seq2Seq deep learning respiratory motion simulator (RMSim) that learns from 4D-CT images and predicts future breathing phases given a static CT image. The predicted respiratory patterns, represented by time-varying displacement vector fields (DVFs) at different breathing phases, are modulated through auxiliary inputs of 1D breathing traces so that a larger amplitude in the trace results in more significant predicted deformation. Stacked 3D-ConvLSTMs are used to capture the spatial-temporal respiration patterns. Training loss includes a smoothness loss in the DVF and mean-squared error between the predicted and ground truth phase images. A spatial transformer deforms the static CT with the predicted DVF to generate the predicted phase image. 10-phase 4D-CTs of 140 internal patients were used to train and test RMSim. The trained RMSim was then used to augment a public DIR challenge dataset for training VoxelMorph to show the effectiveness of RMSim-generated deformation augmentation. We validated our RMSim output with both private and public benchmark datasets (healthy and cancer patients). The proposed approach can be used for validating DIR algorithms as well as for patient-specific augmentations to improve deep learning DIR algorithms. The code, pretrained models, and augmented DIR validation datasets will be released at https://github.com/nadeemlab/SeqX2Y.


title: Bounds for Permutation Arrays under Kendall Tau Metric

abstract: Permutation arrays under the Kendall-$\tau$ metric have been considered for error-correcting codes. Given $n$ and $d\in [1..\binom{n}{2}]$, the task is to find a large permutation array of permutations on $n$ symbols with pairwise Kendall-$\tau$ distance at least $d$. Let $P(n,d)$ denote the maximum size of any permutation array of permutations on $n$ symbols with pairwise Kendall-$\tau$ distance $d$. New algorithms and several theorems are presented, giving improved lower bounds for $P(n,d)$. Also, $(n,m,d)$-arrays are defined, which are permutation arrays on n symbols with Kendall-$\tau$ distance d, with the restriction that symbols {1...(n-m)} appear in increasing order. Let $P(n,m,d)$ denote the maximum size of any $(n,m,d)$-array. For example, (n,m,d)-arrays are useful for recursively computing lower bounds for $P(n,d)$. Lower and upper bounds are given for $P(n.m,d)$.


title: An inductive model structure for strict $\infty$-categories

abstract: We construct a left semi-model category of ``marked strict $\infty$-categories'' for which the fibrant objects are those whose marked arrows satisfy natural closure properties and are weakly invertible. The canonical model structure on strict $\infty$-categories can be recovered as a left Bousfield localization of this model structure. We show that an appropriate extension of the Street nerve to the marked setting produces a Quillen adjunction between our model category and the Verity model structure for complicial sets, generalizing previous results by the second named author. Finally, we use this model structure to study, in the setting of strict $\infty$-categories, the idea that there are several non-equivalent notions of weak $(\infty,\infty)$-categories - depending on what tower of $(\infty,n)$-categories is used. We show that there ought to be at least three different notions of $(\infty,\infty)$-categories.


title: Anomalously high supercurrent density in a two-dimensional topological material

abstract: Ongoing advances in superconductors continue to revolutionize technology thanks to the increasingly versatile and robust availability of lossless supercurrent. In particular high supercurrent density can lead to more efficient and compact power transmission lines, high-field magnets, as well as high-performance nanoscale radiation detectors and superconducting spintronics. Here, we report the discovery of an unprecedentedly high superconducting critical current density (17 MA/cm2 at 0 T and 7 MA/cm2 at 8 T) in 1T'-WS2, exceeding those of all reported two-dimensional superconductors to date. 1T'-WS2 features a strongly anisotropic (both in- and out-of-plane) superconducting state that violates the Pauli paramagnetic limit signaling the presence of unconventional superconductivity. Spectroscopic imaging of the vortices further substantiates the anisotropic nature of the superconducting state. More intriguingly, the normal state of 1T'-WS2 carries topological properties. The band structure obtained via angle-resolved photoemission spectroscopy and first-principles calculations points to a Z2 topological invariant. The concomitance of topology and superconductivity in 1T'-WS2 establishes it as a topological superconductor candidate, which is promising for the development of quantum computing technology.

categories: cond-mat.supr-con cond-mat.mes-hall cond-mat.str-el physics.app-ph


title: Model-based Offline Reinforcement Learning with Local Misspecification

abstract: We present a model-based offline reinforcement learning policy performance lower bound that explicitly captures dynamics model misspecification and distribution mismatch and we propose an empirical algorithm for optimal offline policy selection. Theoretically, we prove a novel safe policy improvement theorem by establishing pessimism approximations to the value function. Our key insight is to jointly consider selecting over dynamics models and policies: as long as a dynamics model can accurately represent the dynamics of the state-action pairs visited by a given policy, it is possible to approximate the value of that particular policy. We analyze our lower bound in the LQR setting and also show competitive performance to previous lower bounds on policy selection across a set of D4RL tasks.


title: Resetting induced multimodality

abstract: Properties of stochastic systems are defined by the noise type and deterministic forces acting on the system. In out-of-equilibrium setups, e.g., for motions under action of L\'evy noises, the existence of the stationary state is not only determined by the potential but also by the noise. Potential wells need to be steeper than parabolic in order to assure existence of stationary states. The existence of stationary states, in sub-harmonic potential wells, can be restored by stochastic resetting, which is the protocol of starting over at random times. Herein we demonstrate that the combined action of L\'evy noise and Poissonian stochastic resetting can result in the phase transition between non-equilibrium stationary states of various multimodality in the overdamped system in super-harmonic potentials. Fine-tuned resetting rates can increase the modality of stationary states, while for high resetting rates the multimodality is destroyed as the stochastic resetting limits the spread of particles.


title: Recursive characterisations of random matrix ensembles and associated combinatorial objects

abstract: We give an overview of the recursive characterisations of random matrix ensembles that are currently at the forefront of random matrix theory by way of studying two classes of ensembles using two different types of recursive schemes: Established theory on Selberg correlation integrals is used to derive linear differential equations on the eigenvalue densities and resolvents of the classical matrix ensembles, which lead to $1$-point recursions, understood to be analogues of the Harer-Zagier recursion, for the expansion coefficients of the associated $1$-point cumulants, while loop equation analysis is used to recursively compute some leading order correlator expansion coefficients pertaining to certain products of random matrices that have recently come into interest due to their connections to Muttalib-Borodin ensembles and integrals of Harish-Chandra-Itzykson-Zuber type. We also show how the aforementioned differential equations can be used to characterise the large $N$ limiting statistics of the classical matrix ensembles' eigenvalue densities in the global and edge scaling regimes and moreover give a comprehensive review of how the Isserlis-Wick theorem implies ribbon graph interpretations for mixed moments and cumulants of some of the ensembles at hand. It is expected that this thesis will serve as a valuable resource for readers wanting a well-rounded introduction to classical matrix ensembles, matrix product ensembles, $1$-point recursions, loop equations, Selberg correlation integrals, and ribbon graphs.


title: Just Another Day on Twitter: A Complete 24 Hours of Twitter Data

abstract: At the end of October 2022, Elon Musk concluded his acquisition of Twitter. In the weeks and months before that, several questions were publicly discussed that were not only of interest to the platform's future buyers, but also of high relevance to the Computational Social Science research community. For example, how many active users does the platform have? What percentage of accounts on the site are bots? And, what are the dominating topics and sub-topical spheres on the platform? In a globally coordinated effort of 80 scholars to shed light on these questions, and to offer a dataset that will equip other researchers to do the same, we have collected all 375 million tweets published within a 24-hour time period starting on September 21, 2022. To the best of our knowledge, this is the first complete 24-hour Twitter dataset that is available for the research community. With it, the present work aims to accomplish two goals. First, we seek to answer the aforementioned questions and provide descriptive metrics about Twitter that can serve as references for other researchers. Second, we create a baseline dataset for future research that can be used to study the potential impact of the platform's ownership change.


title: Vortex sheet solutions for the Ginzburg-Landau system in cylinders: symmetry and global minimality

abstract: We consider the Ginzburg-Landau energy $E_\epsilon$ for $\mathbb{R}^M$-valued maps defined in a cylinder shape domain $B^N\times (0,1)^n$ satisfying a degree-one vortex boundary condition on $\partial B^N\times (0,1)^n$ in dimensions $M\geq N\geq 2$ and $n\geq 1$. The aim is to study the radial symmetry of global minimizers of this variational problem. We prove the following: if $N\geq 7$, then for every $\epsilon>0$, there exists a unique global minimizer which is given by the non-escaping radially symmetric vortex sheet solution $u_\epsilon(x,z)=(f_\epsilon(|x|) \frac{x}{|x|}, 0_{\mathbb{R}^{M-N}})$, $\forall x\in B^N$ that is invariant in $z\in (0,1)^n$. If $2\leq N \leq 6$ and $M\geq N+1$, the following dichotomy occurs between escaping and non-escaping solutions: there exists $\epsilon_N>0$ such that $\bullet$ if $\epsilon\in (0, \epsilon_N)$, then every global minimizer is an escaping radially symmetric vortex sheet solution of the form $R \tilde u_\epsilon$ where $\tilde{u}_{\epsilon}(x,z)=(\tilde{f}_{\epsilon}(|x|) \frac{x}{|x|}, 0_{\mathbb{R}^{M-N-1}}, g_{\epsilon}(|x|))$ is invariant in $z$-direction with $g_\epsilon>0$ in $(0,1)$ and $R\in O(M)$ is an orthogonal transformation keeping invariant the space $\mathbb{R}^N\times \{0_{\mathbb{R}^{M-N}}\}$; $\bullet$ if $\epsilon\geq \epsilon_N$, then the non-escaping radially symmetric vortex sheet solution $u_\epsilon(x,z)=(f_\epsilon(|x|) \frac{x}{|x|}, 0_{\mathbb{R}^{M-N}})$, $\forall x\in B^N, z\in (0,1)^n$ is the unique global minimizer; moreover, there are no bounded escaping solutions in this case. We also discuss the problem of vortex sheet $\mathbb{S}^{M-1}$-valued harmonic maps.


title: Semidefinite Relaxations for Robust Multiview Triangulation

abstract: We propose an approach based on convex relaxations for certifiably optimal robust multiview triangulation. To this end, we extend existing relaxation approaches to non-robust multiview triangulation by incorporating a truncated least squares cost function. We propose two formulations, one based on epipolar constraints and one based on fractional reprojection constraints. The first is lower dimensional and remains tight under moderate noise and outlier levels, while the second is higher dimensional and therefore slower but remains tight even under extreme noise and outlier levels. We demonstrate through extensive experiments that the proposed approaches allow us to compute provably optimal reconstructions even under significant noise and a large percentage of outliers.


title: Emotional Interaction Qualities: Vocabulary, Modalities, Actions, And Mapping

abstract: Have you ever typed particularly powerful on your keyboard, maybe even harsh, to write and send a message with some emphasis of your emotional state or message? Did it work? Probably not. It didn't affect how you typed or interacted with your mouse. But what if you had other, connected devices, with other modalities for inputs and outputs? Which would you have chosen, and how would you characterize your interactions with them? We researched with our multisensory and multimodal tool, the Loaded Dice, in co-design workshops the design space of IoT usage scenarios: what interaction qualities users want, characterized using an interaction vocabulary, and how they might map them to a selection of sensors and actuators. We discuss based on our experience some thoughts of such a mapping.


title: Strong Variability in AzV 493, an Extreme Oe-Type Star in the SMC

abstract: We present 18 years of OGLE photometry together with spectra obtained over 12 years, revealing that the early Oe star AzV 493 shows strong photometric (Delta I < 1.2 mag) and spectroscopic variability with a dominant, 14.6-year pattern and ~40-day oscillations. We estimate stellar parameters T_eff = 42000 K, log L/L_sun = 5.83 +/- 0.15, M/M_sun = 50 +/- 9, and vsini = 370 +/- 40 km/s. Direct spectroscopic evidence shows episodes of both gas ejection and infall. There is no X-ray detection, and it is likely a runaway star. AzV 493 may have an unseen companion on a highly eccentric (e > 0.93) orbit. We propose that close interaction at periastron excites ejection of the decretion disk, whose variable emission-line spectrum suggests separate inner and outer components, with an optically thick outer component obscuring both the stellar photosphere and the emission-line spectrum of the inner disk at early phases in the photometric cycle. It is plausible that AzV 493's mass and rotation have been enhanced by binary interaction followed by the core-collapse supernova explosion of the companion, which now could be either a black hole or neutron star. This system in the Small Magellanic Cloud can potentially shed light on OBe decretion disk formation and evolution, massive binary evolution, and compact binary progenitors.


title: The Field Structure of Free Photons

abstract: Using a quantum field theoretic description of the photon it is shown that, as intuitively expected but not before theoretically proven, the vector potential of a photon has a likely amplitude associated with a discrete frequency and therefore energy, and momentum. In particular, by finding the wave-functional for the vector potential, it is shown that the likely absolute amplitude spectrum has delta function at a given frequency. This analysis is extended to n-photon systems. It shows that such systems have a vector potential distribution whose most likely element has a strong sinusoidal component which has an amplitude corresponding to n-fold more energy than a single photon system. An analogous result for photons of different energy is also derived. Through the use of Parseval's theorem for stochastic systems, the calculations and associated analyses introduces a simple tool for exploring the nature of QFT Schrodinger wave-functional generally.


title: Learning Modulo Theories

abstract: Recent techniques that integrate \emph{solver layers} into Deep Neural Networks (DNNs) have shown promise in bridging a long-standing gap between inductive learning and symbolic reasoning techniques. In this paper we present a set of techniques for integrating \emph{Satisfiability Modulo Theories} (SMT) solvers into the forward and backward passes of a deep network layer, called SMTLayer. Using this approach, one can encode rich domain knowledge into the network in the form of mathematical formulas. In the forward pass, the solver uses symbols produced by prior layers, along with these formulas, to construct inferences; in the backward pass, the solver informs updates to the network, driving it towards representations that are compatible with the solver's theory. Notably, the solver need not be differentiable. We implement \layername as a Pytorch module, and our empirical results show that it leads to models that \emph{1)} require fewer training samples than conventional models, \emph{2)} that are robust to certain types of covariate shift, and \emph{3)} that ultimately learn representations that are consistent with symbolic knowledge, and thus naturally interpretable.


title: Synesthetic Dice: Sensors, Actuators, And Mappings

abstract: How bright can you cry? How loud does the sun shine? We developed a multisensory and multimodal tool, the Loaded Dice, for use in co-design workshops to research the design space of IoT usage scenarios. The Loaded Dice incorporate the principle of a technical synesthesia, being able to map any of the included sensors to any of the included actuators. With just a turn of one of the cubical devices it is possible to create a new combination. We discuss the core principles of the Loaded Dice, what sensors and actuators are included, how they relate to human senses, and how we realized a meaningful mapping between sensors and actuators. We further discuss where we see additional potential in the Loaded Dice to support synesthetic exploration - as Synesthetic Dice - so that you can eventually find out who cries brighter.


title: Densities for Elliptic Curves over Global Function Fields

abstract: Let $K$ be a global function field. Using Haar measures, we compute the densities of the Kodiara types and Tamagawa numbers of elliptic curves over a completion of $K$. Also, we prove results about the number of iterations of Tate's algorithm that are completed when the algorithm is used on elliptic curves over a completion of $K$.


title: Study of the Uncertainties of the Galactic Radio Background as a Calibration Source for Radio Arrays

abstract: The indirect detection of cosmic rays via the radio signal of extensive air showers is gaining a lot of ground. Many new arrays of radio antennas are under construction or in the phase of development. Calibrating these arrays is important for the reconstruction of observed events and for the comparability between observatories. Using reference antennas in calibration campaigns is not ideal because of uncertainties on their signal output strength that are large or difficult to assess. In a different approach the arrays can be calibrated against the Galactic radio emission as the dominant source of background. This so-called Galactic Calibration relies on predictions of the diffuse Galactic radio emission, for which models are publicly available. We present a comparison of these models in the frequency range from 10 to 408 MHz in order to estimate the systematic uncertainties on the strength of the Galactic background. We do this comparison on a global level as well as adapted for selected radio arrays and discuss implications for applying the Galactic calibration method. Furthermore we study the influence of the quiet Sun as an additional source of radio emission in the sky.


title: Primordial black holes in loop quantum cosmology: The effect on the threshold

abstract: Primordial black holes form in the early Universe and constitute one of the most viable candidates for dark matter. The study of their formation process requires the determination of a critical energy density perturbation threshold $\delta_\mathrm{c}$, which in general depends on the underlying gravity theory. Up to now, the majority of analytic and numerical techniques calculate $\delta_\mathrm{c}$ within the framework of general relativity. In this work, using simple physical arguments we estimate semi-analytically the PBH formation threshold within the framework of quantum gravity, working for concreteness within loop quantum cosmology (LQC). In particular, for low mass PBHs formed close to the quantum bounce, we find a reduction in the value of $\delta_\mathrm{c}$ up to $50\%$ compared to the general relativistic regime quantifying for the first time to the best of our knowledge how quantum effects can influence PBH formation within a quantum gravity framework. Finally, by varying the Barbero-Immirzi parameter $\gamma$ of loop quantum gravity (LQG) we show its effect on the value of $\delta_\mathrm{c}$ while using the observational/phenomenological signatures associated to ultra-light PBHs, namely the ones affected by LQG effects, we propose the PBH portal as a novel probe to constrain the potential quantum nature of gravity.


title: Secure synchronization of artificial neural networks used to correct errors in quantum cryptography

abstract: Quantum cryptography can provide a very high level of data security. However, a big challenge of this technique is errors in quantum channels. Therefore, error correction methods must be applied in real implementations. An example is error correction based on artificial neural networks. This paper considers the practical aspects of this recently proposed method and analyzes elements which influence security and efficiency. The synchronization process based on mutual learning processes is analyzed in detail. The results allowed us to determine the impact of various parameters. Additionally, the paper describes the recommended number of iterations for different structures of artificial neural networks and various error rates. All this aims to support users in choosing a suitable configuration of neural networks used to correct errors in a secure and efficient way.


title: Metaverse for Wireless Systems: Architecture, Advances, Standardization, and Open Challenges

abstract: The growing landscape of emerging wireless applications is a key driver toward the development of novel wireless system designs. Such a design can be based on the metaverse that uses a virtual model of the physical world systems along with other schemes/technologies (e.g., optimization theory, machine learning, and blockchain). A metaverse using a virtual model performs proactive intelligent analytics prior to a user request for efficient management of the wireless system resources. Additionally, a metaverse will enable self-sustainability to operate wireless systems with the least possible intervention from network operators. Although the metaverse can offer many benefits, it faces some challenges as well. Therefore, in this tutorial, we discuss the role of a metaverse in enabling wireless applications. We present an overview, key enablers, design aspects (i.e., metaverse for wireless and wireless for metaverse), and a novel high-level architecture of metaverse-based wireless systems. We discuss metaverse management, reliability, and security of the metaverse-based system. Furthermore, we discuss recent advances and standardization of metaverse-enabled wireless system. Finally, we outline open challenges and present possible solutions.


title: Communication-Efficient Collaborative Regret Minimization in Multi-Armed Bandits

abstract: In this paper, we study the collaborative learning model, which concerns the tradeoff between parallelism and communication overhead in multi-agent multi-armed bandits. For regret minimization in multi-armed bandits, we present the first set of tradeoffs between the number of rounds of communication among the agents and the regret of the collaborative learning process.


title: Limitless stability for Graph Convolutional Networks

abstract: This work establishes rigorous, novel and widely applicable stability guarantees and transferability bounds for graph convolutional networks -- without reference to any underlying limit object or statistical distribution. Crucially, utilized graph-shift operators (GSOs) are not necessarily assumed to be normal, allowing for the treatment of networks on both undirected- and for the first time also directed graphs. Stability to node-level perturbations is related to an 'adequate (spectral) covering' property of the filters in each layer. Stability to edge-level perturbations is related to Lipschitz constants and newly introduced semi-norms of filters. Results on stability to topological perturbations are obtained through recently developed mathematical-physics based tools. As an important and novel example, it is showcased that graph convolutional networks are stable under graph-coarse-graining procedures (replacing strongly-connected sub-graphs by single nodes) precisely if the GSO is the graph Laplacian and filters are regular at infinity. These new theoretical results are supported by corresponding numerical investigations.


title: Dynamical vs Supernova Acceleration of OB Stars in the Small Magellanic Cloud

abstract: We use the RIOTS4 sample of SMC field OB stars to determine the origin of massive runaways in this low-metallicity galaxy using Gaia proper motions, together with stellar masses obtained from RIOTS4 data. These data allow us to estimate the relative contributions of stars accelerated by the dynamical ejection vs binary supernova mechanisms, since dynamical ejection favors faster, more massive runaways, while SN ejection favors the opposite trend. In addition, we use the frequencies of classical OBe stars, high-mass X-ray binaries, and non-compact binaries to discriminate between the mechanisms. Our results show that the dynamical mechanism dominates by a factor of 2 - 3. This also implies a significant contribution from two-step acceleration that occurs when dynamically ejected binaries are followed by SN kicks. We update our published quantitative results from Gaia DR2 proper motions with new data from DR3.


title: 3DShape2VecSet: A 3D Shape Representation for Neural Fields and Generative Diffusion Models

abstract: We introduce 3DShape2VecSet, a novel shape representation for neural fields designed for generative diffusion models. Our shape representation can encode 3D shapes given as surface models or point clouds, and represents them as neural fields. The concept of neural fields has previously been combined with a global latent vector, a regular grid of latent vectors, or an irregular grid of latent vectors. Our new representation encodes neural fields on top of a set of vectors. We draw from multiple concepts, such as the radial basis function representation and the cross attention and self-attention function, to design a learnable representation that is especially suitable for processing with transformers. Our results show improved performance in 3D shape encoding and 3D shape generative modeling tasks. We demonstrate a wide variety of generative applications: unconditioned generation, category-conditioned generation, text-conditioned generation, point-cloud completion, and image-conditioned generation.


title: On granularity of prosodic representations in expressive text-to-speech

abstract: In expressive speech synthesis it is widely adopted to use latent prosody representations to deal with variability of the data during training. Same text may correspond to various acoustic realizations, which is known as a one-to-many mapping problem in text-to-speech. Utterance, word, or phoneme-level representations are extracted from target signal in an auto-encoding setup, to complement phonetic input and simplify that mapping. This paper compares prosodic embeddings at different levels of granularity and examines their prediction from text. We show that utterance-level embeddings have insufficient capacity and phoneme-level tend to introduce instabilities when predicted from text. Word-level representations impose balance between capacity and predictability. As a result, we close the gap in naturalness by 90% between synthetic speech and recordings on LibriTTS dataset, without sacrificing intelligibility.


title: Personalised Federated Learning On Heterogeneous Feature Spaces

abstract: Most personalised federated learning (FL) approaches assume that raw data of all clients are defined in a common subspace i.e. all clients store their data according to the same schema. For real-world applications, this assumption is restrictive as clients, having their own systems to collect and then store data, may use heterogeneous data representations. We aim at filling this gap. To this end, we propose a general framework coined FLIC that maps client's data onto a common feature space via local embedding functions. The common feature space is learnt in a federated manner using Wasserstein barycenters while the local embedding functions are trained on each client via distribution alignment. We integrate this distribution alignement mechanism into a federated learning approach and provide the algorithmics of FLIC. We compare its performances against FL benchmarks involving heterogeneous input features spaces. In addition, we provide theoretical insights supporting the relevance of our methodology.


title: Far-Infrared Signatures for a Two-Steps Pressure-Driven Metallization in Transition Metal Dichalcogenides

abstract: We present a high-pressure investigation of the semiconductor-to-metal transition in MoS2 and WS2 carried out by synchrotron-based far-infrared spectroscopy, to reconcile the controversial estimates of the metallization pressure found in the literature and gain new insight into the mechanisms ruling this electronic transition. Two spectral descriptors are found indicative of the onset of metallicity and of the origin of the free carriers in the metallic state: the absorbance spectral weight, whose abrupt increase defines the metallization pressure threshold, and the asymmetric lineshape of the E1u peak, whose pressure evolution, interpreted within the Fano model, suggests the electrons in the metallic state originate from n-type doping levels. Combining our results with those reported in the literature, we hypothesize a two-step mechanism is at work in the metallization process, in which the pressure-induced hybridization between doping and conduction band states drives an early metallic behaviour, while the band-gap closes at higher pressures.


title: A realization of poset associahedra

abstract: Given any connected poset $P$, we give a simple realization of Galashin's poset associahedron $\mathscr{A}(P)$ as a convex polytope in $\mathbb{R}^P.$ The realization is inspired by the description of $\mathscr{A}(P)$ as a compactification of the configuration space of order-preserving maps $P \to \mathbb{R}.$ In addition, we give an analogous realization for Galashin's affine poset cyclohedra.


title: SMILE: Robust Network Localization via Sparse and Low-Rank Matrix Decomposition

abstract: Motivated by collaborative localization in robotic sensor networks, we consider the problem of large-scale network localization where location estimates are derived from inter-node radio signals. Well-established methods for network localization commonly assume that all radio links are line-of-sight and subject to Gaussian noise. However, the presence of obstacles which cause non-line-of-sight attenuation present distinct challenges. To enable robust network localization, we present Sparse Matrix Inference and Linear Embedding (SMILE), a novel approach which draws on both the well-known Locally Linear Embedding (LLE) algorithm and recent advances in sparse plus low-rank matrix decomposition. We demonstrate that our approach is robust to noisy signal propagation, severe attenuation due to non-line-of-sight, and missing pairwise measurements. Our experiments include simulated large-scale networks, an 11-node sensor network, and an 18-node network of mobile robots and static anchor radios in a GPS-denied limestone mine. Our findings indicate that SMILE outperforms classical multidimensional scaling (MDS) which ignores the effect of non-line of sight (NLOS), as well as outperforming state-of-the-art robust network localization algorithms that do account for NLOS attenuation including a graph convolutional network-based approach. We demonstrate that this improved accuracy is not at the cost of complexity, as SMILE sees reduced computation time for very large networks which is important for position estimation updates in a dynamic setting, e.g for mobile robots.


title: Switchable anomalous Hall effects in polar-stacked 2D antiferromagnet MnBi2Te4

abstract: Van der Waals (vdW) assembly allows controlling symmetry of two-dimensional (2D) materials that determines their physical properties. Especially interesting is the recently demonstrated breaking inversion symmetry by polar layer stacking to realize novel electronic, magnetic, and transport properties of 2D vdW materials switchable by induced electric polarization. Here, based on symmetry analyses and density-functional calculations, we explore the emergence of the anomalous Hall effect (AHE) in antiferromagnetic MnBi2Te4 films assembled by polar layer stacking. We demonstrate that breaking PT symmetry in an MnBi2Te4 bilayer makes this 2D material magnetoelectric and produces a spontaneous AHE switchable by electric polarization. We find that reversable polarization at one of the interfaces in a three-layer MnBi2Te4 film drives a metal-insulator transition, as well as switching between an AHE and quantum AHE (QAHE). Finally, we predict that engineering an interlayer polarization in a three-layer MnBi2Te4 film allows converting MnBi2Te4 from a trivial insulator to a Chern insulator. Overall, our work emphasizes the emergence of quantum-transport phenomena in 2D vdW antiferromagnets by polar layer stacking, which do not exist in this material in the bulk or bulk-like thin-film forms.


title: Transient Foreshock Structures Upstream of Mars: Implications of the Small Martian Bow Shock

abstract: We characterize the nature of magnetic structures in the foreshock region of Mars associated with discontinuities in the solar wind. The structures form at the upstream edge of moving foreshocks caused by slow rotations in the interplanetary magnetic field (IMF). The solar wind plasma density and the IMF strength noticeably decrease inside the structures' core, and a compressional shock layer is present at their sunward side, making them consistent with foreshock bubbles (FBs). Ion populations responsible for these structures include backstreaming ions that only appear within the moving foreshock, and accelerated reflected ions from the quasi-perpendicular bow shock. Both ion populations accumulate near the upstream edge of the moving foreshock which facilitates FB formation. Reflected ions with hybrid trajectories that straddle between the quasi-perpendicular and quasi-parallel bow shocks during slow IMF rotations contribute to formation of foreshock transients.


title: Coupled continuum modeling of size-segregation driven by shear-strain-rate gradients and flow in dense, bidisperse granular media

abstract: Dense granular systems that consist of particles of disparate sizes segregate based on size during flow, resulting in complex, coupled segregation and flow patterns. The ability to predict how granular mixtures segregate is important in the design of industrial processes and the understanding of geophysical phenomena. The two primary drivers of size-segregation are pressure gradients and shear-strain-rate gradients. In this work, we isolate size-segregation driven by shear-strain-rate gradients by studying two dense granular flow geometries with constant pressure fields: gravity-driven flow down a long vertical chute with rough parallel walls and annular shear flow with rough inner and outer walls. We perform discrete element method (DEM) simulations of dense flow of bidisperse granular systems in both flow geometries, while varying system parameters, such as the flow rate, flow configuration size, fraction of large/small grains, and grain-size ratio, and use DEM data to inform continuum constitutive equations for the relative flux of large and small particles. When the resulting continuum model for the dynamics of size-segregation is coupled with the nonlocal granular fluidity model (a nonlocal continuum model for dense granular flow rheology), we show that both flow fields and segregation dynamics may be simultaneously captured using this coupled, continuum system of equations.


title: Boundary Aware U-Net for Glacier Segmentation

abstract: Large-scale study of glaciers improves our understanding of global glacier change and is imperative for monitoring the ecological environment, preventing disasters, and studying the effects of global climate change. Glaciers in the Hindu Kush Himalaya (HKH) are particularly interesting as the HKH is one of the world's most sensitive regions for climate change. In this work, we: (1) propose a modified version of the U-Net for large-scale, spatially non-overlapping, clean glacial ice, and debris-covered glacial ice segmentation; (2) introduce a novel self-learning boundary-aware loss to improve debris-covered glacial ice segmentation performance; and (3) propose a feature-wise saliency score to understand the contribution of each feature in the multispectral Landsat 7 imagery for glacier mapping. Our results show that the debris-covered glacial ice segmentation model trained using self-learning boundary-aware loss outperformed the model trained using dice loss. Furthermore, we conclude that red, shortwave infrared, and near-infrared bands have the highest contribution toward debris-covered glacial ice segmentation from Landsat 7 images.


title: First Observations of the Brown Dwarf HD 19467 B with JWST

abstract: We observed HD 19467 B with JWST's NIRCam in six filters spanning 2.5-4.6 $\mu m$ with the Long Wavelength Bar coronagraph. The brown dwarf HD 19467 B was initially identified through a long-period trend in the radial velocity of G3V star HD 19467. HD 19467 B was subsequently detected via coronagraphic imaging and spectroscopy, and characterized as a late-T type brown dwarf with approximate temperature $\sim1000$K. We observed HD 19467 B as a part of the NIRCam GTO science program, demonstrating the first use of the NIRCam Long Wavelength Bar coronagraphic mask. The object was detected in all 6 filters (contrast levels of $2\times10^{-4}$ to $2\times10^{-5}$) at a separation of 1.6 arcsec using Angular Differential Imaging (ADI) and Synthetic Reference Differential Imaging (SynRDI). Due to a guidestar failure during acquisition of a pre-selected reference star, no reference star data was available for post-processing. However, RDI was successfully applied using synthetic Point Spread Functions (PSFs) developed from contemporaneous maps of the telescope's optical configuration. Additional radial velocity data (from Keck/HIRES) are used to constrain the orbit of HD 19467 B. Photometric data from TESS are used to constrain the properties of the host star, particularly its age. NIRCam photometry, spectra and photometry from literature, and improved stellar parameters are used in conjunction with recent spectral and evolutionary substellar models to derive physical properties for HD 19467 B. Using an age of 9.4$\pm$0.9 Gyr inferred from spectroscopy, Gaia astrometry, and TESS asteroseismology, we obtain a model-derived mass of 62$\pm 1M_{J}$, which is consistent within 2-$\sigma$ with the dynamically derived mass of 81$^{+14}_{-12}M_{J}$.


title: Graph Scattering beyond Wavelet Shackles

abstract: This work develops a flexible and mathematically sound framework for the design and analysis of graph scattering networks with variable branching ratios and generic functional calculus filters. Spectrally-agnostic stability guarantees for node- and graph-level perturbations are derived; the vertex-set non-preserving case is treated by utilizing recently developed mathematical-physics based tools. Energy propagation through the network layers is investigated and related to truncation stability. New methods of graph-level feature aggregation are introduced and stability of the resulting composite scattering architectures is established. Finally, scattering transforms are extended to edge- and higher order tensorial input. Theoretical results are complemented by numerical investigations: Suitably chosen cattering networks conforming to the developed theory perform better than traditional graph-wavelet based scattering approaches in social network graph classification tasks and significantly outperform other graph-based learning approaches to regression of quantum-chemical energies on QM7.


title: Attacking Important Pixels for Anchor-free Detectors

abstract: Deep neural networks have been demonstrated to be vulnerable to adversarial attacks: subtle perturbation can completely change the prediction result. Existing adversarial attacks on object detection focus on attacking anchor-based detectors, which may not work well for anchor-free detectors. In this paper, we propose the first adversarial attack dedicated to anchor-free detectors. It is a category-wise attack that attacks important pixels of all instances of a category simultaneously. Our attack manifests in two forms, sparse category-wise attack (SCA) and dense category-wise attack (DCA), that minimize the $L_0$ and $L_\infty$ norm-based perturbations, respectively. For DCA, we present three variants, DCA-G, DCA-L, and DCA-S, that select a global region, a local region, and a semantic region, respectively, to attack. Our experiments on large-scale benchmark datasets including PascalVOC, MS-COCO, and MS-COCO Keypoints indicate that our proposed methods achieve state-of-the-art attack performance and transferability on both object detection and human pose estimation tasks.


title: Superconductivity of anomalous pseudospin

abstract: Spin-orbit coupling driven by broken inversion symmetry ($I$) is known to lead to unusual magnetic response of superconductors, including extremely large critical fields for spin-singlet superconductors. This unusual response is also known to appear in materials that have $I$, provided there is local $I$-breaking: fermions participating in superconductivity reside on crystal sites that lack $I$. Here we show that this unusual response exists even when the crystal sites preserve $I$. Indeed, we argue that the symmetry of Kramers degenerate fermionic pseudospin is more relevant than the local crystal site symmetry. We examine and classify non-symmorphic materials with momentum space spin-textures that exhibit an anomalous pseudospin with different symmetry properties than usual spin-1/2. We find that this anomalous pseudospin does not depend on the existence of local $I$ breaking crystal sites and it optimizes the unusual magnetic response traditionally associated with locally noncentrosymmetric superconductors, dramatically extending the range of relevant materials. We further show this anomalous pseudospin leads to fully gapped `nodal' superconductors and provides additional insight into the breakdown of Blount's theorem for pseudospin triplet superconductors. We apply our results to UPt$_3$, BiS$_2$-based superconductors, Fe-based superconductors, and paramagnetic UCoGe.


title: Neural-Symbolic Inference for Robust Autoregressive Graph Parsing via Compositional Uncertainty Quantification

abstract: Pre-trained seq2seq models excel at graph semantic parsing with rich annotated data, but generalize worse to out-of-distribution (OOD) and long-tail examples. In comparison, symbolic parsers under-perform on population-level metrics, but exhibit unique strength in OOD and tail generalization. In this work, we study compositionality-aware approach to neural-symbolic inference informed by model confidence, performing fine-grained neural-symbolic reasoning at subgraph level (i.e., nodes and edges) and precisely targeting subgraph components with high uncertainty in the neural parser. As a result, the method combines the distinct strength of the neural and symbolic approaches in capturing different aspects of the graph prediction, leading to well-rounded generalization performance both across domains and in the tail. We empirically investigate the approach in the English Resource Grammar (ERG) parsing problem on a diverse suite of standard in-domain and seven OOD corpora. Our approach leads to 35.26% and 35.60% error reduction in aggregated Smatch score over neural and symbolic approaches respectively, and 14% absolute accuracy gain in key tail linguistic categories over the neural model, outperforming prior state-of-art methods that do not account for compositionality or uncertainty.


title: Evidencing the squeezed dark nuclear spin state in lead halide perovskites

abstract: Coherent many-body states are highly promising for robust and scalable quantum information processing. While far-reaching theoretical predictions have been made for various implementations, direct experimental evidence of their appealing properties can be challenging. Here, we demonstrate coherent optical manipulation of the nuclear spin ensemble in the lead halide perovskite semiconductor FAPbBr$_3$ (FA=formamidinium), targeting a long-postulated collective dark state that is insensitive to optical pumping. Via optical orientation of localized hole spins we drive the nuclear many-body system into an entangled state, requiring a weak magnetic field of only a few Millitesla strength at cryogenic temperatures. During its fast build-up, the nuclear polarization along the optical axis remains small, while the transverse nuclear spin fluctuations are strongly reduced, corresponding to spin squeezing as evidenced by a strong violation of the generalized nuclear squeezing-inequality with $\xi_s < 0.3$. The dark state evidenced in this process corresponds to an approximately 750-body entanglement between the nuclei. Dark nuclear spin states can be exploited to store quantum information benefiting from their long-lived many-body coherence and to perform quantum measurements with a precision beyond the standard limit.


title: Learning to Generate All Feasible Actions

abstract: Several machine learning (ML) applications are characterized by searching for an optimal solution to a complex task. The search space for this optimal solution is often very large, so large in fact that this optimal solution is often not computable. Part of the problem is that many candidate solutions found via ML are actually infeasible and have to be discarded. Restricting the search space to only the feasible solution candidates simplifies finding an optimal solution for the tasks. Further, the set of feasible solutions could be re-used in multiple problems characterized by different tasks. In particular, we observe that complex tasks can be decomposed into subtasks and corresponding skills. We propose to learn a reusable and transferable skill by training an actor to generate all feasible actions. The trained actor can then propose feasible actions, among which an optimal one can be chosen according to a specific task. The actor is trained by interpreting the feasibility of each action as a target distribution. The training procedure minimizes a divergence of the actor's output distribution to this target. We derive the general optimization target for arbitrary f-divergences using a combination of kernel density estimates, resampling, and importance sampling. We further utilize an auxiliary critic to reduce the interactions with the environment. A preliminary comparison to related strategies shows that our approach learns to visit all the modes in the feasible action space, demonstrating the framework's potential for learning skills that can be used in various downstream tasks.


title: How poor is the stimulus? Evaluating hierarchical generalization in neural networks trained on child-directed speech

abstract: When acquiring syntax, children consistently choose hierarchical rules over competing non-hierarchical possibilities. Is this preference due to a learning bias for hierarchical structure, or due to more general biases that interact with hierarchical cues in children's linguistic input? We explore these possibilities by training LSTMs and Transformers - two types of neural networks without a hierarchical bias - on data similar in quantity and content to children's linguistic input: text from the CHILDES corpus. We then evaluate what these models have learned about English yes/no questions, a phenomenon for which hierarchical structure is crucial. We find that, though they perform well at capturing the surface statistics of child-directed speech (as measured by perplexity), both model types generalize in a way more consistent with an incorrect linear rule than the correct hierarchical rule. These results suggest that human-like generalization from text alone requires stronger biases than the general sequence-processing biases of standard neural network architectures.


title: Nik Defense: An Artificial Intelligence Based Defense Mechanism against Selfish Mining in Bitcoin

abstract: The Bitcoin cryptocurrency has received much attention recently. In the network of Bitcoin, transactions are recorded in a ledger. In this network, the process of recording transactions depends on some nodes called miners that execute a protocol known as mining protocol. One of the significant aspects of mining protocol is incentive compatibility. However, literature has shown that Bitcoin mining's protocol is not incentive-compatible. Some nodes with high computational power can obtain more revenue than their fair share by adopting a type of attack called the selfish mining attack. In this paper, we propose an artificial intelligence-based defense against selfish mining attacks by applying the theory of learning automata. The proposed defense mechanism ignores private blocks by assigning weight based on block discovery time and changes current Bitcoin's fork resolving policy by evaluating branches' height difference in a self-adaptive manner utilizing learning automata. To the best of our knowledge, the proposed protocol is the literature's first learning-based defense mechanism. Simulation results have shown the superiority of the proposed mechanism against tie-breaking mechanism, which is a well-known defense. The simulation results have shown that the suggested defense mechanism increases the profit threshold up to 40\% and decreases the revenue of selfish attackers.


title: The formation of supermassive black holes from Population III.1 seeds. II. Evolution to the local universe

abstract: We present predictions for cosmic evolution of populations of supermassive black holes (SMBHs) forming from Population III.1 seeds, i.e., early, metal-free dark matter minihalos forming far from other sources, parameterized by isolation distance, $d_{\rm{iso}}$. Extending previous work that explored this scenario to $z=10$, we follow evolution of a $(60\:{\rm{Mpc}})^3$ volume to $z=0$. We focus on evolution of SMBH comoving number densities, halo occupation fractions, angular clustering and 3D clustering, exploring a range of $d_{\rm{iso}}$ constrained by observed local number densities of SMBHs. We also compute synthetic projected observational fields, in particular a case comparable to the Hubble Ultra Deep Field. We compare Pop III.1 seeding to a simple halo mass threshold model, commonly adopted in cosmological simulations of galaxy formation. Major predictions of the Pop III.1 model include that all SMBHs form by $z\sim25$, after which their comoving number densities are near-constant, with low merger rates. Occupation fractions evolve to concentrate SMBHs in the most massive halos by $z=0$, but with rare cases in halos down to $\sim10^8\:M_\odot$. The $d_{\rm{iso}}$ scale at epoch of formation, e.g., $100\:$kpc-proper at $z\sim30$, i.e., $\sim3\:$Mpc-comoving, is imprinted in the SMBH two-point angular correlation function, remaining discernible as a low-amplitude feature to $z\sim1$. The SMBH 3D two-point correlation function at $z=0$ also shows lower amplitude compared to equivalently massive halos. We discuss prospects for testing these predictions with observational surveys of SMBH populations.


title: Steinberg quotients, Weyl Characters, and Kazhdan-Lusztig Polynomials

abstract: Let $G$ be a reductive group over a field of prime characteristic. An indecomposable tilting module for $G$ whose highest weight lies above the Steinberg weight has a character that is divisible by the Steinberg character. The resulting "Steinberg quotient" carries important information about $G$-modules, and in previous work we studied patterns in the weight multiplicities of these characters. In this paper we broaden our scope to include quantum Steinberg quotients, and show how the multiplicities in these characters relate to algebraic Steinberg quotients, Weyl characters, and evaluations of Kazhdan-Lusztig polynomials. We give an explicit algorithm for computing minimal characters that possess a key attribute of Steinberg quotients. We provide computations which show that these minimal characters are not always equal to quantum Steinberg quotients, but are close in several nontrivial cases.


title: Hydrogen atom/molecule adsorption on 2D metallic porphyrin: A first-principles study

abstract: Hydrogen is a promising element for applications in new energy sources like fuel cells. One key issue for such applications is storing hydrogen. And, to improve storage capacity, understanding the interaction mechanism between hydrogen and possible storage materials is critical. This work uses DFT simulations to comprehensively investigate the adsorption mechanism of H/H$_2$ on the 2D metallic porphyrins with one transition metal in its center. Our results suggest that the mechanism for adsorption of H (H$_2$) is chemisorption (physisorption). The maximum adsorption energy for atomic hydrogen was $-3.7$ eV for 2D porphyrins embedded with vanadium or chromium atoms. Our results also revealed charge transfer of up $-0.43$ e to chemisorbed H atoms. In contrast, the maximum adsorption energy calculated for molecular hydrogen was $-122.5$ meV for 2D porphyrins embedded with scandium atoms. Furthermore, charge transfer was minimal for physisorption. Finally, we also determined that uniaxial strain has a minimal effect on the adsorption properties of 2D metallic porphyrins.


title: Cross-domain recommendation via user interest alignment

abstract: Cross-domain recommendation aims to leverage knowledge from multiple domains to alleviate the data sparsity and cold-start problems in traditional recommender systems. One popular paradigm is to employ overlapping user representations to establish domain connections, thereby improving recommendation performance in all scenarios. Nevertheless, the general practice of this approach is to train user embeddings in each domain separately and then aggregate them in a plain manner, often ignoring potential cross-domain similarities between users and items. Furthermore, considering that their training objective is recommendation task-oriented without specific regularizations, the optimized embeddings disregard the interest alignment among user's views, and even violate the user's original interest distribution. To address these challenges, we propose a novel cross-domain recommendation framework, namely COAST, to improve recommendation performance on dual domains by perceiving the cross-domain similarity between entities and aligning user interests. Specifically, we first construct a unified cross-domain heterogeneous graph and redefine the message passing mechanism of graph convolutional networks to capture high-order similarity of users and items across domains. Targeted at user interest alignment, we develop deep insights from two more fine-grained perspectives of user-user and user-item interest invariance across domains by virtue of affluent unsupervised and semantic signals. We conduct intensive experiments on multiple tasks, constructed from two large recommendation data sets. Extensive results show COAST consistently and significantly outperforms state-of-the-art cross-domain recommendation algorithms as well as classic single-domain recommendation methods.


title: Multi-limb Split Learning for Tumor Classification on Vertically Distributed Data

abstract: Brain tumors are one of the life-threatening forms of cancer. Previous studies have classified brain tumors using deep neural networks. In this paper, we perform the later task using a collaborative deep learning technique, more specifically split learning. Split learning allows collaborative learning via neural networks splitting into two (or more) parts, a client-side network and a server-side network. The client-side is trained to a certain layer called the cut layer. Then, the rest of the training is resumed on the server-side network. Vertical distribution, a method for distributing data among organizations, was implemented where several hospitals hold different attributes of information for the same set of patients. To the best of our knowledge this paper will be the first paper to implement both split learning and vertical distribution for brain tumor classification. Using both techniques, we were able to achieve train and test accuracy greater than 90\% and 70\%, respectively.


title: Field-free switching of perpendicular magnetization in an ultrathin epitaxial magnetic insulator

abstract: For energy efficient and fast magnetic memories, switching of perpendicular magnetization by the spin-orbit torque (SOT) appears as a very promising solution, even more using magnetic insulators that suppress electrical shunting. This SOT switching generally requires the assistance of an in-plane magnetic field to break the symmetry. Here, we present experiments demonstrating the field-free SOT switching of perpendicularly magnetized layers of the thulium iron garnet (Tm_{3}Fe_{5}O_{12}) magnetic insulator. The polarity of the switching loops, clockwise (CW) or counter-clockwise (CCW), is determined by the direction of the initial current pulses, in contrast with field-assisted switchings in which this polarity is controlled by the direction of the field. After an independent determination of Dzyaloshinskii-Moriya interaction (DMI), we relate the field free switching to the interplay of SOT and DMI and the polarity of the loops to the imprint of a N\'eel domain wall induced by the first pulse, in agreement with Kerr imaging. Our observation and interpretation of field-free electrical switching of a magnetic insulator is an important milestone for future low power spintronic devices.


title: Photons from dark photon solitons via parametric resonance

abstract: Wave-like dark matter made of spin-1 particles (dark photons) is expected to form ground state clumps called "vector solitons", which can have different polarizations. In this work, we consider the interaction of dark photons with photons, expressed as dimension-6 operators, and study the electromagnetic radiation that arises from an isolated vector soliton due to parametric resonant amplification of the ambient electromagnetic field. We characterize the directional dependence and polarization of the outgoing radiation, which depends on the operator as well as the polarization state of the underlying vector soliton. We discuss the implications of this radiation for the stability of solitons and as a possible channel for detecting mergers of vector solitons through astrophysical observations.


title: Multi-channel Medium Access Control Protocols for Wireless Networks within Computing Packages

abstract: Wireless communications at the chip scale emerge as a interesting complement to traditional wire-based approaches thanks to their low latency, inherent broadcast nature, and capacity to bypass pin constraints. However, as current trends push towards massive and bandwidth-hungry processor architectures, there is a need for wireless chip-scale networks that exploit and share as many channels as possible. In this context, this work addresses the issue of channel sharing by exploring the design space of multi-channel Medium Access Control (MAC) protocols for chip-scale networks. Distinct channel assignment strategies for both random access and token passing are presented and evaluated under realistic traffic patterns. It is shown that, even with the improvements enabled by the multiple channels, both protocols maintain their intrinsic advantages and disadvantages.


title: Fast Bayesian inference for spatial mean-parameterized Conway-Maxwell-Poisson models

abstract: Count data with complex features arise in many disciplines, including ecology, agriculture, criminology, medicine, and public health. Zero inflation, spatial dependence, and non-equidispersion are common features in count data. There are two classes of models that allow for these features -- he mode-parameterized Conway--Maxwell--Poisson (COMP) distribution and the generalized Poisson model. However both require the use of either constraints on the parameter space or a parameterization that leads to challenges in interpretability. We propose a spatial mean-parameterized COMP model that retains the flexibility of these models while resolving the above issues. We use a Bayesian spatial filtering approach in order to efficiently handle high-dimensional spatial data and we use reversible-jump MCMC to automatically choose the basis vectors for spatial filtering. The COMP distribution poses two additional computational challenges -- an intractable normalizing function in the likelihood and no closed-form expression for the mean. We propose a fast computational approach that addresses these challenges by, respectively, introducing an efficient auxiliary variable algorithm and pre-computing key approximations for fast likelihood evaluation. We illustrate the application of our methodology to simulated and real datasets, including Texas HPV-cancer data and US vaccine refusal data.


title: Proof of a conjecture of Krawchuk and Rampersad

abstract: We prove a 2018 conjecture of Krawchuk and Rampersad on the extremal behavior of $c(n)$, where $c(n)$ counts the number of length-$n$ factors of the Thue-Morse word $\mathbf{t}$, up to cyclic rotation.


title: Surface heating steers planetary-scale ocean circulation

abstract: Gyres are central features of large-scale ocean circulation and are involved in transporting tracers such as heat, nutrients, and carbon-dioxide within and across ocean basins. Traditionally, the gyre circulation is thought to be driven by surface winds and quantified via Sverdrup balance, but it has been proposed that surface buoyancy fluxes may also contribute to gyre forcing. Through a series of eddy-permitting global ocean model simulations with perturbed surface forcing, the relative contribution of wind stress and surface heat flux forcing to the large-scale ocean circulation is investigated, focusing on the subtropical gyres. In addition to gyre strength being linearly proportional to wind stress, it is shown that the gyre circulation is strongly impacted by variations in the surface heat flux (specifically, its meridional gradient) through a rearrangement of the ocean's buoyancy structure. On shorter timescales ($\sim$ decade), the gyre circulation anomalies are proportional to the magnitude of the surface heat flux gradient perturbation, with up to $\sim 0.15\,\mathrm{Sv}$ anomaly induced per $\mathrm{W}\,\mathrm{m}^{-2}$ change in the surface heat flux. On timescales longer than a decade, the gyre response to surface buoyancy flux gradient perturbations becomes non-linear as ocean circulation anomalies feed back onto the buoyancy structure induced by the surface buoyancy fluxes. These interactions complicate the development of a buoyancy-driven theory for the gyres to complement the Sverdrup relation. The flux-forced simulations underscore the importance of surface buoyancy forcing in steering the large-scale ocean circulation.


title: The effect of primary school education on preventive behaviours during COVID-19 in Japan

abstract: Education plays a critical role on promoting preventive behaviours against the spread of pandemics. In Japan, hand-washing education in primary schools was positively correlated with preventive behaviours against COVID-19 transmission for adults in 2020 during the early stages of COVID-19 [1]. The following year, the Tokyo Olympics were held in Japan, and a state of emergency was declared several times. Public perceptions of and risks associated with the pandemic changed drastically with the emergence of COVID-19 vaccines. We re-examine whether effect of hand-washing education on preventive behaviours persisted by covering a longer period of the COVID-19 pandemic than previous studies. 26 surveys were conducted nearly once a month for 30 months from March 2020 (the early stage of COVID-19) to September 2022 in Japan. By corresponding with the same individuals across surveys, we comprehensively gathered data on preventive behaviours during this period. In addition, we asked about hand-washing education they had received in their primary school. We used the data to investigate how and the degree to which school education is associated with pandemic mitigating preventive behaviours. We found that hand-washing education in primary school is positively associated with behaviours such as hand washing and mask wearing as a COVID-19 preventive measure, but not related to staying at home. We observed a statistically significant difference in hand washing between adults who received childhood hand-washing education and those who did not. This difference persisted throughout the study period. In comparison, the difference in mask wearing between the two groups was smaller, but still statistically significant. Furthermore, there was no difference in staying at home between them.


title: Generalized Munchausen Reinforcement Learning using Tsallis KL Divergence

abstract: Many policy optimization approaches in reinforcement learning incorporate a Kullback-Leilbler (KL) divergence to the previous policy, to prevent the policy from changing too quickly. This idea was initially proposed in a seminal paper on Conservative Policy Iteration, with approximations given by algorithms like TRPO and Munchausen Value Iteration (MVI). We continue this line of work by investigating a generalized KL divergence -- called the Tsallis KL divergence -- which use the $q$-logarithm in the definition. The approach is a strict generalization, as $q = 1$ corresponds to the standard KL divergence; $q > 1$ provides a range of new options. We characterize the types of policies learned under the Tsallis KL, and motivate when $q >1$ could be beneficial. To obtain a practical algorithm that incorporates Tsallis KL regularization, we extend MVI, which is one of the simplest approaches to incorporate KL regularization. We show that this generalized MVI($q$) obtains significant improvements over the standard MVI($q = 1$) across 35 Atari games.