title: Teaching Unknown Objects by Leveraging Human Gaze and Augmented Reality in Human-Robot Interaction

abstract: Robots are becoming increasingly popular in a wide range of environments due to their exceptional work capacity, precision, efficiency, and scalability. This development has been further encouraged by advances in Artificial Intelligence, particularly Machine Learning. By employing sophisticated neural networks, robots are given the ability to detect and interact with objects in their vicinity. However, a significant drawback arises from the underlying dependency on extensive datasets and the availability of substantial amounts of training data for these object detection models. This issue becomes particularly problematic when the specific deployment location of the robot and the surroundings, are not known in advance. The vast and ever-expanding array of objects makes it virtually impossible to comprehensively cover the entire spectrum of existing objects using preexisting datasets alone. The goal of this dissertation was to teach a robot unknown objects in the context of Human-Robot Interaction (HRI) in order to liberate it from its data dependency, unleashing it from predefined scenarios. In this context, the combination of eye tracking and Augmented Reality created a powerful synergy that empowered the human teacher to communicate with the robot and effortlessly point out objects by means of human gaze. This holistic approach led to the development of a multimodal HRI system that enabled the robot to identify and visually segment the Objects of Interest in 3D space. Through the class information provided by the human, the robot was able to learn the objects and redetect them at a later stage. Due to the knowledge gained from this HRI based teaching, the robot's object detection capabilities exhibited comparable performance to state-of-the-art object detectors trained on extensive datasets, without being restricted to predefined classes, showcasing its versatility and adaptability.


title: Grade restriction and D-brane transport for a non-abelian GLSM of an elliptic curve

abstract: We discuss a simple model for D-brane transport in non-abelian GLSMs. The model is the elliptic curve version of a non-abelian GLSM introduced by Hori and Tong and has gauge group U(2). It has two geometric phases, both of which describe the same elliptic curve, once realised as a codimension five complete intersection in G(2,5) and once as a determinantal variety. The determinantal phase is strongly coupled with unbroken SU(2). There are two singular points in the moduli space where the theory has a Coulomb branch. Using grade restriction rules, we show how to transport B-branes between the two phases along paths avoiding the singular points. With the help of the GLSM hemisphere partition function we compute analytic continuation matrices and monodromy matrices, confirming results obtained by different methods.


title: Multi Armed Bandit based Resource Allocation in Near Memory Processing Architectures

abstract: Recent advances in 3D fabrication have allowed handling the memory bottlenecks for modern data-intensive applications by bringing the computation closer to the memory, enabling Near Memory Processing (NMP). Memory Centric Networks (MCN) are advanced memory architectures that use NMP architectures, where multiple stacks of the 3D memory units are equipped with simple processing cores, allowing numerous threads to execute concurrently. The performance of the NMP is crucially dependent upon the efficient task offloading and task-to-NMP allocation. Our work presents a multi-armed bandit (MAB) based approach in formulating an efficient resource allocation strategy for MCN. Most existing literature concentrates only on one application domain and optimizing only one metric, i.e., either execution time or power. However, our solution is more generic and can be applied to diverse application domains. In our approach, we deploy Upper Confidence Bound (UCB) policy to collect rewards and eventually use it for regret optimization. We study the following metrics: instructions per cycle, execution times, NMP core cache misses, packet latencies, and power consumption. Our study covers various applications from PARSEC and SPLASH2 benchmarks suite. The evaluation shows that the system's performance improves by ~11% on average and an average reduction in total power consumption by ~12%.


title: Understanding Crypto-Ransomware

abstract: Crypto-Ransomware has been increasing in sophistication since it first appeared in September 2013, leveraging new attack vectors, incorporating advanced encryption algorithms, and expanding the number of file types it targets. In this report, we dissect nearly 30 samples of ransomware variants that have been encountered since September 2013, revealing a trend of increasing sophistication.


title: An example related to Whitney's extension problem for $L^{2,p}(\mathbb{R}^2)$ when $1<p<2$

abstract: In this paper, we prove the existence of a bounded linear extension operator $T: L^{2,p}(E) \rightarrow L^{2,p}(\mathbb{R}^2)$ when $1<p<2$, where $E \subset \mathbb{R}^2$ is a certain discrete set with fractal structure. Our proof makes use of a theorem of Fefferman-Klartag on the existence of linear extension operators for radially symmetric binary trees.


title: Dynamics and fluctuations of minimally-structured glass formers

abstract: The mean-field theory (MFT) of simple glasses, which is exact in the limit of infinite spatial dimensions, $d\rightarrow\infty$, offers theoretical insight as well as quantitative predictions about certain features of $d=3$ systems. In order to more systematically relate the behavior of physical systems to MFT, however, various finite-$d$ corrections need to be accounted for. Although some efforts along this direction have already been undertaken, theoretical and technical challenges hinder progress. A general approach to sidestep many of these difficulties consists of simulating minimally-structured models whose behavior smoothly converges to that described by the MFT as $d$ increases, so as to permit a dimensional extrapolation. Using this approach, we here extract the small fluctuations around the dynamical MFT captured by a standard liquid-state observable the non-Gaussian parameter $\alpha_2$. The results provide insight into the physical origin of these fluctuations as well as a quantitative reference with which to compare observations in more realistic glass formers.


title: Approximation algorithms for k-median problems on complex networks: theory and practice

abstract: Finding the k-medianin a network involves identifying a subset of k vertices that minimize the total distance to all other vertices in a graph. This problem has been extensively studied in computer science, graph theory, operations research, and numerous areas due to its significance in a wide range of applications. While known to be computationally challenging (NP-hard) several approximation algorithms have been proposed, most with high-order polynomial-time complexity. However, the graph topology of complex networks with heavy-tailed degree distributions present characteristics that can be exploited to yield custom-tailored algorithms. We compare eight algorithms specifically designed for complex networks and evaluate their performance based on accuracy and efficiency for problems of varying sizes and application areas. Rather than relying on a small number of problems, we conduct over 16,000 experiments covering a wide range of network sizes and k-median{} values. While individual results vary, a few methods provide consistently good results. We draw general conclusions about how algorithms perform in practice and provide general guidelines for solutions.


title: PolarVis: Towards Web-based Polarimetric Analysis

abstract: Astronomers performing polarimetric analysis on astronomical images often have to manually identify locations on their objects of interest, such as galaxies, which exhibit the influence of magnetic forces due to interaction with their environments or inherent processes. These locations are known as Lines of Sight (LoS). Analysing the various lines of sight can provide insight into the electromagnetic nature of the astrophysical object in question and its surroundings. For each LoS, astronomers generate diagnostic plots to map out the variation of the corresponding electromagnetic field, such as those of fractional polarisation and Faraday spectra. However, associating the different LoS diagnostic plots to their positions on an astronomical image requires alternating between the plots and the images. As a result, determining whether the location of the LoS influences its magnetic field variation by analysing its diagnostic plots becomes arduous due to the absence of a direct way of linking the two. PolarVis is an effort towards allowing an almost instant view of the interactive diagnostic plots corresponding to a given line of sight at the click of a button on that line of sight on the image, using an interactive web-based FITS viewer -- JS9.


title: A covariant regulator for entanglement entropy: proofs of the Bekenstein bound and QNEC

abstract: While von Neumann entropies for subregions in quantum field theory universally contain ultraviolet divergences, differences between von Neumann entropies are finite and well-defined in many physically relevant scenarios. We demonstrate that such a notion of entropy differences can be rigorously defined in quantum field theory in a general curved spacetime by introducing a novel, covariant regulator for the entropy based on the modular crossed product. This regulator associates a type II von Neumann algebra to each spacetime subregion, resulting in well-defined renormalized entropies. This prescription reproduces formulas for entropy differences that coincide with heuristic formulas widely used in the literature, and we prove that it satisfies desirable properties such as unitary invariance and concavity. As an application, we provide proofs of the Bekenstein bound and the quantum null energy condition, formulated directly in terms of vacuum-subtracted von Neumann entropies.


title: Is there a black hole in the center of the Sun?

abstract: There is probably not a black hole in the center of the sun. Despite this detail, our goal in this work to convince the reader that this question is interesting and that work studying stars with central black holes is well motivated. If primordial black holes exist then they may exist in sufficiently large numbers to explain the dark matter in the universe. While primordial black holes may form at almost any mass, the asteroid-mass window between $10^{-16} - 10^{-10} ~ \textrm{M}_\odot$ remains a viable dark matter candidate and these black holes could be captured by stars upon formation. Such a star, partially powered by accretion luminosity from a microscopic black hole in its core, has been called a `Hawking star.' Stellar evolution of Hawking stars is highly nontrivial and requires detailed stellar evolution models, which were developed in our recent work. We present here full evolutionary models of solar mass Hawking stars using two accretion schemes: one with a constant radiative efficiency, and one that is new in this work that uses an adaptive radiative efficiency to model the effects of photon trapping.

categories: astro-ph.SR astro-ph.CO astro-ph.HE gr-qc


title: Outcomes of Sub-Neptune Collisions

abstract: Observed high multiplicity planetary systems are often tightly packed. Numerical studies indicate that such systems are susceptible to dynamical instabilities. Dynamical instabilities in close-in tightly packed systems, similar to those found in abundance by Kepler, often lead to planet-planet collisions. For sub-Neptunes, the dominant type of observed exoplanets, the planetary mass is concentrated in a rocky core, but the volume is dominated by a low-density gaseous envelope. For these, using the traditional `sticky-sphere' assumption is questionable. Using both N-body integration and smoothed-particle hydrodynamics, we have simulated sub-Neptune collisions for a wide range in realistic kinematic properties such as impact parameters ($b^{\prime}$) and impact velocities ($v_{\rm{im}}$) to study the possible outcomes in detail. We find that the majority ($\sim 76\%$) of the collisions with kinematic properties similar to what is expected from dynamical instabilities in multiplanet systems may not lead to mergers of sub-Neptunes. Instead, the sub-Neptunes separate from each other, often with significant atmosphere loss. When mergers do occur, they can involve significant mass loss even from the core and can sometimes lead to complete disruption of one or both planets. Sub-Neptunes merge or disrupt if $b^{\prime}<b_{\rm{c}}$, a critical value dependent on $ v_{\rm{im}}/v_{\rm{esc}}$, where $v_{\rm{esc}}$ is the escape velocity from the surface of the hypothetical merged planet assuming sticky-sphere. For $ v_{\rm{im}}/v_{\rm{esc}}\lesssim2.5$, $b_{\rm{c}}\propto (v_{\rm{im}}/v_{\rm{esc}})^{-2}$, and collisions with $b^{\prime}<b_{\rm{c}}$ typically leads to mergers. On the other hand, for $ v_{\rm{im}}/v_{\rm{esc}}\gtrsim2.5$, $b_{\rm{c}}\propto v_{\rm{im}}/v_{\rm{esc}}$, and the collisions with $b^{\prime}<b_{\rm{c}}$ can result in complete destruction of one or both sub-Neptunes.


title: Identifying Energy-Dependent Flavor Transitions in High-Energy Astrophysical Neutrino Measurements

abstract: The flavor composition of TeV--PeV astrophysical neutrinos, i.e., the proportion of neutrinos of different flavors in their flux, is a versatile probe of high-energy astrophysics and fundamental physics. Because flavor identification is challenging and the number of detected high-energy astrophysical neutrinos is limited, so far measurements of the flavor composition have represented an average over the range of observed neutrino energies. Yet, this washes out the potential existence of changes in the flavor composition with energy and weakens our sensitivity to the many models that posit them. For the first time, we measure the energy dependence of the flavor composition, looking for a transition from low to high energies. Our present-day measurements, based on the 7.5-year public sample of IceCube High-Energy Starting Events (HESE), find no evidence of a flavor transition. The observation of HESE and through-going muons jointly by next-generation neutrino telescopes Baikal-GVD, IceCube-Gen2, KM3NeT, P-ONE, TAMBO, and TRIDENT may identify a flavor transition around 200TeV by 2030. By 2040, we could infer the flavor composition with which neutrinos are produced with enough precision to establish the transition from neutrino production via the full pion decay chain at low energies to muon-damped pion decay at high energies.


title: The Quality/Cosmology Tension for a Post-Inflation QCD Axion

abstract: It is difficult to construct a post-inflation QCD axion model that solves the axion quality problem (and hence the Strong CP problem) without introducing a cosmological disaster. In a post-inflation axion model, the axion field value is randomized during the Peccei-Quinn phase transition, and axion domain walls form at the QCD phase transition. We emphasize that the gauge equivalence of all minima of the axion potential (i.e., domain wall number one) is insufficient to solve the cosmological domain wall problem. The axion string on which a domain wall ends must exist as an individual object (as opposed to a multi-string state), and it must be produced in the early universe. These conditions are often not satisfied in concrete models. Post-inflation axion models also face a potential problem from fractionally charged relics; solving this problem often leads to low-energy Landau poles for Standard Model gauge couplings, reintroducing the quality problem. We study several examples, finding that models that solve the quality problem face cosmological problems, and vice versa. This is not a no-go theorem; nonetheless, we argue that it is much more difficult than generally appreciated to find a viable post-inflation QCD axion model. Successful examples may have a nonstandard cosmological history (e.g., multiple types of cosmic axion strings of different tensions), undermining the widespread expectation that the post-inflation QCD axion scenario predicts a unique mass for axion dark matter.


title: Bounds on EFT's in an expanding Universe

abstract: We find bounds on the Wilson coefficients of effective field theories (EFTs) living in a Universe undergoing expansion by requiring that its modes do not propagate further than a minimally coupled photon by a resolvable amount. To do so, we compute the spatial shift suffered by the EFT modes at a fixed time slice within the WKB approximation and the regime of validity of the EFT. We analyze the bounds arising on shift-symmetric scalars and curved space generalizations of Galileons.


title: Multiparticle Factorization and the Rigidity of String Theory

abstract: Is string theory uniquely determined by self-consistency? Causality and unitarity seemingly permit a multitude of putative deformations, at least at the level of two-to-two scattering. Motivated by this question, we initiate a systematic exploration of the constraints on scattering from higher-point factorization, which imposes extraordinarily restrictive sum rules on the residues and spectra defined by a given amplitude. These bounds handily exclude several proposed deformations of the string: the simplest "bespoke" amplitudes with tunable masses and a family of modified string integrands from "binary geometry." While the string itself passes all tests, our formalism directly extracts the three-point amplitudes for the low-lying string modes without the aid of worldsheet vertex operators.


title: Van-Hove tuning of Fermi surface instabilities through compensated metallicity

abstract: Van-Hove (vH) singularities in the vicinity of the Fermi level facilitate the emergence of electronically mediated Fermi surface instabilities. This is because they provide a momentum-localized enhancement of density of states enhancing selective electronic scattering channels. High-temperature topological superconductivity has been argued for in graphene at vH filling which, however, has so far proven inaccessible due to the demanded large doping from pristine half filling. We propose compensated metallicity as a path to unlock vH-driven pairing close to half filling in an electronic honeycomb lattice model. It is enabled through the strong breaking of chiral symmetry from intra-sublattice hybridization, leading to the emergence of a hole pocket (hp) nearby the van-Hove points $M$ at the Brillouin zone boundary and an electron pocket (ep) around the zone center $\Gamma$. While the ep is radially symmetric and barely contributing to the electronic ordering selection, the hp is dominated by its vH signature and yields electronic order at elevated scales.


title: Quantum Simulation of Realistic Materials in First Quantization Using Non-local Pseudopotentials

abstract: This paper improves and demonstrates the usefulness of the first quantized plane-wave algorithms for the quantum simulation of electronic structure, developed by Babbush et al. and Su et al. We describe the first quantum algorithm for first quantized simulation that accurately includes pseudopotentials. We focus on the Goedecker-Tetter-Hutter (GTH) pseudopotential, which is among the most accurate and widely used norm-conserving pseudopotentials enabling the removal of core electrons from the simulation. The resultant screened nuclear potential regularizes cusps in the electronic wavefunction so that orders of magnitude fewer plane waves are required for a chemically accurate basis. Despite the complicated form of the GTH pseudopotential, we are able to block encode the associated operator without significantly increasing the overall cost of quantum simulation. This is surprising since simulating the nuclear potential is much simpler without pseudopotentials, yet is still the bottleneck. We also generalize prior methods to enable the simulation of materials with non-cubic unit cells, which requires nontrivial modifications. Finally, we combine these techniques to estimate the block-encoding costs for commercially relevant instances of heterogeneous catalysis (e.g. carbon monoxide adsorption on transition metals) and compare to the quantum resources needed to simulate materials in second quantization. We conclude that for computational cells with many particles, first quantization often requires meaningfully less spacetime volume.


title: Bjorken $x$ weighted Energy-Energy Correlators from the Target Fragmentation Region to the Current Fragmentation Region

abstract: We present the complete spectrum for the Bjorken $x$ weighted Energy-Energy Correlation in the deep inelastic scattering (DIS) process, from the target fragmentation region to the current fragmentation region, in the Breit frame. The corresponding collinear and transverse momentum-dependent logarithms are resummed to all orders with the accuracy of NLL and N$^3$LL, respectively. And the results in the full region are matched with ${\cal O}(\alpha^2_s)$ fixed-order calculation. The final numerical predictions are presented for both EIC and CEBAF kinematics.


title: Optimization of an Optical Testbed for Characterization of EXCLAIM u-Spec Integrated Spectrometers

abstract: We describe a testbed to characterize the optical response of compact superconducting on-chip spectrometers in development for the Experiment for Cryogenic Large-Aperture Intensity Mapping (EXCLAIM) mission. EXCLAIM is a balloonborne far-infrared experiment to probe the CO and CII emission lines in galaxies from redshift 3.5 to the present. The spectrometer, called u-Spec, comprises a diffraction grating on a silicon chip coupled to kinetic inductance detectors (KIDs) read out via a single microwave feedline. We use a prototype spectrometer for EXCLAIM to demonstrate our ability to characterize the spectrometers spectral response using a photomixer source. We utilize an on-chip reference detector to normalize relative to spectral structure from the off-chip optics and a silicon etalon to calibrate the absolute frequency.


title: Supermassive black hole spin evolution in cosmological simulations with OpenGadget3

abstract: Mass and spin of massive black holes (BHs) at the centre of galaxies evolve due to gas accretion and mergers with other BHs. Besides affecting e.g. the evolution of relativistic jets, the BH spin determines the efficiency with which the BH radiates energy. Using cosmological, hydrodynamical simulations, we investigate the evolution of the BH spin across cosmic time and its role in controlling the joint growth of supermassive BHs and their host galaxies. We implement a sub-resolution prescription that models the BH spin, accounting for both BH coalescence and misaligned accretion through a geometrically thin, optically thick disc. We investigate how BH spin evolves in two idealised setups, in zoomed-in simulations, and in a cosmological volume. The latter simulation allows us to retrieve statistically robust results as for the evolution and distribution of BH spins as a function of BH properties. We find that BHs with $M_{\rm BH}\lesssim 2 \times 10^{7}\;{\rm M}_{\odot}$ grow through gas accretion, occurring mostly in a coherent fashion that favours spin-up. Above $M_{\rm BH}\gtrsim 2 \times 10^{7}~{\rm M}_{\odot}$ the gas angular momentum directions of subsequent accretion episodes are often uncorrelated with each other. The probability of counter-rotating accretion and hence spin-down increases with BH mass. In the latter mass regime, BH coalescence plays an important role. The spin magnitude displays a wide variety of histories, depending on the dynamical state of the gas feeding the BH and the relative contribution of mergers and gas accretion. As a result of their combined effect, we observe a broad range of values of the spin magnitude at the high-mass end. Our predictions for the distributions of BH spin and spin-dependent radiative efficiency as a function of BH mass are in very good agreement with observations.


title: The hardness of quantum spin dynamics

abstract: Recent experiments demonstrated quantum computational advantage in random circuit sampling and Gaussian boson sampling. However, it is unclear whether these experiments can lead to practical applications even after considerable research effort. On the other hand, simulating the quantum coherent dynamics of interacting spins has been considered as a potential first useful application of quantum computers, providing a possible quantum advantage. Despite evidence that simulating the dynamics of hundreds of interacting spins is challenging for classical computers, concrete proof is yet to emerge. We address this problem by proving that sampling from the output distribution generated by a wide class of quantum spin Hamiltonians is a hard problem for classical computers. Our proof is based on the Taylor series of the output probability, which contains the permanent of a matrix as a coefficient when bipartite spin interactions are considered. We devise a classical algorithm that extracts the coefficient using an oracle estimating the output probability. Since calculating the permanent is #P-hard, such an oracle does not exist unless the polynomial hierarchy collapses. With an anticoncentration conjecture, the hardness of the sampling task is also proven. Based on our proof, we estimate that an instance involving about 200 spins will be challenging for classical devices but feasible for intermediate-scale quantum computers with fault-tolerant qubits.


title: A non-perturbative second law of black hole mechanics in effective field theory

abstract: We describe a method for defining dynamical black hole entropy in gravitational effective field theories (EFTs). The entropy is constructed order by order in derivatives. For any fixed number of derivatives, the entropy satisfies a non-perturbative second law of black hole mechanics if the black hole remains within the regime of validity of EFT. In equilibrium the entropy reduces to the Wald entropy. It reduces to the entropy defined by Hollands et al in theories of vacuum gravity with up to 10 derivatives.


title: Warm Surprises from Cold Duets: N-Body Simulations with Two-Component Dark Matter

abstract: We explore extensive N-body simulations with two-component cold dark matter candidates. We delve into the temperature evolution, power spectrum, density perturbation, and maximum circular velocity functions. We find that the substantial mass difference between the two candidates and the annihilation of the heavier components to the lighter ones effectively endow the latter with warm dark matter-like behavior, taking advantage of all distinct features that warm dark matter candidates offer, without observational bounds on the warm dark matter mass. Moreover, we demonstrate that the two-component dark matter model aligns well with observational data, providing valuable insights into where and how to search for the elusive dark matter candidates in terrestrial experiments.


title: CLIP as RNN: Segment Countless Visual Concepts without Training Endeavor

abstract: Existing open-vocabulary image segmentation methods require a fine-tuning step on mask labels and/or image-text datasets. Mask labels are labor-intensive, which limits the number of categories in segmentation datasets. Consequently, the vocabulary capacity of pre-trained VLMs is severely reduced after fine-tuning. However, without fine-tuning, VLMs trained under weak image-text supervision tend to make suboptimal mask predictions. To alleviate these issues, we introduce a novel recurrent framework that progressively filters out irrelevant texts and enhances mask quality without training efforts. The recurrent unit is a two-stage segmenter built upon a frozen VLM. Thus, our model retains the VLM's broad vocabulary space and equips it with segmentation ability. Experiments show that our method outperforms not only the training-free counterparts, but also those fine-tuned with millions of data samples, and sets the new state-of-the-art records for both zero-shot semantic and referring segmentation. Concretely, we improve the current record by 28.8, 16.0, and 6.9 mIoU on Pascal VOC, COCO Object, and Pascal Context.


title: The Global Asteroseismology Project Proof of Concept: Asteroseismology of Massive Stars with Continuous Ground-Based Observations

abstract: Massive (>~ 8 solar masses) stars are the progenitors of many astrophysical systems, yet key aspects of their structure and evolution are poorly understood. Asteroseismology has the potential to solve these open puzzles, however, sampling both the short period pulsations and long period beat patterns of massive stars poses many observational challenges. Ground-based single-site observations require years or decades to discern the main oscillation modes. Multi-site campaigns were able to shorten this time span, but have not been able to scale up to population studies on samples of objects. Space-based observations can achieve both continuous sampling and observe large numbers of objects, however, most lack the multi-band data that is often necessary for mode identification and removing model degeneracies. Here, we develop and test a new ground-based observational strategy for discerning and identifying the main oscillation modes of a massive star in a few months, in a way that can be scaled to large samples. We do so using the Las Cumbres Observatory - a unique facility consisting of robotic, homogeneous telescopes operating as a global network, overcoming most of the challenges of previous multi-site efforts, but presenting new challenges which we tailor our strategy to address. This work serves as the proof of concept for the Global Asteroseismology Project, which aims to move massive star asteroseismology from single-objects to bulk studies, unleashing its full potential in constraining stellar structure and evolution models. This work also demonstrates the ability of the Las Cumbres Observatory to perform multi-site continuous observations for various science goals.


title: Interpreting the long-term variability of the changing-look AGN Mrk 1018

abstract: We present a thorough study of the Changing-Look Active Galactic Nucleus (CL-AGN) Mrk 1018, utilizing an extensive dataset spanning optical, UV, and X-ray spectro-photometric data from 2005 to 2019. We analysed X-ray spectra and broad-band photometry, and performed optical-to-X-ray spectral energy distribution (SED) fitting to comprehend the observed changing-look behaviour. We found that over the 14 years in analysis, significant changes in X-ray spectra occurred, as the hardness ratio increases by a factor of ~2. We validated also the broad-band dimming, with optical, UV, and X-ray luminosities decreasing by factors of >7, >24 and ~9, respectively. These dims are attributed to the declining UV emission. We described the X-ray spectra with a two-Comptonization model, revealing a consistent hot comptonizing medium but a cooling warm component. This cooling, linked to the weakening of the magnetic fields in the accretion disk, explains the UV dimming. We propose that the weakening is caused by the formation of a jet, in turn originated from the change of state of the inner accretion flow. Our optical-to-X-ray SED fitting supports this conclusion, as the normalised accretion rate is super-critical ($\mu=$0.06>0.02) in the bright state and sub-critical ($\mu=$0.01<0.02) in the faint state. Instabilities arising at the interface of the state-transition are able to reduce the viscous timescale to the observed ~10 years of Mrk 1018 variability. We explored a possible triggering mechanism for this state transition, involving gaseous clouds pushed onto the AGN sub-pc regions by a recent merging event or by cold chaotic accretion. This scenario, if validated by future simulations, could enhance our understanding of CL-AGN and raises questions about an accretion rate of ~0.02, coupled with minor disturbances in the accretion disk, being the primary factor in the changing-look phenomenon.


title: Orbital Torus Imaging: Acceleration, density, and dark matter in the Galactic disk measured with element abundance gradients

abstract: Under the assumption of a simple and time-invariant gravitational potential, many Galactic dynamics techniques infer the Milky Way's mass and dark matter distribution from stellar kinematic observations. These methods typically rely on parameterized potential models of the Galaxy and must take into account non-trivial survey selection effects, because they make use of the density of stars in phase space. Large-scale spectroscopic surveys now supply information beyond kinematics in the form of precise stellar label measurements (especially element abundances). These element abundances are known to correlate with orbital actions or other dynamical invariants. Here, we use the Orbital Torus Imaging (OTI) framework that uses abundance gradients in phase space to map orbits. In many cases these gradients can be measured without detailed knowledge of the selection function. We use stellar surface abundances from the APOGEE survey combined with kinematic data from the Gaia mission. Our method reveals the vertical ($z$-direction) orbit structure in the Galaxy and enables empirical measurements of the vertical acceleration field and orbital frequencies in the disk. From these measurements, we infer the total surface mass density, $\Sigma$, and midplane volume density, $\rho_0$, as a function of Galactocentric radius and height. Around the Sun, we find $\Sigma_{\odot}(z=1.1$ kpc)$=72^{+6}_{-9}$M$_{\odot}$pc$^{-2}$ and $\rho_{\odot}(z=0)=0.081^{+0.015}_{-0.009}$ M$_{\odot}$pc$^{-3}$ using the most constraining abundance ratio, [Mg/Fe]. This corresponds to a dark matter contribution in surface density of $\Sigma_{\odot,\mathrm{DM}}(z=1.1$ kpc)$=24\pm4$ M$_{\odot}$pc$^{-2}$, and in total volume mass density of $\rho_{\odot,\mathrm{DM}}(z=0)=0.011\pm0.002$ M$_{\odot}$pc$^{-3}$. Moreover, using these mass density values we estimate the scale length of the low-$\alpha$ disc to be $h_R=2.24\pm0.06$kpc.


title: Combining QED and Approximate N${}^3$LO QCD Corrections in a Global PDF Fit: MSHT20qed_an3lo PDFs

abstract: We present the MSHT20qed_an3lo parton distribution functions (PDFs). These result from the first global PDF analysis to combine QED and approximate N${}^3$LO (aN${}^3$LO) QCD corrections in the theoretical calculation of the PDF evolution and cross sections entering the fit. We examine the PDF impact, and find that the effect of QED is relatively mild in comparison to the aN${}^3$LO corrections, although it should still be accounted for at the level of precision now required. These QED corrections are in addition found to roughly factorise from the QCD corrections; that is, their relative impact on the PDFs is roughly the same at NNLO and aN${}^3$LO. The fit quality exhibits a very small deterioration at aN${}^3$LO upon the inclusion of QED corrections, which is rather smaller than the deterioration observed at NNLO in QCD. The impact on several cross sections at N${}^3$LO is also examined, including the Higgs cross section via gluon fusion at N${}^3$LO. Finally, a LO in QCD fit that includes QED corrections is also presented: the MSHT20qed_lo set.


title: Universal definition of the non-conformal trace anomaly

abstract: We show that there exists a generalized, universal notion of the trace anomaly for theories which are not conformally invariant at the classical level. The definition is suitable for any regularization scheme and clearly states to what extent the classical equations of motion should be used, thus resolving existing controversies surrounding previous proposals. Additionally, we exhibit the link between our definition of the anomaly and the functional Jacobian arising from a Weyl transformation.


title: Mirror dualities with four supercharges

abstract: We consider 3d N=2 non-abelian Hanany-Witten brane setups with chiral flavor symmetry. We propose that the associated field theories are quivers with improved bifundamentals, instead of standard bifundamentals. The improved bifundamental is a strongly coupled SCFT that carries one more U(1) global symmetry than the standard bifundamental. As a consequence, our proposal overcomes the long standing challenge of associating to each N=2 brane setup a gauge theory with the full rank global symmetry, allowing the study of all the usual supersymmetric observables, such as superconformal index, sphere partition function, chiral ring and moduli space. The construction passes many non-trivial tests, for instance we algorithmically prove that any two improved quivers associated to S-dual brane setups are infrared dual. The 3d N=2 mirror dualities can be uplifted to 4d dualities with 4d improved bifundamentals connecting USp(2N) nodes.


title: Strongly interacting photons in 2D waveguide QED

abstract: One dimensional confinement in waveguide Quantum Electrodynamics (QED) plays a crucial role to enhance light-matter interactions and to induce a strong quantum nonlinear optical response. In two or higher dimensional settings, this response is reduced since photons can be emitted within a larger phase space, opening the question whether strong photon-photon interaction can be still achieved. In this study, we positively answer this question for the case of a 2D square array of atoms coupled to the light confined into a two-dimensional waveguide. More specifically, we demonstrate the occurrence of long-lived two-photon repulsive and bound states with genuine 2D features. Furthermore, we observe signatures of these effects also in free-space atomic arrays in the form of weakly-subradiant in-band scattering resonances. Our findings provide a paradigmatic signature of the presence of strong photon-photon interactions in 2D waveguide QED.


title: GMTalker: Gaussian Mixture-based Audio-Driven Emotional talking video Portraits

abstract: Synthesizing high-fidelity and emotion-controllable talking video portraits, with audio-lip sync, vivid expressions, realistic head poses, and eye blinks, has been an important and challenging task in recent years. Most existing methods suffer in achieving personalized and precise emotion control, smooth transitions between different emotion states, and the generation of diverse motions. To tackle these challenges, we present GMTalker, a Gaussian mixture-based emotional talking portraits generation framework. Specifically, we propose a Gaussian mixture-based expression generator that can construct a continuous and disentangled latent space, achieving more flexible emotion manipulation. Furthermore, we introduce a normalizing flow-based motion generator pretrained on a large dataset with a wide-range motion to generate diverse head poses, blinks, and eyeball movements. Finally, we propose a personalized emotion-guided head generator with an emotion mapping network that can synthesize high-fidelity and faithful emotional video portraits. Both quantitative and qualitative experiments demonstrate our method outperforms previous methods in image quality, photo-realism, emotion accuracy, and motion diversity.


title: Electroweak Scattering at the Muon Shot

abstract: It has long been recognized that the scattering of electroweak particles at very high energies is dominated by vector boson fusion, which probes the origin of electroweak symmetry breaking and offers a unique window into the ultraviolet regime of the SM. Previous studies assume SM-like couplings and rely on the effective $W$ approximation (or electroweak parton distribution), whose validity is well-established within the SM but not yet studied in the presence of anomalous Higgs couplings. In this work, we critically examine the electroweak production of two Higgs bosons in the presence of anomalous $VVh$ and $VVhh$ couplings. We compute the corresponding helicity amplitudes and compare the cross section results in the effective $W$ approximation with the full fixed-order calculation. In particular, we identify two distinct classes of anomalous Higgs couplings, whose effects are not captured by vector boson fusion and effective $W$ approximation. Such very high energy electroweak scatterings can be probed at the Muon Shot, a multi-TeV muon collider upon which we base our study, although similar considerations apply to other high energy colliders.


title: Reacting like Humans: Incorporating Intrinsic Human Behaviors into NAO through Sound-Based Reactions to Fearful and Shocking Events for Enhanced Sociability

abstract: Robots' acceptability among humans and their sociability can be significantly enhanced by incorporating human-like reactions. Humans can react to environmental events very quickly and without thinking. An instance where humans show natural reactions is when they encounter a sudden and loud sound that startles or frightens them. During such moments, individuals may instinctively move their hands, turn toward the origin of the sound, and try to determine the event's cause. This inherent behavior motivated us to explore this less-studied part of social robotics. In this work, a multi-modal system composed of an action generator, sound classifier, and YOLO object detector was designed to sense the environment and, in the presence of sudden loud sounds, show natural human fear reactions; and finally, locate the fear-causing sound source in the environment. These valid generated motions and inferences could imitate intrinsic human reactions and enhance the sociability of robots. For motion generation, a model based on LSTM and MDN networks was proposed to synthesize various motions. Also, in the case of sound detection, a transfer learning model was preferred that used the spectrogram of the sound signals as its input. After developing individual models for sound detection, motion generation, and image recognition, they were integrated into a comprehensive "fear" module implemented on the NAO robot. Finally, the fear module was tested in practical application and two groups of experts and non-experts (in the robotics area) filled out a questionnaire to evaluate the performance of the robot. We indicated that the proposed module could convince the participants that the Nao robot acts and reasons like a human when a sudden and loud sound is in the robot's peripheral environment, and additionally showed that non-experts have higher expectations about social robots and their performance.


title: Topological Signal Processing on Quantum Computers for Higher-Order Network Analysis

abstract: Predicting and analyzing global behaviour of complex systems is challenging due to the intricate nature of their component interactions. Recent work has started modelling complex systems using networks endowed with multiway interactions among nodes, known as higher-order networks. In this context, simplicial complexes are a class of higher-order networks that have received significant attention due to their topological structure and connections to Hodge theory. Topological signal processing utilizes these connections to analyze and manipulate signals defined on non-Euclidean domains such as simplicial complexes. In this work, we present a general quantum algorithm for implementing filtering processes in TSP and describe its application to extracting network data based on the Hodge decomposition. We leverage pre-existing tools introduced in recent quantum algorithms for topological data analysis and combine them with spectral filtering techniques using the quantum singular value transformation framework. While this paper serves as a proof-of-concept, we obtain a super-polynomial improvement over the best known classical algorithms for TSP filtering processes, modulo some important caveats about encoding and retrieving the data from a quantum state. The proposed algorithm generalizes the applicability of tools from quantum topological data analysis to novel applications in analyzing high-dimensional complex systems.


title: A Multi-Precision Quadratic Regularization Method for Unconstrained Optimization with Rounding Error Analysis

abstract: We propose a multi-precision extension of the Quadratic Regularization (R2) algorithm that enables it to take advantage of low-precision computations, and by extension to decrease energy consumption during the solve. The lower the precision in which computations occur, the larger the errors induced in the objective value and gradient, as well as in all other computations that occur in the course of the iterations. The Multi-Precision R2 (MPR2) algorithm monitors the accumulation of rounding errors with two aims: to provide guarantees on the result of all computations, and to permit evaluations of the objective and its gradient in the lowest precision possible while preserving convergence properties. MPR2's numerical results show that single precision offers enough accuracy for more than half of objective evaluations and most of gradient evaluations during the algorithm's execution. However, MPR2 fails to converge on several problems of the test set for which double precision does not offer enough precision to ensure the convergence conditions before reaching a first order critical point. That is why we propose a practical version of MPR2 with relaxed conditions which converges for almost as many problems as R2 and potentially enables to save about 50 % time and 60% energy for objective evaluation and 50 % time and 70% energy for gradient evaluation.


title: Topology change and non-geometry at infinite distance

abstract: The distance conjecture diagnoses viable low-energy effective realisation of consistent theories of quantum gravity by examining their breakdown at infinite distance in their parameter space. At the same time, infinite distance points in parameter space are naturally intertwined with string dualities. We explore the implications of the distance conjecture when T-duality is applied to curved compact manifolds and in presence of (non-)geometric fluxes. We provide evidence of how divergent potentials signal pathological infinite distance points in the scalar field space where towers of light states cannot be sustained by the curved background. This leads us to suggest an extension to the current statement of the Swampland distance conjecture in curved spaces or in presence of non-trivial fluxes supporting the background.


title: Non-standard Wigner doublets

abstract: Guided by a conservative formulation in investigating the physical content of quantum fields, we explore non-standard Wigner classes of particles that could provide the basis for self-interaction models to dark matter. We critically contrast the analysis with long-standing constraints to non-standard Wigner classes in the literature to discuss the model's viability.


title: Extremely metal-poor stars in the Fornax and Carina dwarf spheroidal galaxies

abstract: We present our analysis of VLT/UVES and X-shooter observations of six very metal-poor stars, including four stars at [Fe/H]$\approx$$-3$ in the Fornax and Carina dwarf spheroidal (dSph) galaxies. To date, this metallicity range in these two galaxies has not yet been investigated fully, or at all in some cases. The chemical abundances of 25 elements are presented, based on 1D and local thermodynamic equilibrium (LTE) model atmospheres. We discuss the different elemental groups, and find that $\alpha$- and iron-peak elements in these two systems are generally in good agreement with the Milky Way halo at the same metallicity. Our analysis reveals that none of the six stars we studied exhibits carbon enhancement, which is noteworthy given the prevalence of carbon-enhanced metal-poor stars without s-process enhancement (CEMP-no) in the Galaxy at similarly low metallicities. Our compilation of literature data shows that the fraction of CEMP-no stars in dSph galaxies is significantly lower than in the Milky Way, and than in ultra-faint dwarf galaxies. Furthermore, we report the discovery of the lowest metallicity, [Fe/H]=$-2.92$, r-process rich (r-I) star in a dSph galaxy. This star, fnx_06_019, has $\rm[Eu/Fe]=+0.8$, and also shows enhancement of La, Nd, and Dy, $\rm[X/Fe]>+0.5$. Our new data in Carina and Fornax help populate the extremely low metallicity range in dSph galaxies, and add to the evidence of a low fraction of CEMP-no stars in these systems.


title: Schr\"odinger evolution of a scalar field in Riemannian and pseudoRiemannian expanding metrics

abstract: We study the quantum field theory (QFT) of a scalar field in the Schr\"odinger picture in the functional formulation. We derive a formula for the evolution kernel in a flat expanding metric. We discuss a transition between Riemannian and pseudoRiemannian metrics (signature inversion). We express the real time Schr\"odinger evolution by the Brownian motion . We discuss the Feynman integral for a scalar field in a radiation background. We show that the unitary Schr\"odinger evolution for positive time can go over for negative time into a dissipative evolution described by diffusive paths.


title: Dzyaloshinskii-Moriya interaction inducing weak ferromagnetism in centrosymmetric altermagnets and weak ferrimagnetism in noncentrosymmetric altermagnets

abstract: The Dzyaloshinskii-Moriya interaction (DMI) has explained successfully the weak ferromagnetism in some centrosymmetric antiferromagnets. However, in the last years, it was generally claimed that the DMI is not effective in centrosymmetric systems. We reconciled these views by separating the conventional antiferromagnets and altermagnets. Altermagnets represent collinear antiferromagnetic compounds with spin-up and spin-down sublattices connected only by mirror and roto-translational symmetries. Consequently, the system shows even-parity wave spin order in the k-space lifting the Kramer's degeneracy in the non-relativistic band structure. We emphasize that the DMI can create weak ferromagnetism in centrosymmetric altermagnets while it is not effective in centrosymmetric conventional antiferromagnets. Additionally, DMI can create weak ferromagnetism or weak ferrimagnetism in noncentrosymmetric altermagnets. Once the spin-orbit coupling is included in an altermagnetic system without time-reversal symmetry, the components of spin moments of the two sublattices along the Neel vector are antiparallel but the other two spin components orthogonal to the N\'eel vector can be either parallel or antiparallel for centrosymetric systems. For noncentrosymmetric systems, we can have different bands showing parallel or antiparallel spin components resulting in weak ferrimagnetism. We can divide the altermagnetic compounds into classes based on the weak ferromagnetism or weak ferrimagnetism properties. The weak ferromagnetism and weak ferrimagnetism induced by DMI is a property exclusively of the altermagnets, not present in either ferromagnets or conventional antiferromagnets so we propose that altermagnets should be classified based on this property.


title: Bayesian Online Learning for Consensus Prediction

abstract: Given a pre-trained classifier and multiple human experts, we investigate the task of online classification where model predictions are provided for free but querying humans incurs a cost. In this practical but under-explored setting, oracle ground truth is not available. Instead, the prediction target is defined as the consensus vote of all experts. Given that querying full consensus can be costly, we propose a general framework for online Bayesian consensus estimation, leveraging properties of the multivariate hypergeometric distribution. Based on this framework, we propose a family of methods that dynamically estimate expert consensus from partial feedback by producing a posterior over expert and model beliefs. Analyzing this posterior induces an interpretable trade-off between querying cost and classification performance. We demonstrate the efficacy of our framework against a variety of baselines on CIFAR-10H and ImageNet-16H, two large-scale crowdsourced datasets.


title: I Open at the Close: A Deep Reinforcement Learning Evaluation of Open Streets Initiatives

abstract: The open streets initiative "opens" streets to pedestrians and bicyclists by closing them to cars and trucks. The initiative, adopted by many cities across North America, increases community space in urban environments. But could open streets also make cities safer and less congested? We study this question by framing the choice of which streets to open as a reinforcement learning problem. In order to simulate the impact of opening streets, we first compare models for predicting vehicle collisions given network and temporal data. We find that a recurrent graph neural network, leveraging the graph structure and the short-term temporal dependence of the data, gives the best predictive performance. Then, with the ability to simulate collisions and traffic, we frame a reinforcement learning problem to find which streets to open. We compare the streets in the NYC Open Streets program to those proposed by a Q-learning algorithm. We find that the streets proposed by the Q-learning algorithm have reliably better outcomes, while streets in the program have similar outcomes to randomly selected streets. We present our work as a step toward principally choosing which streets to open for safer and less congested cities. All our code and data are available on Github.


title: A Note on Local Convergence of Iterative Processes for Pipe Network Analysis

abstract: Analysis of pipe networks involves computing flow rates and pressure differences on pipe segments in the network, given the external inflow/outflow values. This analysis can be conducted using iterative methods, among which the algorithms of Hardy Cross and Newton-Raphson have historically been applied in practice. In this note, we address the mathematical analysis of the local convergence of these algorithms. The loop-based Newton-Raphson algorithm converges quadratically fast, and we provide estimates for its convergence radius to correct some estimates in the previous literature. In contrast, we show that the convergence of the Hardy Cross algorithm is only linear. This provides theoretical confirmation of experimental observations reported earlier in the literature.


title: A Conditioned Unsupervised Regression Framework Attuned to the Dynamic Nature of Data Streams

abstract: In scenarios where obtaining real-time labels proves challenging, conventional approaches may result in sub-optimal performance. This paper presents an optimal strategy for streaming contexts with limited labeled data, introducing an adaptive technique for unsupervised regression. The proposed method leverages a sparse set of initial labels and introduces an innovative drift detection mechanism to enable dynamic model adaptations in response to evolving patterns in the data. To enhance adaptability, we integrate the ADWIN (ADaptive WINdowing) algorithm with error generalization based on Root Mean Square Error (RMSE). ADWIN facilitates real-time drift detection, while RMSE provides a robust measure of model prediction accuracy. This combination enables our multivariate method to effectively navigate the challenges of streaming data, continuously adapting to changing patterns while maintaining a high level of predictive precision. We evaluate the performance of our multivariate method across various public datasets, comparing it to non-adapting baselines. Through comprehensive assessments, we demonstrate the superior efficacy of our adaptive regression technique for tasks where obtaining labels in real-time is a significant challenge. The results underscore the method's capacity to outperform traditional approaches and highlight its potential in scenarios characterized by label scarcity and evolving data patterns.


title: On Rosenbaum's Rank-based Matching Estimator

abstract: In two influential contributions, Rosenbaum (2005, 2020) advocated for using the distances between component-wise ranks, instead of the original data values, to measure covariate similarity when constructing matching estimators of average treatment effects. While the intuitive benefits of using covariate ranks for matching estimation are apparent, there is no theoretical understanding of such procedures in the literature. We fill this gap by demonstrating that Rosenbaum's rank-based matching estimator, when coupled with a regression adjustment, enjoys the properties of double robustness and semiparametric efficiency without the need to enforce restrictive covariate moment assumptions. Our theoretical findings further emphasize the statistical virtues of employing ranks for estimation and inference, more broadly aligning with the insights put forth by Peter Bickel in his 2004 Rietz lecture (Bickel, 2004).


title: Dynamical friction in self-interacting ultralight dark matter

abstract: We explore how dynamical friction in an ultralight dark matter (ULDM) background is affected by dark matter self-interactions. We calculate the force of dynamical friction on a point mass moving through a uniform ULDM background with self-interactions, finding that the force of dynamical friction vanishes for sufficiently strong repulsive self-interactions. Using the pseudospectral solver $\texttt{UltraDark.jl}$, we show with simulations that reasonable values of the ULDM self-interaction strength and particle mass cause $\mathcal{O}(1)$ differences in the acceleration of an object like a supermassive black hole (SMBH) traveling near the center of a soliton, relative to the case with no self-interactions. For example, repulsive self-interactions with $\lambda = 10^{-90}$ yield a deceleration due to dynamical friction $\approx70\%$ smaller than a model with no self-interactions. We discuss the observational implications of our results for SMBHs near soliton centers and for massive satellite galaxies falling into ultralight axion halos and show that outcomes are dependent on whether a self-interaction is present or not.


title: A Perspective of Q-value Estimation on Offline-to-Online Reinforcement Learning

abstract: Offline-to-online Reinforcement Learning (O2O RL) aims to improve the performance of offline pretrained policy using only a few online samples. Built on offline RL algorithms, most O2O methods focus on the balance between RL objective and pessimism, or the utilization of offline and online samples. In this paper, from a novel perspective, we systematically study the challenges that remain in O2O RL and identify that the reason behind the slow improvement of the performance and the instability of online finetuning lies in the inaccurate Q-value estimation inherited from offline pretraining. Specifically, we demonstrate that the estimation bias and the inaccurate rank of Q-value cause a misleading signal for the policy update, making the standard offline RL algorithms, such as CQL and TD3-BC, ineffective in the online finetuning. Based on this observation, we address the problem of Q-value estimation by two techniques: (1) perturbed value update and (2) increased frequency of Q-value updates. The first technique smooths out biased Q-value estimation with sharp peaks, preventing early-stage policy exploitation of sub-optimal actions. The second one alleviates the estimation bias inherited from offline pretraining by accelerating learning. Extensive experiments on the MuJoco and Adroit environments demonstrate that the proposed method, named SO2, significantly alleviates Q-value estimation issues, and consistently improves the performance against the state-of-the-art methods by up to 83.1%.


title: Adaptive Phase Estimation with Squeezed Vacuum Approaching the Quantum Limit

abstract: Phase estimation plays a central role in communications, sensing, and information processing. Quantum correlated states, such as squeezed states, enable phase estimation beyond the shot-noise limit, and in principle approach the ultimate quantum limit in precision, when paired with optimal quantum measurements. However, physical realizations of optimal quantum measurements for optical phase estimation with quantum correlated states are still unknown. Here we address this problem by introducing an adaptive Gaussian measurement strategy for optical phase estimation with squeezed vacuum states that, by construction, approaches the quantum limit in precision. This strategy builds from a comprehensive set of locally optimal POVMs through rotations and homodyne measurements and uses the Adaptive Quantum State Estimation framework for optimizing the adaptive measurement process, which, under certain regularity conditions, guarantees asymptotic optimality for this quantum parameter estimation problem. As a result, the adaptive phase estimation strategy based on locally-optimal homodyne measurements achieves the quantum limit within the phase interval of $[0, \pi/2)$. Furthermore, we generalize this strategy by including heterodyne measurements, enabling phase estimation across the full range of phases from $[0, \pi)$, where squeezed vacuum allows for unambiguous phase encoding. Remarkably, for this phase interval, which is the maximum range of phases that can be encoded in squeezed vacuum, this estimation strategy maintains an asymptotic quantum-optimal performance, representing a significant advancement in quantum metrology.


title: Multi-compartmental modeling for Gd-EOB-DTPA using sparse human DCE-MRI data

abstract: Purpose: To derive kinetic equations for multi-compartmental contrast agent distribution from first principles and apply it to two and three compartments for Gd-EOB-DTPA using low time resolution human liver DCE-MRI data. Methods: The continuity and diffusion equation were combined and used to derive a general form for differential equations governing multi-compartmental particle exchange. They were applied to two (equivalent to the Tofts model) and three compartments. Both models were fit to human DCE-MRI data with low temporal resolution and three compartment model's parameters' implications are discussed. Results: The model derived for two compartments is shown to be equivalent with the Tofts model. Using reasonable biological and physical assumptions an analytical solution for the three compartment model is obtained. The three compartment model was able to fit all Gd-EOB-DTPA DCE-MRI data whereas the Tofts model did not. Differences were increased for cases with large as well as rapid uptake of contrast. Conclusion: We demonstrated the ability to fit sparse Gd-EOB-DTPA DCE-MRI data using a three compartment model derived from first principles whose parameters can be used to help quantify overall and regional liver function.


title: Radiosensitization beyond DNA damage: Monte Carlo simulations of realistic nanomaterial biodistributions

abstract: We investigated cellular distribution of a tumor-specific gadolinium chelate in 4T1 and U87 cancer cells with the goal to generate more realistic geometries for Monte Carlo simulations of radiation interaction with nanoparticles in cells. Cells were exposed to the agent in-vitro for 30 minutes to 72 hours before being fixed and imaged using transmission electron microscopy. Initially, electron-dense areas consistent with gadolinium were observable throughout the cytoplasm. At six hours those areas were restricted to endosomes and at 24 hours or longer electron dense areas were only found in lysosomes. Lysosomes were on average larger in 4T1 cells, which were exposed to 1 mM concentration compared to U87 cells, exposed to 1 uM. Based on this information we built an extended cell model for Monte Carlo simulations that includes lysosomes with discrete nanoparticle enclaves in addition to mitochondria and the nucleus.


title: 1-D Isentropic Euler flows: Self-similar Vacuum Solutions

abstract: We consider one-dimensional self-similar solutions to the isentropic Euler system when the initial data are at vacuum to the left of the origin. For $x>0$ the initial velocity and sound speed are of form $u_0(x)=u_+x^{1-\lambda}$ and $c_0(x)=c_+x^{1-\lambda}$, for constants $u_+\in\RR$, $c_+>0$, $\lambda\in\RR$. We analyze the resulting solutions in terms of the similarity parameter $\lambda$, the adiabatic exponent $\gamma$, and the initial (signed) Mach number $\text{Ma}=u_+/c_+$. Restricting attention to locally bounded data, we find that when the sound speed initially decays to zero in a H\"older manner ($0<\lambda<1$), the resulting flow is always defined globally. Furthermore, there are three regimes depending on $\text{Ma}$: for sufficiently large positive $\text{Ma}$-values, the solution is continuous and the initial H\"older decay is immediately replaced by $C^1$-decay to vacuum along a stationary vacuum interface; for moderate values of $\text{Ma}$, the solution is again continuous and with an accelerating vacuum interface along which $c^2$ decays linearly to zero (i.e., a "physical singularity''); for sufficiently large negative $\text{Ma}$-values, the solution contains a shock wave emanating from the initial vacuum interface and propagating into the fluid, together with a physical singularity along an accelerating vacuum interface. In contrast, when the sound speed initially decays to zero in a $C^1$ manner ($\lambda<0$), a global flow exists only for sufficiently large positive values of $\text{Ma}$. Non-existence of global solutions for smaller $\text{Ma}$-values is due to rapid growth of the data at infinity and is unrelated to the presence of a vacuum.


title: The discrete adiabatic quantum linear system solver has lower constant factors than the randomized adiabatic solver

abstract: The solution of linear systems of equations is the basis of many other quantum algorithms, and recent results provided an algorithm with optimal scaling in both the condition number $\kappa$ and the allowable error $\epsilon$ [PRX Quantum \textbf{3}, 0403003 (2022)]. That work was based on the discrete adiabatic theorem, and worked out an explicit constant factor for an upper bound on the complexity. Here we show via numerical testing on random matrices that the constant factor is in practice about 1,500 times smaller than the upper bound found numerically in the previous results. That means that this approach is far more efficient than might naively be expected from the upper bound. In particular, it is over an order of magnitude more efficient than using a randomised approach from [arXiv:2305.11352] that claimed to be more efficient.


title: Unleashed from Constrained Optimization: Quantum Computing for Quantum Chemistry Employing Generator Coordinate Method

abstract: Hybrid quantum-classical approaches offer potential solutions for quantum chemistry problems, but they also introduce challenges such as the barren plateau and the exactness of the ansatze. These challenges often manifest as constrained optimization problems without a guarantee of identifying global minima. In this work, we highlight the interconnection between constrained optimization and generalized eigenvalue problems, using a unique class of non-orthogonal and overcomplete basis sets generated by Givens rotation-type canonical transformations on a reference state. Employing the generator coordinate approach, we represent the wave function in terms of these basis sets. The ensuing generalized eigenvalue problem yields rigorous lower bounds on energy, outperforming the conventional variational quantum eigensolver (VQE) that employs the same canonical transformations in its ansatze. Our approach effectively tackles the barren plateau issue and the heuristic nature of numerical minimizers in the standard VQE, making it ideal for intricate quantum chemical challenges. For real-world applications, we propose an adaptive scheme for selecting these transformations, emphasizing the linear expansion of the non-orthogonal basis sets. This ensures a harmonious balance between accuracy and efficiency in hybrid quantum-classical simulations. Our analysis and suggested methodology further broaden the applications of quantum computing in quantum chemistry. Notably, they pave the way for alternative strategies in excited state computation and Hamiltonian downfolding, laying the groundwork for sophisticated quantum simulations in chemistry.


title: Bootstrapping Boundary QED Part I

abstract: We use the numerical conformal bootstrap to study boundary quantum electrodynamics, the theory of a four dimensional photon in a half space coupled to charged conformal matter on the boundary. This system is believed to be a boundary conformal field theory with an exactly marginal coupling corresponding to the strength of the interaction between the photon and the matter degrees of freedom. In part one of this project, we present three results. We show how the Maxwell equations put severe constraints on boundary three-point functions involving two currents and a symmetric traceless tensor. We use semi-definite programming to show that any three dimensional conformal field theory with a global U(1) symmetry must have a spin two gap less than about 1.05. Finally, combining a numerical bound on an OPE coefficient and some Ward identities involving the current and the displacement operator, we bound the displacement operator two-point function above. This upper bound also constrains a boundary contribution to the anomaly in the trace of the stress tensor for these types of theories.


title: Scaling Culture in Blockchain Gaming: Generative AI and Pseudonymous Engagement

abstract: Managing rapidly growing decentralized gaming communities brings unique challenges at the nexus of cultural economics and technology. This paper introduces a streamlined analytical framework that utilizes Large Language Models (LLMs), in this instance open-access generative pre-trained transformer (GPT) models, offering an efficient solution with deeper insights into community dynamics. The framework aids moderators in identifying pseudonymous actor intent, moderating toxic behavior, rewarding desired actions to avoid unintended consequences of blockchain-based gaming, and gauging community sentiment as communities venture into metaverse platforms and plan for hypergrowth. This framework strengthens community controls, eases onboarding, and promotes a common moral mission across communities while reducing agency costs by 95 pct. Highlighting the transformative role of generative AI, the paper emphasizes its potential to redefine the cost of cultural production. It showcases the utility of GPTs in digital community management, expanding their implications in cultural economics and transmedia storytelling.


title: GP+: A Python Library for Kernel-based learning via Gaussian Processes

abstract: In this paper we introduce GP+, an open-source library for kernel-based learning via Gaussian processes (GPs) which are powerful statistical models that are completely characterized by their parametric covariance and mean functions. GP+ is built on PyTorch and provides a user-friendly and object-oriented tool for probabilistic learning and inference. As we demonstrate with a host of examples, GP+ has a few unique advantages over other GP modeling libraries. We achieve these advantages primarily by integrating nonlinear manifold learning techniques with GPs' covariance and mean functions. As part of introducing GP+, in this paper we also make methodological contributions that (1) enable probabilistic data fusion and inverse parameter estimation, and (2) equip GPs with parsimonious parametric mean functions which span mixed feature spaces that have both categorical and quantitative variables. We demonstrate the impact of these contributions in the context of Bayesian optimization, multi-fidelity modeling, sensitivity analysis, and calibration of computer models.


title: Non-mechanical steering of the optical beam in spectral-domain Optical Coherence Tomography

abstract: We demonstrate in a {\it proof-of-concept} experiment spectral-domain optical coherence tomography where steering of the optical beam that probes the sample in a transverse scan does not make use of any mechanical element. Steering is done with the help of a spatial light modulator, that introduces a spatially-dependent phase between the two orthogonal polarization components of an optical beam, and some optical elements that control the polarization of light. We demonstrate that making use of the non-mechanical steering system considered here, we can reproduce the main traits of imaging with standard OCT that makes use of mechanical-assisted optical beam steering.


title: Real-time Network Intrusion Detection via Decision Transformers

abstract: Many cybersecurity problems that require real-time decision-making based on temporal observations can be abstracted as a sequence modeling problem, e.g., network intrusion detection from a sequence of arriving packets. Existing approaches like reinforcement learning may not be suitable for such cybersecurity decision problems, since the Markovian property may not necessarily hold and the underlying network states are often not observable. In this paper, we cast the problem of real-time network intrusion detection as casual sequence modeling and draw upon the power of the transformer architecture for real-time decision-making. By conditioning a causal decision transformer on past trajectories, consisting of the rewards, network packets, and detection decisions, our proposed framework will generate future detection decisions to achieve the desired return. It enables decision transformers to be applied to real-time network intrusion detection, as well as a novel tradeoff between the accuracy and timeliness of detection. The proposed solution is evaluated on public network intrusion detection datasets and outperforms several baseline algorithms using reinforcement learning and sequence modeling, in terms of detection accuracy and timeliness.


title: A Class of Computational Methods to Reduce Selection Bias when Designing Phase 3 Clinical Trials

abstract: When designing confirmatory Phase 3 studies, one usually evaluates one or more efficacious and safe treatment option(s) based on data from previous studies. However, several retrospective research articles reported the phenomenon of ``diminished treatment effect in Phase 3'' based on many case studies. Even under basic assumptions, it was shown that the commonly used estimator could substantially overestimate the efficacy of selected group(s). As alternatives, we propose a class of computational methods to reduce estimation bias and mean squared error (MSE) with a broader scope of multiple treatment groups and flexibility to accommodate summary results by group as input. Based on simulation studies and a real data example, we provide practical implementation guidance for this class of methods under different scenarios. For more complicated problems, our framework can serve as a starting point with additional layers built in. Proposed methods can also be widely applied to other selection problems.


title: Machine Learning and Citizen Science Approaches for Monitoring the Changing Environment

abstract: This dissertation will combine new tools and methodologies to answer pressing questions regarding inundation area and hurricane events in complex, heterogeneous changing environments. In addition to remote sensing approaches, citizen science and machine learning are both emerging fields that harness advancing technology to answer environmental management and disaster response questions.


title: Science research from the Instituto Argentino de Radioastronomia

abstract: In this article, I will present some figures and milestones of the written production of the Instituto Argentino de Radioastronomia (IAR), as well as a personal review of the scientific achievements carried out in recent years by the researchers working at the IAR. I will also briefly describe the scientific objectives of the IAR's flagship project, the Multipurpose Interferometric Array (MIA), in the context of the instrumental projects that have lately been or are being installed on Argentine soil.


title: Analysis of Drop Collisions

abstract: When do droplets merge and when do they bounce? Over the last 10 years, advances in experimental techniques, such as high-speed cameras, have enabled us to make important discoveries on how the dynamics of thin gas films can influence the behaviour of liquid droplets. However, the exact physical mechanisms governing these phenomena are a subject of debate, leading to a variety of theoretical, experimental and computational approaches. In this document, I will study two problems: the head-on collision of two droplets, and the impact of a drop on a solid surface. I will outline the conventional fluid mechanical description, how it fails, and various candidates for the ``missing physics''.


title: Axion String Source Modelling

abstract: In this paper, we perform an investigation into the effect of the string radius of curvature $R_\mathrm{\,Gaussian}$ on the magnitude and relative magnitude of the massive and massless radiation from axion (global) string configurations, motivated by qualitative observations from string network simulations. We construct initial conditions from travelling wave solutions on a global string for two colliding Gaussians, performing parameter scans over amplitude $A$ and standard deviation $\sigma_\mathrm{d}$. We show that the energy emitted via massless radiation obeys a power law $E_\mathrm{massless} \propto A^{\gamma}$, where the coefficient $\gamma$ depends on the curvature regime. Massive radiation is exponentially suppressed approximately as $E_{\mathrm{massive}} \propto e^{-\zeta R_\mathrm{\,Gaussian}}$ in the quasi-linear regime $\sigma_\mathrm{d} \gg \delta$ and exhibits power-law decay $E_{\mathrm{massive}} \propto (R_\mathrm{\,Gaussian})^{-\gamma}$ in the nonlinear regime where $\sigma_\mathrm{d} \lesssim 2\delta$, with different $\gamma$ in different regimes of $R_\mathrm{\,Gaussian}$. In certain regions of the nonlinear regime, massive particle radiation comprises up to 50% of the total energy emitted. Drawing on a known parallel between axion radiation from global strings and gravitational radiation from Abelian-Higgs strings, this suggests that massive particle radiation channel may become of equal significance to the massless (gravitational) channel for nonlinear burst signals where $R < \sigma_\mathrm{d}$, unless we are in the regime where additional loops are generated. We also estimate the spectral index $q$ of the axion radiation for different amplitudes, showing that a higher proportion of radiation is emitted in high frequency modes as the curvature increases, bounded by $q \gtrsim 1$ for the configurations studied.


title: Powers averaging for actions on $C(X)$-algebras

abstract: Given a unital $C(X)$-algebra $A$ discrete group $\Gamma$ and an action $\alpha: \Gamma\to \text{aut}(A)$ which leaves $C(X)$ invariant and such that $C(X)\rtimes_{\alpha,r} \Gamma$ is simple, and a $2$-cocycle $\omega$, we obtain a bijective correspondence between maximal $\Gamma$-invariant ideals of $A$ and maximal ideals in $A\rtimes_{\alpha,\omega,r} \Gamma$. In particular, $A\rtimes_{\alpha,\omega,r} \Gamma$ is simple if and only if $A$ has no $\Gamma$-invariant ideals.


title: Nash equilibria for dividend distribution with competition

abstract: We construct a Nash equilibrium in feedback form for a class of two-person stochastic games with absorption arising from corporate finance. More precisely, the paper focusses on a strategic dynamic game in which two financially-constrained firms operate in the same market. The firms distribute dividends and are faced with default risk. The strategic interaction arises from the fact that if one firm defaults, the other one becomes a monopolist and increases its profitability. To determine a Nash equilibrium in feedback form, we develop two different concepts depending on the initial endowment of each firm. If one firm is richer than the other one, then we use a notion of control vs.\ strategy equilibrium. If the two firms have the same initial endowment (hence they are symmetric in our setup) then we need mixed strategies in order to construct a symmetric equilibrium.


title: Distribution of the elemental regression weights with t-distributed co-variate measurement errors

abstract: In this article, a heuristic approach is used to determined the best approximate distribution of $\dfrac{Y_1}{Y_1 + Y_2}$, given that $Y_1,Y_2$ are independent, and each of $Y_1$ and $Y$ is distributed as the $\mathcal{F}$-distribution with common denominator degrees of freedom. The proposed approximate distribution is subject to graphical comparisons and distributional tests. The proposed distribution is used to derive the distribution of the elemental regression weight $\omega_E$, where $E$ is the elemental regression set.


title: Brain-optimized inference improves reconstructions of fMRI brain activity

abstract: The release of large datasets and developments in AI have led to dramatic improvements in decoding methods that reconstruct seen images from human brain activity. We evaluate the prospect of further improving recent decoding methods by optimizing for consistency between reconstructions and brain activity during inference. We sample seed reconstructions from a base decoding method, then iteratively refine these reconstructions using a brain-optimized encoding model that maps images to brain activity. At each iteration, we sample a small library of images from an image distribution (a diffusion model) conditioned on a seed reconstruction from the previous iteration. We select those that best approximate the measured brain activity when passed through our encoding model, and use these images for structural guidance during the generation of the small library in the next iteration. We reduce the stochasticity of the image distribution at each iteration, and stop when a criterion on the "width" of the image distribution is met. We show that when this process is applied to recent decoding methods, it outperforms the base decoding method as measured by human raters, a variety of image feature metrics, and alignment to brain activity. These results demonstrate that reconstruction quality can be significantly improved by explicitly aligning decoding distributions to brain activity distributions, even when the seed reconstruction is output from a state-of-the-art decoding algorithm. Interestingly, the rate of refinement varies systematically across visual cortex, with earlier visual areas generally converging more slowly and preferring narrower image distributions, relative to higher-level brain areas. Brain-optimized inference thus offers a succinct and novel method for improving reconstructions and exploring the diversity of representations across visual brain areas.


title: Near-Optimal Differentially Private k-Core Decomposition

abstract: Recent work by Dhulipala et al. \cite{DLRSSY22} initiated the study of the $k$-core decomposition problem under differential privacy via a connection between low round/depth distributed/parallel graph algorithms and private algorithms with small error bounds. They showed that one can output differentially private approximate $k$-core numbers, while only incurring a multiplicative error of $(2 +\eta)$ (for any constant $\eta >0$) and additive error of $\poly(\log(n))/\eps$. In this paper, we revisit this problem. Our main result is an $\eps$-edge differentially private algorithm for $k$-core decomposition which outputs the core numbers with no multiplicative error and $O(\text{log}(n)/\eps)$ additive error. This improves upon previous work by a factor of 2 in the multiplicative error, while giving near-optimal additive error. Our result relies on a novel generalized form of the sparse vector technique, which is especially well-suited for threshold-based graph algorithms; thus, we further strengthen the connection between distributed/parallel graph algorithms and differentially private algorithms.


title: Learning for System Identification of NDAE-modeled Power Systems

abstract: System identification through learning approaches is emerging as a promising strategy for understanding and simulating dynamical systems, which nevertheless faces considerable difficulty when confronted with power systems modeled by differential-algebraic equations (DAEs). This paper introduces a neural network (NN) framework for effectively learning and simulating solution trajectories of DAEs. The proposed framework leverages the synergy between Implicit Runge-Kutta (IRK) time-stepping schemes tailored for DAEs and NNs (including a differential NN (DNN)). The framework enforces an NN to cooperate with the algebraic equation of DAEs as hard constraints and is suitable for the identification of the ordinary differential equation (ODE)-modeled dynamic equation of DAEs using an existing penalty-based algorithm. Finally, the paper demonstrates the efficacy and precision of the proposed NN through the identification and simulation of solution trajectories for the considered DAE-modeled power system.


title: Bose-Einstein condensation by polarization gradient laser cooling

abstract: Attempts to create quantum degenerate gases without evaporative cooling have been pursued since the early days of laser cooling, with the consensus that polarization gradient cooling (PGC, also known as "optical molasses") alone cannot reach condensation. In the present work, we report that simple PGC can generate a small Bose-Einstein condensate (BEC) inside a corrugated micrometer-sized optical dipole trap. The experimental parameters enabling BEC creation were found by machine learning, which increased the atom number by a factor of 5 and decreased the temperature by a factor of 2.5, corresponding to almost two orders of magnitude gain in phase space density. When the trapping light is slightly misaligned through a microscopic objective lens, a BEC of $\sim 250$ $^{87}$Rb atoms is formed inside a local dimple within 40 ms of PGC.


title: Majority is Not Required: A Rational Analysis of the Private Double-Spend Attack from a Sub-Majority Adversary

abstract: We study the incentives behind double-spend attacks on Nakamoto-style Proof-of-Work cryptocurrencies. In these systems, miners are allowed to choose which transactions to reference with their block, and a common strategy for selecting transactions is to simply choose those with the highest fees. This can be problematic if these transactions originate from an adversary with substantial (but less than 50\%) computational power, as high-value transactions can present an incentive for a rational adversary to attempt a double-spend attack if they expect to profit. The most common mechanism for deterring double-spend attacks is for the recipients of large transactions to wait for additional block confirmations (i.e., to increase the attack cost). We argue that this defense mechanism is not satisfactory, as the security of the system is contingent on the actions of its users. Instead, we propose that defending against double-spend attacks should be the responsibility of the miners; specifically, miners should limit the amount of transaction value they include in a block (i.e., reduce the attack reward). To this end, we model cryptocurrency mining as a mean-field game in which we augment the standard mining reward function to simulate the presence of a rational, double-spending adversary. We design and implement an algorithm which characterizes the behavior of miners at equilibrium, and we show that miners who use the adversary-aware reward function accumulate more wealth than those who do not. We show that the optimal strategy for honest miners is to limit the amount of value transferred by each block such that the adversary's expected profit is 0. Additionally, we examine Bitcoin's resilience to double-spend attacks. Assuming a 6 block confirmation time, we find that an attacker with at least 25% of the network mining power can expect to profit from a double-spend attack.


title: The classifying element for quotients of Fermat curves

abstract: Suppose $C$ is a cyclic Galois cover of the projective line branched at the three points $0$, $1$, and $\infty$. Under a mild condition on the ramification, we determine the structure of the graded Lie algebra of the lower central series of the fundamental group of $C$ in terms of a basis which is well-suited to studying the action of the absolute Galois group of $\mathbb{Q}$.


title: Leveraging Large Language Models to Build and Execute Computational Workflows

abstract: The recent development of large language models (LLMs) with multi-billion parameters, coupled with the creation of user-friendly application programming interfaces (APIs), has paved the way for automatically generating and executing code in response to straightforward human queries. This paper explores how these emerging capabilities can be harnessed to facilitate complex scientific workflows, eliminating the need for traditional coding methods. We present initial findings from our attempt to integrate Phyloflow with OpenAI's function-calling API, and outline a strategy for developing a comprehensive workflow management system based on these concepts.


title: The DEAP-3600 liquid argon optical model and NEST updates

abstract: As liquid argon (LAr) detectors are made at progressively larger sizes, accurate models of LAr optical properties become increasingly important for simulating light transport, understanding signals, and developing analyses. The refractive index, group velocity, and Rayleigh scattering length are particularly important for vacuum ultraviolet (VUV) and visible photons in detectors with diameters much greater than one meter. While optical measurements in the VUV are sparse, recent measurements of the group velocity of 128 nm photons in LAr provide valuable constraints on these parameters. These calculations are further complicated by the dependence of optical parameters on thermodynamic properties that might fluctuate or vary throughout the argon volume. This manuscript presents the model used by DEAP-3600, a dark matter direct detection experiment at SNOLAB using a 3.3 tonne LAr scintillation counter. Existing data and thermodynamic models are synthesized to estimate the wavelength-dependent refractive index, group velocity, and Rayleigh scattering length within the detector, and parameters' uncertainties are estimated. This model, along with in situ measurements of LAr scintillation properties, is benchmarked against data collected in DEAP-3600, providing a method for modeling optical properties in large LAr detectors and for propagating their uncertainties through downstream simulations. Updates are also presented of the Noble Element Simulation Technique (NEST) software, widely used to model scintillation and ionization signals in argon- and xenon-based detectors.


title: Lagrangian Liouville models of multiphase flows with randomly forced inertial particles

abstract: Eulerian-Lagrangian models of particle-laden (multiphase) flows describe fluid flow and particle dynamics in the Eulerian and Lagrangian frameworks respectively. Regardless of whether the flow is turbulent or laminar, the particle dynamics is stochastic because the suspended particles are subjected to random forces. We use a polynomial chaos expansion (PCE), rather than a postulated constitutive law, to capture structural and parametric uncertainties in the particles' forcing. The stochastic particle dynamics is described by a joint probability density function (PDF) of a particle's position and velocity and random coefficients in the PCE. We deploy the method of distributions (MoD) to derive a deterministic (Liouville-type) partial-differential equation for this PDF. We reformulate this PDF equation in a Lagrangian form, obtaining PDF flow maps and tracing events and their probability in the phase space. That is accomplished via a new high-order spectral scheme, which traces, marginalizes and computes moments of the high-dimensional joint PDF and comports with high-order carrier-phase solvers. Our approach has lower computational cost than either high-order Eulerian solvers or Monte Carlo methods, is not subjected to a CFL condition, does not suffer from Gibbs oscillations and does not require (order-reducing) filtering and regularization techniques. These features are demonstrated on several test cases.


title: Internal structure and analytical representation of preference relations defined on infinite-dimensional real vector spaces

abstract: The paper deals with partial and weak preference relations defined on infinite-dimensional vector spaces and compatible with algebraic operations. By a partial preference we mean an asymmetric and transitive binary relation, while a weak preference is such a partial preference for which the indifference relation corresponding it is transitive (an indifference relation is the complement to the union of a partial preference and the reverse to it). Our first aim is to study the internal structure of compatible partial and weak preferences. Using these results we then prove that a compatible weak preference admits an analytical representation by means of a step-linear function, while a compatible partial preference can be analytically represent by the family of step-linear functions.


title: Confined magnon dispersion in ferromagnetic and antiferromagnetic thin films in a second quantization approach: the case of Fe and NiO

abstract: We present a methodology based on the calculation of the inelastic scattering from magnons via the spin scattering function in confined geometries such as thin films using a second quantization formalism, for both ferromagnetic and antiferromagnetic materials. The case studies are chosen with an aim to demonstrate the effects of film thickness and crystal orientation on magnon modes, using bcc Fe(100) and NiO with (100) and (111) crystallographic orientations as prototypical systems. Due to the quantization of the quasi-momentum we observe a granularity in the inelastic spectra in the reciprocal space path reflecting the orientation of the thin film. This approach also allows to capture softer modes that appear due to the partial interaction of magnetic moments close to the surface in a thin film geometry, in addition to bulk modes. The softer modes are also affected by crystallographic orientations as illustrated by the different surface-related peaks of NiO magnon density of states at approximately ~ 65 meV for (100) and ~ 42 meV for (111). Additionally, we explore the role of anisotropy on magnon modes, revealing that introducing anisotropy to both Fe and NiO films increases the overall hardness of the magnon modes. The introduction of a surface anisotropy produces a shift of the surface-related magnon DOS peak to higher energies with increased surface anisotropy, and in some cases leading to surface confined mode.


title: Permutations avoiding bipartite partially ordered patterns have a regular insertion encoding

abstract: We prove that any class of permutations defined by avoiding a partially ordered pattern (POP) with height at most two has a regular insertion encoding and thus has a rational generating function. Then, we use Combinatorial Exploration to find combinatorial specifications and generating functions for hundreds of other permutation classes defined by avoiding a size 5 POP, allowing us to resolve several conjectures of Gao and Kitaev and of Chen and Lin.


title: Understanding disorder in Silicon quantum computing platforms: Scattering mechanisms in Si/SiGe quantum wells

abstract: Motivated by recent experiments on Si/SiGe quantum wells with a co-design of high electron mobility and large valley splitting [B. Paquelet Wuetz, et al., Nature Communications 14, 1385 (2023); D. D. Esposti, et al., arXiv:2309.02832], suitable for a Si-based spin qubit quantum computing platform, we examine the role of disorder by theoretically calculating mobility and quantum mobility from various scattering mechanisms and their dependence on the electron density. At low electron densities $n_e < 4 \times 10^{11}$ cm$^{-2}$, we find that mobility is limited by remote Coulomb impurities in the capping layer, whereas interface roughness becomes the significant limiting factor at higher densities. We also find that alloy disorder scattering is not a limiting mechanism in the reported high-mobility structures. We estimate the critical density of the disorder-driven low-density metal-insulator transition using the Anderson-Ioffe-Regel localization criterion and qualitatively explain the breakdown of the Boltzmann-Born theory at low densities. We also estimate the critical density by considering inhomogeneous density fluctuations induced by long-range Coulomb disorder in the system, and find a larger critical density compared to the one obtained from the Anderson-Ioffe-Regel criterion. For quantum mobility, our calculation suggests that remote and distant background impurities are likely the limiting scattering sources across all density. Future measurements of quantum mobility should provide more information on the distribution of background impurities inside the SiGe barriers. Moreover, we extend our theoretical analysis to the effect of quantum degeneracy on transport properties and predict the mobility and the critical density for the metal-insulator transition in spin-polarized high-mobility structures under an external parallel magnetic field.


title: CaVE: A Cone-Aligned Approach for Fast Predict-then-optimize with Binary Linear Programs

abstract: The end-to-end predict-then-optimize framework, also known as decision-focused learning, has gained popularity for its ability to integrate optimization into the training procedure of machine learning models that predict the unknown cost (objective function) coefficients of optimization problems from contextual instance information. Naturally, most of the problems of interest in this space can be cast as integer linear programs. In this work, we focus on binary linear programs (BLPs) and propose a new end-to-end training method to predict-then-optimize. Our method, Cone-aligned Vector Estimation (CaVE), aligns the predicted cost vectors with the normal cone corresponding to the true optimal solution of a training instance. When the predicted cost vector lies inside the cone, the optimal solution to the linear relaxation of the binary problem is optimal. This alignment not only produces decision-aware learning models but also dramatically reduces training time as it circumvents the need to solve BLPs to compute a loss function with its gradients. Experiments across multiple datasets show that our method exhibits a favorable trade-off between training time and solution quality, particularly with large-scale optimization problems such as vehicle routing, a hard BLP that has yet to benefit from predict-then-optimize methods in the literature due to its difficulty.