title: Imaging the inner astronomical unit of Herbig Be star HD 190073

abstract: Inner regions of protoplanetary disks host many complex physical processes such as star-disk interactions, magnetic fields, planet formation, and the migration of new planets. To directly study this region requires milli-arcsecond angular resolution, beyond the diffraction limit of the world's largest optical telescopes and even too small for the mm-wave interferometer ALMA. However, we can use infrared interferometers to image the inner astronomical unit. Here, we present new results from the CHARA and VLTI arrays for the young and luminous Herbig Be star HD 190073. We detect a sub-AU cavity surrounded by a ring-like structure that we interpret as the dust destruction front. We model the shape with 6 radial profiles, 3 symmetric and 3 asymmetric, and present a model-free image reconstruction. All the models are consistent with a near face-on disk with inclination $\lesssim 20^\circ$, and we measure an average ring radius of 1.4 $\pm 0.2$ mas (1.14 AU). Around $48\%$ of the total flux comes from the disk with ~$15\%$ of that emission appearing to emerge from inside the inner rim. The cause of emission is still unclear, perhaps due to different dust grain compositions or gas emission. The skewed models and the imaging point to an off-center star, possibly due to binarity. Our image shows a sub-AU structure, which seems to move between the two epochs inconsistently with Keplerian motion and we discuss possible explanations for this apparent change.


title: Continuous phase transitions between fractional quantum Hall states and symmetry-protected topological states

abstract: We study quantum phase transitions in Bose-Fermi mixtures driven by interspecies interaction in the quantum Hall regime. In the absence of such an interaction, the bosons and fermions form their respective fractional quantum Hall (FQH) states at certain filling factors. A symmetry-protected topological (SPT) state is identified as the ground state for strong interspecies interaction. The phase transitions between them are proposed to be described by Chern-Simons-Higgs field theories. For a simple microscopic Hamiltonian, we present numerical evidence for the existence of the SPT state and a continuous transition to the FQH state. It is also found that the entanglement entropy between the bosons and fermions exhibits scaling behavior in the vicinity of this transition.


title: Universality of the Collins-Soper kernel in lattice calculations

abstract: The Collins-Soper (CS) kernel is a nonperturbative function that characterizes the rapidity evolution of transverse-momentum-dependent parton distribution functions (TMDPDFs) and wave functions. In this Letter, we calculate the CS kernel for pion and proton targets and for quasi-TMDPDFs of leading and next-to-leading power. The calculations are carried out on the CLS ensemble H101 with dynamical $N_f=2+1$ clover-improved Wilson fermions. Our analyses demonstrate the consistency of different lattice extractions of the CS kernel for mesons and baryons, as well as for twist-two and twist-three operators, even though lattice artifacts could be significant. This consistency corroborates the universality of the lattice-determined CS kernel and suggests that a high-precision determination of it is in reach.


title: Ground(less) Truth: A Causal Framework for Proxy Labels in Human-Algorithm Decision-Making

abstract: A growing literature on human-AI decision-making investigates strategies for combining human judgment with statistical models to improve decision-making. Research in this area often evaluates proposed improvements to models, interfaces, or workflows by demonstrating improved predictive performance on "ground truth" labels. However, this practice overlooks a key difference between human judgments and model predictions. Whereas humans reason about broader phenomena of interest in a decision -- including latent constructs that are not directly observable, such as disease status, the "toxicity" of online comments, or future "job performance" -- predictive models target proxy labels that are readily available in existing datasets. Predictive models' reliance on simplistic proxies makes them vulnerable to various sources of statistical bias. In this paper, we identify five sources of target variable bias that can impact the validity of proxy labels in human-AI decision-making tasks. We develop a causal framework to disentangle the relationship between each bias and clarify which are of concern in specific human-AI decision-making tasks. We demonstrate how our framework can be used to articulate implicit assumptions made in prior modeling work, and we recommend evaluation strategies for verifying whether these assumptions hold in practice. We then leverage our framework to re-examine the designs of prior human subjects experiments that investigate human-AI decision-making, finding that only a small fraction of studies examine factors related to target variable bias. We conclude by discussing opportunities to better address target variable bias in future research.


title: Preconditioned Score-based Generative Models

abstract: Score-based generative models (SGMs) have recently emerged as a promising class of generative models. However, a fundamental limitation is that their sampling process is slow due to a need for many (\eg, $2000$) iterations of sequential computations. An intuitive acceleration method is to reduce the sampling iterations which however causes severe performance degradation. We assault this problem to the ill-conditioned issues of the Langevin dynamics and reverse diffusion in the sampling process. Under this insight, we propose a model-agnostic {\bf\em preconditioned diffusion sampling} (PDS) method that leverages matrix preconditioning to alleviate the aforementioned problem. PDS alters the sampling process of a vanilla SGM at marginal extra computation cost, and without model retraining. Theoretically, we prove that PDS preserves the output distribution of the SGM, no risk of inducing systematical bias to the original sampling process. We further theoretically reveal a relation between the parameter of PDS and the sampling iterations,easing the parameter estimation under varying sampling iterations. Extensive experiments on various image datasets with a variety of resolutions and diversity validate that our PDS consistently accelerates off-the-shelf SGMs whilst maintaining the synthesis quality. In particular, PDS can accelerate by up to $29\times$ on more challenging high resolution (1024$\times$1024) image generation. Compared with the latest generative models (\eg, CLD-SGM, DDIM, and Analytic-DDIM), PDS can achieve the best sampling quality on CIFAR-10 at a FID score of 1.99. Our code is made publicly available to foster any further research https://github.com/fudan-zvg/PDS.


title: ASFAP Working Groups Activity Summary: Biophysics, Light Sources, Atomic and Molecular Physics, Condensed Matter and Materials Physics, and Earth Sciences

abstract: Various panel sessions were organized to highlight the activities of the African Strategy for Fundamental and Applied Physics (ASFAP) Working Groups during the second African Conference of Fundamental and Applied Physics (ACP2021) that was held in March 7-11, 2022. A joint session was devoted to highlight the activities assigned to the Light Sources, Accelerators, Biophysics, Earth Sciences, Atomic and Molecular Physics, and Condensed Matter and Materials Physics Working Groups. Major outcomes and recommendations are demonstrated and deliberated in this contribution.


title: A Myhill-Nerode Theorem for Generalized Automata, with Applications to Pattern Matching and Compression

abstract: The model of generalized automata, introduced by Eilenberg in 1974, allows representing a regular language more concisely than conventional automata by allowing edges to be labeled not only with characters, but also strings. Giammaresi and Montalbano introduced a notion of determinism for generalized automata [STACS 1995]. While generalized deterministic automata retain many properties of conventional deterministic automata, the uniqueness of a minimal generalized deterministic automaton is lost. In the first part of the paper, we show that the lack of uniqueness can be explained by introducing a set $ \mathcal{W(A)} $ associated with a generalized automaton $ \mathcal{A} $. By fixing $ \mathcal{W(A)} $, we are able to derive for the first time a full Myhill-Nerode theorem for generalized automata, which contains the textbook Myhill-Nerode theorem for conventional automata as a degenerate case. In the second part of the paper, we show that the set $ \mathcal{W(A)} $ leads to applications for pattern matching and data compression. Wheeler automata [TCS 2017, SODA 2020] are a popular class of automata that can be compactly stored using $ e \log \sigma (1 + o(1)) + O(e) $ bits ($ e $ being the number of edges, $ \sigma $ being the size of the alphabet) in such a way that pattern matching queries can be solved in $ \tilde{O}(m) $ time ($ m $ being the length of the pattern). In the paper, we show how to extend these results to generalized automata. More precisely, a Wheeler generalized automata can be stored using $ \mathfrak{e} \log \sigma (1 + o(1)) + O(e + rn) $ bits so that pattern matching queries can be solved in $ \tilde{O}(r m) $ time, where $ \mathfrak{e} $ is the total length of all edge labels, $ r $ is the maximum length of an edge label and $ n $ is the number of states.


title: N/MEMS Biosensors: An Introduction

abstract: In the 21st century, biosensors have gathered much wider attention than ever before, irrespective of the technology that promises to bring them forward. With the recent COVID-19 outbreak, the concern and efforts to restore global health and well-being are rising at an unprecedented rate. A requirement to develop precise, fast, point-of-care, reliable, easily disposable/reproducible and low-cost diagnostic tools have ascended. Biosensors form a primary element of hand-held medical kits, tools, products, and/or instruments. They have a very wide range of applications such as nearby environmental checks, detecting the onset of a disease, food quality, drug discovery, medicine dose control, and many more. Thischapter explains how Nano/Micro-Electro-Mechanical Systems (N/MEMS) can be enabling technology towards a sustainable, scalable, ultra-miniaturized, easy-to-use, energy efficient, and integrated bio/chemical sensing system. This study provides a deeper insight into the fundamentals, recent advances, and potential end applications of N/MEMS sensors and integrated systems to detect and measure the concentration of biological and/or chemical analytes. Transduction principle/s, materials, efficient designs including readout technique, and sensor performance are explained. This is followed by a discussion on how N/MEMS biosensors continue to evolve. The challenges and possible opportunities are also discussed.


title: Scattered polarized radiation of extrasolar circumplanetary rings

abstract: We have investigated the impact of circumplanetary rings consisting of spherical micrometer-sized particles on the net scattered light polarization of extrasolar gas giants. Using the three-dimensional Monte Carlo radiative transfer code POLARIS, we studied the impact of the macroscopic parameters that define the ring, such as its radius and inclination, and the chemical composition of the ring particles on the net scattered polarization. For the spherical ring particles, we applied the Mie scattering theory. We studied the flux and polarization of the scattered stellar radiation as a function of planetary phase angle and wavelength from the optical to the near-infrared. For the chosen grain size distribution, the dust particles in the ring show strong forward scattering at the considered wavelengths. Thus, the reflected flux of the planet dominates the total reflected and polarized flux at small phase angles. However, the scattered and polarized flux of the ring increase at large phase angles and exceeds the total reflected planetary flux. For large rings that contain silicate particles, the total reflected flux is dominated by the radiation scattered by the dust in the ring at all phase angles. As a result, the orientation of polarization is parallel to the scattering plane at small phase angles. In contrast, for a ring that contains water ice particles, the orientation of polarization is parallel to the scattering plane at large phase angles. Depending on the ring inclination and orientation, the total reflected and polarized flux show a specific distribution as well. Large particles show a strong polarization at large phase angles compared to smaller particles. For a Jupiter-like atmosphere that contains methane and aerosols, methane absorption features are missing in the spectrum of a ringed planet.


title: Blow-up of dynamically restricted critical norms near a potential Navier-Stokes singularity

abstract: In this paper we develop new methods to obtain regularity criteria for the three-dimensional Navier-Stokes equations in terms of dynamically restricted endpoint critical norms: the critical Lebesgue norm in general or the critical weak Lebesgue norm in the axisymmetric case. This type of results is inspired in particular by a work of Neustupa (2014), which handles certain non endpoint critical norms. Our work enables to have a better understanding of the nonlocal effect of the pressure on the regularity of the solutions.


title: Nonparametric estimation of multivariate hidden Markov models using tensor-product B-splines

abstract: For multivariate time series driven by underlying states, hidden Markov models (HMMs) constitute a powerful framework which can be flexibly tailored to the situation at hand. However, in practice it can be challenging to choose an adequate emission distribution for multivariate observation vectors. For example, the marginal data distribution may not immediately reveal the within-state distributional form, and also the different data streams may operate on different supports, rendering the common approach of using a multivariate normal distribution inadequate. Here we explore a nonparametric estimation of the emission distributions within a multivariate HMM based on tensor-product B-splines. In two simulation studies, we show the feasibility of our modelling approach and demonstrate potential pitfalls of inappropriate choices of parametric distributions. To illustrate the practical applicability, we present a case study where we use an HMM to model the bivariate time series comprising the lengths and angles of goalkeeper passes during the UEFA EURO 2020, investigating the effect of match dynamics on the teams' tactics.


title: Bi-objective risk-averse facility location using a subset-based representation of the conditional value-at-risk

abstract: For many real-world decision-making problems subject to uncertainty, it may be essential to deal with multiple and often conflicting objectives while taking the decision-makers' risk preferences into account. Conditional value-at-risk (CVaR) is a widely applied risk measure to address risk-averseness of the decision-makers. In this paper, we use the subset-based polyhedral representation of the CVaR to reformulate the bi-objective two-stage stochastic facility location problem presented in Nazemi et al. (2021). We propose an approximate cutting-plane method to deal with this more computationally challenging subset-based formulation. Then, the cutting plane method is embedded into the epsilon-constraint method, the balanced-box method, and a recently developed matheuristic method to address the bi-objective nature of the problem. Our computational results show the effectiveness of the proposed method. Finally, we discuss how incorporating an approximation of the subset-based polyhedral formulation affects the obtained solutions.


title: Near-Optimal Cryptographic Hardness of Agnostically Learning Halfspaces and ReLU Regression under Gaussian Marginals

abstract: We study the task of agnostically learning halfspaces under the Gaussian distribution. Specifically, given labeled examples $(\mathbf{x},y)$ from an unknown distribution on $\mathbb{R}^n \times \{ \pm 1\}$, whose marginal distribution on $\mathbf{x}$ is the standard Gaussian and the labels $y$ can be arbitrary, the goal is to output a hypothesis with 0-1 loss $\mathrm{OPT}+\epsilon$, where $\mathrm{OPT}$ is the 0-1 loss of the best-fitting halfspace. We prove a near-optimal computational hardness result for this task, under the widely believed sub-exponential time hardness of the Learning with Errors (LWE) problem. Prior hardness results are either qualitatively suboptimal or apply to restricted families of algorithms. Our techniques extend to yield near-optimal lower bounds for related problems, including ReLU regression.


title: DEPAS: De-novo Pathology Semantic Masks using a Generative Model

abstract: The integration of artificial intelligence into digital pathology has the potential to automate and improve various tasks, such as image analysis and diagnostic decision-making. Yet, the inherent variability of tissues, together with the need for image labeling, lead to biased datasets that limit the generalizability of algorithms trained on them. One of the emerging solutions for this challenge is synthetic histological images. However, debiasing real datasets require not only generating photorealistic images but also the ability to control the features within them. A common approach is to use generative methods that perform image translation between semantic masks that reflect prior knowledge of the tissue and a histological image. However, unlike other image domains, the complex structure of the tissue prevents a simple creation of histology semantic masks that are required as input to the image translation model, while semantic masks extracted from real images reduce the process's scalability. In this work, we introduce a scalable generative model, coined as DEPAS, that captures tissue structure and generates high-resolution semantic masks with state-of-the-art quality. We demonstrate the ability of DEPAS to generate realistic semantic maps of tissue for three types of organs: skin, prostate, and lung. Moreover, we show that these masks can be processed using a generative image translation model to produce photorealistic histology images of two types of cancer with two different types of staining techniques. Finally, we harness DEPAS to generate multi-label semantic masks that capture different cell types distributions and use them to produce histological images with on-demand cellular features. Overall, our work provides a state-of-the-art solution for the challenging task of generating synthetic histological images while controlling their semantic information in a scalable way.


title: Multiple Appropriate Facial Reaction Generation in Dyadic Interaction Settings: What, Why and How?

abstract: According to the Stimulus Organism Response (SOR) theory, all human behavioral reactions are stimulated by context, where people will process the received stimulus and produce an appropriate reaction. This implies that in a specific context for a given input stimulus, a person can react differently according to their internal state and other contextual factors. Analogously, in dyadic interactions, humans communicate using verbal and nonverbal cues, where a broad spectrum of listeners' non-verbal reactions might be appropriate for responding to a specific speaker behaviour. There already exists a body of work that investigated the problem of automatically generating an appropriate reaction for a given input. However, none attempted to automatically generate multiple appropriate reactions in the context of dyadic interactions and evaluate the appropriateness of those reactions using objective measures. This paper starts by defining the facial Multiple Appropriate Reaction Generation (fMARG) task for the first time in the literature and proposes a new set of objective evaluation metrics to evaluate the appropriateness of the generated reactions. The paper subsequently introduces a framework to predict, generate, and evaluate multiple appropriate facial reactions.


title: Indexed profunctors over 2-categories

abstract: We define the notion of an indexed profunctor over a 2-category, and use it to develop an abstract theory of limits. The theory subsumes (conical) limits, weighted limits, ends and Kan extensions. Results include an abstract version of the theorem that right adjoint functors preserve limits, and an abstract account of the phenomenon that the comparison arrow implicated in the preservation of a limit by a functor is natural in any functorial variable the functor happens to depend on. These results make an extensive use of the data and axioms of an indexed profunctor over a 2-category.


title: Two-point measurement energy statistics from particle scattering

abstract: We show that the energy statistics resulting from a two-point measurement of an isolated quantum system subject to a time-dependent driving protocol can be probed by subjecting the same system to a collision with a suitably prepared incoming particle. This means that the particle acts both as an external drive and as an energy measurement device for the quantum system and that energy fluctuations can be defined within a fully autonomous setting.


title: Learning a quantum channel from its steady-state

abstract: We present a scalable method for learning local non-unital quantum channels using only local expectation values measured on their steady state. Our method can be seen as the channel equivalent of the algorithms for learning local Hamiltonians from their ground states. Such channels can be efficiently implemented on the quantum computer using mid-circuit measurements or RESET gates. Using engineered dissipative dynamics, our method provides a simple way to check how well a given noise model describes the actual noisy dynamics within the hardware, when all qubits are actively used for a non-negligible amount of time. Moreover, given a parameterized noise model, our method can be applied to learn the underlying noise parameters for the entire system. We demonstrate our method numerically and experimentally on an IBMQ machine, and show that a full noise model can be verified and learned from Pauli measurements on a single state.


title: SelectionBias: An R Package for Bounding Selection Bias

abstract: Selection bias can occur when subjects are included or excluded in the analysis based upon some selection criteria for the study population. The bias can jeopardize the validity of the study and sensitivity analyses for assessing the effect of the selection are desired. One method of sensitivity analysis is to construct bounds for the bias. In this work, we present an R package that can be used to calculate two previously proposed bounds for selection bias for the causal relative risk and causal risk difference for both the total and the selected population. The first bound, derived by Smith and VanderWeele (SV), is based on values of sensitivity parameters that describe parts of the joint distribution of the outcome, treatment, selection indicator and unobserved variables. The second bound is based solely on the observed data, and is therefore referred to as an assumption free (AF) bound. In addition to a tutorial for the bounds and the R package, we derive additional properties for the SV bound. We show that the sensitivity parameters are variation independent and derive feasible regions for them. Furthermore, a bound is sharp if it is a priori known that the bias can be equal to the value of the bound, given the values of the selected sensitivity parameters. Conditions for the SV bound to be sharp in the selected subpopulation are provided based on the observed data. We illustrate both the R package and the properties of the bound with a simulated dataset that emulates a study where the effect of the zika virus on microcephaly in Brazil is investigated.


title: Functional imaging through scattering medium via fluorescence speckle demixing and localization

abstract: Recently, fluorescence-based optical techniques have emerged as a powerful tool to probe information in the mammalian brain. However, tissue heterogeneities prevent clear imaging of deep neuron bodies due to light scattering. While several up-to-date approaches based on ballistic light allow to retrieve information at shallow depths inside the brain, non-invasive localization and functional imaging at depth still remains a challenge. It was recently shown that functional signals from time-varying fluorescent emitters located behind scattering samples could be retrieved by using a matrix factorization algorithm. Here we show that the seemingly information-less, low-contrast fluorescent speckle patterns recovered by the algorithm can be used to locate each individual emitter, even in the presence of background fluorescence. We test our approach by imaging the temporal activity of large groups of fluorescent sources behind different scattering phantoms mimicking biological tissues, and through a brain slice with a thickness of ~200 micron.


title: Releasing Memory with Optimistic Access: A Hybrid Approach to Memory Reclamation and Allocation in Lock-Free Programs

abstract: Lock-free data structures are an important tool for the development of concurrent programs as they provide scalability, low latency and avoid deadlocks, livelocks and priority inversion. However, they require some sort of additional support to guarantee memory reclamation. The Optimistic Access (OA) method has most of the desired properties for memory reclamation, but since it allows memory to be accessed after being reclaimed, it is incompatible with the traditional memory management model. This renders it unable to release memory to the memory allocator/operating system, and, as such, it requires a complex memory recycling mechanism. In this paper, we extend the lock-free general purpose memory allocator LRMalloc to support the OA method. By doing so, we are able to simplify the memory reclamation method implementation and also allow memory to be reused by other parts of the same process. We further exploit the virtual memory system provided by the operating system and hardware in order to make it possible to release reclaimed memory to the operating system.


title: Density jump as a function of magnetic field strength for parallel collisionless shocks with anisotropic upstream pressure

abstract: The properties of collisionless shocks are frequently assessed in the magnetohydrodynamics (MHD) model. Yet, in a collisionless plasma, an ambient magnetic field can sustain a stable anisotropy in the upstream or the downstream, resulting in a departure from the MHD predicted behavior. We present a model allowing to derive the downstream anisotropy, hence the shock density jump, in terms of the upstream quantities. For simplicity, the case of a parallel shock in pair plasma is considered. Contrary to previous works where the upstream was assumed isotropic, here the upstream anisotropy $A=T_\perp/T_\parallel$ is a free parameter. The strong sonic shock regime is formally identical to the isotropic upstream case. Yet, for intermediate sonic Mach numbers, a variety of behaviors appear as a result of the anisotropy of the upstream.

categories: physics.plasm-ph astro-ph.HE astro-ph.SR physics.geo-ph


title: Certification of genuine time-bin and energy-time entanglement with integrated photonics

abstract: Time-bin (TB) and energy-time (ET) entanglement are crucial resources for long-distance quantum information processing. Recently, major efforts have been made to produce compact high-quality sources of TB/ET entangled photons based on solid-state integrated technologies. However, these attempts failed to close the so-called "post-selection loophole". Here, we present an integrated photonic general Bell-test chip for genuine (i.e., free of the post-selection loophole) TB and ET entanglement certification. We report a violation of a Bell inequality by more than 10 standard deviations using our device based on the "hug" interferometric scheme. The experiment also demonstrates that the hug scheme, previously exploited for ET entanglement, can also be used for genuine TB entanglement.


title: Transferable Deep Metric Learning for Clustering

abstract: Clustering in high dimension spaces is a difficult task; the usual distance metrics may no longer be appropriate under the curse of dimensionality. Indeed, the choice of the metric is crucial, and it is highly dependent on the dataset characteristics. However a single metric could be used to correctly perform clustering on multiple datasets of different domains. We propose to do so, providing a framework for learning a transferable metric. We show that we can learn a metric on a labelled dataset, then apply it to cluster a different dataset, using an embedding space that characterises a desired clustering in the generic sense. We learn and test such metrics on several datasets of variable complexity (synthetic, MNIST, SVHN, omniglot) and achieve results competitive with the state-of-the-art while using only a small number of labelled training datasets and shallow networks.


title: On cosmological Inflation In Palatini $F(R,\phi)$ Gravity

abstract: Single field inflationary models are investigated within Palatini quadratic gravity represented by $R+\alpha R^2$ along with a non-minimal coupling of the form $f(\phi) R$ between the inflaton field $\phi$ and the gravity. The treatment is performed in the Einstein frame, where the minimal coupling to gravity is recovered through conformal transformation. We consider various limits of the model with different inflationary scenarios characterized as canonical slow-roll inflation in the limit $\alpha \dot{\phi}^2\ll (1+f(\phi)) $, constant-roll k-inflation for $\alpha \ll 1$, and slow-roll K-inflation for$ \alpha \gg 1$ . A cosine and exponential potential are examined with the limits mentioned above and different well-motivated non-minimal couplings to gravity. We compare the theoretical results, exemplified by the tensor-to-scalar $r$ ratio and spectral index $n_s$, with the recent observational results of Planck 2018 $\&$ BICEP/Keck . Furthermore, we include the results of a new study forecast precision with which $n_s$ and $r$ can be constrained by currently envisaged observations, including CMB (Simons Observatory, CMB-S4, and LiteBIRD)


title: On finite generation in magnitude (co)homology, and its torsion

abstract: The aim of this paper is to apply the framework, which was developed by Sam and Snowden, to study structural properties of graph homologies, in the spirit of Ramos, Miyata and Proudfoot. Our main results concern the magnitude homology of graphs introduced by Hepworth and Willerton. More precisely, for graphs of bounded genus, we prove that magnitude cohomology, in each homological degree, has rank which grows at most polynomially in the number of vertices, and that its torsion is bounded. As a consequence, we obtain analogous results for path homology of (undirected) graphs. We complement the work with a proof that the category of planar graphs of bounded genus and marked edges, with contractions, is quasi-Gr\"obner.


title: Nonlocal-interaction vortices

abstract: We consider sequences of quadratic non-local functionals, depending on a small parameter $\e$, that approximate the Dirichlet integral by a well-known result by Bourgain, Brezis and Mironescu. Similarly to what is done for hard-core approximations to vortex energies in the case of the Dirichlet integral, we further scale such energies by $|\log\e|^{-1}$ and restrict them to $S^1$-valued functions. We introduce a notion of convergence of functions to integral currents with respect to which such energies are equi-coercive, and show the converge to a vortex energy, similarly to the limit behaviour of Ginzburg-Landau energies at the vortex scaling.


title: An Empirical Evaluation of Using Large Language Models for Automated Unit Test Generation

abstract: Unit tests play a key role in ensuring the correctness of software. However, manually creating unit tests is a laborious task, motivating the need for automation. Large Language Models (LLMs) have recently been applied to this problem, utilizing additional training or few-shot learning on examples of existing tests. This paper presents a large-scale empirical evaluation on the effectiveness of LLMs for automated unit test generation without additional training or manual effort, providing the LLM with the signature and implementation of the function under test, along with usage examples extracted from documentation. We also attempt to repair failed generated tests by re-prompting the model with the failing test and error message. We implement our approach in TestPilot, a test generation tool for JavaScript that automatically generates unit tests for all API functions in an npm package. We evaluate TestPilot using OpenAI's gpt3.5-turbo LLM on 25 npm packages with a total of 1,684 API functions. The generated tests achieve a median statement coverage of 70.2% and branch coverage of 52.8%, significantly improving on Nessie, a recent feedback-directed JavaScript test generation technique, which achieves only 51.3% statement coverage and 25.6% branch coverage. We also find that 92.8% of TestPilot's generated tests have no more than 50% similarity with existing tests (as measured by normalized edit distance), with none of them being exact copies. Finally, we run TestPilot with two additional LLMs, OpenAI's older code-cushman-002 LLM and the open LLM StarCoder. Overall, we observed similar results with the former (68.2% median statement coverage), and somewhat worse results with the latter (54.0% median statement coverage), suggesting that the effectiveness of the approach is influenced by the size and training set of the LLM, but does not fundamentally depend on the specific model.


title: Low-dimensional Data-based Surrogate Model of a Continuum-mechanical Musculoskeletal System Based on Non-intrusive Model Order Reduction

abstract: In recent decades, the main focus of computer modeling has been on supporting the design and development of engineering prototyes, but it is now ubiquitous in non-traditional areas such as medical rehabilitation. Conventional modeling approaches like the finite element~(FE) method are computationally costly when dealing with complex models, making them of limited use for purposes like real-time simulation or deployment on low-end hardware, if the model at hand cannot be simplified in a useful manner. Consequently, non-traditional approaches such as surrogate modeling using data-driven model order reduction are used to make complex high-fidelity models more widely available anyway. They often involve a dimensionality reduction step, in which the high-dimensional system state is transformed onto a low-dimensional subspace or manifold, and a regression approach to capture the reduced system behavior. While most publications focus on one dimensionality reduction, such as principal component analysis~(PCA) (linear) or autoencoder (nonlinear), we consider and compare PCA, kernel PCA, autoencoders, as well as variational autoencoders for the approximation of a structural dynamical system. In detail, we demonstrate the benefits of the surrogate modeling approach on a complex FE model of a human upper-arm. We consider both the models deformation and the internal stress as the two main quantities of interest in a FE context. By doing so we are able to create a computationally low cost surrogate model which captures the system behavior with high approximation quality and fast evaluations.


title: Unleashing the Power of Electrocardiograms: A novel approach for Patient Identification in Healthcare Systems with ECG Signals

abstract: Over the course of the past two decades, a substantial body of research has substantiated the viability of utilising cardiac signals as a biometric modality. This paper presents a novel approach for patient identification in healthcare systems using electrocardiogram signals. A convolutional neural network is used to classify users based on images extracted from ECG signals. The proposed identification system is evaluated in multiple databases, providing a comprehensive understanding of its potential in real-world scenarios. The impact of Cardiovascular Diseases on generic user identification has been largely overlooked in previous studies. The presented method takes into account the cardiovascular condition of the patients, ensuring that the results obtained are not biased or limited. Furthermore, the results obtained are consistent and reliable, with lower error rates and higher accuracy metrics, as demonstrated through extensive experimentation. All these features make the proposed method a valuable contribution to the field of patient identification in healthcare systems, and make it a strong contender for practical applications.


title: Reanalysis of the spin direction distribution of Galaxy Zoo SDSS spiral galaxies

abstract: The distribution of the spin directions of spiral galaxies in the Sloan Digital Sky Survey has been a topic of debate in the past two decades, with conflicting conclusions reported even in cases where the same data was used. Here we follow one of the previous experiments by applying the SpArcFiRe algorithm to annotate the spin directions in original dataset of Galaxy Zoo 1. The annotation of the galaxy spin directions is done after a first step of selecting the spiral galaxies in three different manners: manual analysis by Galaxy Zoo classifications, by a model-driven computer analysis, and with no selection of spiral galaxies. The results show that when spiral galaxies are selected by Galaxy Zoo volunteers, the distribution of their spin directions as determined by SpArcFiRe is not random, which agrees with previous reports. When selecting the spiral galaxies using a model-driven computer analysis or without selecting the spiral galaxies at all, the distribution is also not random. Simple binomial distribution analysis shows that the probability of the parity violation to occur by chance is lower than 0.01. Fitting the spin directions as observed from Earth to cosine dependence exhibits a dipole axis with statistical strength of 2.33$\sigma$ to 3.97$\sigma$. These experiments show that regardless of the selection mechanism and the analysis method, all experiments show similar conclusions. These results are aligned with previous reports using other methods and telescopes, suggesting that the spin directions of spiral galaxies as observed from Earth exhibit a dipole axis formed by their spin directions. Possible explanations can be related to the large-scale structure of the Universe, or to internal structure of galaxies. The catalogs of annotated galaxies generated as part of this study are available.


title: Generalized Weak Galerkin Finite Element Methods for Biharmonic Equations

abstract: The generalized weak Galerkin (gWG) finite element method is proposed and analyzed for the biharmonic equation. A new generalized discrete weak second order partial derivative is introduced in the gWG scheme to allow arbitrary combinations of piecewise polynomial functions defined in the interior and on the boundary of general polygonal or polyhedral elements. The error estimates are established for the numerical approximation in a discrete H^2 norm and a L^2 norm. The numerical results are reported to demonstrate the accuracy and flexibility of our proposed gWG method for the biharmonic equation.


title: Improper ferroelectricity in ultrathin hexagonal ferrites film

abstract: The suppression of ferroelectricity in ultrathin films of improper ferroelectric hexagonal ferrites or manganites has been attributed to the effect of interfacial clamping, however, the quantitative understanding and related phenomenological model are still lacking. In this work, we report the paraelectric-to-ferroelectric phase transition of epitaxial h-ScFeO3 films with different thickness through in-situ reflection high-energy electron diffraction (RHEED). Based on the interfacial clamping model and the Landau theory, we show that the thickness-dependence of the ferroelectric Curie temperature can be understood in terms of the characteristic length of interfacial clamping layer and the bulk Curie temperature. Furthermore, we found that the critical thickness of improper ferroelectricity is proportional to the characteristic length of interfacial clamping layer. These results reveal the essential role of mechanical clamping from interface on the improper ferroelectricity of hexagonal ferrites or manganites, and could serve as the guidance to achieve robust improper ferroelectricity in ultrathin films.


title: INSPIRE: INvestigating Stellar Population In RElics IV. The Initial Mass Function slope in relics

abstract: In the last decade, growing evidence has emerged supporting a non-universal stellar Initial Mass Function (IMF) in massive galaxies, with a larger number of dwarf stars with respect to the Milky-Way (bottom-heavy IMF). However, a consensus about the mechanisms that cause IMF variations is yet to be reached. Recently, it has been suggested that stars formed early-on in cosmic time, via a star formation burst, could be characterised by a bottom-heavy IMF. A promising way to confirm this is to use relics, ultra-compact massive galaxies, almost entirely composed by these "pristine" stars. The INSPIRE Project aims at assembling a large sample of confirmed relics, that can serve as laboratory to investigate on the conditions of star formation in the first 1-3 Gyr of the Universe. In this third INSPIRE paper, we build a high signal-to-noise spectrum from five relics and one from five galaxies with similar sizes, masses, and kinematical properties, but characterised by a more extended star formation history (non-relics). Our detailed stellar population analysis suggests a systematically bottom-heavier IMF slope for relics than for non-relics, adding new observational evidence for the non-universality of the IMF at various redshifts and further supporting the above proposed physical scenario.


title: Fourier-RNNs for Modelling Noisy Physics Data

abstract: Classical sequential models employed in time-series prediction rely on learning the mappings from the past to the future instances by way of a hidden state. The Hidden states characterise the historical information and encode the required temporal dependencies. However, most existing sequential models operate within finite-dimensional Euclidean spaces which offer limited functionality when employed in modelling physics relevant data. Alternatively recent work with neural operator learning within the Fourier space has shown efficient strategies for parameterising Partial Differential Equations (PDE). In this work, we propose a novel sequential model, built to handle Physics relevant data by way of amalgamating the conventional RNN architecture with that of the Fourier Neural Operators (FNO). The Fourier-RNN allows for learning the mappings from the input to the output as well as to the hidden state within the Fourier space associated with the temporal data. While the Fourier-RNN performs identical to the FNO when handling PDE data, it outperforms the FNO and the conventional RNN when deployed in modelling noisy, non-Markovian data.


title: Dynamical properties of coarse-grained linear SDEs

abstract: Coarse-graining or model reduction is a term describing a range of approaches used to extend the time-scale of molecular simulations by reducing the number of degrees of freedom. In the context of molecular simulation, standard coarse-graining approaches approximate the potential of mean force and use this to drive an effective Markovian model. To gain insight into this process, the simple case of a quadratic energy is studied in an overdamped setting. A hierarchy of reduced models is derived and analysed, and the merits of these different coarse-graining approaches are discussed. In particular, while standard recipes for model reduction accurately capture static equilibrium statistics, it is shown that dynamical statistics such as the mean-squared displacement display systematic error, even when a system exhibits large time-scale separation. In the linear setting studied, it is demonstrated both analytically and numerically that such models can be augmented in a simple way to better capture dynamical statistics.


title: Poisson points, resetting, universality and the role of the last item

abstract: For a stochastic process reset at random times, we discuss to what extent the probabilities of some orderings of observables associated with the intervals of time between resetting events are universal, i.e., independent of the choice of the observables, and in particular, to what extent universality depends on the choice of the distribution of these intervals. For Poissonian resetting, universality relies only on a combinatorial argument and on the statistical properties of Poisson points. For a generic distribution of time intervals between resets, universality no longer holds in general.


title: Harnessing the Power of Long-Range Entanglement for Clifford Circuit Synthesis

abstract: In superconducting architectures, limited connectivity remains a significant challenge for the synthesis and compilation of quantum circuits. We consider models of entanglement-assisted computation where long-range operations are achieved through injections of large GHZ states. These are prepared using ancillary qubits acting as an ``entanglement bus,'' unlocking global operation primitives such as multi-qubit Pauli rotations and fan out gates. We derive bounds on the circuit size for several well studied problems, such as CZ circuit, CX circuit, and Clifford circuit synthesis. In particular, in an architecture using one such entanglement bus, we give an $O(n^3)$-complexity synthesis scheme for arbitrary Clifford operations requiring at most $2n + 1$ layers of entangled-state-injections. In a square-lattice architecture with two entanglement buses, we show that a graph state can be synthesized using at most $\lceil \frac{1}{2}n \rceil + 1$ layers of GHZ state injections, and Clifford operations require only $\lceil\frac{3}{2} n \rceil + O(\sqrt n)$ layers of GHZ state injections.


title: Compact anisotropic stable hypersurfaces with free boundary in convex solid cones

abstract: We consider a convex solid cone $\mathcal{C}\subset\mathbb{R}^{n+1}$ with vertex at the origin and boundary $\partial\mathcal{C}$ smooth away from $0$. Our main result shows that a compact two-sided hypersurface $\Sigma$ immersed in $\mathcal{C}$ with free boundary in $\partial\mathcal{C}\setminus\{0\}$ and minimizing, up to second order, an anisotropic area functional under a volume constraint is contained in a Wulff-shape.


title: On the minimal faithful degree of Rhodes semisimple semigroups

abstract: In this paper we compute the minimum degree of a faithful representation by partial transformations of a finite semigroup admitting a faithful completely reducible matrix representation over the field of complex numbers. This includes all inverse semigroups, and hence our results generalize earlier results of Easdown and Schein on the minimal faithful degree of an inverse semigroup. It also includes well-studied monoids like full matrix monoids over finite fields and the monoid of binary relations (i.e., matrices over the Boolean semiring). Our answer reduces the computation to considerations of permutation representations of maximal subgroups that are faithful when restricted to distinguished normal subgroups. This is analogous to (and inspired by) recent results of the second author on the minimal number of irreducible constituents in a faithful completely reducible complex matrix representation of a finite semigroup that admits one. To illustrate what happens when a finite semigroup does not admit a faithful completely reducible representation, we compute the minimum faithful degree of the opposite monoid of the full transformation monoid.


title: Imitation from Observation With Bootstrapped Contrastive Learning

abstract: Imitation from observation (IfO) is a learning paradigm that consists of training autonomous agents in a Markov Decision Process (MDP) by observing expert demonstrations without access to its actions. These demonstrations could be sequences of environment states or raw visual observations of the environment. Recent work in IfO has focused on this problem in the case of observations of low-dimensional environment states, however, access to these highly-specific observations is unlikely in practice. In this paper, we adopt a challenging, but more realistic problem formulation, learning control policies that operate on a learned latent space with access only to visual demonstrations of an expert completing a task. We present BootIfOL, an IfO algorithm that aims to learn a reward function that takes an agent trajectory and compares it to an expert, providing rewards based on similarity to agent behavior and implicit goal. We consider this reward function to be a distance metric between trajectories of agent behavior and learn it via contrastive learning. The contrastive learning objective aims to closely represent expert trajectories and to distance them from non-expert trajectories. The set of non-expert trajectories used in contrastive learning is made progressively more complex by bootstrapping from roll-outs of the agent learned through RL using the current reward function. We evaluate our approach on a variety of control tasks showing that we can train effective policies using a limited number of demonstrative trajectories, greatly improving on prior approaches that consider raw observations.


title: Towards Agile Text Classifiers for Everyone

abstract: Text-based safety classifiers are widely used for content moderation and increasingly to tune generative language model behavior - a topic of growing concern for the safety of digital assistants and chatbots. However, different policies require different classifiers, and safety policies themselves improve from iteration and adaptation. This paper introduces and evaluates methods for agile text classification, whereby classifiers are trained using small, targeted datasets that can be quickly developed for a particular policy. Experimenting with 7 datasets from three safety-related domains, comprising 15 annotation schemes, led to our key finding: prompt-tuning large language models, like PaLM 62B, with a labeled dataset of as few as 80 examples can achieve state-of-the-art performance. We argue that this enables a paradigm shift for text classification, especially for models supporting safer online discourse. Instead of collecting millions of examples to attempt to create universal safety classifiers over months or years, classifiers could be tuned using small datasets, created by individuals or small organizations, tailored for specific use cases, and iterated on and adapted in the time-span of a day.


title: Fourier Neural Operator for Plasma Modelling

abstract: Predicting plasma evolution within a Tokamak is crucial to building a sustainable fusion reactor. Whether in the simulation space or within the experimental domain, the capability to forecast the spatio-temporal evolution of plasma field variables rapidly and accurately could improve active control methods on current tokamak devices and future fusion reactors. In this work, we demonstrate the utility of using Fourier Neural Operator (FNO) to model the plasma evolution in simulations and experiments. Our work shows that the FNO is capable of predicting magnetohydrodynamic models governing the plasma dynamics, 6 orders of magnitude faster than the traditional numerical solver, while maintaining considerable accuracy (NMSE $\sim 10^{-5})$. Our work also benchmarks the performance of the FNO against other standard surrogate models such as Conv-LSTM and U-Net and demonstrate that the FNO takes significantly less time to train, requires less parameters and outperforms other models. We extend the FNO approach to model the plasma evolution observed by the cameras positioned within the MAST spherical tokamak. We illustrate its capability in forecasting the formation of filaments within the plasma as well as the heat deposits. The FNO deployed to model the camera is capable of forecasting the full length of the plasma shot within half the time of the shot duration.


title: Direct electrocaloric characterization of ceramic films

abstract: Reliable and accurate characterization of the electrocaloric effect is necessary to understand the intrinsic properties of materials. To date, several methods have been developed to directly measure the electrocaloric effect. However, each of them has some limitations, making them less suitable for characterizing ceramic films, which rely almost exclusively on less accurate indirect methods. Here, a new approach is proposed to address the process of rapid heat dissipation in ceramic films and to detect the electrically induced temperature change before it thermally bonds with the surrounding elements. By using a polymer substrate that slows heat dissipation to the substrate and fast infrared imaging, a substantial part of the adiabatic electrocaloric effect in Pb(Mg1/3Nb2/3)O3-based ceramic films is captured. Infrared imaging provides a robust technique to reduce the ratio between the adiabatic and the measured electrocaloric temperature change in micrometer-sized ceramic films to a single-digit number, ~3.5. The obtained results are validated with another direct thermometric method and compared with the results obtained with an indirect approach. Despite different measurement principles, the results obtained with the two direct methods agree well. The proposed approach is timely and could open a door to verify the predicted giant electrocaloric effects in ceramic films, thus accelerating the process of their integration into functional devices.


title: Probabilistic Circuits That Know What They Don't Know

abstract: Probabilistic circuits (PCs) are models that allow exact and tractable probabilistic inference. In contrast to neural networks, they are often assumed to be well-calibrated and robust to out-of-distribution (OOD) data. In this paper, we show that PCs are in fact not robust to OOD data, i.e., they don't know what they don't know. We then show how this challenge can be overcome by model uncertainty quantification. To this end, we propose tractable dropout inference (TDI), an inference procedure to estimate uncertainty by deriving an analytical solution to Monte Carlo dropout (MCD) through variance propagation. Unlike MCD in neural networks, which comes at the cost of multiple network evaluations, TDI provides tractable sampling-free uncertainty estimates in a single forward pass. TDI improves the robustness of PCs to distribution shift and OOD data, demonstrated through a series of experiments evaluating the classification confidence and uncertainty estimates on real-world data.


title: 5D Elko spinor field non-minimally coupled to nonmetricity in $f(Q)$ gravity

abstract: This paper aims to investigate the localization of the five-dimensional spinor field known as Elko (dual-helicity eigenspinors of the charge conjugation operator) by employing a Yukawa-like geometrical coupling in which the Elko field is non-minimally coupled to nonmetricity scalar $Q$. We adopt the braneworld scenarios in which the first-order formalism with sine-Gordon and linear superpotentials is employed to obtain the warp factors. A linear function supports the zero-mode trapping within the geometric coupling, leading to the same effective potential as the scalar field. Moreover, an exotic term must be added to obtain real-valued massive modes. Such modes are investigated through the Schr\"{o}dinger-like approach.


title: Molecular hydrodynamic theory of the velocity autocorrelation function

abstract: The velocity autocorrelation function (VACF) encapsulates extensive information about a fluid's molecular-structural and hydrodynamic properties. We address the following fundamental question: How well can a purely hydrodynamic description recover the molecular features of a fluid as exhibited by the VACF? To this end, we formulate a bona fide hydrodynamic theory of the tagged-particle VACF for simple fluids. Our approach is distinguished from previous efforts in two key ways: collective hydrodynamic modes are modeled by \emph{linear} hydrodynamic equations; the fluid's static kinetic energy spectrum is identified as a necessary initial condition for the momentum current correlation. Our formulation leads to a natural physical interpretation of the hydrodynamic VACF as a superposition of quasinormal hydrodynamic modes weighted commensurately with the static kinetic energy spectrum, which appears to be essential to bridging continuum hydrodynamical behavior and discrete-particle kinetics. Our methodology yields VACF calculations quantitatively on par with existing approaches for liquid noble gases and alkali metals; moreover, our hydrodynamic model for the self-intermediate scattering function extends the applicable domain to low densities where the Schmidt number is of order unity, enabling calculations for gases and supercritical fluids.

categories: cond-mat.soft cond-mat.stat-mech physics.chem-ph physics.flu-dyn


title: Multi-Agent Path Integral Control for Interaction-Aware Motion Planning in Urban Canals

abstract: Autonomous vehicles that operate in urban environments shall comply with existing rules and reason about the interactions with other decision-making agents. In this paper, we introduce a decentralized and communication-free interaction-aware motion planner and apply it to Autonomous Surface Vessels (ASVs) in urban canals. We build upon a sampling-based method, namely Model Predictive Path Integral control (MPPI), and employ it to, in each time instance, compute both a collision-free trajectory for the vehicle and a prediction of other agents' trajectories, thus modeling interactions. To improve the method's efficiency in multi-agent scenarios, we introduce a two-stage sample evaluation strategy and define an appropriate cost function to achieve rule compliance. We evaluate this decentralized approach in simulations with multiple vessels in real scenarios extracted from Amsterdam's canals, showing superior performance than a state-of-the-art trajectory optimization framework and robustness when encountering different types of agents.


title: Automatic Noise Filtering with Dynamic Sparse Training in Deep Reinforcement Learning

abstract: Tomorrow's robots will need to distinguish useful information from noise when performing different tasks. A household robot for instance may continuously receive a plethora of information about the home, but needs to focus on just a small subset to successfully execute its current chore. Filtering distracting inputs that contain irrelevant data has received little attention in the reinforcement learning literature. To start resolving this, we formulate a problem setting in reinforcement learning called the $\textit{extremely noisy environment}$ (ENE), where up to $99\%$ of the input features are pure noise. Agents need to detect which features provide task-relevant information about the state of the environment. Consequently, we propose a new method termed $\textit{Automatic Noise Filtering}$ (ANF), which uses the principles of dynamic sparse training in synergy with various deep reinforcement learning algorithms. The sparse input layer learns to focus its connectivity on task-relevant features, such that ANF-SAC and ANF-TD3 outperform standard SAC and TD3 by a large margin, while using up to $95\%$ fewer weights. Furthermore, we devise a transfer learning setting for ENEs, by permuting all features of the environment after 1M timesteps to simulate the fact that other information sources can become relevant as the world evolves. Again, ANF surpasses the baselines in final performance and sample complexity. Our code is available at https://github.com/bramgrooten/automatic-noise-filtering


title: Between Generating Noise and Generating Images: Noise in the Correct Frequency Improves the Quality of Synthetic Histopathology Images for Digital Pathology

abstract: Artificial intelligence and machine learning techniques have the promise to revolutionize the field of digital pathology. However, these models demand considerable amounts of data, while the availability of unbiased training data is limited. Synthetic images can augment existing datasets, to improve and validate AI algorithms. Yet, controlling the exact distribution of cellular features within them is still challenging. One of the solutions is harnessing conditional generative adversarial networks that take a semantic mask as an input rather than a random noise. Unlike other domains, outlining the exact cellular structure of tissues is hard, and most of the input masks depict regions of cell types. However, using polygon-based masks introduce inherent artifacts within the synthetic images - due to the mismatch between the polygon size and the single-cell size. In this work, we show that introducing random single-pixel noise with the appropriate spatial frequency into a polygon semantic mask can dramatically improve the quality of the synthetic images. We used our platform to generate synthetic images of immunohistochemistry-treated lung biopsies. We test the quality of the images using a three-fold validation procedure. First, we show that adding the appropriate noise frequency yields 87% of the similarity metrics improvement that is obtained by adding the actual single-cell features. Second, we show that the synthetic images pass the Turing test. Finally, we show that adding these synthetic images to the train set improves AI performance in terms of PD-L1 semantic segmentation performances. Our work suggests a simple and powerful approach for generating synthetic data on demand to unbias limited datasets to improve the algorithms' accuracy and validate their robustness.


title: Towards fully bayesian analyses in Lattice QCD

abstract: We present a promising method to learn physical parameters from a bayesian inference, using modern tools to replace both our traditional fits and the way errors are computed and propagated. A few models are built as illustrations for a realistic case with Lattice QCD data, and appear to extract a lot of information with good stability. We discuss the evaluation of these models with either a fully bayesian approach or information criteria, as well as the model-building challenges which remain to be solved.


title: Primes in tuples and Romanoff's theorem

abstract: We obtain a lower bound for a number of primes in tuples. As applications, we obtain a lower bound for the Romanoff type representation functions.


title: The Ungar Games

abstract: Let $L$ be a finite lattice. An Ungar move sends an element $x\in L$ to the meet of $\{x\}\cup T$, where $T$ is a subset of the set of elements covered by $x$. We introduce the following Ungar game. Starting at the top element of $L$, two players -- Atniss and Eeta -- take turns making nontrivial Ungar moves; the first player who cannot do so loses the game. Atniss plays first. We say $L$ is an Atniss win (respectively, Eeta win) if Atniss (respectively, Eeta) has a winning strategy in the Ungar game on $L$. We first prove that the number of principal order ideals in the weak order on $S_n$ that are Eeta wins is $O(0.95586^nn!)$. We then consider a broad class of intervals in Young's lattice that includes all principal order ideals, and we characterize the Eeta wins in this class; we deduce precise enumerative results concerning order ideals in rectangles and type-$A$ root posets. We also characterize and enumerate principal order ideals in Tamari lattices that are Eeta wins. Finally, we conclude with some open problems and a short discussion of the computational complexity of Ungar games.


title: On the Iwasawa invariants of BDP Selmer groups and BDP P-adic L-fucntions

abstract: Let $p$ be an odd prime. Let $f_1$ and $f_2$ be weight-two Hecke eigen-cuspforms with isomorphic residual Galois representations at $p$. Greenberg--Vatsal and Emerton--Pollack--Weston showed that if $p$ is a good ordinary prime for the two forms, the Iwasawa invariants of their $p$-primary Selmer groups and $p$-adic $L$-functions over the cyclotomic $\mathbb{Z}_p$-extension of $\mathbb{Q}$ are closely related. The goal of this article is to generalize these results to the anticyclotomic setting. More precisely, let $K$ be an imaginary quadratic field where $p$ splits. Suppose that the generalized Heegner hypothesis holds with respect to both $(f_1,K)$ and $(f_2,K)$. We study relations between the Iwasawa invariants of the BDP Selmer groups and the BDP $p$-adic $L$-functions of $f_1$ and $f_2$.


title: Improving robot navigation in crowded environments using intrinsic rewards

abstract: Autonomous navigation in crowded environments is an open problem with many applications, essential for the coexistence of robots and humans in the smart cities of the future. In recent years, deep reinforcement learning approaches have proven to outperform model-based algorithms. Nevertheless, even though the results provided are promising, the works are not able to take advantage of the capabilities that their models offer. They usually get trapped in local optima in the training process, that prevent them from learning the optimal policy. They are not able to visit and interact with every possible state appropriately, such as with the states near the goal or near the dynamic obstacles. In this work, we propose using intrinsic rewards to balance between exploration and exploitation and explore depending on the uncertainty of the states instead of on the time the agent has been trained, encouraging the agent to get more curious about unknown states. We explain the benefits of the approach and compare it with other exploration algorithms that may be used for crowd navigation. Many simulation experiments are performed modifying several algorithms of the state-of-the-art, showing that the use of intrinsic rewards makes the robot learn faster and reach higher rewards and success rates (fewer collisions) in shorter navigation times, outperforming the state-of-the-art.


title: Implications of the Convergence of Language and Vision Model Geometries

abstract: Large-scale pretrained language models (LMs) are said to ``lack the ability to connect [their] utterances to the world'' (Bender and Koller, 2020). If so, we would expect LM representations to be unrelated to representations in computer vision models. To investigate this, we present an empirical evaluation across three different LMs (BERT, GPT2, and OPT) and three computer vision models (VMs, including ResNet, SegFormer, and MAE). Our experiments show that LMs converge towards representations that are partially isomorphic to those of VMs, with dispersion, and polysemy both factoring into the alignability of vision and language spaces. We discuss the implications of this finding.


title: VA-DepthNet: A Variational Approach to Single Image Depth Prediction

abstract: We introduce VA-DepthNet, a simple, effective, and accurate deep neural network approach for the single-image depth prediction (SIDP) problem. The proposed approach advocates using classical first-order variational constraints for this problem. While state-of-the-art deep neural network methods for SIDP learn the scene depth from images in a supervised setting, they often overlook the invaluable invariances and priors in the rigid scene space, such as the regularity of the scene. The paper's main contribution is to reveal the benefit of classical and well-founded variational constraints in the neural network design for the SIDP task. It is shown that imposing first-order variational constraints in the scene space together with popular encoder-decoder-based network architecture design provides excellent results for the supervised SIDP task. The imposed first-order variational constraint makes the network aware of the depth gradient in the scene space, i.e., regularity. The paper demonstrates the usefulness of the proposed approach via extensive evaluation and ablation analysis over several benchmark datasets, such as KITTI, NYU Depth V2, and SUN RGB-D. The VA-DepthNet at test time shows considerable improvements in depth prediction accuracy compared to the prior art and is accurate also at high-frequency regions in the scene space. At the time of writing this paper, our method -- labeled as VA-DepthNet, when tested on the KITTI depth-prediction evaluation set benchmarks, shows state-of-the-art results, and is the top-performing published approach.


title: Transient Hemodynamics Prediction Using an Efficient Octree-Based Deep Learning Model

abstract: Patient-specific hemodynamics assessment could support diagnosis and treatment of neurovascular diseases. Currently, conventional medical imaging modalities are not able to accurately acquire high-resolution hemodynamic information that would be required to assess complex neurovascular pathologies. Therefore, computational fluid dynamics (CFD) simulations can be applied to tomographic reconstructions to obtain clinically relevant information. However, three-dimensional (3D) CFD simulations require enormous computational resources and simulation-related expert knowledge that are usually not available in clinical environments. Recently, deep-learning-based methods have been proposed as CFD surrogates to improve computational efficiency. Nevertheless, the prediction of high-resolution transient CFD simulations for complex vascular geometries poses a challenge to conventional deep learning models. In this work, we present an architecture that is tailored to predict high-resolution (spatial and temporal) velocity fields for complex synthetic vascular geometries. For this, an octree-based spatial discretization is combined with an implicit neural function representation to efficiently handle the prediction of the 3D velocity field for each time step. The presented method is evaluated for the task of cerebral hemodynamics prediction before and during the injection of contrast agent in the internal carotid artery (ICA). Compared to CFD simulations, the velocity field can be estimated with a mean absolute error of 0.024 m/s, whereas the run time reduces from several hours on a high-performance cluster to a few seconds on a consumer graphical processing unit.


title: Boundedness of log Fano pairs with certain K-stability

abstract: We prove several boundedness results for log Fano pairs with certain K-stability. In particular, we prove that K-semistable log Fano pairs of Maeda type form a log bounded family. We also compute K-semistable domains for some examples.


title: Recommending to Strategic Users

abstract: Recommendation systems are pervasive in the digital economy. An important assumption in many deployed systems is that user consumption reflects user preferences in a static sense: users consume the content they like with no other considerations in mind. However, as we document in a large-scale online survey, users do choose content strategically to influence the types of content they get recommended in the future. We model this user behavior as a two-stage noisy signalling game between the recommendation system and users: the recommendation system initially commits to a recommendation policy, presents content to the users during a cold start phase which the users choose to strategically consume in order to affect the types of content they will be recommended in a recommendation phase. We show that in equilibrium, users engage in behaviors that accentuate their differences to users of different preference profiles. In addition, (statistical) minorities out of fear of losing their minority content exposition may not consume content that is liked by mainstream users. We next propose three interventions that may improve recommendation quality (both on average and for minorities) when taking into account strategic consumption: (1) Adopting a recommendation system policy that uses preferences from a prior, (2) Communicating to users that universally liked ("mainstream") content will not be used as basis of recommendation, and (3) Serving content that is personalized-enough yet expected to be liked in the beginning. Finally, we describe a methodology to inform applied theory modeling with survey results.


title: Large Scale Multi-Lingual Multi-Modal Summarization Dataset

abstract: Significant developments in techniques such as encoder-decoder models have enabled us to represent information comprising multiple modalities. This information can further enhance many downstream tasks in the field of information retrieval and natural language processing; however, improvements in multi-modal techniques and their performance evaluation require large-scale multi-modal data which offers sufficient diversity. Multi-lingual modeling for a variety of tasks like multi-modal summarization, text generation, and translation leverages information derived from high-quality multi-lingual annotated data. In this work, we present the current largest multi-lingual multi-modal summarization dataset (M3LS), and it consists of over a million instances of document-image pairs along with a professionally annotated multi-modal summary for each pair. It is derived from news articles published by British Broadcasting Corporation(BBC) over a decade and spans 20 languages, targeting diversity across five language roots, it is also the largest summarization dataset for 13 languages and consists of cross-lingual summarization data for 2 languages. We formally define the multi-lingual multi-modal summarization task utilizing our dataset and report baseline scores from various state-of-the-art summarization techniques in a multi-lingual setting. We also compare it with many similar datasets to analyze the uniqueness and difficulty of M3LS.


title: Gait design for limbless obstacle aided locomotion using geometric mechanics

abstract: Limbless robots have the potential to maneuver through cluttered environments that conventional robots cannot traverse. As illustrated in their biological counterparts such as snakes and nematodes, limbless locomotors can benefit from interactions with obstacles, yet such obstacle-aided locomotion (OAL) requires properly coordinated high-level self-deformation patterns (gait templates) as well as low-level body adaptation to environments. Most prior work on OAL utilized stereotyped traveling-wave gait templates and relied on local body deformations (e.g., passive body mechanics or decentralized controller parameter adaptation based on force feedback) for obstacle navigation, while gait template design for OAL remains less studied. In this paper, we explore novel gait templates for OAL based on tools derived from geometric mechanics (GM), which thus far has been limited to homogeneous environments. Here, we expand the scope of GM to obstacle-rich environments. Specifically, we establish a model that maps the presence of an obstacle to directional constraints in optimization. In doing so, we identify novel gait templates suitable for sparsely and densely distributed obstacle-rich environments respectively. Open-loop robophysical experiments verify the effectiveness of our identified OAL gaits in obstacle-rich environments. We posit that when such OAL gait templates are augmented with appropriate sensing and feedback controls, limbless locomotors will gain robust function in obstacle rich environments.


title: Uniform volume heating of mixed fuels within the ICF paradigm

abstract: The paper investigates the feasibility of achieving uniform high-power volume heating for a fusion reactor concept employing a mixed fuel composition involving $\text{pBDT}$. The realm of mixed fuel fusion concepts remains relatively unexplored. The pursuit of uniform high-power volume heating presents a technological challenge, yet it bears ramifications for fusion reactor designs. In this study, we introduce the proposition of employing embedded nano-structures that represent structured foams. These structured foams interact with short-pulse lasers, thereby achieving ultra-high power volume heating both within the fuel and the adjacent hohlraums. Notably, structured foams exhibit superior efficiency compared to unstructured foams, plasma or surfaces when it comes to absorbing high-power, short-pulse lasers. The suggested incorporation of these embedded structured foams interacting with an array of ultra-short laser pulses offers a high laser absorption power density, along with meticulous control over energy and power distribution within the fuel, both in spatial and temporal dimensions. This holds the potential for the realization of fusion reactors characterized by straight-forward designs and low complexity, where $Q_F \approx Q_T > 1$ is expected for the fuel and target gains. Depending on the fuel composition they can be strong neutron sources.


title: Efficient Systematic Deletions/Insertions of $0$'s Error Control Codes and the $L_{1}$ Metric (Extended version)

abstract: This paper gives some theory and efficient design of binary block systematic codes capable of controlling the deletions of the symbol ``$0$'' (referred to as $0$-deletions) and/or the insertions of the symbol ``$0$'' (referred to as $0$-insertions). The problem of controlling $0$-deletions and/or $0$-insertions (referred to as $0$-errors) is known to be equivalent to the efficient design of $L_{1}$ metric asymmetric error control codes over the natural alphabet, $\mathbb{N}$. So, $t$ $0$-insertion correcting codes can actually correct $t$ $0$-errors, detect $(t+1)$ $0$-errors and, simultaneously, detect all occurrences of only $0$-deletions or only $0$-insertions in every received word (briefly, they are $t$-Symmetric $0$-Error Correcting/$(t+1)$-Symmetric $0$-Error Detecting/All Unidirectional $0$-Error Detecting ($t$-Sy$0$EC/$(t+1)$-Sy$0$ED/AU$0$ED) codes). From the relations with the $L_{1}$ distance, optimal systematic code designs are given. In general, for all $t,k\in\mathbb{N}$, a recursive method is presented to encode $k$ information bits into efficient systematic $t$-Sy$0$EC/$(t+1)$-Sy$0$ED/AU$0$ED codes of length $$ n\leq k+t\log_{2}k+o(t\log n) $$ as $n\in\mathbb{N}$ increases. Decoding can be efficiently performed by algebraic means using the Extended Euclidean Algorithm (EEA).


title: A Domain Decomposition-Based CNN-DNN Architecture for Model Parallel Training Applied to Image Recognition Problems

abstract: Deep neural networks (DNNs) and, in particular, convolutional neural networks (CNNs) have brought significant advances in a wide range of modern computer application problems. However, the increasing availability of large amounts of datasets as well as the increasing available computational power of modern computers lead to a steady growth in the complexity and size of DNN and CNN models, and thus, to longer training times. Hence, various methods and attempts have been developed to accelerate and parallelize the training of complex network architectures. In this work, a novel CNN-DNN architecture is proposed that naturally supports a model parallel training strategy and that is loosely inspired by two-level domain decomposition methods (DDM). First, local CNN models, that is, subnetworks, are defined that operate on overlapping or nonoverlapping parts of the input data, for example, sub-images. The subnetworks can be trained completely in parallel. Each subnetwork outputs a local decision for the given machine learning problem which is exclusively based on the respective local input data. Subsequently, an additional DNN model is trained which evaluates the local decisions of the local subnetworks and generates a final, global decision. With respect to the analogy to DDM, the DNN can be interpreted as a coarse problem and hence, the new approach can be interpreted as a two-level domain decomposition. In this paper, solely image classification problems using CNNs are considered. Experimental results for different 2D image classification problems are provided as well as a face recognition problem, and a classification problem for 3D computer tomography (CT) scans. The results show that the proposed approach can significantly accelerate the required training time compared to the global model and, additionally, can also help to improve the accuracy of the underlying classification problem.


title: The Burst Observer and Optical Transient Exploring System in the multi-messenger astronomy era

abstract: The Burst Observer and Optical Transient Exploring System (BOOTES) was first designed as an asset of autonomous telescopes that started to be deployed in 1998, taking 24 years to be fully developed around the Earth. Nowadays BOOTES has became a global network of robotic telescopes, being the first one present in all continents, as of 2022. Here we present the details of the network and review its achievements over the last two decades regarding follow-up observations of high-energy transient events. Moreover, considering the recent operations of neutrino and gravitational wave detectors, some hot-topic expectations related to robotic astronomy are discussed within the framework of multi-wavelength astrophysics.


title: Characterizing the VPN Ecosystem in the Wild

abstract: With the shift to working remotely after the COVID-19 pandemic, the use of Virtual Private Networks (VPNs) around the world has nearly doubled. Therefore, measuring the traffic and security aspects of the VPN ecosystem is more important now than ever. It is, however, challenging to detect and characterize VPN traffic since some VPN protocols use the same port number as web traffic and port-based traffic classification will not help. VPN users are also concerned about the vulnerabilities of their VPN connections due to privacy issues. In this paper, we aim at detecting and characterizing VPN servers in the wild, which facilitates detecting the VPN traffic. To this end, we perform Internet-wide active measurements to find VPN servers in the wild, and characterize them based on their vulnerabilities, certificates, locations, and fingerprinting. We find 9.8M VPN servers distributed around the world using OpenVPN, SSTP, PPTP, and IPsec, and analyze their vulnerability. We find SSTP to be the most vulnerable protocol with more than 90% of detected servers being vulnerable to TLS downgrade attacks. Of all the servers that respond to our VPN probes, 2% also respond to HTTP probes and therefore are classified as Web servers. We apply our list of VPN servers to the traffic from a large European ISP and observe that 2.6% of all traffic is related to these VPN servers.


title: Poisson Cluster Process Models for Detecting Ultra-Diffuse Galaxies

abstract: We propose a novel set of Poisson Cluster Process (PCP) models to detect Ultra-Diffuse Galaxies (UDGs), a class of extremely faint, enigmatic galaxies of substantial interest in modern astrophysics. We model the unobserved UDG locations as parent points in a PCP, and infer their positions based on the observed spatial point patterns of their old star cluster systems. Many UDGs have somewhere from a few to hundreds of these old star clusters, which we treat as offspring points in our models. We also present a new framework to construct a marked PCP model using the marks of star clusters. The marked PCP model may enhance the detection of UDGs and offers broad applicability to problems in other disciplines. To assess the overall model performance, we design an innovative assessment tool for spatial prediction problems where only point-referenced ground truth is available, overcoming the limitation of standard ROC analyses where spatial Boolean reference maps are required. We construct a bespoke blocked Gibbs adaptive spatial birth-death-move MCMC algorithm to infer the locations of UDGs using real data from a \textit{Hubble Space Telescope} imaging survey. Based on our performance assessment tool, our novel models significantly outperform existing approaches using the Log-Gaussian Cox Process. We also obtained preliminary evidence that the marked PCP model improves UDG detection performance compared to the model without marks. Furthermore, we find evidence of a potential new ``dark galaxy'' that was not detected by previous methods.


title: Comp2Comp: Open-Source Body Composition Assessment on Computed Tomography

abstract: Computed tomography (CT) is routinely used in clinical practice to evaluate a wide variety of medical conditions. While CT scans provide diagnoses, they also offer the ability to extract quantitative body composition metrics to analyze tissue volume and quality. Extracting quantitative body composition measures manually from CT scans is a cumbersome and time-consuming task. Proprietary software has been developed recently to automate this process, but the closed-source nature impedes widespread use. There is a growing need for fully automated body composition software that is more accessible and easier to use, especially for clinicians and researchers who are not experts in medical image processing. To this end, we have built Comp2Comp, an open-source Python package for rapid and automated body composition analysis of CT scans. This package offers models, post-processing heuristics, body composition metrics, automated batching, and polychromatic visualizations. Comp2Comp currently computes body composition measures for bone, skeletal muscle, visceral adipose tissue, and subcutaneous adipose tissue on CT scans of the abdomen. We have created two pipelines for this purpose. The first pipeline computes vertebral measures, as well as muscle and adipose tissue measures, at the T12 - L5 vertebral levels from abdominal CT scans. The second pipeline computes muscle and adipose tissue measures on user-specified 2D axial slices. In this guide, we discuss the architecture of the Comp2Comp pipelines, provide usage instructions, and report internal and external validation results to measure the quality of segmentations and body composition measures. Comp2Comp can be found at https://github.com/StanfordMIMI/Comp2Comp.


title: Inferring Player Location in Sports Matches: Multi-Agent Spatial Imputation from Limited Observations

abstract: Understanding agent behaviour in Multi-Agent Systems (MAS) is an important problem in domains such as autonomous driving, disaster response, and sports analytics. Existing MAS problems typically use uniform timesteps with observations for all agents. In this work, we analyse the problem of agent location imputation, specifically posed in environments with non-uniform timesteps and limited agent observability (~95% missing values). Our approach uses Long Short-Term Memory and Graph Neural Network components to learn temporal and inter-agent patterns to predict the location of all agents at every timestep. We apply this to the domain of football (soccer) by imputing the location of all players in a game from sparse event data (e.g., shots and passes). Our model estimates player locations to within ~6.9m; a ~62% reduction in error from the best performing baseline. This approach facilitates downstream analysis tasks such as player physical metrics, player coverage, and team pitch control. Existing solutions to these tasks often require optical tracking data, which is expensive to obtain and only available to elite clubs. By imputing player locations from easy to obtain event data, we increase the accessibility of downstream tasks.


title: Beyond Uniform Smoothness: A Stopped Analysis of Adaptive SGD

abstract: This work considers the problem of finding a first-order stationary point of a non-convex function with potentially unbounded smoothness constant using a stochastic gradient oracle. We focus on the class of $(L_0,L_1)$-smooth functions proposed by Zhang et al. (ICLR'20). Empirical evidence suggests that these functions more closely captures practical machine learning problems as compared to the pervasive $L_0$-smoothness. This class is rich enough to include highly non-smooth functions, such as $\exp(L_1 x)$ which is $(0,\mathcal{O}(L_1))$-smooth. Despite the richness, an emerging line of works achieves the $\widetilde{\mathcal{O}}(\frac{1}{\sqrt{T}})$ rate of convergence when the noise of the stochastic gradients is deterministically and uniformly bounded. This noise restriction is not required in the $L_0$-smooth setting, and in many practical settings is either not satisfied, or results in weaker convergence rates with respect to the noise scaling of the convergence rate. We develop a technique that allows us to prove $\mathcal{O}(\frac{\mathrm{poly}\log(T)}{\sqrt{T}})$ convergence rates for $(L_0,L_1)$-smooth functions without assuming uniform bounds on the noise support. The key innovation behind our results is a carefully constructed stopping time $\tau$ which is simultaneously "large" on average, yet also allows us to treat the adaptive step sizes before $\tau$ as (roughly) independent of the gradients. For general $(L_0,L_1)$-smooth functions, our analysis requires the mild restriction that the multiplicative noise parameter $\sigma_1 < 1$. For a broad subclass of $(L_0,L_1)$-smooth functions, our convergence rate continues to hold when $\sigma_1 \geq 1$. By contrast, we prove that many algorithms analyzed by prior works on $(L_0,L_1)$-smooth optimization diverge with constant probability even for smooth and strongly-convex functions when $\sigma_1 > 1$.


title: Hamilton--Jacobi equations for controlled gradient flows: cylindrical test functions

abstract: This work is the second part of a program initiated in arXiv:2111.13258 aiming at the development of an intrinsic geometric well-posedness theory for Hamilton-Jacobi equations related to controlled gradient flow problems in metric spaces. Our main contribution is that of showing that the comparison principle proven therein implies a comparison principle for viscosity solutions relative to smoother Hamiltonians, acting on test functions that are mere cylindrical functions of the underling squared metric distance and whose rigorous definition is achieved from the Evolutional Variational Inequality formulation of gradient flows (EVI). In particular, the new Hamiltonians no longer require to work with test functions containing Tataru's distance. This substantial simplification paves the way for the development of a comprehensive existence theory.


title: EnergyShield: Provably-Safe Offloading of Neural Network Controllers for Energy Efficiency

abstract: To mitigate the high energy demand of Neural Network (NN) based Autonomous Driving Systems (ADSs), we consider the problem of offloading NN controllers from the ADS to nearby edge-computing infrastructure, but in such a way that formal vehicle safety properties are guaranteed. In particular, we propose the EnergyShield framework, which repurposes a controller ''shield'' as a low-power runtime safety monitor for the ADS vehicle. Specifically, the shield in EnergyShield provides not only safety interventions but also a formal, state-based quantification of the tolerable edge response time before vehicle safety is compromised. Using EnergyShield, an ADS can then save energy by wirelessly offloading NN computations to edge computers, while still maintaining a formal guarantee of safety until it receives a response (on-vehicle hardware provides a just-in-time fail safe). To validate the benefits of EnergyShield, we implemented and tested it in the Carla simulation environment. Our results show that EnergyShield maintains safe vehicle operation while providing significant energy savings compared to on-vehicle NN evaluation: from 24% to 54% less energy across a range of wireless conditions and edge delays.


title: Propagation of weakly advantageous mutations in cancer cell population

abstract: Research into somatic mutations in cancer cell DNA and their role in tumour growth and progression between successive stages is crucial for improving our understanding of cancer evolution. Mathematical and computer modelling can provide valuable insights into the scenarios of cancer growth, the roles of somatic mutations, and the types and strengths of evolutionary forces they introduce. Previous studies have developed mathematical models of cancer evolution, incorporating driver and passenger somatic mutations. Driver mutations were assumed to have a strong advantageous effect on the growth of the cancer cell population, while passenger mutations were considered fully neutral or mildly deleterious. However, according to several studies, passenger mutations may have a weakly advantageous effect on tumour growth. In this paper, we develop models of cancer evolution with somatic mutations that introduce a weakly advantageous force to the evolution of cancer cells. The models used in this study can be classified into two categories: deterministic and stochastic. Deterministic models are based on systems of differential equations that balance the average number of cells and mutations during evolution. To verify the results of our deterministic modelling, we use a stochastic model based on the Gillespie algorithm. We compare the predictions of our modelling with some observational data on cancer evolution.


title: Hierarchies among Genuine Multipartite Entangling Capabilities of Quantum Gates

abstract: We categorize quantum gates according to their capability to generate genuine multipartite entanglement based on the hierarchy of multipartite separable states. In particular, when a fixed unitary operator acts on the set of k-separable states, the maximal (average) genuine multipartite entanglement (GME) content produced via that particular unitary operator is determined after maximizing over the set of k-separable input states. We identify unitary operators that are beneficial for generating high GME when the input states are entangled in some bipartition, although the picture can also be reversed in which entanglement in inputs does not help. We characterize maximum entangling power of a variety of unitary operators including special classes of quantum gates, diagonal, permutation and Haar uniformly generated unitary operators by computing generalized geometric measure (GGM) as GME quantifier. We determine the unitary operators and their corresponding inputs which can create the resulting states having maximum GGM.


title: Superconducting diode effect in topological hybrid structures

abstract: Currently, the superconducting diode effect (SDE) is actively discussed due to large application potential in superconducting electronics. In particular, the superconducting hybrid structures based on three-dimensional (3D) topological insulators are among the best candidates due to the strongest spin-orbit coupling (SOC). Most of the theoretical studies of the SDE focus either on full numerical calculation, which is often rather complicated or on the phenomenological approach. In the present paper we perform a comparison of the linearized and nonlinear microscopic approaches in the superconductor/ ferromagnet/ 3D topological insulator (S/F/TI) hybrid structure. Employing the quasiclassical Green's function formalism we solve the problem self-consistently. We show that the results obtained by the linearized approximation are not qualitatively different from the nonlinear solution. Main distinction in the results between the two methods is quantitative, i. e. they yield different supercurrent amplitudes. However, when calculating the so-called diode quality factor the quantitative difference is eliminated and both approaches can result in a good agreement.


title: GFlowNet-EM for learning compositional latent variable models

abstract: Latent variable models (LVMs) with discrete compositional latents are an important but challenging setting due to a combinatorially large number of possible configurations of the latents. A key tradeoff in modeling the posteriors over latents is between expressivity and tractable optimization. For algorithms based on expectation-maximization (EM), the E-step is often intractable without restrictive approximations to the posterior. We propose the use of GFlowNets, algorithms for sampling from an unnormalized density by learning a stochastic policy for sequential construction of samples, for this intractable E-step. By training GFlowNets to sample from the posterior over latents, we take advantage of their strengths as amortized variational inference algorithms for complex distributions over discrete structures. Our approach, GFlowNet-EM, enables the training of expressive LVMs with discrete compositional latents, as shown by experiments on non-context-free grammar induction and on images using discrete variational autoencoders (VAEs) without conditional independence enforced in the encoder.


title: Post-Experiment Forensics and Human-in-the-Loop Interventions in Explainable Autonomous Scanning Probe Microscopy

abstract: The broad adoption of machine learning (ML)-based automated and autonomous experiments (AE) in physical characterization and synthesis requires development of strategies for understanding and intervention in the experimental workflow. Here, we introduce and realize strategies for post-acquisition forensic analysis applied to the deep kernel learning based AE scanning probe microscopy. This approach yields real-time and post-acquisition indicators of the progression of an active learning process interacting with an experimental system. We further illustrate that this approach can be extended towards human-in-the-loop autonomous experiments, where human operators make high-level decisions at high latencies setting the policies for AE, and the ML algorithm performs low-level fast decisions. The proposed approach is universal and can be extended to other physical and chemical imaging techniques and applications such as combinatorial library analysis. The full forensic analysis notebook is publicly available on GitHub at https://github.com/yongtaoliu/Forensics-DKL-BEPS.