Search

We study quasi-periodic eigenvalue problems that arise in the stability analysis of periodic traveling wave solutions to Hamiltonian PDEs. We establish bounds on regions in the complex plane when the eigenvalues may deviate from the imaginary axis, and estimates for the number of such off-axis eigenvalues. These relations hold when the dispersion relation grows sufficiently rapidly in the wavenumber. The proofs involve a Gershgorin disk argument together with the Hamiltonian symmetry of the spectrum. The results are applicable to a broad class of nonlinear dispersive equations including the generalized Korteweg--de Vries, Benjamin--Bona--Mahoney, and Kawanhara equations., Comment: 20 pages; 3 figures

We report here the results obtained from a systematic optical photometric survey aimed at finding new compact binary millisecond pulsars (also known as "spiders"): the COmpact BInary PULsar SEarch (COBIPULSE). We acquired multi-band optical images over one year around $33$ unidentified Fermi-LAT sources, selected as pulsar candidates based on their curved GeV spectra and steady $\gamma$-ray emission. We present the discovery of four optical variables coinciding with the Fermi sources 3FGL J0737.2$-$3233, 3FGL J2117.6$+$3725 (two systems in this field) and 3FGL J2221.6$+$6507, which we propose as new candidate spider systems. Indeed, they all show optical flux modulation consistent with orbital periods of $0.3548(5) \ \mathrm{d}$, $0.25328(6) \ \mathrm{d}$, $0.441961(2) \ \mathrm{d}$, and $0.165(4) \ \mathrm{d}$, respectively, with amplitudes $\gtrsim 0.3 \ \mathrm{mag}$ and colors compatible with companion star temperatures of $5000$--$6000 \ \mathrm{K}$. These properties are consistent with the "redback" sub-class of spider pulsars. If confirmed as a millisecond pulsar, 3FGL J0737.2$-$3233 will be the closest known spider to Earth ($D=659_{-20}^{+16} \ \mathrm{pc}$, from Gaia-DR3 parallax). We searched and did not find any X-ray sources matching our four candidates, placing $3\sigma$ upper limits of $\sim10^{31}$--$10^{32} \ \mathrm{erg} \ \mathrm{s}^{-1}$ ($0.3$--$10 \ \mathrm{keV}$) on their soft X-ray luminosities. We also present and discuss other multi-wavelength information on our spider candidates, from infrared to X-rays., Comment: 38 pages, 93 figures, 4 tables. This paper has been accepted for publication in ApJ

The Zirconium (Z = 40) isotopic chain has attracted interest for more than four decades. The abrupt lowering of the energy of the first $2^+$ state and the increase in the transition strength B(E2; $2_1^\rightarrow 0_1^+$ going from $^{98}$Zr to $^{100}$Zr has been the first example of "quantum phase transition" in nuclear shapes, which has few equivalents in the nuclear chart. Although a multitude of experiments have been performed to measure nuclear properties related to nuclear shapes and collectivity in the region, none of the measured lifetimes were obtained using the Recoil Distance Doppler Shift method in the $\gamma\gamma$-coincidence mode where a gate on the direct feeding transition of the state of interest allows a strict control of systematical errors. This work reports the results of lifetime measurements for the first yrast excited states in $^{98-104}$Zr carried out to extract reduced transition probabilities. The new lifetime values in $\gamma\gamma$-coincidence and $\gamma$-single mode are compared with the results of former experiments. Recent predictions of the Interacting Boson Model with Configuration Mixing, the Symmetry Conserving Configuration Mixing model based on the Hartree-Fock-Bogoliubov approach and the Monte Carlo Shell Model are presented and compared with the experimental data.

Let $W$ be a finite group generated by reflections of a lattice $M$. If a lattice polytope $P \subset M \otimes_{\mathbb Z}\mathbb R$ is preserved by $W$, then we show that the quotient of the projective toric variety $X_P$ by $W$ is isomorphic to the toric variety $X_{P \cap D}$, where $D$ is a fundamental domain for the action of $W$. This answers a question of Horiguchi-Masuda-Shareshian-Song, and recovers results of Blume, of the second author, and of Gui-Hu-Liu., Comment: comments welcome!

Laser-driven free-electron lasers (LDFELs) replace magnetostatic undulators with the electromagnetic fields of a laser pulse. Because the undulator period is half the wavelength of the laser pulse, LDFELs can amplify x rays using lower electron energies and over shorter interaction lengths than a conventional free-electron laser. Here we show that a flying-focus pulse substantially reduces the energy required to reach high gain in an LDFEL by providing a highly uniform, high-intensity field over the entire interaction length. The flying-focus pulse features an intensity peak that travels in the opposite direction of its phase fronts. This enables an LDFEL configuration where an electron beam collides head-on with the phase fronts and experiences a near-constant undulator strength as it co-propagates with the intensity peak. Three-dimensional simulations of this configuration demonstrate the generation of megawatts of coherent x-ray radiation with 20 times less energy than a conventional laser pulse.

The lattice of flats $\mathcal L_M$ of a matroid $M$ is combinatorially well-behaved and, when $M$ is realizable, admits a geometric model in the form of a "matroid Schubert variety". In contrast, the lattice of flats of a polymatroid exhibits many combinatorial pathologies and admits no similar geometric model. We address this situation by defining the lattice $\mathcal L_P$ of "combinatorial flats" of a polymatroid $P$. Combinatorially, $\mathcal L_P$ exhibits good behavior analogous to that of $\mathcal L_M$: it is graded, determines $P$ when $P$ is simple, and is top-heavy. When $P$ is realizable over a field of characteristic 0, we show that $\mathcal L_P$ is modelled by a "polymatroid Schubert variety". Our work generalizes a number of results of Ardila-Boocher and Huh-Wang on matroid Schubert varieties; however, the geometry of polymatroid Schubert varieties is noticeably more complicated than that of matroid Schubert varieties. Many natural questions remain open., Comment: 24 pages, 3 figures, comments welcome!

In a hierarchically formed Universe, galaxies accrete smaller systems that tidally disrupt as they evolve in the host's potential. We present a complete catalogue of disrupting galaxies accreted onto Milky Way-mass haloes from the Auriga suite of cosmological magnetohydrodynamic zoom-in simulations. We classify accretion events as intact satellites, stellar streams, or phase-mixed systems based on automated criteria calibrated to a visually classified sample, and match accretions to their counterparts in haloes re-simulated at higher resolution. Most satellites with a bound progenitor at the present day have lost substantial amounts of stellar mass -- 67 per cent have $f_\text{bound} < 0.97$ (our threshold to no longer be considered intact), while 53 per cent satisfy a more stringent $f_\text{bound} < 0.8$. Streams typically outnumber intact systems, contribute a smaller fraction of overall accreted stars, and are substantial contributors at intermediate distances from the host centre ($\sim$0.1 to $\sim$0.7$R_\text{200m}$, or $\sim$35 to $\sim$250 kpc for the Milky Way). We also identify accretion events that disrupt to form streams around massive intact satellites instead of the main host. Streams are more likely than intact or phase-mixed systems to have experienced preprocessing, suggesting this mechanism is important for setting disruption rates around Milky Way-mass haloes. All of these results are preserved across different simulation resolutions, though we do find some hints that satellites disrupt more readily at lower resolution. The Auriga haloes suggest that disrupting satellites surrounding Milky Way-mass galaxies are the norm and that a wealth of tidal features waits to be uncovered in upcoming surveys., Comment: 16+4 pages, 13+2 figures, 1+1 tables. Submitted to MNRAS, comments are welcome

Galaxies like the Milky Way are surrounded by complex populations of satellites at all stages of tidal disruption. In this paper, we present a dynamical study of the disrupting satellite galaxies in the Auriga simulations that are orbiting 28 distinct Milky Way-mass hosts across three resolutions. We find that the satellite galaxy populations are highly disrupted. The majority of satellites that remain fully intact at present day were accreted recently without experiencing more than one pericentre ($n_{\rm peri} \lesssim 1$) and have large apocentres ($r_{\rm apo} \gtrsim 200$ kpc) and pericentres ($r_{\rm peri} \gtrsim 50$ kpc). The remaining satellites have experienced significant tidal disruption and, given full knowledge of the system, would be classified as stellar streams. We find stellar streams in Auriga across the range of pericentres and apocentres of the known Milky Way dwarf galaxy streams and, interestingly, overlapping significantly with the Milky Way intact satellite population. We find no significant change in satellite orbital distributions across resolution. However, we do see substantial halo-to-halo variance of $(r_\text{peri}, r_\text{apo})$ distributions across host galaxies, as well as a dependence of satellite orbits on host halo mass - systems disrupt at larger pericentres and apocentres in more massive hosts. Our results suggest that either cosmological simulations (including, but not limited to, Auriga) are disrupting satellites far too readily, or that the Milky Way's satellites are more disrupted than current imaging surveys have revealed. Future observing facilities and careful mock observations of these systems will be key to revealing the nature of this apparent discrepancy., Comment: 16 pages, 10 figures, 1 table, submitted to MNRAS

The novel series of hole-doped quadruple manganite perovskites Hg$_{1-x}$Na$_x$Mn$_3$Mn$_4$O$_{12}$ (HNMO) has been synthesized and its charge and orbital order behavior investigated through high-resolution synchrotron powder x-ray diffraction techniques. Through careful Rietveld refinements of structural models $via$ symmetry-motivated approaches, we show that the ground state of HNMO compositions adopts a polar $G$-type charge and orbital ordered state, which is rare in manganite perovskites, and is robust as a sole phase up to a critical doping level. Upon this critical doping, coincident with that in which colossal magnetoresistance (CMR) is maximal in canonical manganite perovskites, electronic phase separation occurs between $G$-type and orbital order with charge disorder-type states. The latter state has recently been identified in Ca$_{1-x}$Na$_x$Mn$_3$Mn$_4$O$_{12}$ perovskites, and proposed to be the competing insulating state from which CMR phenomena emerges. We show the mechanism for the formation of the $G$-type state is due to charge transfer processes which may occur through a coupling of distortions involving structural and charge and orbital degrees of freedom, ultimately driving the polar ground state through an improper-like ferroelectric polarization mechanism. These results will act as an important recipe for designing novel ferroelectric-active materials, in addition to expanding the richness of charge and orbital ordered states in manganite perovskites., Comment: 11 pages, 8 figures

Importance: Many individuals with drug-resistant epilepsy continue to have seizures after resective surgery. Accurate identification of focal brain abnormalities is essential for successful neurosurgical intervention. Current clinical approaches to identify structural abnormalities for surgical targeting in epilepsy do not use diffusion-weighted MRI (dMRI), despite evidence that dMRI abnormalities are present in epilepsy and may relate to the epileptogenic zone. Objective: To investigate whether surgical resection of diffusion abnormalities relates to post-operative seizure freedom. Design: This retrospective case-control study was conducted between 2009 and 2022. Data were acquired at the National Hospital for Neurology and Neurosurgery, UK. Study participants included 200 individuals with drug-resistant focal epilepsy, who underwent resective surgery, and 97 healthy controls used as a normative baseline. Main Outcomes: Spatial overlap between diffusion abnormality clusters and surgical resection masks, and relation to post-surgical outcome. Results: Surgical resections overlapping with the largest abnormal cluster significantly correlated with sustained seizure freedom at 12 months (83% vs 55%; p<0.0001) and over five years (p<0.0001). Notably, resecting only a small proportion of the largest cluster was associated with better seizure outcomes than cases with no resection of this cluster (p=0.008). Furthermore, sparing the largest cluster but resecting other large clusters still improved seizure freedom rates compared to no overlap (p=0.03). Conclusions: Our results suggest that abnormal clusters, identified using dMRI, are integral to the epileptogenic network, and even a partial removal of such an abnormal cluster is sufficient to achieve seizure freedom. The study highlights the potential of incorporating dMRI into pre-surgical planning to improve outcomes in focal epilepsy.

The stars of the Milky Way carry the chemical history of our Galaxy in their atmospheres as they journey through its vast expanse. Like barcodes, we can extract the chemical fingerprints of stars from high-resolution spectroscopy. The fourth data release (DR4) of the Galactic Archaeology with HERMES (GALAH) Survey, based on a decade of observations, provides the chemical abundances of up to 32 elements for 917 588 stars that also have exquisite astrometric data from the $Gaia$ satellite. For the first time, these elements include life-essential nitrogen to complement carbon, and oxygen as well as more measurements of rare-earth elements critical to modern-life electronics, offering unparalleled insights into the chemical composition of the Milky Way. For this release, we use neural networks to simultaneously fit stellar parameters and abundances across the full spectrum, leveraging synthetic grids computed with Spectroscopy Made Easy. These grids account for atomic line formation in non-local thermodynamic equilibrium for 14 elements. In a two-iteration process, we first fit stellar labels for all 1 085 520 spectra, then co-add repeated observations and refine these labels using astrometric data from $Gaia$ and 2MASS photometry, improving the accuracy and precision of stellar parameters and abundances. Our validation thoroughly assesses the reliability of spectroscopic measurements and highlights key caveats for catalogue users. GALAH DR4 represents yet another milestone in Galactic archaeology, combining detailed chemical compositions from multiple nucleosynthetic channels with kinematic information and age estimates. The resulting dataset, covering nearly a million stars, opens new avenues for understanding not only the chemical and dynamical history of the Milky Way, but also the broader questions of the origin of elements and the evolution of planets, stars, and galaxies., Comment: 43 pages, 38 figures to be submitted to PASA. Accompanying the GALAH Data Release 4, see https://www.galah-survey.org and https://cloud.datacentral.org.au/teamdata/GALAH/public/GALAH_DR4/. All code available on http://github.com/svenbuder/GALAH_DR4/ and https://github.com/svenbuder/galah_dr4_paper. Comments welcome

Medical imaging is spearheading the AI transformation of healthcare. Performance reporting is key to determine which methods should be translated into clinical practice. Frequently, broad conclusions are simply derived from mean performance values. In this paper, we argue that this common practice is often a misleading simplification as it ignores performance variability. Our contribution is threefold. (1) Analyzing all MICCAI segmentation papers (n = 221) published in 2023, we first observe that more than 50% of papers do not assess performance variability at all. Moreover, only one (0.5%) paper reported confidence intervals (CIs) for model performance. (2) To address the reporting bottleneck, we show that the unreported standard deviation (SD) in segmentation papers can be approximated by a second-order polynomial function of the mean Dice similarity coefficient (DSC). Based on external validation data from 56 previous MICCAI challenges, we demonstrate that this approximation can accurately reconstruct the CI of a method using information provided in publications. (3) Finally, we reconstructed 95% CIs around the mean DSC of MICCAI 2023 segmentation papers. The median CI width was 0.03 which is three times larger than the median performance gap between the first and second ranked method. For more than 60% of papers, the mean performance of the second-ranked method was within the CI of the first-ranked method. We conclude that current publications typically do not provide sufficient evidence to support which models could potentially be translated into clinical practice., Comment: Paper accepted at MICCAI 2024 conference

It is no secret that statistical modelling often involves making simplifying assumptions when attempting to study complex stochastic phenomena. Spatial modelling of extreme values is no exception, with one of the most common such assumptions being stationarity in the marginal and/or dependence features. If non-stationarity has been detected in the marginal distributions, it is tempting to try to model this while assuming stationarity in the dependence, without necessarily putting this latter assumption through thorough testing. However, margins and dependence are often intricately connected and the detection of non-stationarity in one feature might affect the detection of non-stationarity in the other. This work is an in-depth case study of this interrelationship, with a particular focus on a spatio-temporal environmental application exhibiting well-documented marginal non-stationarity. Specifically, we compare and contrast four different marginal detrending approaches in terms of our post-detrending ability to detect temporal non-stationarity in the spatial extremal dependence structure of a sea surface temperature dataset from the Red Sea., Comment: 18 pages, 7 figures

Cutting-edge abstractive summarisers generate fluent summaries, but the factuality of the generated text is not guaranteed. Early summary factuality evaluation metrics are usually based on n-gram overlap and embedding similarity, but are reported fail to align with human annotations. Therefore, many techniques for detecting factual inconsistencies build pipelines around natural language inference (NLI) or question-answering (QA) models with additional supervised learning steps. In this paper, we revisit similarity-based metrics, showing that this failure stems from the comparison text selection and its granularity. We propose a new zero-shot factuality evaluation metric, Sentence-BERT Score (SBERTScore), which compares sentences between the summary and the source document. It outperforms widely-used word-word metrics including BERTScore and can compete with existing NLI and QA-based factuality metrics on the benchmark without needing any fine-tuning. Our experiments indicate that each technique has different strengths, with SBERTScore particularly effective in identifying correct summaries. We demonstrate how a combination of techniques is more effective in detecting various types of error.

Band gap energy of an organic molecule can be reduced by intermolecular interaction enhancement, and thus, certain polycyclic aromatic hydrocarbons (PAHs), which are insulators with wide band gaps, are expected to undergo insulator-metal transitions by simple compression. Such a pressure-induced electronic transition can be exploited to transform non-metallic organic materials into states featuring intriguing electronic characteristics such as high-temperature superconductivity. Numerous attempts have been made to metalize various small PAHs, but so far only pressure-induced amorphization well below the megabar region was observed. The wide band gap energy of the small PAHs and low chemical stability under simple compression are the bottlenecks. We have investigated the band gap energy evolution and the crystal structural compression of the large PAH molecules, where the band gap energy is significantly reduced by increasing the number of {\pi}-electrons and improved chemical stability with fully benzenoid molecular structure. Herein, we present a pressure-induced transition in dicoronylene, C48H20, an insulator at ambient conditions that transforms into a semi-metallic state above 23.0 GPa with a three-order-of-magnitude reduction in resistivity. In-situ UV-visible absorption, transport property measurement, Raman spectroscopy, X-ray diffraction and density functional theory calculations were performed to provide tentative explanations to the alterations in its electronic structure at high pressure. The discovery of an electronic transition at pressures well below the megabar is a promising step towards realization of a single component purely hydrocarbon molecular metal in the near future., Comment: 22 pages, 5 figures

We introduce Clerical, a programming language for exact real-number computation that combines first-order imperative-style programming with a limit operator for computation of real numbers as limits of Cauchy sequences. We address the semidecidability of the linear ordering of the reals by incorporating nondeterministic guarded choice, through which decisions based on partial comparison operations on reals can be patched together to give total programs. The interplay between mutable state, nondeterminism, and computation of limits is controlled by the requirement that expressions computing limits and guards modify only local state. We devise a domain-theoretic denotational semantics that uses a variant of Plotkin powerdomain construction tailored to our specific version of nondeterminism. We formulate a Hoare-style specification logic, show that it is sound for the denotational semantics, and illustrate the setup by implementing and proving correct a program for computation of $\pi$ as the least positive zero of $\sin$. The modular character of Clerical allows us to compose the program from smaller parts, each of which is shown to be correct on its own. We provide a proof-of-concept OCaml implementation of Clerical, and formally verify parts of the development, notably the soundness of specification logic, in the Coq proof assistant.

The Gemini South telescope is now equipped with a new high-resolution spectrograph called GHOST (the Gemini High-resolution Optical SpecTrograph). This instrument provides high-efficiency, high-resolution spectra covering 347-1060 nm in a single exposure of either one or two targets simultaneously, along with precision radial velocity spectroscopy utilizing an internal calibration source. It can operate at a spectral element resolving power of either 76000 or 56000, and can reach a SNR$\sim$5 in a 1hr exposure on a V$\sim$20.8 mag target in median site seeing, and dark skies (per resolution element). GHOST was installed on-site in June 2022, and we report performance after full integration to queue operations in November 2023, in addition to scientific results enabled by the integration observing runs. These results demonstrate the ability to observe a wide variety of bright and faint targets with high efficiency and precision. With GHOST, new avenues to explore high-resolution spectroscopy have opened up to the astronomical community. These are described, along with the planned and potential upgrades to the instrument., Comment: Accepted in the Astronomical Journal; 26 pages, 24 figures and 4 tables

We consider the problem of direct data-driven predictive control for unknown stochastic linear time-invariant (LTI) systems with partial state observation. Building upon our previous research on data-driven stochastic control, this paper (i) relaxes the assumption of Gaussian process and measurement noise, and (ii) enables optimization of the gain matrix within the affine feedback policy. Output safety constraints are modelled using conditional value-at-risk, and enforced in a distributionally robust sense. Under idealized assumptions, we prove that our proposed data-driven control method yields control inputs identical to those produced by an equivalent model-based stochastic predictive controller. A simulation study illustrates the enhanced performance of our approach over previous designs., Comment: 8 pages, 1 figure, 2 tables, the extended version of an accepted paper in Conference on Decision and Control (CDC). arXiv admin note: text overlap with arXiv:2312.15177

Thermoelectric materials can be designed to support optical resonances across multiple spectral ranges to enable ultra-wide band photodetection. For instance, antimony telluride (Sb2Te3) chalcogenide exhibits interband plasmonic resonances in the visible range and Mie resonances in the mid-infrared (mid-IR) range, while simultaneously possessing large thermoelectric Seebeck coefficients. In this paper, we designed and fabricated Sb2Te3 metasurface devices to achieve resonant absorption for enabling photodetectors operating across an ultra-wideband spectrum, from visible to mid-IR. Furthermore, relying on asymmetric Sb2Te3 metasurface, we demonstrated the thermoelectric photodetectors with polarization-selectivity. This work provides a potential platform towards the portable ultrawide band spectrometers at room temperature, for environmental sensing applications.

Detecting out-of-distribution (OOD) data is crucial in machine learning applications to mitigate the risk of model overconfidence, thereby enhancing the reliability and safety of deployed systems. The majority of existing OOD detection methods predominantly address uni-modal inputs, such as images or texts. In the context of multi-modal documents, there is a notable lack of extensive research on the performance of these methods, which have primarily been developed with a focus on computer vision tasks. We propose a novel methodology termed as attention head masking (AHM) for multi-modal OOD tasks in document classification systems. Our empirical results demonstrate that the proposed AHM method outperforms all state-of-the-art approaches and significantly decreases the false positive rate (FPR) compared to existing solutions up to 7.5\%. This methodology generalizes well to multi-modal data, such as documents, where visual and textual information are modeled under the same Transformer architecture. To address the scarcity of high-quality publicly available document datasets and encourage further research on OOD detection for documents, we introduce FinanceDocs, a new document AI dataset. Our code and dataset are publicly available.

We present GTC-OSIRIS phase-resolved optical spectroscopy of three compact binary MSPs, or `spiders': PSR J1048+2339, PSR J1810+1744, and (for the first time) PSR J1908+2105. For the companion in each system, the temperature is traced throughout its orbit, and radial velocities are measured. The radial velocities are found to vary with the absorption features used when measuring them, resulting in a lower radial velocity curve semi-amplitude measured from the day side of two of the systems when compared to the night: for J1048 ($K_\mathrm{day} = 344 \pm 4$ km s$^{-1}$, $K_\mathrm{night} = 372 \pm 3$ km s$^{-1}$) and, tentatively, for J1810 ($K_\mathrm{day} = 448 \pm 19$ km s$^{-1}$, $K_\mathrm{night} = 491 \pm 32$ km s$^{-1}$). With existing inclination constraints, this gives the neutron star (NS) and companion masses $M_\mathrm{NS} = 1.50 - 2.04$ $M_\odot$ and $M_2 = 0.32 - 0.40$ $M_\odot$ for J1048, and $M_\mathrm{NS} > 1.7$ $M_\odot$ and $M_2 = 0.05 - 0.08$ $M_\odot$ for J1810. For J1908, we find an upper limit of $K_2 < 32$ km s$^{-1}$, which constrains its mass ratio $q = M_2 / M_\mathrm{NS} > 0.55$ and inclination $i < 6.0^\circ$, revealing the previously misunderstood system to be the highest mass ratio, lowest inclination redback yet. This raises questions for the origins of its substantial radio eclipses. Additionally, we find evidence of asymmetric heating in J1048 and J1810, and signs of metal enrichment in J1908. We also explore the impact of inclination on spectroscopic temperatures, and demonstrate that the temperature measured at quadrature ($\phi = 0.25, 0.75$) is essentially independent of inclination, and thus can provide additional constraints on photometric modelling., Comment: Submitted to MNRAS. 18 pages, 19 figures

While many capabilities of language models (LMs) improve with increased training budget, the influence of scale on hallucinations is not yet fully understood. Hallucinations come in many forms, and there is no universally accepted definition. We thus focus on studying only those hallucinations where a correct answer appears verbatim in the training set. To fully control the training data content, we construct a knowledge graph (KG)-based dataset, and use it to train a set of increasingly large LMs. We find that for a fixed dataset, larger and longer-trained LMs hallucinate less. However, hallucinating on $\leq5$% of the training data requires an order of magnitude larger model, and thus an order of magnitude more compute, than Hoffmann et al. (2022) reported was optimal. Given this costliness, we study how hallucination detectors depend on scale. While we see detector size improves performance on fixed LM's outputs, we find an inverse relationship between the scale of the LM and the detectability of its hallucinations., Comment: Published at COLM 2024. 16 pages, 11 figures

Certain linear matrix operators arise naturally in systems analysis and design problems involving cascade interconnections of linear time-invariant systems, including problems of stabilization, estimation, and model order reduction. We conduct here a comprehensive study of these operators and their relevant system-theoretic properties. The general theory is then leveraged to delineate both known and new design methodologies for control, estimation, and model reduction. Several entirely new designs arise from this systematic categorization, including new recursive and low-gain design frameworks for observation of cascaded systems. The benefits of the results beyond the linear time-invariant setting are demonstrated through preliminary extensions for nonlinear systems, with an outlook towards the development of a similarly comprehensive nonlinear theory., Comment: 16 pages, 5 figures, submitted for publication

Status epilepticus (SE) carries risks of morbidity and mortality. Experimental studies have implicated the entorhinal cortex in prolonged seizures; however, studies in large human cohorts are limited. We hypothesised that individuals with temporal lobe epilepsy (TLE) and a history of SE would have more severe entorhinal atrophy compared to others with TLE and no history of SE. 357 individuals with drug resistant temporal lobe epilepsy (TLE) and 100 healthy controls were scanned on a 3T MRI. For all subjects the cortex was segmented, parcellated, and the thickness calculated from the T1-weighted anatomical scan. Subcortical volumes were derived similarly. Cohen's d and Wilcoxon rank-sum tests respectively were used to capture effect sizes and significance. Individuals with TLE and SE had reduced entorhinal thickness compared to those with TLE and no history of SE. The entorhinal cortex was more atrophic ipsilaterally (d=0.51, p<0.001) than contralaterally (d=0.37, p=0.01). Reductions in ipsilateral entorhinal thickness were present in both left TLE (n=22:176, d=0.78, p<0.001), and right TLE (n=19:140, d=0.31, p=0.04), albeit with a smaller effect size in right TLE. Several other regions exhibited atrophy in individuals with TLE, but these did not relate to a history of SE. These findings suggest potential involvement or susceptibility of the entorhinal cortex in prolonged seizures.

The border-collision normal form is a piecewise-linear family of continuous maps that describe the dynamics near border-collision bifurcations. Most prior studies assume each piece of the normal form is invertible, as is generic from an abstract viewpoint, but in applied problems one piece of the map often has degenerate range, corresponding to a zero determinant. This provides simplification, yet even in two dimensions the dynamics can be incredibly rich. The purpose of this paper is to determine broadly how the dynamics of the two-dimensional border-collision normal form with a zero determinant differs for different values of its parameters. We identify parameter regions of period-adding, period-incrementing, mode-locking, and component doubling of chaotic attractors, and characterise the dominant bifurcation boundaries. The intention is for the results to enable border-collision bifurcations in mathematical models to be analysed more easily and effectively, and we illustrate this with a flu epidemic model and two stick-slip friction oscillator models. We also describe three novel bifurcation structures that remain to be explored.

CT is a main modality for imaging liver diseases, valuable in detecting and localizing liver tumors. Traditional anomaly detection methods analyze reconstructed images to identify pathological structures. However, these methods may produce suboptimal results, overlooking subtle differences among various tissue types. To address this challenge, here we employ generative diffusion prior to inpaint the liver as the reference facilitating anomaly detection. Specifically, we use an adaptive threshold to extract a mask of abnormal regions, which are then inpainted using a diffusion prior to calculating an anomaly score based on the discrepancy between the original CT image and the inpainted counterpart. Our methodology has been tested on two liver CT datasets, demonstrating a significant improvement in detection accuracy, with a 7.9% boost in the area under the curve (AUC) compared to the state-of-the-art. This performance gain underscores the potential of our approach to refine the radiological assessment of liver diseases., Comment: 10 pages, 11 figures

We report the detection of 239 trans-Neptunian Objects discovered through the on-going New Horizons survey for distant minor bodies being performed with the Hyper Suprime-Cam mosaic imager on the Subaru Telescope. These objects were discovered in images acquired with either the r2 or the recently commissioned EB-gri filter using shift and stack routines. Due to the extremely high stellar density of the search region down stream of the spacecraft, new machine learning techniques had to be developed to manage the extremely high false positive rate of bogus candidates produced from the shift and stack routines. We report discoveries as faint as r2$\sim26.5$. We highlight an overabundance of objects found at heliocentric distances $R\gtrsim70$~au compared to expectations from modelling of the known outer Solar System. If confirmed, these objects betray the presence of a heretofore unrecognized abundance of distant objects that can help explain a number of other observations that otherwise remain at odds with the known Kuiper Belt, including detections of serendipitous stellar occultations, and recent results from the Student Dust Counter on-board the New Horizons spacecraft., Comment: Accepted for publication in the Planetary Science Journal, 28 pages, 7 figures, 3 tables

Reduced Order Models (ROMs) form essential tools across engineering domains by virtue of their function as surrogates for computationally intensive digital twinning simulators. Although purely data-driven methods are available for ROM construction, schemes that allow to retain a portion of the physics tend to enhance the interpretability and generalization of ROMs. However, physics-based techniques can adversely scale when dealing with nonlinear systems that feature parametric dependencies. This study introduces a generative physics-based ROM that is suited for nonlinear systems with parametric dependencies and is additionally able to quantify the confidence associated with the respective estimates. A main contribution of this work is the conditioning of these parametric ROMs to features that can be derived from monitoring measurements, feasibly in an online fashion. This is contrary to most existing ROM schemes, which remain restricted to the prescription of the physics-based, and usually a priori unknown, system parameters. Our work utilizes conditional Variational Autoencoders to continuously map the required reduction bases to a feature vector extracted from limited output measurements, while additionally allowing for a probabilistic assessment of the ROM-estimated Quantities of Interest. An auxiliary task using a neural network-based parametrization of suitable probability distributions is introduced to re-establish the link with physical model parameters. We verify the proposed scheme on a series of simulated case studies incorporating effects of geometric and material nonlinearity under parametric dependencies related to system properties and input load characteristics.

Adversarial attacks on explainability models have drastic consequences when explanations are used to understand the reasoning of neural networks in safety critical systems. Path methods are one such class of attribution methods susceptible to adversarial attacks. Adversarial learning is typically phrased as a constrained optimisation problem. In this work, we propose algebraic adversarial examples and study the conditions under which one can generate adversarial examples for integrated gradients. Algebraic adversarial examples provide a mathematically tractable approach to adversarial examples.

The border-collision normal form describes the local dynamics in continuous systems with switches when a fixed point intersects a switching surface. For one-dimensional cases where the bifurcation creates or destroys only fixed points and period-two orbits, we show that the standard local equivalence of normal forms, topological conjugacy, can be replaced by differentiable conjugacy provided an extra term is added to the normal form. In these cases topological conjugacy is so weak that a range of values can be used for the coefficients in the normal form. The extension to differentiable conjugacy explains why the usual choice of slopes in the standard normal form is privileged. This highlights the importance of differentiable conjugacies and the need for extended normal forms.

As the number of scientific publications and preprints is growing exponentially, several attempts have been made to navigate this complex and increasingly detailed landscape. These have almost exclusively taken unsupervised approaches that fail to incorporate domain knowledge and lack the structural organisation required for intuitive interactive human exploration and discovery. Especially in highly interdisciplinary fields, a deep understanding of the connectedness of research works across topics is essential for generating insights. We have developed a unique approach to data navigation that leans on geographical visualisation and uses hierarchically structured domain knowledge to enable end-users to explore knowledge spaces grounded in their desired domains of interest. This can take advantage of existing ontologies, proprietary intelligence schemata, or be directly derived from the underlying data through hierarchical topic modelling. Our approach uses natural language processing techniques to extract named entities from the underlying data and normalise them against relevant domain references and navigational structures. The knowledge is integrated by first calculating similarities between entities based on their shared extracted feature space and then by alignment to the navigational structures. The result is a knowledge graph that allows for full text and semantic graph query and structured topic driven navigation. This allows end-users to identify entities relevant to their needs and access extensive graph analytics. The user interface facilitates graphical interaction with the underlying knowledge graph and mimics a cartographic map to maximise ease of use and widen adoption. We demonstrate an exemplar project using our generalisable and scalable infrastructure for an academic biomedical literature corpus that is grounded against hundreds of different named domain entities.

We say a finite word $x$ is a palindromic periodicity if there exist two palindromes $p$ and $s$ such that $|x| \geq |ps|$ and $x$ is a prefix of the word $(ps)^\omega = pspsps\cdots$. In this paper we examine the palindromic periodicities occurring in some classical infinite words, such as Sturmian words, episturmian words, the Thue-Morse word, the period-doubling word, the Rudin-Shapiro word, the paperfolding word, and the Tribonacci word, and prove a number of results about them. We also prove results about words with the smallest number of palindromic periodicities., Comment: revised version with minor corrections and more details of proofs

With high equilibrium temperatures and tidally locked rotation, ultra-hot Jupiters (UHJs) are unique laboratories within which to probe extreme atmospheric physics and chemistry. In this paper, we present high-resolution dayside spectroscopy of the UHJ WASP-189b obtained with the new Gemini High-resolution Optical SpecTrograph (GHOST) at the Gemini South Observatory. The observations, which cover three hours of post-eclipse orbital phases, were obtained during the instrument's System Verification run. We detect the planet's atmosphere via the Doppler cross-correlation technique, and recover a detection of neutral iron in the planet's dayside atmosphere at a significance of 7.5$\sigma$ in the red-arm of the data, verifying the presence of a thermal inversion. We also investigate the presence of other species in the atmosphere and discuss the implications of model injection/recovery tests. These results represent the first atmospheric characterization of an exoplanet with GHOST's high-resolution mode, and demonstrate the potential of this new instrument in detecting and studying ultra-hot exoplanet atmospheres., Comment: 21 pages, 12 figures, 3 tables. Accepted for publication in The Astronomical Journal

Automated Essay Scoring (AES) holds significant promise in the field of education, helping educators to mark larger volumes of essays and provide timely feedback. However, Arabic AES research has been limited by the lack of publicly available essay data. This study introduces AR-AES, an Arabic AES benchmark dataset comprising 2046 undergraduate essays, including gender information, scores, and transparent rubric-based evaluation guidelines, providing comprehensive insights into the scoring process. These essays come from four diverse courses, covering both traditional and online exams. Additionally, we pioneer the use of AraBERT for AES, exploring its performance on different question types. We find encouraging results, particularly for Environmental Chemistry and source-dependent essay questions. For the first time, we examine the scale of errors made by a BERT-based AES system, observing that 96.15 percent of the errors are within one point of the first human marker's prediction, on a scale of one to five, with 79.49 percent of predictions matching exactly. In contrast, additional human markers did not exceed 30 percent exact matches with the first marker, with 62.9 percent within one mark. These findings highlight the subjectivity inherent in essay grading, and underscore the potential for current AES technology to assist human markers to grade consistently across large classes.

The central regions of the Milky Way constitute a unique laboratory for a wide swath of astrophysical studies, consequently the inner $\sim$400 pc has been the target of numerous large surveys at all accessible wavelengths. In this paper we present a catalog of sources at 25 and 37 $\mu$m located within all of the regions observed with the SOFIA/FORCAST instrument in the inner $\sim$200 pc of the Galaxy. The majority of the observations were obtained as part of the SOFIA Cycle 7 Galactic Center Legacy program survey, which was designed to complement the Spitzer/MIPS 24 $\mu$m catalog in regions saturated in the MIPS observations. Due to the wide variety of source types captured by our observations at 25 and 37 $\mu$m, we do not limit the FORCAST source catalog to unresolved point sources, or treat all sources as if they are point-like sources. The catalog includes all detectable sources in the regions, resulting in a catalog of 950 sources, including point sources, compact sources, and extended sources. We also provide the user with metrics to discriminate between the source types., Comment: 29 pages, 13 figures, Accepted to ApJ

In our study, we first constructed a dataset from the tweets of the top 100 medical influencers with the highest Influencer Score during the COVID-19 pandemic. This dataset was then used to construct a socio-semantic network, mapping both their identities and key topics, which are crucial for understanding their impact on public health discourse. To achieve this, we developed a few-shot multi-label classifier to identify influencers and their network actors' identities, employed BERTopic for extracting thematic content, and integrated these components into a network model to analyze their impact on health discourse. To ensure the reproducibility of our results, we have made the code available at https://github.com/ZhijinGuo/Medinfluencer., Comment: ACM SIGKDD 2024 Workshop epiDAMIK 2024: The 7th International Workshop on Epidemiology meets Data Mining and Knowledge Discovery

In a post-industrial society, the workplace is dominated primarily by Knowledge Work, which is achieved mostly through human cognitive processing, such as analysis, comprehension, evaluation, and decision-making. Many of these processes have limited support from technology in the same way that physical tasks have been enabled through a host of tools from hammers to shovels and hydraulic lifts. To develop a suite of cognitive tools, we first need to understand which processes humans use to complete work tasks. In the past century several classifications (e.g., Blooms) of cognitive processes have emerged, and we assessed their viability as the basis for designing tools that support cognitive work. This study re-used an existing data set composed of interviews of environmental scientists about their core work. While the classification uncovered many instances of cognitive process, the results showed that the existing cognitive process classifications do not provide a sufficiently comprehensive deconstruction of the human cognitive processes; the work is quite simply too abstract to be operational.

Coral reefs are increasingly subjected to major disturbances threatening the health of marine ecosystems. Substantial research underway to develop intervention strategies that assist reefs in recovery from, and resistance to, inevitable future climate and weather extremes. To assess potential benefits of interventions, mechanistic understanding of coral reef recovery and resistance patterns is essential. Recent evidence suggests that more than half of the reefs surveyed across the Great Barrier Reef (GBR) exhibit deviations from standard recovery modelling assumptions when the initial coral cover is low ($\leq 10$\%). New modelling is necessary to account for these observed patterns to better inform management strategies. We consider a new model for reef recovery at the coral cover scale that accounts for biphasic recovery patterns. The model is based on a multispecies Richards' growth model that includes a change point in the recovery patterns. Bayesian inference is applied for uncertainty quantification of key parameters for assessing reef health and recovery patterns. This analysis is applied to benthic survey data from the Australian Institute of Marine Sciences (AIMS). We demonstrate agreement between model predictions and data across every recorded recovery trajectory with at least two years of observations following disturbance events occurring between 1992--2020. This new approach will enable new insights into the biological, ecological and environmental factors that contribute to the duration and severity of biphasic coral recovery patterns across the GBR. These new insights will help to inform managements and monitoring practice to mitigate the impacts of climate change on coral reefs.

We propose a simple mathematical model to describe the mechanical relaxation of cells within a curved epithelial tissue layer represented by an arbitrary curve in two-dimensional space. The model represents the mechanics of the cell body either by straight springs between points of the curve, or by curved springs whose shape follows the curve. To understand the collective behaviour of these discrete models of cells at the broader tissue scale, we devise an appropriate continuum limit in which the number of cells is constant but the number of springs tends to infinity. The continuum limit shows that (i)~the straight spring model and the curved spring model converge to the same dynamics; and (ii)~the density of cells becomes governed by a diffusion equation in arc length space with second-order accuracy, where diffusion may be linear or nonlinear depending on the choice of the spring restoring force law. Our derivation of the continuum limit justifies that to reach consistent dynamics as the number of springs increases, the spring restoring force laws must be rescaled appropriately. Despite mechanical relaxation occurring within a curved tissue layer, we find that the curvature of the tissue does not affect tangential stress nor the mechanics-induced redistribution of cells within the layer in the continuum limit. However, the cell's normal stress does depend on curvature due to surface tension induced by the tangential forces. By characterising the full stress state of a cell, these models provide a basis to represent further mechanobiological processes., Comment: 34 pages, 13 figures; v2: minor edits

Generalized additive models (GAMs, Hastie & Tibshirani, 1990; Wood, 2017) are an extension of the generalized linear model that allows the effects of covariates to be modelled as smooth functions. GAMs are increasingly used in many areas of science (e.g. Pedersen, Miller, Simpson, & Ross, 2019; Simpson, 2018) because the smooth functions allow nonlinear relationships between covariates and the response to be learned from the data through the use of penalized splines. Within the R (R Core Team, 2024) ecosystem, Simon Wood's mgcv package (Wood, 2017) is widely used to fit GAMs and is a Recommended package that ships with R as part of the default install. A growing number of other R packages build upon mgcv, for example as an engine to fit specialised models not handled by mgcv itself (e.g. GJMR, Marra & Radice, 2023), or to make use of the wide range of splines available in mgcv (e.g. brms, B\"urkner, 2017). The gratia package builds upon mgcv by providing functions that make working with GAMs easier. gratia takes a tidy approach (Wickham, 2014) providing ggplot2 (Wickham, 2016) replacements for mgcv's base graphics-based plots, functions for model diagnostics and exploration of fitted models, and a family of functions for drawing samples from the posterior distribution of a fitted GAM. Additional functionality is provided to facilitate the teaching and understanding of GAMs., Comment: 9 pages, 4 figures, submitted to Journal of Open Source Software

In vitro cell biology experiments are routinely used to characterize cell migration properties under various experimental conditions. These experiments can be interpreted using lattice-based random walk models to provide insight into underlying biological mechanisms, and continuum limit partial differential equation (PDE) descriptions of the stochastic models can be used to efficiently explore model properties instead of relying on repeated stochastic simulations. Working with efficient PDE models is of high interest for parameter estimation algorithms that typically require a large number of forward model simulations. Quantitative data from cell biology experiments usually involves non-negative cell counts in different regions of the experimental images, and it is not obvious how to relate finite, noisy count data to the solutions of continuous PDE models that correspond to noise-free density profiles. In this work we illustrate how to develop and implement likelihood-based methods for parameter estimation, parameter identifiability and model prediction for lattice-based models describing collective migration with an arbitrary number of interacting subpopulations. We implement a standard additive Gaussian measurement error model as well as a new physically-motivated multinomial measurement error model that relates noisy count data with the solution of continuous PDE models. Both measurement error models lead to similar outcomes for parameter estimation and parameter identifiability, whereas the standard additive Gaussian measurement error model leads to non-physical prediction outcomes. In contrast, the new multinomial measurement error model involves a lower computational overhead for parameter estimation and identifiability analysis, as well as leading to physically meaningful model predictions., Comment: 34 pages, 7 figures

Optical photothermal microscopy is a powerful, emerging method that overcomes the diffraction limit in infrared hyperspectral imaging by utilizing a visible probe laser beam to detect local temperature-induced modulation at the visible diffraction limit. However, the spectral range of this technique has been limited by the tuning range of infrared sources, which is generally restricted to the fingerprint window with commercially available quantum cascade lasers. In this work, ultra-broadband synchrotron infrared radiation was used for infrared photothermal imaging and spectroscopy, spanning the entire mid-infrared range. Both optical- and fluorescence-detected photothermal modalities were performed, demonstrating improved spectral range when compared to optical photothermal microscopy using commercial sources and improved spatial resolution when compared to synchrotron micro-spectroscopy measurements. Following these initial validation studies, synchrotron Fourier-transform optical photothermal infrared spectroscopy (FT-OPTIR) in combination with synchrotron micro-spectroscopy measurements were used to differentiate cells in mouse brain tissue sections.

Chemical functionalization of diamond surfaces by hydrogen is an important method for controlling the charge state of near-surface fluorescent color centers, an essential process in fabricating devices such as diamond field-effect transistors and chemical sensors, and a required first step for realizing families of more complex terminations through subsequent chemical processing. In all these cases, termination is typically achieved using hydrogen plasma sources which can etch or damage the diamond as well as deposited materials or embedded colour centers. This work explores alternative methods for lower-damage hydrogenation of diamond surfaces, specifically the annealing of diamond samples in high-purity, non-explosive mixtures of nitrogen and hydrogen gas, and the exposure of samples to microwave hydrogen plasmas in the absence of intentional stage heating. The effectiveness of these methods are characterized by x-ray photoelectron spectroscopy, and comparison of the results to density-functional modelling of the surface hydrogenation energetics implicates surface oxygen ligands as the primary factor limiting the termination quality of annealed samples. Finally, photoluminescence spectroscopy is used to verify that both the annealing and reduced sample temperature plasma methods are non-destructive to near-surface ensembles of nitrogen-vacancy centers, in stark contrast to plasma treatments which use heated sample stages., Comment: 24 pages, 4 figures

This research explores the use of History to enhance education in the field of ICT For Sustainability ICT4S in response to a challenge from the ICT4S 2023 conference. No previous studies were found in ICT4S but the literature on History and Education for Sustainable Development is reviewed. An ICT4S lecturer collaborated with History lecturers to add an historic parallel to each weeks teaching on a Sustainable Business and Computing unit for final year undergraduate BSc Business Computing students. A list of the topics and rationale is provided. Student perceptions were surveyed before and after the teaching and semi-structured interviews carried out. A majority of students saw relevance to their degree and career. There was an increase in the proportion of students with interest in History. The paper explores the lessons learned from the interdisciplinary collaboration, including topic choice, format and perceived value. The project has enhanced the way we approach our subjects as computing and history educators. We believe this is the first empirical, survey-based study of the use of history to enhance ICT4S education. The team will extend the research to a larger unit covering a wider range of computing degrees.

Magnetic resonance imaging (MRI) is a crucial tool to identify brain abnormalities in a wide range of neurological disorders. In focal epilepsy MRI is used to identify structural cerebral abnormalities. For covert lesions, machine learning and artificial intelligence algorithms may improve lesion detection if abnormalities are not evident on visual inspection. The success of this approach depends on the volume and quality of training data. Herein, we release an open-source dataset of preprocessed MRI scans from 442 individuals with drug-refractory focal epilepsy who had neurosurgical resections, and detailed demographic information. The MRI scan data includes the preoperative 3D T1 and where available 3D FLAIR, as well as a manually inspected complete surface reconstruction and volumetric parcellations. Demographic information includes age, sex, age of onset of epilepsy, location of surgery, histopathology of resected specimen, occurrence and frequency of focal seizures with and without impairment of awareness, focal to bilateral tonic-clonic seizures, number of anti-seizure medications (ASMs) at time of surgery, and a total of 1764 patient years of post-surgical follow up. Crucially, we also include resection masks delineated from post-surgical imaging. To demonstrate the veracity of our data, we successfully replicated previous studies showing long-term outcomes of seizure freedom in the range of around 50%. Our imaging data replicates findings of group level atrophy in patients compared to controls. Resection locations in the cohort were predominantly in the temporal and frontal lobes. We envisage our dataset, shared openly with the community, will catalyse the development and application of computational methods in clinical neurology.

A laser pulse composed of a fundamental and an appropriately phased second harmonic can drive a time-dependent current of photoionized electrons that generates broadband THz radiation. Over the propagation distances relevant to many experiments, dispersion causes the relative phase between the harmonics to evolve. This "dephasing" slows the accumulation of THz energy and results in a multi-cycle THz pulse with significant angular dispersion. Here, we introduce a novel optical configuration that compensates the relative phase evolution, allowing for the formation of a half-cycle THz pulse with almost no angular dispersion. The configuration uses the spherical aberration of an axilens to map a prescribed radial phase variation in the near field to a desired longitudinal phase variation in the far field. Simulations that combine this configuration with an ultrashort flying focus demonstrate the formation of a half-cycle THz pulse with a controlled emission angle and 1/4 the angular divergence of the multi-cycle pulse created by a conventional optical configuration.

A nonsmooth fold is where an equilibrium or limit cycle of a nonsmooth dynamical system hits a switching manifold and collides and annihilates with another solution of the same type. We show that beyond the bifurcation the leading-order truncation to the system in general has no bounded invariant set. This is proved for boundary equilibrium bifurcations of Filippov systems, hybrid systems, and continuous piecewise-smooth ODEs, and grazing-type events for which the truncated form is a continuous piecewise-linear map. The omitted higher-order terms are expected to be incapable of altering the local dynamics qualitatively, implying the system has no local invariant set on one side of a nonsmooth fold, and we demonstrate this with an example. Thus if the equilibrium or limit cycle is attracting the bifurcation causes the local attractor of the system to tip to a new state. The results also help explain global aspects of the bifurcation structures of the truncated systems.