title: The brightest GRB ever detected: GRB 221009A as a highly luminous event at z = 0.151

abstract: Context: The extreme luminosity of gamma-ray bursts (GRBs) makes them powerful beacons for studies of the distant Universe. The most luminous bursts are typically detected at moderate/high redshift, where the volume for seeing such rare events is maximized and the star-formation activity is greater than at z = 0. For distant events, not all observations are feasible, such as at TeV energies. Aims: Here we present a spectroscopic redshift measurement for the exceptional GRB 221009A, the brightest GRB observed to date with emission extending well into the TeV regime. Methods: We used the X-shooter spectrograph at the ESO Very Large Telescope (VLT) to obtain simultaneous optical to near-IR spectroscopy of the burst afterglow 0.5 days after the explosion. Results: The spectra exhibit both absorption and emission lines from material in a host galaxy at z = 0.151. Thus GRB 221009A was a relatively nearby burst with a luminosity distance of 745 Mpc. Its host galaxy properties (star-formation rate and metallicity) are consistent with those of LGRB hosts at low redshift. This redshift measurement yields information on the energy of the burst. The inferred isotropic energy release, $E_{\rm iso} > 5 \times 10^{54}$ erg, lies at the high end of the distribution, making GRB 221009A one of the nearest and also most energetic GRBs observed to date. We estimate that such a combination (nearby as well as intrinsically bright) occurs between once every few decades to once per millennium.


title: Holograms In Our World

abstract: In AdS/CFT, the entanglement wedge EW$(B)$ is the portion of the bulk geometry that can be reconstructed from a boundary region $B$; in other words, EW$(B)$ is the hologram of $B$. We extend this notion to arbitrary spacetimes. Given any gravitating region $a$, we define a max- and a min-entanglement wedge, $e_{\rm max}(a)$ and $e_{\rm min}(a)$, such that $e_{\rm min}(a)\supset e_{\rm max}(a)\supset a$. Unlike their analogues in AdS/CFT, these two spacetime regions can differ already at the classical level, when the generalized entropy is approximated by the area. All information outside $a$ in $e_{\rm max}(a)$ can flow inwards towards $a$, through quantum channels whose capacity is controlled by the areas of intermediate homology surfaces. In contrast, all information outside $e_{\rm min}(a)$ can flow outwards. The generalized entropies of appropriate entanglement wedges obey strong subadditivity, suggesting that they represent the von Neumann entropies of ordinary quantum systems. The entanglement wedges of suitably independent regions satisfy a no-cloning relation. This suggests that it may be possible for an observer in $a$ to summon information from spacelike related points in $e_{\rm max}(a)$, using resources that transcend the semiclassical description of $a$.


title: Quantum control of Rydberg atoms for mesoscopic-scale quantum state and circuit preparation

abstract: Individually trapped Rydberg atoms show significant promise as a platform for scalable quantum simulation and for development of programmable quantum computers. In particular, the Rydberg blockade effect can be used to facilitate both fast qubit-qubit interactions and long coherence times via low-lying electronic states encoding the physical qubits. To bring existing Rydberg-atom-based platforms a step closer to fault-tolerant quantum computation, we demonstrate high-fidelity state and circuit preparation in a system of five atoms. We specifically show that quantum control can be used to reliably generate fully connected cluster states and to simulate the error-correction encoding circuit based on the 'Perfect Quantum Error Correcting Code' by Laflamme et al. [Phys. Rev. Lett. 77, 198 (1996)]. Our results make these ideas and their implementation directly accessible to experiments and demonstrate a promising level of noise tolerance with respect to experimental errors. With this approach, we motivate the application of quantum control in small subsystems in combination with the standard gate-based quantum circuits for direct and high-fidelity implementation of few-qubit modules.


title: Tensor network study of the Shastry-Sutherland model with weak interlayer coupling

abstract: The layered material SrCu$_2$(BO$_3$)$_2$ has long been studied because of its fascinating physics in a magnetic field and under pressure. Many of its properties are remarkably well described by the Shastry-Sutherland model (SSM) - a two-dimensional frustrated spin system. However, the extent of the intermediate plaquette phase discovered in SrCu$_2$(BO$_3$)$_2$ under pressure is significantly smaller than predicted in theory, which is likely due to the weak interlayer coupling that is present in the material but neglected in the model. Using state-of-the-art tensor network methods we study the SSM with a weak interlayer coupling and show that the intermediate plaquette phase is destabilized already at a smaller value around $J''/J\sim0.04-0.05$ than previously predicted from series expansion. Based on our phase diagram we estimate the effective interlayer coupling in SrCu$_2$(BO$_3$)$_2$ to be around $J''/J\sim0.027$ at ambient pressure.


title: Phase transition in Stabilizer Entropy and efficient purity estimation

abstract: Stabilizer Entropy (SE) quantifies the spread of a state in the basis of Pauli operators. It is a computationally tractable measure of non-stabilizerness and thus a useful resource for quantum computation. SE can be moved around a quantum system, effectively purifying a subsystem from its complex features. We show that there is a phase transition in the residual subsystem SE as a function of the density of non-Clifford resources. This phase transition has important operational consequences: it marks the onset of a subsystem purity estimation protocol that requires $poly(n)exp(t)$ many queries to a circuit containing $t$ non-Clifford gates that prepares the state from a stabilizer state. Then, for $t=O(\log_2 n)$, it estimates the purity with polynomial resources and, for highly entangled states, attains an exponential speed-up over the known state-of-the-art algorithms.


title: The Kormendy relation of early-type galaxies as a function of wavelength in Abell S1063, MACS J0416.1-2403, and MACS J1149.5+2223

abstract: The wavelength dependence of the Kormendy relation (KR) is well characterised at low redshift but poorly studied at intermediate redshifts. The KR provides information on the evolution of the population of early-type galaxies (ETGs). Therefore, by studying it, we may shed light on the assembly processes of these objects and their size evolution. As studies at different redshifts are generally conducted in different rest-frame wavebands, it is important to investigate whether the KR is dependent on wavelength. Knowledge of such a dependence is fundamental to correctly interpreting the conclusions we might draw from these studies. We analyse the KRs of the three Hubble Frontier Fields clusters, Abell S1063 (z = 0.348), MACSJ0416.1-2403 (z = 0.396), and MACS J1149.5+2223 (z = 0.542), as a function of wavelength. This is the first time the KR of ETGs has been explored consistently over such a large range of wavelengths at intermediate redshifts. We exploit very deep HST photometry, ranging from the observed B-band to the H-band, and MUSE integral field spectroscopy. We improve the structural parameter estimation we performed in a previous work by means of a newly developed Python package called morphofit. With its use on cluster ETGs, we find that the KR slopes increase smoothly with wavelength from the optical to the near-infrared (NIR) bands in all three clusters, with the intercepts becoming fainter at lower redshifts due to the passive ageing of the ETG stellar populations. The slope trend is consistent with previous findings at lower redshifts. The slope increase with wavelength implies that smaller ETGs are more centrally concentrated than larger ETGs in the NIR with respect to the optical regime. As different bands probe different stellar populations in galaxies, the slope increase also implies that smaller ETGs have stronger internal gradients with respect to larger ETGs.


title: Comment on "Traversable wormhole dynamics on a quantum processor"

abstract: A recent article [Nature 612, 51-55 (2022)] claims to observe traversable wormhole dynamics in an experiment. This claim is based upon performing a teleportation protocol using a Hamiltonian that consists of seven Majorana fermions with five fully-commuting terms. The Hamiltonian is generated via a machine-learning procedure designed to replicate the teleportation behavior of the Sachdev-Ye-Kitaev (SYK) model. The authors claim that the learned Hamiltonian reproduces gravitational dynamics of the SYK model and demonstrates gravitational teleportation through an emergent wormhole. We find: (i) in contrast to these claims, the learned Hamiltonian does not exhibit thermalization; (ii) the teleportation signal only resembles the SYK model for operators that were used in the machine-learning training; (iii) the observed perfect size winding is in fact a generic feature of small-size, fully-commuting models, and does not appear to persist in larger-size fully-commuting models or in non-commuting models at equivalent system sizes

categories: quant-ph cond-mat.quant-gas gr-qc hep-th


title: Can LIGO Detect Non-Annihilating Dark Matter?

abstract: Dark matter from the galactic halo can accumulate in neutron stars and transmute them into sub-2.5 $M_{\odot}$ black holes if the dark matter particles are heavy, stable, and have interactions with nucleons. We show that non-detection of gravitational waves from mergers of such low-mass black holes can constrain the interactions of non-annihilating dark matter particles with nucleons. We find benchmark constraints with LIGO O3 data, viz., $\sigma_{\chi n} \geq {\cal O}(10^{-47})$ cm$^2$ for bosonic DM with $m_\chi\sim$ PeV (or $m_\chi\sim$ GeV, if they can Bose-condense) and $\geq {\cal O}(10^{-46})$ cm$^2$ for fermionic DM with $m_\chi \sim 10^3$ PeV. These bounds depend on the priors on DM parameters and on the currently uncertain binary neutron star merger rate density. However, with increased exposure by the end of this decade, LIGO will probe cross-sections that are many orders of magnitude below the neutrino floor and completely test the dark matter solution to missing pulsars in the Galactic center, demonstrating a windfall science-case for gravitational wave detectors as probes of particle dark matter.


title: Revisiting oldest stars as cosmological probes: new constraints on the Hubble constant

abstract: Despite the tremendous advance of observational cosmology, the value of the Hubble constant ($H_0$) is still controversial (the so called ``Hubble tension'') because of the inconsistency between local/late-time measurements and those derived from the cosmic microwave background. As the age of the Universe is very sensitive to $H_0$, we explored whether the present-day oldest stars could place independent constraints on the Hubble constant. To this purpose, we selected from the literature the oldest objects (globular clusters, stars, white dwarfs, ultra-faint and dwarf spheroidal galaxies) with accurate age estimates. Adopting a conservative prior on their formation redshifts ($11 \leq z_{\rm f} \leq 30$) and assuming $\Omega_{\rm M} = 0.3 \pm 0.02$, we developed a method based on Bayesian statistics to estimate the Hubble constant. We selected the oldest objects ($>13.3$ Gyr) and estimated $H_0$ both for each of them individually and for the average ages of homogeneous subsamples. Statistical and systematic uncertainties were properly taken into account. The constraints based on individual ages indicate that $H_0<70.6$ km/s/Mpc when selecting the most accurate estimates. If the ages are averaged and analyzed independently for each subsample, the most stringent constraints imply $H_0<73.0$ with a probability of 90.3% and errors around 2.5 km/s/Mpc. We also constructed an ``accuracy matrix'' to assess how the constraints on $H_0$ become more stringent with further improvements in the accuracy of stellar ages and $\Omega_{\rm M}$. The results show the high potential of the oldest stars as independent and competitive cosmological probes not only limited to the Hubble constant.


title: Non-perturbative investigation of low eccentricity exterior mean motion resonances

abstract: Mean motion resonances are important in the analysis and understanding of the dynamics of planetary systems. While perturbative approaches have been dominant in many previous studies, recent non-perturbative approaches have revealed novel properties in the low eccentricity regime for interior mean motion resonances of Jupiter in the fundamental model of the circular planar restricted three body model. Here we extend the non-perturbative investigation to exterior mean motion resonances in the low eccentricity regime (up to about 0.1) and for perturber mass in the range 5e-5 to 1e-3 (in units of the central mass). Our results demonstrate that first order exterior resonances have two branches at low eccentricity as well as low-eccentricity bridges connecting neighboring first order resonances. With increasing perturber mass, higher order resonances dissolve into chaos whereas low order resonances persist with larger widths in their radial extent but smaller azimuthal widths. For low order resonances, we also detect secondary resonances arising from small integer commensurabilities between resonant librations and the synodic frequency. These secondary resonances contribute significantly to generating the chaotic sea that typically occurs near mean motion resonances of higher-mass perturbers.


title: Search for scalar induced gravitational waves in the International Pulsar Timing Array Data Release 2 and NANOgrav 12.5 years datasets

abstract: We perform a Bayesian search in the latest Pulsar Timing Array (PTA) datasets for a stochastic gravitational wave (GW) background sourced by curvature perturbations at scales $10^5~\text{Mpc}^{-1}\lesssim k\lesssim 10^8~\text{Mpc}^{-1}$. These re-enter the Hubble horizon at temperatures around and below the QCD crossover phase transition in the early Universe. We include a stochastic background of astrophysical origin in our search and properly account for constraints on the curvature power spectrum from the overproduction of primordial black holes (PBHs). We find that the International PTA Data Release 2 significantly favors the astrophysical model for its reported common-spectrum process, over the curvature-induced background. On the other hand, the two interpretations fit the NANOgrav 12.5 years dataset equally well. We then set new upper limits on the amplitude of the curvature power spectrum at small scales. These are independent from, and competitive with, indirect astrophysical bounds from the abundance of PBH dark matter. Upcoming PTA data releases will provide the strongest probe of the curvature power spectrum around the QCD epoch.


title: Accounting for non-ideal mixing effects in the hydrogen-helium equation of state

abstract: The equation of state for hydrogen and helium is fundamental for studying stars and giant planets. It has been shown that because of interactions at atomic and molecular levels, the behaviour of a mixture of hydrogen and helium cannot be accurately represented by considering these elements separately. This paper aims at providing a simple method to account for interactions between hydrogen and helium in interior and evolution models of giant planets. Using on the one hand ab initio simulations that involve a system of interacting hydrogen and helium particles and pure equations of state for hydrogen and helium on the other, we derived the contributions in density and entropy of the interactions between hydrogen and helium particles. We show that relative variations of up to 15% in density and entropy arise when non-ideal mixing is accounted for. These non-ideal mixing effects must be considered in interior models of giant planets based on accurate gravity field measurements, particularly in the context of variations in the helium-to-hydrogen ratio. They also affect the mass-radius relation of exoplanets. We provide a table that contains the volume and entropy of mixing as a function of pressure and temperature. This table is to be combined with pure hydrogen and pure helium equations of state to obtain an equation of state that self-consistently includes mixing effects for any hydrogen and helium mixing ratio and may be used to model the interior structure and evolution of giant planets to brown dwarfs. Non-linear mixing must be included in accurate calculations of the equations of state of hydrogen and helium. Uncertainties on the equation of state still exist, however. Ab initio calculations of the behaviour of the hydrogen-helium mixture in the megabar regime for various compositions should be performed in order to gain accuracy.


title: Beyond semiclassical time: dynamics in quantum cosmology

abstract: We review two approaches to the definition of the Hilbert space and evolution in mechanical theories with local time-reparametrization invariance, which are often used as toy models of quantum gravity. The first approach is based on the definition of invariant relational observables, whereas the second formalism consists of a perturbative construction of the Hilbert space and a weak-coupling expansion of the Hamiltonian constraint, which is frequently performed as part of the Born-Oppenheimer treatment in quantum cosmology. We discuss in which sense both approaches exhibit an inner product that is gauge-fixed via an operator version of the usual Faddeev-Popov procedure, and, in the second approach, how the unitarity of the effective Schr\"{o}dinger evolution is established perturbatively. We note that a conditional probability interpretation of the physical states is possible, so that both formalisms are examples of quantum mechanics with a relational dynamics.


title: Flaring activity from magnetic reconnection in BL Lacertae

abstract: The evolution of the spectral energy distribution during flares constrains models of particle acceleration in blazar jets. The archetypical blazar BL Lac provided a unique opportunity to study spectral variations during an extended strong flaring episode from 2020-2021. During its brightest $\gamma$-ray state, the observed flux (0.1-300 GeV) reached up to $2.15\,\times\,10^{-5}\,\rm{ph\,cm^{-2}\,s^{-1}}$, with sub-hour scale variability. The synchrotron hump extended into the X-ray regime showing a minute-scale flare with an associated peak shift of inverse-Compton hump in gamma-rays. In shock acceleration models, a high Doppler factor value $>$100 is required to explain the observed rapid variability, change of state, and $\gamma$-ray peak shift. Assuming particle acceleration in mini-jets produced by magnetic reconnection during flares, on the other hand, alleviates the constraint on required bulk Doppler factor. In such jet-in-jet models, observed spectral shift to higher energies (towards TeV regime) and simultaneous rapid variability arises from the accidental alignment of a magnetic plasmoid with the direction of the line of sight. We infer a magnetic field of $\sim0.6\,\rm{G}$ in a reconnection region located at the edge of BLR ($\sim0.02\,\rm{pc}$). The scenario is further supported by log-normal flux distribution arising from merging of plasmoids in reconnection region.


title: Isometric tensor network representations of two-dimensional thermal states

abstract: Tensor networks provide a useful tool to describe low-dimensional complex many-body systems. Finding efficient algorithms to use these methods for finite-temperature simulations in two dimensions is a continuing challenge. Here, we use the class of recently introduced isometric tensor network states, which can also be directly realized with unitary gates on a quantum computer. We utilize a purification ansatz to efficiently represent thermal states of the transverse field Ising model. By performing an imaginary-time evolution starting from infinite temperature, we find that this approach offers a different way with low computational complexity to represent states at finite temperatures.


title: A structured jet explains the extreme GRB 221009A

abstract: Long duration gamma-ray bursts (GRBs) are powerful cosmic explosions, signaling the death of massive stars. Among them, GRB 221009A is by far the brightest burst ever observed. Due to its enormous energy ($E_\textrm{iso}\!\approx$10$^{55}$ erg) and proximity ($z\!\approx$0.15), GRB 221009A is an exceptionally rare event that pushes the limits of our theories. We present multi-wavelength observations covering the first three months of its afterglow evolution. The X-ray brightness decays as a power-law with slope $\approx\!t^{-1.66}$, which is not consistent with standard predictions for jetted emission. We attribute this behavior to a shallow energy profile of the relativistic jet. A similar trend is observed in other energetic GRBs, suggesting that the most extreme explosions may be powered by structured jets launched by a common central engine.


title: Momentum-independent magnetic excitation continuum in the honeycomb iridate H$_3$LiIr$_2$O$_6$

abstract: In the search for realizations of Quantum Spin Liquids (QSL), it is essential to understand the interplay between inherent disorder and the correlated fluctuating spin ground state. H$_3$LiIr$_2$O$_6$ is regarded as a spin liquid proximate to the Kitaev-limit (KQSL) in which H zero-point motion and stacking faults are known to be present. Bond disorder has been invoked to account for the existence of unexpected low-energy spin excitations. Controversy remains about the nature of the underlying correlated state and if any KQSL physics survives. Here, we use resonant X-ray spectroscopies to map the collective excitations in H$_3$LiIr$_2$O$_6$ and characterize its magnetic state. We uncover a broad bandwidth and momentum-independent continuum of magnetic excitations at low temperatures that are distinct from the paramagnetic state. The center energy and high-energy tail of the continuum are consistent with expectations for dominant ferromagnetic Kitaev interactions between dynamically fluctuating spins. The absence of a momentum dependence to these excitations indicates a broken translational invariance. Our data support an interpretation of H$_3$LiIr$_2$O$_6$ as a disordered topological spin liquid in close proximity to bond-disordered versions of the KQSL. Our results shed light on how random disorder affects topological magnetic states and have implications for future experimental and theoretical works toward realizing the Kitaev model in condensed matter systems


title: Linear optical logical Bell state measurements with optimal loss-tolerance threshold

abstract: Quantum threshold theorems impose hard limits on the hardware capabilities to process quantum information. We derive tight and fundamental upper bounds to loss-tolerance thresholds in different linear-optical quantum information processing settings through an adversarial framework, taking into account the intrinsically probabilistic nature of linear optical Bell measurements. For logical Bell state measurements - ubiquitous operations in photonic quantum information - we demonstrate analytically that linear optics can achieve the fundamental loss threshold imposed by the no-cloning theorem even though, following the work of Lee et al., (Phys. Rev. A 100, 052303 (2019)), the constraint was widely assumed to be stricter. We spotlight the assumptions of the latter publication and find their bound holds for a logical Bell measurement built from adaptive physical linear-optical Bell measurements. We also give an explicit even stricter bound for non-adaptive Bell measurements.


title: MAGIC: Manipulating Avatars and Gestures to Improve Remote Collaboration

abstract: Remote collaborative work has become pervasive in many settings, from engineering to medical professions. Users are immersed in virtual environments and communicate through life-sized avatars that enable face-to-face collaboration. Within this context, users often collaboratively view and interact with virtual 3D models, for example, to assist in designing new devices such as customized prosthetics, vehicles, or buildings. However, discussing shared 3D content face-to-face has various challenges, such as ambiguities, occlusions, and different viewpoints that all decrease mutual awareness, leading to decreased task performance and increased errors. To address this challenge, we introduce MAGIC, a novel approach for understanding pointing gestures in a face-to-face shared 3D space, improving mutual understanding and awareness. Our approach distorts the remote user\'s gestures to correctly reflect them in the local user\'s reference space when face-to-face. We introduce a novel metric called pointing agreement to measure what two users perceive in common when using pointing gestures in a shared 3D space. Results from a user study suggest that MAGIC significantly improves pointing agreement in face-to-face collaboration settings, improving co-presence and awareness of interactions performed in the shared space. We believe that MAGIC improves remote collaboration by enabling simpler communication mechanisms and better mutual awareness.


title: Dimensionality dependent electronic structure of the exfoliated van der Waals antiferromagnet NiPS$_3$

abstract: Resonant Inelastic X-ray Scattering (RIXS) was used to measure the local electronic structure in few-layer exfoliated flakes of the van der Waals antiferromagnet NiPS$_3$. The resulting spectra show a systematic softening and broadening of $NiS_6$ multiplet excitations with decreasing layer count from the bulk to three atomic layers (3L). These trends are driven by a decrease in the transition metal-ligand and ligand-ligand hopping integrals, and in the charge-transfer energy: $\Delta$ = 0.60 eV in the bulk and 0.22 eV in 3L NiPS$_3$. Relevant intralayer magnetic exchange integrals computed from the electronic parameters exhibit a systematic decrease in the average interaction strength with thickness and place 2D NiPS$_3$ close to the phase boundary between stripy and spiral antiferromagnetic order, which may explain the apparent vanishing of long-range order in the 2D limit. This study explicitly demonstrates the influence of $inter$layer electronic interactions on $intra$layer ones in insulating magnets. As a consequence, the magnetic Hamiltonian in few-layer insulating magnets can be significantly different from that in the bulk.


title: From Reality Keys to Oraclize. A Deep Dive into the History of Bitcoin Oracles

abstract: Before the advent of alternative blockchains such as Ethereum, the future of decentralization was all in the hands of Bitcoin. Together with Nakamoto itself, early developers were trying to leverage Bitcoin potential to decentralize traditionally centralized applications. However, being Bitcoin a decentralized machine, available non-trustless oracles were considered unsuitable. Therefore, strategies had to be elaborated to solve the so-called oracle problem in the newborn scenario. By interviewing early developers and crawling early forums and repositories, this paper aims to retrace and reconstruct the chain of events and contributions that gave birth to oracles on Bitcoin. The evolution of early trust models and approaches to solving the oracle problem is also outlined. Analyzing technical and social barriers to building oracles on Bitcoin, the transition to Ethereum will also be discussed.

categories: cs.CR cs.CY cs.GT cs.IT econ.GN math.IT q-fin.EC


title: Meeting the Needs of Low-Resource Languages: The Value of Automatic Alignments via Pretrained Models

abstract: Large multilingual models have inspired a new class of word alignment methods, which work well for the model's pretraining languages. However, the languages most in need of automatic alignment are low-resource and, thus, not typically included in the pretraining data. In this work, we ask: How do modern aligners perform on unseen languages, and are they better than traditional methods? We contribute gold-standard alignments for Bribri--Spanish, Guarani--Spanish, Quechua--Spanish, and Shipibo-Konibo--Spanish. With these, we evaluate state-of-the-art aligners with and without model adaptation to the target language. Finally, we also evaluate the resulting alignments extrinsically through two downstream tasks: named entity recognition and part-of-speech tagging. We find that although transformer-based methods generally outperform traditional models, the two classes of approach remain competitive with each other.


title: Heyting frames and Esakia duality

abstract: We introduce the category of Heyting frames and show that it is equivalent to the category of Heyting algebras and dually equivalent to the category of Esakia spaces. This provides a frame-theoretic perspective on Esakia duality for Heyting algebras. We also generalize these results to the setting of Brouwerian algebras and Brouwerian semilattices by introducing the corresponding categories of Brouwerian frames and extending the above equivalences and dual equivalences. This provides a frame-theoretic perspective on generalized Esakia duality for Brouwerian algebras and Brouwerian semilattices.


title: On the Automorphism Group of Token Graphs of Complete Bipartite Graphs

abstract: Let $G$ be a graph of order $n$ and let $k\in \{1,2,\ldots,n-1\}$. The $k$-token graph of $G$ is the graph, whose vertices are all the $k$-subsets of vertices of $G$, where two such $k$-sets are adjacent whenever their symmetric difference is an edge of $G$. In this paper, we determine the automorphism group of the $k$-token graph of the complete bipartite graph $K_{m,n}$.


title: Coupling of ferromagnetic and antiferromagnetic spin dynamics in Mn$_{2}$Au/NiFe thin-film bilayers

abstract: We investigate magnetization dynamics of Mn$_{2}$Au/Py (Ni$_{80}$Fe$_{20}$) thin film bilayers using broadband ferromagnetic resonance (FMR) and Brillouin light scattering spectroscopy. Our bilayers exhibit two resonant modes with zero-field frequencies up to almost 40 GHz, far above the single-layer Py FMR. Our model calculations attribute these modes to the coupling of the Py FMR and the two antiferromagnetic resonance (AFMR) modes of Mn2Au. The coupling-strength is in the order of 1.6 T$\cdot$nm at room temperature for nm-thick Py. Our model reveals the dependence of the hybrid modes on the AFMR frequencies and interfacial coupling as well as the evanescent character of the spin waves that extend across the Mn$_{2}$Au/Py interface


title: The CARMENES search for exoplanets around M dwarfs. Line-by-line sensitivity to activity in M dwarfs

abstract: Radial velocities (RVs) measured from high-resolution stellar spectra are routinely used to detect and characterise orbiting exoplanet companions. The different lines present in stellar spectra are created by several species, which are non-uniformly affected by stellar variability features such as spots or faculae. Stellar variability distorts the shape of the spectral absorption lines from which precise RVs are measured, posing one of the main problems in the study of exoplanets. In this work we aim to study how the spectral lines present in M dwarfs are independently impacted by stellar activity. We used CARMENES optical spectra of six active early- and mid-type M dwarfs to compute line-by-line RVs and study their correlation with several well-studied proxies of stellar activity. We are able to classify spectral lines based on their sensitivity to activity in five M dwarfs displaying high levels of stellar activity. We further used this line classification to compute RVs with activity-sensitive lines and less sensitive lines, enhancing or mitigating stellar activity effects in the RV time series. For specific sets of the least activity-sensitive lines, the RV scatter decreases by ~ 2 to 5 times the initial one, depending on the star. Finally, we compare these lines in the different stars analysed, finding the sensitivity to activity to vary from star to star. Despite the high density of lines and blends present in M dwarf stellar spectra, we find that a line-by-line approach is able to deliver precise RVs. Line-by-line RVs are also sensitive to stellar activity effects, and they allow for an accurate selection of activity-insensitive lines to mitigate activity effects in RV. However, we find stellar activity effects to vary in the same insensitive lines from star to star.


title: Evaluating Trade-offs in Computer Vision Between Attribute Privacy, Fairness and Utility

abstract: This paper investigates to what degree and magnitude tradeoffs exist between utility, fairness and attribute privacy in computer vision. Regarding privacy, we look at this important problem specifically in the context of attribute inference attacks, a less addressed form of privacy. To create a variety of models with different preferences, we use adversarial methods to intervene on attributes relating to fairness and privacy. We see that that certain tradeoffs exist between fairness and utility, privacy and utility, and between privacy and fairness. The results also show that those tradeoffs and interactions are more complex and nonlinear between the three goals than intuition would suggest.


title: A universal sheaf of algebras governing representations of vector fields on quasi-projective varieties

abstract: We construct a quasi-coherent sheaf of associative algebras which controls a category of $AV$-modules over a smooth quasi-projective variety. We establish a local structure theorem, proving that in \'etale charts these associative algebras decompose into a tensor product of the algebra of differential operators and the universal enveloping algebra of the Lie algebra of power series vector fields vanishing at the origin.


title: COVID-VTS: Fact Extraction and Verification on Short Video Platforms

abstract: We introduce a new benchmark, COVID-VTS, for fact-checking multi-modal information involving short-duration videos with COVID19- focused information from both the real world and machine generation. We propose, TwtrDetective, an effective model incorporating cross-media consistency checking to detect token-level malicious tampering in different modalities, and generate explanations. Due to the scarcity of training data, we also develop an efficient and scalable approach to automatically generate misleading video posts by event manipulation or adversarial matching. We investigate several state-of-the-art models and demonstrate the superiority of TwtrDetective.


title: InfoNCE Loss Provably Learns Cluster-Preserving Representations

abstract: The goal of contrasting learning is to learn a representation that preserves underlying clusters by keeping samples with similar content, e.g. the ``dogness'' of a dog, close to each other in the space generated by the representation. A common and successful approach for tackling this unsupervised learning problem is minimizing the InfoNCE loss associated with the training samples, where each sample is associated with their augmentations (positive samples such as rotation, crop) and a batch of negative samples (unrelated samples). To the best of our knowledge, it was unanswered if the representation learned by minimizing the InfoNCE loss preserves the underlying data clusters, as it only promotes learning a representation that is faithful to augmentations, i.e., an image and its augmentations have the same representation. Our main result is to show that the representation learned by InfoNCE with a finite number of negative samples is also consistent with respect to clusters in the data, under the condition that the augmentation sets within clusters may be non-overlapping but are close and intertwined, relative to the complexity of the learning function class.


title: A Deep Learning Technique to Control the Non-linear Dynamics of a Gravitational-wave Interferometer

abstract: In this work we developed a deep learning technique that successfully solves a non-linear dynamic control problem. Instead of directly tackling the control problem, we combined methods in probabilistic neural networks and a Kalman-Filter-inspired model to build a non-linear state estimator for the system. We then used the estimated states to implement a trivial controller for the now fully observable system. We applied this technique to a crucial non-linear control problem that arises in the operation of the LIGO system, an interferometric gravitational-wave observatory. We demonstrated in simulation that our approach can learn from data to estimate the state of the system, allowing a successful control of the interferometer's mirror . We also developed a computationally efficient model that can run in real time at high sampling rate on a single modern CPU core, one of the key requirements for the implementation of our solution in the LIGO digital control system. We believe these techniques could be used to help tackle similar non-linear control problems in other applications.


title: Abundant lattice instability in kagome metal ScV$_6$Sn$_6$

abstract: Kagome materials are emerging platforms for studying charge and spin orders. In this work, we have revealed a rich lattice instability in a $\mathbb{Z}_2$ kagome metal ScV$_6$Sn$_6$ by first-principles calculations. Beyond verifying the $\sqrt3 \times \sqrt3 \times 3$ charge density wave (CDW) order observed by the recent experiment, we further identified three more possible CDW structures, i.e., $\sqrt3 \times \sqrt3 \times 2$ CDW with $P6/mmm$ symmetry, $2 \times 2 \times 2$ CDW with $Immm$ symmetry, and $2 \times 2 \times 2$ CDW with $P6/mmm$ symmetry. The former two are more energetically favored than the $\sqrt3 \times \sqrt3 \times 3$ phase, while the third one is comparable in energy. These CDW distortions involve mainly out-of-plane motions of Sc and Sn atoms, while V atoms constituting the kagome net are almost unchanged. We attribute the lattice instability to the smallness of Sc atomic radius. In contrast, such an instability disappears in its sister compounds $R$V$_6$Sn$_6$ ($R$ is Y, or a rare-earth element) because $R$ has a larger radius. Our work indicates that ScV$_6$Sn$_6$ might exhibit varied CDW phases in different experimental conditions and provides insights to explore rich charge orders in kagome materials.


title: Possible Existence of Extremely Neutron-Rich Superheavy Nuclei in Neutron Star Crusts Under a Superstrong Magnetic Field

abstract: We investigate outer crust compositions for a wide range of magnetic field strengths, up to $B\simeq 4\times10^{18}$ G, employing the latest experimental nuclear masses supplemented with various mass models. The essential effects of the magnetic field are twoholds: 1) Enhancement of electron fraction, which is connected to that of protons via the charge neutrality condition, due to the Landau-Rabi quantization of electron motion perpendicular to the field, namely, neutron-richness is suppressed for a given pressure. As a result, 2) nuclei can exist at higher pressures without dripping out neutrons. By exploring optimal outer-crust compositions from all possible nuclei predicted by theoretical models, we find that neutron-rich heavy nuclei with neutron magic numbers 50, 82, 126, as well as 184, with various proton numbers emerge for $B\gtrsim 10^{17}$ G. Moreover, we show that superheavy nuclei with proton numbers $Z\ge104$, including unknown elements such as $Z=119, 120, 122$, and/or $124$, depending on mass models, may emerge as an equilibrium composition at bottom layers of the outer crust for $B\gtrsim 10^{18}$ G, which are extremely neutron-rich ($N\approx 260$-$287$, i.e., $N/Z\approx 2.2$-$2.4$). We point out that those extremely neutron-rich superheavy nuclei locate around the next neutron magic number $N=258$ after $N=184$, underlining importance of nuclear structure calculations under such really exotic, extreme conditions. We demonstrate how the superstrong magnetic field substantially alters crustal properties of neutron stars, which may have detectable consequences.


title: A Spherical Shells Model of Atmospheric Absorption for Instrument Calibration

abstract: We present a model for atmospheric absorption of solar ultraviolet (UV) radiation. The initial motivation for this work is to predict this effect and correct it in Sounding Rocket (SR) experiments. In particular, the Full-sun Ultraviolet Rocket Spectrograph (FURST) is anticipated to launch in mid-2023. FURST has the potential to observe UV absorption while imaging solar spectra between 120-181 nm, at a resolution of R > 2x10$^4$ ($\Delta$ V < $\pm$ 15 km/s), and at altitudes of between 110-255 km. This model uses estimates for density and temperature, as well as laboratory measurements of the absorption cross-section, to predict the UV absorption of solar radiation at high altitudes. Refraction correction is discussed and partially implemented but is negligible for the results presented. Absorption by molecular Oxygen is the primary driver within the UV spectral range of our interest. The model is built with a wide range of applications in mind. The primary result is a method for inversion of the absorption cross-section from images obtained during an instrument flight, even if atmospheric observations were not initially intended. The potential to obtain measurements of atmospheric properties is an exciting prospect, especially since sounding rockets are the only method currently available for probing this altitude in situ. Simulation of noisy spectral images along the FURST flight profile is performed using data from the High-Resolution Telescope and Spectrograph (HRTS) SR and the FISM2 model for comparison. Analysis of these simulated signals allows us to capture the Signal-to-Noise Ratio (SNR) of FURST and the capability to measure atmospheric absorption properties as a function of altitude. Based on the prevalence of distinct spectral features, our calculations demonstrate that atmospheric absorption may be used to perform wavelength calibration from in-flight data.


title: AGN Accretion Discs

abstract: I present the elements of accretion-disc physics applied to active galactic nuclei. The accretion driving mechanisms are discussed, then models of geometrically-thin discs, both stationary and time-dependent, are addressed. Disc's self-gravitation in the AGN context is presented. The shapes of spectral line from accretion discs are explained in both newtonian and relativistic cases. The physics of disc's thermal and viscous instabilities is dicussed in detail and models are applied to AGN discs. Finally, the problems of thick accretion flows (ADAFs and slim discs), disc coronae, winds and jets are discussed in some detail. This is a graduate-student level lecture, not a review article


title: Commonsense Reasoning for Conversational AI: A Survey of the State of the Art

abstract: Large, transformer-based pretrained language models like BERT, GPT, and T5 have demonstrated a deep understanding of contextual semantics and language syntax. Their success has enabled significant advances in conversational AI, including the development of open-dialogue systems capable of coherent, salient conversations which can answer questions, chat casually, and complete tasks. However, state-of-the-art models still struggle with tasks that involve higher levels of reasoning - including commonsense reasoning that humans find trivial. This paper presents a survey of recent conversational AI research focused on commonsense reasoning. The paper lists relevant training datasets and describes the primary approaches to include commonsense in conversational AI. The paper also discusses benchmarks used for evaluating commonsense in conversational AI problems. Finally, the paper presents preliminary observations of the limited commonsense capabilities of two state-of-the-art open dialogue models, BlenderBot3 and LaMDA, and its negative effect on natural interactions. These observations further motivate research on commonsense reasoning in conversational AI.


title: Less is different: why sparse networks with inhibition differ from complete graphs

abstract: In neuronal systems, inhibition contributes to stabilizing dynamics and regulating pattern formation. Through developing mean field theories of neuronal models, using complete graph networks, inhibition is commonly viewed as one ``control parameter'' of the system, promoting an absorbing phase transition. Here, we show that for low connectivity sparse networks, inhibition weight is not a control parameter of the transition. We present analytical and simulation results using generic stochastic integrate-and-fire neurons that, under specific restrictions, become other simpler stochastic neuron models common in literature, which allow us to show that our results are valid for those models as well. We also give a simple explanation about why the inhibition role depends on topology, even when the topology has a dimensionality greater than the critical one. The absorbing transition independence of the inhibitory weight may be an important feature of a sparse network, as it will allow the network to maintain a near-critical regime, self-tuning average excitation, but at the same time, have the freedom to adjust inhibitory weights for computation, learning, and memory, exploiting the benefits of criticality.


title: Multi-Channel Target Speaker Extraction with Refinement: The WavLab Submission to the Second Clarity Enhancement Challenge

abstract: This paper describes our submission to the Second Clarity Enhancement Challenge (CEC2), which consists of target speech enhancement for hearing-aid (HA) devices in noisy-reverberant environments with multiple interferers such as music and competing speakers. Our approach builds upon the powerful iterative neural/beamforming enhancement (iNeuBe) framework introduced in our recent work, and this paper extends it for target speaker extraction. We therefore name the proposed approach as iNeuBe-X, where the X stands for extraction. To address the challenges encountered in the CEC2 setting, we introduce four major novelties: (1) we extend the state-of-the-art TF-GridNet model, originally designed for monaural speaker separation, for multi-channel, causal speech enhancement, and large improvements are observed by replacing the TCNDenseNet used in iNeuBe with this new architecture; (2) we leverage a recent dual window size approach with future-frame prediction to ensure that iNueBe-X satisfies the 5 ms constraint on algorithmic latency required by CEC2; (3) we introduce a novel speaker-conditioning branch for TF-GridNet to achieve target speaker extraction; (4) we propose a fine-tuning step, where we compute an additional loss with respect to the target speaker signal compensated with the listener audiogram. Without using external data, on the official development set our best model reaches a hearing-aid speech perception index (HASPI) score of 0.942 and a scale-invariant signal-to-distortion ratio improvement (SI-SDRi) of 18.8 dB. These results are promising given the fact that the CEC2 data is extremely challenging (e.g., on the development set the mixture SI-SDR is -12.3 dB). A demo of our submitted system is available at WAVLab CEC2 demo.


title: Searching for TeV Dark Matter in Irregular dwarf galaxies with HAWC Observatory

abstract: We present the results of dark matter (DM) searches in a sample of 31 dwarf irregular (dIrr) galaxies within the field of view of the HAWC Observatory. dIrr galaxies are DM dominated objects, which astrophysical gamma-ray emission is estimated to be negligible with respect to the secondary gamma-ray flux expected by annihilation or decay of Weakly Interacting Massive Particles (WIMPs). While we do not see any statistically significant DM signal in dIrr galaxies, we present the exclusion limits ($95\%~\text{C.L.}$) for annihilation cross-section and decay lifetime for WIMP candidates with masses between $1$ and $100~\text{TeV}$. Exclusion limits from dIrr galaxies are relevant and complementary to benchmark dwarf Spheroidal (dSph) galaxies. In fact, dIrr galaxies are targets kinematically different from benchmark dSph, preserving the footprints of different evolution histories. We compare the limits from dIrr galaxies to those from ultrafaint and classical dSph galaxies previously observed with HAWC. We find that the contraints are comparable to the limits from classical dSph galaxies and $\thicksim2$ orders of magnitude weaker than the ultrafaint dSph limits.


title: Interpretable Deep Learning Methods for Multiview Learning

abstract: Technological advances have enabled the generation of unique and complementary types of data or views (e.g. genomics, proteomics, metabolomics) and opened up a new era in multiview learning research with the potential to lead to new biomedical discoveries. We propose iDeepViewLearn (Interpretable Deep Learning Method for Multiview Learning) for learning nonlinear relationships in data from multiple views while achieving feature selection. iDeepViewLearn combines deep learning flexibility with the statistical benefits of data and knowledge-driven feature selection, giving interpretable results. Deep neural networks are used to learn view-independent low-dimensional embedding through an optimization problem that minimizes the difference between observed and reconstructed data, while imposing a regularization penalty on the reconstructed data. The normalized Laplacian of a graph is used to model bilateral relationships between variables in each view, therefore, encouraging selection of related variables. iDeepViewLearn is tested on simulated and two real-world data, including breast cancer-related gene expression and methylation data. iDeepViewLearn had competitive classification results and identified genes and CpG sites that differentiated between individuals who died from breast cancer and those who did not. The results of our real data application and simulations with small to moderate sample sizes suggest that iDeepViewLearn may be a useful method for small-sample-size problems compared to other deep learning methods for multiview learning.


title: ANSEL Photobot: A Robot Event Photographer with Semantic Intelligence

abstract: Our work examines the way in which large language models can be used for robotic planning and sampling, specifically the context of automated photographic documentation. Specifically, we illustrate how to produce a photo-taking robot with an exceptional level of semantic awareness by leveraging recent advances in general purpose language (LM) and vision-language (VLM) models. Given a high-level description of an event we use an LM to generate a natural-language list of photo descriptions that one would expect a photographer to capture at the event. We then use a VLM to identify the best matches to these descriptions in the robot's video stream. The photo portfolios generated by our method are consistently rated as more appropriate to the event by human evaluators than those generated by existing methods.


title: The muon deficit problem: a new method to calculate the muon rescaling factors and the Heitler-Matthews beta exponent

abstract: Simulations of extensive air showers using current hadronic interaction models predict too small numbers of muons compared to events observed in the air-shower experiments, which is known as the muon-deficit problem. In this work, we present a new method to calculate the factor by which the muon signal obtained via Monte-Carlo simulations must be rescaled to match the data, as well as the beta exponent from the Heitler-Matthews model which governs the number of muons found in an extensive air shower as a function of the mass and the energy of the primary cosmic ray. This method uses the so-called z variable (difference between the total reconstructed and the simulated signals), which is connected to the muon signal and is roughly independent of the zenith angle, but depends on the mass of the primary cosmic ray. Using a mock dataset built from QGSJetII-04, we show that such a method allows us to reproduce the average muon signal from this dataset using Monte-Carlo events generated with the EPOS-LHC hadronic model, with accuracy better than 6%. As a consequence of the good recovery of the muon signal for each primary included in the analysis, also the beta exponent can be obtained with accuracy of less than 1% for the studied system. Detailed simulations show a dependence of the beta exponent on hadronic interaction properties, thus the determination of this parameter is important for understanding the muon deficit problem.


title: RADIUS: Risk-Aware, Real-Time, Reachability-Based Motion Planning

abstract: Deterministic methods for motion planning guarantee safety amidst uncertainty in obstacle locations by trying to restrict the robot from operating in any possible location that an obstacle could be in. Unfortunately, this can result in overly conservative behavior. Chance-constrained optimization can be applied to improve the performance of motion planning algorithms by allowing for a user-specified amount of bounded constraint violation. However, state-of-the-art methods rely either on moment-based inequalities, which can be overly conservative, or make it difficult to satisfy assumptions about the class of probability distributions used to model uncertainty. To address these challenges, this work proposes a real-time, risk-aware reachability-based motion planning framework called RADIUS. The method first generates a reachable set of parameterized trajectories for the robot offline. At run time, RADIUS computes a closed-form over-approximation of the risk of a collision with an obstacle. This is done without restricting the probability distribution used to model uncertainty to a simple class (e.g., Gaussian). Then, RADIUS performs real-time optimization to construct a trajectory that can be followed by the robot in a manner that is certified to have a risk of collision that is less than or equal to a user-specified threshold. The proposed algorithm is compared to several state-of-the-art chance-constrained and deterministic methods in simulation, and is shown to consistently outperform them in a variety of driving scenarios. A demonstration of the proposed framework on hardware is also provided.


title: Profiling Cold New Early Dark Energy

abstract: Recent interest in New Early Dark Energy (NEDE), a cosmological model with a vacuum energy component decaying in a triggered phase transition around recombination, has been sparked by its impact on the Hubble tension. Previous constraints on the model parameters were derived in a Bayesian framework with Markov-chain Monte Carlo (MCMC) methods. In this work, we instead perform a frequentist analysis using the profile likelihood in order to assess the impact of prior volume effects on the constraints. We constrain the maximal fraction of NEDE $f_\mathrm{NEDE}$, finding $f_\mathrm{NEDE}=0.076^{+0.040}_{-0.035}$ at $68 \%$ CL with our baseline dataset and similar constraints using either data from SPT-3G, ACT or full-shape large-scale structure, showing a preference over $\Lambda$CDM even in the absence of a SH0ES prior on $H_0$. While this is stronger evidence for NEDE than obtained with the corresponding Bayesian analysis, our constraints broadly match those obtained by fixing the NEDE trigger mass. Including the SH0ES prior on $H_0$, we obtain $f_\mathrm{NEDE}=0.136^{+0.024}_{-0.026}$ at $68 \%$ CL. Furthermore, we compare NEDE with the Early Dark Energy (EDE) model, finding similar constraints on the maximal energy density fractions and $H_0$ in the two models. At $68 \%$ CL in the NEDE model, we find $H_0 = 69.56^{+1.16}_{-1.29} \text{ km s}^{-1}\text{ Mpc}^{-1}$ with our baseline and $H_0 = 71.62^{+0.78}_{-0.76} \text{ km s}^{-1}\text{ Mpc}^{-1}$ when including the SH0ES measurement of $H_0$, thus corroborating previous conclusions that the NEDE model provides a considerable alleviation of the $H_0$ tension.


title: Market-Based Probability of Stock Returns

abstract: Markets possess all available information on stock returns. The randomness of market trade determines the statistics of stock returns. This paper describes the dependence of the first four market-based statistical moments of stock returns on statistical moments and correlations of current and past trade values. The mean return of trades during the averaging period coincides with Markowitz's definition of portfolio value weighted return. We derive the market-based volatility of return and return-value correlations. We present approximations of the characteristic functions and probability measures of stock return by a finite number of market-based statistical moments. To forecast market-based average return or volatility of return, one should predict the statistical moments and correlations of current and past market trade values at the same time horizon.


title: Galaxy cluster rotation revealed in the MACSIS simulations with the kinetic Sunyaev-Zeldovich effect

abstract: The kinetic Sunyaev-Zeldovich (kSZ) effect has now become a clear target for ongoing and future studies of the cosmic microwave background (CMB) and cosmology. Aside from the bulk cluster motion, internal motions also lead to a kSZ signal. In this work, we study the rotational kSZ effect caused by coherent large-scale motions of the cluster medium using cluster hydrodynamic cosmological simulations. To utilise the rotational kSZ as a cosmological probe, simulations offer some of the most comprehensive data sets that can inform the modelling of this signal. In this work, we use the MACSIS data set to investigate the rotational kSZ effect in massive clusters specifically. Based on these models, we test stacking approaches and estimate the amplitude of the combined signal with varying mass, dynamical state, redshift and map-alignment geometry. We find that the dark matter, galaxy and gas spins are generally misaligned, an effect that can cause a sub-optimal estimation of the rotational kSZ effect when based on galaxy motions. Furthermore, we provide halo-spin-mass scaling relations that can be used to build a statistical model of the rotational kSZ. The rotational kSZ contribution, which is largest in massive unrelaxed clusters ($\gtrsim$100 $\mu$K), could be relevant to studies of higher-order CMB temperature signals, such as the moving lens effect. The limited mass range of the MACSIS sample strongly motivates an extended investigation of the rotational kSZ effect in large-volume simulations to refine the modelling, particularly towards lower mass and higher redshift, and provide forecasts for upcoming cosmological CMB experiments (e.g. Simons Observatory, SKA-2) and X-ray observations (e.g. Athena/X-IFU).


title: The Expressive Power of Tuning Only the Normalization Layers

abstract: Feature normalization transforms such as Batch and Layer-Normalization have become indispensable ingredients of state-of-the-art deep neural networks. Recent studies on fine-tuning large pretrained models indicate that just tuning the parameters of these affine transforms can achieve high accuracy for downstream tasks. These findings open the questions about the expressive power of tuning the normalization layers of frozen networks. In this work, we take the first step towards this question and show that for random ReLU networks, fine-tuning only its normalization layers can reconstruct any target network that is $O(\sqrt{\text{width}})$ times smaller. We show that this holds even for randomly sparsified networks, under sufficient overparameterization, in agreement with prior empirical work.


title: Scalable Multi-Agent Reinforcement Learning with General Utilities

abstract: We study the scalable multi-agent reinforcement learning (MARL) with general utilities, defined as nonlinear functions of the team's long-term state-action occupancy measure. The objective is to find a localized policy that maximizes the average of the team's local utility functions without the full observability of each agent in the team. By exploiting the spatial correlation decay property of the network structure, we propose a scalable distributed policy gradient algorithm with shadow reward and localized policy that consists of three steps: (1) shadow reward estimation, (2) truncated shadow Q-function estimation, and (3) truncated policy gradient estimation and policy update. Our algorithm converges, with high probability, to $\epsilon$-stationarity with $\widetilde{\mathcal{O}}(\epsilon^{-2})$ samples up to some approximation error that decreases exponentially in the communication radius. This is the first result in the literature on multi-agent RL with general utilities that does not require the full observability.


title: Retrieval study of cool, directly imaged exoplanet 51 Eri b

abstract: Retrieval methods are a powerful analysis technique for modelling exoplanetary atmospheres by estimating the bulk physical and chemical properties that combine in a forward model to best-fit an observed spectrum, and they are increasingly being applied to observations of directly-imaged exoplanets. We have adapted TauREx3, the Bayesian retrieval suite, for the analysis of near-infrared spectrophotometry from directly-imaged gas giant exoplanets and brown dwarfs. We demonstrate TauREx3's applicability to sub-stellar atmospheres by presenting results for brown dwarf benchmark GJ 570D which are consistent with previous retrieval studies, whilst also exhibiting systematic biases associated with the presence of alkali lines. We also present results for the cool exoplanet 51 Eri b, the first application of a free chemistry retrieval analysis to this object, using spectroscopic observations from GPI and SPHERE. While our retrieval analysis is able to explain spectroscopic and photometric observations without employing cloud extinction, we conclude this may be a result of employing a flexible temperature-pressure profile which is able to mimic the presence of clouds. We present Bayesian evidence for an ammonia detection with a 2.7$\sigma$ confidence, the first indication of ammonia in an exoplanetary atmosphere. This is consistent with this molecule being present in brown dwarfs of a similar spectral type. We demonstrate the chemical similarities between 51 Eri b and GJ 570D in relation to their retrieved molecular abundances. Finally, we show that overall retrieval conclusions for 51 Eri b can vary when employing different spectral data and modelling components, such as temperature-pressure and cloud structures.


title: Online Tool Selection with Learned Grasp Prediction Models

abstract: Deep learning-based grasp prediction models have become an industry standard for robotic bin-picking systems. To maximize pick success, production environments are often equipped with several end-effector tools that can be swapped on-the-fly, based on the target object. Tool-change, however, takes time. Choosing the order of grasps to perform, and corresponding tool-change actions, can improve system throughput; this is the topic of our work. The main challenge in planning tool change is uncertainty - we typically cannot see objects in the bin that are currently occluded. Inspired by queuing and admission control problems, we model the problem as a Markov Decision Process (MDP), where the goal is to maximize expected throughput, and we pursue an approximate solution based on model predictive control, where at each time step we plan based only on the currently visible objects. Special to our method is the idea of void zones, which are geometrical boundaries in which an unknown object will be present, and therefore cannot be accounted for during planning. Our planning problem can be solved using integer linear programming (ILP). However, we find that an approximate solution based on sparse tree search yields near optimal performance at a fraction of the time. Another question that we explore is how to measure the performance of tool-change planning: we find that throughput alone can fail to capture delicate and smooth behavior, and propose a principled alternative. Finally, we demonstrate our algorithms on both synthetic and real world bin picking tasks.


title: An Experimentation Infrastructure for Quantitative Measurements of Cyber Resilience

abstract: The vulnerability of cyber-physical systems to cyber attack is well known, and the requirement to build cyber resilience into these systems has been firmly established. The key challenge this paper addresses is that maturing this discipline requires the development of techniques, tools, and processes for objectively, rigorously, and quantitatively measuring the attributes of cyber resilience. Researchers and program managers need to be able to determine if the implementation of a resilience solution actually increases the resilience of the system. In previous work, a table top exercise was conducted using a notional heavy vehicle on a fictitious military mission while under a cyber attack. While this exercise provided some useful data, more and higher fidelity data is required to refine the measurement methodology. This paper details the efforts made to construct a cost-effective experimentation infrastructure to provide such data. It also presents a case study using some of the data generated by the infrastructure.


title: Enhancing Deep Knowledge Tracing with Auxiliary Tasks

abstract: Knowledge tracing (KT) is the problem of predicting students' future performance based on their historical interactions with intelligent tutoring systems. Recent studies have applied multiple types of deep neural networks to solve the KT problem. However, there are two important factors in real-world educational data that are not well represented. First, most existing works augment input representations with the co-occurrence matrix of questions and knowledge components\footnote{\label{ft:kc}A KC is a generalization of everyday terms like concept, principle, fact, or skill.} (KCs) but fail to explicitly integrate such intrinsic relations into the final response prediction task. Second, the individualized historical performance of students has not been well captured. In this paper, we proposed \emph{AT-DKT} to improve the prediction performance of the original deep knowledge tracing model with two auxiliary learning tasks, i.e., \emph{question tagging (QT) prediction task} and \emph{individualized prior knowledge (IK) prediction task}. Specifically, the QT task helps learn better question representations by predicting whether questions contain specific KCs. The IK task captures students' global historical performance by progressively predicting student-level prior knowledge that is hidden in students' historical learning interactions. We conduct comprehensive experiments on three real-world educational datasets and compare the proposed approach to both deep sequential KT models and non-sequential models. Experimental results show that \emph{AT-DKT} outperforms all sequential models with more than 0.9\% improvements of AUC for all datasets, and is almost the second best compared to non-sequential models. Furthermore, we conduct both ablation studies and quantitative analysis to show the effectiveness of auxiliary tasks and the superior prediction outcomes of \emph{AT-DKT}.


title: A unified gas-kinetic particle method for frequency-dependent radiative transfer equations with isotropic scattering process on unstructured mesh

abstract: In this paper, we extend the unified kinetic particle (UGKP) method to the frequency-dependent radiative transfer equation with both absorption-emission and scattering processes. The extended UGKP method could not only capture the diffusion and free transport limit, but also provide a smooth transition in the physical and frequency space in the regime between the above two limits. The proposed scheme has the properties of asymptotic-preserving, regime-adaptive, and entropy-preserving, which make it an accurate and efficient scheme in the simulation of multiscale photon transport problems. The methodology of scheme construction is a coupled evolution of macroscopic energy equation and the microscopic radiant intensity equation, where the numerical flux in macroscopic energy equation and the closure in microscopic radiant intensity equation are constructed based on the integral solution. Both numerical dissipation and computational complexity are well controlled especially in the optical thick regime. A 2D multi-thread code on a general unstructured mesh has been developed. Several numerical tests have been simulated to verify the numerical scheme and code, covering a wide range of flow regimes. The numerical scheme and code that we developed are highly demanded and widely applicable in the high energy density engineering applications.


title: Effective Data Augmentation With Diffusion Models

abstract: Data augmentation is one of the most prevalent tools in deep learning, underpinning many recent advances, including those from classification, generative models, and representation learning. The standard approach to data augmentation combines simple transformations like rotations and flips to generate new images from existing ones. However, these new images lack diversity along key semantic axes present in the data. Current augmentations cannot alter the high-level semantic attributes, such as animal species present in a scene, to enhance the diversity of data. We address the lack of diversity in data augmentation with image-to-image transformations parameterized by pre-trained text-to-image diffusion models. Our method edits images to change their semantics using an off-the-shelf diffusion model, and generalizes to novel visual concepts from a few labelled examples. We evaluate our approach on few-shot image classification tasks, and on a real-world weed recognition task, and observe an improvement in accuracy in tested domains.


title: An implicit unified gas-kinetic wave-particle method for radiative transport process

abstract: The unified gas-kinetic wave-particle method (UGKWP) has been developed for the multiscale gas, plasma, and multiphase flow transport processes for the past years. In this work, we propose an implicit unified gas-kinetic wave-particle (IUGKWP) method to remove the CFL time step constraint. Based on the local integral solution of the radiative transfer equation (RTE), the particle transport processes are categorized into the long-$\lambda$ streaming process and the short-$\lambda$ streaming process comparing to a local physical characteristic time $t_p$. In the construction of the IUGKWP method, the long-$\lambda$ streaming process is tracked by the implicit Monte Carlo (IMC) method; the short-$\lambda$ streaming process is evolved by solving the implicit moments equations; and the photon distribution is closed by a local integral solution of RTE. In the IUGKWP method, the multiscale flux of radiation energy and the multiscale closure of photon distribution are constructed based on the local integral solution. The IUGKWP method preserves the second-order asymptotic expansion of RTE in the optically thick regime and adapts its computational complexity to the flow regime. The numerical dissipation is well controlled, and the teleportation error is significantly reduced in the optically thick regime. The computational complexity of the IUGKWP method decreases exponentially as the Knudsen number approaches zero, and the computational efficiency is remarkably improved in the optically thick regime. The IUGKWP is formulated on a generalized unstructured mesh, and multidimensional 2D and 3D algorithms are developed. Numerical tests are presented to validate the capability of IUGKWP in capturing the multiscale photon transport process. The algorithm and code will apply in the engineering applications of inertial confinement fusion (ICF).


title: Experimenting with Emerging RISC-V Systems for Decentralised Machine Learning

abstract: Decentralised Machine Learning (DML) enables collaborative machine learning without centralised input data. Federated Learning (FL) and Edge Inference are examples of DML. While tools for DML (especially FL) are starting to flourish, many are not flexible and portable enough to experiment with novel processors (e.g., RISC-V), non-fully connected network topologies, and asynchronous collaboration schemes. We overcome these limitations via a domain-specific language allowing us to map DML schemes to an underlying middleware, i.e. the FastFlow parallel programming library. We experiment with it by generating different working DML schemes on x86-64 and ARM platforms and an emerging RISC-V one. We characterise the performance and energy efficiency of the presented schemes and systems. As a byproduct, we introduce a RISC-V porting of the PyTorch framework, the first publicly available to our knowledge.


title: A Sequential Quadratic Programming Method for Optimization with Stochastic Objective Functions, Deterministic Inequality Constraints and Robust Subproblems

abstract: In this paper, a robust sequential quadratic programming method of [1] for constrained optimization is generalized to problem with stochastic objective function, deterministic equality and inequality constraints. A stochastic line search scheme in [2] is employed to globalize the steps. We show that in the case where the algorithm fails to terminate in finite number of iterations, the sequence of iterates will converge almost surely to a Karush-Kuhn-Tucker point under the assumption of extended Mangasarian-Fromowitz constraint qualification. We also show that, with a specific sampling method, the probability of the penalty parameter approaching infinity is 0. Encouraging numerical results are reported.


title: Effect of dynamical screening in the Bethe-Salpeter framework: Excitons in crystalline naphthalene

abstract: Solving the Bethe-Salpeter equation (BSE) for the optical polarization functions is a first principles means to model optical properties of materials including excitonic effects. One almost ubiquitously used approximation neglects the frequency dependence of the screened electron-hole interaction. This is commonly justified by the large difference in magnitude of electronic plasma frequency and exciton binding energy. We incorporated dynamical effects into the screening of the electron-hole interaction in the BSE using two different approximations as well as exact diagonalization of the exciton Hamiltonian. We compare these approaches for a naphthalene organic crystal, for which the difference between exciton binding energy and plasma frequency is only about a factor of ten. Our results show that in this case, corrections due to dynamical screening are about 15\,\% of the exciton binding energy. We analyze the effect of screening dynamics on optical absorption across the visible spectral range and use our data to establish an \emph{effective} screening model as a computationally efficient approach to approximate dynamical effects in complex materials in the future.


title: Second Quantization: Gating a Quantum Dot Through the Sequential Removal of Single Electrons from a Nanoscale Floating Gate

abstract: We use the tip of an atomic force microscope (AFM) to charge floating metallic gates defined on the surface of a Si/SiGe heterostructure. The AFM tip serves as an ideal and movable cryogenic switch, allowing us to bias a floating gate to a specific voltage and then lock the charge on the gate by withdrawing the tip. Biasing with an AFM tip allows us to reduce the size of a quantum dot floating gate electrode down to $\sim100~\mathrm{nm}$. Measurements of the conductance through a quantum dot formed beneath the floating gate indicate that its charge changes in discrete steps. From the statistics of the single-electron leakage events, we determine the floating gate leakage resistance $R \sim 10^{19}~ \mathrm{Ohm}$ - a value immeasurable by conventional means.


title: Topological Neural Discrete Representation Learning \`a la Kohonen

abstract: Unsupervised learning of discrete representations from continuous ones in neural networks (NNs) is the cornerstone of several applications today. Vector Quantisation (VQ) has become a popular method to achieve such representations, in particular in the context of generative models such as Variational Auto-Encoders (VAEs). For example, the exponential moving average-based VQ (EMA-VQ) algorithm is often used. Here we study an alternative VQ algorithm based on the learning rule of Kohonen Self-Organising Maps (KSOMs; 1982) of which EMA-VQ is a special case. In fact, KSOM is a classic VQ algorithm which is known to offer two potential benefits over the latter: empirically, KSOM is known to perform faster VQ, and discrete representations learned by KSOM form a topological structure on the grid whose nodes are the discrete symbols, resulting in an artificial version of the topographic map in the brain. We revisit these properties by using KSOM in VQ-VAEs for image processing. In particular, our experiments show that, while the speed-up compared to well-configured EMA-VQ is only observable at the beginning of training, KSOM is generally much more robust than EMA-VQ, e.g., w.r.t. the choice of initialisation schemes. Our code is public.


title: Analyzing the Engagement of Social Relationships During Life Event Shocks in Social Media

abstract: Individuals experiencing unexpected distressing events, shocks, often rely on their social network for support. While prior work has shown how social networks respond to shocks, these studies usually treat all ties equally, despite differences in the support provided by different social relationships. Here, we conduct a computational analysis on Twitter that examines how responses to online shocks differ by the relationship type of a user dyad. We introduce a new dataset of over 13K instances of individuals' self-reporting shock events on Twitter and construct networks of relationship-labeled dyadic interactions around these events. By examining behaviors across 110K replies to shocked users in a pseudo-causal analysis, we demonstrate relationship-specific patterns in response levels and topic shifts. We also show that while well-established social dimensions of closeness such as tie strength and structural embeddedness contribute to shock responsiveness, the degree of impact is highly dependent on relationship and shock types. Our findings indicate that social relationships contain highly distinctive characteristics in network interactions and that relationship-specific behaviors in online shock responses are unique from those of offline settings.


title: Hyperbolic Axisymmetric Shallow Water Moment Equations

abstract: Models for shallow water flow often assume that the lateral velocity is constant over the water height. The recently derived shallow water moment equations are an extension of these standard shallow water equations. The extended models allow for vertical changes in the lateral velocity, resulting in a system that is more accurate in situations where the horizontal velocity varies considerably over the height of the fluid. Unfortunately, already the one-dimensional models lack global hyperbolicity, an important property of partial differential equations that ensures that disturbances have a finite propagation speed. In this paper we show that the loss of hyperbolicity also occurs in two-dimensional axisymmetric systems. First, a cylindrical moment model is obtained by starting from the cylindrical incompressible Navier-Stokes equations. We derive two-dimensional axisymmetric Shallow Water Moment Equations by imposing axisymmetry in the cylindrical model. The loss of hyperbolicity is observed by directly evaluating the propagation speeds and plotting the hyperbolicity region. A hyperbolic axisymmetric moment model is then obtained by modifying the system matrix in analogy to the one-dimensional case, for which the hyperbolicity problem has already been observed and overcome. The new model is written in analytical form and we prove that hyperbolicity can be guaranteed. To test the new model, numerical simulations with both discontinuous and continuous initial data are performed. We show that the hyperbolic model leads to less oscillations than the non-hyperbolic model in the case of a discontinuous initial height profile. The hyperbolic model is preferred over the non-hyperbolic model in this specific scenario for shorter times. Further, we observe that the error decreases when the order of the model is increased in a test case with smooth initial data.


title: AI Security Threats against Pervasive Robotic Systems: A Course for Next Generation Cybersecurity Workforce

abstract: Robotics, automation, and related Artificial Intelligence (AI) systems have become pervasive bringing in concerns related to security, safety, accuracy, and trust. With growing dependency on physical robots that work in close proximity to humans, the security of these systems is becoming increasingly important to prevent cyber-attacks that could lead to privacy invasion, critical operations sabotage, and bodily harm. The current shortfall of professionals who can defend such systems demands development and integration of such a curriculum. This course description includes details about seven self-contained and adaptive modules on "AI security threats against pervasive robotic systems". Topics include: 1) Introduction, examples of attacks, and motivation; 2) - Robotic AI attack surfaces and penetration testing; 3) - Attack patterns and security strategies for input sensors; 4) - Training attacks and associated security strategies; 5) - Inference attacks and associated security strategies; 6) - Actuator attacks and associated security strategies; and 7) - Ethics of AI, robotics, and cybersecurity.


title: ATLAS-TEIDE: The next generations of ATLAS units for the Teide Observatory

abstract: In this work we present the design of the ATLAS unit (Asteroid Terrestrial-impact Last Alert System) that will be installed at Teide Observatory in Tenerife island (Spain). ATLAS-Teide will be built by the Instituto de Astrofisica de Canarias (IAC) and will be operated as part of the ATLAS network in the framework of an operation and science exploitation agreement between the IAC and the ATLAS team at University of Hawaii. ATLAS-Teide will be the first ATLAS unit based on commercial on the shelf (COTS) components. Its design is modular, each module (building block) consist of four Celestron RASA 11 telescopes that point to the same sky field, equipped with QHY600PRO CMOS cameras on an equatorial Direct Drive mount. Each module is equivalent to a 56cm effective diameter telescope and provides a 7.3 deg^2 field of view and a 1.26 arcsec/pix plate scale. ATLAS-Teide will consist of four ATLAS modules in a roll-off roof building. This configuration allows to cover the same sky area of the actual ATLAS telescopes. The first ATLAS module was installed in November 2022 in an existing clamshell at the TO. This module (ATLAS-P) is being used as a prototype to test the system capabilities, develop the needed software (control, image processing, etc.) and complete the fully integration of ATLAS-Teide in the ATLAS network. The preliminary results of the tests are presented here, and the benefits of the new ATLAS design are discussed.


title: Mathematical models for order of mutation problem in myeloproliferative neoplasm: non-additivity and non-commutativity

abstract: In some patients of myeloproliferative neoplasm, two genetic mutations can be found: JAK2 V617F and TET2. When one mutation is present or not, the other mutation has different effects on regulating gene expressions. Besides, when both mutations are present, the order of occurrence might make a difference. In this paper, we build nonlinear ordinary differential equation models and Markov chain models to explain such phenomena.


title: Tight Auditing of Differentially Private Machine Learning

abstract: Auditing mechanisms for differential privacy use probabilistic means to empirically estimate the privacy level of an algorithm. For private machine learning, existing auditing mechanisms are tight: the empirical privacy estimate (nearly) matches the algorithm's provable privacy guarantee. But these auditing techniques suffer from two limitations. First, they only give tight estimates under implausible worst-case assumptions (e.g., a fully adversarial dataset). Second, they require thousands or millions of training runs to produce non-trivial statistical estimates of the privacy leakage. This work addresses both issues. We design an improved auditing scheme that yields tight privacy estimates for natural (not adversarially crafted) datasets -- if the adversary can see all model updates during training. Prior auditing works rely on the same assumption, which is permitted under the standard differential privacy threat model. This threat model is also applicable, e.g., in federated learning settings. Moreover, our auditing scheme requires only two training runs (instead of thousands) to produce tight privacy estimates, by adapting recent advances in tight composition theorems for differential privacy. We demonstrate the utility of our improved auditing schemes by surfacing implementation bugs in private machine learning code that eluded prior auditing techniques.


title: ColibriES: A Milliwatts RISC-V Based Embedded System Leveraging Neuromorphic and Neural Networks Hardware Accelerators for Low-Latency Closed-loop Control Applications

abstract: End-to-end event-based computation has the potential to push the envelope in latency and energy efficiency for edge AI applications. Unfortunately, event-based sensors (e.g., DVS cameras) and neuromorphic spike-based processors (e.g., Loihi) have been designed in a decoupled fashion, thereby missing major streamlining opportunities. This paper presents ColibriES, the first-ever neuromorphic hardware embedded system platform with dedicated event-sensor interfaces and full processing pipelines. ColibriES includes event and frame interfaces and data processing, aiming at efficient and long-life embedded systems in edge scenarios. ColibriES is based on the Kraken system-on-chip and contains a heterogeneous parallel ultra-low power (PULP) processor, frame-based and event-based camera interfaces, and two hardware accelerators for the computation of both event-based spiking neural networks and frame-based ternary convolutional neural networks. This paper explores and accurately evaluates the performance of event data processing on the example of gesture recognition on ColibriES, as the first step of full-system evaluation. In our experiments, we demonstrate a chip energy consumption of 7.7 \si{\milli\joule} and latency of 164.5 \si{\milli\second} of each inference with the DVS Gesture event data set as an example for closed-loop data processing, showcasing the potential of ColibriES for battery-powered applications such as wearable devices and UAVs that require low-latency closed-loop control.


title: Meta-Reinforcement Learning via Exploratory Task Clustering

abstract: Meta-reinforcement learning (meta-RL) aims to quickly solve new tasks by leveraging knowledge from prior tasks. However, previous studies often assume a single mode homogeneous task distribution, ignoring possible structured heterogeneity among tasks. Leveraging such structures can better facilitate knowledge sharing among related tasks and thus improve sample efficiency. In this paper, we explore the structured heterogeneity among tasks via clustering to improve meta-RL. We develop a dedicated exploratory policy to discover task structures via divide-and-conquer. The knowledge of the identified clusters helps to narrow the search space of task-specific information, leading to more sample efficient policy adaptation. Experiments on various MuJoCo tasks showed the proposed method can unravel cluster structures effectively in both rewards and state dynamics, proving strong advantages against a set of state-of-the-art baselines.


title: Optimal Distributed Voltage Control via Primal Dual Gradient Dynamics

abstract: The rapidly increasing penetration of inverter-based resources into a power transmission network requires more sophisticated voltage control strategies considering their inherent output variabilities. In addition, faults and load variations affect the voltage profile over the power network. This paper proposes a Primal Dual Gradient Dynamics based optimal distributed voltage control approach that optimizes outputs of distributed reactive power sources to maintain an acceptable voltage profile while preserving operational limits. Case studies of this new approach on IEEE test systems have verified its effectiveness.


title: Learning to Substitute Ingredients in Recipes

abstract: Recipe personalization through ingredient substitution has the potential to help people meet their dietary needs and preferences, avoid potential allergens, and ease culinary exploration in everyone's kitchen. To address ingredient substitution, we build a benchmark, composed of a dataset of substitution pairs with standardized splits, evaluation metrics, and baselines. We further introduce Graph-based Ingredient Substitution Module (GISMo), a novel model that leverages the context of a recipe as well as generic ingredient relational information encoded within a graph to rank plausible substitutions. We show through comprehensive experimental validation that GISMo surpasses the best performing baseline by a large margin in terms of mean reciprocal rank. Finally, we highlight the benefits of GISMo by integrating it in an improved image-to-recipe generation pipeline, enabling recipe personalization through user intervention. Quantitative and qualitative results show the efficacy of our proposed system, paving the road towards truly personalized cooking and tasting experiences.


title: Functional Connectivity Dynamics show Resting-State Instability and Rightward Parietal Dysfunction in ADHD

abstract: Attention Deficit/Hyperactivity Disorder (ADHD) is one of the most common neurodevelopmental disorders in children and is characterised by inattention, impulsiveness and hyperactivity. While several studies have analysed the static functional connectivity in the resting-state functional MRI (rs-fMRI) of ADHD patients, detailed investigations are required to characterize the connectivity dynamics in the brain. In an attempt to establish a link between attention instability and the dynamic properties of Functional Connectivity (FC), we investigated the differences in temporal variability of FC between 40 children with ADHD and 40 Typically Developing (TD) children. Using a sliding-window method to segment the rs-fMRI scans in time, we employed seed-to-voxel correlation analysis for each window to obtain time-evolving seed connectivity maps for seeds placed in the posterior cingulate cortex (PCC) and the medial prefrontal cortex (mPFC). For each subject, the standard deviation of the voxel connectivity time series was used as a measure of the temporal variability of FC. Results showed that ADHD patients exhibited significantly higher variability in dFC than TD children in the cingulo-temporal, cingulo-parietal, fronto-temporal, and fronto-parietal networks ($p_{FWE} < 0.05$). Atypical temporal variability was observed in the left and right temporal gyri, the anterior cingulate cortex, and lateral regions of the right parietal cortex. The observations are consistent with visual attention issues, executive control deficit, and rightward parietal dysfunction reported in ADHD, respectively. These results help in understanding the disorder with a fresh perspective linking behavioural inattention with instability in FC in the brain.


title: \texttt{NuHamil}: A numerical code to generate nuclear two- and three-body matrix elements from chiral effective field theory

abstract: The applicability of nuclear {\it ab initio} calculations has rapidly extended over the past decades. However, starting research projects is still challenging due to the required numerical expertise in the generation of underlying nuclear interaction matrix elements and many-body calculations. To ease the first issue, in this paper we introduce the numerical code \texttt{NuHamil} to generate the nucleon-nucleon (NN) and three-nucleon (3N) matrix elements expressed in a spherical harmonic-oscillator basis, inputs of many-body calculations. The ground-state energies for the selected doubly closed shell nuclei are calculated with the no-core shell-model (NCSM) and in-medium similarity renormalization group (IMSRG). The code is written in modern Fortran, and OpenMP+MPI hybrid parallelization is available for the 3N matrix-element calculations.


title: WALLABY Pilot Survey: Hydra Cluster Galaxies UV and HI morphometrics

abstract: Galaxy morphology in atomic hydrogen (HI) and in the ultra-violet (UV) are closely linked. This has motivated their combined use to quantify morphology over the full H i disk for both H i and UV imaging. We apply galaxy morphometrics: Concentration, Asymmetry, Gini, M20 and Multimode-Intensity-Deviation statistics to the first moment-0 maps of the WALLABY survey of galaxies in the Hydra cluster center. Taking advantage of this new HI survey, we apply the same morphometrics over the full HI extent on archival GALEX FUV and NUV data to explore how well HI truncated, extended ultraviolet disk (XUV) and other morphological phenomena can be captured using pipeline WALLABY data products. Extended HI and UV disks can be identified relatively straightforward from their respective concentration. Combined with WALLABY HI, even the shallowest GALEX data is sufficient to identify XUV disks. Our second goal is to isolate galaxies undergoing ram-pressure stripping in the H i morphometric space. We employ four different machine learning techniques, a decision tree, a k-nearest neighbour, a support-vector machine, and a random forest. Up to 80% precision and recall are possible with the Random Forest giving the most robust results.


title: On Rank Energy Statistics via Optimal Transport: Continuity, Convergence, and Change Point Detection

abstract: This paper considers the use of recently proposed optimal transport-based multivariate test statistics, namely rank energy and its variant the soft rank energy derived from entropically regularized optimal transport, for the unsupervised nonparametric change point detection (CPD) problem. We show that the soft rank energy enjoys both fast rates of statistical convergence and robust continuity properties which lead to strong performance on real datasets. Our theoretical analyses remove the need for resampling and out-of-sample extensions previously required to obtain such rates. In contrast the rank energy suffers from the curse of dimensionality in statistical estimation and moreover can signal a change point from arbitrarily small perturbations, which leads to a high rate of false alarms in CPD. Additionally, under mild regularity conditions, we quantify the discrepancy between soft rank energy and rank energy in terms of the regularization parameter. Finally, we show our approach performs favorably in numerical experiments compared to several other optimal transport-based methods as well as maximum mean discrepancy.


title: A property of $\gamma$ curves derived from topology of relative trisected 4-manifolds

abstract: Any two $(g; 0, 0, k)$-trisected closed 4-manifolds are related to each other by cutting and regluing in regard to $\gamma$ curves by an element of the Torelli group. Moreover, Lambert-Cole showed that this element can be replaced by an element of the Johnson kernel. We consider similar situation in the relative case. Furthermore, we relate this situation to cork twist.


title: The qudit Pauli group: non-commuting pairs, non-commuting sets, and structure theorems

abstract: Qudits with local dimension $d>2$ can have unique structure and uses that qubits ($d=2$) cannot. Qudit Pauli operators provide a very useful basis of the space of qudit states and operators. We study the structure of the qudit Pauli group for any, including composite, $d$ in several ways. To cover composite values of $d$, we work with modules over commutative rings, which generalize the notion of vector spaces over fields. For any specified set of commutation relations, we construct a set of qudit Paulis satisfying those relations. We also study the maximum size of sets of Paulis that mutually non-commute and sets that non-commute in pairs. Finally, we give methods to find near minimal generating sets of Pauli subgroups, calculate the sizes of Pauli subgroups, and find bases of logical operators for qudit stabilizer codes. Useful tools in this study are normal forms from linear algebra over commutative rings, including the Smith normal form, alternating Smith normal form, and Howell normal form of matrices. Possible applications of this work include the construction and analysis of qudit stabilizer codes, entanglement assisted codes, parafermion codes, and fermionic Hamiltonian simulation.


title: Self-supervised Registration and Segmentation of the Ossicles with A Single Ground Truth Label

abstract: AI-assisted surgeries have drawn the attention of the medical image research community due to their real-world impact on improving surgery success rates. For image-guided surgeries, such as Cochlear Implants (CIs), accurate object segmentation can provide useful information for surgeons before an operation. Recently published image segmentation methods that leverage machine learning usually rely on a large number of manually predefined ground truth labels. However, it is a laborious and time-consuming task to prepare the dataset. This paper presents a novel technique using a self-supervised 3D-UNet that produces a dense deformation field between an atlas and a target image that can be used for atlas-based segmentation of the ossicles. Our results show that our method outperforms traditional image segmentation methods and generates a more accurate boundary around the ossicles based on Dice similarity coefficient and point-to-point error comparison. The mean Dice coefficient is improved by 8.51% with our proposed method.