title: HD 222925: a New Opportunity to Explore the Astrophysical and Nuclear Conditions of r-process Sites

abstract: With the most trans-iron elements detected of any star outside the Solar System, HD 222925 represents the most complete chemical inventory among r-process-enhanced, metal-poor stars. While the abundance pattern of the heaviest elements identified in HD 222925 agrees with the scaled Solar r-process residuals, as is characteristic of its r-process-enhanced classification, the newly measured lighter r-process elements display marked differences from their Solar counterparts. In this work, we explore which single astrophysical site (if any) produced the entire range of elements ($34\leq Z\leq 90$) present in HD 222925. We find that the abundance pattern of lighter r-process elements newly identified in HD 222925 presents a challenge for our existing nucleosynthesis models to reproduce. The most likely astrophysical explanation for the elemental pattern of HD 222925 is that its light r-process elements were created in rapidly expanding ejecta (e.g., from shocked, dynamical ejecta of compact object merger binaries). However, we find that the light r-process-element pattern can also be successfully reproduced by employing different nuclear mass models, indicating a need for a fresh investigation of nuclear input data for elements with $46\lesssim Z\lesssim 52$ by experimental methods. Either way, the new elemental abundance pattern of HD 222925 -- particularly the abundances obtained from space-based, ultraviolet (UV) data -- call for a deeper understanding of both astrophysical r-process sites and nuclear data. Similar UV observations of additional r-process-enhanced stars will be required to determine whether the elemental abundance pattern of HD 222925 is indeed a canonical template for the r-process at low metallicity.


title: Interpolated SelectionConv for Spherical Images and Surfaces

abstract: We present a new and general framework for convolutional neural network operations on spherical (or omnidirectional) images. Our approach represents the surface as a graph of connected points that doesn't rely on a particular sampling strategy. Additionally, by using an interpolated version of SelectionConv, we can operate on the sphere while using existing 2D CNNs and their weights. Since our method leverages existing graph implementations, it is also fast and can be fine-tuned efficiently. Our method is also general enough to be applied to any surface type, even those that are topologically non-simple. We demonstrate the effectiveness of our technique on the tasks of style transfer and segmentation for spheres as well as stylization for 3D meshes. We provide a thorough ablation study of the performance of various spherical sampling strategies.


title: Effective Constitutive Relations for Simulating CO2 Capillary Trapping in Heterogeneous Reservoirs with Fluvial Sedimentary Architecture

abstract: Carbon dioxide (CO2) storage reservoirs commonly exhibit sedimentary architecture that reflects fluvial deposition. The heterogeneity in petrophysical properties arising from this architecture influences the dynamics of injected CO2. We previously used a geocellular modeling approach to represent this heterogeneity, including heterogeneity in constitutive saturation relationships. The dynamics of CO2 plumes in fluvial reservoirs was investigated during and after injection. It was shown that small-scale (centimeter to meter) features play a critical role in capillary trapping processes and have a primary effect on physical- and dissolution-trapping of CO2, and on the ultimate distribution of CO2 in the reservoir. Heterogeneity in saturation functions at that small scale enhances capillary trapping (snap off), creates capillary pinning, and increases the surface area of the plume. The understanding of these small-scale trapping processes from previous work is here used to develop effective saturation relationships that represent, at a larger scale, the integral effect of these processes. While it is generally not computationally feasible to represent the small-scale heterogeneity directly in a typical reservoir simulation, the effective saturation relationships for capillary pressure and relative permeability presented here, along with an effective intrinsic permeability, allow better representation of the total physical trapping at the scale of larger model grid cells, as typically used in reservoir simulations, and thus the approach diminishes limits on cell size and decreases simulation time in reservoir simulations.


title: Proximal Learning With Opponent-Learning Awareness

abstract: Learning With Opponent-Learning Awareness (LOLA) (Foerster et al. [2018a]) is a multi-agent reinforcement learning algorithm that typically learns reciprocity-based cooperation in partially competitive environments. However, LOLA often fails to learn such behaviour on more complex policy spaces parameterized by neural networks, partly because the update rule is sensitive to the policy parameterization. This problem is especially pronounced in the opponent modeling setting, where the opponent's policy is unknown and must be inferred from observations; in such settings, LOLA is ill-specified because behaviorally equivalent opponent policies can result in non-equivalent updates. To address this shortcoming, we reinterpret LOLA as approximating a proximal operator, and then derive a new algorithm, proximal LOLA (POLA), which uses the proximal formulation directly. Unlike LOLA, the POLA updates are parameterization invariant, in the sense that when the proximal objective has a unique optimum, behaviorally equivalent policies result in behaviorally equivalent updates. We then present practical approximations to the ideal POLA update, which we evaluate in several partially competitive environments with function approximation and opponent modeling. This empirically demonstrates that POLA achieves reciprocity-based cooperation more reliably than LOLA.


title: Division algebraic symmetry breaking

abstract: Reframing certain well-known particle models in terms of normed division algebras leads to two new results for BSM physics. (1) We identify a sequence of complex structures which induces a cascade of breaking symmetries: Spin(10) $\mapsto$ Pati-Salam $\mapsto$ Left-Right symmetric $\mapsto$ Standard model + B-L (both pre- and post-Higgs-mechanism). These complex structures derive from the octonions, then from the quaternions, then from the complex numbers. (2) We provide, also for the first time we believe, an explicit demonstration of left-right symmetric Higgs representations stemming from quaternionic triality, tri($\mathbb{H}$). Upon the breaking of $su(2)_R$, our Higgs reduces to the familiar standard model Higgs.


title: Output Feedback Tube MPC-Guided Data Augmentation for Robust, Efficient Sensorimotor Policy Learning

abstract: Imitation learning (IL) can generate computationally efficient sensorimotor policies from demonstrations provided by computationally expensive model-based sensing and control algorithms. However, commonly employed IL methods are often data-inefficient, requiring the collection of a large number of demonstrations and producing policies with limited robustness to uncertainties. In this work, we combine IL with an output feedback robust tube model predictive controller (RTMPC) to co-generate demonstrations and a data augmentation strategy to efficiently learn neural network-based sensorimotor policies. Thanks to the augmented data, we reduce the computation time and the number of demonstrations needed by IL, while providing robustness to sensing and process uncertainty. We tailor our approach to the task of learning a trajectory tracking visuomotor policy for an aerial robot, leveraging a 3D mesh of the environment as part of the data augmentation process. We numerically demonstrate that our method can learn a robust visuomotor policy from a single demonstration--a two-orders of magnitude improvement in demonstration efficiency compared to existing IL methods.


title: Distributed MPC for Self-Organized Cooperation of Multiagent Systems -- Extended Version

abstract: We present a sequential distributed model predictive control (MPC) scheme for cooperative control of multi-agent systems with dynamically decoupled heterogeneous nonlinear agents subject to individual constraints. In the scheme, we explore the idea of using tracking MPC with artificial references to let agents coordinate their cooperation without external guidance. Each agent combines a tracking MPC with artificial references, the latter penalized by a suitable coupling cost. They solve an individual optimization problem for this artificial reference and an input that tracks it, only communicating the former to its neighbors in a communication graph. This puts the cooperative problem on a different layer than the handling of the dynamics and constraints, loosening the connection between the two. We provide sufficient conditions on the formulation of the cooperative problem and the coupling cost for the closed-loop system to asymptotically achieve it. Since the dynamics and the cooperative problem are only loosely connected, classical results from distributed optimization can be used to this end. We illustrate the scheme's application to consensus and formation control.


title: Absence of localization in two-dimensional Clifford circuits

abstract: We analyze a Floquet circuit with random Clifford gates in one and two spatial dimensions. By using random graphs and methods from percolation theory, we prove in the two dimensional setting that some local operators grow at ballistic rate, which implies the absence of localization. In contrast, the one-dimensional model displays a strong form of localization characterized by the emergence of left and right-blocking walls in random locations. We provide additional insights by complementing our analytical results with numerical simulations of operator spreading and entanglement growth, which show the absence (presence) of localization in two-dimension (one-dimension). Furthermore, we unveil that the spectral form factor of the Floquet unitary in two-dimensional circuits behaves like that of quasi-free fermions with chaotic single particle dynamics, with an exponential ramp that persists till times scaling linearly with the size of the system. Our work sheds light on the nature of disordered, Floquet Clifford dynamics and its relationship to fully chaotic quantum dynamics.


title: PERI: Part Aware Emotion Recognition In The Wild

abstract: Emotion recognition aims to interpret the emotional states of a person based on various inputs including audio, visual, and textual cues. This paper focuses on emotion recognition using visual features. To leverage the correlation between facial expression and the emotional state of a person, pioneering methods rely primarily on facial features. However, facial features are often unreliable in natural unconstrained scenarios, such as in crowded scenes, as the face lacks pixel resolution and contains artifacts due to occlusion and blur. To address this, in the wild emotion recognition exploits full-body person crops as well as the surrounding scene context. In a bid to use body pose for emotion recognition, such methods fail to realize the potential that facial expressions, when available, offer. Thus, the aim of this paper is two-fold. First, we demonstrate our method, PERI, to leverage both body pose and facial landmarks. We create part aware spatial (PAS) images by extracting key regions from the input image using a mask generated from both body pose and facial landmarks. This allows us to exploit body pose in addition to facial context whenever available. Second, to reason from the PAS images, we introduce context infusion (Cont-In) blocks. These blocks attend to part-specific information, and pass them onto the intermediate features of an emotion recognition network. Our approach is conceptually simple and can be applied to any existing emotion recognition method. We provide our results on the publicly available in the wild EMOTIC dataset. Compared to existing methods, PERI achieves superior performance and leads to significant improvements in the mAP of emotion categories, while decreasing Valence, Arousal and Dominance errors. Importantly, we observe that our method improves performance in both images with fully visible faces as well as in images with occluded or blurred faces.


title: Symmetric Homology is Representation Homology

abstract: Symmetric homology is a natural generalization of cyclic homology, in which symmetric groups play the role of cyclic groups. In the case of associative algebras, the symmetric homology theory was introduced by Z. Fiedorowicz \cite{F} and was further developed in the work of S. Ault \cite{Au1, Au2}. In this paper, we show that, for algebras defined over a field of characteristic $0$, the symmetric homology theory is naturally equivalent to the (one-dimensional) representation homology theory introduced by the authors (jointly with G. Khachatryan) in \cite{BKR}. Using known results on representation homology, we compute symmetric homology explicitly for basic algebras, such as polynomial algebras and universal enveloping algebras of (DG) Lie algebras. As an application, we prove two conjectures of Ault and Fiedorowicz, including the main conjecture of \cite{AF07} on topological interpretation of symmetric homology of polynomial algebras.


title: Impact of Electronic Correlations on High-Pressure Iron: Insights from Time-Dependent Density Functional Theory

abstract: We present a comprehensive investigation of the electrical and thermal conductivity of iron under high pressures at ambient temperature, employing the real-time formulation of time-dependent density functional theory (RT-TDDFT). Specifically, we examine the influence of a Hubbard correction (+\textit{U}) to account for strong electron correlations. Our calculations based on RT-TDDFT demonstrate that the evaluated electrical conductivity for both high-pressure body-centered cubic (BCC) and hexagonal close-packed (HCP) iron phases agrees well with experimental data. Furthermore, we explore the anisotropy in the thermal conductivity of HCP iron under high pressure, and our findings are consistent with experimental observations. Interestingly, we find that the incorporation of the +\textit{U} correction significantly impacts the ground state and linear response properties of iron at pressures below 50 GPa, with its influence diminishing as pressure increases. This study offers valuable insights into the influence of electronic correlations on the electronic transport properties of iron under extreme conditions.


title: Efficient Privacy-Preserving Machine Learning with Lightweight Trusted Hardware

abstract: In this paper, we propose a new secure machine learning inference platform assisted by a small dedicated security processor, which will be easier to protect and deploy compared to today's TEEs integrated into high-performance processors. Our platform provides three main advantages over the state-of-the-art: (i) We achieve significant performance improvements compared to state-of-the-art distributed Privacy-Preserving Machine Learning (PPML) protocols, with only a small security processor that is comparable to a discrete security chip such as the Trusted Platform Module (TPM) or on-chip security subsystems in SoCs similar to the Apple enclave processor. In the semi-honest setting with WAN/GPU, our scheme is 4X-63X faster than Falcon (PoPETs'21) and AriaNN (PoPETs'22) and 3.8X-12X more communication efficient. We achieve even higher performance improvements in the malicious setting. (ii) Our platform guarantees security with abort against malicious adversaries under honest majority assumption. (iii) Our technique is not limited by the size of secure memory in a TEE and can support high-capacity modern neural networks like ResNet18 and Transformer. While previous work investigated the use of high-performance TEEs in PPML, this work represents the first to show that even tiny secure hardware with really limited performance can be leveraged to significantly speed-up distributed PPML protocols if the protocol can be carefully designed for lightweight trusted hardware.


title: Intrinsic ion migration dynamics in a one-dimensional organic metal halide hybrid

abstract: Metal halide perovskites possess many physical properties amenable to optoelectronic applications, whereas the realization of these potentials has been hampered by their environmental and electronic instabilities. The morphological and molecular low dimensional perovskites and perovskite related materials have shown much promise in enhancing the chemical stability due to their unique molecular structures. Here we report on robust and reproducible four-terminal (4T) electrical measurements in a one-dimensional (1D) organic metal halide hybrid, (R)-{\alpha}-methylbenzylammonium lead triiodide ((R-{\alpha}-MBA)PbI3) made possible by its chemical stability. The results reveal a distinct intrinsic ion migration dynamic of single exponential, which underlies the unique 4T I-V characteristics. The dynamic is directly verified by real-time measurements of the transient ionic current. Our observations are consistent with photo-activation and field-assisted ion migration. The elucidated intrinsic ion dynamics may provide the physical basis for understanding and modelling the ubiquitous hysteresis in metal halides based electronic devices and new insights into the dynamics of ion migration in metal halide perovskites and hybrids in general.


title: Monochromatic spanning trees and matchings in ordered complete graphs

abstract: Two independent edges in ordered graphs can be nested, crossing or separated. These relations define six types of subgraphs, depending on which relations are forbidden. We refine a remark by Erd\H{o}s and Rado that every 2-coloring of the edges of a complete graph contains a monochromatic spanning tree. We show that forbidding one relation we always have a monochromatic (non-nested, non-crossing, non-separated) spanning tree in a 2-edge-colored ordered complete graph. On the other hand, if two relations are forbidden, then it is possible that we have monochromatic (nested, separated, crossing) subtrees of size only half the number of vertices. The existence of a monochromatic non-nested spanning tree in 2-colorings of ordered complete graphs verifies a more general conjecture for twisted drawings. Our second subject is to refine the Ramsey number of matchings for ordered complete graphs. Cockayne and Lorimer proved that for given positive integers t, n, m=(t-1)(n-1)+2n is the smallest integer with the following property: every t-coloring of the edges of a complete graph Km contains a monochromatic matching with n edges. We conjecture a strengthening: t-colored ordered complete graphs on m vertices contain monochromatic non-nested and also non-separated matchings with n edges. We prove this conjecture for some special cases. (i) Every t-colored ordered complete graph on t+3 vertices contains a monochromatic non-nested matching of size two. (ii) Every 2-colored ordered complete graph on 3n-1 vertices contains a monochromatic non-separated matching of size n. For nested, separated, and crossing matchings the situation is different. The smallest m ensuring a monochromatic matching of size n in every t-coloring is 2(t(n-1))+1) in the first two cases and one less in the third case.


title: Discipline Reputation Evaluation Based on PhD Exchange Network

abstract: When reputation evaluation indicators become targets, existing indicators will lose the role of indicating the true quality; At present, the evaluation of discipline reputation mostly focuses on subjective evaluation based on objective data, and there is a dispute about reliability and validity; Due to different indicators and weight settings, it is difficult to make horizontal comparison among disciplines; The evaluation also has a certain time lag. In order to solve the above four problems, this study explores a new method of discipline reputation evaluation. Taking the business administration discipline as an example, it collects data of 5848 doctoral graduates who first entered teaching posts, establishes a directed adjacency matrix from the employment unit to the doctoral degree awarding unit, and uses the theory and method of social network analysis to conduct quantitative analysis on the doctoral mutual employment network. The results show that: (1) PhD exchange network can explain discipline reputation and is a new indicator to measure discipline reputation; (2) From the perspective of employment behavior among colleges and universities, there is horizontal flow and downward flow between the head colleges and universities, and downward flow is mainly among the middle and lower colleges. There is a time lag between college talent recruitment and academic achievement output. Therefore, the mining of the structural characteristics and network evolution trend of the PhD exchange network based on the "foot voting" of doctoral graduates is faster than the discipline ranking based on the follow-up achievement indicators to reflect the changes in the discipline quality, which can be used to warn the changes in the discipline quality.


title: Testing the data framework for an AI algorithm in preparation for high data rate X-ray facilities

abstract: The advent of next-generation X-ray free electron lasers will be capable of delivering X-rays at a repetition rate approaching 1 MHz continuously. This will require the development of data systems to handle experiments at these type of facilities, especially for high throughput applications, such as femtosecond X-ray crystallography and X-ray photon fluctuation spectroscopy. Here, we demonstrate a framework which captures single shot X-ray data at the LCLS and implements a machine-learning algorithm to automatically extract the contrast parameter from the collected data. We measure the time required to return the results and assess the feasibility of using this framework at high data volume. We use this experiment to determine the feasibility of solutions for `live' data analysis at the MHz repetition rate.


title: Class-Level Confidence Based 3D Semi-Supervised Learning

abstract: Recent state-of-the-art method FlexMatch firstly demonstrated that correctly estimating learning status is crucial for semi-supervised learning (SSL). However, the estimation method proposed by FlexMatch does not take into account imbalanced data, which is the common case for 3D semi-supervised learning. To address this problem, we practically demonstrate that unlabeled data class-level confidence can represent the learning status in the 3D imbalanced dataset. Based on this finding, we present a novel class-level confidence based 3D SSL method. Firstly, a dynamic thresholding strategy is proposed to utilize more unlabeled data, especially for low learning status classes. Then, a re-sampling strategy is designed to avoid biasing toward high learning status classes, which dynamically changes the sampling probability of each class. To show the effectiveness of our method in 3D SSL tasks, we conduct extensive experiments on 3D SSL classification and detection tasks. Our method significantly outperforms state-of-the-art counterparts for both 3D SSL classification and detection tasks in all datasets.


title: Correlation between slip precursors and topological length scales at the onset of frictional sliding

abstract: Understanding the interplay between concurrent length scales is a fundamental issue in many problems involving friction between sliding interfaces, from tribology to the study of earthquakes and seismic faults. On the one hand, a macroscopic sliding event is preceded by slip precursors with a characteristic propagation length scale. On the other hand, the emergent frictional properties can be modified by surface patterning depending on their geometric length scale. This suggests that macroscopic sliding of structured surfaces is governed by the interplay between the length scale of the slip precursors and those characterizing the geometric features. In this paper, we investigate these aspects by means of numerical simulations using a two-dimensional spring-block model. We discuss the influence of the geometric features on the occurrence and localization of slip precursors, extending the study to interfaces characterized by two geometric length scales. We find that different types of detachment sequences are triggered by specific surface structures, depending on their scales and relation to sliding direction, leading to a macroscopically smooth transition to sliding in the case of hierarchical and/or anisotropic features. These concepts could be exploited in devices switching from static to dynamic sliding, and can contribute to an improvement in the understanding and interpretation of seismic data


title: Turbulent flows over porous lattices: alteration of near-wall turbulence and pore-flow amplitude modulation

abstract: Turbulent flows over porous lattices consisting of rectangular cuboid pores are investigated using scale-resolving direct numerical simulations. Beyond a certain threshold which is primarily determined by the wall-normal Darcy permeability, ${K_y}^+$, near-wall turbulence transitions from its canonical regime, marked by the presence of streak-like structures, to another marked by the presence of spanwise coherent structures reminiscent of the Kelvin-Helmholtz (K-H) type of instability. This permeability threshold agrees well with that previously established in studies where permeable-wall boundary conditions had been used as surrogates for a porous substrate. None of the substrates investigated demonstrate any drag reduction relative to smooth-wall turbulent flow. At the permeable surface, a significant component of the flow is that which adheres to the pore geometry and undergoes amplitude modulation (AM). This pore-coherent flow remains notable within the substrates, highlighting the importance of the porous substrate's microstructure when the overlying flow is turbulent, an aspect which cannot be accounted for when using continuum-based approaches to model porous media flows or effective representations such as wall boundary conditions. The severity of the AM is enhanced in the K-H-like regime, which has implications when designing porous substrates for transport processes. This suggests that the surface of the substrate can have a geometry which is different than the rest of it and tailored to influence the overlying flow in a particular way.


title: Analytical Study of Hybrid Surface Phonon-Plasmon-Polaritons (HSP3) in Symmetric Nonlinear-Graphene-hBN Heterostructures

abstract: The hybridization of hyperbolic Polaritonics with THz Plasmonics has attracted immense attention due to its fascinating applications in this region. However, to effectively enhance the performance of these coupled modes, one way is the usage of smart materials in the heterostructures. Here, we introduce a symmetric graphene-based structure containing hexagonal Boron Nitride (hBN) layers, where two nonlinear layers have been utilized as claddings to improve the confinement of Hybrid Surface Phonon-Plasmon-Polaritons. A new analytical model is derived. A set of nonlinear equations are solved numerically and the obtained results are reported. Harnessing a nonlinear medium together with hybrid graphene-hBN layers allows one to tune the propagation properties of the structure by changing the chemical potential, the relaxation time, and the nonlinear factor. The study has been done for two frequency regions: the upper Reststrahlen and the lower Reststrahlen bands. The numerical results show a long propagation length (Lprop=300 {\mu}m) and a large Figure of Merit (FOM=98) for HSP3 propagating in the upper Reststrahlen band at the frequency of 45 THz. The proposed structure and its analytical model can open new insight into the design of tunable devices in mid-infrared plasmonics.


title: Graph Attention Networks Unveil Determinants of Intra- and Inter-city Health Disparity

abstract: Understanding the determinants underlying variations in urban health status is important for informing urban design and planning, as well as public health policies. Multiple heterogeneous urban features could modulate the prevalence of diseases across different neighborhoods in cities and across different cities. This study examines heterogeneous features related to socio-demographics, population activity, mobility, and the built environment and their non-linear interactions to examine intra- and inter-city disparity in prevalence of four disease types: obesity, diabetes, cancer, and heart disease. Features related to population activity, mobility, and facility density are obtained from large-scale anonymized mobility data. These features are used in training and testing graph attention network (GAT) models to capture non-linear feature interactions as well as spatial interdependence among neighborhoods. We tested the models in five U.S. cities across the four disease types. The results show that the GAT model can predict the health status of people in neighborhoods based on the top five determinant features. The findings unveil that population activity and built-environment features along with socio-demographic features differentiate the health status of neighborhoods to such a great extent that a GAT model could predict the health status using these features with high accuracy. The results also show that the model trained on one city can predict health status in another city with high accuracy, allowing us to quantify the inter-city similarity and discrepancy in health status. The model and findings provide novel approaches and insights for urban designers, planners, and public health officials to better understand and improve health disparities in cities by considering the significant determinant features and their interactions.


title: Lattice-Based Quantum Advantage from Rotated Measurements

abstract: Trapdoor claw-free functions (TCFs) are immensely valuable in cryptographic interactions between a classical client and a quantum server. Typically, a protocol has the quantum server prepare a superposition of two-bit strings of a claw and then measure it using Pauli-$X$ or $Z$ measurements. In this paper, we demonstrate a new technique that uses the entire range of qubit measurements from the $XY$-plane. We show the advantage of this approach in two applications. First, building on (Brakerski et al. 2018, Kalai et al. 2022), we show an optimized two-round proof of quantumness whose security can be expressed directly in terms of the hardness of the LWE (learning with errors) problem. Second, we construct a one-round protocol for blind remote preparation of an arbitrary state on the $XY$-plane up to a Pauli-$Z$ correction.


title: Cross-Domain Aspect Extraction using Transformers Augmented with Knowledge Graphs

abstract: The extraction of aspect terms is a critical step in fine-grained sentiment analysis of text. Existing approaches for this task have yielded impressive results when the training and testing data are from the same domain. However, these methods show a drastic decrease in performance when applied to cross-domain settings where the domain of the testing data differs from that of the training data. To address this lack of extensibility and robustness, we propose a novel approach for automatically constructing domain-specific knowledge graphs that contain information relevant to the identification of aspect terms. We introduce a methodology for injecting information from these knowledge graphs into Transformer models, including two alternative mechanisms for knowledge insertion: via query enrichment and via manipulation of attention patterns. We demonstrate state-of-the-art performance on benchmark datasets for cross-domain aspect term extraction using our approach and investigate how the amount of external knowledge available to the Transformer impacts model performance.


title: Majorana coupling and Kondo screening of localized spins

abstract: We perform a theoretical analysis of the fate of local magnetic moment of a quantum dot coupled to a normal metallic lead and a topological superconducting wire hosting Majorana modes at the ends. By means of simple analytical tools and numerical renormalization group calculations we show that the proximity of Majorana mode reduces the magnetic moment from $1/4$, characteristic of a free spin $1/2$, to $1/16$. The coupling to the normal lead then causes the Kondo effect, such that the magnetic moment is fully screened below the Kondo temperature. The latter is vastly increased for strong coupling to Majorana mode.


title: Housing Forecasts via Stock Market Indicators

abstract: Through the reinterpretation of housing data as candlesticks, we extend Nature Scientific Reports' article by Liang and Unwin [LU22] on stock market indicators for COVID-19 data, and utilize some of the most prominent technical indicators from the stock market to estimate future changes in the housing market, comparing the findings to those one would obtain from studying real estate ETF's. By providing an analysis of MACD, RSI, and Candlestick indicators (Bullish Engulfing, Bearish Engulfing, Hanging Man, and Hammer), we exhibit their statistical significance in making predictions for USA data sets (using Zillow Housing data) and also consider their applications within three different scenarios: a stable housing market, a volatile housing market, and a saturated market. In particular, we show that bearish indicators have a much higher statistical significance then bullish indicators, and we further illustrate how in less stable or more populated countries, bearish trends are only slightly more statistically present compared to bullish trends.


title: TEFL: Turbo Explainable Federated Learning for 6G Trustworthy Zero-Touch Network Slicing

abstract: Sixth-generation (6G) networks anticipate intelligently supporting a massive number of coexisting and heterogeneous slices associated with various vertical use cases. Such a context urges the adoption of artificial intelligence (AI)-driven zero-touch management and orchestration (MANO) of the end-to-end (E2E) slices under stringent service level agreements (SLAs). Specifically, the trustworthiness of the AI black-boxes in real deployment can be achieved by explainable AI (XAI) tools to build transparency between the interacting actors in the slicing ecosystem, such as tenants, infrastructure providers and operators. Inspired by the turbo principle, this paper presents a novel iterative explainable federated learning (FL) approach where a constrained resource allocation model and an \emph{explainer} exchange -- in a closed loop (CL) fashion -- soft attributions of the features as well as inference predictions to achieve a transparent and SLA-aware zero-touch service management (ZSM) of 6G network slices at RAN-Edge setup under non-independent identically distributed (non-IID) datasets. In particular, we quantitatively validate the faithfulness of the explanations via the so-called attribution-based \emph{confidence metric} that is included as a constraint in the run-time FL optimization task. In this respect, Integrated-Gradient (IG) as well as Input $\times$ Gradient and SHAP are used to generate the attributions for the turbo explainable FL (TEFL), wherefore simulation results under different methods confirm its superiority over an unconstrained Integrated-Gradient \emph{post-hoc} FL baseline.


title: Bidiagonal Decompositions of Vandermonde-Type Matrices of Arbitrary Rank

abstract: We present a method to derive new explicit expressions for bidiagonal decompositions of Vandermonde and related matrices such as the (q-, h-) Bernstein-Vandermonde ones, among others. These results generalize the existing expressions for nonsingular matrices to matrices of arbitrary rank. For totally nonnegative matrices of the above classes, the new decompositions can be computed efficiently and to high relative accuracy componentwise in floating point arithmetic. In turn, matrix computations (e.g., eigenvalue computation) can also be performed efficiently and to high relative accuracy.


title: The weak null condition on Kerr backgrounds

abstract: We study a system of semilinear wave equations on Kerr backgrounds that satisfies the weak null condition. Under the assumption of small initial data, we prove global existence and pointwise decay estimates.


title: Joint User Localization and Location Calibration of A Hybrid Reconfigurable Intelligent Surface

abstract: The recent research in the emerging technology of reconfigurable intelligent surfaces (RISs) has identified its high potential for localization and sensing. However, to accurately localize a user placed in the area of influence of an RIS, the RIS location needs to be known a priori and its phase profile is required to be optimized for localization. In this paper, we study the problem of the joint localization of a hybrid RIS (HRIS) and a user, considering that the former is equipped with a single reception radio-frequency (RF) chain enabling simultaneous tunable reflections and sensing via power splitting. Focusing on the downlink of a multi-antenna base station, we present a multi-stage approach for the estimation of the HRIS position and orientation as well as the user position. Our simulation results, including comparisons with the Cram\'er-Rao lower bounds, demonstrate the efficiency of the proposed localization approach, while showcasing that there exists an optimal HRIS power splitting ratio for the desired multi-parameter estimation problem.


title: Dialogue system with humanoid robot

abstract: Today, as seen in smart speakers, spoken dialogue technology is rapidly advancing to enable human-like interaction. However, current dialogue systems cannot pay attention not only to the content of speech, but also to the way of speaking and eye contact and facial expressions, while watching the facial expressions of the person with whom one is speaking. Therefore, this study participated in a Japanese competition called the "Dialogue Robot Competition" and attempted to develop a dialogue system that includes control of not only the content of speech but also the robot's facial expressions and gaze in order to realize a humanoid robot that can naturally interact with humans.


title: Galois representations ramified at one prime and with suitably large image

abstract: Let $p\geq 7$ be a prime and $n>1$ be a natural number. We show that there exist infinitely many Galois representations $\varrho:Gal(\bar{\mathbb{Q}}/\mathbb{Q})\rightarrow GL_{n}(\mathbb{Z}_p)$ which are unramified outside $\{p, \infty\}$ with large image. More precisely, the Galois representations constructed have image containing the kernel of the mod-$p^t$ reduction map $SL_n(\mathbb{Z}_p)\rightarrow SL_n(\mathbb{Z}/p^t\mathbb{Z})$, where $t:=8(n^2-n)\left(3+\lfloor log_p(2^n+1)\rfloor\right)+8$. The results are proven via a purely Galois theoretic lifting construction. When $p\equiv 1\mod{4}$, our results are conditional since in this case, we assume a very weak version of Vandiver's conjecture.


title: Long-time behaviour of interaction models on Riemannian manifolds with bounded curvature

abstract: We investigate the long-time behaviour of solutions to a nonlocal partial differential equation on smooth Riemannian manifolds of bounded sectional curvature. The equation models self-collective behaviour with intrinsic interactions that are modelled by an interaction potential. We consider attractive interaction potentials and establish sufficient conditions for a consensus state to form asymptotically. In addition, we quantify the approach to consensus, by deriving a convergence rate for the diameter of the solution's support. The analytical results are supported by numerical simulations for the equation set up on the rotation group.


title: A tale of two theories of gravity in asymptotically Anti-de Sitter spacetime

abstract: We consider two BF formulations of the theory of gravity with a negative cosmological constant, of Plebanski and of MacDowell-Mansouri. Both give the standard Einstein equations in the bulk but differ in expressions of edge charges. We compute the asymptotic charges explicitly in both theories for AdS-Schwarzschild, AdS-Kerr, and AdS-Taub--NUT solutions. We find that while in the case of the Plebanski theory the charges are divergent, they are finite for MacDowell-Mansouri theory. Furthermore, we show that in the case of the arbitrary asymptotically AdS spacetimes, MacDowell--Mansouri asymptotic charges, action, and symplectic form are all finite. Therefore MacDowell-Mansouri theory of gravity in asymptotically AdS spaces does not need any counterterms.


title: Spatial Heterogeneity Localizes Turing Patterns in Reaction-Cross-Diffusion Systems

abstract: Motivated by bacterial chemotaxis and multi-species ecological interactions in heterogeneous environments, we study a general one-dimensional reaction-cross-diffusion system in the presence of spatial heterogeneity in both transport and reaction terms. Under a suitable asymptotic assumption that the transport is slow over the domain, while gradients in the reaction heterogeneity are not too sharp, we study the stability of a heterogeneous steady state approximated by the system in the absence of transport. Using a WKB ansatz, we find that this steady state can undergo a Turing-type instability in subsets of the domain, leading to the formation of localized patterns. The boundaries of the pattern-forming regions are given asymptotically by `local' Turing conditions corresponding to a spatially homogeneous analysis parameterized by the spatial variable. We developed a general open-source code which is freely available, and show numerical examples of this localized pattern formation in a Schnakenberg cross-diffusion system, a Keller-Segel chemotaxis model, and the Shigesada-Kawasaki-Teramoto model with heterogeneous parameters. We numerically show that the patterns may undergo secondary instabilities leading to spatiotemporal movement of spikes, though these remain approximately within the asymptotically predicted localized regions. This theory can elegantly differentiate between spatial structure due to background heterogeneity, from spatial patterns emergent from Turing-type instabilities.


title: An intergalactic medium temperature from a giant radio galaxy

abstract: The warm-hot intergalactic medium (warm-hot IGM, or WHIM) pervades the filaments of the Cosmic Web and harbours half of the Universe's baryons. The WHIM's thermodynamic properties are notoriously hard to measure. Here we estimate a galaxy group - WHIM boundary temperature using a new method. In particular, we use a radio image of the giant radio galaxy (giant RG, or GRG) created by NGC 6185, a massive nearby spiral. We analyse this extraordinary object with a Bayesian 3D lobe model and deduce an equipartition pressure $P_\mathrm{eq} = 6 \cdot 10^{-16}\ \mathrm{Pa}$ -- among the lowest found in RGs yet. Using an X-ray-based statistical conversion for Fanaroff-Riley II RGs, we find a true lobe pressure $P = 1.5\substack{+1.7\\-0.4}\cdot 10^{-15}\ \mathrm{Pa}$. Cosmic Web reconstructions, group catalogues, and MHD simulations furthermore imply an $\mathrm{Mpc}$-scale IGM density $1 + \delta_\mathrm{IGM} = 40\substack{+30\\-10}$. The buoyantly rising lobes are crushed by the IGM at their inner side, where an approximate balance between IGM and lobe pressure occurs: $P_\mathrm{IGM} \approx P$. The ideal gas law then suggests an IGM temperature $T_\mathrm{IGM} = 11\substack{+12\\-5} \cdot 10^6\ \mathrm{K}$, or $k_\mathrm{B}T_\mathrm{IGM} = 0.9\substack{+1.0\\-0.4}\ \mathrm{keV}$, at the virial radius -- consistent with X-ray-derived temperatures of similarly massive groups. Interestingly, the method is not performing at its limit: in principle, estimates $T_\mathrm{IGM} \sim 4 \cdot 10^6\ \mathrm{K}$ are already possible -- rivalling the lowest X-ray measurements available. The technique's future scope extends from galaxy group outskirts to the WHIM. In conclusion, we demonstrate that observations of GRGs in Cosmic Web filaments are finally sensitive enough to probe the thermodynamics of galaxy groups and beyond.


title: Optimized parametric inference for the inner loop of the Multigrid Ensemble Kalman Filter

abstract: Essential features of the Multigrid Ensemble Kalman Filter (G. Moldovan, G. Lehnasch, L. Cordier, M. Meldi, A multigrid/ensemble Kalman filter strategy for assimilation of unsteady flows, Journal of Computational Physics 443-110481) recently proposed for Data Assimilation of fluid flows are investigated and assessed in this article. The analysis is focused on the improvement in performance due to the inner loop. In this step, data from solutions calculated on the higher resolution levels of the multigrid approach are used as surrogate observations to improve the model prediction on the coarsest levels of the grid. The latter represents the level of resolution used to run the ensemble members for global Data Assimilation. The method is tested over two classical one-dimensional problems, namely the linear advection problem and the Burgers' equation. The analyses encompass a number of different aspects, such as different grid resolutions. The results indicate that the contribution of the inner loop is essential in obtaining accurate flow reconstruction and global parametric optimization. These findings open exciting perspectives of application to grid-dependent reduced-order models extensively used in fluid mechanics applications for complex flows, such as Large Eddy Simulation (LES).


title: Degrees of the stretched Kostka quasi-polynomials

abstract: We provide a type-uniform formula for the degree of the stretched Kostka quasi-polynomial $K_{\lambda,\mu}(N)$ in all classical types, improving a previous result by McAllister in $\mathfrak{sl}_r(\mathbb{C})$. Our proof relies on a combinatorial model for the weight multiplicity by Berenstein and Zelevinsky.


title: Note on the density of ISE and a related diffusion

abstract: The integrated super-Brownian excursion (ISE) is the occupation measure of the spatial component of the head of the Brownian snake with lifetime process the normalized Brownian excursion. It is a random probability measure on $\mathbb{R}$, and it is known to describe the continuum limit of the distribution of labels in various models of random discrete labelled trees. We show that $f_{ISE}$, its (random) density has a.s. a derivative $f'_{ISE}$ which is continuous and $\left(\frac{1}{2}-a\right)$-H\"older for any $a >0$ but for no $a<0$ (proving a conjecture of Bousquet-M\'elou and Janson). We conjecture that $f_{ISE}$ can be represented as a second-order diffusion of the form $$df'_{ISE}(t) = \sqrt{2f_{ISE}(t)}\, dB_t + g\left(f'_{ISE}(t), f_{ISE}(t),\int_{-\infty}^t f_{ISE}(s)ds\right)dt,$$ for some continuous function $g$, for $t>0$, and we give a number of remarks and questions in that direction. The proof of regularity is based on a moment estimate coming from a discrete model of trees, while the heuristic of the diffusion comes from an analogous statement in the discrete setting, which is a reformulation of explicit product formulas of Bousquet-M\'elou and the first author (2012).


title: Uncertainty in Extreme Multi-label Classification

abstract: Uncertainty quantification is one of the most crucial tasks to obtain trustworthy and reliable machine learning models for decision making. However, most research in this domain has only focused on problems with small label spaces and ignored eXtreme Multi-label Classification (XMC), which is an essential task in the era of big data for web-scale machine learning applications. Moreover, enormous label spaces could also lead to noisy retrieval results and intractable computational challenges for uncertainty quantification. In this paper, we aim to investigate general uncertainty quantification approaches for tree-based XMC models with a probabilistic ensemble-based framework. In particular, we analyze label-level and instance-level uncertainty in XMC, and propose a general approximation framework based on beam search to efficiently estimate the uncertainty with a theoretical guarantee under long-tail XMC predictions. Empirical studies on six large-scale real-world datasets show that our framework not only outperforms single models in predictive performance, but also can serve as strong uncertainty-based baselines for label misclassification and out-of-distribution detection, with significant speedup. Besides, our framework can further yield better state-of-the-art results based on deep XMC models with uncertainty quantification.


title: Nonparametric Quantile Regression: Non-Crossing Constraints and Conformal Prediction

abstract: We propose a nonparametric quantile regression method using deep neural networks with a rectified linear unit penalty function to avoid quantile crossing. This penalty function is computationally feasible for enforcing non-crossing constraints in multi-dimensional nonparametric quantile regression. We establish non-asymptotic upper bounds for the excess risk of the proposed nonparametric quantile regression function estimators. Our error bounds achieve optimal minimax rate of convergence for the Holder class, and the prefactors of the error bounds depend polynomially on the dimension of the predictor, instead of exponentially. Based on the proposed non-crossing penalized deep quantile regression, we construct conformal prediction intervals that are fully adaptive to heterogeneity. The proposed prediction interval is shown to have good properties in terms of validity and accuracy under reasonable conditions. We also derive non-asymptotic upper bounds for the difference of the lengths between the proposed non-crossing conformal prediction interval and the theoretically oracle prediction interval. Numerical experiments including simulation studies and a real data example are conducted to demonstrate the effectiveness of the proposed method.


title: Computational pathology in renal disease: a comprehensive perspective

abstract: Computational pathology is a field that has complemented various subspecialties of diagnostic pathology over the last few years. In this article a brief analyzis the different applications in nephrology is developed. To begin, an overview of the different forms of image production is provided. To continue, the most frequent applications of computer vision models, the salient features of the different clinical applications, and the data protection considerations encountered are described. To finish the development, I delve into the interpretability of these applications, expanding in depth on the three dimensions of this area.


title: MedCLIP: Contrastive Learning from Unpaired Medical Images and Text

abstract: Existing vision-text contrastive learning like CLIP aims to match the paired image and caption embeddings while pushing others apart, which improves representation transferability and supports zero-shot prediction. However, medical image-text datasets are orders of magnitude below the general images and captions from the internet. Moreover, previous methods encounter many false negatives, i.e., images and reports from separate patients probably carry the same semantics but are wrongly treated as negatives. In this paper, we decouple images and texts for multimodal contrastive learning thus scaling the usable training data in a combinatorial magnitude with low cost. We also propose to replace the InfoNCE loss with semantic matching loss based on medical knowledge to eliminate false negatives in contrastive learning. We prove that MedCLIP is a simple yet effective framework: it outperforms state-of-the-art methods on zero-shot prediction, supervised classification, and image-text retrieval. Surprisingly, we observe that with only 20K pre-training data, MedCLIP wins over the state-of-the-art method (using around 200K data). Our code is available at https://github.com/RyanWangZf/MedCLIP.


title: Optimized Telecloning Circuits: Theory and Practice of Nine NISQ Clones

abstract: Although perfect copying of an unknown quantum state is not possible, approximate cloning is possible in quantum mechanics. Quantum telecloning is a variant of approximate quantum cloning which uses quantum teleportation to allow for the use of classical communication to create physically separate clones of a quantum state. We present results of a of $1 \rightarrow 9$ universal, symmetric, optimal quantum telecloning implementation on a cloud accessible quantum computer - the Quantinuum H1-1 device. The H1-1 device allows direct creation of the telecloning protocol due to real time classical if-statements that are conditional on the mid-circuit measurement outcome of a Bell measurement. In this implementation, we also provide an improvement over previous work for the circuit model description of quantum telecloning, which reduces the required gate depth and gate count for an all-to-all connectivity. The demonstration of creating $9$ approximate clones on a quantum processor is the largest number of clones that has been generated, telecloning or otherwise.


title: Episodic accretion and mergers during growth of massive protostars

abstract: 3D simulations of high mass young stellar object (HMYSO) growth show that their circumstellar discs fragment onto multiple self-gravitating objects. Accretion of these by HMYSO may explain episodic accretion bursts discovered recently. We post-process results of a previous 3D simulation of a HMYSO disc with a 1D code that resolves the disc and object dynamics down to the stellar surface. We find that burst-like deposition of material into the inner disc seen in 3D simulations by itself does not always signify powerful accretion bursts. Only high density post-collapse clumps crossing the inner computational boundary may result in observable bursts. The rich physics of the inner disc has a significant impact on the expected accretion bursts: (1) In the standard turbulent viscosity discs, migrating objects can stall at a migration trap at the distance of a few au from the star. However, in discs powered by magnetised winds, the objects are able to cross the trap and produce bursts akin to those observed so far. (2) Migrating objects may interact with and modify the thermal (hydrogen ionisation) instability of the inner disc, which can be responsible for longer duration and lower luminosity bursts in HMYSOs. (3) If the central star is bloated to a fraction of an au by a previous episode of high accretion rate, or if the migrating object is particularly dense, a merger rather than a disc-mediated accretion burst results; (4) Object disruption bursts may be super-Eddington, leading to episodic feedback on HMYSO surroundings via powerful outflows.


title: Informed Neutrality in Minimalist Market Design: A Case Study on a Constitutional Crisis in India

abstract: In a 3-2 split verdict, the Supreme Court approved the exclusion of India's socially and economically backward classes from its affirmative action measures to address economic deprivation. Dissenting justices, including the Chief Justice of India, protested the Majority Opinion for sanctioning "an avowedly exclusionary and discriminatory principle." To justify their controversial decision, the majority justices relied on technical arguments that are categorically false. The confusion of the justices is due to a subtle technical aspect of the affirmative action system in India: the significance of overlaps between members of protected groups. Conventionally, protected classes were determined by the caste system, which meant they did not overlap. The addition of a new protected class defined by economic criteria alters this structure, unless it is artificially enforced. The majority justices failed to appreciate the significance of this critical change and inaccurately argued that the controversial exclusion is a technical necessity to provide benefits to previously unprotected members of a new class. We show that this case could have been resolved with three competing policies that each avoid the controversial exclusion. One of these policies aligns with the core arguments in the Majority Opinion, whereas a second aligns with those in the Dissenting Opinion.


title: Energy of boundary of spacetime

abstract: We consider how the energy can be stored in the boundary of spacetime, in particular in a spherical bubble that can be made by a quantum gravitational process. Our calculation is performed within the framework of classical Einstein gravity by identifying the Gibbons-Hawking-York term as the membrane action. We show that the energy of the bubble can be given consistently with the Schwarzschild metric. The solution of the consistency condition suggests positive membrane tension, which explains why we do not observe such topological defects in ordinary experiences and also gives a mechanism for suppressing the spacetime with boundary in a dynamical way.


title: Granger causal inference on DAGs identifies genomic loci regulating transcription

abstract: When a dynamical system can be modeled as a sequence of observations, Granger causality is a powerful approach for detecting predictive interactions between its variables. However, traditional Granger causal inference has limited utility in domains where the dynamics need to be represented as directed acyclic graphs (DAGs) rather than as a linear sequence, such as with cell differentiation trajectories. Here, we present GrID-Net, a framework based on graph neural networks with lagged message passing for Granger causal inference on DAG-structured systems. Our motivating application is the analysis of single-cell multimodal data to identify genomic loci that mediate the regulation of specific genes. To our knowledge, GrID-Net is the first single-cell analysis tool that accounts for the temporal lag between a genomic locus becoming accessible and its downstream effect on a target gene's expression. We applied GrID-Net on multimodal single-cell assays that profile chromatin accessibility (ATAC-seq) and gene expression (RNA-seq) in the same cell and show that it dramatically outperforms existing methods for inferring regulatory locus-gene links, achieving up to 71% greater agreement with independent population genetics-based estimates. By extending Granger causality to DAG-structured dynamical systems, our work unlocks new domains for causal analyses and, more specifically, opens a path towards elucidating gene regulatory interactions relevant to cellular differentiation and complex human diseases at unprecedented scale and resolution.


title: Expectations Formation with Fat-tailed Processes: Evidence from Sales Forecasts

abstract: We empirically analyze a large sample of firm sales growth expectations. We find that the relationship between forecast errors and lagged revision is non-linear. Forecasters underreact to typical (positive or negative) news about future sales, but overreact to very significant news. To account for this non-linearity, we propose a simple framework, where (1) sales growth dynamics have a fat-tailed high frequency component and (2) forecasters use a simple linear rule. This framework qualitatively fits several additional features of data on sales growth dynamics, forecast errors, and stock returns.


title: Second-order self-consistent field algorithms: from classical to quantum nuclei

abstract: This work presents a general framework for deriving exact and approximate Newton self-consistent field (SCF) orbital optimization algorithms by leveraging concepts borrowed from differential geometry. Within this framework, we extend the augmented Roothaan--Hall (ARH) algorithm to unrestricted electronic and nuclear-electronic calculations. We demonstrate that ARH yields an excellent compromise between stability and computational cost for SCF problems that are hard to converge with conventional first-order optimization strategies. In the electronic case, we show that ARH overcomes the slow convergence of orbitals in strongly-correlated molecules with the example of several iron-sulfur clusters. For nuclear-electronic calculations, ARH significantly enhances the convergence already for small molecules, as demonstrated for a series of protonated water clusters.


title: Doubly-robust and heteroscedasticity-aware sample trimming for causal inference

abstract: A popular method for variance reduction in observational causal inference is propensity-based trimming, the practice of removing units with extreme propensities from the sample. This practice has theoretical grounding when the data are homoscedastic and the propensity model is parametric (Yang and Ding, 2018; Crump et al. 2009), but in modern settings where heteroscedastic data are analyzed with non-parametric models, existing theory fails to support current practice. In this work, we address this challenge by developing new methods and theory for sample trimming. Our contributions are three-fold: first, we describe novel procedures for selecting which units to trim. Our procedures differ from previous work in that we trim not only units with small propensities, but also units with extreme conditional variances. Second, we give new theoretical guarantees for inference after trimming. In particular, we show how to perform inference on the trimmed subpopulation without requiring that our regressions converge at parametric rates. Instead, we make only fourth-root rate assumptions like those in the double machine learning literature. This result applies to conventional propensity-based trimming as well and thus may be of independent interest. Finally, we propose a bootstrap-based method for constructing simultaneously valid confidence intervals for multiple trimmed sub-populations, which are valuable for navigating the trade-off between sample size and variance reduction inherent in trimming. We validate our methods in simulation, on the 2007-2008 National Health and Nutrition Examination Survey, and on a semi-synthetic Medicare dataset and find promising results in all settings.


title: Simplex Range Searching Revisited: How to Shave Logs in Multi-Level Data Structures

abstract: We revisit the classic problem of simplex range searching and related problems in computational geometry. We present a collection of new results which improve previous bounds by multiple logarithmic factors that were caused by the use of multi-level data structures. Highlights include the following: $\bullet$ For a set of $n$ points in a constant dimension $d$, we give data structures with $O(n^d)$ (or slightly better) space that can answer simplex range counting queries in optimal $O(\log n)$ time and simplex range reporting queries in optimal $O(\log n + k)$ time, where $k$ denotes the output size. For semigroup range searching, we obtain $O(\log n)$ query time with $O(n^d\mathop{\rm polylog}n)$ space. Previous data structures with similar space bounds by Matou\v{s}ek from nearly three decades ago had $O(\log^{d+1}n)$ or $O(\log^{d+1}n + k)$ query time. $\bullet$ For a set of $n$ simplices in a constant dimension $d$, we give data structures with $O(n)$ space that can answer stabbing counting queries (counting the number of simplices containing a query point) in $O(n^{1-1/d})$ time, and stabbing reporting queries in $O(n^{1-1/d}+k)$ time. Previous data structures had extra $\log^d n$ factors in space and query time. $\bullet$ For a set of $n$ (possibly intersecting) line segments in 2D, we give a data structure with $O(n)$ space that can answer ray shooting queries in $O(\sqrt{n})$ time. This improves Wang's recent data structure [SoCG'20] with $O(n\log n)$ space and $O(\sqrt{n}\log n)$ query time.


title: Faster Matrix Multiplication via Asymmetric Hashing

abstract: Fast matrix multiplication is one of the most fundamental problems in algorithm research. The exponent of the optimal time complexity of matrix multiplication is usually denoted by $\omega$. This paper discusses new ideas for improving the laser method for fast matrix multiplication. We observe that the analysis of higher powers of the Coppersmith-Winograd tensor [Coppersmith & Winograd 1990] incurs a "combination loss", and we partially compensate for it using an asymmetric version of CW's hashing method. By analyzing the eighth power of the CW tensor, we give a new bound of $\omega<2.371866$, which improves the previous best bound of $\omega<2.372860$ [Alman & Vassilevska Williams 2020]. Our result breaks the lower bound of $2.3725$ in [Ambainis, Filmus & Le Gall 2015] because of the new method for analyzing component (constituent) tensors.


title: Bifurcation results for nonlinear eigenvalue problems involving the (p,q)-Laplace operator

abstract: In this paper, we analyze an eigenvalue problem for nonlinear elliptic operators involving homogeneous Dirichlet boundary conditions in a open smooth bounded domain. We prove bifurcation results from trivial solutions and from infinity for the considered nonlinear eigenvalue problem. We also show the existence of multiple solutions of the nonlinear problem using variational methods.


title: Intra-Source Style Augmentation for Improved Domain Generalization

abstract: The generalization with respect to domain shifts, as they frequently appear in applications such as autonomous driving, is one of the remaining big challenges for deep learning models. Therefore, we propose an intra-source style augmentation (ISSA) method to improve domain generalization in semantic segmentation. Our method is based on a novel masked noise encoder for StyleGAN2 inversion. The model learns to faithfully reconstruct the image preserving its semantic layout through noise prediction. Random masking of the estimated noise enables the style mixing capability of our model, i.e. it allows to alter the global appearance without affecting the semantic layout of an image. Using the proposed masked noise encoder to randomize style and content combinations in the training set, ISSA effectively increases the diversity of training data and reduces spurious correlation. As a result, we achieve up to $12.4\%$ mIoU improvements on driving-scene semantic segmentation under different types of data shifts, i.e., changing geographic locations, adverse weather conditions, and day to night. ISSA is model-agnostic and straightforwardly applicable with CNNs and Transformers. It is also complementary to other domain generalization techniques, e.g., it improves the recent state-of-the-art solution RobustNet by $3\%$ mIoU in Cityscapes to Dark Z\"urich.


title: Entity-Focused Dense Passage Retrieval for Outside-Knowledge Visual Question Answering

abstract: Most Outside-Knowledge Visual Question Answering (OK-VQA) systems employ a two-stage framework that first retrieves external knowledge given the visual question and then predicts the answer based on the retrieved content. However, the retrieved knowledge is often inadequate. Retrievals are frequently too general and fail to cover specific knowledge needed to answer the question. Also, the naturally available supervision (whether the passage contains the correct answer) is weak and does not guarantee question relevancy. To address these issues, we propose an Entity-Focused Retrieval (EnFoRe) model that provides stronger supervision during training and recognizes question-relevant entities to help retrieve more specific knowledge. Experiments show that our EnFoRe model achieves superior retrieval performance on OK-VQA, the currently largest outside-knowledge VQA dataset. We also combine the retrieved knowledge with state-of-the-art VQA models, and achieve a new state-of-the-art performance on OK-VQA.


title: Polynomial bounds on torsion from a fixed geometric isogeny class of elliptic curves

abstract: We show there exist polynomial bounds on torsion of elliptic curves which come from a fixed geometric isogeny class. More precisely, for an elliptic curve $E_0$ defined over a number field $F_0$, for each $\epsilon>0$ there exist constants $c_\epsilon:=c_\epsilon(E_0,F_0),C_\epsilon:=C_\epsilon(E_0,F_0)>0$ such that for any elliptic curve $E_{/F}$ geometrically isogenous to $E_0$, if $E(F)$ has a point of order $N$ then \[ N\leq c_\epsilon\cdot [F:\mathbb{Q}]^{1/2+\epsilon}, \] and one also has \[ \# E(F)[\textrm{tors}] \leq C_\epsilon\cdot [F:\mathbb{Q}]^{1+\epsilon}. \]


title: Phelps Property U and $C(K)$ spaces

abstract: A subspace $X$ of a Banach space $Y$ has $\textit{Property U}$ whenever every continuous linear functional on $X$ has a unique norm-preserving (i.e., Hahn$-$Banach) extension to $Y$ (Phelps, 1960). Throughout this document we introduce and develop a systematic study of the existence of $\textit{U-embeddings}$ between Banach spaces $X$ and $Y$, that is, isometric embeddings of $X$ into $Y$ whose ranges have property U. In particular, we are interested in the case that $Y=C(K)$, where $K$ is a compact Hausdorff topological space. We provide results for general Banach spaces and for some specific set-ups, such as $X$ being a finite-dimensional space or a $C(K)$-space.


title: Inference in conditioned dynamics through causality restoration

abstract: Computing observables from conditioned dynamics is typically computationally hard, because, although obtaining independent samples efficiently from the unconditioned dynamics is usually feasible, generally most of the samples must be discarded (in a form of importance sampling) because they do not satisfy the imposed conditions. Sampling directly from the conditioned distribution is non-trivial, as conditioning breaks the causal properties of the dynamics which ultimately renders the sampling procedure efficient. One standard way of achieving it is through a Metropolis Monte-Carlo procedure, but this procedure is normally slow and a very large number of Monte-Carlo steps is needed to obtain a small number of statistically independent samples. In this work, we propose an alternative method to produce independent samples from a conditioned distribution. The method learns the parameters of a generalized dynamical model that optimally describe the conditioned distribution in a variational sense. The outcome is an effective, unconditioned, dynamical model, from which one can trivially obtain independent samples, effectively restoring causality of the conditioned distribution. The consequences are twofold: on the one hand, it allows us to efficiently compute observables from the conditioned dynamics by simply averaging over independent samples. On the other hand, the method gives an effective unconditioned distribution which is easier to interpret. The method is flexible and can be applied virtually to any dynamics. We discuss an important application of the method, namely the problem of epidemic risk assessment from (imperfect) clinical tests, for a large family of time-continuous epidemic models endowed with a Gillespie-like sampler. We show that the method compares favorably against the state of the art, including the soft-margin approach and mean-field methods.

categories: physics.data-an cond-mat.stat-mech cs.LG q-bio.PE


title: Domain Adaptation in 3D Object Detection with Gradual Batch Alternation Training

abstract: We consider the problem of domain adaptation in LiDAR-based 3D object detection. Towards this, we propose a simple yet effective training strategy called Gradual Batch Alternation that can adapt from a large labeled source domain to an insufficiently labeled target domain. The idea is to initiate the training with the batch of samples from the source and target domain data in an alternate fashion, but then gradually reduce the amount of the source domain data over time as the training progresses. This way the model slowly shifts towards the target domain and eventually better adapt to it. The domain adaptation experiments for 3D object detection on four benchmark autonomous driving datasets, namely ONCE, PandaSet, Waymo, and nuScenes, demonstrate significant performance gains over prior arts and strong baselines.


title: Comparing Embedded Graphs Using Average Branching Distance

abstract: Graphs drawn in the plane are ubiquitous, arising from data sets through a variety of methods ranging from GIS analysis to image classification to shape analysis. A fundamental problem in this type of data is comparison: given a set of such graphs, can we rank how similar they are, in such a way that we capture their geometric "shape" in the plane? In this paper we explore a method to compare two such embedded graphs, via a simplified combinatorial representation called a tail-less merge tree which encodes the structure based on a fixed direction. First, we examine the properties of a distance designed to compare merge trees called the branching distance, and show that the distance as defined in previous work fails to satisfy some of the requirements of a metric. We incorporate this into a new distance function called average branching distance to compare graphs by looking at the branching distance for merge trees defined over many directions. Despite the theoretical issues, we show that the definition is still quite useful in practice by using our open-source code to cluster data sets of embedded graphs.


title: Landmark Enforcement and Style Manipulation for Generative Morphing

abstract: Morph images threaten Facial Recognition Systems (FRS) by presenting as multiple individuals, allowing an adversary to swap identities with another subject. Morph generation using generative adversarial networks (GANs) results in high-quality morphs unaffected by the spatial artifacts caused by landmark-based methods, but there is an apparent loss in identity with standard GAN-based morphing methods. In this paper, we propose a novel StyleGAN morph generation technique by introducing a landmark enforcement method to resolve this issue. Considering this method, we aim to enforce the landmarks of the morph image to represent the spatial average of the landmarks of the bona fide faces and subsequently the morph images to inherit the geometric identity of both bona fide faces. Exploration of the latent space of our model is conducted using Principal Component Analysis (PCA) to accentuate the effect of both the bona fide faces on the morphed latent representation and address the identity loss issue with latent domain averaging. Additionally, to improve high frequency reconstruction in the morphs, we study the train-ability of the noise input for the StyleGAN2 model.


title: Inflationary and phase-transitional primordial magnetic fields in galaxy clusters

abstract: Primordial magnetic fields (PMFs) are possible candidates for explaining the observed magnetic fields in galaxy clusters. Two competing scenarios of primordial magnetogenesis have been discussed in the literature: inflationary and phase-transitional. We study the amplification of both large- and small-scale correlated magnetic fields, corresponding to inflation- and phase transition-generated PMFs, in a massive galaxy cluster. We employ high-resolution magnetohydrodynamic cosmological zoom-in simulations to resolve the turbulent motions in the intracluster medium. We find that the turbulent amplification is more efficient for the large-scale inflationary models, while the phase transition-generated seed fields show moderate growth. The differences between the models are imprinted on the spectral characteristics of the field (such as the amplitude and the shape of the magnetic power spectrum) and therefore on the final correlation length. We find a one order of magnitude difference between the final strengths of the inflation- and phase transition-generated magnetic fields, and a factor of 1.5 difference between their final coherence scales. Thus, the final configuration of the magnetic field retains information about the PMF generation scenarios. Our findings have implications for future extragalactic Faraday rotation surveys with the possibility of distinguishing between different magnetogenesis scenarios.


title: Application Performance Modeling via Tensor Completion

abstract: Performance tuning, software/hardware co-design, and job scheduling are among the many tasks that rely on models to predict application performance. We propose and evaluate low-rank tensor decomposition for modeling application performance. We discretize the input and configuration domains of an application using regular grids. Application execution times mapped within grid-cells are averaged and represented by tensor elements. We show that low-rank canonical-polyadic (CP) tensor decomposition is effective in approximating these tensors. We further show that this decomposition enables accurate extrapolation of unobserved regions of an application's parameter space. We then employ tensor completion to optimize a CP decomposition given a sparse set of observed execution times. We consider alternative piecewise/grid-based models and supervised learning models for six applications and demonstrate that CP decomposition optimized using tensor completion offers higher prediction accuracy and memory-efficiency for high-dimensional performance modeling.


title: An Adaptive Hybrid Control Algorithm for Sender-Receiver Clock Synchronization

abstract: This paper presents an innovative hybrid systems approach to the sender-receiver synchronization of timers. Via the hybrid systems framework, we unite the traditional sender-receiver algorithm for clock synchronization with an online, adaptive strategy to achieve synchronization of the clock rates to exponentially synchronize a pair of clocks connected over a network. Following the conventions of the algorithm, clock measurements of the nodes are given at periodic time instants, and each node uses these measurements to achieve synchronization. For this purpose, we introduce a hybrid system model of a network with continuous and impulsive dynamics that captures the sender-receiver algorithm as a state-feedback controller to synchronize the network clocks. Moreover, we provide sufficient design conditions that ensure attractivity of the synchronization set.


title: Mermin polytopes in quantum computation and foundations

abstract: Mermin square scenario provides a simple proof for state-independent contextuality. In this paper, we study polytopes $\text{MP}_\beta$ obtained from the Mermin scenario, parametrized by a function $\beta$ on the set of contexts. Up to combinatorial isomorphism, there are two types of polytopes $\text{MP}_0$ and $\text{MP}_1$ depending on the parity of $\beta$. Our main result is the classification of the vertices of these two polytopes. In addition, we describe the graph associated with the polytopes. All the vertices of $\text{MP}_0$ turn out to be deterministic. This result provides a new topological proof of a celebrated result of Fine characterizing noncontextual distributions on the CHSH scenario. $\text{MP}_1$ can be seen as a nonlocal toy version of $\Lambda$-polytopes, a class of polytopes introduced for the simulation of universal quantum computation. In the $2$-qubit case, we provide a decomposition of the $\Lambda$-polytope using $\text{MP}_1$, whose vertices are classified, and the nonsignaling polytope of the $(2,3,2)$ Bell scenario, whose vertices are well-known.


title: Extreme mass loss in low-mass type Ib/c supernova progenitors

abstract: Many core collapse supernovae (SNe) with hydrogen-poor and low-mass ejecta, such as ultra-stripped SNe and type Ibn SNe, are observed to interact with dense circumstellar material (CSM). These events likely arise from the core-collapse of helium stars which have been heavily stripped by a binary companion and ejected significant mass during the last weeks to years of their lives. In helium star models run to days before core-collapse, we identify a range of helium core masses $\approx 2.5 -3 M_{\odot}$ whose envelopes expand substantially due to helium shell burning while the core undergoes neon and oxygen burning. When modeled in binary systems, the rapid expansion of these helium stars induces extremely high rates of late-stage mass transfer ($\dot{M} \gtrsim 10^{-2} M_\odot/{\rm yr}$) beginning weeks to decades before core-collapse. We consider two scenarios for producing CSM in these systems: either mass transfer remains stable and mass loss is driven from the system in the vicinity of the accreting companion, or mass transfer becomes unstable and causes a common envelope event (CEE) through which the helium envelope is unbound. The ensuing CSM properties are consistent with the CSM masses ($\sim 10^{-2}-1 M_\odot$) and radii ($\sim 10^{13}-10^{16} {\rm cm}$) inferred for ultra-stripped SNe and several type Ibn SNe. Furthermore, systems that undergo a CEE could produce short-period NS binaries that merge in less than 100 Myr.


title: On Hitting Times for General Quantum Markov Processes

abstract: Random walks (or Markov chains) are models extensively used in theoretical computer science. Several tools, including analysis of quantities such as hitting and mixing times, are helpful for devising randomized algorithms. A notable example is Sch\"oning's algorithm for the satisfiability (SAT) problem. In this work, we use the density-matrix formalism to define a quantum Markov chain model which directly generalizes classical walks, and we show that a common tools such as hitting times can be computed with a similar formula as the one found in the classical theory, which we then apply to known quantum settings such as Grover's algorithm.


title: Assessment of various Hamiltonian partitionings for the electronic structure problem on a quantum computer using the Trotter approximation

abstract: Solving the electronic structure problem via unitary evolution of the electronic Hamiltonian is one of the promising applications of digital quantum computers. One of the practical strategies to implement the unitary evolution is via Trotterization, where a sequence of short-time evolutions of fast-forwardable (i.e. efficiently diagonalizable) Hamiltonian fragments is used. Given multiple choices of possible Hamiltonian decompositions to fast-forwardable fragments, the accuracy of the Hamiltonian evolution depends on the choice of the fragments. We assess efficiency of multiple Hamiltonian partitioning techniques using fermionic and qubit algebras for the Trotterization. Use of symmetries of the electronic Hamiltonian and its fragments significantly reduces the Trotter error. This reduction makes fermionic-based partitioning Trotter errors lower compared to those in qubit-based techniques. However, from the simulation-cost standpoint, fermionic methods tend to introduce quantum circuits with a greater number of T-gates at each Trotter step and thus are more computationally expensive compared to their qubit counterparts.


title: Kindling the First Stars: I. Dependence of Detectability of the First Stars with JWST on the Pop III Stellar Masses

abstract: The first Pop III stars formed out of primordial, metal free gas, in minihalos at z>20, and kickstarted the cosmic processes of reionizaton and enrichment. While these stars are likely more massive than their enriched counterparts, the current unknowns of their astrophysics include; when the first Pop III stars ignited, how massive they were, and when and how the era of the first stars ended. Investigating these questions requires an exploration of a multi-dimensional parameter space, including the slope of the Pop III stellar initial mass function (IMF) and the strength of the non-ionizing UV background. In this work, we present a novel model which treats both the slope and maximum mass of Pop III stars as truly free parameters while including the physics of the fragmentation of primordial gas. Our results also hint at a non-universal Pop III IMF which is dependent on the efficiency of primordial gas fragmentation. Our relatively simple model reproduces the results from hydrodynamic simulations, but with a computational efficiency which allows us to investigate the observable differences between a wide range of potential Pop III IMFs. In addition, the evolution of the number density of Pop III stars may provide insight into the evolution of the H2 dissociating background. While the slope of the Pop III IMF does not significantly affect the predicted number density of the first stars, more top heavy IMFs produce Pop III star clusters which are 2-3 magnitudes brighter than their more bottom heavy counterparts. While the Pop III star clusters are too dim for direct detection by JWST, we find they are within the reach of gravitational lensing.


title: Simple and Effective Unsupervised Speech Translation

abstract: The amount of labeled data to train models for speech tasks is limited for most languages, however, the data scarcity is exacerbated for speech translation which requires labeled data covering two different languages. To address this issue, we study a simple and effective approach to build speech translation systems without labeled data by leveraging recent advances in unsupervised speech recognition, machine translation and speech synthesis, either in a pipeline approach, or to generate pseudo-labels for training end-to-end speech translation models. Furthermore, we present an unsupervised domain adaptation technique for pre-trained speech models which improves the performance of downstream unsupervised speech recognition, especially for low-resource settings. Experiments show that unsupervised speech-to-text translation outperforms the previous unsupervised state of the art by 3.2 BLEU on the Libri-Trans benchmark, on CoVoST 2, our best systems outperform the best supervised end-to-end models (without pre-training) from only two years ago by an average of 5.0 BLEU over five X-En directions. We also report competitive results on MuST-C and CVSS benchmarks.


title: Mixed Isogeometric Methods for Linear Elasticity with Weakly Imposed Symmetry

abstract: We consider and discretize a mixed formulation for linear elasticity with weakly imposed symmetry in two and three dimensions. Whereas existing methods mainly deal with simplicial or polygonal meshes, we take advantage of isogeometric analysis (IGA) and consequently allow for shapes with curved boundaries. To introduce the discrete spaces we use isogeometric discrete differential forms defined by proper B-spline spaces. For the proposed schemes a proof of well-posedness and an error estimate are given. Further we discuss our ansatz by means of different numerical examples.


title: $\lambda$-MIMO: Massive MIMO via Modulo Sampling

abstract: Massive multiple-input multiple-output (M-MIMO) architecture is the workhorse of modern communication systems. Currently, two fundamental bottlenecks, namely, power consumption and receiver saturation, limit the full potential achievement of this technology. These bottlenecks are intricately linked with the analog-to-digital converter (ADC) used in each radio frequency (RF) chain. The power consumption in MIMO systems grows exponentially with the ADC's bit budget while ADC saturation causes permanent loss of information. This motivates the need for a solution that can simultaneously tackle the above-mentioned bottlenecks while offering advantages over existing alternatives such as low-resolution ADCs. Taking a radically different approach to this problem, we propose $\lambda$-MIMO architecture which uses modulo ADCs ($M_\lambda$-ADC) instead of a conventional ADC. Our work is inspired by the Unlimited Sampling Framework. $M_\lambda$-ADC in the RF chain folds high dynamic range signals into low dynamic range modulo samples, thus alleviating the ADC saturation problem. At the same time, digitization of modulo signal results in high resolution quantization. In the novel $\lambda$-MIMO context, we discuss baseband signal reconstruction, detection and uplink achievable sum-rate performance. The key takeaways of our work include, (a) leveraging higher signal-to-quantization noise ratio (SQNR), (b) detection and average uplink sum-rate performances comparable to a conventional, infinite-resolution ADC when using a $1$-$2$ bit $M_\lambda$-ADC. This enables higher order modulation schemes e.g. $1024$ QAM that seemed previously impossible, (c) superior trade-off between energy efficiency and bit budget, thus resulting in higher power efficiency. Numerical simulations and modulo ADC based hardware experiments corroborate our theory and reinforce the clear benefits of $\lambda$-MIMO approach.


title: Rethinking Prototypical Contrastive Learning through Alignment, Uniformity and Correlation

abstract: Contrastive self-supervised learning (CSL) with a prototypical regularization has been introduced in learning meaningful representations for downstream tasks that require strong semantic information. However, to optimize CSL with a loss that performs the prototypical regularization aggressively, e.g., the ProtoNCE loss, might cause the "coagulation" of examples in the embedding space. That is, the intra-prototype diversity of samples collapses to trivial solutions for their prototype being well-separated from others. Motivated by previous works, we propose to mitigate this phenomenon by learning Prototypical representation through Alignment, Uniformity and Correlation (PAUC). Specifically, the ordinary ProtoNCE loss is revised with: (1) an alignment loss that pulls embeddings from positive prototypes together; (2) a uniformity loss that distributes the prototypical level features uniformly; (3) a correlation loss that increases the diversity and discriminability between prototypical level features. We conduct extensive experiments on various benchmarks where the results demonstrate the effectiveness of our method in improving the quality of prototypical contrastive representations. Particularly, in the classification down-stream tasks with linear probes, our proposed method outperforms the state-of-the-art instance-wise and prototypical contrastive learning methods on the ImageNet-100 dataset by 2.96% and the ImageNet-1K dataset by 2.46% under the same settings of batch size and epochs.


title: Curriculum Reinforcement Learning using Optimal Transport via Gradual Domain Adaptation

abstract: Curriculum Reinforcement Learning (CRL) aims to create a sequence of tasks, starting from easy ones and gradually learning towards difficult tasks. In this work, we focus on the idea of framing CRL as interpolations between a source (auxiliary) and a target task distribution. Although existing studies have shown the great potential of this idea, it remains unclear how to formally quantify and generate the movement between task distributions. Inspired by the insights from gradual domain adaptation in semi-supervised learning, we create a natural curriculum by breaking down the potentially large task distributional shift in CRL into smaller shifts. We propose GRADIENT, which formulates CRL as an optimal transport problem with a tailored distance metric between tasks. Specifically, we generate a sequence of task distributions as a geodesic interpolation (i.e., Wasserstein barycenter) between the source and target distributions. Different from many existing methods, our algorithm considers a task-dependent contextual distance metric and is capable of handling nonparametric distributions in both continuous and discrete context settings. In addition, we theoretically show that GRADIENT enables smooth transfer between subsequent stages in the curriculum under certain conditions. We conduct extensive experiments in locomotion and manipulation tasks and show that our proposed GRADIENT achieves higher performance than baselines in terms of learning efficiency and asymptotic performance.


title: BirdSoundsDenoising: Deep Visual Audio Denoising for Bird Sounds

abstract: Audio denoising has been explored for decades using both traditional and deep learning-based methods. However, these methods are still limited to either manually added artificial noise or lower denoised audio quality. To overcome these challenges, we collect a large-scale natural noise bird sound dataset. We are the first to transfer the audio denoising problem into an image segmentation problem and propose a deep visual audio denoising (DVAD) model. With a total of 14,120 audio images, we develop an audio ImageMask tool and propose to use a few-shot generalization strategy to label these images. Extensive experimental results demonstrate that the proposed model achieves state-of-the-art performance. We also show that our method can be easily generalized to speech denoising, audio separation, audio enhancement, and noise estimation.


title: 1/f noise and anomalous scaling in L\'evy noise-driven on-off intermittency

abstract: On-off intermittency occurs in nonequilibrium physical systems close to bifurcation points and is characterised by an aperiodic switching between a large-amplitude "on" state and a small-amplitude "off" state. L\'evy on-off intermittency is a recently introduced generalisation of on-off intermittency to multiplicative L\'evy noise, which depends on a stability parameter $\alpha$ and a skewness parameter $\beta$. Here, we derive two novel results on L\'evy on-off intermittency by leveraging known exact results on the first-passage time statistics of L\'evy flights. First, we compute anomalous critical exponents explicitly as a function of arbitrary L\'evy noise parameters $(\alpha,\beta)$ for the first time, by a heuristic method, complementing previous results. The predictions are verified using numerical solutions of the fractional Fokker-Planck equation. Second, we derive the power spectrum $S(f)$ of L\'evy on-off intermittency and show that it displays a power law $S(f)\propto f^\kappa$ at low frequencies $f$, where $\kappa\in (-1,0)$ depends on the noise parameters $\alpha,\beta$. An explicit expression for $\kappa$ is obtained in terms of $(\alpha,\beta)$. The predictions are verified using long time series realisations of L\'evy on-off intermittency. Our findings help shed light on instabilities subject to non-equilibrium, power-law-distributed fluctuations, emphasizing that their properties can differ starkly from the case of Gaussian fluctuations.


title: Radiation Transport Two-Temperature GRMHD Simulations of Warped Accretion Disks

abstract: In many black hole systems, the accretion disk is expected to be misaligned with respect to the black hole spin axis. If the scale height of the disk is much smaller than the misalignment angle, the spin of the black hole can tear the disk into multiple, independently precessing `sub-disks'. This is most likely to happen during outbursts in black hole X-Ray binaries (BHXRBs) and in active galactic nuclei (AGN) accreting above a few percent of the Eddington limit, because the disk becomes razor-thin. Disk tearing has the potential to explain variability phenomena including quasi-periodic oscillations (QPOs) in BHXRBs and changing-look phenomena in AGN. Here, we present the first radiative two-temperature GRMHD simulation of a strongly tilted ($65^{\circ}$) accretion disk around a $M_{BH}=10M_{\odot}$ black hole, which tears and precesses. This leads to luminosity swings between a few percent and $50 \%$ of the Eddington limit on sub-viscous timescales. Surprisingly, even where the disk is radiation pressure dominated, the accretion disk is thermally stable over $t \gtrsim 21,000 r_g/c$. This suggests warps play an important role in stabilizing the disk against thermal collapse. The disk forms two nozzle shocks perpendicular to the line of nodes where the scale height of the disk decreases $10$-fold and the electron temperature reaches $T_e \sim 10^8-10^9 K$. In addition, optically thin gas crossing the tear between the inner and outer disk gets heated to $T_e \sim 10^8 K$. This suggests that warped disks may emit a Comptonized spectrum that deviates substantially from idealized models.


title: Bayesian Optimization over Discrete and Mixed Spaces via Probabilistic Reparameterization

abstract: Optimizing expensive-to-evaluate black-box functions of discrete (and potentially continuous) design parameters is a ubiquitous problem in scientific and engineering applications. Bayesian optimization (BO) is a popular, sample-efficient method that leverages a probabilistic surrogate model and an acquisition function (AF) to select promising designs to evaluate. However, maximizing the AF over mixed or high-cardinality discrete search spaces is challenging standard gradient-based methods cannot be used directly or evaluating the AF at every point in the search space would be computationally prohibitive. To address this issue, we propose using probabilistic reparameterization (PR). Instead of directly optimizing the AF over the search space containing discrete parameters, we instead maximize the expectation of the AF over a probability distribution defined by continuous parameters. We prove that under suitable reparameterizations, the BO policy that maximizes the probabilistic objective is the same as that which maximizes the AF, and therefore, PR enjoys the same regret bounds as the original BO policy using the underlying AF. Moreover, our approach provably converges to a stationary point of the probabilistic objective under gradient ascent using scalable, unbiased estimators of both the probabilistic objective and its gradient. Therefore, as the number of starting points and gradient steps increase, our approach will recover of a maximizer of the AF (an often-neglected requisite for commonly used BO regret bounds). We validate our approach empirically and demonstrate state-of-the-art optimization performance on a wide range of real-world applications. PR is complementary to (and benefits) recent work and naturally generalizes to settings with multiple objectives and black-box constraints.


title: Helpful Neighbors: Leveraging Neighbors in Geographic Feature Pronunciation

abstract: If one sees the place name Houston Mercer Dog Run in New York, how does one know how to pronounce it? Assuming one knows that Houston in New York is pronounced "how-ston" and not like the Texas city, then one can probably guess that "how-ston" is also used in the name of the dog park. We present a novel architecture that learns to use the pronunciations of neighboring names in order to guess the pronunciation of a given target feature. Applied to Japanese place names, we demonstrate the utility of the model to finding and proposing corrections for errors in Google Maps. To demonstrate the utility of this approach to structurally similar problems, we also report on an application to a totally different task: Cognate reflex prediction in comparative historical linguistics. A version of the code has been open-sourced (https://github.com/google-research/google-research/tree/master/cognate_inpaint_neighbors).