title: Time evolution of the galactic $B- \rho$ relation: the impact of the magnetic field morphology

abstract: One of the most frequently used indicators to characterize the magnetic field's influence on star formation is the relation between magnetic field strength and gas density ($B-\rho$ relation), usually expressed as $B \propto \rho^{\kappa}$. The value of $\kappa$ is an indication of the dynamical importance of the magnetic field during gas compression. Investigating the global magnetic field's impact on this relation and its evolution, we conduct MHD simulations of Milky-Way-like galaxies including gravity, star formation, and supernova feedback along with non-equilibrium chemistry up to $H_2$ formation fueling star formation. Two initial magnetic field morphologies are studied: one completely ordered (toroidal) and the other completely random. In these models, we study the dynamical importance of the magnetic field through the plasma $\beta$ and the $B-\rho$ relation. For both magnetic morphologies, low-density regions are thermally supported, while high-density regions are magnetically dominated. Equipartition is reached earlier and at lower densities in the toroidal model. However, the $B-\rho$ relation is not unique even within the same galaxy, as it consistently includes two different branches for a given density, with $\kappa$ ranging from about 0.2 to 0.8. The mean value of $\kappa$ for each model also displays significant variations over time, which supersede the differences between the two models. While our findings suggest that the magnetic field morphology does influence the galactic $B-\rho$ relation, its impact is transient, since time-averaged differences between the models fall within the large temporal scatter. The context and time-dependent nature of the $B-\rho$ relation underscore the need for comprehensive research and observations to understand the intricate role of magnetic fields in star formation processes across diverse galactic environments.


title: On the universal exceptional structure of world-volume theories in string and M-theory

abstract: A universal structure of world-volume theories of half-BPS branes in string and M-theory in terms of exceptional generalised geometry is observed. Previous constructions are extended in two ways: from internal $d$-dimensional space to full 11- or 10-dimensional spacetimes, by coupling to the tensor hierarchy gauge fields, and from $E_{d(d)}$ with $d\leq 6$ to $E_{7(7)}$ and $E_{8(8)}$. This is done by a clarification of the role of the tensor hierarchy and its end in the context of brane world-volume theories. The exceptional structure of the gauged $\sigma$-model of the Kaluza-Klein monopole is provided as a new example.


title: CLusters in the Uv as EngineS (CLUES). II. Sub-kpc scale outflows driven by stellar feedback

abstract: We analyze the far-ultraviolet spectroscopy of 20 young and massive star clusters (YSCs) in 11 nearby star-forming galaxies. We probe the interstellar gas intervening along the line of sight, detecting several metal absorption lines of a wide range of ionization potentials, from 6.0 eV to 77.5 eV. Multiple-component Voigt fits to the absorption lines are used to study the kinematics of the gas. We find that nearly all targets in the sample feature gas outflowing from 30 up to 190 km per second, often both in the neutral and ionized phase. The outflow velocities correlate with the underlying stellar population properties directly linked to the feedback: the mass of the YSCs, the photon production rate and the instantaneous mechanical luminosity produced by stellar winds and SNe. We detect a neutral inflow in 4 targets, which we interpret as likely not associated with the star cluster but tracing larger scale gas kinematics. A comparison between the outflows energy and that produced by the associated young stellar populations suggests an average coupling efficiency of 10 per cent with a broad scatter. Our results extend the relation found in previous works between galactic outflows and the host galaxy star-formation rate to smaller scales, pointing towards the key role that clustered star formation and feedback play in regulating galaxy growth.


title: Investigating finite-size effects in random matrices by counting resonances

abstract: Resonance counting is an intuitive and widely used tool in Random Matrix Theory and Anderson Localization. Its undoubted advantage is its simplicity: in principle, it is easily applicable to any random matrix ensemble. On the downside, the notion of resonance is ill-defined, and the `number of resonances' does not have a direct mapping to any commonly used physical observable like the participation entropy, the fractal dimensions, or the gap ratios (r-parameter), restricting the method's predictive power to the thermodynamic limit only where it can be used for locating the Anderson localization transition. In this work, we reevaluate the notion of resonances and relate it to measurable quantities, building a foundation for the future application of the method to finite-size systems. To access the HTML version of the paper & discuss it with the authors, visit https://enabla.com/pub/558.


title: Operational Methods Applied to the Spherical Mean and X-Ray Transform

abstract: We use the framework of operational calculus to derive the operators associated with the spherical mean and related means of a function in $n$-dimensional space. These operators are then used to simplify the treatment and discussion of the means, i.e., to determine their power series expansion, to find the series solution to the corresponding inversion problem, to derive the associated differential equations that the means satisfy, or give meaning to fractional spherical means. Finally, we use the framework to derive the inversion formula of the X-ray transform in an operational manner.


title: Variational Autoencoders for Regression: Recovering Fully Leptonic $b\bar{b}W^+W^-$ in Di-Higgs Searches

abstract: The search for double Higgs production in $b\bar{b}W^+W^-$, where both $W$ bosons decay to leptons, has been rehabilitated as a good option to look for that key process to the Standard Model scalar sector study in the LHC. The missing neutrinos, however, hinder the reconstruction of useful information like the Higgs pair mass, which is very sensitive to the trilinear Higgs self-coupling. We present a solution to that problem using a Variational Autoencoder for Regression (VAER) to reconstruct the Higgs and top pairs decays $hh,t\bar{t}\to b\bar{b}W^+W^-\to b\bar{b}\ell^+\ell^{\prime -}\nu_\ell\bar{\nu}_{\ell^\prime}$. The algorithm predicts the invariant mass of non-resonant $hh$ irrespective of the trilinear coupling, even for events whose Higgs self-couplings were never presented to it. VAER is also able to identify a new Higgs resonance in an unsupervised way, showing generalization power for events not presented in its training phase. Finally, we demonstrate that VAER prediction is as useful to statistical inference as ground truth simulated distributions by computing a $\chi^2$ between trilinear coupling hypotheses based on binned invariant mass distributions of $b\bar{b}\ell^+\ell^{\prime -}\nu_\ell\bar{\nu}_{\ell^\prime}$.


title: Additively manufactured ceramics for compact quantum technologies

abstract: Quantum technologies are advancing from fundamental research in specialized laboratories to practical applications in the field, driving the demand for robust, scalable, and reproducible system integration techniques. Ceramic components can be pivotal thanks to high stiffness, low thermal expansion, and excellent dimensional stability under thermal stress. We explore lithography-based additive manufacturing of technical ceramics especially for miniaturized physics packages and electro-optical systems. This approach enables functional systems with precisely manufactured, intricate structures and high mechanical stability while minimizing size and weight. It facilitates rapid prototyping, simplifies fabrication and leads to highly integrated, reliable devices. As an electrical insulator with low outgassing and high temperature stability, printed technical ceramics such as Al2O3 and AlN bridge a technology gap in quantum technology and offer advantages over other printable materials. We demonstrate this potential with CerAMRef, a micro-integrated rubidium D2 line optical frequency reference on a printed Al2O3 micro-optical bench and housing. The frequency instability of the reference is comparable to laboratory setups while the volume of the integrated spectroscopy setup is only 6 ml. We identify potential for future applications in compact atomic magnetometers, miniaturized optical atom traps, and vacuum system integration.


title: Quantum optics with giant atoms in a structured photonic bath

abstract: We present a general framework to tackle quantum optics problems with giant atoms, i.e. quantum emitters each coupled {\it non-locally} to a structured photonic bath (typically a lattice) of any dimension. The theory encompasses the calculation and general properties of Green's functions, atom-photon bound states (BSs), collective master equations and decoherence-free Hamiltonians (DFHs), and is underpinned by a formalism where a giant atom is formally viewed as a normal atom lying at a fictitious location. As a major application, we provide for the first time a general criterion to predict/engineer DFHs of giant atoms, which can be applied both in and out of the photonic continuum and regardless of the structure or dimensionality of the photonic bath. This is used to show novel DFHs in 2D baths such as a square lattice and photonic graphene.


title: Extracting the current-phase-relation of a monolithic three-dimensional nano-constriction using a DC-current-tunable superconducting microwave cavity

abstract: Superconducting circuits with nonlinear elements such as Josephson tunnel junctions or kinetic inductance nanowires are the workhorse for microwave quantum and superconducting sensing technologies. For devices, which can be operated at high temperatures and large magnetic fields, nano-constrictions as nonlinear elements are recently under intense investigation. Constrictions, however, are far less understood than conventional Josephson tunnel junctions, and their current-phase-relationships (CPRs) -- although highly important for device design -- are hard to predict. Here, we present a niobium microwave cavity with a monolithically integrated, neon-ion-beam patterned three-dimensional (3D) nano-constriction. By design, we obtain a DC-current-tunable microwave circuit and characterize how the bias-current-dependent constriction properties impact the cavity resonance. Based on the results of these experiments, we reconstruct the CPR of the nanoconstriction. Finally, we discuss the Kerr nonlinearity of the device, a parameter important for many high-dynamic-range applications and an experimental probe for the second and third derivatives of the CPR. Our platform provides a useful method to comprehensively characterize nonlinear elements integrated in microwave circuits and could be of interest for current sensors, hybrid quantum systems and parametric amplifiers. Our findings furthermore contribute to a better understanding of nano-fabricated 3D constrictions.


title: Exploiting symmetries in nuclear Hamiltonians for ground state preparation

abstract: The Lipkin and Agassi models are simplified nuclear models that provide natural test beds for quantum simulation methods. Prior work has investigated the suitability of the Variational Quantum Eigensolver (VQE) to find the ground state of these models. There is a growing awareness that if VQE is to prove viable, we will need problem inspired ans\"{a}tze that take into account the symmetry properties of the problem and use clever initialization strategies. Here, by focusing on the Lipkin and Agassi models, we investigate how to do this in the context of nuclear physics ground state problems. We further use our observations to discus the potential of new classical, but quantum-inspired, approaches to learning ground states in nuclear problems.


title: The Kormendy relation of cluster galaxies in PPS regions

abstract: We study a sample of 936 early-type galaxies located in 48 low-z regular galaxy clusters with $M_{200}\geq 10^{14}~ M_\odot$ at $z< 0.1$. We examine variations in the Kormendy relation (KR) according to their location in the projected phase space (PPS) of the clusters. We have used a combination of Bayesian statistical methods to identify possible differences between the fitted relations. Our results indicate that the overall KR is better fitted when we take into account the information about PPS regions. We also find that objects with time since infall $\geq 6.5$ Gyr have a significant statistical difference of the KR coefficients relative to objects that are more recent in the cluster environment. We show that giant central ellipticals are responsible for tilting the KR relation towards smaller slopes. These galaxies present a late growth probably due to cumulative preprocessing during infall, plus cannibalism and accretion of smaller stripped objects near the center of the clusters.


title: Geometric logarithmic-Hardy and Hardy-Poincar\'e inequalities on stratified groups

abstract: We develop a unified strategy to obtain the geometric logarithmic Hardy inequality on any open set M of a stratified group, provided the validity of the Hardy inequality in this setting, where the so-called "weight" is regarded to be any measurable non-negative function on M . Provided the legitimacy of the latter for some open set and for some weight, we also show an inequality that is an extension of the "generalised Poincar\'e inequality" introduced by Beckner with the addition of a weight, and this is referred to as the "geometric Hardy-Poincar\'e inequality". The aforesaid inequalities become explicit in the case where M is the half-space of the group and the weight is the distance function from the boundary, and in the case where M is just the whole group (or any open set in the group), in which case the weight is the "horizontal norm" on the first stratum of the group. For the second case, the semi-Gaussian analogue of the derived inequalities is proved, when the Gaussian measure is regarded with respect to the first stratum of the group. Applying our results to the case where the group is just the (abelian) Euclidean space we generalise the classical probabilistic Poincar\'e inequality by adding weights.


title: SusFL: Energy-Aware Federated Learning-based Monitoring for Sustainable Smart Farms

abstract: We propose a novel energy-aware federated learning (FL)-based system, namely SusFL, for sustainable smart farming to address the challenge of inconsistent health monitoring due to fluctuating energy levels of solar sensors. This system equips animals, such as cattle, with solar sensors with computational capabilities, including Raspberry Pis, to train a local deep-learning model on health data. These sensors periodically update Long Range (LoRa) gateways, forming a wireless sensor network (WSN) to detect diseases like mastitis. Our proposed SusFL system incorporates mechanism design, a game theory concept, for intelligent client selection to optimize monitoring quality while minimizing energy use. This strategy ensures the system's sustainability and resilience against adversarial attacks, including data poisoning and privacy threats, that could disrupt FL operations. Through extensive comparative analysis using real-time datasets, we demonstrate that our FL-based monitoring system significantly outperforms existing methods in prediction accuracy, operational efficiency, system reliability (i.e., mean time between failures or MTBF), and social welfare maximization by the mechanism designer. Our findings validate the superiority of our system for effective and sustainable animal health monitoring in smart farms. The experimental results show that SusFL significantly improves system performance, including a $10\%$ reduction in energy consumption, a $15\%$ increase in social welfare, and a $34\%$ rise in Mean Time Between Failures (MTBF), alongside a marginal increase in the global model's prediction accuracy.


title: Comparison of variational discretizations for a convection-diffusion problem

abstract: For a model convection-diffusion problem, we obtain new error estimates for a general upwinding finite element discretization based on bubble modification of the test space. The key analysis tool is based on finding representations of the optimal norms on the trial spaces at the continuous and discrete levels. We analyze and compare the standard linear discretization, the saddle point least square and upwinding Petrov-Galerkin methods. We conclude that the bubble upwinding Petrov-Galerkin method is the most performant discretization for the one dimensional model. Our results for the model convection-diffusion problem can be extended for creating new and efficient discretizations for the multidimensional cases.


title: Information Capacity Regret Bounds for Bandits with Mediator Feedback

abstract: This work addresses the mediator feedback problem, a bandit game where the decision set consists of a number of policies, each associated with a probability distribution over a common space of outcomes. Upon choosing a policy, the learner observes an outcome sampled from its distribution and incurs the loss assigned to this outcome in the present round. We introduce the policy set capacity as an information-theoretic measure for the complexity of the policy set. Adopting the classical EXP4 algorithm, we provide new regret bounds depending on the policy set capacity in both the adversarial and the stochastic settings. For a selection of policy set families, we prove nearly-matching lower bounds, scaling similarly with the capacity. We also consider the case when the policies' distributions can vary between rounds, thus addressing the related bandits with expert advice problem, which we improve upon its prior results. Additionally, we prove a lower bound showing that exploiting the similarity between the policies is not possible in general under linear bandit feedback. Finally, for a full-information variant, we provide a regret bound scaling with the information radius of the policy set.


title: Backdoor Attack against One-Class Sequential Anomaly Detection Models

abstract: Deep anomaly detection on sequential data has garnered significant attention due to the wide application scenarios. However, deep learning-based models face a critical security threat - their vulnerability to backdoor attacks. In this paper, we explore compromising deep sequential anomaly detection models by proposing a novel backdoor attack strategy. The attack approach comprises two primary steps, trigger generation and backdoor injection. Trigger generation is to derive imperceptible triggers by crafting perturbed samples from the benign normal data, of which the perturbed samples are still normal. The backdoor injection is to properly inject the backdoor triggers to comprise the model only for the samples with triggers. The experimental results demonstrate the effectiveness of our proposed attack strategy by injecting backdoors on two well-established one-class anomaly detection models.


title: The intermediate neutron capture process. V. The i-process in AGB stars with overshoot

abstract: The intermediate neutron capture process (i-process) can develop during proton ingestion events (PIE), potentially during the early stages of low-mass low-metallicity asymptotic giant branch (AGB) stars. We examine the impact of overshoot mixing on the triggering and development of i-process nucleosynthesis in AGB stars of various initial masses and metallicities. We computed AGB stellar models, with initial masses of 1, 2, 3, and 4 M$_{\odot}$ and metallicities in the $-2.5 \le $ [Fe/H] $\le 0$ range, using the stellar evolution code STAREVOL with a network of 1160 nuclei coupled to the transport equations. We considered different overshooting profiles below and above the thermal pulses, and below the convective envelope. The occurrence of PIEs is found to be primarily governed by the amount of overshooting at the top of pulse ($f_{\rm top}$) and to increase with rising $f_{\rm top}$. For $f_{\rm top} =$ 0, 0.02, 0.04, and 0.1, we find that 0 %, 6 %, 24 %, and 86 % of our 21 AGB models with $-2<$ [Fe/H] $<0$ experience a PIE, respectively. We also find that PIEs leave a $^{13}$C-pocket at the bottom of the pulse that can give rise to an additional radiative s-process nucleosynthesis, and ultimately produce a noticeable mixed i+s chemical signature at the surface. Finally, the chemical abundance patterns of 22 observed r/s-stars candidates with $-2<$ [Fe/H] $<-1$ are found to be in reasonable agreement with our AGB model predictions. The binary status of the dwarfs/giants being unclear, we suggest that these stars have acquired their chemical pattern either from the mass transfer of a now-extinct AGB companion or from an early generation AGB star that polluted the natal cloud. Stricter constraints from multi-dimensional hydrodynamical models on overshoot coefficients could deliver new insights into the contribution of AGB stars to heavy elements in the Universe.


title: The affect of Some Meteorological Parameters on Particulate Matters Concentration Over Iraq using Remote Sensing dataset

abstract: Numerous countries have built urban stations for monitoring the amount of PM2.5 in the atmosphere. In Iraq, there aren't enough stations to monitor PM2.5 pollution levels across all governorates. As a result, satellite remote sensing data is used in the majority of studies aimed at monitoring PM2.5 and the impact of other factors on it. The current study aimed to analyze the spatial and temporal distribution of (PM2.5) and its relationship with the meteorological parameters.(Air temperature, Relative humidity, Precipitation and wind speed) in Iraq during two periods (2001 and 2022). The dataset adopted in the study were downloaded from the Giovanni user interface which is based on satellite remote sensing data and reanalysis by MERRA-2model which produce by NASA. The output results shows that, the seasonal and annual PM2.5 concentration values increased from 2001 to 2022 due especially in the center and south of Iraq with the highest values of PM2.5 concentration recorded in the summers of 2001 and 2022 being 172.41 micro.g/m3 and 190.06 micro.g/m3 (increased 10.24%), respectively. Because of the low average temperature and the influence of northeasterly winds bringing continental air from Central Asia, PM2.5 values in northern and northeastern Iraq are lower than those in the center and southern regions. in 2001, they ranged from 8.41 to 12.6 micro.g/m3, whereas in 2022, they ranged from 9.02 to 15.98 micro.g/m3 throughout the year. Rainfall during the cold months in the north and northeast is an essential factor in cleaning the air of PM2.5. Also, study results indicate that the max. of PM 2.5 values have consistently exceeded the upper limits of PM2.5 quarterly standards set by both the US and Iraqi regulations, for the years 2001 and 2022, but the min. PM2.5 values are within both standards.


title: NiOx passivation in perovskite solar cells: from surface reactivity to device performance

abstract: Non-stoichiometric nickel oxide (NiOx) is the only metal oxide successfully used as hole transport material in p-i-n type perovskite solar cells (PSCs). Its favorable opto-electronic properties and facile large-scale preparation methods are potentially relevant for future commercialization of PSCs, though currently low operational stability of PSCs containing NiOx hole transport layers are reported. Poorly understood degradation reactions at the interface to the perovskite are seen as cause for the inferior stability and a variety of interface passivation approaches have been shown to be effective in improving the overall solar cell performance. To gain a better understanding of the processes happening at this interface, we systematically passivated possible specific defects on NiOx with three different categories of organic/inorganic compounds. The effects on the NiOx and the perovskite (MAPbI3) were investigated using x-ray photoelectron spectroscopy (XPS), X-ray diffraction (XRD), and scanning electron microscopy (SEM) where we find that the structural stability and film formation can be significantly affected. In combination with Density Functional Theory (DFT) calculations, a likely origin of NiOx-perovskite degradation interactions is proposed. The surface passivated NiOx was incorporated into MAPbI3 based PSCs and its influence on overall performance, particularly operational stability, was investigated by current-voltage (J-V), impedance spectroscopy (IS), and open circuit voltage decay (OCVD) measurements. Interestingly, we find that a superior structural stability due to an interface passivation must not relate to high operational stability. The discrepancy comes from the formation of excess ions at the interface which negatively impacts all solar cell parameters.


title: Random minimal surfaces in spheres

abstract: For any dimension $n>1$, we construct a branched minimal immersion $\psi_n$ from a closed Riemann surface $\Sigma_n$ to the round $n$-sphere of radius $\sqrt{8}$, such that if $\Sigma_n$ is endowed with the pullback metric and if $K$ is its Gaussian curvature, then $\Sigma_n$ is almost hyperbolic in the sense that $$\lim_{n\to \infty} \frac{1}{\mathrm{Area}(\Sigma_n)}\int_{\Sigma_n} |K+1|=0$$ and $\Sigma_n$ Benjamini-Schramm converges to the hyperbolic plane. Our proof is based on a connection between minimal surface theory and random matrix theory. The maps $\psi_n$ are obtained by applying the spherical Plateau problem to random unitary representations $\rho_N$ of the free group $F_2$.


title: Quantum effects in gravity beyond the Newton potential from a delocalised quantum source

abstract: Recent progress in table-top experiments offers the opportunity to show for the first time that gravity is not compatible with a classical description. In all current experimental proposals, such as the generation of gravitationally induced entanglement between two quantum sources of gravity, gravitational effects can be explained with the Newton potential, namely in a regime that is consistent with the weak-field limit of general relativity and does not probe the field nature of gravity. Hence, the Newtonian origin of the effects is a limitation to the conclusions on the nature of gravity that can be drawn from these experiments. Here, we identify two effects that overcome this limitation: they cannot be reproduced using the Newton potential and are independent of graviton emission. First, we show that the interaction between a generic quantum source of gravity, e.g. in a wide Gaussian state, and a test particle cannot be reproduced with the Newton potential nor with a known classical theory or gravity. Hence, observing the form of this interaction would require either a modification to classical gravity or its quantum description. Second, we show that the quantum commutator between the gravitational field and its canonically conjugate momentum appears as an additional term in the relative phase of a generic quantum source interacting with a test particle. Observing this term in the phase would be a test of the gravitational field as a quantum mediator. Identifying stronger quantum aspects of gravity than those reproducible with the Newton potential is crucial to prove the nonclassicality of the gravitational field and to plan a new generation of experiments testing quantum aspects of gravity in a broader sense than what proposed so far.


title: Thompson Sampling in Partially Observable Contextual Bandits

abstract: Contextual bandits constitute a classical framework for decision-making under uncertainty. In this setting, the goal is to learn the arms of highest reward subject to contextual information, while the unknown reward parameters of each arm need to be learned by experimenting that specific arm. Accordingly, a fundamental problem is that of balancing exploration (i.e., pulling different arms to learn their parameters), versus exploitation (i.e., pulling the best arms to gain reward). To study this problem, the existing literature mostly considers perfectly observed contexts. However, the setting of partial context observations remains unexplored to date, despite being theoretically more general and practically more versatile. We study bandit policies for learning to select optimal arms based on the data of observations, which are noisy linear functions of the unobserved context vectors. Our theoretical analysis shows that the Thompson sampling policy successfully balances exploration and exploitation. Specifically, we establish the followings: (i) regret bounds that grow poly-logarithmically with time, (ii) square-root consistency of parameter estimation, and (iii) scaling of the regret with other quantities including dimensions and number of arms. Extensive numerical experiments with both real and synthetic data are presented as well, corroborating the efficacy of Thompson sampling. To establish the results, we introduce novel martingale techniques and concentration inequalities to address partially observed dependent random variables generated from unspecified distributions, and also leverage problem-dependent information to sharpen probabilistic bounds for time-varying suboptimality gaps. These techniques pave the road towards studying other decision-making problems with contextual information as well as partial observations.


title: Experiments with Encoding Structured Data for Neural Networks

abstract: The project's aim is to create an AI agent capable of selecting good actions in a game-playing domain called Battlespace. Sequential domains like Battlespace are important testbeds for planning problems, as such, the Department of Defense uses such domains for wargaming exercises. The agents we developed combine Monte Carlo Tree Search (MCTS) and Deep Q-Network (DQN) techniques in an effort to navigate the game environment, avoid obstacles, interact with adversaries, and capture the flag. This paper will focus on the encoding techniques we explored to present complex structured data stored in a Python class, a necessary precursor to an agent.


title: An Evaluation of Real-time Adaptive Sampling Change Point Detection Algorithm using KCUSUM

abstract: Detecting abrupt changes in real-time data streams from scientific simulations presents a challenging task, demanding the deployment of accurate and efficient algorithms. Identifying change points in live data stream involves continuous scrutiny of incoming observations for deviations in their statistical characteristics, particularly in high-volume data scenarios. Maintaining a balance between sudden change detection and minimizing false alarms is vital. Many existing algorithms for this purpose rely on known probability distributions, limiting their feasibility. In this study, we introduce the Kernel-based Cumulative Sum (KCUSUM) algorithm, a non-parametric extension of the traditional Cumulative Sum (CUSUM) method, which has gained prominence for its efficacy in online change point detection under less restrictive conditions. KCUSUM splits itself by comparing incoming samples directly with reference samples and computes a statistic grounded in the Maximum Mean Discrepancy (MMD) non-parametric framework. This approach extends KCUSUM's pertinence to scenarios where only reference samples are available, such as atomic trajectories of proteins in vacuum, facilitating the detection of deviations from the reference sample without prior knowledge of the data's underlying distribution. Furthermore, by harnessing MMD's inherent random-walk structure, we can theoretically analyze KCUSUM's performance across various use cases, including metrics like expected delay and mean runtime to false alarms. Finally, we discuss real-world use cases from scientific simulations such as NWChem CODAR and protein folding data, demonstrating KCUSUM's practical effectiveness in online change point detection.


title: Standing spin waves in Permalloy-NiO bilayers as a probe of the interfacial exchange coupling

abstract: Ferromagnetic/Antiferromagnetic (FM/AFM) bilayers dynamics have been a recent topic of interest due to the interaction occurring at the interface, where the magnetic moments of the AFM can be imprinted into the FM, and the exchange bias field can affect these dynamics. Here, we investigate Permalloy (Py) and NiO (Py/NiO) hybrids and for comparison single Py films in the broad Py thickness range varied from few nm to 200 nm by using static (Kerr effect) and dynamic (spin waves) measurements along with micromagnetic simulations. We observe hybrid modes between uniform (ferromagnetic resonance FMR, n=0) and perpendicular standing spin waves (PSSWs, n=1, 2) and a clear enhancement of the PSSWs modes frequencies upon interfacing Py with NiO both from experiments and simulations. This enhancement becomes less pronounced as the thickness of the film increases, demonstrating its interfacial origin rooted in the exchange coupling between the FM and AFM layers. Besides, through micromagnetic simulations, we investigate and correlate changes in spatial profiles of the PSSWs with the interfacial exchange coupling. As the thickness is increased, we see that the n=1 and n=2 modes begin to couple with the fundamental FMR mode, resulting in asymmetric (with respect the Py layer center) modes. Our results suggest that PSSWs detection in a ferromagnet offers a means of probing the interfacial exchange coupling with the adjacent AFM layer. Furthermore, the controlled spatial symmetry breaking by the AFM layer enables engineering of PSSWs with different spatial profiles in the FM.


title: Parallel Play Saves Quantifiers

abstract: The number of quantifiers needed to express first-order properties is captured by two-player combinatorial games called multi-structural (MS) games. We play these games on linear orders and strings, and introduce a technique we call "parallel play", that dramatically reduces the number of quantifiers needed in many cases. Linear orders and strings are the most basic representatives of ordered structures -- a class of structures that has historically been notoriously difficult to analyze. Yet, in this paper, we provide upper bounds on the number of quantifiers needed to characterize different-sized subsets of these structures, and prove that they are tight up to constant factors, including, in some cases, up to a factor of $1+\varepsilon$, for arbitrarily small $\varepsilon$.


title: LAVE: LLM-Powered Agent Assistance and Language Augmentation for Video Editing

abstract: Video creation has become increasingly popular, yet the expertise and effort required for editing often pose barriers to beginners. In this paper, we explore the integration of large language models (LLMs) into the video editing workflow to reduce these barriers. Our design vision is embodied in LAVE, a novel system that provides LLM-powered agent assistance and language-augmented editing features. LAVE automatically generates language descriptions for the user's footage, serving as the foundation for enabling the LLM to process videos and assist in editing tasks. When the user provides editing objectives, the agent plans and executes relevant actions to fulfill them. Moreover, LAVE allows users to edit videos through either the agent or direct UI manipulation, providing flexibility and enabling manual refinement of agent actions. Our user study, which included eight participants ranging from novices to proficient editors, demonstrated LAVE's effectiveness. The results also shed light on user perceptions of the proposed LLM-assisted editing paradigm and its impact on users' creativity and sense of co-creation. Based on these findings, we propose design implications to inform the future development of agent-assisted content editing.


title: Quasicrystals as an intermediate form of matter between crystalline and amorphous solids

abstract: Quasicrystals have been observed in a variety of materials ranging from metal alloys to block copolymers and represent an "intermediate" form of matter between crystals and amorphous materials (glasses and liquids) in that their structural and dynamical properties can not readily described in terms of conventional solid-state models of liquids and solids. In the present work, we present a comprehensive analysis of basic thermodynamic and dynamic properties of quasicrystals to better understand the nature of the atomic motion underlying diffusion and structural relaxation in these materials. As our model system, we investigate a dodecagonal quasicrystal using molecular dynamics (MD) simulations in two dimensions (2D), subject to periodic boundary conditions. We observe a two-stage relaxation dynamics in the self-intermediate scattering function $F_s(k,t)$ of our quasicrystal material involving a fast $\beta$-relaxation on a ps timescale and $\alpha$ relaxation process having a highly temperature dependent relaxation time whose activation energy varies in concert with the extent of string-like collective motion. Multi-step relaxation of the intermediate scattering function and string-like collective atomic motion have similarly been observed ubiquitously in glass-forming liquids at low temperatures and in crystalline materials at elevated temperatures where structural relaxation and diffusion are both non-Arrhenius. After examining the dynamics of our quasi-crystalline material in great detail, we conclude that its dynamics more closely resemble observations on metallic glass-forming liquids, in qualitative accord with previous neutron scattering studies.


title: On very regular representations in presence of index

abstract: We discuss the possibility of very regular subgroups of a Lie group, in presence of an index figure. Further, representations that reduce action to a very regular boundary.


title: Spinorial Yamabe-type equations and the B\"ar-Hijazi-Lott invariant

abstract: We consider on a closed Riemannian spin manifold $(M^n,g,\sigma)$ the spinorial Yamabe type equation $D_g\varphi=\lambda|\varphi|^{\frac{2}{n-1}}\varphi$, where $\varphi$ is a spinor field and $\lambda$ is a positive constant. For a normalized solution $\varphi$ of this equation we find a positive lower bound for $\lambda^2$. As an application we obtain an explicit lower bound of the B\"ar-Hijazi-Lott invariant for some spin manifolds with positive scalar curvature.


title: Streaming algorithm for balance gain and cost with cardinality constraint on the integer lattice

abstract: Team formation problem is a very important problem in the labor market, and it is proved to be NP-hard. In this paper, we design an efficient bicriteria streaming algorithms to construct a balance between gain and cost in a team formation problem with cardinality constraint on the integer lattice. To solve this problem, we establish a model for maximizing the difference between a nonnegative normalized monotone submodule function and a nonnegative linear function. Further, we discuss the case where the first function of the object function is $\alpha$--weakly submodular. Combining the lattice binary search with the threshold method, we present an online algorithm called bicriteria streaming algorithms. Meanwhile, we give detailed analysis for both of these models.


title: Stability for the 3D Riemannian Penrose inequality

abstract: We show that Schwarzschild manifold is stable for the 3-dimensional Riemannian Penrose inequality in the pointed measured Gromov-Hausdorff topology modulo negligible spikes and boundary area perturbations.


title: Early Warning Signals for Bifurcations Embedded in High Dimensions

abstract: Recent work has highlighted the utility of methods for early warning signal detection in dynamic systems approaching critical tipping thresholds. Often these tipping points resemble local bifurcations, whose low dimensional dynamics can play out on a manifold embedded in a much higher dimensional state space. In many cases of practical relevance, the form of this embedding is poorly understood or entirely unknown. This paper explores how measurement of the critical phenomena that generically precede such bifurcations can be used to make inferences about the properties of their embeddings, and, conversely, how prior knowledge about the mechanism of bifurcation can robustify predictions of an oncoming tipping event. These modes of analysis are first demonstrated on a simple fluid flow system undergoing a Hopf bifurcation. The same approach is then applied to data associated with the West African monsoon shift, with results corroborated by existing models of the same system. This example highlights the effectiveness of the methodology even when applied to complex climate data, and demonstrates how a well-resolved spatial structure associated with the onset of atmospheric instability can be inferred purely from time series measurements.


title: Mutation of {\tau}-exceptional pairs and sequences

abstract: We introduce a notion of mutation for $\tau$-exceptional sequences of modules over arbitrary finite dimensional algebras. For hereditary algebras, we show that this coincides with the classical mutation of exceptional sequences. For rank two algebras, we show that mutation of $\tau$-exceptional sequences is transitive if and only if mutation of support $\tau$-tilting modules in the sense of Adachi-Iyama-Reiten is transitive.


title: How to Discern Important Urgent News?

abstract: We found that a simple property of clusters in a clustered dataset of news correlate strongly with importance and urgency of news (IUN) as assessed by LLM. We verified our finding across different news datasets, dataset sizes, clustering algorithms and embeddings. The found correlation should allow using clustering (as an alternative to LLM) for identifying the most important urgent news, or for filtering out unimportant articles.


title: Spontaneous Emission in the presence of Quantum Mirrors

abstract: Arrays of atoms coupled to waveguides can behave as mirrors. We consider an array of $\Lambda$-type three-level atoms wherein preparing the atoms in one ground state or another leads to reflection or transmission of the guided electromagnetic field; a superposition of the two ground states thus corresponds to a coherent superposition of mirror-like and transparent boundary conditions. We analyze the spontaneous emission of an excited two-level atom in the presence of such a quantum mirror, and inside a cavity formed by quantum mirrors, demonstrating that the resulting dynamics of the excited atom can exhibit exotic features, e.g., a superposition of Rabi cycle and exponential decay. Our results pave the way for exploring quantum electrodynamics (QED) phenomena in a paradigm wherein boundary conditions can exhibit quantum superpositions and correlations.


title: Uncertainty analysis of the plasma impedance probe

abstract: A plasma impedance probe (PIP) is a type of in-situ, radio-frequency (RF) probe that is traditionally used to measure plasma properties (e.g. density) in low-density environments such as the Earth's ionosphere. We believe that PIPs are underrepresented in laboratory settings, in part because PIP operation and analysis has not been optimized for signal-to-noise ratio (SNR), reducing the probe's accuracy, upper density limit, and acquisition rate. This work presents our efforts in streamlining and simplifying the PIP design, model, calibration, and analysis for unmagnetized laboratory plasmas, in both continuous and pulsed PIP operation. The focus of this work is a Monte Carlo uncertainty analysis, which identifies operational and analysis procedures that improve SNR by multiple orders of magnitude. Additionally, this analysis provides evidence that the sheath resonance (and not the plasma frequency as previously believed) sets the PIP's upper density limit, which likely provides an additional method for extending the PIP's density limit.


title: Effective module lattices and their shortest vectors

abstract: We prove tight probabilistic bounds for the shortest vectors in module lattices over number fields using the results of arXiv:2308.15275. Moreover, establishing asymptotic formulae for counts of fixed rank matrices with algebraic integer entries and bounded Euclidean length, we prove an approximate Rogers integral formula for discrete sets of module lattices obtained from lifts of algebraic codes. This in turn implies that the moment estimates of arXiv:2308.15275 as well as the aforementioned bounds on the shortest vector also carry through for large enough discrete sets of module lattices.


title: Linear and nonlinear receptivity mechanisms in boundary layers subject to free-stream turbulence

abstract: Large-eddy simulations of a flat-plate boundary layer, without a leading edge, subject to multiple levels of incoming free stream turbulence are considered in the present work. Within an input-output model where non-linear terms of the incompressible Navier-Stokes equations are treated as an external forcing, we manage to separate inputs related to perturbations coming through the intake of the numerical domain, whose evolution represents a linear mechanism, and the volumetric non-linear forcing due to triadic interactions. With these, we perform the full reconstruction of the statistics of the flow, as measured in the simulations, to quantify pairs of wavenumbers and frequencies more affected by either linear or non-linear receptivity mechanisms. Inside the boundary layer, different wavenumbers at near-zero frequency reveal streaky structures. Those that are amplified predominantly via linear interactions with the incoming vorticity occur upstream and display transient growth, while those generated by the non-linear forcing are the most energetic and appear in more downstream positions. The latter feature vortices growing proportionally to the laminar boundary layer thickness, along with a velocity profile that agrees with the optimal amplification obtained by linear transient growth theory. The numerical approach presented is general and could potentially be extended to any simulation for which receptivity to incoming perturbations needs to be assessed.


title: A New Radio to Overcome Critical Link Budgets

abstract: We propose Multi-Antenna (MA) Towards Inband Shift Keying (TISK): a new multi-carrier radio concept to cope with critical link budgets. In contrast to common proposals that rely on analog beamforming at both transmitter and receiver, MA-TISK does not require beam alignment. The transmitted signals have all constant envelope in continuous time, which allows for efficient, low-cost power amplification and up-conversion. The concept is compatible with any linear PSK-modulation as well as pulse position modulation. Each sub-carrier is sent over a separate antenna that is equipped with a voltage-controlled oscillator. The phases of these oscillators are controlled by digital baseband. Temporal signal combining makes up for the lack of beamforming gain at the transmitter. A common message may be broadcast to many receivers, simultaneously. Demodulation can be efficiently implemented by means of fast Fourier transform. MA-TISK does not suffer from spectral re-growth issues plaguing other constant envelope modulations like GMSK. Almost rectangular signal spectra similar to those for linear modulation with root-raised-cosine pulse shaping are possible. For the 100 MHz-wide spectral mask of 5G downlink, QPSK-modulation allows for 160 MBit/s with 5.74 MHz subcarrier spacing when using 16 transmit antennas. The wide carrier spacing makes the signals insensitive to Doppler effects. There is no loss in link budget gain compared to spatial beamforming at the transmitter.


title: Annihilation dynamics during spiral defect chaos revealed by particle models

abstract: Pair-annihilation events are ubiquitous in a variety of spatially extended systems and are often studied using computationally expensive simulations. Here we develop an approach in which we simulate the pair-annihilation of spiral wave tips in cardiac models using a computationally efficient particle model. Spiral wave tips are represented as particles with dynamics governed by diffusive behavior and short-ranged attraction. The parameters for diffusion and attraction are obtained by comparing particle motion to the trajectories of spiral wave tips in cardiac models during spiral defect chaos. The particle model reproduces the annihilation rates of the cardiac models and can determine the statistics of spiral wave dynamics, including its mean termination time. We show that increasing the attraction coefficient sharply decreases the mean termination time, making it a possible target for pharmaceutical intervention


title: Discrete Probabilistic Inference as Control in Multi-path Environments

abstract: We consider the problem of sampling from a discrete and structured distribution as a sequential decision problem, where the objective is to find a stochastic policy such that objects are sampled at the end of this sequential process proportionally to some predefined reward. While we could use maximum entropy Reinforcement Learning (MaxEnt RL) to solve this problem for some distributions, it has been shown that in general, the distribution over states induced by the optimal policy may be biased in cases where there are multiple ways to generate the same object. To address this issue, Generative Flow Networks (GFlowNets) learn a stochastic policy that samples objects proportionally to their reward by approximately enforcing a conservation of flows across the whole Markov Decision Process (MDP). In this paper, we extend recent methods correcting the reward in order to guarantee that the marginal distribution induced by the optimal MaxEnt RL policy is proportional to the original reward, regardless of the structure of the underlying MDP. We also prove that some flow-matching objectives found in the GFlowNet literature are in fact equivalent to well-established MaxEnt RL algorithms with a corrected reward. Finally, we study empirically the performance of multiple MaxEnt RL and GFlowNet algorithms on multiple problems involving sampling from discrete distributions.


title: Interpretable Generative Adversarial Imitation Learning

abstract: Imitation learning methods have demonstrated considerable success in teaching autonomous systems complex tasks through expert demonstrations. However, a limitation of these methods is their lack of interpretability, particularly in understanding the specific task the learning agent aims to accomplish. In this paper, we propose a novel imitation learning method that combines Signal Temporal Logic (STL) inference and control synthesis, enabling the explicit representation of the task as an STL formula. This approach not only provides a clear understanding of the task but also allows for the incorporation of human knowledge and adaptation to new scenarios through manual adjustments of the STL formulae. Additionally, we employ a Generative Adversarial Network (GAN)-inspired training approach for both the inference and the control policy, effectively narrowing the gap between the expert and learned policies. The effectiveness of our algorithm is demonstrated through two case studies, showcasing its practical applicability and adaptability.


title: The optimal placement of the head in the noun phrase. The case of demonstrative, numeral, adjective and noun

abstract: The word order of a sentence is shaped by multiple principles. The principle of syntactic dependency distance minimization is in conflict with the principle of surprisal minimization (or predictability maximization) in single head syntactic dependency structures: while the former predicts that the head should be placed at the center of the linear arrangement, the latter predicts that the head should be placed at one of the ends (either first or last). A critical question is when surprisal minimization (or predictability maximization) should surpass syntactic dependency distance minimization. In the context of single head structures, it has been predicted that this is more likely to happen when two conditions are met, i.e. (a) fewer words are involved and (b) words are shorter. Here we test the prediction on the noun phrase when it is composed of a demonstrative, a numeral, an adjective and a noun. We find that, across preferred orders in languages, the noun tends to be placed at one of the ends, confirming the theoretical prediction. We also show evidence of anti locality effects: syntactic dependency distances in preferred orders are longer than expected by chance.


title: Towards Tight Convex Relaxations for Contact-Rich Manipulation

abstract: We present a method for global motion planning of robotic systems that interact with the environment through contacts. Our method directly handles the hybrid nature of such tasks using tools from convex optimization. We formulate the motion-planning problem as a shortest-path problem in a graph of convex sets, where a path in the graph corresponds to a contact sequence and a convex set models the quasi-static dynamics within a fixed contact mode. For each contact mode, we use semidefinite programming to relax the nonconvex dynamics that results from the simultaneous optimization of the object's pose, contact locations, and contact forces. The result is a tight convex relaxation of the overall planning problem, that can be efficiently solved and quickly rounded to find a feasible contact-rich trajectory. As a first application of this technique, we focus on the task of planar pushing. Exhaustive experiments show that our convex-optimization method generates plans that are consistently within a small percentage of the global optimum. We demonstrate the quality of these plans on a real robotic system.


title: Gravitational signatures of ALP dark matter fragmentation

abstract: The misalignment mechanism for axion-like particles (ALPs) is a leading explanation for dark matter. In this work we investigate ALPs with non-periodic potentials, which allow for large misalignment of the field from the minimum and make it possible for ALPs to match the relic density of dark matter in a large part of the parameter space. Such potentials give rise to self-interactions which can trigger an exponential growth of fluctuations in the ALP field via parametric resonance, leading to the fragmentation of the field. The fluctuations later collapse to halos that can be dense enough to produce observable gravitational effects. These effects would provide a probe of dark matter even if it does not couple to the Standard Model (or too feebly). We determine the relevant regions of parameter space in the (ALP mass, decay constant)-plane and compare predictions in different axion fragmentation models. These proceedings are a short version of arXiv:2305.03756


title: Weighted Brunn-Minkowski Theory II: Inequalities for Mixed Measures and Applications

abstract: In "Weighted Brunn-Minkowski Theory I", the prequel to this work, we discussed how recent developments on concavity of measures have laid the foundations of a nascent weighted Brunn-Minkowski theory. In particular, we defined the mixed measures of three convex bodies and obtained its integral representation. In this work, we obtain inequalities for mixed measures, such as a generalization of Fenchel's inequality; this provides a new, simpler proof of the classical volume case. Moreover, we show that mixed measures are connected to the study of log-submodularity and supermodularity of the measure of Minkowski sums of convex bodies. This elaborates on the recent investigations of these properties for the Lebesgue measure. We conclude by establishing that the only Radon measures that are supermodular over the class of compact, convex sets are multiples of the Lebesgue measure. Motivated by this result, we then discuss weaker forms of supermodularity by restricting the class of convex sets.


title: A variable ionized disk wind in MAXI J1803-298 revealed by NICER

abstract: We present the results from the NICER observation data of MAXI J1803-298 across the entire 2021 outburst. In the intermediate and soft state, we detect significant absorption lines at $\sim 7.0$ keV and $\sim 6.7$ keV, arising from the X-ray disk wind outflowing with a velocity of hundreds of km per second along our line of sight. The fitting results from photoionized model suggest that the wind is driven by thermal pressure and the mass-loss rate is low. We find a clear transition for iron from predominantly H-like to predominantly He-like during the intermediate-to-soft state transition. Our results indicate this transition for iron is caused by the evolution of the illuminating spectrum and the slow change of the geometric properties of the disk wind together. The coexistence of disk wind and QPOs features in intermediate state is also reported. Our study makes MAXI J1803-298 the first source in which a transition from optical wind to X-ray wind is detected, offering new insights into the evolution of disk winds across an entire outburst and long-term coupling of accretion disks and mass outflows around accreting black holes.


title: Higher relative index theorems for foliations

abstract: In this paper we solve the general case of the cohomological relative index problem for foliations of non-compact manifolds. In particular, we significantly generalize the groundbreaking results of Gromov and Lawson, [GL83], to Dirac operators defined along the leaves of foliations of non-compact complete Riemannian manifolds, by involving all the terms of the Connes-Chern character, especially the higher order terms in Haefliger cohomology. The zero-th order term corresponding to holonomy invariant measures was carried out in [BH21] and becomes a special case of our main results here. In particular, for two leafwise Dirac operators on two foliated manifolds which agree near infinity, we define a relative topological index and the Connes-Chern character of a relative analytic index, both being in relative Haefliger cohomology. We show that these are equal. This invariant can be paired with closed holonomy invariant currents (which agree near infinity) to produce higher relative scalar invariants. When we relate these invariants to the leafwise index bundles, we restrict to Riemannian foliations on manifolds of sub-exponential growth. This allows us {to prove a higher relative index bundle theorem}, extending the classical index bundle theorem of [BH08]. Finally, we construct examples of foliations and use these invariants to prove that their spaces of leafwise positive scalar curvature metrics have infinitely many path-connected components, completely new results which are not available from [BH21]. In particular, these results confirm the well-known idea that important geometric information of foliations is embodied in the higher terms of the A-hat genus.


title: Two-component integrable extension of general heavenly equation

abstract: We introduce an integrable two-component extension of the general heavenly equation and prove that the solutions of this extension are in one-to-one correspondence with 4-dimensional hyper-para-Hermitian metrics. Furthermore, we demonstrate that if the metrics in question are hyper-para-K\"ahler, then our system reduces to the general heavenly equation. We also present an infinite hierarchy of nonlocal symmetries and a recursion operator for the system under study.


title: Multi-Antenna Towards Inband Shift Keying

abstract: We propose a new continuous phase frequency shift keying that is particularly suited for multi-antenna communications when the link budget is critical and beam alignment is problematic. It combines the constant envelope of frequency modulation with low-rate repetition coding in order to compensate for the absence of transmit beamforming. Although it is a frequency modulation, its transmit signal shows close to rectangular spectral shape. Similar to GSM's Gaussian minimum shift keying, it can be well approximated by linear modulation, when combined with differential precoding. This allows for easy coherent demodulation by means of a windowed fast Fourier transform.


title: Modeling Blood Alcohol Concentration Using Fractional Differential Equations Based on the $\psi$-Caputo Derivative

abstract: We propose a novel dynamical model for blood alcohol concentration that incorporates $\psi$-Caputo fractional derivatives. Using the generalized Laplace transform technique, we successfully derive an analytic solution for both the alcohol concentration in the stomach and the alcohol concentration in the blood of an individual. These analytical formulas provide us a straightforward numerical scheme, which demonstrates the efficacy of the $\psi$-Caputo derivative operator in achieving a better fit to real experimental data on blood alcohol levels available in the literature. In comparison to existing classical and fractional models found in the literature, our model outperforms them significantly. Indeed, by employing a simple yet non-standard kernel function $\psi(t)$, we are able to reduce the error by more than half, resulting in an impressive gain improvement of 59 percent.


title: Entanglement degradation under local dissipative Landau-Zener noise

abstract: We study entanglement degradation when noise on one share of an entangled pair is described by the dissipative Landau-Zener model. We show that spin-coupling direction to the environment significantly affects entanglement dynamics. In particular, for zero bath temperature in the slow-driving regime with transversal coupling, entanglement remains intact and in the fast-driving regime transversal noise have less destructive affects on entanglement compared to the longitudinal noise. Furthermore, we show that non-adiabatic dynamic is more in favour of preserving entanglement compared to adiabatic evolution.


title: LaserSAM: Zero-Shot Change Detection Using Visual Segmentation of Spinning LiDAR

abstract: This paper presents an approach for applying camera perception techniques to spinning LiDAR data. To improve the robustness of long-term change detection from a 3D LiDAR, range and intensity information are rendered into virtual perspectives using a pinhole camera model. Hue-saturation-value image encoding is used to colourize the images by range and near-IR intensity. The LiDAR's active scene illumination makes it invariant to ambient brightness, which enables night-to-day change detection without additional processing. Using the range-colourized, perspective image allows existing foundation models to detect semantic regions. Specifically, the Segment Anything Model detects semantically similar regions in both a previously acquired map and live view from a path-repeating robot. By comparing the masks in both views, changes in the live scan are detected. Results indicate that the Segment Anything Model accurately captures the shape of arbitrary changes introduced into scenes. The proposed method achieves a segmentation intersection over union of 73.3% when evaluated in unstructured environments and 80.4% when evaluated within the planning corridor. Changes can be detected reliably through day-to-night illumination variations. After pixel-level masks are generated, the one-to-one correspondence with 3D points means that the 2D masks can be used directly to recover the 3D location of the changes. The detected 3D changes are avoided in a closed loop by treating them as obstacles in a local motion planner. Experiments on an unmanned ground vehicle demonstrate the performance of the method.


title: Maximum Likelihood Degrees of Brownian Motion Tree Models: Star Trees and Root Invariance

abstract: A Brownian motion tree (BMT) model is a Gaussian model whose associated set of covariance matrices is linearly constrained according to common ancestry in a phylogenetic tree. We study the complexity of inferring the maximum likelihood (ML) estimator for a BMT model by computing its ML-degree. Our main result is that the ML-degree of the BMT model on a star tree with $n + 1$ leaves is $2^{n+1}-2n-3$, which was previously conjectured by Am\'endola and Zwiernik. We also prove that the ML-degree of a BMT model is independent of the choice of the root. The proofs rely on the toric geometry of concentration matrices in a BMT model. Toward this end, we produce a combinatorial formula for the determinant of the concentration matrix of a BMT model, which generalizes the Cayley-Pr\"ufer theorem to complete graphs with weights given by a tree.


title: A Computationally Efficient Learning-Based Model Predictive Control for Multirotors under Aerodynamic Disturbances

abstract: Neglecting complex aerodynamic effects hinders high-speed yet high-precision multirotor autonomy. In this paper, we present a computationally efficient learning-based model predictive controller that simultaneously optimizes a trajectory that can be tracked within the physical limits (on thrust and orientation) of the multirotor system despite unknown aerodynamic forces and adapts the control input. To do this, we leverage the well-known differential flatness property of multirotors, which allows us to transform their nonlinear dynamics into a linear model. The main limitation of current flatness-based planning and control approaches is that they often neglect dynamic feasibility. This is because these constraints are nonlinear as a result of the mapping between the input, i.e., multirotor thrust, and the flat state. In our approach, we learn a novel representation of the drag forces by learning the mapping from the flat state to the multirotor thrust vector (in a world frame) as a Gaussian Process (GP). Our proposed approach leverages the properties of GPs to develop a convex optimal controller that can be iteratively solved as a second-order cone program (SOCP). In simulation experiments, our proposed approach outperforms related model predictive controllers that do not account for aerodynamic effects on trajectory feasibility, leading to a reduction of up to 55% in absolute tracking error.


title: Hacktivism Goes Orbital: Investigating NB65's Breach of ROSCOSMOS

abstract: In March of 2022, Network battalion 65 (NB65), a hacktivist affiliate of Anonymous, publicly asserted its successful breach of ROSCOSMOS's satellite imaging capabilities in response to Russia's invasion of Ukraine. NB65 disseminated a series of primary sources as substantiation, proclaiming the incapacitation of ROSCOSMOS's space-based vehicle monitoring system and doxing of related proprietary documentation. Despite the profound implications of hacktivist incursions into the space sector, the event has garnered limited attention due to the obscurity of technical attack vectors and ROCOSMOS's denial of NB65's allegations. Through analysis of NB65's released primary sources of evidence, this paper uncovers the probable vulnerabilities and exploits that enabled the alleged breach into ROSCOSMOS's ground and space segment. Additionally, we highlight lessons learned and the consequences this event has for the global aerospace community.


title: The E302 instability-versus-efficiency experiment at FACET-II

abstract: We discuss plans for the E302 instability-efficiency experiment, starting in 2024 at the recently upgraded FACET-II facility at SLAC National Accelerator Laboratory. The beam-breakup instability will be the main area of study for the E302 experiment. With the imaging spectrometer at FACET-II, we introduce a novel technique for observing and quantifying the amplitude growth of the trailing bunch due to the transverse instability. Using the transverse position on the spectrometer screen and the transfer matrix of the magnetic lattice used for the spectrometer, we aim to extract a $x'$-$E$ charge distribution that can be used to quantify the amplitude of the beam. By varying the trailing bunch's charge and, hence, the beam loading of the accelerating field, we aim to adjust the wake-to-beam power transfer efficiency in the E302 experiment. We plan to quantify the amplitude for different configurations of the beam charge and, hence, investigate the relationship between the beam-breakup instability and efficiency. We use a combination of particle-in-cell (PIC) codes to simulate a beam-driven plasma wakefield accelerator from start-to-end with a FACET-II-like spectrometer and demonstrate the methodology that will be used for the instability studies at the E302 experiment.


abstract: The report describes the discussions from the Workshop on Mathematical Opportunities in Digital Twins (MATH-DT) from December 11-13, 2023, George Mason University. It illustrates that foundational Mathematical advances are required for Digital Twins (DTs) that are different from traditional approaches. A traditional model, in biology, physics, engineering or medicine, starts with a generic physical law (e.g., equations) and is often a simplification of reality. A DT starts with a specific ecosystem, object or person (e.g., personalized care) representing reality, requiring multi -scale, -physics modeling and coupling. Thus, these processes begin at opposite ends of the simulation and modeling pipeline, requiring different reliability criteria and uncertainty assessments. Additionally, unlike existing approaches, a DT assists humans to make decisions for the physical system, which (via sensors) in turn feeds data into the DT, and operates for the life of the physical system. While some of the foundational mathematical research can be done without a specific application context, one must also keep specific applications in mind for DTs. E.g., modeling a bridge or a biological system (a patient), or a socio-technical system (a city) is very different. The models range from differential equations (deterministic/uncertain) in engineering, to stochastic in biology, including agent-based. These are multi-scale hybrid models or large scale (multi-objective) optimization problems under uncertainty. There are no universal models or approaches. For e.g., Kalman filters for forecasting might work in engineering, but can fail in biomedical domain. Ad hoc studies, with limited systematic work, have shown that AI/ML methods can fail for simple engineering systems and can work well for biomedical problems. A list of `Mathematical Opportunities and Challenges' concludes the report.


title: Experimental study of cavitation development in liquid in pulsed non-uniform electric field under the action of ponderomotive forces

abstract: In this paper, the Rayleigh scattering method is used to study the formation of cavitation in water in a pulsed inhomogeneous electric field when a nanosecond high voltage pulse is applied to a needle-like electrode. The observational results confirm the theoretical picture of cavitation development under the action of electrostrictive ponderomotive forces. The values of the negative pressure, the average size of the cavitation nanovoids and their concentration were obtained from these measurements, which are in agreement with theoretical estimates.


title: Turbulent activity in the near-wall region of adverse pressure gradient turbulent boundary layers

abstract: Two direct numerical simulation (DNS) databases are investigated to understand the effect of the outer-layer turbulence on the inner layer's structures and energy transfer mechanisms. The first DNS database is the non-equilibrium adverse-pressure-gradient (APG) turbulence boundary layer (TBL) of Gungor et al. (2022). Its Reynolds number and the inner-layer pressure gradient parameter reach above 8000 and 10, respectively. The shape factor spans between 1.4 and 3.3, which indicates the flow has various velocity defect situations. The second database is the same flow as the first one but the outer layer turbulence is artificially eliminated in this flow. Turbulence is removed above 0.15 local boundary layer thickness. For the analysis, we chose four streamwise positions with small, moderate, large, and very-large velocity defect. We compare the wall-normal distribution of Reynolds stresses, two-point correlations and spectral distributions of energy, production and pressure strain. The results show that the inner layer turbulence can sustain itself when the outer-layer turbulence does not exist regardless of the velocity defect or the pressure gradient. The two-point correlations of both cases show that outer large-scale structures affect the inner layer structures significantly. The streamwise extent of the correlation contours scales with pressure-viscous units. This shows the importance of the pressure gradient's effect on the inner-layer structures. The spectral distributions demonstrate that the energy transfer mechanisms are probably the same in the inner layer regardless of the velocity defect, which suggests the near-wall cycle may exist even in very-large defect APG TBLs where the mean shear in the inner layer is considerably lower than small-defect APG TBLs.


title: Universal Manipulation Interface: In-The-Wild Robot Teaching Without In-The-Wild Robots

abstract: We present Universal Manipulation Interface (UMI) -- a data collection and policy learning framework that allows direct skill transfer from in-the-wild human demonstrations to deployable robot policies. UMI employs hand-held grippers coupled with careful interface design to enable portable, low-cost, and information-rich data collection for challenging bimanual and dynamic manipulation demonstrations. To facilitate deployable policy learning, UMI incorporates a carefully designed policy interface with inference-time latency matching and a relative-trajectory action representation. The resulting learned policies are hardware-agnostic and deployable across multiple robot platforms. Equipped with these features, UMI framework unlocks new robot manipulation capabilities, allowing zero-shot generalizable dynamic, bimanual, precise, and long-horizon behaviors, by only changing the training data for each task. We demonstrate UMI's versatility and efficacy with comprehensive real-world experiments, where policies learned via UMI zero-shot generalize to novel environments and objects when trained on diverse human demonstrations. UMI's hardware and software system is open-sourced at https://umi-gripper.github.io.


title: Quantum dynamics of superconductor-quantum dot-superconductor Josephson junctions

abstract: Josephson junctions constructed from superconductor-quantum dot-superconductor (S-QD-S) heterostructures have been used to realize a variety of voltage-tunable superconducting quantum devices, including qubits and parametric amplifiers. In such devices, the interplay between the charge degree of freedom associated with the quantum dot and its environment must be considered for faithful modeling of circuit dynamics. Here we describe the self-consistent quantization of a capacitively-shunted S-QD-S junction via path-integral formulation. In the effective Hamiltonian, the Josephson potential for the Andreev bound states reproduces earlier results for static phase bias, whereas the charging energy term has new features: (i) the system's capacitance is renormalized by the junction gate voltage, an effect which depends on the strength of the tunneling rates between the dot and its superconducting leads as well, and (ii) an additional charge offset appears for asymmetric junctions. These results are important to understand future experiments and quantum devices incorporating S-QD-S junctions in arbitrary impedance environments.


title: Rough path theory

abstract: The theory of rough paths arose from a desire to establish continuity properties of ordinary differential equations involving terms of low regularity. While essentially an analytic theory, its main motivation and applications are in stochastic analysis, where it has given a new perspective on It\^o calculus and a meaning to stochastic differential equations driven by irregular paths outside the setting of semi-martingales. In this survey, we present some of the main ideas that enter rough path theory. We discuss complementary notions of solutions for rough differential equations and the related notion of path signature, and give several applications and generalisations of the theory.


title: Jones-Wenzl projectors and the Khovanov homotopy of the infinite twist

abstract: We construct and study a lift of Jones-Wenzl projectors to the setting of Khovanov spectra, and provide a realization of such lifted projectors via a Cooper-Krushkal-like sequence of maps. We also give a polynomial action on the 3-strand spectral projector allowing a complete computation of the 3-colored Khovanov spectrum of the unknot, proving a conjecture of Lobb-Orson-Sch\"{u}tz. As a byproduct, we disprove a conjecture of Lawson-Lipshitz-Sarkar on the topological Hochschild homology of tangle spectra.


title: Chromatic symmetric functions and polynomial invariants of trees

abstract: Stanley asked whether a tree is determined up to isomorphism by its chromatic symmetric function. We approach Stanley's problem by studying the relationship between the chromatic symmetric function and other invariants. First, we prove Crew's conjecture that the chromatic symmetric function of a tree determines its generalized degree sequence, which enumerates vertex subsets by cardinality and the numbers of internal and external edges. Second, we prove that the restriction of the generalized degree sequence to subtrees contains exactly the same information as the subtree polynomial, which enumerates subtrees by cardinality and number of leaves. Third, we construct arbitrarily large families of trees sharing the same subtree polynomial, proving and generalizing a conjecture of Eisenstat and Gordon.


title: HI-GAN: Hierarchical Inpainting GAN with Auxiliary Inputs for Combined RGB and Depth Inpainting

abstract: Inpainting involves filling in missing pixels or areas in an image, a crucial technique employed in Mixed Reality environments for various applications, particularly in Diminished Reality (DR) where content is removed from a user's visual environment. Existing methods rely on digital replacement techniques which necessitate multiple cameras and incur high costs. AR devices and smartphones use ToF depth sensors to capture scene depth maps aligned with RGB images. Despite speed and affordability, ToF cameras create imperfect depth maps with missing pixels. To address the above challenges, we propose Hierarchical Inpainting GAN (HI-GAN), a novel approach comprising three GANs in a hierarchical fashion for RGBD inpainting. EdgeGAN and LabelGAN inpaint masked edge and segmentation label images respectively, while CombinedRGBD-GAN combines their latent representation outputs and performs RGB and Depth inpainting. Edge images and particularly segmentation label images as auxiliary inputs significantly enhance inpainting performance by complementary context and hierarchical optimization. We believe we make the first attempt to incorporate label images into inpainting process.Unlike previous approaches requiring multiple sequential models and separate outputs, our work operates in an end-to-end manner, training all three models simultaneously and hierarchically. Specifically, EdgeGAN and LabelGAN are first optimized separately and further optimized inside CombinedRGBD-GAN to enhance inpainting quality. Experiments demonstrate that HI-GAN works seamlessly and achieves overall superior performance compared with existing approaches.


title: Correlation Clustering with Vertex Splitting

abstract: We explore Cluster Editing and its generalization Correlation Clustering with a new operation called permissive vertex splitting which addresses finding overlapping clusters in the face of uncertain information. We determine that both problems are NP-hard, yet they exhibit significant differences in parameterized complexity and approximability. For Cluster Editing with Permissive Vertex Splitting, we show a polynomial kernel when parameterized by the solution size and develop a polynomial-time algorithm with approximation factor 7. In the case of Correlation Clustering, we establish para-NP-hardness when parameterized by solution size and demonstrate that computing an $n^{1-\epsilon}$-approximation is NP-hard for any constant $\epsilon > 0$. Additionally, we extend the established link between Correlation Clustering and Multicut to the setting with permissive vertex splitting.


title: Ultrafast Photochemistry and Electron Diffraction for Cyclobutanone in the S2 State: Surface Hopping with Time-Dependent Density Functional Theory

abstract: We simulate the photodynamics of gas-phase cyclobutanone excited to the S$_2$ state using fewest switches surface hopping (FSSH) dynamics powered by time-dependent density functional theory (TDDFT). We predict a total C3+C2 photoproduct yield of 9%, with a C3:C2 product ratio of 1:8. Two primary S$_2$$\rightarrow$S$_1$ conical intersections are identified: $\beta$ stretch and CCH bend, with the higher energy $\beta$ stretch being associated with sub-picosecond S$_2$ decay. Excited state lifetimes computed with respect to electronic state populations were found to be 7.0 ps (S$_2$$\rightarrow$S$_1$) and 550 fs (S$_1$$\rightarrow$S$_0$). We also generate time-resolved difference pair distribution functions ($\Delta$PDFs) from our TDDFT-FSSH dynamics results in order to generate direct comparisons to ultrafast electron diffraction experiment observables. Global and target analysis of time-resolved $\Delta$PDFs produced a distinct set of lifetimes: i) a 0.462 ps decay, and ii) a 16.8 ps decay that both resemble the S$_2$ minimum, as well as iii) a long ($>$ nanosecond) decay that resembles the S$_1$ minimum geometry and the fully separated C3/C2 products. Finally, we contextualize our results by considering the impact of the most likely sources of significant errors.


title: LoVoCCS -- II. Weak Lensing Mass Distributions, Red-Sequence Galaxy Distributions, and Their Alignment with the Brightest Cluster Galaxy in 58 Nearby X-ray-Luminous Galaxy Clusters

abstract: The Local Volume Complete Cluster Survey (LoVoCCS) is an on-going program to observe nearly a hundred low-redshift X-ray-luminous galaxy clusters (redshifts $0.03<z<0.12$ and X-ray luminosities in the 0.1-2.4 keV band $L_{X500c}>10^{44}$ erg/s) with the Dark Energy Camera (DECam), capturing data in $u,g,r,i,z$ bands with a $5\sigma$ point source depth of approximately 25-26th AB magnitudes. Here, we map the aperture masses in 58 galaxy cluster fields using weak gravitational lensing. These clusters span a variety of dynamical states, from nearly relaxed to merging systems, and approximately half of them have not been subject to detailed weak lensing analysis before. In each cluster field, we analyze the alignment between the 2D mass distribution described by the aperture mass map, the 2D red-sequence (RS) galaxy distribution, and the brightest cluster galaxy (BCG). We find that the orientations of the BCG and the RS distribution are strongly aligned throughout the interiors of the clusters: the median misalignment angle is 19 deg within 2 Mpc. We also observe the alignment between the orientations of the RS distribution and the overall cluster mass distribution (by a median difference of 32 deg within 1 Mpc), although this is constrained by galaxy shape noise and the limitations of our cluster sample size. These types of alignment suggest long-term dynamical evolution within the clusters over cosmic timescales.


title: Athermal granular creep in a quenched sandpile

abstract: Creep is a generic descriptor of slow motions -- in the context of materials, it describes quasi-static deformation of a solid when subjected to stresses below the global yield, at which all rigidity collapses and the material flows. Here, we experimentally investigate creep, flow, and the transition between the two states in a granular heap flow. Within the surface flowing layer the dimensionless strain rate diminishes with depth, there is an absence of spatial correlations, and there is no aging dynamics. Beneath this layer, the bulk creeps via localized avalanches of plasticity, and there is significant aging. The transition between fast surface flow and slow bulk creep and aging is observed to be in the vicinity of a critical inertial number of $I = 10^{-5}$. Surprisingly, at the cessation of surface flow and the `quenching' of the pile, creep persists in the absence of the flowing layer; albeit with significant differences for a pile that experiences a long duration of surface flow (strongly annealed) and one where flow during preparation does not last long (weakly annealed). Our results contribute to an emerging view of athermal granular creep, showing similarities across dry and submerged systems. Quenched quiescent heaps that creep indefinitely, however, present a challenge to granular rheology, and open new possibilities for interpreting and casting creep and deformation of soils in nature.


title: What to Do When Your Discrete Optimization Is the Size of a Neural Network?

abstract: Oftentimes, machine learning applications using neural networks involve solving discrete optimization problems, such as in pruning, parameter-isolation-based continual learning and training of binary networks. Still, these discrete problems are combinatorial in nature and are also not amenable to gradient-based optimization. Additionally, classical approaches used in discrete settings do not scale well to large neural networks, forcing scientists and empiricists to rely on alternative methods. Among these, two main distinct sources of top-down information can be used to lead the model to good solutions: (1) extrapolating gradient information from points outside of the solution set (2) comparing evaluations between members of a subset of the valid solutions. We take continuation path (CP) methods to represent using purely the former and Monte Carlo (MC) methods to represent the latter, while also noting that some hybrid methods combine the two. The main goal of this work is to compare both approaches. For that purpose, we first overview the two classes while also discussing some of their drawbacks analytically. Then, on the experimental section, we compare their performance, starting with smaller microworld experiments, which allow more fine-grained control of problem variables, and gradually moving towards larger problems, including neural network regression and neural network pruning for image classification, where we additionally compare against magnitude-based pruning.


title: Highlighting the Safety Concerns of Deploying LLMs/VLMs in Robotics

abstract: In this paper, we highlight the critical issues of robustness and safety associated with integrating large language models (LLMs) and vision-language models (VLMs) into robotics applications. Recent works focus on using LLMs and VLMs to improve the performance of robotics tasks, such as manipulation and navigation. Despite these improvements, analyzing the safety of such systems remains underexplored yet extremely critical. LLMs and VLMs are highly susceptible to adversarial inputs, prompting a significant inquiry into the safety of robotic systems. This concern is important because robotics operate in the physical world where erroneous actions can result in severe consequences. This paper explores this issue thoroughly, presenting a mathematical formulation of potential attacks on LLM/VLM-based robotic systems and offering experimental evidence of the safety challenges. Our empirical findings highlight a significant vulnerability: simple modifications to the input can drastically reduce system effectiveness. Specifically, our results demonstrate an average performance deterioration of 19.4% under minor input prompt modifications and a more alarming 29.1% under slight perceptual changes. These findings underscore the urgent need for robust countermeasures to ensure the safe and reliable deployment of advanced LLM/VLM-based robotic systems.


title: Morphology and Mach Number Dependence of Subsonic Bondi-Hoyle Accretion

abstract: We carry out three-dimensional computations of the accretion rate onto an object (of size $R_{\rm sink}$ and mass $m$) as it moves through a uniform medium at a subsonic speed $v_{\infty}$. The object is treated as a fully-absorbing boundary (e.g. a black hole). In contrast to early conjectures, we show that when $R_{\rm sink}\ll R_{A}=2Gm/v^2$ the accretion rate is independent of $v_{\infty}$ and only depends on the entropy of the ambient medium, its adiabatic index, and $m$. Our numerical simulations are conducted using two different numerical schemes via the Athena++ and Arepo hydrodynamics solvers, which reach nearly identical steady-state solutions. We find that pressure gradients generated by the isentropic compression of the flow near the accretor are sufficient to suspend much of the surrounding gas in a near-hydrostatic equilibrium, just as predicted from the spherical Bondi-Hoyle calculation. Indeed, the accretion rates for steady flow match the Bondi-Hoyle rate, and are indicative of isentropic flow for subsonic motion where no shocks occur. We also find that the accretion drag may be predicted using the Safronov number, $\Theta=R_{A}/R_{\rm sink}$, and is much less than the dynamical friction for sufficiently small accretors ($R_{\rm sink}\ll R_{A}$).


title: Exploration-Driven Policy Optimization in RLHF: Theoretical Insights on Efficient Data Utilization

abstract: Reinforcement Learning from Human Feedback (RLHF) has achieved impressive empirical successes while relying on a small amount of human feedback. However, there is limited theoretical justification for this phenomenon. Additionally, most recent studies focus on value-based algorithms despite the recent empirical successes of policy-based algorithms. In this work, we consider an RLHF algorithm based on policy optimization (PO-RLHF). The algorithm is based on the popular Policy Cover-Policy Gradient (PC-PG) algorithm, which assumes knowledge of the reward function. In PO-RLHF, knowledge of the reward function is not assumed and the algorithm relies on trajectory-based comparison feedback to infer the reward function. We provide performance bounds for PO-RLHF with low query complexity, which provides insight into why a small amount of human feedback may be sufficient to get good performance with RLHF. A key novelty is our trajectory-level elliptical potential analysis technique used to infer reward function parameters when comparison queries rather than reward observations are used. We provide and analyze algorithms in two settings: linear and neural function approximation, PG-RLHF and NN-PG-RLHF, respectively.


title: Non-adaptive Bellman-Ford: Yen's improvement is optimal

abstract: The Bellman-Ford algorithm for single-source shortest paths repeatedly updates tentative distances in an operation called relaxing an edge. In several important applications a non-adaptive (oblivious) implementation is preferred, which means fixing the entire sequence of relaxations upfront, independent of the edge-weights. In a dense graph on $n$ vertices, the algorithm in its standard form performs $(1 + o(1))n^3$ relaxations. An improvement by Yen from 1970 reduces the number of relaxations by a factor of two. We show that no further constant-factor improvements are possible, and every non-adaptive deterministic algorithm based on relaxations must perform $(\frac{1}{2} - o(1))n^3$ steps. This improves an earlier lower bound of Eppstein of $(\frac{1}{6} - o(1))n^3$. Given that a non-adaptive randomized variant of Bellman-Ford with at most $(\frac{1}{3} + o(1))n^3$ relaxations (with high probability) is known, our result implies a strict separation between deterministic and randomized strategies, answering an open question of Eppstein.


title: Evaluating NeRFs for 3D Plant Geometry Reconstruction in Field Conditions

abstract: We evaluate different Neural Radiance Fields (NeRFs) techniques for reconstructing (3D) plants in varied environments, from indoor settings to outdoor fields. Traditional techniques often struggle to capture the complex details of plants, which is crucial for botanical and agricultural understanding. We evaluate three scenarios with increasing complexity and compare the results with the point cloud obtained using LiDAR as ground truth data. In the most realistic field scenario, the NeRF models achieve a 74.65% F1 score with 30 minutes of training on the GPU, highlighting the efficiency and accuracy of NeRFs in challenging environments. These findings not only demonstrate the potential of NeRF in detailed and realistic 3D plant modeling but also suggest practical approaches for enhancing the speed and efficiency of the 3D reconstruction process.


title: Exploiting spacetime symmetry in dissipative nonlinear multimode amplifiers for output control

abstract: Time-reversal symmetry enables shaping input waves to control output waves in many linear and nonlinear systems; however energy dissipation violates such symmetry. We consider a saturated multimode fiber amplifier in which light generates heat flow and suffers nonlinear thermo-optical scattering, breaking time-reversal symmetry. We identify a spacetime symmetry which maps the target output back to an input field. This mapping employs phase conjugation, gain and absorption substitution but not time reversal, and holds in steady-state and for slowly varying inputs. Our results open the possibility of output control of a saturated multimode fiber amplifier.


title: Wavelet compressed, modified Hilbert transform in the space-time discretization of the heat equation

abstract: On a finite time interval $(0,T)$, we consider the multiresolution Galerkin discretization of a modified Hilbert transform $\mathcal H_T$ which arises in the space-time Galerkin discretization of the linear diffusion equation. To this end, we design spline-wavelet systems in $(0,T)$ consisting of piecewise polynomials of degree $\geq 1$ with sufficiently many vanishing moments which constitute Riesz bases in the Sobolev spaces $ H^{s}_{0,}(0,T)$ and $ H^{s}_{,0}(0,T)$. These bases provide multilevel splittings of the temporal discretization spaces into "increment" or "detail" spaces of direct sum type. Via algebraic tensor-products of these temporal multilevel discretizations with standard, hierarchic finite element spaces in the spatial domain (with standard Lagrangian FE bases), sparse space-time tensor-product spaces are obtained, which afford a substantial reduction in the number of the degrees of freedom as compared to time-marching discretizations. In addition, temporal spline-wavelet bases allow to compress certain nonlocal integrodifferential operators which appear in stable space-time variational formulations of initial-boundary value problems, such as the heat equation and the acoustic wave equation. An efficient preconditioner is proposed that affords linear complexity solves of the linear system of equations which results from the sparse space-time Galerkin discretization.


title: Classifying topology in photonic crystal slabs with radiative environments

abstract: In the recent years, photonic Chern materials have attracted substantial interest as they feature topological edge states that are robust against disorder, promising to realize defect-agnostic integrated photonic crystal slab devices. However, the out-of-plane radiative losses in those photonic Chern slabs has been previously neglected, yielding limited accuracy for predictions of these systems' topological protection. Here, we develop a general framework for measuring the topological protection in photonic systems, such as in photonic crystal slabs, while accounting for in-plane and out-of-plane radiative losses. Our approach relies on the spectral localizer that combines the position and Hamiltonian matrices of the system to draw a real-picture of the system's topology. This operator-based approach to topology allows us to use an effective Hamiltonian directly derived from the full-wave Maxwell equations after discretization via finite-elements method (FEM), resulting in the full account of all the system's physical processes. As the spectral FEM-localizer is constructed solely from FEM discretization of the system's master equation, the proposed framework is applicable to any physical system and is compatible with commonly used FEM software. Moving forward, we anticipate the generality of the method to aid in the topological classification of a broad range of complex physical systems.

categories: physics.optics cond-mat.mes-hall math-ph math.MP


title: A 3D phase-field based Eulerian variational framework for multiphase fluid-structure interaction with contact dynamics

abstract: Using a fixed Eulerian mesh, the phase-field method has been successfully utilized for a broad range of moving boundary problems involving multiphase fluids and single-phase fluid-structure interaction. Nevertheless, multiphase fluids interacting with multiple solids are rarely explored, especially for large-scale finite element simulations with contact dynamics. In this work, we introduce a novel parallelized three-dimensional fully Eulerian variational framework for simulating multiphase fluids interacting with multiple deformable solids subjected to contact dynamics. In the framework, each solid or fluid phase is identified by a standalone phase indicator. Moreover the phase indicators are initialized by the grid cell method, which restricts the calculation to several grid cells. A diffuse interface description is employed for a smooth interpolation of the physical properties across the phases, yielding unified mass and momentum conservation equations for the coupled dynamical interactions. For each solid object, temporal integration is carried out to track the strain evolution in an Eulerian frame of reference. The coupled differential equations are solved in a partitioned iterative manner. We first verify the framework against reference numerical data in a two-dimensional case of a rotational disk in a lid-driven cavity flow. The case is generalized to a rotational sphere in a lid-driven cavity flow to showcase large deformation and rotational motion of solids and examine the convergence in three dimensions. We then simulate the falling of an immersed solid sphere on an elastic block under gravitational force to demonstrate the translational motion and the solid-to-solid contact in a fluid environment. Finally, we demonstrate the framework for a ship-ice interaction problem involving multiphase fluids with an air-water interface and contact between a floating ship and ice floes.