Search

It is well-known that the conditional mutual information of a quantum state is zero if, and only if, the quantum state is a quantum Markov chain. Replacing the Umegaki relative entropy in the definition of the conditional mutual information by the Belavkin-Staszewski (BS) relative entropy, we obtain the BS-conditional mutual information, and we call the states with zero BS-conditional mutual information Belavkin-Staszewski quantum Markov chains. In this article, we establish a correspondence which relates quantum Markov chains and BS-quantum Markov chains. This correspondence allows us to find a recovery map for the BS-relative entropy in the spirit of the Petz recovery map. Moreover, we prove a structural decomposition of the Belavkin-Staszewski quantum Markov chains and also study states for which the BS-conditional mutual information is only approximately zero. As an application, we show that the conditional mutual information of the state associated to the Gibbs state of a quantum spin chain with local, finite-range, translation-invariant interactions at any positive temperature via this new correspondence decays superexponentially fast with the size of the middle system., Comment: 19 pages, 1 figure

A growing number of directly-imaged companions have been recently characterised, with robust constraints on carbon-to-oxygen ratios and even isotopic ratios. Many companions and isolated targets have also shown spectral variability. In this work we observed the super-Jupiter AB~Pictoris~b across four consecutive nights using VLT/CRIRES+ as part of the ESO SupJup survey, exploring how the constraints on chemical composition and temperature profile change over time using spectral line shape variations between nights. We performed atmospheric retrievals of the high-resolution observations and found broadly consistent results across all four nights, but there were differences for some parameters. We clearly detect H$_2$O, $^{12}$CO and $^{13}$CO in each night, but abundances varied by $\sim2\sigma$, which was correlated to the deep atmosphere temperature profiles. We also found differences in the $^{12}$C$/^{13}$C ratios in each night by up to $\sim3\sigma$, which seemed to be correlated with the cloud deck pressure. Our combined retrieval simultaneously analysing all nights together constrained broadly the average of each night individually, with the C/O$=0.59\pm0.01$, consistent with solar composition, and $^{12}$C$/^{13}$C~$ = 102\pm8$, slightly higher than the ISM and Solar System values. We also find a low projected rotational velocity, suggesting that AB~Pictoris~b is either intrinsically a slow rotator due to its young age or that the spin axis is observed pole-on with a $\sim90^\circ$ misalignment with its orbit inclination. Future observations will be able to further explore the variability and orbit of AB~Pictoris~b as well as for other companions., Comment: 22 pages, 17 figures, 4 tables, accepted for publication in MNRAS

Context. The supernova remnant (SNR) W44 and its surroundings are a prime target for studying the acceleration of cosmic rays (CRs). Several previous studies established an extended gamma-ray emission that is set apart from the radio shell of W44. This emission is thought to originate from escaped high-energy CRs that interact with a surrounding dense molecular cloud complex. Aims. We present a detailed analysis of Fermi-LAT data with an emphasis on the spatial and spectral properties of W44 and its surroundings. We also report the results of the observations performed with the MAGIC telescopes of the northwestern region of W44. Finally, we present an interpretation model to explain the gamma-ray emission of the SNR and its surroundings. Methods. We first performed a detailed spatial analysis of 12 years of Fermi-LAT data at energies above 1 GeV, in order to exploit the better angular resolution, while we set a threshold of 100MeV for the spectral analysis. We performed a likelihood analysis of 174 hours of MAGIC data above 130 GeV using the spatial information obtained with Fermi-LAT. Results. The combined spectra of Fermi-LAT and MAGIC, extending from 100MeV to several TeV, were used to derive constraints on the escape of CRs. Using a time-dependent model to describe the particle acceleration and escape from the SNR, we show that the maximum energy of the accelerated particles has to be ' 40 GeV. However, our gamma-ray data suggest that a small number of lower-energy particles also needs to escape. We propose a novel model, the broken-shock scenario, to account for this effect and explain the gamma-ray emission.

Mrk 421 was in its most active state around early 2010, which led to the highest TeV gamma-ray flux ever recorded from any active galactic nuclei. We aim to characterize the multiwavelength behavior during this exceptional year for Mrk 421, and evaluate whether it is consistent with the picture derived with data from other less exceptional years. We investigated the period from November 5, 2009, (MJD 55140) until July 3, 2010, (MJD 55380) with extensive coverage from very-high-energy (VHE; E$\,>\,$100$\,$GeV) gamma rays to radio with MAGIC, VERITAS, Fermi-LAT, RXTE, Swift, GASP-WEBT, VLBA, and a variety of additional optical and radio telescopes. We investigated the variability and correlation behavior among different energy bands in great detail. We find the strongest variability in X-rays and VHE gamma rays, and PSDs compatible with power-law functions. We observe strong correlations between X-rays and VHE gamma rays. We also report a marginally significant positive correlation between high-energy (HE; E$\,>\,$100$\,$MeV) gamma rays and the ultraviolet band. We detected marginally significant correlations between the HE and VHE gamma rays, and between HE gamma rays and the X-ray, that disappear when the large flare in February 2010 is excluded from the correlation study. The activity of Mrk 421 also yielded the first ejection of features in the VLBA images of the jet of Mrk 421. Yet the large uncertainties in the ejection times of these radio features prevent us from firmly associating them to the specific flares recorded during the campaign. We also show that the collected multi-instrument data are consistent with a scenario where the emission is dominated by two regions, a compact and extended zone, which could be considered as a simplified implementation of an energy-stratified jet as suggested by recent IXPE observations., Comment: Accepted for publication in Astronomy & Astrophysics. Corresponding authors: Felix Schmuckermaier, David Paneque, Axel Arbet Engels

The trimmed mean of $n$ scalar random variables from a distribution $P$ is the variant of the standard sample mean where the $k$ smallest and $k$ largest values in the sample are discarded for some parameter $k$. In this paper, we look at the finite-sample properties of the trimmed mean as an estimator for the mean of $P$. Assuming finite variance, we prove that the trimmed mean is ``sub-Gaussian'' in the sense of achieving Gaussian-type concentration around the mean. Under slightly stronger assumptions, we show the left and right tails of the trimmed mean satisfy a strong ratio-type approximation by the corresponding Gaussian tail, even for very small probabilities of the order $e^{-n^c}$ for some $c>0$. In the more challenging setting of weaker moment assumptions and adversarial sample contamination, we prove that the trimmed mean is minimax-optimal up to constants., Comment: 37 pages

Strong gravitational lenses come in many forms, but are typically divided into two populations: galaxies, and groups and clusters of galaxies. When calculating the properties of the images we expect to see from these lenses, it is typically assumed that each lens is roughly a singular isothermal sphere. In reality, the largest objects in the Universe (i.e. galaxy clusters) are highly irregular and composed of many components due to a history of (or active) hierarchical mergers. In this work, we analyze the discrepancies in the observables of strongly lensed transients in both scenarios, namely relative magnifications, time delays, and image multiplicities. Focusing on gravitational waves, we compare the detection rates between the single spherical dark matter halo models found in the literature, and publicly available state-of-the-art cluster lens models. We find there to be approximately an order of magnitude fewer detection of strongly lensed transients in the realistic model case, likely caused by their loss of overall strong lensing optical depth. We also report detection rates in the weak lensing or single-image regime. Additionally, we find a systemic shift towards lower time delays between the brightest image pairs in the cases of the realistic models, as well as higher fractions of positive versus negative parity images, as seen elsewhere in the literature. This significant deviation in the joint relative magnification factor-time delay distribution will hinder the feasibility of the reconstruction of lenses through time domain transients alone, but can still provide a lower limit on the lens mass., Comment: 13 pages, 10 figures

Jetted Active Galactic Nuclei (AGN) exhibit variability across a wide range of time scales. Traditionally, this variability can often be modeled well as a stochastic process. However, in certain cases, jetted AGN variability displays regular patterns, enabling us to conduct investigations aimed at understanding its origins. Additionally, a novel type of variability has emerged in jetted AGN lightcurves, specifically, the observation of a long-term trend characterized by a linear increase of the flux with time in blazars such as PG 1553+113, which is among the objects most likely to display periodic behavior. In this paper, we present the results of a systematic search for long-term trends, spanning $\approx$10\, years, utilizing 12 years of Fermi-LAT observations. The study is focused on detecting the presence of linear or quadratic long-term trends in a sample of 3308 jetted AGN. Our analysis has identified 40 jetted AGN that exhibit long-term trends, each with distinct properties, which we also characterize in this study. These long-term trends may originate from the dynamics of a supermassive black hole binary system, or they could be the result of intrinsic phenomena within the jet itself. Our findings can help in addressing questions pertaining to the astrophysical origins of variability and periodicity within jetted AGN., Comment: 22 pages, 4 tables, 6 fgures

It is imperative to develop an intrusion prevention system (IPS), specifically designed for autonomous robotic systems. This is due to the unique nature of these cyber-physical systems (CPS), which are not merely typical distributed systems. These systems employ their own systems software (i.e. robotic middleware and frameworks) and execute distinct components to facilitate interaction with various sensors and actuators, and other robotic components (e.g. cognitive subsystems). Furthermore, as cyber-physical systems, they engage in interactions with humans and their physical environment, as exemplified by social robots. These interactions can potentially lead to serious consequences, including physical damage. In response to this need, we have designed and implemented RIPS, an intrusion prevention system tailored for robotic applications based on ROS 2, the framework that has established itself as the de facto standard for developing robotic applications. This manuscript provides a comprehensive exposition of the issue, the security aspects of ROS 2 applications, and the key points of the threat model we created for our robotic environment. It also describes the architecture and the implementation of our initial research prototype and a language specifically designed for defining detection and prevention rules for diverse, real-world robotic scenarios. Moreover, the manuscript provides a comprehensive evaluation of the approach, that includes a set of experiments with a real social robot executing a well known testbed used in international robotic competitions., Comment: 45 pages, 5 figures

The BL Lacertae object VER J0521+211 underwent a notable flaring episode in February 2020. A short-term monitoring campaign, led by the MAGIC (Major Atmospheric Gamma Imaging Cherenkov) collaboration, covering a wide energy range from radio to very-high-energy (VHE, 100 GeV < E < 100 TeV) gamma rays was organised to study its evolution. These observations resulted in a consistent detection of the source over six consecutive nights in the VHE gamma-ray domain. Combining these nightly observations with an extensive set of multiwavelength data made modelling of the blazar's spectral energy distribution (SED) possible during the flare. This modelling was performed with a focus on two plausible emission mechanisms: i) a leptonic two-zone synchrotron-self-Compton scenario, and ii) a lepto-hadronic one-zone scenario. Both models effectively replicated the observed SED from radio to the VHE gamma-ray band. Furthermore, by introducing a set of evolving parameters, both models were successful in reproducing the evolution of the fluxes measured in different bands throughout the observing campaign. Notably, the lepto-hadronic model predicts enhanced photon and neutrino fluxes at ultra-high energies (E > 100 TeV). While the photon component, generated via decay of neutral pions, is not directly observable as it is subject to intense pair production (and therefore extinction) through interactions with the cosmic microwave background photons, neutrino detectors (e.g. IceCube) can probe the predicted neutrino component. Finally, the analysis of the gamma-ray spectra, as observed by MAGIC and the Fermi-LAT telescopes, yielded a conservative 95\% confidence upper limit of z \leq 0.244 for the redshift of this blazar., Comment: Accepted for publication in A&A

For gravitationally lensed type II signals, the phase of the dominant (2, 2) mode and the higher order (3, 3) mode is offset by $-\pi/12$. Using this, we develop a test for type II imagery by allowing the phases of the (2,2) and (3,3) modes to vary separately and introducing a new waveform parameter to represent the phase offset between the two. We use simulated, low mass ratio, precessing signals to show that the test can reproduce the $-\pi/12$ phase offset when detected by three detectors for H-L optimal SNR $\gtrsim$ 40 and $\mathcal{M} \leq 30$. We analyze GW190412 and GW190814 using this parameterization, measuring the offset to be $0.13^{+0.22}_{-0.17}$ for GW190412 and $-0.05^{+0.20}_{-0.22}$ for GW190814. We also measure the Bayes Factor in support of no phase offset, $\log_{10} \mathcal{B}_{\Delta \varphi = 0}$, to be $-0.14$ for GW190412 and $0.21$ for GW190814., Comment: 10 pages, 7 figures

Observations of strongly gravitationally lensed gravitational wave (GW) sources provide a unique opportunity for constraining their transverse motion, which otherwise is exceedingly hard for GW mergers in general. Strong lensing makes this possible when two or more images of the lensed GW source are observed, as each image essentially allows the observer to see the GW source from different directional lines-of-sight. If the GW source is moving relative to the lens and observer, the observed GW signal from one image will therefore generally appear blue- or redshifted compared to GW signal from the other image. This velocity induced differential Doppler shift gives rise to an observable GW phase shift between the GW signals from the different images, which provides a rare glimpse into the relative motion of GW sources and their host environment across redshift. We illustrate that detecting such GW phase shifts is within reach of next-generation ground-based detectors such as Einstein Telescope, that is expected to detect $\sim$hundreds of lensed GW mergers per year. This opens up completely new ways of inferring the environment of GW sources, as well as studying cosmological velocity flows across redshift., Comment: comments welcome

Although a good portion of elementary linear algebra concerns itself with matrices over a field such as $\mathbb{R}$ or $\mathbb{C}$, many combinatorial problems naturally surface when we instead work with matrices over a finite field. As some recent work has been done in these areas, we turn our attention to the problem of enumerating the square matrices with entries in $\mathbb{Z}_{p^k}$ that are diagonalizable over $\mathbb{Z}_{p^k}$. This turns out to be significantly more nontrivial than its finite field counterpart due to the presence of zero divisors in $\mathbb{Z}_{p^k}$., Comment: 17 pages, 4 figures, 3 tables, published in Involve, a Journal of Mathematics

A majority of blazars exhibit variable emission across the entire electromagnetic spectrum, observed over various time scales. In particular, discernible periodic patterns are detected in the {\gamma}-ray light curves of a few blazars, such as PG 1553+113, S5 1044+71, and PKS 0426-380. The presence of trends, flares, and noise complicates the detection of periodicity, requiring careful analysis to determine whether these patterns are related to emission mechanisms within the source or occur by chance. We employ Singular Spectrum Analysis (SSA) for the first time on data from the Large Area Telescope (LAT) aboard the Fermi Gamma-ray Space Telescope to systematically search for periodicity in the time domain, using 28-day binned light curves. Our aim is to isolate any potential periodic nature of the emission from trends and noise, thereby reducing uncertainties in revealing periodicity. Additionally, we aim to characterize long-term trends and develop a forecasting algorithm based on SSA, enabling accurate predictions of future emission behavior. We apply SSA to analyze 494 sources detected by Fermi-LAT, focusing on identifying and isolating periodic components from trends and noise in their {\gamma}-ray light curves. We calculate the Lomb-Scargle Periodogram for the periodic components extracted by SSA to determine the most significant periods. The local and global significance of these periods is then assessed to validate their authenticity. Our analysis identifies 46 blazars as potential candidates for quasi-periodic {\gamma}-ray emissions, each with a local significance level >= 2{\sigma}. Notably, 33 of these candidates exhibit a local significance of >= 4{\sigma} (corresponding to a global significance of >= 2.2{\sigma}). Our findings introduce 25 new {\gamma}-ray candidates, effectively doubling the number of potentially periodic sources., Comment: 23 pages, 8 figures, 4 tables; Submitted to A&A

This paper examines the role of 5G communication in the Mobilities for EU project, a collaborative initiative involving 29 partners and 11 pilots aimed at revolutionizing urban mobility through electrification, automation, and connectivity. Focusing on Dresden as a Lead City, we explore the integration of 27 innovative solutions, including autonomous freight transport, eBuses, and charging robots, using a 5G communication network as the central framework. We analyze how 5G enables seamless connectivity and real-time data processing across diverse technologies, fostering interdependencies and synergies. This approach not only provides a cohesive understanding of the project's scope but also demonstrates 5G's critical role in smart city infrastructure. We evaluate the anticipated impact on sustainability metrics such as air quality, noise levels, CO2 emissions, and traffic congestion. The paper concludes by discussing challenges and strategies in leveraging 5G for comprehensive urban mobility solutions and its potential impact on future smart city developments., Comment: 6 pages

Magnetron sputtering is an essential technique in combinatorial materials science, enabling the efficient synthesis of thin-film materials libraries with continuous compositional gradients. For exploring multidimensional search spaces, minimizing preliminary experiments is essen-tial, as numerous materials libraries are required to adequately cover the space, making it crucial to fabricate only those libraries that are absolutely necessary. This can be achieved by Monte Carlo particle simulations to model the deposition profile, e.g. by SIMTRA, which is an established package mainly designed for single cathode simulations. A strong enhance-ment of its capabilities is the development of a Python-based wrapper, designed to simulate multi-cathode sputter processes through parallel Monte Carlo simulations. By modeling a sputter chamber and determining the relationship between deposition power and rate for an exemplary quaternary system Ni-Pd-Pt-Ru, we achieve a match between simulated and measured compositions, with a mean Euclidean distance of 3.5%. The object-oriented design of the package allows easy customization and enables the definition of complex sputter sys-tems. Due to parallelization, simulating multiple cathodes results in no additional simulation time. These additions extend the capabilities of SIMTRA making it applicable in combinatorial materials research.

A core-corona decomposition of compact (neutron) star models is compared to recent NICER data of masses and radii. It is in particular interesting to capture the outlier XTE~J1814-338. Instead of integrating the TOV equations from the center to surface, we follow here another pathway by accommodating all uncertainties of the equation(s) of state (EoS) at supra-nuclear density or/and an unknown dark matter admixture in a parameterization of the core by its radius $r_x$, the included mass $m_x$ and the pressure $p_x$ at $r_x$. The corona, which may be dubbed also envelope or halo or outer crust, is assumed to be of standard-model matter where the EoS is supposed to be faithfully known., Comment: arXiv admin note: substantial text overlap with arXiv:2302.01389

We study a flaring activity of the HSP Mrk421 that was characterized from radio to very-high-energy (VHE; E $>0.1$TeV) gamma rays with MAGIC, Fermi-LAT, Swift, XMM-Newton and several optical and radio telescopes. These observations included, for the first time for a gamma-ray flare of a blazar, simultaneous X-ray polarization measurements with IXPE. We find substantial variability in both X-rays and VHE gamma rays throughout the campaign, with the highest VHE flux above 0.2 TeV occurring during the IXPE observing window, and exceeding twice the flux of the Crab Nebula. However, the VHE and X-ray spectra are on average softer, and the correlation between these two bands weaker that those reported in previous flares of Mrk421. IXPE reveals an X-ray polarization degree significantly higher than that at radio and optical frequencies. The X-ray polarization angle varies by $\sim$100$^\circ$ on timescales of days, and the polarization degree changes by more than a factor 4. The highest X-ray polarization degree reaches 26%, around which a X-ray counter-clockwise hysteresis loop is measured with XMM-Newton. It suggests that the X-ray emission comes from particles close to the high-energy cutoff, hence possibly probing an extreme case of the Turbulent Extreme Multi-Zone model. We model the broadband emission with a simplified stratified jet model throughout the flare. The polarization measurements imply an electron distribution in the X-ray emitting region with a very high minimum Lorentz factor, which is expected in electron-ion plasma, as well as a variation of the emitting region size up to a factor of three during the flaring activity. We find no correlation between the fluxes and the evolution of the model parameters, which indicates a stochastic nature of the underlying physical mechanism. Such behaviour would be expected in a highly turbulent electron-ion plasma crossing a shock front., Comment: Submitted to Astronomy and Astrophysics. Corresponding authors: Axel Arbet-Engels, Lea Heckmann, David Paneque

OT 081 is a well-known, luminous blazar that is remarkably variable in many energy bands. We present the first broadband study of the source which includes very-high-energy (VHE, $E>$100\,GeV) $\gamma$-ray data taken by the MAGIC and H.E.S.S. imaging Cherenkov telescopes. The discovery of VHE $\gamma$-ray emission happened during a high state of $\gamma$-ray activity in July 2016, observed by many instruments from radio to VHE $\gamma$-rays. We identify four states of activity of the source, one of which includes VHE $\gamma$-ray emission. Variability in the VHE domain is found on daily timescales. The intrinsic VHE spectrum can be described by a power-law with index $3.27\pm0.44_{\rm stat}\pm0.15_{\rm sys}$ (MAGIC) and $3.39\pm0.58_{\rm stat}\pm0.64_{\rm sys}$ (H.E.S.S.) in the energy range of 55--300\,GeV and 120--500\,GeV, respectively. The broadband emission cannot be sucessfully reproduced by a simple one-zone synchrotron self-Compton model. Instead, an additional external Compton component is required. We test a lepto-hadronic model that reproduces the dataset well and a proton-synchrotron dominated model that requires an extreme proton luminosity. Emission models that are able to successfully represent the data place the emitting region well outside of the Broad Line Region (BLR) to a location at which the radiative environment is dominated by the infrared thermal radiation field of the dusty torus. In the scenario described by this flaring activity, the source appears to be an FSRQ, in contrast with past categorizations. This suggests that the source can be considered to be a transitional blazar, intermediate between BL~Lac and FSRQ objects., Comment: Accepted on MNRAS Corresponding authors: M. Manganaro, J. Becerra Gonz\'alez, M. Seglar-Arroyo, D. A. Sanchez

The integration of non-terrestrial networks (NTN) into 5G new radio (NR) has opened up the possibility of developing a new positioning infrastructure using NR signals from Low-Earth Orbit (LEO) satellites. LEO-based cellular positioning offers several advantages, such as a superior link budget, higher operating bandwidth, and large forthcoming constellations. Due to these factors, LEO-based positioning, navigation, and timing (PNT) is a potential enhancement for NTN in 6G cellular networks. However, extending the existing terrestrial cellular positioning methods to LEO-based NTN positioning requires considering key fundamental enhancements. These include creating broad positioning beams orthogonal to conventional communication beams, time-domain processing at the user equipment (UE) to resolve large delay and Doppler uncertainties, and efficiently accommodating positioning reference signals (PRS) from multiple satellites within the communication resource grid. In this paper, we present the first set of design insights by incorporating these enhancements and thoroughly evaluating LEO-based positioning, considering the constraints and capabilities of the NR-NTN physical layer. To evaluate the performance of LEO-based NTN positioning, we develop a comprehensive NR-compliant simulation framework, including LEO orbit simulation, transmission (Tx) and receiver (Rx) architectures, and a positioning engine incorporating the necessary enhancements. Our findings suggest that LEO-based NTN positioning could serve as a complementary infrastructure to existing Global Navigation Satellite Systems (GNSS) and, with appropriate enhancements, may also offer a viable alternative., Comment: 7 pages, 6 figures, submitted to IEEE Communications Magazine

We present the experimental apparatus enabling the observation of the heteronuclear Efimov effect in an optically trapped ultracold mixture of $^6$Li-$^{133}$Cs with high-resolution control of the interactions. A compact double-species Zeeman slower consisting of four interleaving helical coils allows for a fast-switching between two optimized configurations for either Li or Cs and provides an efficient sequential loading into their respective MOTs. By means of a bichromatic optical trapping scheme based on species-selective trapping we prepare mixtures down to 100 nK of 1$\times$ 10$^4$ Cs atoms and 7$\times$ 10$^3$ Li atoms. Highly stable magnetic fields allow high-resolution atom-loss spectroscopy and enable to resolve splitting in the loss feature of a few tens of milligauss. These features allowed for a detailed study of the Efimov effect.

Recent efforts in the development of autonomous driving technology have induced great advancements in perception, planning and control systems. Model predictive control is one of the most popular advanced control methods, but its application to nonlinear systems still depends on the development of computationally efficient methods. This work presents a nonlinear model predictive control formulation based on real-time iteration using an implicit discretization of the system's dynamics, with the objective of achieving greater prediction accuracy and lower computational cost when dealing with stiff dynamical systems, as is the case for vehicle dynamics. The proposed method is described and later evaluated on a simulation scenario considering modeling errors and external disturbances. The presented results demonstrate the effectiveness of the method when it comes to tracking a given trajectory and its low computational burden, measured in terms of execution time., Comment: This works was submitted, accepted and presented at the 2024 Simp\'osio Brasileiro de Automa\c{c}\~ao Inteligente - SBAI

The notion of defects in crystalline phases of matter has been extremely powerful for understanding crystal growth, deformation and melting. These discontinuities in the periodic order of crystals are mathematically described by the Burgers vector, derived from the particle displacements, which encapsulates the direction and magnitude of slip relative to the undeformed state. Since the reference structure of the crystal is known a priori, the Burgers vector can be determined experimentally using both imaging and diffraction methods to measure the lattice distortion, and thus infer the particle displacements. Glasses have structures that lack the periodicity of crystals, and thus a well-defined reference state. Yet, measurable structural parameters can still be obtained from diffraction from a glass. Here we examine the usefulness of these parameters to probe deformation in glasses. We find that coordinated transformations in the centrosymmetry of local particle arrangements are a strong marker of plastic events. Moreover, we investigate two geometric indicators that can be derived from distortions in local diffraction patterns, namely the continuous Burgers vector and the quadrupolar strain. We find that the Burgers vector again emerges as a robust and sensitive metric for understanding local structural transformations due to mechanical deformation, even in structural glasses., Comment: 13 pages 5 figures

The Duke Robotics Club is proud to present our robot for the 2023 RoboSub Competition: Oogway. Oogway marks one of the largest design overhauls in club history. Beyond a revamped formfactor, some of Oogway's notable features include all-new computer vision software, advanced sonar integration, novel acoustics hardware processing, and upgraded stereoscopic cameras. Oogway was built on the principle of independent, well-integrated, and reliable subsystems. Individual components and subsystems were tested and designed separately. Oogway's most advanced capabilities are a result of the tight integration between these subsystems. Such examples include sonar-assisted computer vision algorithms and robot-agnostic controls configured in part through the robot's 3D model. The success of constructing and testing Oogway in under 2 year's time can be attributed to 20+ contributing club members, supporters within Duke's Pratt School of Engineering, and outside sponsors., Comment: arXiv admin note: text overlap with arXiv:2410.09684

The Duke Robotics Club is proud to present our robot for the 2024 RoboSub Competition: Oogway. Now in its second year, Oogway has been dramatically upgraded in both its capabilities and reliability. Oogway was built on the principle of independent, well-integrated, and reliable subsystems. Individual components and subsystems were tested and designed separately. Oogway's most advanced capabilities are a result of the tight integration between these subsystems. Such examples include a re-envisioned controls system, an entirely new electrical stack, advanced sonar integration, additional cameras and system monitoring, a new marker dropper, and a watertight capsule mechanism. These additions enabled Oogway to prequalify for Robosub 2024.

In the presence of ring exchange interactions, bosons in a ladder-like lattice may form the bosonic analogon of a correlated metal, known as the d-wave Bose liquid (DBL). In this Letter, we show that a chain of trapped ions with three internal levels can be used to mimic a ladder-like system constrained to a maximum occupation of one boson per rung. The setup enables tunable ring exchange interactions, giving rise to a transition between a polarized regime with all bosons confined to one leg and the DBL regime. The latter state is characterized by a splitting of the peak in the momentum distribution and an oscillating pair correlation function., Comment: 7 pages, 6 figures

Instruments for gamma-ray astronomy at Very High Energies ($E>100\,{\rm GeV}$) have traditionally derived their scientific results through proprietary data and software. Data standardisation has become a prominent issue in this field both as a requirement for the dissemination of data from the next generation of gamma-ray observatories and as an effective solution to realise public data legacies of current-generation instruments. Specifications for a standardised gamma-ray data format have been proposed as a community effort and have already been successfully adopted by several instruments. We present the first production of standardised data from the Major Atmospheric Gamma-ray Imaging Cherenkov (MAGIC) telescopes. We converted $166\,{\rm h}$ of observations from different sources and validated their analysis with the open-source software Gammapy. We consider six data sets representing different scientific and technical analysis cases and compare the results obtained analysing the standardised data with open-source software against those produced with the MAGIC proprietary data and software. Aiming at a systematic production of MAGIC data in this standardised format, we also present the implementation of a database-driven pipeline automatically performing the MAGIC data reduction from the calibrated down to the standardised data level. In all the cases selected for the validation, we obtain results compatible with the MAGIC proprietary software, both for the manual and for the automatic data productions. Part of the validation data set is also made publicly available, thus representing the first large public release of MAGIC data. This effort and this first data release represent a technical milestone toward the realisation of a public MAGIC data legacy., Comment: Accepted for publication in the Journal of High Energy Astrophysics

In our recent investigation, we utilized a century's worth of archival optical data to search for a decade-long periodicity from the blazar PG 1553+113, finding a hint of a 22-year period. Building on this foundation, the current study extends our analysis to include 10 blazars from the Fermi Large Area Telescope 2FHL catalog to uncover similar long-term periodic behavior. To ensure the reliability of our findings, we consider the impact of observational limitations, such as temporal gaps and uneven sampling, which could potentially introduce artifacts or false periodic signals. Our analysis reveals that 4 of these blazars (AP Librae, MKN 421, MKN 501, PG 1246+586) exhibit decade-long periods in their optical light curves, albeit 3 of them may be influenced by noise. However, a likely genuine period of approximately 51 $\pm$ 9 yr is identified for MKN 421., Comment: 14 pages, 13 figures, Submitted to The Astrophysical Journal

High-resolution spectroscopic characterization of young super-Jovian planets enables precise constraints on elemental and isotopic abundances of their atmospheres. As part of the ESO SupJup Survey, we present high-resolution spectral observations of two wide-orbit super-Jupiters in YSES 1 (or TYC 8998-760-1) using the upgraded VLT/CRIRES+ (R~100,000) in K-band. We carry out free atmospheric retrieval analyses to constrain chemical and isotopic abundances, temperature structures, rotation velocities, and radial velocities. We confirm the previous detection of 13CO in YSES 1 b at a higher significance of 12.6{\sigma}, but point to a higher 12CO/13CO ratio of 88+/-13 (1{\sigma} confidence interval), consistent with the primary's isotope ratio 66+/-5. We retrieve a solar-like composition in YSES 1 b with a C/O=0.57+/-0.01, indicating a formation via gravitational instability or core accretion beyond the CO iceline. Additionally, the observations lead to detections of H2O and CO in the outer planet YSES 1 c at 7.3{\sigma} and 5.7{\sigma}, respectively. We constrain the atmospheric C/O ratio of YSES 1 c to be either solar or subsolar (C/O=0.36+/-0.15), indicating the accretion of oxygen-rich solids. The two companions have distinct vsini, 5.34+/-0.14 km/s for YSES 1 b and 11.3+/-2.1 km/s for YSES 1 c, despite their similar natal environments. This may indicate different spin axis inclinations or effective magnetic braking by the long-lived circumplanetary disk around YSES 1 b. YSES 1 represents an intriguing system for comparative studies of super-Jovian companions and linking present atmospheres to formation histories., Comment: 25 pages, 11 figure, accepted for publication in AJ. The extracted CRIRES+ spectra of the YSES-1 system can be found at https://doi.org/10.5281/zenodo.13664032

Predicting the emergence of multiple chronic conditions (MCC) is crucial for early intervention and personalized healthcare, as MCC significantly impacts patient outcomes and healthcare costs. Graph neural networks (GNNs) are effective methods for modeling complex graph data, such as those found in MCC. However, a significant challenge with GNNs is their reliance on an existing graph structure, which is not readily available for MCC. To address this challenge, we propose a novel generative framework for GNNs that constructs a representative underlying graph structure by utilizing the distribution of the data to enhance predictive analytics for MCC. Our framework employs a graph variational autoencoder (GVAE) to capture the complex relationships in patient data. This allows for a comprehensive understanding of individual health trajectories and facilitates the creation of diverse patient stochastic similarity graphs while preserving the original feature set. These variations of patient stochastic similarity graphs, generated from the GVAE decoder, are then processed by a GNN using a novel Laplacian regularization technique to refine the graph structure over time and improves the prediction accuracy of MCC. A contextual Bandit is designed to evaluate the stochastically generated graphs and identify the best-performing graph for the GNN model iteratively until model convergence. We validate the performance of the proposed contextual Bandit algorithm against $\varepsilon$-Greedy and multi-armed Bandit algorithms on a large cohort (n = 1,592) of patients with MCC. These advancements highlight the potential of the proposed approach to transform predictive healthcare analytics, enabling a more personalized and proactive approach to MCC management.

On the 29th of May 2023, the LIGO-Virgo-KAGRA Collaboration observed a compact binary coalescence event consistent with a neutron star-black hole merger, though the heavier object of mass 2.5-4.5 $M_\odot$ would fall into the purported lower mass gap. An alternative explanation for apparent observations of events in this mass range has been suggested as strongly gravitationally lensed binary neutron stars. In this scenario, magnification would lead to the source appearing closer and heavier than it really is. Here, we investigate the chances and possible consequences for the GW230529 event to be gravitationally lensed. We find this would require high magnifications and we obtain low rates for observing such an event, with a relative fraction of lensed versus unlensed observed events of $2 \times 10^{-3}$ at most. When comparing the lensed and unlensed hypotheses accounting for the latest rates and population model, we find a 1/58 chance of lensing, disfavoring this option. Moreover, when the magnification is assumed to be strong enough to bring the mass of the heavier binary component below the standard limits on neutron star masses, we find high probability for the lighter object to have a sub-solar mass, making the binary even more exotic than a mass-gap neutron star-black hole system. Even when the secondary is not sub-solar, its tidal deformability would likely be measurable, which is not the case for GW230529. Finally, we do not find evidence for extra lensing signatures such as the arrival of additional lensed images, type-II image dephasing, or microlensing. Therefore, we conclude it is unlikely for GW230529 to be a strongly gravitationally lensed binary neutron star signal., Comment: 15 pages, 11 figures

Interfacial ferroelectricity, prevalent in various parallel-stacked layered materials, allows switching of out-of-plane ferroelectric order by in-plane sliding of adjacent layers. Its resilience against doping potentially enables next-generation storage and logic devices. However, studies have been limited to indirect sensing or visualization of ferroelectricity. For transition metal dichalcogenides, there is little knowledge about the influence of ferroelectric order on their intrinsic valley and excitonic properties. Here, we report direct probing of ferroelectricity in few-layer 3R-MoS$_2$ using reflectance contrast spectroscopy. Contrary to a simple electrostatic perception, layer-hybridized excitons with out-of-plane electric dipole moment remain decoupled from ferroelectric ordering, while intralayer excitons with in-plane dipole orientation are sensitive to it. Ab initio calculations identify stacking-specific interlayer hybridization leading to this asymmetric response. Exploiting this sensitivity, we demonstrate optical readout and control of multi-state polarization with hysteretic switching in a field-effect device. Time-resolved Kerr ellipticity reveals a direct correspondence between spin-valley dynamics and stacking order.

We study the regularity of weak solutions to nonlocal in time subdiffusion equations for a wide class of weakly singular kernels appearing in the generalised fractional derivative operator. We prove a weak Harnack inequality for nonnegative weak supersolutions and Holder continuity of weak solutions to such problems. Our results substantially extend the results from our previous work [12] by leaving the framework of distributed order fractional time derivatives and considering a general PC kernel and by also allowing for an inhomogeneity in the PDE from a Lebesgue space of mixed type., Comment: 52 pages

While $U$(1) spin liquids have been extensively studied in both quantum and classical regimes, exact classical $\mathbb{Z}_2$ spin liquids arising from models with nearest-neighbor, bilinear spin interactions are still rare. In this Letter, we explore the four-color Kitaev model as a minimal model for stabilizing classical $\mathbb{Z}_2$ spin liquids across a broad family of tricoordinated lattices. By formulating a $\mathbb{Z}_2$ lattice gauge theory, we identify this spin liquid as being described by an emergent Gauss's law with effective charge-2 condensation, and deconfined fractionalized bond-charge excitations. We complement our findings with Monte Carlo simulations, revealing a crossover from a high-temperature paramagnet to a low-temperature liquid phase characterized by residual entropy, classical $\mathbb{Z}_2$ flux order, and diffuse spin structure factors., Comment: 6 pages, 4 figures

Machine-learned interatomic potentials (MILPs) are rapidly gaining interest for molecular modeling, as they provide a balance between quantum-mechanical level descriptions of atomic interactions and reasonable computational efficiency. However, questions remain regarding the stability of simulations using these potentials, as well as the extent to which the learned potential energy function can be extrapolated safely. Past studies have reported challenges encountered when MILPs are applied to classical benchmark systems. In this work, we show that some of these challenges are related to the characteristics of the training datasets, particularly the inclusion of rigid constraints. We demonstrate that long stability in simulations with MILPs can be achieved by generating unconstrained datasets using unbiased classical simulations if the fast modes are correctly sampled. Additionally, we emphasize that in order to achieve precise energy predictions, it is important to resort to enhanced sampling techniques for dataset generation, and we demonstrate that safe extrapolation of MILPs depends on judicious choices related to the system's underlying free energy landscape and the symmetry features embedded within the machine learning models.

The implementation of realistic quantum devices requires a solid understanding of the nonlocal resources present in quantum channels, and the effects of decoherence on them. Here we quantify nonlocality of bipartite quantum channels and identify its component resisting the effects of dephasing noise. Despite its classical nature, we demonstrate that the latter plays a relevant role in performing quantum protocols, such as state transformations and quantum coding for noisy communication. In the converse direction, we show that simulating certain stochastic processes with quantum channels undergoing decoherence has a communication advantage with respect to their classical simulation.

Schenkerian Analysis (SchA) is a uniquely expressive method of music analysis, combining elements of melody, harmony, counterpoint, and form to describe the hierarchical structure supporting a work of music. However, despite its powerful analytical utility and potential to improve music understanding and generation, SchA has rarely been utilized by the computer music community. This is in large part due to the paucity of available high-quality data in a computer-readable format. With a larger corpus of Schenkerian data, it may be possible to infuse machine learning models with a deeper understanding of musical structure, thus leading to more "human" results. To encourage further research in Schenkerian analysis and its potential benefits for music informatics and generation, this paper presents three main contributions: 1) a new and growing dataset of SchAs, the largest in human- and computer-readable formats to date (>140 excerpts), 2) a novel software for visualization and collection of SchA data, and 3) a novel, flexible representation of SchA as a heterogeneous-edge graph data structure.

Acoustic waves (AW) have recently emerged as an energy-efficient ice removal procedure compatible with functional and industrial-relevant substrates. However, critical aspects at fundamental and experimental levels have yet to be disclosed to optimize their operational conditions. Identifying the processes and mechanisms by which different types of AWs induce de-icing are some of these issues. Herein, using model LiNbO3 systems and two types of interdigitated transducers, we analyze the de-icing and anti-icing efficiencies and mechanisms driven by Rayleigh surface acoustic waves (R-SAW) and Lamb waves with 120 and 510 um wavelengths, respectively. Through the experimental analysis of de-icing and active anti-icing processes and the finite element simulation of the AW generation, propagation, and interaction with small ice aggregates, we disclose that Lamb waves are more favorable than R-SAWs to induce de-icing and/or prevent the freezing of droplets. Prospects for applications of this study are supported by proof of concept experiments, including de-icing in an ice wind tunnel, demonstrating that Lamb waves can efficiently remove ice layers covering large LN substrates. Results indicate that the de-icing mechanism may differ for Lamb waves or R-SAWs and that the wavelength must be considered as an important parameter for controlling the efficiency.

Surrogate models that combine dimensionality reduction and regression techniques are essential to reduce the need for costly high-fidelity computational fluid dynamics data. New approaches using $\beta$-Variational Autoencoder ($\beta$-VAE) architectures have shown promise in obtaining high-quality low-dimensional representations of high-dimensional flow data while enabling physical interpretation of their latent spaces. We propose a surrogate model based on latent space regression to predict pressure distributions on a transonic wing given the flight conditions: Mach number and angle of attack. The $\beta$-VAE model, enhanced with Principal Component Analysis (PCA), maps high-dimensional data to a low-dimensional latent space, showing a direct correlation with flight conditions. Regularization through $\beta$ requires careful tuning to improve overall performance, while PCA preprocessing helps to construct an effective latent space, improving autoencoder training and performance. Gaussian Process Regression is used to predict latent space variables from flight conditions, showing robust behavior independent of $\beta$, and the decoder reconstructs the high-dimensional pressure field data. This pipeline provides insight into unexplored flight conditions. Furthermore, a fine-tuning process of the decoder further refines the model, reducing the dependence on $\beta$ and enhancing accuracy. Structured latent space, robust regression performance, and significant improvements in fine-tuning collectively create a highly accurate and efficient surrogate model. Our methodology demonstrates the effectiveness of $\beta$-VAEs for aerodynamic surrogate modeling, offering a rapid, cost-effective, and reliable alternative for aerodynamic data prediction., Comment: 19 pages, 12 figures

Symmetry-protected topological phases have fundamentally changed our understanding of quantum matter. An archetypal example of such a quantum phase of matter is the Haldane phase, containing the spin-1 Heisenberg chain. The intrinsic quantum nature of such phases, however, often makes it challenging to study them using classical means. Here, we use trapped-ion qutrits to natively engineer spin-1 chains within the Haldane phase. Using a scalable, deterministic procedure to prepare the Affleck-Kennedy-Lieb-Tasaki (AKLT) state within the Haldane phase, we study the topological features of this system on a qudit quantum processor. Notably, we verify the long-range string order of the state, despite its short-range correlations, and observe spin fractionalization of the physical spin-1 particles into effective qubits at the chain edges, a defining feature of this system. The native realization of Haldane physics on a qudit quantum processor and the scalable preparation procedures open the door to the efficient exploration of a wide range of systems beyond spin-1/2, Comment: 12 pages, 6 figures

Quantum phase transitions reveal deep insights into the behavior of many-body quantum systems, but identifying these transitions without well-defined order parameters remains a significant challenge. In this work, we introduce a novel approach to constructing phase diagrams using the vector field of the reduced fidelity susceptibility (RFS). This method maps quantum phases and formulates an optimization problem to discover observables corresponding to order parameters. We demonstrate the effectiveness of our approach by applying it to well-established models, including the Axial Next Nearest Neighbour Interaction (ANNNI) model, a cluster state model, and a chain of Rydberg atoms. By analyzing observable decompositions into eigen-projectors and finite-size scaling, our method successfully identifies order parameters and characterizes quantum phase transitions with high precision. Our results provide a powerful tool for exploring quantum phases in systems where conventional order parameters are not readily available.

Children of undocumented immigrants exist in permanent limbo. They are fearful every day that their parents will be deported and wondering if they have any future in the country. Many of these students identify as Latinx. Instead of focusing on their education, these students struggle with this uncertainty and as a result are often absent from school or inattentive. Their teachers also struggle to motivate them and sometimes to protect them. The so-called broken immigration system hurts schools and creates victims across the spectrum of race and ethnicity in the United States, but it is especially acute for these students. This collaborative research brief updates and builds on a previous analysis of a 2017-2018 survey examining the harmful impacts of immigration enforcement on Latinx children of undocumented immigrants. It summarizes key findings from "Schools Under Siege: The Impact of Immigration Enforcement on Educational Equity," but updates the analysis with more recent data. [Additional collaborators of the research brief are the Center for the Transformation of Schools, School of Education & Information Studies, UCLA and the Latino Policy & Politics Institute, Meyer and Renee Luskin School of Public Affairs, UCLA.]