title: Revisiting Weighted Aggregation in Federated Learning with Neural Networks

abstract: In federated learning (FL), weighted aggregation of local models is conducted to generate a global model, and the aggregation weights are normalized (the sum of weights is 1) and proportional to the local data sizes. In this paper, we revisit the weighted aggregation process and gain new insights into the training dynamics of FL. First, we find that the sum of weights can be smaller than 1, causing global weight shrinking effect (analogous to weight decay) and improving generalization. We explore how the optimal shrinking factor is affected by clients' data heterogeneity and local epochs. Second, we dive into the relative aggregation weights among clients to depict the clients' importance. We develop client coherence to study the learning dynamics and find a critical point that exists. Before entering the critical point, more coherent clients play more essential roles in generalization. Based on the above insights, we propose an effective method for Federated Learning with Learnable Aggregation Weights, named as FedLAW. Extensive experiments verify that our method can improve the generalization of the global model by a large margin on different datasets and models.


title: Balanced Audiovisual Dataset for Imbalance Analysis

abstract: The imbalance problem is widespread in the field of machine learning, which also exists in multimodal learning areas caused by the intrinsic discrepancy between modalities of samples. Recent works have attempted to solve the modality imbalance problem from algorithm perspective, however, they do not fully analyze the influence of modality bias in datasets. Concretely, existing multimodal datasets are usually collected under specific tasks, where one modality tends to perform better than other ones in most conditions. In this work, to comprehensively explore the influence of modality bias, we first split existing datasets into different subsets by estimating sample-wise modality discrepancy. We surprisingly find that: the multimodal models with existing imbalance algorithms consistently perform worse than the unimodal one on specific subsets, in accordance with the modality bias. To further explore the influence of modality bias and analyze the effectiveness of existing imbalance algorithms, we build a balanced audiovisual dataset, with uniformly distributed modality discrepancy over the whole dataset. We then conduct extensive experiments to re-evaluate existing imbalance algorithms and draw some interesting findings: existing algorithms only provide a compromise between modalities and suffer from the large modality discrepancy of samples. We hope that these findings could facilitate future research on the modality imbalance problem.


title: From paintbrush to pixel: A review of deep neural networks in AI-generated art

abstract: This paper delves into the fascinating field of AI-generated art and explores the various deep neural network architectures and models that have been utilized to create it. From the classic convolutional networks to the cutting-edge diffusion models, we examine the key players in the field. We explain the general structures and working principles of these neural networks. Then, we showcase examples of milestones, starting with the dreamy landscapes of DeepDream and moving on to the most recent developments, including Stable Diffusion and DALL-E 2, which produce mesmerizing images. A detailed comparison of these models is provided, highlighting their strengths and limitations. Thus, we examine the remarkable progress that deep neural networks have made so far in a short period of time. With a unique blend of technical explanations and insights into the current state of AI-generated art, this paper exemplifies how art and computer science interact.


title: GLUECons: A Generic Benchmark for Learning Under Constraints

abstract: Recent research has shown that integrating domain knowledge into deep learning architectures is effective -- it helps reduce the amount of required data, improves the accuracy of the models' decisions, and improves the interpretability of models. However, the research community is missing a convened benchmark for systematically evaluating knowledge integration methods. In this work, we create a benchmark that is a collection of nine tasks in the domains of natural language processing and computer vision. In all cases, we model external knowledge as constraints, specify the sources of the constraints for each task, and implement various models that use these constraints. We report the results of these models using a new set of extended evaluation criteria in addition to the task performances for a more in-depth analysis. This effort provides a framework for a more comprehensive and systematic comparison of constraint integration techniques and for identifying related research challenges. It will facilitate further research for alleviating some problems of state-of-the-art neural models.


title: Conformers are All You Need for Visual Speech Recognition

abstract: Visual speech recognition models extract visual features in a hierarchical manner. At the lower level, there is a visual front-end with a limited temporal receptive field that processes the raw pixels depicting the lips or faces. At the higher level, there is an encoder that attends to the embeddings produced by the front-end over a large temporal receptive field. Previous work has focused on improving the visual front-end of the model to extract more useful features for speech recognition. Surprisingly, our work shows that complex visual front-ends are not necessary. Instead of allocating resources to a sophisticated visual front-end, we find that a linear visual front-end paired with a larger Conformer encoder results in lower latency, more efficient memory usage, and improved WER performance. We achieve a new state-of-the-art of 12.8% WER for visual speech recognition on the TED LRS3 dataset, which rivals the performance of audio-only models from just four years ago.


title: How Generative AI models such as ChatGPT can be (Mis)Used in SPC Practice, Education, and Research? An Exploratory Study

abstract: Generative Artificial Intelligence (AI) models such as OpenAI's ChatGPT have the potential to revolutionize Statistical Process Control (SPC) practice, learning, and research. However, these tools are in the early stages of development and can be easily misused or misunderstood. In this paper, we give an overview of the development of Generative AI. Specifically, we explore ChatGPT's ability to provide code, explain basic concepts, and create knowledge related to SPC practice, learning, and research. By investigating responses to structured prompts, we highlight the benefits and limitations of the results. Our study indicates that the current version of ChatGPT performs well for structured tasks, such as translating code from one language to another and explaining well-known concepts but struggles with more nuanced tasks, such as explaining less widely known terms and creating code from scratch. We find that using new AI tools may help practitioners, educators, and researchers to be more efficient and productive. However, in their current stages of development, some results are misleading and wrong. Overall, the use of generative AI models in SPC must be properly validated and used in conjunction with other methods to ensure accurate results.


title: A matrix-free macro-element variant of the hybridized discontinuous Galerkin method

abstract: We investigate a macro-element variant of the hybridized discontinuous Galerkin (HDG) method, using patches of standard simplicial elements that can have non-matching interfaces. Coupled via the HDG technique, our method enables local refinement by uniform simplicial subdivision of each macro-element. By enforcing one spatial discretization for all macro-elements, we arrive at local problems per macro-element that are embarrassingly parallel, yet well balanced. Therefore, our macro-element variant scales efficiently to n-node clusters and can be tailored to available hardware by adjusting the local problem size to the capacity of a single node, while still using moderate polynomial orders such as quadratics or cubics. Increasing the local problem size means simultaneously decreasing, in relative terms, the global problem size, hence effectively limiting the proliferation of degrees of freedom. The global problem is solved via a matrix-free iterative technique that also heavily relies on macro-element local operations. We investigate and discuss the advantages and limitations of the macro-element HDG method via an advection-diffusion model problem.


title: Level set-based multiscale topology optimization for a thermal cloak design problem using the homogenization method

abstract: Artificially designed composite materials consist of microstructures, that exhibit various thermal properties depending on their shapes, such as anisotropic thermal conductivity. One of the representative applications of such composite materials for thermal control is the thermal cloak. This study proposed a topology optimization method based on a level set method for a heat conduction problem to optimally design composite materials that achieve a thermal cloak. The homogenization method was introduced to evaluate its effective thermal conductivity coefficient. Then, we formulated a multiscale topology optimization method for the composite materials in the framework of the homogenization method, where the microstructures were optimized to minimize objective functions defined using the macroscopic temperature field. We presented examples of optimal structures in a two-dimensional problem and discussed the validity of the obtained structures.


title: Multiverse Predictions for Habitability: Element Abundances

abstract: We investigate the dependence of elemental abundances on physical constants, and the implications this has for the distribution of complex life for various proposed habitability criteria. We consider three main sources of abundance variation: differing supernova rates, alpha burning in massive stars, and isotopic stability, and how each affects the metal-to-rock ratio and the abundances of carbon, oxygen, nitrogen, phosphorus, sulfur, silicon, magnesium, and iron. Our analysis leads to several predictions for which habitability criteria are correct by determining which ones make our observations of the physical constants, as well as a few other observed features of our universe, most likely. Our results indicate that carbon-rich or carbon-poor planets are uninhabitable, slightly magnesium-rich planets are habitable, and life does not depend on nitrogen abundance too sensitively. We also find suggestive but inconclusive evidence that metal-rich planets and phosphorus-poor planets are habitable. These predictions can then be checked by probing regions of our universe that closely resemble normal environments in other universes. If any of these predictions are found to be wrong, the multiverse scenario would predict that the majority of observers are born in universes differing substantially from ours, and so can be ruled out, to varying degrees of statistical significance.


title: 'The Taurus': Cattle Breeds & Diseases Identification Mobile Application using Machine Learning

abstract: Dairy farming plays an important role in agriculture for thousands of years not only in Sri Lanka but also in so many other countries. When it comes to dairy farming cattle is an indispensable animal. According to the literature surveys almost 3.9 million cattle and calves die in a year due to different types of diseases. The causes of diseases are mainly bacteria, parasites, fungi, chemical poisons and etc. Infectious diseases can be a greatest threat to livestock health. The mortality rate of cattle causes a huge impact on social, economic and environmental damage. In order to decrease this negative impact, the proposal implements a cross-platform mobile application to easily analyze and identify the diseases which cattle suffer from and give them a solution and also to identify the cattle breeds. The mobile application is designed to identify the breeds by analyzing the images of the cattle and identify diseases after analyzing the videos and the images of affected areas. Then make a model to identify the weight and the age of a particular cow and suggest the best dose of the medicine to the identified disease. This will be a huge advantage to farmers as well as to dairy industry. The name of the proposed mobile application is 'The Taurus' and this paper address the selected machine learning and image processing models and the approaches taken to identify the diseases, breeds and suggest the prevention methods and medicine to the identified disease.


title: Effect of rigid top surface on the onset of phototactic bioconvection

abstract: Phototaxis is the movement of microorganisms towards or away from light. In presence of dim (intense) light, microorganisms move towards (away from) the light sources, which is known as positive (negative) phototaxis. In this study, we present a model in which isotropic scattering algae suspension is illuminated by both diffuse and collimated irradiation, and the top surface is assumed to be rigid. Bimodal steady states result from the scattering Effect (for almost purely scattering suspension), which transits into a unimodal steady state as diffuse light intensity increases. The linear stability of the same suspension is also examined using the linear perturbation theory, which predicts the solutions' stable and oscillatory nature for specific parameter values. Here, we also observe that suspension becomes more stable in a rigid top surface and diffuse irradiation.


title: Gas Kinetic Schemes for Solving the Magnetohydrodynamic Equations with Pressure Anisotropy

abstract: In many astrophysical plasmas, the Coulomb collision is insufficient to maintain an isotropic temperature, and the system is driven to the anisotropic regime. In this case, magnetohydrodynamic (MHD) models with anisotropic pressure are needed to describe such a plasma system. To solve the anisotropic MHD equation numerically, we develop a robust Gas-Kinetic flux scheme for non-linear MHD flows. Using anisotropic velocity distribution functions, the numerical flux functions are derived for updating the macroscopic plasma variables. The schemes is suitable for finite-volume solvers which utilize a conservative form of the mass, momentum and total energy equations, and can be easily applied to multi-fluid problems and extended to more generalized double polytropic plasma systems. Test results show that the numerical scheme is very robust and performs well for both linear wave and non-linear MHD problems.


title: PointFISH -- learning point cloud representations for RNA localization patterns

abstract: Subcellular RNA localization is a critical mechanism for the spatial control of gene expression. Its mechanism and precise functional role is not yet very well understood. Single Molecule Fluorescence in Situ Hybridization (smFISH) images allow for the detection of individual RNA molecules with subcellular accuracy. In return, smFISH requires robust methods to quantify and classify RNA spatial distribution. Here, we present PointFISH, a novel computational approach for the recognition of RNA localization patterns. PointFISH is an attention-based network for computing continuous vector representations of RNA point clouds. Trained on simulations only, it can directly process extracted coordinates from experimental smFISH images. The resulting embedding allows scalable and flexible spatial transcriptomics analysis and matches performance of hand-crafted pipelines.


title: A Reinforcement Learning Framework for Online Speaker Diarization

abstract: Speaker diarization is a task to label an audio or video recording with the identity of the speaker at each given time stamp. In this work, we propose a novel machine learning framework to conduct real-time multi-speaker diarization and recognition without prior registration and pretraining in a fully online and reinforcement learning setting. Our framework combines embedding extraction, clustering, and resegmentation into the same problem as an online decision-making problem. We discuss practical considerations and advanced techniques such as the offline reinforcement learning, semi-supervision, and domain adaptation to address the challenges of limited training data and out-of-distribution environments. Our approach considers speaker diarization as a fully online learning problem of the speaker recognition task, where the agent receives no pretraining from any training set before deployment, and learns to detect speaker identity on the fly through reward feedbacks. The paradigm of the reinforcement learning approach to speaker diarization presents an adaptive, lightweight, and generalizable system that is useful for multi-user teleconferences, where many people might come and go without extensive pre-registration ahead of time. Lastly, we provide a desktop application that uses our proposed approach as a proof of concept. To the best of our knowledge, this is the first approach to apply a reinforcement learning approach to the speaker diarization task.


title: Effect of Geomagnetic Storms on a Power Network at Mid Latitudes

abstract: Solar activities may disturb the geomagnetic field and impact the power grid via geomagnetically induced currents. We study active and reactive powers as well as the power factor of Iran's power grid transformers (230 kV and 400 kV) and their correlations with geomagnetic disturbances indices (SYM-H < -30 nT and sizable horizontal geomagnetic field fluctuation) from 19 March 2018 to 20 March 2020. Out of 128,627 cases with a transformer power factor of less than 0.7, we observe that 12,112 samples correlated with SYM-H. Our investigation shows that about 4 percent of two years, the SYM-H has values less than -30 nT. Analysis of high-performance transformers (a power factor greater than 0.7 at 95 percent of working time) shows at least a 55 percent correlation of power factor less than 0.7 and SYM-H less than -30 nT. We observe that the transformers' power factor of Rafsanjan-Kerman and Sarcheshmeh-Kerman (wye connection on 230 kV side) substations decreased to less than 0.7 and correlated with SYM-H. We show that the reactive power of the Sefidrood-Guilan and Shahid Beheshti-Guilan transformers (wye configurations) increased considerably on 9 January 2020 and positively correlated with SYM-H may produce large GICs at this part of the grid. We observe that the increase in reactive power at the Shahid Beheshti-Guilan substation correlated with the sizable changes in the horizontal field recorded by the Jaipur station. However, more details (temperature and current of transformers) need records to estimate the impact of geomagnetic induction current on the power grid.


title: Edge coalitions in graphs

abstract: An edge coalition in a graph $G=(V,E)$ consists of two disjoint sets of edges $E_1$ and $E_2$, neither of which is an edge dominating set but whose union $E_1\cup E_2$ is an edge dominating set. An edge coalition partition in a graph $G$ of order $n=|V|$ and size $m$ is an edge partition $\pi=\{E_1,\cdots,E_k\}$ so that every set $E_i$ of $\pi$ either is a singleton edge dominating set, or is not an edge dominating set but forms an edge coalition with another set $E_j$ which is not an edge dominating set. In this paper we introduce the concept of edge coalition and show that there exists edge coalition for some graphs and trees. The graphs $G$ with small and size number of edge coalition are characterized. Finally, coalition graphs of special graphs are studied.


title: Spatial gradient consistency for unsupervised learning of hyperspectral demosaicking: Application to surgical imaging

abstract: Hyperspectral imaging has the potential to improve intraoperative decision making if tissue characterisation is performed in real-time and with high-resolution. Hyperspectral snapshot mosaic sensors offer a promising approach due to their fast acquisition speed and compact size. However, a demosaicking algorithm is required to fully recover the spatial and spectral information of the snapshot images. Most state-of-the-art demosaicking algorithms require ground-truth training data with paired snapshot and high-resolution hyperspectral images, but such imagery pairs with the exact same scene are physically impossible to acquire in intraoperative settings. In this work, we present a fully unsupervised hyperspectral image demosaicking algorithm which only requires exemplar snapshot images for training purposes. We regard hyperspectral demosaicking as an ill-posed linear inverse problem which we solve using a deep neural network. We take advantage of the spectral correlation occurring in natural scenes to design a novel inter spectral band regularisation term based on spatial gradient consistency. By combining our proposed term with standard regularisation techniques and exploiting a standard data fidelity term, we obtain an unsupervised loss function for training deep neural networks, which allows us to achieve real-time hyperspectral image demosaicking. Quantitative results on hyperspetral image datasets show that our unsupervised demosaicking approach can achieve similar performance to its supervised counter-part, and significantly outperform linear demosaicking. A qualitative user study on real snapshot hyperspectral surgical images confirms the results from the quantitative analysis. Our results suggest that the proposed unsupervised algorithm can achieve promising hyperspectral demosaicking in real-time thus advancing the suitability of the modality for intraoperative use.


title: End-to-end kilonova models of neutron-star mergers with delayed black-hole formation

abstract: We investigate the nucleosynthesis and kilonova properties of binary neutron-star (NS) merger models which lead to intermediate remnant lifetimes of ~0.1-1seconds until black-hole (BH) formation and describe all components of material ejected during the dynamical merger phase, NS-remnant evolution, and final viscous disintegration of the BH torus after gravitational collapse. To this end we employ a combination of hydrodynamics, nucleosynthesis, and radiative-transfer tools to achieve a consistent end-to-end modeling of the system and its observables. We adopt a novel version of the Shakura-Sunyaev scheme allowing to vary the approximate turbulent viscosity inside the NS remnant independently of the surrounding disk. We find that asymmetric progenitors lead to shorter remnant lifetimes and enhanced ejecta masses, although the viscosity affects the absolute values of these characteristics. The integrated production of lanthanides and heavier elements in such binary systems is sub-solar, suggesting that the considered scenarios contribute in a sub-dominant fashion to r-process enrichment. One reason is that BH-tori formed after delayed collapse exhibit less neutron-rich conditions than typically found, and often assumed in previous BH-torus models, for early BH formation. The outflows in our models feature strong anisotropy as a result of the lanthanide-poor polar neutrino-driven wind pushing aside lanthanide-rich dynamical ejecta. Considering the complexity of the models, the estimated kilonova light curves show promising agreement with AT2017gfo after times of several days, while the remaining inconsistencies at early times could possibly be overcome in binary configurations with a more dominant neutrino-driven wind relative to the dynamical ejecta.


title: The Evryscope Fast Transient Engine: Real-Time Detection for Rapidly Evolving Transients

abstract: Astrophysical transients with rapid development on sub-hour timescales are intrinsically rare. Due to their short durations, events like stellar superflares, optical flashes from gamma-ray bursts, and shock breakouts from young supernovae are difficult to identify on timescales that enable spectroscopic followup. This paper presents the Evryscope Fast Transient Engine (EFTE), a new data reduction pipeline designed to provide low-latency transient alerts from the Evryscopes, a North-South pair of ultra-wide-field telescopes with an instantaneous footprint covering 38% of the entire sky, and tools for building long-term light curves from Evryscope data. EFTE leverages the optical stability of the Evryscopes by using a simple direct image subtraction routine suited to continuously monitoring the transient sky at minute cadence. Candidates are produced within the base Evryscope two-minute cadence for 98.5% of images, and internally filtered using VetNet, a convolutional neural network real-bogus classifier. EFTE provides an extensible, robust architecture for transient surveys probing similar timescales, and serves as the software testbed for the real-time analysis pipelines and public data distribution systems for the Argus Array, a next generation all-sky observatory with a data rate 62x higher than Evryscope.


title: Active galactic nuclei with high-resolution X-ray spectroscopy

abstract: The imminent launch of XRISM will usher in an era of high-resolution X-ray spectroscopy. For active galactic nuclei (AGN) this is an exciting epoch that is full of massive potential for uncovering the ins and outs of supermassive black hole accretion. In this work, we review AGN research topics that are certain to advance in the coming years with XRISM and prognosticate the possibilities with Athena and Arcus. Specifically, our discussion focuses on: (i) the relatively slow moving ionised winds known as warm absorbers and obscurers; (ii) the iron emitting from different regions of the inner and outer disc, broad line region, and torus; and (iii) the ultrafast outflows that may be the key to understanding AGN feedback.


title: The Concentration-Mass Relation of Massive, Dynamically Relaxed Galaxy Clusters: Agreement Between Observations and $\Lambda$CDM Simulations

abstract: The relationship linking a galaxy cluster's total mass with the concentration of its mass profile and its redshift is a fundamental prediction of the Cold Dark Matter (CDM) paradigm of cosmic structure formation. However, confronting those predictions with observations is complicated by the fact that simulated clusters are not representative of observed samples where detailed mass profile constraints are possible. In this work, we calculate the Symmetry-Peakiness-Alignment (SPA) morphology metrics for maps of X-ray emissivity from THE THREE HUNDRED project hydrodynamical simulations of galaxy clusters at four redshifts, and thereby select a sample of morphologically relaxed, simulated clusters, using observational criteria. These clusters have on average earlier formation times than the full sample, confirming that they are both morphologically and dynamically more relaxed than typical. We constrain the concentration-mass-redshift relation of both the relaxed and complete sample of simulated clusters, assuming power-law dependences on mass ($\kappa_m$) and $1+z$ ($\kappa_\zeta$), finding $\kappa_m = -0.12 \pm 0.07$ and $\kappa_\zeta = -0.27 \pm 0.19$ for the relaxed subsample. From an equivalently selected sample of massive, relaxed clusters observed with ${\it Chandra}$, we find $\kappa_m = -0.12 \pm 0.08$ and $\kappa_\zeta = -0.48 \pm 0.19$, in good agreement with the simulation predictions. The simulated and observed samples also agree well on the average concentration at a pivot mass and redshift providing further validation of the $\Lambda$CDM paradigm in the properties of the largest gravitationally collapsed structures observed. This also represents the first clear detection of decreasing concentration with redshift, a longstanding prediction of simulations, in data.


title: Scary Barbie: An Extremely Energetic, Long-Duration Tidal Disruption Event Candidate Without a Detected Host Galaxy at z = 0.995

abstract: We report multi-wavelength observations and characterization of the ultraluminous transient AT 2021lwx (ZTF20abrbeie; aka ``Barbie'') identified in the alert stream of the Zwicky Transient Facility (ZTF) using a Recommender Engine For Intelligent Transient Tracking (REFITT) filter on the ANTARES alert broker. From a spectroscopically measured redshift of 0.995, we estimate a peak observed pseudo-bolometric luminosity of log (L$_{\text{max}} / [\text{erg}/\text{s}]$) = 45.7 from slowly fading ztf-$\it{g}$ and ztf-$r$ light curves spanning over 1000 observer-frame days. The host galaxy is not detected in archival Pan-STARRS observations ($g > 23.3$ mag), implying a lower limit to the outburst amplitude of more than 5 mag relative to the quiescent host galaxy. Optical spectra from Lick and Keck Observatories exhibit strong emission lines with narrow cores from the H Balmer series and ultraviolet semi-forbidden lines of Si III] $\lambda$1892, C III] $\lambda$1909, and C II] $\lambda$2325. Typical nebular lines in AGN spectra from ions such as [O II] and [O III] are not detected. These spectral features, along with the smooth light curve that is unlike most AGN flaring activity, and the luminosity that exceeds any observed or theorized supernova, lead us to conclude that AT 2021lwx is most likely an extreme tidal disruption event (TDE). Modeling of ZTF photometry with MOSFiT suggests that the TDE was between a $\approx 14 M_{\odot}$ star and a supermassive black hole of mass $M_{\text{BH}} \sim$ $10^{8} M_{\odot}$. Continued monitoring of the still-evolving light curve along with deep imaging of the field once AT 2021lwx has faded can test this hypothesis and potentially detect the host galaxy.


title: Quench dynamics of the Schwinger model via variational quantum algorithms

abstract: We investigate the real-time dynamics of the $(1+1)$-dimensional U(1) gauge theory known as the Schwinger model via variational quantum algorithms. Specifically, we simulate quench dynamics in the presence of an external electric field. First, we use a variational quantum eigensolver to obtain the ground state of the system in the absence of an external field. With this as the initial state, we perform real-time evolution under an external field via a fixed-depth, parameterized circuit whose parameters are updated using McLachlan's variational principle. We use the same Ansatz for initial state preparation and time evolution, by which we are able to reduce the overall circuit depth. We test our method with a classical simulator and confirm that the results agree well with exact diagonalization.

categories: hep-ph cond-mat.str-el hep-lat hep-th quant-ph


title: COSMOS2020: Exploring the dawn of quenching for massive galaxies at 3 < z < 5 with a new colour selection method

abstract: We select and characterise a sample of massive (log(M$_{*}/$M$_{\odot})>10.6$) quiescent galaxies (QGs) at $3<z<5$ in the latest COSMOS2020 catalogue. QGs are selected using a new rest-frame colour selection method, based on their probability of belonging to the quiescent group defined by a Gaussian Mixture Model (GMM) trained on rest-frame colours ($NUV-U, U-V, V-J$) of similarly massive galaxies at $2<z<3$. We calculate the quiescent probability threshold above which a galaxy is classified as quiescent using simulated galaxies from the SHARK semi-analytical model. We find that at $z\geq3$ in SHARK, the GMM/$NUVU-VJ$ method out-performs classical rest-frame $UVJ$ selection and is a viable alternative. We select galaxies as quiescent based on their probability in COSMOS2020 at $3<z<5$, and compare the selected sample to both $UVJ$ and $NUVrJ$ selected samples. We find that although the new selection matches $UVJ$ and $NUVrJ$ in number, the overlap between colour selections is only $\sim50-80\%$, implying that rest-frame colour commonly used at lower redshifts selections cannot be equivalently used at $z>3$. We compute median rest-frame SEDs for our sample and find the median quiescent galaxy at $3<z<5$ has a strong Balmer/4000 Angstrom break, and residual $NUV$ flux indicating recent quenching. We find the number densities of the entire quiescent population (including post-starbursts) more than doubles from $3.5\pm2.2\times10^{-6}$ Mpc$^{-3}$ at $4<z<5$ to $1.4\pm0.4\times10^{-5}$ Mpc$^{-3}$ at $3<z<4$, confirming that the onset of massive galaxy quenching occurs as early as $3<z<5$.


title: Complex path simulations of geometrically frustrated ladders

abstract: Quantum systems with geometrical frustration remain an outstanding challenge for numerical simulations due to the infamous numerical sign problem. Here, we overcome this obstruction via complex path integration in a geometrically frustrated ladder of interacting bosons at finite density. This enables studies of the many-body ground state properties, otherwise inaccessible with standard quantum Monte Carlo methods. Specifically, we study a chemical potential tuned quantum phase transition, along which we track the emergence of quasi long range order and critical softening of the single particle gap. We chart future methodological improvements and applications in generalized geometrically frustrated lattice models.


title: An Atlas of Color-selected Quiescent Galaxies at $z>3$ in Public $JWST$ Fields

abstract: We present the results of a systematic search for candidate quiescent galaxies in the distant Universe in eleven $JWST$ fields with publicly available observations collected during the first three months of operations and covering an effective sky area of $\sim145$ arcmin$^2$. We homogeneously reduce the new $JWST$ data and combine them with existing observations from the $Hubble\,Space\,Telescope$. We select a robust sample of $\sim80$ candidate quiescent and quenching galaxies at $3 < z < 5$ using two methods: (1) based on their rest-frame $UVJ$ colors, and (2) a novel quantitative approach based on Gaussian Mixture Modeling of the $NUV-U$, $U-V$, and $V-J$ rest-frame color space, which is more sensitive to recently quenched objects. We measure comoving number densities of massive ($M_\star\geq 10^{10.6} M_\odot$) quiescent galaxies consistent with previous estimates relying on ground-based observations, after homogenizing the results in the literature with our mass and redshift intervals. However, we find significant field-to-field variations of the number densities up to a factor of $2-3$, highlighting the effect of cosmic variance and suggesting the presence of overdensities of red quiescent galaxies at $z>3$, as it could be expected for highly clustered massive systems. Importantly, $JWST$ enables the robust identification of quenching/quiescent galaxy candidates at lower masses and higher redshifts than before, challenging standard formation scenarios. All data products, including the literature compilation, are made publicly available.


title: The coevolution of supermassive black holes and galaxies in luminous AGN over a wide range of redshift

abstract: It is well known that supermassive black holes (SMBHs) and their host galaxies co-evolve. A manifestation of this co-evolution is the correlation that has been found between the SMBH mass, M$_{BH}$, and the galaxy bulge or stellar mass, M$_*$. The cosmic evolution of this relation, though, is still a matter of debate. In this work, we examine the M$_{BH}-$M$_*$ relation, using 687 X-ray luminous (median $\rm log\,[L_{X,2-10keV}(ergs^{-1})]=44.3$), broad line AGN, at $\rm 0.2<z<4.0$ (median $\rm z\approx 1.4$) that lie in the XMM-{\it{XXL}} field. Their M$_{BH}$ and M$_*$ range from $\rm 7.5<log\,[M_{BH}\,(M_\odot)]<9.5$ and $\rm 10<log\,[M_*(M_\odot)]<12$, respectively. Most of the AGN live in star-forming galaxies and their Eddington ratios range from 0.01 to 1, with a median value of 0.06. Our results show that M$_{BH}$ and M$_*$ are correlated ($\rm r=0.47\pm0.21$, averaged over different redshift intervals). Our analysis also shows that the mean ratio of the M$_{BH}$ and M$_*$ does not evolve with redshift, at least up to $\rm z=2$ and has a value of $\rm log($M$_{BH}/$M$_*)=-2.44$. The majority of the AGN ($75\%$) are in a SMBH mass growth dominant phase. In these systems, the M$_{BH}-$M$_*$ correlation is weaker and their M$_*$ tends to be lower (for the same M$_{BH}$) compared to systems that are in a galaxy mass growth phase. Our findings suggest that the growth of black hole mass occurs first, while the early stellar mass assembly may not be so efficient.


title: Comment on: "Possible Relation between the Cosmological Constant and Standard Model Parameters"

abstract: It has recently been proposed [M. P. Hertzberg and A. Loeb, arXiv:2302.09090] that the observed value of the cosmological constant can be related to the physical parameters of the Standard Model. In the present compact note we point out that the derivation of the claimed cosmological-constant-standard-model relation presented in \cite{HL} is, unfortunately, erroneous.


title: Clumpy wind studies and the non-detection of cyclotron line in OAO 1657-415

abstract: Winds of massive stars are suspected to be inhomogeneous (or clumpy), which biases the measures of their mass loss rates. In High Mass X-ray Binaries (HMXBs), the compact object can be used as an orbiting X-ray point source to probe the wind and constrain its clumpiness. We perform spectro-timing analysis of the HMXB OAO 1657-415 with non-simultaneous NuSTAR and NICER observations. We compute the hardness ratio from the energy-resolved light curves, and using an adaptive rebinning technique, we thus select appropriate time segments to search for rapid spectral variations on timescales of a few hundreds to thousands of seconds. Column density and intensity of Iron K$\alpha$ line were strongly correlated, and the recorded spectral variations were consistent with accretion from a clumpy wind. We also illustrate a novel framework to measure clump sizes, masses in HMXBs more accurately based on absorption measurements and orbital parameters of the source. We then discuss the limitations posed by current X-ray spacecrafts in such measurements and present prospects with future X-ray missions. We find that the source pulse profiles show a moderate dependence on energy. We identify a previously undetected dip in the pulse profile visible throughout the NuSTAR observation near spin phase 0.15 possibly caused by intrinsic changes in accretion geometry close to the neutron star. We do not find any evidence for the debated cyclotron line at $\sim$ 36\,keV in the time-averaged or the phase-resolved spectra with NuSTAR.


title: Dynamical mean-field theory for R\'{e}nyi entanglement entropy and mutual Information in Hubbard Model

abstract: Quantum entanglement, lacking any classical counterpart, provides a fundamental new route to characterize the quantum nature of many-body states. In this work, we discuss an implementation of a new path integral method [Phys. Rev. Res. 2, 033505 (2020)] for fermions to compute entanglement for extended subsystems in the Hubbard model within dynamical mean field theory (DMFT) in one and two dimensions. The new path integral formulation measures entanglement by applying a ``kick" to the underlying interacting fermions. We show that the R\'{e}nyi entanglement entropy can be extracted efficiently within the DMFT framework by integrating over the strength of the kick term. Using this method, we compute the second R\'{e}nyi entropy as a function of subsystem size for metallic and Mott insulating phases of the Hubbard model. We explore the thermal entropy to entanglement crossover in the subsystem R\'{e}nyi entropy in the correlated metallic phase. We show that the subsystem-size scaling of second R\'{e}nyi entropy is well described by the crossover formula which interpolates between the volume-law thermal R\'{e}nyi entropy and the universal boundary-law R\'{e}nyi entanglement entropy with logarithmic violation, as predicted by conformal field theory. We also study the mutual information across the Mott metal-insulator transition.

categories: cond-mat.str-el cond-mat.dis-nn cond-mat.stat-mech hep-th quant-ph


title: A lensed radio jet at milli-arcsecond resolution II: Constraints on fuzzy dark matter from an extended gravitational arc

abstract: Using a single gravitational lens system observed at $\lesssim5$ milli-arcsecond resolution with very long baseline interferometry (VLBI), we place a lower bound on the mass of the fuzzy dark matter (FDM) particle, ruling out $m_\chi \leq 4.4\times10^{-21}~\mathrm{eV}$ with a 20:1 posterior odds ratio relative to a smooth lens model. We generalize our result to non-scalar and multiple-field models, such as vector FDM, with $m_{\chi,\mathrm{vec}} > 1.4 \times 10^{-21}~\mathrm{eV}$. Due to the extended source structure and high angular resolution of the observation, our analysis is directly sensitive to the presence of granule structures in the main dark matter halo of the lens, which is the most generic prediction of FDM theories. A model based on well-understood physics of ultra-light dark matter fields in a gravitational potential well makes our result robust to a wide range of assumed dark matter fractions and velocity dispersions in the lens galaxy. Our result is competitive with other lower bounds on $m_\chi$ from past analyses, which rely on intermediate modelling of structure formation and/or baryonic effects. Higher resolution observations taken at 10 to 100 GHz could improve our constraints by up to 2 orders of magnitude in the future.


title: Black Hole Perturbations and Electric-Magnetic Duality

abstract: Black holes can be electromagnetically charged, or carry vector charge from new fundamental fields. Their response to small fluctuations is of paramount importance to study gravitational wave generation. However, the usual even and odd sectors of gravitoelectromagnetic waves couple if the black hole is magnetically charged, a fact that complicates significantly the perturbative approach. In this paper, perturbation theory based on harmonic expansion is extended to have manifest invariance under electric-magnetic duality. As a result, the equations decouple into two generalised even and odd sectors, each governed by master wave equations that include the most general coupling to a dyonic source. These can be used to compute, in a simple manner, the gravitational and electromagnetic radiation emitted in the interaction of the most general spherically symmetric black holes of the Einstein--Maxwell theory with electromagnetically charged matter.


title: The growth of brightest cluster galaxies in the TNG300 simulation:dissecting the contributions from mergers and in situ star formation

abstract: We investigate the formation of brightest cluster galaxies (BCGs) in the TNG300 cosmological simulation of the IllustrisTNG project. Our cluster sample consists of 700 haloes with $M_{200} \geq 5 \times 10^{13} \, \mathrm{M}_{\odot}$ at $z=0$, along with their progenitors at earlier epochs. This includes 280 systems with $M_{200} \geq 10^{14} \, \mathrm{M}_{\odot}$ at $z=0$, as well as three haloes with $M_{200} \geq 10^{15} \, \mathrm{M}_{\odot}$. We find that the stellar masses and star formation rates of our simulated BCGs are in good agreement with observations at $z \lesssim 0.4$, and that they have experienced, on average, $\sim$2 ($\sim$3) major mergers since $z=1$ ($z=2$). Separating the BCG from the intracluster light (ICL) by means of a fixed 30 kpc aperture, we find that the fraction of stellar mass contributed by ex situ (i.e. accreted) stars at $z=0$ is approximately 70, 80, and 90 per cent for the BCG, BCG+ICL, and ICL, respectively. Tracking our simulated BCGs back in time using the merger trees, we find that they became dominated by ex situ stars at $z \sim $1-2, and that half of the stars that are part of the BCG at $z=0$ formed early ($z \sim 3$) in other galaxies, but `assembled' onto the BCG until later times ($z \approx 0.8$ for the whole sample, $z \approx 0.5$ for BCGs in $M_{200} \geq 5 \times 10^{14} \, \mathrm{M}_{\odot}$ haloes). Finally, we show that the stellar mass profiles of BCGs are often dominated by ex situ stars at all radii, with stars from major mergers being found closer to the centre, while stars that were tidally stripped from other galaxies dominate the outer regions.


title: Renormalisable Non-local Quark-Gluon Interaction: Mass Gap, Chiral Symmetry Breaking & Scale Invariance

abstract: We derive a Nambu--Jona-Lasinio (NJL) model from a non-local gauge theory and show that it has confining properties at low energies. In particular, we present an extended approach to non-local QCD and a complete revision of the technique of Bender, Milton and Savage applied to non-local theories, providing a set of Dyson--Schwinger equations in differential form. In the local case, we obtain closed form solutions in the simplest case of the scalar field and extend it to the Yang--Mills field. In general, for non-local theories, we use a perturbative technique and a Fourier series and show how higher-order harmonics are heavily damped due to the presence of the non-local factor. The spectrum of the theory is analysed for the non-local Yang--Mills sector and found to be in agreement with the local results on the lattice in the limit of the non-locality mass parameter running to infinity. In the non-local case, we confine ourselves to a non-locality mass that is sufficiently large compared to the mass scale arising from the integration of the Dyson--Schwinger equations. Such a choice results in good agreement, in the proper limit, with the spectrum of the local theory. We derive the gap equation for the fermions in the theory that gives some indication of quark confinement in the non-local NJL case as well. Confinement seems to be a rather ubiquitous effect that removes some degrees of freedom in the original action, favouring the appearance of new observable states, as seen, e.g.,\ for quantum chromodynamics at lower energies.


title: Effect of spacetime geometry on neutrino oscillations

abstract: Propagation of fermions in spacetime requires a spin connection, which can be split into a universal gravitational part and a non-universal ``contorsion'' part. The latter is non-dynamical and can be eliminated from the theory, leaving an effective four-fermion interaction with unknown coupling constants. The generic form of the contorsion-fermion coupling, and thus the four-fermion interaction, breaks chiral symmetry. This interaction affects all fermions -- in particular neutrinos passing through matter will notice a contribution to their effective Hamiltonian, just like the MSW effect coming from weak interactions, but diagonal in the mass basis rather than the flavor basis. Then there is a possibility that this geometrical contribution is not negligible for neutrinos passing through normal matter, provided the coupling constants are not too small. We calculate the matter potential due to this interaction and thus write the conversion and survival probabilities including its effect. We plot conversion probabilities of $\nu_\mu$ to $\nu_e$ and $\nu_\tau$ for a baseline of 1300 km, as well as their dependence on the CP phase, with and without the geometrical interaction. We also plot the survival probability of reactor $\bar{\nu}_e$ for different baselines.


title: Enhanced C$_2$H$_2$ absorption within Jupiter's southern auroral oval from Juno UVS observations

abstract: Reflected sunlight observations from the Ultraviolet Spectrograph (UVS) on the Juno spacecraft were used to study the distribution of acetylene (C$_2$H$_2$) at Jupiter's south pole. We find that the shape of the C$_2$H$_2$ absorption feature varies significantly across the polar region, and this can be used to infer spatial variability in the C$_2$H$_2$ abundance. There is a localized region of enhanced C$_2$H$_2$ absorption which coincides with the location of Jupiter's southern polar aurora; the C$_2$H$_2$ abundance poleward of the auroral oval is a factor of 3 higher than adjacent quiescent, non-auroral longitudes. This builds on previous infrared studies which found enhanced C$_2$H$_2$ abundances within the northern auroral oval. This suggests that Jupiter's upper-atmosphere chemistry is being strongly influenced by the influx of charged auroral particles and demonstrates the necessity of developing ion-neutral photochemical models of Jupiter's polar regions.


title: Using helium 10830 {\AA} transits to constrain planetary magnetic fields

abstract: Planetary magnetic fields can affect the predicted mass loss rate for close-in planets that experience large amounts of UV irradiation. In this work, we present a method to detect the magnetic fields of close-in exoplanets undergoing atmospheric escape using transit spectroscopy at the 10830 Angstrom line of helium. Motivated by previous work on hydrodynamic and magneto-hydrodynamic photoevaporation, we suggest that planets with magnetic fields that are too weak to control the outflow's topology lead to blue-shifted transits due to day-to-night-side flows. In contrast, strong magnetic fields prevent this day-to-night flow, as the gas is forced to follow the magnetic field's roughly dipolar topology. We post-process existing 2D photoevaporation simulations to test this concept, computing synthetic transit profiles in helium. As expected, we find that hydrodynamically dominated outflows lead to blue-shifted transits on the order of the sound speed of the gas. Strong surface magnetic fields lead to unshifted or slightly red-shifted transit profiles. High-resolution observations can distinguish between these profiles; however, eccentricity uncertainties generally mean that we cannot conclusively say velocity shifts are due to the outflow for individual planets. The majority of helium observations are blue-shifted, which could be a tentative indication that close-in planets generally have surface dipole magnetic field strengths $\lesssim 0.1$ gauss. More 3D hydrodynamic and magneto-hydrodynamic are needed to confirm this conclusion robustly.


title: The pride of lions around Messier 105

abstract: We undertook a search for new dwarf galaxies in the Leo-I group using the data from the DECaLS digital sky survey. Five new presumed members of this group have been found in a wide vicinity of ${\rm M}\,105 ({\rm NGC}\,3379$). Currently, the group has a population of $83$ galaxies, $33$ of which have measured radial velocities. More than half of the group members belong to early types with no signs of ongoing star formation. About a quarter of the galaxies are outside the group's virial radius, $R_v = 385$~kpc. The presence of multiple systems with a size of about 15~kpc is evident in the group, but there are no noticeable global flat or filamentary substructures. The luminosity function of the group looks to be deficient in galaxies with absolute magnitudes in the interval $M_B = [-18, -15]$ mag. The ${\rm M}\,105$ group is characterized by a radial velocity dispersion of $136$~km~s$^{-1}$, orbital mass estimate $(5.76\pm1.32)\times 10^{12}~M_{\odot}$, and the total mass-to-K-band-luminosity ratio $(17.8\pm4.1) M_\odot/L_\odot$. The neighboring group of galaxies around ${\rm M}\,66 ({\rm NGC}\,3627$) has a similar virial radius, $390$~kpc, velocity dispersion, $135$~km~s$^{-1}$, and total mass-to-luminosity ratio, $(15.6\pm3.9) M_\odot/L_\odot$. Both groups in the Leo constellation are approaching the Local Group with a velocity of about 100~km~s$^{-1}$. In the background of the ${\rm M}\,105$ group, we noted a group of 6 galaxies with an unusually low virial mass-to-luminosity ratio, $M_T/L_K = (4.1\pm2.2) M_\odot/L_\odot$.


title: Block-encoding structured matrices for data input in quantum computing

abstract: The cost of data input can dominate the run-time of quantum algorithms. Here, we consider data input of arithmetically structured matrices via block encoding circuits, the input model for the quantum singular value transform and related algorithms. We demonstrate how to construct block encoding circuits based on an arithmetic description of the sparsity and pattern of repeated values of a matrix. We present schemes yielding different subnormalisations of the block encoding; a comparison shows that the best choice depends on the specific matrix. The resulting circuits reduce flag qubit number according to sparsity, and data loading cost according to repeated values, leading to an exponential improvement for certain matrices. We give examples of applying our block encoding schemes to a few families of matrices, including Toeplitz and tridiagonal matrices.


title: Test beam measurement of ATLAS ITk Short Strip module at warm and cold operational temperature

abstract: This study is focused on an investigation of the performance of the Short Strip module developed by the ATLAS Inner Tracker (ITk) strip collaboration using electron beams of energy 5.4 GeV and 5.8 GeV at the DESY-II Testbeam Facility. The noise at +30 C and -30 C was measured. The ratio of the two measurements is compared with a circuit-model calculation. The measured noise at -30 C is compared with the maximum noise that would correspond to an acceptable signal-to-noise ratio after the expected radiation damage from operation at HL-LHC. The measured charge distributions at +30 C and -30 C are compared with GEANT4 simulations. The detection efficiency and noise-occupancy were measured as a function of threshold at both +30 C and -30 C. The average cluster width was measured as a function of threshold. Scans of detection efficiency versus threshold at different delay settings were used to reconstruct the pulse shape in time. The resulting pulse shape was compared with a circuit model calculation.


title: Twin Theories, Polytope Mutations and Quivers for GTPs

abstract: We propose a unified perspective on two sets of objects that usually arise in the study of bipartite field theories. Each of the sets consists of a polytope, or equivalently a toric Calabi-Yau, and a quiver theory. We refer to the two sets of objects as original and twin. In the simplest cases, the two sides of the correspondence are connected by the graph operation known as untwisting. The democratic treatment that we advocate raises new questions regarding the connections between these objects, some of which we explore. With this motivation in mind, we establish a correspondence between the mutations of the original polytope and the twin quiver. This leads us to propose that non-toric twin quivers are naturally associated to generalized toric polygons (GTPs) and we explore various aspects of this idea. Supporting evidence includes global symmetries, the ability of twin quivers to encode the generalized $s$-rule, and the connection between the mutations of polytopes and of configurations of webs of 5-branes suspended from 7-branes. We introduce three methods for constructing twin quivers for GTPs. We also investigate the connection between twin quivers obtained using different toric phases. Twin quivers provide a powerful new perspective on GTPs. The ideas presented in this paper may represent a step towards the generalization of brane tilings to GTPs.


title: Machine learning for the prediction of safe and biologically active organophosphorus molecules

abstract: Drug discovery is a complex process with a large molecular space to be considered. By constraining the search space, the fragment-based drug design is an approach that can effectively sample the chemical space of interest. Here we propose a framework of Recurrent Neural Networks (RNN) with an attention model to sample the chemical space of organophosphorus molecules using the fragment-based approach. The framework is trained with a ZINC dataset that is screened for high druglikeness scores. The goal is to predict molecules with similar biological action modes as organophosphorus pesticides or chemical warfare agents yet less toxic to humans. The generated molecules contain a starting fragment of PO2F but have a bulky hydrocarbon side chain limiting its binding effectiveness to the targeted protein.


title: A NICER viewing angle on the accretion stream of Vela X-1

abstract: Vela X-1 is the archetypical eclipsing high-mass X-ray binary, composed of a neutron star (NS) accreting the B-star wind. It was observed by nearly all X-ray observatories, often multiple times, featuring a rich spectrum of variable emission lines. Yet, the precise origin of these lines in the binary system remains uncertain. We perform a systematic, orbital phase-dependent analysis of the reflected Fe K$\alpha$ fluorescence line at 6.4 keV using over 100 NICER observations. We resolve the line variability into 500s time bins and find that it is predominantly due to variation in the ionizing flux, with a moderate underlying phase dependence over the 9-day orbital period. Our analysis reveals a significant reflection component that cannot originate from the companion B-star alone. We also find that an appreciable portion of the B-star surface is obscured opposite the eclipse, and this obscuration is not symmetric around the mid-point (phase=0.5). We argue that an accretion stream, from the B-star to the NS and distorted by the orbital motion, is responsible for both the additional fluorescence emission component and for obscuring the B-star.


title: Polymer Complexation: Partially Ionizable Asymmetric Polyelectrolytes

abstract: Studies of the thermodynamics of complex coacervation of pairs of symmetric, strongly ionizable, oppositely charged polyelectrolyte chains are abundant. To generalize such understanding to asymmetric chain lengths and variable ionizability (chemical charge density), frequently observed in experiments, we present a theoretical framework to analyze the effective charge and size of the complex and the thermodynamics of complexation of two polyions as a function of such asymmetries. The free energy ensuing from the Edwards' Hamiltonian undergoes variational extremization, and explicitly accounts for the screened Coulomb and non-electrostatic interactions among monomers within individual polyions and between two polyions. Assuming maximal ion-pair formation of the complexed part, the system free energy comprising configurational entropy of the polyions and free-ion entropy of the small ions is minimized. The thermodynamic drive for complexation is found to increase with the ionizability of the symmetric polyions and to be maximum for symmetric chain lengths for equally ionizable polyions. The effective charge and size of the complex increase with asymmetry in charge density, where the size can be substantially larger than a collapsed globule found for symmetric chains. The regimes of enthalpy- and entropy-driven complexation are found, respectively, for low and high Coulomb strengths. The crossover strength is found to be strongly dependent on the dielectric environment and salt, but marginally dependent on the charge density, thus implying an entropy-driven process at moderate strengths. The key results match the trends in simulations and experiments, and are expected to provide insight for asymmetric complexation in real systems.


title: Duality and Fluxes in the Sen Formulation of Self-Dual Fields

abstract: In Sen's formulation of self-dual fields one finds two closed forms: $H^{(g)}$ and $H^{(s)}$. Only the former couples to sources and the spacetime metric. The latter has the wrong sign kinetic term but decouples and hence might be regarded as an unphysical artifact. In this letter we illustrate how an electromagnetic duality associated to the potential for $H^{(s)}$ gives rise to a T-like duality in the partition function for $H^{(g)}$. We then compute the partition function on a (4k+2)-dimensional torus highlighting its dependence on the choice of flux of $H^{(s)}$. Lastly we compute the two-point function of Wilson Surface operators.


title: A Search for OH 18-cm Emission from Intermediate-Velocity Gas at High Galactic Latitudes

abstract: We present search results of 22 high latitude (b > 25 deg.) sightlines for OH 18-cm emission using the 305-m radio telescope at the Arecibo Observatory. These sightlines appear in neutral hydrogen emission at intermediate velocities (V_lsr values ranging from -90 to -20 km/s) and are predicted to have a sufficient molecular composition so as to be detectable in molecular emission. Such objects, known as Intermediate-Velocity Molecular Clouds (IVMCs), have historically been detected through 12CO emission. Recent studies indicate that IVMCs may be widespread in the Galaxy and have important implications for models of the interstellar medium and star formation. However, we report non-detections of OH emission toward the 22 sightlines and provide stringent upper limits on the OH column density. Using available HI and Av data in combination with existing state-of-the-art PDR models, we estimate H2 column densities and find that they are more than an order of magnitude lower than the predicted values. We also find that the hydrogen volume density of these clouds is less than roughly 25 per cubic centimeter. In addition, we discuss the known IVMCs with previous 12CO detections in the context of the PDR models. Our analysis of these clouds indicates that the structure of molecular material in IVMCs is morphologically clumpy. These results motivate the need for future sensitive, on-the-fly searches (rather than targeted searches) for CO emission from IVMCs with of order roughly 1' resolution. High angular resolution (1') HI and Av data will also be helpful to better constrain the structure and composition of IVMCs.


title: The Near-Infrared Spectrograph (NIRSpec) on the James Webb Space Telescope V. Optimal algorithms for planning multi-object spectroscopic observations

abstract: We present an overview of the capabilities and key algorithms employed in the so-called eMPT software suite developed for planning scientifically optimized, multi-object spectroscopic (MOS) observations with the Micro-Shutter Array (MSA) of the Near-Infrared Spectrograph (NIRSpec) instrument on board the James Webb Space Telescope (JWST), the first multi-object spectrograph to operate in space. NIRSpec MOS mode is enabled by a programmable MSA, a regular grid of ~250,000 individual apertures that projects to a static, semi-regular pattern of available slits on the sky and makes the planning and optimization of an MSA observation a rather complex task. As such, the eMPT package is offered to the NIRSpec user community as a supplement to the MSA Planning Tool (MPT) included in the STScI Astronomer's Proposal Tool (APT) to assist in the planning of NIRSpec MOS proposals requiring advanced functionality to meet ambitious science goals. The eMPT produces output that can readily be imported and incorporated into the user's observing program within the APT to generate a customized MPT MOS observation. Furthermore, its novel algorithms and modular approach make it highly flexible and customizable, providing users the option to finely control the workflow and even insert their own software modules to tune their MSA slit masks to the particular scientific objectives at hand.


title: Adiabatic magnon spectra with and without constraining field: Benchmark against an exact magnon spectrum

abstract: The spectrum of magnon excitations in magnetic materials can be obtained exactly from the transverse dynamic magnetic susceptibility, which is however in practice numerically expensive. Many ab initio approaches therefore consider instead the adiabatic magnon spectrum, which assumes a separation of time scales of magnons and electronic excitations. There exist two alternative implementations for adiabatic magnon spectra: one based on the magnetic force theorem (MFT) and the other with a constraining field that enforces static non-collinear spin configurations. We benchmark both implementations against the exact magnon spectrum of an exactly solvable mean-field model. While both adiabatic methods are equally valid in the low magnon energy and strong Stoner coupling limits, we find that the constraining field method performs better than the MFT in both the cases of strong Stoner coupling and high magnon energies,while the MFT performs better for combined weak coupling and low magnon energies.


title: Dealing with Collinearity in Large-Scale Linear System Identification Using Gaussian Regression

abstract: Many problems arising in control require the determination of a mathematical model of the application. This has often to be performed starting from input-output data, leading to a task known as system identification in the engineering literature. One emerging topic in this field is estimation of networks consisting of several interconnected dynamic systems. We consider the linear setting assuming that system outputs are the result of many correlated inputs, hence making system identification severely ill-conditioned. This is a scenario often encountered when modeling complex cybernetics systems composed by many sub-units with feedback and algebraic loops. We develop a strategy cast in a Bayesian regularization framework where any impulse response is seen as realization of a zero-mean Gaussian process. Any covariance is defined by the so called stable spline kernel which includes information on smooth exponential decay. We design a novel Markov chain Monte Carlo scheme able to reconstruct the impulse responses posterior by efficiently dealing with collinearity. Our scheme relies on a variation of the Gibbs sampling technique: beyond considering blocks forming a partition of the parameter space, some other (overlapping) blocks are also updated on the basis of the level of collinearity of the system inputs. Theoretical properties of the algorithm are studied obtaining its convergence rate. Numerical experiments are included using systems containing hundreds of impulse responses and highly correlated inputs.


title: Simulating the LOcal Web (SLOW): I. Anomalies in the local density field

abstract: Context: Several observations of the local Universe (LU) point towards the existence of very prominent structures. The presence of massive galaxy clusters and local super clusters on the one hand, but also large local voids and under-densities on the other hand. However, it is highly non trivial to connect such different observational selected tracers to the underlying dark matter (DM) distribution. Methods (abridged): We used a 500 Mpc/h large constrained simulation of the LU with initial conditions based on peculiar velocities derived from the CosmicFlows-2 catalogue and follow galaxy formation physics directly in the hydro-dynamical simulations to base the comparison on stellar masses of galaxies or X-ray luminosity of clusters. We also used the 2668 Mpc/h large cosmological box from the Magneticum simulations to evaluate the frequency of finding such anomalies in random patches within simulations. Results: We demonstrate that haloes and galaxies in our constrained simulation trace the local DM density field very differently. Thereby, this simulation reproduces the observed 50% under-density of galaxy clusters and groups within the sphere of ~100 Mpc when applying the same mass or X-ray luminosity limit used in the observed cluster sample (CLASSIX), which is consistent with a ~1.5$\sigma$ feature. At the same time, the simulation reproduces the observed over-density of massive galaxy clusters within the same sphere, which on its own also corresponds to a ~1.5$\sigma$ feature. Interestingly, we find that only 44 out of 15635 random realizations (i.e. 0.28%) are matching both anomalies, making the LU to be a ~3$\sigma$ environment. We finally compared a mock galaxy catalogue with the observed distribution of galaxies in the LU, finding also a match to the observed factor of two over-density at ~16 Mpc as well as the observed 15% under-density at ~40 Mpc distance.


title: Imaging of exocomets with infrared interferometry

abstract: Active comets have been detected in several exoplanetary systems, although so far only indirectly, when the dust or gas in the extended coma has transited in front of the stellar disk. The large optical surface and relatively high temperature of an active cometary coma also makes it suitable to study with direct imaging, but the angular separation is generally too small to be reachable with present-day facilities. However, future imaging facilities with the ability to detect terrestrial planets in the habitable zones of nearby systems will also be sensitive to exocomets in such systems. Here we examine several aspects of exocomet imaging, particularly in the context of the Large Interferometer for Exoplanets (LIFE), which is a proposed space mission for infrared imaging and spectroscopy through nulling interferometry. We study what capabilities LIFE would have for acquiring imaging and spectroscopy of exocomets, based on simulations of the LIFE performance as well as statistical properties of exocomets that have recently been deduced from transit surveys. We find that for systems with extreme cometary activities such as beta Pictoris, sufficiently bright comets may be so abundant that they overcrowd the LIFE inner field of view. More nearby and moderately active systems such as epsilon Eridani or Fomalhaut may turn out to be optimal targets. If the exocomets have strong silicate emission features, such as in comet Hale-Bopp, it may become possible to study the mineralogy of individual exocometary bodies. We also discuss the possibility of exocomets as false positives for planets, with recent deep imaging of alpha Centauri as one hypothetical example. Such contaminants could be common, primarily among young debris disk stars, but should be rare among the main sequence population. We discuss strategies to mitigate the risk of any such false positives.


title: V4142 Sgr: a Double Periodic Variable with an accretor surrounded by the accretion-disk's atmosphere

abstract: Context: A detailed study of the close interacting binary V4142\,Sgr based on photometric and spectroscopic analysis is presented.This system belongs to the enigmatic class of Algol-like variables showing a long photometric cycle of unknown nature. Aims: Performing photometric data-mining and spectroscopic observations covering the orbital cycle, we obtain the orbital parameters and the stellar properties of the binary system, along with the physical properties of the accretion disk located around the hot star. Insights on the evolutive path of the system are obtained. Methods: The light curve was modeled through an inverse modeling method using a theoretical light curve of the binary system, considering the light curve contribution of both stars and the accretion disk of the hot star to obtain the fundamental parameters. To constrain the main stellar parameters the mass ratio was fixed, as well as the donor temperature using the obtained values from our spectroscopic analysis including deblending methods to isolate the spectral lines of the stellar components. The system parameters were compared with a grid of binary star evolutive models in order to get insights on the evolutionary history of the system. Results: The orbital period and the long cycle were re-calculated and found to be of $30.633 \pm 0.002 ~\mathrm{days}$ and $1201 \pm 14 ~\mathrm{days}$. The spectral analysis reveals H$\alpha$ double emission with a persistent $V \leq R$ asymmetry which is considered evidence of a possible wind emergin from the hotspot region...


title: Can Learning Be Explained By Local Optimality In Low-rank Matrix Recovery?

abstract: We explore the local landscape of low-rank matrix recovery, aiming to reconstruct a $d_1\times d_2$ matrix with rank $r$ from $m$ linear measurements, some potentially noisy. When the true rank is unknown, overestimation is common, yielding an over-parameterized model with rank $k\geq r$. Recent findings suggest that first-order methods with the robust $\ell_1$-loss can recover the true low-rank solution even when the rank is overestimated and measurements are noisy, implying that true solutions might emerge as local or global minima. Our paper challenges this notion, demonstrating that, under mild conditions, true solutions manifest as \textit{strict saddle points}. We study two categories of low-rank matrix recovery, matrix completion and matrix sensing, both with the robust $\ell_1$-loss. For matrix sensing, we uncover two critical transitions. With $m$ in the range of $\max\{d_1,d_2\}r\lesssim m\lesssim \max\{d_1,d_2\}k$, none of the true solutions are local or global minima, but some become strict saddle points. As $m$ surpasses $\max\{d_1,d_2\}k$, all true solutions become unequivocal global minima. In matrix completion, even with slight rank overestimation and mild noise, true solutions either emerge as non-critical or strict saddle points.


title: Sedition Hunters: A Quantitative Study of the Crowdsourced Investigation into the 2021 U.S. Capitol Attack

abstract: Social media platforms have enabled extremists to organize violent events, such as the 2021 U.S. Capitol Attack. Simultaneously, these platforms enable professional investigators and amateur sleuths to collaboratively collect and identify imagery of suspects with the goal of holding them accountable for their actions. Through a case study of Sedition Hunters, a Twitter community whose goal is to identify individuals who participated in the 2021 U.S. Capitol Attack, we explore what are the main topics or targets of the community, who participates in the community, and how. Using topic modeling, we find that information sharing is the main focus of the community. We also note an increase in awareness of privacy concerns. Furthermore, using social network analysis, we show how some participants played important roles in the community. Finally, we discuss implications for the content and structure of online crowdsourced investigations.


title: The Neutrino Magnetic Moment Portal and Supernovae: New Constraints and Multimessenger Opportunities

abstract: We scrutinize the hypothesis that gauge singlet fermions -- sterile neutrinos -- interact with Standard Model particles through the transition magnetic moment portal. These interactions lead to the production of sterile neutrinos in supernovae followed by their decay into photons and active neutrinos which can be detected at $\gamma$-ray telescopes and neutrino detectors, respectively. We find that the non-observation of active neutrinos and photons from sterile-neutrino decay associated to SN1987A yields the strongest constraints to date on magnetic-moment-coupled sterile neutrinos if their masses are inside a $0.1-100$ MeV window. Assuming a near-future galactic supernova explosion, we estimate the sensitivity of several present and near-future experiments, including Fermi-LAT, e-ASTROGAM, DUNE, and Hyper-Kamiokande, to magnetic-moment-coupled sterile neutrinos. We also study the diffuse photon and neutrino fluxes produced in the decay of magnetic-moment coupled sterile neutrinos produced in all past supernova explosions and find that the absence of these decay daughters yields the strongest constraints to date for sterile neutrino masses inside a $1-100$ keV window.


title: Dirac bracket and time dependent constraints

abstract: We provide a compact derivation of the Dirac bracket and of the equations of motion for second class constrained systems when the constraints are time dependent. The examples of Parameterized Mechanics and of General Relativity after gauge fixing are given, and the need for the use of time dependent gauge fixing conditions in these examples is illustrated geometrically.


title: Counting rational points on weighted projective spaces over number fields

abstract: Deng (arXiv:math/9812082) gave an asymptotic formula for the number of rational points on a weighted projective space over a number field with respect to a certain height function. We prove a generalization of Deng's result involving a morphism between weighted projective spaces, allowing us to count rational points whose image under this morphism has bounded height. This method provides a more general and simpler proof for a result of the first-named author and Najman on counting elliptic curves with prescribed level structures over number fields. We further include some examples of applications to modular curves.


title: Two-dimensional inversion symmetry as the fundamental symmetry of incompressible quantum Hall fluids

abstract: Two dimensional inversion symmetry ($180^{\circ}$ rotations in the ``Hall plane'' that hosts the incompressible electron fluid that exhibits the quantized Hall effect) is identified as its fundamental unbroken symmetry. A consequence is that the integers $p$ and $q$ which define both the Landau level filling factor $\nu$ = $p/q$ and the elementary fractional charge $\pm e/q$ of topological excitations, cannot have a common divisor greater than 2.


title: Measuring city-scale green infrastructure drawdown dynamics using internet-connected sensors in Detroit

abstract: The impact of green infrastructure (GI) on the urban drainage landscape remains largely unmeasured at high temporal and spatial scales. To that end, a data toolchain is introduced, underpinned by a novel wireless sensor network for continuously measuring real-time water levels in GI. The internet-connected sensors enable the collection of high-resolution data across large regions. A case study in Detroit (MI, US) is presented, where the water levels of 14 GI sites were measured in-situ from June to September 2021. The large dataset is analyzed using an automated storm segmentation methodology, which automatically extracts and analyzes individual storms from measurement time series. Storms are used to parameterize a dynamical system model of GI drawdown dynamics. The model is completely described by the decay constant {\alpha}, which is directly proportional to the drawdown rate. The parameter is analyzed across storms to compare GI dynamics between sites and to determine the major design and physiographic features that drive drawdown dynamics. A correlation analysis using Spearman's rank correlation coefficient reveals that depth to groundwater, imperviousness, longitude, and drainage area to surface area ratio are the most important features explaining GI drawdown dynamics in Detroit. A discussion is provided to contextualize these finding and explore the implications of data-driven strategies for GI design and placement.


title: Differentiable Rendering with Reparameterized Volume Sampling

abstract: In view synthesis, a neural radiance field approximates underlying density and radiance fields based on a sparse set of scene pictures. To generate a pixel of a novel view, it marches a ray through the pixel and computes a weighted sum of radiance emitted from a dense set of ray points. This rendering algorithm is fully differentiable and facilitates gradient-based optimization of the fields. However, in practice, only a tiny opaque portion of the ray contributes most of the radiance to the sum. We propose a simple end-to-end differentiable sampling algorithm based on inverse transform sampling. It generates samples according to the probability distribution induced by the density field and picks non-transparent points on the ray. We utilize the algorithm in two ways. First, we propose a novel rendering approach based on Monte Carlo estimates. This approach allows for evaluating and optimizing a neural radiance field with just a few radiance field calls per ray. Second, we use the sampling algorithm to modify the hierarchical scheme proposed in the original NeRF work. We show that our modification improves reconstruction quality of hierarchical models, at the same time simplifying the training procedure by removing the need for auxiliary proposal network losses.


title: Sidelobe Modeling and Mitigation for a Three Mirror Anastigmat Cosmic Microwave Background Telescope

abstract: Telescopes measuring cosmic microwave background (CMB) polarization on large angular scales require exquisite control of systematic errors to ensure the fidelity of the cosmological results. In particular, far-sidelobe contamination from wide angle scattering is a potentially prominent source of systematic error for large aperture microwave telescopes. Here we describe and demonstrate a ray-tracing-based modeling technique to predict far sidelobes for a Three Mirror Anistigmat (TMA) telescope designed to observe the CMB from the South Pole. Those sidelobes are produced by light scattered in the receiver optics subsequently interacting with the walls of the surrounding telescope enclosure. After comparing simulated sidelobe maps and angular power spectra for different enclosure wall treatments, we propose a highly scattering surface that would provide more than an order of magnitude reduction in the degree-scale far-sidelobe contrast compared to a typical reflective surface. We conclude by discussing the fabrication of a prototype scattering wall panel and presenting measurements of its angular scattering profile.


title: Complexity of Maker-Breaker Games on Edge Sets of Graphs

abstract: We initiate the study of the algorithmic complexity of Maker-Breaker games played on the edge sets of general graphs. We mainly consider the perfect matching game and the $H$-game. Maker wins if she claims the edges of a perfect matching in the first, and a copy of a fixed graph $H$ in the second. We prove that deciding who wins the perfect matching game and the $H$-game is PSPACE-complete, even for the latter in small-diameter graphs if $H$ is a tree. Toward finding the smallest graph $H$ for which the $H$-game is PSPACE-complete, we also prove that such an $H$ of order 51 and size 57 exists. We then give several positive results for the $H$-game. As the $H$-game is already PSPACE-complete when $H$ is a tree, we mainly consider the case where $H$ belongs to a subclass of trees. In particular, we design two linear-time algorithms, both based on structural characterizations, to decide the winners of the $P_4$-game in general graphs and the $K_{1,\ell}$-game in trees. Then, we prove that the $K_{1,\ell}$-game in any graph, and the $H$-game in trees are both FPT parameterized by the length of the game, notably adding to the short list of games with this property, which is of independent interest. Another natural direction to take is to consider the $H$-game when $H$ is a cycle. While we were unable to resolve this case, we prove that the related arboricity-$k$ game is polynomial-time solvable. In particular, when $k=2$, Maker wins this game if she claims the edges of any cycle.


title: Detecting virtual photons in ultrastrongly coupled superconducting quantum circuits

abstract: Light-matter interaction and understanding the fundamental physics behind is essential for emerging quantum technologies. Solid-state devices may explore new regimes where coupling strengths are "ultrastrong", i.e., comparable to the energies of the subsystems. New exotic phenomena occur the common root of many of them being the fact that the entangled vacuum contains virtual photons. They herald the lack of conservation of the number of excitations which is the witness of ultrastrong coupling breaking the U(1) symmetry. Despite more than a decade of research, the detection of ground-state virtual photons still awaits demonstration. In this work, we recognize the "conspiring" set of experimental challenges and show how to overcome them, thus providing a solution to this long-standing problem. We find that combining a superinductor-based unconventional "light fluxonium" qudit and coherent control yields a highly efficient, faithful, and selective conversion of virtual photons into real ones. This enables their detection with resources available to present-day quantum technologies.


title: Electron transport and scattering mechanisms in ferromagnetic monolayer Fe$_3$GeTe$_2$

abstract: We study intrinsic charge-carrier scattering mechanisms and determine their contribution to the transport properties of the two-dimensional ferromagnet Fe$_3$GeTe$_2$. We use state-of-the-art first-principles calculations combined with the model approaches to elucidate the role of the electron-phonon and electron-magnon interactions in the electronic transport. Our findings show that the charge carrier scattering in Fe$_3$GeTe$_2$ is dominated by the electron-phonon interaction, while the role of magnetic excitations is marginal. At the same time, the magnetic ordering is shown to effect essentially on the electron-phonon coupling and its temperature dependence. This leads to a sublinear temperature dependence of the electrical resistivity near the Curie temperature, which is in line with experimental observations. The room temperature resistivity is estimated to be $\sim$35 $\mu \Omega \cdot$cm which may be considered as an intrinsic limit for monolayer Fe$_3$GeTe$_2$.


title: Improved uncertainty quantification for neural networks with Bayesian last layer

abstract: Uncertainty quantification is an important task in machine learning - a task in which standardneural networks (NNs) have traditionally not excelled. This can be a limitation for safety-critical applications, where uncertainty-aware methods like Gaussian processes or Bayesian linear regression are often preferred. Bayesian neural networks are an approach to address this limitation. They assume probability distributions for all parameters and yield distributed predictions. However, training and inference are typically intractable and approximations must be employed. A promising approximation is NNs with Bayesian last layer (BLL). They assume distributed weights only in the linear output layer and yield a normally distributed prediction. To approximate the intractable Bayesian neural network, point estimates of the distributed weights in all but the last layer should be obtained by maximizing the marginal likelihood. This has previously been challenging, as the marginal likelihood is expensive to evaluate in this setting. We present a reformulation of the log-marginal likelihood of a NN with BLL which allows for efficient training using backpropagation. Furthermore, we address the challenge of uncertainty quantification for extrapolation points. We provide a metric to quantify the degree of extrapolation and derive a method to improve the uncertainty quantification for these points. Our methods are derived for the multivariate case and demonstrated in a simulation study. In comparison to Bayesian linear regression with fixed features, and a Bayesian neural network trained with variational inference, our proposed method achieves the highest log-predictive density on test data.


title: Fiber-coupled plug-and-play heralded single photon source based on Ti:LiNbO$_3$ and polymer technology

abstract: A reliable, but cost-effective generation of single-photon states is key for practical quantum communication systems. For real-world deployment, waveguide sources offer optimum compatibility with fiber networks and can be embedded in hybrid integrated modules. Here, we present the first chip-size fully integrated fiber-coupled Heralded Single Photon Source (HSPS) module based on a hybrid integration of a nonlinear lithium niobate waveguide into a polymer board. Photon pairs at 810nm (signal) and 1550nm (idler) are generated via parametric down-conversion pumped at 532nm in the $\mathrm{LiNbO_3}$ waveguide. The pairs are splitted in the polymer board and routed to separate output ports. The module has a size of $(2 \times 1)\, \mathrm{cm^2}$ and is fully fiber-coupled with one pump input fiber and two output fibers. We measure a heralded second-order correlation function of $g_h^{(2)}=0.05$ with a heralding efficiency of $\eta_h=4.5\, \mathrm{\%}$ at low pump powers.


title: HLSDataset: Open-Source Dataset for ML-Assisted FPGA Design using High Level Synthesis

abstract: Machine Learning (ML) has been widely adopted in design exploration using high level synthesis (HLS) to give a better and faster performance, and resource and power estimation at very early stages for FPGA-based design. To perform prediction accurately, high-quality and large-volume datasets are required for training ML models.This paper presents a dataset for ML-assisted FPGA design using HLS, called HLSDataset. The dataset is generated from widely used HLS C benchmarks including Polybench, Machsuite, CHStone and Rossetta. The Verilog samples are generated with a variety of directives including loop unroll, loop pipeline and array partition to make sure optimized and realistic designs are covered. The total number of generated Verilog samples is nearly 9,000 per FPGA type. To demonstrate the effectiveness of our dataset, we undertake case studies to perform power estimation and resource usage estimation with ML models trained with our dataset. All the codes and dataset are public at the github repo.We believe that HLSDataset can save valuable time for researchers by avoiding the tedious process of running tools, scripting and parsing files to generate the dataset, and enable them to spend more time where it counts, that is, in training ML models.


title: Learning to Retrieve Engaging Follow-Up Queries

abstract: Open domain conversational agents can answer a broad range of targeted queries. However, the sequential nature of interaction with these systems makes knowledge exploration a lengthy task which burdens the user with asking a chain of well phrased questions. In this paper, we present a retrieval based system and associated dataset for predicting the next questions that the user might have. Such a system can proactively assist users in knowledge exploration leading to a more engaging dialog. The retrieval system is trained on a dataset which contains ~14K multi-turn information-seeking conversations with a valid follow-up question and a set of invalid candidates. The invalid candidates are generated to simulate various syntactic and semantic confounders such as paraphrases, partial entity match, irrelevant entity, and ASR errors. We use confounder specific techniques to simulate these negative examples on the OR-QuAC dataset and develop a dataset called the Follow-up Query Bank (FQ-Bank). Then, we train ranking models on FQ-Bank and present results comparing supervised and unsupervised approaches. The results suggest that we can retrieve the valid follow-ups by ranking them in higher positions compared to confounders, but further knowledge grounding can improve ranking performance.


title: Anchor-based optimization of energy density functionals

abstract: A new anchor-based optimization method of defining the energy density functionals (EDFs) is proposed. In this approach, the optimization of the parameters of EDF is carried out for the selected set of spherical anchor nuclei the physical observables of which are modified by the correction function which takes into account the global performance of EDF. It is shown that the use of this approach leads to a substantial improvement in global description of binding energies for several classes of covariant EDFs. The computational cost of defining a new functional within this approach is drastically lower as compared with the one for the optimization which includes the global experimental data on spherical, transitional and deformed nuclei into the fitting protocol.


title: MultiRobustBench: Benchmarking Robustness Against Multiple Attacks

abstract: The bulk of existing research in defending against adversarial examples focuses on defending against a single (typically bounded Lp-norm) attack, but for a practical setting, machine learning (ML) models should be robust to a wide variety of attacks. In this paper, we present the first unified framework for considering multiple attacks against ML models. Our framework is able to model different levels of learner's knowledge about the test-time adversary, allowing us to model robustness against unforeseen attacks and robustness against unions of attacks. Using our framework, we present the first leaderboard, MultiRobustBench, for benchmarking multiattack evaluation which captures performance across attack types and attack strengths. We evaluate the performance of 16 defended models for robustness against a set of 9 different attack types, including Lp-based threat models, spatial transformations, and color changes, at 20 different attack strengths (180 attacks total). Additionally, we analyze the state of current defenses against multiple attacks. Our analysis shows that while existing defenses have made progress in terms of average robustness across the set of attacks used, robustness against the worst-case attack is still a big open problem as all existing models perform worse than random guessing.


title: Earthquake magnitude distribution and aftershocks: a statistical geometry explanation

abstract: The emergence of a power-law distribution for the energy released during an earthquake is investigated in several models. Generic features are identified which are based on the self-affine behavior of the stress field prior to an event. This field behaves at large scale as a random trajectory in 1 dimension of space and a random surface in 2 dimensions. Using concepts of statistical mechanics and results on the properties of these random objects, several predictions are obtained and verified, in particular the value of the power-law exponent of the earthquake energy distribution (the Gutenberg-Richter law) as well as a mechanism for the existence of aftershocks after a large earthquake (the Omori law).


title: Robust Majorana bound states in magnetic topological insulator nanoribbons with fragile chiral edge channels

abstract: Magnetic topological insulators in the quantum anomalous Hall regime host ballistic chiral edge channels. When proximitized by an $s$-wave superconductor, these edge states offer the potential for realizing topological superconductivity and Majorana bound states without the detrimental effect of large externally-applied magnetic fields on superconductivity. Realizing well-separated unpaired Majorana bound states requires magnetic topological insulator ribbons with a width of the order of the transverse extent of the edge state, however, which is expected to bring the required ribbon width down to around 100 nm. In this regime, it is known to be extremely difficult to retain the ballistic nature of chiral edge channels and realize a quantized Hall conductance. In this paper, we study the impact of disorder in such magnetic topological insulator nanoribbons and compare the fragility of ballistic chiral edge channels with the stability of Majorana bound states when the ribbon is covered by a superconducting film. We find that the Majorana bound states exhibit greater robustness against disorder than the underlying chiral edge channels.


title: Do Orcas Have Semantic Language? Machine Learning to Predict Orca Behaviors Using Partially Labeled Vocalization Data

abstract: Orcinus orca (killer whales) exhibit complex calls. They last about a second. In a call, an orca typically uses multiple frequencies simultaneously, varies the frequencies, and varies their volumes. Behavior data is hard to obtain because orcas live under water and travel quickly. Sound data is relatively easy to capture. As a science goal, we would like to know whether orca vocalizations constitute a semantic language. We do this by studying whether machine learning can predict behavior from vocalizations. Such prediction would also help scientific research and safety applications because one would like to predict behavior while only having to capture sound. A significant challenge in this process is lack of labeled data. We work with recent recordings of McMurdo Sound orcas [Wellard et al. 2020] where each recording is labeled with the behaviors observed during the recording. This yields a dataset where sound segments - continuous vocalizations that can be thought of as call sequences or more general structures - within the recordings are labeled with superfluous behaviors. Despite that, with a careful combination of recent machine learning techniques, we achieve 96.4% classification accuracy. This suggests that orcas do use a semantic language. It is also promising for research and applications.


title: Discovering Process Models that Support Desired Behavior and Avoid Undesired Behavior

abstract: Process discovery is one of the primary process mining tasks and starting point for process improvements using event data. Existing process discovery techniques aim to find process models that best describe the observed behavior. The focus can be on recall (i.e., replay fitness) or precision. Here, we take a different perspective. We aim to discover a process model that allows for the good behavior observed, and does not allow for the bad behavior. In order to do this, we assume that we have a desirable event log ($L^+$) and an undesirable event log ($L^-$). For example, the desirable event log consists of the cases that were handled within two weeks, and the undesirable event log consists of the cases that took longer. Our discovery approach explores the tradeoff between supporting the cases in the desirable event log and avoiding the cases in the undesirable event log. The proposed framework uses a new inductive mining approach that has been implemented and tested on several real-life event logs. Experimental results show that our approach outperforms other approaches that use only the desirable event log ($L^+$). This supports the intuitive understanding that problematic cases can and should be used to improve processes.


title: Molecular Beam Epitaxy of Twin-Free Bi2Se3 and Sb2Te3 on In2Se3/InP(111)B Virtual Substrates

abstract: Three-dimensional topological insulators (3D-TIs) are a new generation of materials with insulating bulk and exotic metallic surface states that facilitate a wide variety of ground-breaking applications. However, utilization of the surface channels is often hampered by the presence of crystal defects, such as antisites, vacancies and twin domains. For terahertz device applications, twinning is shown to be highly deleterious. Previous attempts to reduce twins using technologically important InP(111) substrates have been promising, but have failed to completely suppress twin domains while preserving high structural quality. Here we report growth of twin-free molecular beam epitaxial Bi2Se3 and Sb2Te3 structures on ultra-thin In2Se3 layers formed by a novel selenium passivation technique during the oxide desorption of smooth, non-vicinal InP(111)B substrates, without the use of an indium source. The formation of un-twinned In2Se3 provides a favorable template to fully suppress twin domains in 3D-TIs, greatly broadening novel device applications in the terahertz regime.