Search

This study explores the potential of large language models (LLMs) for identifying and examining intertextual relationships within biblical, Koine Greek texts. By evaluating the performance of LLMs on various intertextuality scenarios the study demonstrates that these models can detect direct quotations, allusions, and echoes between texts. The LLM's ability to generate novel intertextual observations and connections highlights its potential to uncover new insights. However, the model also struggles with long query passages and the inclusion of false intertextual dependences, emphasizing the importance of expert evaluation. The expert-in-the-loop methodology presented offers a scalable approach for intertextual research into the complex web of intertextuality within and beyond the biblical corpus.

When developing machine learning models, image quality assessment (IQA) measures are a crucial component for evaluation. However, commonly used IQA measures have been primarily developed and optimized for natural images. In many specialized settings, such as medical images, this poses an often-overlooked problem regarding suitability. In previous studies, the IQA measure HaarPSI showed promising behavior for natural and medical images. HaarPSI is based on Haar wavelet representations and the framework allows optimization of two parameters. So far, these parameters have been aligned for natural images. Here, we optimize these parameters for two annotated medical data sets, a photoacoustic and a chest X-Ray data set. We observe that they are more sensitive to the parameter choices than the employed natural images, and on the other hand both medical data sets lead to similar parameter values when optimized. We denote the optimized setting, which improves the performance for the medical images notably, by HaarPSI$_{MED}$. The results suggest that adapting common IQA measures within their frameworks for medical images can provide a valuable, generalizable addition to the employment of more specific task-based measures., Comment: 5 pages, 3 figures, 2 tables

We introduce Image2Struct, a benchmark to evaluate vision-language models (VLMs) on extracting structure from images. Our benchmark 1) captures real-world use cases, 2) is fully automatic and does not require human judgment, and 3) is based on a renewable stream of fresh data. In Image2Struct, VLMs are prompted to generate the underlying structure (e.g., LaTeX code or HTML) from an input image (e.g., webpage screenshot). The structure is then rendered to produce an output image (e.g., rendered webpage), which is compared against the input image to produce a similarity score. This round-trip evaluation allows us to quantitatively evaluate VLMs on tasks with multiple valid structures. We create a pipeline that downloads fresh data from active online communities upon execution and evaluates the VLMs without human intervention. We introduce three domains (Webpages, LaTeX, and Musical Scores) and use five image metrics (pixel similarity, cosine similarity between the Inception vectors, learned perceptual image patch similarity, structural similarity index measure, and earth mover similarity) that allow efficient and automatic comparison between pairs of images. We evaluate Image2Struct on 14 prominent VLMs and find that scores vary widely, indicating that Image2Struct can differentiate between the performances of different VLMs. Additionally, the best score varies considerably across domains (e.g., 0.402 on sheet music vs. 0.830 on LaTeX equations), indicating that Image2Struct contains tasks of varying difficulty. For transparency, we release the full results at https://crfm.stanford.edu/helm/image2struct/v1.0.1/., Comment: NeurIPS 2024. First three authors contributed equally

Synchrotron X-ray emission has been detected from nearly a dozen young supernova remnants (SNRs). X-rays of synchrotron origin exhibit linear polarization in a regular, non-randomly oriented magnetic field. The significant polarized X-ray emission from four such SNRs has already been reported on the basis of observations with the Imaging X-ray Polarimetry Explorer (IXPE). The magnetic-field structure as derived from IXPE observations is radial for Cassiopeia A, Tycho's SNR, and SN 1006, and tangential for RX J1713.7-3946. The latter together with the recent detection of a tangential magnetic field in SNR 1E 0102.2-7219 by the Australia Telescope Compact Array in the radio band shows that tangential magnetic fields can also be present in young SNRs. Thus, the dichotomy in polarization between young and middle-aged SNRs (radial magnetic fields in young SNRs, but tangential magnetic fields in middle-aged SNRs), previously noticed in the radio band, deserves additional attention. The present analysis of IXPE observations determines, for the first time, a magnetic-field structure in the northwestern rim of Vela Jr, also known as RX J0852.0-4622, and provides a new example of a young SNR with a tangential magnetic field., Comment: Accepted for publication in Astronomy and Astrophysics

X-ray polarization is a unique new probe of the particle acceleration in astrophysical jets made possible through the Imaging X-ray Polarimetry Explorer. Here we report on the first dense X-ray polarization monitoring campaign on the blazar Mrk 421. Our observations were accompanied by an even denser radio and optical polarization campaign. We find significant short-timescale variability in both X-ray polarization degree and angle, including a $\sim90^\circ$ angle rotation about the jet axis. We attribute this to random variations of the magnetic field, consistent with the presence of turbulence but also unlikely to be explained by turbulence alone. At the same time, the degree of lower-energy polarization is significantly lower and shows no more than mild variability. Our campaign provides further evidence for a scenario in which energy-stratified shock-acceleration of relativistic electrons, combined with a turbulent magnetic field, is responsible for optical to X-ray synchrotron emission in blazar jets., Comment: 23 pages, including 8 pages of appendices. 12 figures, 3 tables. Submitted to ApJ

We present JWST/NIRSpec PRISM spectroscopic characterization of GHZ9 at z= 10.145 $\pm$ 0.010, currently the most distant source detected by the Chandra X-ray Observatory. The spectrum reveals several UV high-ionization lines, including CII, SiIV, [NIV], CIV, HeII, OIII], NIII], and CIII]. The prominent rest-frame equivalent widths (EW(CIV)$\simeq$65A, EW(HeII)$\simeq$18A, EW(CIII])$\simeq$48A) show the presence of a hard radiation field, while the analysis of line ratio diagnostics suggest this galaxy hosts both AGN and star-formation activity. GHZ9 is nitrogen-enriched (6--9.5 times solar), carbon-poor (0.2--0.65 times solar), metal-poor (Z = 0.01--0.1 Z$_{\odot}$), and compact ($<$ 106 pc), similarly to GNz11, GHZ2, and recently discovered N-enhanced high redshift objects. We exploited the newly available JWST/NIRSpec and NIRCam dataset to perform an independent analysis of the Chandra data confirming that GHZ9 is the most likely JWST source associated to X-ray emission at 0.5-7 keV. Assuming a spectral index $\Gamma$ = 2.3 (1.8), we estimate a black hole (BH) mass of 1.60 $\pm$ 0.31 (0.48 $\pm$ 0.09) $\times$ 10$^8$M$_{\odot}$, which is consistent either with Eddington-accretion onto heavy ($\geq$ 10$^6$ M$_{\odot}$) BH seeds formed at z=18, or super-Eddington accretion onto a light seed of $\sim$ 10$^2-10^4$ M$_{\odot}$ at z = 25. The corresponding BH-to-stellar mass ratio M$_{BH}$/M$_{star}$= 0.33$\pm$0.22 (0.10$\pm$0.07), with a stringent limit $>$0.02, implies an accelerated growth of the BH mass with respect to the stellar mass. GHZ9 is the ideal target to constrain the early phases of AGN-galaxy coevolution with future multi-frequency observations., Comment: Submitted to ApJ

Scanning objects with a more tightly focused beam (for example of photons or electrons) can provide higher-resolution images. However the stronger localisation of energy deposition can damage tissues in organic samples or may rearrange the chemical structure or physical properties of inorganic materials. Scanning an object with a broad beam can deliver an equivalent probe energy but spreads it over a much wider footprint. Sharp images can be reconstructed from the diffuse implanted signal when a decoding step can recover a delta-like impulse response. Huffman sequences, by design, have the optimal delta-like autocorrelation for aperiodic (non-cyclic) convolution and are well-conditioned. Here we adapt 1D Huffman sequences to design 2D Huffman-like discrete arrays that have spatially broad, relatively thin and uniform intensity profiles that retain excellent aperiodic autocorrelation metrics. Examples of broad shaped diffuse beams were developed for the case of x-ray imaging. A variety of masks were fabricated by the deposition of finely structured layers of tantalum on a silicon oxide wafer. The layers form a pattern of discrete pixels that modify the shape of an incident uniform beam of low energy x-rays as it passes through the mask. The intensity profiles of the x-ray beams after transmission through these masks were validated, first by acquiring direct-detector x-ray images of the masks, and second by raster scanning a pinhole over each mask pattern, pixel-by-pixel, collecting "bucket" signals as applied in traditional ghost imaging. The masks were then used to raster scan the shaped x-ray beam over several simple binary and "gray" test objects, again producing bucket signals, from which sharp reconstructed object images were obtained by deconvolving their bucket images.

The US opioid overdose epidemic has been a major public health concern in recent decades. There has been increasing recognition that its etiology is rooted in part in the social contexts that mediate substance use and access; however, reliable statistical measures of social influence are lacking in the literature. We use Facebook's social connectedness index (SCI) as a proxy for real-life social networks across diverse spatial regions that help quantify social connectivity across different spatial units. This is a measure of the relative probability of connections between localities that offers a unique lens to understand the effects of social networks on health outcomes. We use SCI to develop a variable, called "deaths in social proximity", to measure the influence of social networks on opioid overdose deaths (OODs) in US counties. Our results show a statistically significant effect size for deaths in social proximity on OODs in counties in the United States, controlling for spatial proximity, as well as demographic and clinical covariates. The effect size of standardized deaths in social proximity in our cluster-robust linear regression model indicates that a one-standard-deviation increase, equal to 11.70 more deaths per 100,000 population in the social proximity of ego counties in the contiguous United States, is associated with thirteen more deaths per 100,000 population in ego counties. To further validate our findings, we performed a series of robustness checks using a network autocorrelation model to account for social network effects, a spatial autocorrelation model to capture spatial dependencies, and a two-way fixed-effect model to control for unobserved spatial and time-invariant characteristics. These checks consistently provide statistically robust evidence of positive social influence on OODs in US counties.

Evidence for high-redshift supermassive black holes challenges standard scenarios for how such objects form in the early universe. Here, we entertain the possibility that a fraction of the cosmological dark matter could be ultra-strongly self interacting. This would imply that gravothermal collapse occur at early times in the cores of dark matter halos, followed by accretion. We study under which conditions on the abundance and interaction strength and structure of such ultra self-interacting dark matter the black holes resulting from the end-point of gravothermal core collapse can seed the observed, early-forming supermassive black holes. We find, depending on the velocity dependence of the self-interaction cross section, a bimodal structure in the favored parameter space, where data points to either a small collapsing dark matter fraction with a large cross section, or a large fraction and a relatively small cross section. While self-interaction cross sections with different velocity dependence can explain observations, we find that the best, self-consistent results correspond to a Rutherford-like self-interaction, typical of long-range dark-sector forces with light mediators. We discuss complementary observational probes if this scenario is realized in nature, focusing especially on the expected intermediate mass black holes predicted to exist in smaller galaxies., Comment: 28 pages, 6 Figures

Intels efforts to build a practical quantum computer are focused on developing a scalable spin-qubit platform leveraging industrial high-volume semiconductor manufacturing expertise and 300 mm fabrication infrastructure. Here, we provide an overview of the design, fabrication, and demonstration of a new customized quantum test chip, which contains 12-quantum-dot spin-qubit linear arrays, code named Tunnel Falls. These devices are fabricated using immersion and extreme ultraviolet lithography (EUV), along with other standard high-volume manufacturing (HVM) processes, as well as production-level process control. We present key device features and fabrication details, as well as qubit characterization results confirming device functionality. These results corroborate our fabrication methods and are a crucial step towards scaling of extensible 2D qubit array schemes.

A light-front holographic model is used to illustrate an algebraic scheme for constructing a representation of a hadron's zero-skewness generalised parton distribution (GPD) from its valence-quark distribution function (DF) and electromagnetic form factor, $F_H$, without reference to deeply virtual Compton scattering data. The hadron's mass distribution gravitational form factor, $A_H$, calculated from this GPD is harder than $F_H$; and, for each hadron, the associated mass-density profile is more compact than the analogous charge profile, with each pion near-core density being larger than that of its proton partner. These features are independent of the scheme employed., Comment: 7 pages, 6 figures

Stellar active regions like spots and faculae can distort the shapes of spectral lines, inducing variations in the radial velocities that are often orders of magnitude larger than the signals from Earth-like planets. Efforts to mitigate these activity signals have hitherto focused on either the time or the velocity (wavelength) domains. We present a physics-driven Gaussian process (GP) framework to model activity signals directly in time series of line profiles or Cross-Correlation Functions (CCFs). Unlike existing methods which correct activity signals in line profile time series, our approach exploits the time correlation between velocity (wavelength) bins in the line profile variations, and is based on a simplified but physically motivated model for the origin of these variations. When tested on both synthetic and real data sets with signal-to-noise ratios down to $\sim$ 100, our method was able to separate the planetary signal from the activity signal, even when their periods were identical. We also conducted injection/recovery tests using two years of realistically sampled HARPS-N solar data, demonstrating the ability of the method to accurately recover a signal induced by a 1.5-Earth mass planet with a semi-amplitude of 0.3 m/s and a period of 33 days during high solar activity., Comment: Accepted for publication in MNRAS

Recently there has been a proliferation of intrinsic motivation (IM) reward-shaping methods to learn in complex and sparse-reward environments. These methods can often inadvertently change the set of optimal policies in an environment, leading to suboptimal behavior. Previous work on mitigating the risks of reward shaping, particularly through potential-based reward shaping (PBRS), has not been applicable to many IM methods, as they are often complex, trainable functions themselves, and therefore dependent on a wider set of variables than the traditional reward functions that PBRS was developed for. We present an extension to PBRS that we prove preserves the set of optimal policies under a more general set of functions than has been previously proven. We also present {\em Potential-Based Intrinsic Motivation} (PBIM) and {\em Generalized Reward Matching} (GRM), methods for converting IM rewards into a potential-based form that are useable without altering the set of optimal policies. Testing in the MiniGrid DoorKey and Cliff Walking environments, we demonstrate that PBIM and GRM successfully prevent the agent from converging to a suboptimal policy and can speed up training. Additionally, we prove that GRM is sufficiently general as to encompass all potential-based reward shaping functions. This paper expands on previous work introducing the PBIM method, and provides an extension to the more general method of GRM, as well as additional proofs, experimental results, and discussion., Comment: To be submit to joint AIJ-JAIR special track for award-winning papers. arXiv admin note: substantial text overlap with arXiv:2402.07411

We present JWST/NIRSpec PRISM follow-up of candidate galaxies at z=9-11 selected from deep JWST/NIRCam photometry in GLASS-JWST Early Release Science data. We spectroscopically confirm six sources with secure redshifts at z = 9.52-10.43, each showing multiple emission lines. An additional object is likely at z = 10.66, based on its Lya-break and a single emission feature, while one source is a lower redshift interloper. The sample includes the first JWST-detected candidate at z=10, GHZ1/GLASS-z10, which we confirm at z = 9.875, and the X-ray detected AGN GHZ9 confirmed at z = 10.145. Three objects in our sample, including GHZ9, have EW(CIII])>20A and occupy a region compatible with AGN emission in the EW(CIII]) vs CIV/CIII] diagram. The spectroscopic sample confirms a high abundance of galaxies at z > 9. We measure a number density of z=10 galaxies in the GLASS-JWST ERS field that is a factor of >3 higher than other JWST-based estimates at demagnified rest-frame magnitudes of -21 < Muv < -19. We find that the positions of these galaxies in redshift and angular space are not consistent with all of them being part of a unique progenitor of present-day galaxy clusters. The high density of objects in the GLASS region can be explained either by clustering on large scales or by a superposition of different forming structures of which we observe only the brightest members. By considering all the spectroscopic z=10 sources in the Abell-2744 field, we identify two potential galaxy proto-clusters centered around GHZ9 and JD1, with relative separations between their members of 1-2 pMpc. The potential AGN nature of three of the sources in our sample lends support to a scenario in which the high abundance of bright sources determined by JWST surveys at cosmic dawn may be affected by AGN contribution to their UV luminosity., Comment: Submitted to A&A

A version control system, such as Git, requires a way to integrate changes from different developers or branches. Given a merge scenario, a merge tool either outputs a clean integration of the changes, or it outputs a conflict for manual resolution. A clean integration is correct if it preserves intended program behavior, and is incorrect otherwise (e.g., if it causes a test failure). Manual resolution consumes valuable developer time, and correcting a defect introduced by an incorrect merge is even more costly. New merge tools have been proposed, but they have not yet been evaluated against one another. Prior evaluations do not properly distinguish between correct and incorrect merges, are not evaluated on a realistic set of merge scenarios, and/or do not compare to state-of-the-art tools. We have performed a more realistic evaluation. The results differ significantly from previous claims, setting the record straight and enabling better future research. Our novel experimental methodology combines running test suites, examining merges on deleted branches, and accounting for the cost of incorrect merges. Based on these evaluations, we created a merge tool that out-performs all previous tools under most assumptions. It handles the most common merge scenarios in practice., Comment: ASE 2024

We introduce the open source code PHOEBUS (phifty one ergs blows up a star) for astrophysical general relativistic radiation magnetohydrodynamic simulations. PHOEBUS is designed for, but not limited to, high energy astrophysical environments such as core-collapse supernovae, neutron star mergers, black-hole accretion disks, and similar phenomena. General relativistic magnetohydrodynamics are modeled in the Valencia formulation with conservative finite volume methods. Neutrino radiation transport is included with Monte Carlo and moment methods. PHOEBUS is built on the PARTHENON (Grete et al. 2022) performance portable adaptive mesh refinement framework, uses a GPU first development strategy, and is capable of modeling a large dynamic range in space and time. PHOEBUS utilizes KOKKOS for on-node parallelism and supports both CPU and GPU architectures. We describe the physical model employed in PHOEBUS, the numerical methods used, and demonstrate a suite of test problems to showcase its abilities. We demonstrate weak scaling to over 500 H100 GPUs., Comment: 21 pages, 15 figures. Code available online at https://github.com/lanl/phoebus

We present constraints on the shape of Kepler-51d, which is a super-puff with a mass $\sim6\,M_\oplus$ and a radius $\sim9\,R_\oplus$, based on detailed modeling of the transit light curve from JWST NIRSpec. The projected shape of this extremely low-density planet is consistent with being spherical, and a projected oblateness $f_\perp>0.2$ can be excluded regardless of the spin obliquity angles. If this is taken as the limit on the true shape of the planet, Kepler-51d is rotating at $\lesssim 50\%$ of its break-up spin rate, or its rotation period is $\gtrsim 33\,$hr. In the more plausible situation that the planetary spin is aligned with its orbital direction to within $30^\circ$, then its oblateness is $<0.08$, which corresponds to a dimensionless spin rate $\lesssim30\%$ of the break-up rotation and a dimensional rotation period $\gtrsim 53\,$hr. This seems to contradict the theoretical expectation that planets with such low masses may be spinning near break-up. We point out the usefulness of the stellar mean density and the orbital eccentricity in constraining the shape of the transiting planet, so planets with well-characterized host and orbital parameters are preferred in the detection of planetary oblateness with the JWST transit method., Comment: 11 pages, 4 figures. Submitted to ApJL

One of the most challenging problems in the computational study of localization in quantum manybody systems is to capture the effects of rare events, which requires sampling over exponentially many disorder realizations. We implement an efficient procedure on a quantum processor, leveraging quantum parallelism, to efficiently sample over all disorder realizations. We observe localization without disorder in quantum many-body dynamics in one and two dimensions: perturbations do not diffuse even though both the generator of evolution and the initial states are fully translationally invariant. The disorder strength as well as its density can be readily tuned using the initial state. Furthermore, we demonstrate the versatility of our platform by measuring Renyi entropies. Our method could also be extended to higher moments of the physical observables and disorder learning.

Current benchmarks for assessing vision-language models (VLMs) often focus on their perception or problem-solving capabilities and neglect other critical aspects such as fairness, multilinguality, or toxicity. Furthermore, they differ in their evaluation procedures and the scope of the evaluation, making it difficult to compare models. To address these issues, we extend the HELM framework to VLMs to present the Holistic Evaluation of Vision Language Models (VHELM). VHELM aggregates various datasets to cover one or more of the 9 aspects: visual perception, knowledge, reasoning, bias, fairness, multilinguality, robustness, toxicity, and safety. In doing so, we produce a comprehensive, multi-dimensional view of the capabilities of the VLMs across these important factors. In addition, we standardize the standard inference parameters, methods of prompting, and evaluation metrics to enable fair comparisons across models. Our framework is designed to be lightweight and automatic so that evaluation runs are cheap and fast. Our initial run evaluates 22 VLMs on 21 existing datasets to provide a holistic snapshot of the models. We uncover new key findings, such as the fact that efficiency-focused models (e.g., Claude 3 Haiku or Gemini 1.5 Flash) perform significantly worse than their full models (e.g., Claude 3 Opus or Gemini 1.5 Pro) on the bias benchmark but not when evaluated on the other aspects. For transparency, we release the raw model generations and complete results on our website (https://crfm.stanford.edu/helm/vhelm/v2.0.1). VHELM is intended to be a living benchmark, and we hope to continue adding new datasets and models over time., Comment: NeurIPS 2024. First three authors contributed equally

We introduce a new capability of the Neil Gehrels Swift Observatory, dubbed `continuous commanding,' achieving 10 seconds latency response time on-orbit to unscheduled Target of Opportunity requests. This allows Swift to respond to early warning gravitational-wave detections, rapidly slewing the Burst Alert Telescope (BAT) across the sky to place the GW origin in the BAT field of view at merger time. This will dramatically increase the GW/GRB co-detection rate, and enable prompt arcminute localization of a neutron star merger. We simulate the full Swift response to a GW early warning alert, including input sky maps produced at different warning times, a complete model of Swift's attitude control system, and a full accounting of the latency between the GW detectors and the spacecraft. 60 s of early warning doubles the rate of prompt GRB detections with arcminute position, and 140 s guarantees observation anywhere on the unocculted sky, even with localization areas >> 1000 deg$^2$. While 140 s is beyond current gravitational wave detector sensitivities, 30-70 s is achievable today. We show that the detection yield is now limited by the latency of LIGO/Virgo cyber-infrastructure, and motivate focus on its reduction. Continuous commanding is now a general capability of Swift, significantly increasing its versatility in response to the growing demands of time-domain astrophysics. We demonstrate this potential on an externally triggered Fast Radio Burst, slewing 81 degrees across the sky, and collecting X-ray and UV photons from the source position < 150 s after the trigger was received from the Canadian Hydrogen Intensity Mapping Experiment (CHIME), thereby setting the earliest and deepest such constraints on high energy activity from non-repeating FRBs. The Swift Team invites proposals for novel scientific applications of ultra-low latency UV, X-ray, and gamma-ray observations., Comment: Update to version accepted for publication in ApJL. Comments, questions, and new ideas for applications of ultra-low latency UV, X-ray, and gamma-ray observations are welcome!

Long-range ensemble forecasts are typically verified as anomalies with respect to a lead-time dependent climatology to remove the influence of systematic biases. However, common methods for calculating anomalies result in statistical inconsistencies between forecast and verification anomalies, even for a perfectly reliable ensemble. It is important to account for these systematic effects when evaluating ensemble forecast systems, particularly when tuning a model to improve the reliability of forecast anomalies or when comparing spread-error diagnostics between systems with different reforecast periods. Here, we show that unbiased variances and spread-error ratios can be recovered by deriving estimators that are consistent with the values that would be achieved when calculating anomalies relative to the true, but unknown, climatological mean. An elegant alternative is to construct forecast anomalies separately for each member, which ensures that forecast and verification anomalies are defined relative to reference climatologies with the same sampling uncertainty. This alternative approach has no impact on forecast ensemble means or deterministic scores but systematically modifies the total variance and ensemble spread of forecast anomalies in such a way that anomaly-based spread-error ratios are unbiased without any explicit correction for climatology sample size. Furthermore, the improved statistical consistency of forecast and verification anomalies means that probabilistic forecast skill is optimized when the underlying forecast is also perfectly reliable. Alternative methods for anomaly calculation can thus impact probabilistic forecast skill, especially when anomalies are defined relative to climatologies with a small sample size. Finally, we demonstrate the equivalence of anomalies calculated using different methods after applying an unbiased statistical calibration.

Gaia astrometry of nearby stars is precise enough to detect the tiny displacements induced by substellar companions, but radial velocity data are needed for definitive confirmation. Here we present radial velocity follow-up observations of 28 M and K stars with candidate astrometric substellar companions, which led to the confirmation of two systems, Gaia-4b and Gaia-5b, and the refutation of 21 systems as stellar binaries. Gaia-4b is a massive planet ($M = 11.8 \pm 0.7 \:\mathrm{M_J}$) in a $P = 571.3 \pm 1.4\:\mathrm{day}$ orbit with a projected semi-major axis $a_0=0.312 \pm 0.040\:\mathrm{mas}$ orbiting a $0.644 \pm 0.02 \:\mathrm{M_\odot}$ star. Gaia-5b is a brown dwarf ($M = 20.9 \pm 0.5\:\mathrm{M_J}$) in a $P = 358.58 \pm 0.19\:\mathrm{days}$ eccentric $e=0.6412 \pm 0.0027$ orbit with a projected angular semi-major axis of $a_0 = 0.947 \pm 0.038\:\mathrm{mas}$ around a $0.34 \pm 0.03 \mathrm{M_\odot}$ star. Gaia-4b is one of the first exoplanets discovered via the astrometric technique, and is one of the most massive planets known to orbit a low-mass star., Comment: Submitted to AAS journals. 26 pages, 12 figures, 4 tables

We present a novel class of cloaking in which a region of space is concealed from an ensemble of diffusing particles whose individual trajectories are governed by a stochastic (Langevin) equation. In particular, we simulate how different interpretations of the Langevin equation affect the cloaking performance of an annular single-layer invisibility cloak of smoothly varying diffusivity in two dimensions. Near-perfect cloaking is achieved under the Ito convention, indicated by the cloak preventing particles from accessing an inner core while simultaneously preserving the particle density outside the cloak relative to simulations involving no protected region (and no cloak). Even better cloaking performance can be achieved by regularising the singular behaviour of the cloak -- which we demonstrate through two different approaches. These results establish the foundations of ``stochastic cloaking'', which we believe to be a significant milestone following that of optical and thermal cloaking., Comment: 12 pages (6 pages main), 4 figures (2 figures main)

Parton distribution functions (DFs) have long been recognised as key measures of hadron structure. Today, theoretical prediction of such DFs is becoming possible using diverse methods for continuum and lattice analyses of strong interaction (QCD) matrix elements. Recent developments include a demonstration that continuum and lattice analyses yield mutually consistent results for all pion DFs, with behaviour on the far valence domain of light-front momentum fraction that matches QCD expectations. Theory is also delivering an understanding of the distributions of proton mass and spin amongst its constituents, which varies, of course, with the resolving scale of the measuring probe. Aspects of these developments are sketched herein., Comment: 9 pages, 3 figures, 1 table. Review/perspective commissioned by Chin. Phys. Lett

In this paper we use level set-based topology optimisation to design three-dimensional periodic piezoelectric materials with enhanced properties. Our methodology is fully memory-distributed and written in Julia using the package GridapTopOpt. We compare and assess several existing iterative solvers with respect to their weak scalability and find that an approximate Schur complement preconditioned GMRES method demonstrates the best performance and scalability for solving the piezoelectric homogenisation equations. We use the developed techniques to computationally design high-resolution piezoelectric metamaterials with enhanced stiffness and piezoelectric properties that yield new insights into material design for sensor, hydrophone, and actuator applications. We suggest two robust structures with simple geometric features that exhibit enhanced piezoelectric properties several times larger than those of the base material. We find that level set-based topology optimisation is well suited to problems involving piezoelectricity and has the advantage of avoiding large regions of intermediate density material.

The Perseus cluster is the brightest X-ray cluster in the sky and is known as a cool-core galaxy cluster. Being a very nearby cluster, it has been extensively studied. This has provided a comprehensive view of the physical processes that operate in the intracluster medium (ICM), including feedback from the AGN 3C84 and measurements of ICM turbulence. Additionally, the Perseus cluster contains a central radio mini-halo. This diffuse radio source traces cosmic ray electrons (re-)accelerated in-situ in the ICM. Here we report on LOFAR high-band antenna 120-168 MHz observations of the Perseus cluster that probe a range of four orders of magnitude in angular scales. In our 0.3 arcsec resolution image, we find that the northern extension of the 3C84 lobe consists of several narrow 1.5-3 kpc parallel strands of emission. In addition, we detect steep-spectrum filaments associated with a previous outburst of the central AGN radio emission filling two known X-ray ghost cavities. At 7 arcsec resolution, our images show a complex structured radio mini-halo, with several edges and filaments. At resolutions of 26 arcsec and 80 arcsec, we discover diffuse radio emission with a 1.1 Mpc extent. We classify this emission as a giant radio halo and its properties are distinct from the inner mini-halo. We also detect two diffuse sources at projected cluster centric radii of 0.7 and 1.0 Mpc. Finally, we observe a 0.9 Mpc long trail of radio emission from the cluster member galaxy IC310, connecting it with the giant radio halo. Together with other recent studies of relaxed clusters, our LOFAR observations indicate that cluster-wide radio emission could be (more) common in cool-core clusters. In the case of the Perseus cluster, a past off-axis merger event that preserved the cool core might have generated enough turbulence to produce an extended radio halo observable at low frequencies., Comment: accepted for publication in A&A, 18 pages, 19 figures

We apply spatial dynamical-systems techniques to prove that certain spatiotemporal patterns in reversible reaction-diffusion equations undergo snaking bifurcations. That is, in a narrow region of parameter space, countably many branches of patterned states coexist that connect at towers of saddle-node bifurcations. Our patterns of interest are contact defects, which are 1-dimensional time-periodic patterns with a spatially oscillating core region that at large distances from the origin in space resemble pure temporally oscillatory states and arise as natural analogues of spiral and target waves in one spatial dimension. We show that these solutions lie on snaking branches that have a more complex structure than has been seen in other contexts. In particular, we predict the existence of families of asymmetric travelling defect solutions with arbitrary background phase offsets, in addition to symmetric standing target and spiral patterns. We prove the presence of these additional patterns by reconciling results in classic ODE studies with results from the spatial-dynamics study of patterns in PDEs and using geometrical information contained in the stable and unstable manifolds of the background wavetrains and their natural equivariance structure., Comment: 32 pages, 11 figures

Kepler-51 is a $\lesssim 1\,\mathrm{Gyr}$-old Sun-like star hosting three transiting planets with radii $\approx 6$-$9\,R_\oplus$ and orbital periods $\approx 45$-$130\,\mathrm{days}$. Transit timing variations (TTVs) measured with past Kepler and Hubble Space Telescope (HST) observations have been successfully modeled by considering gravitational interactions between the three transiting planets, yielding low masses and low mean densities ($\lesssim 0.1\,\mathrm{g/cm^3}$) for all three planets. However, the transit time of the outermost transiting planet Kepler-51d recently measured by the James Webb Space Telescope (JWST) 10 years after the Kepler observations is significantly discrepant from the prediction made by the three-planet TTV model, which we confirmed with ground-based and follow-up HST observations. We show that the departure from the three-planet model is explained by including a fourth outer planet, Kepler-51e, in the TTV model. A wide range of masses ($\lesssim M_\mathrm{Jup}$) and orbital periods ($\lesssim 10\,\mathrm{yr}$) are possible for Kepler-51e. Nevertheless, all the coplanar solutions found from our brute-force search imply masses $\lesssim 10\,M_\oplus$ for the inner transiting planets. Thus their densities remain low, though with larger uncertainties than previously estimated. Unlike other possible solutions, the one in which Kepler-51e is around the $2:1$ mean motion resonance with Kepler-51d implies low orbital eccentricities ($\lesssim 0.05$) and comparable masses ($\sim 5\,M_\oplus$) for all four planets, as is seen in other compact multi-planet systems. This work demonstrates the importance of long-term follow-up of TTV systems for probing longer period planets in a system., Comment: 48 pages, 26 figures, accepted for publication in AJ

The prompt spectra of gamma-ray bursts are known to follow broadband continuum behavior over decades in energy. GRB 221009A, given the moniker the brightest of all time (BOAT), is the brightest gamma-ray burst identified in half a century of observations, and was first identified by the Fermi Gamma-ray Burst Monitor (GBM). On behalf of the Fermi-GBM Team, Lesage et al. (2023) described the initial GBM analysis. Ravasio et al. (2024) report the identification of a spectral line in part of the prompt emission of this burst, which they describe as evolving over 80 s from $\sim$12 MeV to 6 MeV. We report a GBM Team analysis on the Ravasio Line: 1) We cannot identify an instrumental effect that could have produced this signal, and 2) our method of calculating the statistical significance of the line shows it easily exceeds the 5$\sigma$ discovery threshold. We additionally comment on the claim of the line beginning at earlier time intervals, up to 37 MeV, as reported in Zhang et al. (2024). We find that it is reasonable to utilize these measurements for characterization of the line evolution, with caution. We encourage theoretical studies exploring this newly discovered gamma-ray burst spectral feature, unless any rigorous alternative explanation unrelated to the emission from GRB 221009A is identified.

We model emissivities of the HCN and CO $J=1-0$ transitions using measured properties of clouds found in normal star forming galaxies and more extreme systems. These models are compared with observations of HCN and CO $J=1-0$ transitions. We combine these model emissivities with predictions of gravoturbulent models of star formation, explore the impact of excitation and optical depth on CO and HCN emission, and assess if observed HCN/CO ratios track the fraction of gravitationally-bound dense gas, $f_\mathrm{grav}$, in molecular clouds. Our modeled HCN/CO ratios and emissivities are consistent with measurements from observations. CO emission shows a range of optical depths across different environments, from optically thick in normal galaxies to moderately optically thin in extreme systems. HCN is only moderately optically thick, with significant subthermal excitation in both normal and extreme galaxies. We find an anticorrelation between HCN/CO and $f_\mathrm{grav}$ as predicted by gravoturbulent models of star formation. Instead this ratio tracks gas at moderate densities ($n>10^{3.5}\ \mathrm{cm}^{-3}$), which is below the standard dense gas threshold of $n>10^{4.5}\ \mathrm{cm}^{-3}$. Variations in CO emissivity depend strongly on optical depth, due to variations in the dynamics of the cloud gas. HCN emissivity depends more strongly on excitation, and thus does not directly track variations in CO emissivity. We conclude that a single line ratio, such as HCN/CO, will not consistently track the fraction of gravitationally-bound, star-forming gas if the critical density for star formation varies in molecular clouds. This work highlights important uncertainties that need to be considered when observationally applying an HCN conversion factor in order to estimate the dense (i.e. $n>10^{4.5}\ \mathrm{cm}^{-3}$) gas content in nearby galaxies., Comment: Accepted for publication in A&A

This paper presents a novel approach to multi-robot planning and collaboration. We demonstrate a cognitive strategy for robots in human-robot teams that incorporates metacognition, natural language communication, and explainability. The system is embodied using the HARMONIC architecture that flexibly integrates cognitive and control capabilities across the team. We evaluate our approach through simulation experiments involving a joint search task by a team of heterogeneous robots (a UGV and a drone) and a human. We detail the system's handling of complex, real-world scenarios, effective action coordination between robots with different capabilities, and natural human-robot communication. This work demonstrates that the robots' ability to reason about plans, goals, and attitudes, and to provide explanations for actions and decisions are essential prerequisites for realistic human-robot teaming., Comment: Submitted to ICRA 2025 Conference, Atlanta, GA, USA

We present HARMONIC, a framework for implementing cognitive robots that transforms general-purpose robots into trusted teammates capable of complex decision-making, natural communication and human-level explanation. The framework supports interoperability between a strategic (cognitive) layer for high-level decision-making and a tactical (robot) layer for low-level control and execution. We describe the core features of the framework and our initial implementation, in which HARMONIC was deployed on a simulated UGV and drone involved in a multi-robot search and retrieval task., Comment: Accepted for presentation at ICRA@40. 23-26 September 2024, Rotterdam, Netherlands

We perform a theoretical study of vacuum polarization corrections to the hyperfine structure in many-electron atoms. Calculations are performed for systems of interest for precision atomic tests of fundamental physics belonging to the alkali-metal atoms and singly-ionized alkaline earths. The vacuum polarization is considered in the Uehling approximation, and we study the many-body effects core relaxation, core polarization, and valence-core correlations in the relativistic framework. We find that for s states, the relative vacuum polarization correction may be well-approximated by that for hydrogen-like ions, though for all other states account of many-body effects -- in particular, the polarization of the core -- is needed to obtain the correct sign and magnitude of the effect.

Lattice gauge theories (LGTs) can be employed to understand a wide range of phenomena, from elementary particle scattering in high-energy physics to effective descriptions of many-body interactions in materials. Studying dynamical properties of emergent phases can be challenging as it requires solving many-body problems that are generally beyond perturbative limits. We investigate the dynamics of local excitations in a $\mathbb{Z}_2$ LGT using a two-dimensional lattice of superconducting qubits. We first construct a simple variational circuit which prepares low-energy states that have a large overlap with the ground state; then we create particles with local gates and simulate their quantum dynamics via a discretized time evolution. As the effective magnetic field is increased, our measurements show signatures of transitioning from deconfined to confined dynamics. For confined excitations, the magnetic field induces a tension in the string connecting them. Our method allows us to experimentally image string dynamics in a (2+1)D LGT from which we uncover two distinct regimes inside the confining phase: for weak confinement the string fluctuates strongly in the transverse direction, while for strong confinement transverse fluctuations are effectively frozen. In addition, we demonstrate a resonance condition at which dynamical string breaking is facilitated. Our LGT implementation on a quantum processor presents a novel set of techniques for investigating emergent particle and string dynamics.

We report on the discovery of a transiting giant planet around the 3500 K M3-dwarf star TOI-6383A located 172 pc from Earth. It was detected by the Transiting Exoplanet Survey Satellite (TESS) and confirmed by a combination of ground-based follow-up photometry and precise radial velocity measurements. This planet has an orbital period of $\sim$1.791 days, mass of 1.040$\pm$0.094 $M_J$ and a radius of 1d.008$^{+0.036}_{-0.033} ~R_J$, resulting in a mean bulk density of 1.26$^{+0.18}_{-0.17}$ g cm$^{-3}$. TOI-6383A has an M-dwarf companion star, TOI-6383B, which has a stellar effective temperature $T_{eff}$ $\sim$ 3100 K and a projected orbital separation of 3100 AU. TOI-6383A is a low-mass dwarf star hosting a giant planet and is an intriguing object for planetary evolution studies due to its high planet-to-star mass ratio. This discovery is part of the \textit{Searching for Giant Exoplanets around M-dwarf Stars (GEMS)} Survey, intending to provide robust and accurate estimates of the occurrence of GEMS and the statistics on their physical and orbital parameters. This paper presents an interesting addition to the small number of confirmed GEMS, particularly notable since its formation necessitates massive, ust-rich protoplanetary discs and high accretion efficiency ($>$ 10\%)., Comment: 20 pages, 8 figures

Though parameter shift rules have drastically improved gradient estimation methods for several types of quantum circuits, leading to improved performance in downstream tasks, so far they have not been transferable to linear optics with single photons. In this work, we derive an analytical formula for the gradients in these circuits with respect to phaseshifters via a generalized parameter shift rule, where the number of parameter shifts depends linearly on the total number of photons. Experimentally, this enables access to derivatives in photonic systems without the need for finite difference approximations. Building on this, we propose two strategies through which one can reduce the number of shifts in the expression, and hence reduce the overall sample complexity. Numerically, we show that this generalized parameter-shift rule can converge to the minimum of a cost function with fewer parameter update steps than alternative techniques. We anticipate that this method will open up new avenues to solving optimization problems with photonic systems, as well as provide new techniques for the experimental characterization and control of linear optical systems., Comment: 13 + 11 pages, 7 figures

Some X-ray binaries containing an energetic pulsar in orbit around a normal star accelerate particles to high energies in the shock cone formed where the pulsar and stellar winds collide. The magnetic field geometry in the acceleration region in such binaries is unknown. We performed the first measurement of the polarization of the X-ray synchrotron emission from a gamma-ray emitting binary system. We observed PSR B1259-63 with the Imaging X-ray Polarimetry Explorer (IXPE) during an X-ray bright phase following the periastron passage in June 2024. X-ray polarization is detected with a polarization degree of $8.3\% \pm 1.5\%$ at a significance of $5.3 \sigma$. The X-ray polarization angle is aligned with the axis of the shock cone at the time of the observation. This indicates that the predominant component of the magnetic field in the acceleration region is oriented perpendicular to the shock cone axis., Comment: 5 pages, to appear in the Astrophysical Journal Letters

A symmetry-preserving analysis of strong interaction quantum field equations is used to complete a unified treatment of pion, kaon, nucleon electromagnetic and gravitational form factors. Findings include a demonstration that the pion near-core pressure is roughly twice that in the proton, so both are significantly greater than that of a neutron star; parton species separations of the nucleon's three gravitational form factors, in which, inter alia, the glue-to-quark ratio for each form factor is seen to take the same constant value, independent of momentum transfer; and a determination of proton radii orderings, with the mechanical (normal force) radius being less than the mass-energy radius, which is less than the proton charge radius. This body of predictions should prove useful in an era of anticipated experiments that will enable them to be tested., Comment: 11 pages, 8 figures, 1 table

Creating a physiologically relevant shear stress in organ-on-a-chip (OOC) devices requires careful tailoring of microfluidic flow parameters. It is currently fairly common to use a simple approximation assuming a constant viscosity, even for serum-based media. Here, we show that a popular nutrient solution (Dulbecco's Modified Eagle Medium supplemented with Fetal Bovine Serum) requires a more complex treatment (i.e., is a non-Newtonian fluid), with observed shear stress values significantly greater than reported in literature. We measure the rheology of the solutions and combine it with a 3-dimensional flow field measurement to derive shear stress at the channel surface. We verify the experiments with numerical simulations, finding good agreement and deriving flow properties. Finally, we provide relevant expressions for shear stress approximation, suitable for development of OOC devices with various geometries.

We present a series of five new broadband X-ray observations of the ultraluminous X-ray source Holmberg IX X-1, performed by $XMM$-$Newton$ and $NuSTAR$ in coordination. The first three of these show high soft X-ray fluxes but a near total collapse of the high-energy ($\gtrsim$15 keV) emission, previously seen to be surprisingly stable across all prior broadband observations of the source. The latter two show a recovery in hard X-rays, remarkably once again respecting the same stable high-energy flux exhibited by all of the archival observations. We also present a joint analysis of all broadband observations of Holmberg IX X-1 to date (encompassing 11 epochs in total) in order to investigate whether it shows the same luminosity-temperature behaviour as NGC 1313 X-1 (which also shows a stable high-energy flux), whereby the hotter disc component in the spectrum exhibits two distinct, positively-correlated tracks in the luminosity-temperature plane. Holmberg IX X-1 may show similar behaviour, but the results depend on whether the highest energy emission is assumed to be an up-scattering corona or an accretion column. The strongest evidence for this behaviour is found in the former case, while in the latter the new 'soft' epochs appear distinct from the other high-flux epochs. We discuss possible explanations for these new 'soft' spectra in the context of the expected structure of super-Eddington accretion flows around black holes and neutron stars, and highlight a potentially interesting analogy with the recent destruction and re-creation of the corona seen in the AGN 1ES 1927+654., Comment: 15 pages, 5 figures, submitted for publication in MNRAS (comments welcome)

Tables, figures, and listings (TFLs) are essential tools for summarizing clinical trial data. Creation of TFLs for reporting activities is often a time-consuming task encountered routinely during the execution of clinical trials. This study explored the use of large language models (LLMs) to automate the generation of TFLs through prompt engineering and few-shot transfer learning. Using public clinical trial data in ADaM format, our results demonstrated that LLMs can efficiently generate TFLs with prompt instructions, showcasing their potential in this domain. Furthermore, we developed a conservational agent named Clinical Trial TFL Generation Agent: An app that matches user queries to predefined prompts that produce customized programs to generate specific predefined TFLs.

Many commonly studied species now have more than one chromosome-scale genome assembly, revealing a large amount of genetic diversity previously missed by approaches that map short reads to a single reference. However, many species still lack multiple reference genomes and correctly aligning references to build pangenomes is challenging, limiting our ability to study this missing genomic variation in population genetics. Here, we argue that $k$-mers are a crucial stepping stone to bridging the reference-focused paradigms of population genetics with the reference-free paradigms of pangenomics. We review current literature on the uses of $k$-mers for performing three core components of most population genetics analyses: identifying, measuring, and explaining patterns of genetic variation. We also demonstrate how different $k$-mer-based measures of genetic variation behave in population genetic simulations according to the choice of $k$, depth of sequencing coverage, and degree of data compression. Overall, we find that $k$-mer-based measures of genetic diversity scale consistently with pairwise nucleotide diversity ($\pi$) up to values of about $\pi = 0.025$ ($R^2 = 0.97$) for neutrally evolving populations. For populations with even more variation, using shorter $k$-mers will maintain the scalability up to at least $\pi = 0.1$. Furthermore, in our simulated populations, $k$-mer dissimilarity values can be reliably approximated from counting bloom filters, highlighting a potential avenue to decreasing the memory burden of $k$-mer based genomic dissimilarity analyses. For future studies, there is a great opportunity to further develop methods to identifying selected loci using $k$-mers., Comment: 6 figures, 34 pages

A symmetry preserving treatment of a vector $\otimes$ vector contact interaction (SCI) is used as the basis for calculations of the two pion transverse momentum dependent parton distribution functions (TMDs); namely, that for unpolarised valence degrees-of-freedom and the analogous Boer-Mulders (BM) function. Amongst other things, the analysis enables the following themes to be addressed: the quark current mass dependence of pion TMDs; the impact of the gauge link model on the positivity constraint that bounds the BM function relative to the unpolarised TMD; the equivalence of direct diagrammatic and light-front wave function TMD calculations; and the size of the BM shift. Interpreted astutely, these SCI results enable one to draw insightful pictures of pion TMDs., Comment: 14 pages, 9 figures, 1 table

Over the Nullarbor Plain in South Australia, the Desert Fireball Network detected a fireball on the night of 1 June 2019 (7:30 pm local time), and six weeks later recovered a single meteorite (42 g) named Arpu Kuilpu. This meteorite was then distributed to a consortium of collaborating institutions to be measured and analyzed by a number of methodologies including: SEM-EDS, EPMA, ICP-MS, gamma-ray spectrometry, ideal gas pycnometry, magnetic susceptibility measurement, {\mu}CT, optical microscopy, and accelerator and noble gas mass spectrometry techniques. These analyses revealed that Arpu Kuilpu is an unbrecciated H5 ordinary chondrite, with minimal weathering (W0-1) and minimal shock (S2). The olivine and pyroxene mineral compositions (in mol%) are Fa: 19.2 +- 0.2, and Fs: 16.8 +- 0.2, further supporting the H5 type and class. The measured oxygen isotopes are also consistent with an H chondrite ({\delta}17O = 2.904 +- 0.177; {\delta}18O = 4.163 +- 0.336; {\Delta}17O = 0.740 +- 0.002). Ideal gas pycnometry measured bulk and grain densities of 3.66 +- 0.02 and 3.77 +- 0.02 g cm-3, respectively, yielding a porosity of 3.0 % +- 0.7. The magnetic susceptibility of this meteorite is log X = 5.16 +- 0.08. The most recent impact-related heating event experienced by Arpu Kuilpu was measured by 40Ar/39Ar chronology to be 4467 +- 16 Ma, while the cosmic ray exposure age is estimated to be between 6-8 Ma. The noble gas isotopes, radionuclides, and fireball observations all indicate that Arpu Kuilpu's meteoroid was quite small (maximum radius of 10 cm, though more likely between 1-5 cm). Although this meteorite is a rather ordinary ordinary chondrite, its prior orbit resembled that of a Jupiter Family Comet (JFC) further lending support to the assertion that many cm- to m-sized objects on JFC orbits are asteroidal rather than cometary in origin.

An essential requirement for universal quantum control of the centre-of-mass motion of levitated objects is the development of a precise readout of all three translational degrees of freedom. Improving that precision presents one key challenge: collecting all the information on the object's position, encoded in the scattered light, equivalent to minimising the measurement imprecision. Here, we propose a new detection technique based on spatial mode decomposition, which addresses this problem using a simple integrated setup, where all of the light back-scattered from a levitated nanoparticle couples into a spatial mode sorter. We observe that each mode of the sorter pairs predominantly to the in-elastically scattered field generated by the object's motion along a particular spatial axis. This results in each translational degree of freedom being selectively encoded in the amplitude of orthogonal information channels. Using this approach, we report measurement efficiencies ($\eta_{\mathrm{tot}}^x$, $\eta_{\mathrm{tot}}^y$, $\eta_{\mathrm{tot}}^z$) = (0.17, 0.15, 0.30), implying that our technique is capable of reaching the 3D quantum ground state.

Recent successes in applying reinforcement learning (RL) for robotics has shown it is a viable approach for constructing robotic controllers. However, RL controllers can produce many collisions in environments where new obstacles appear during execution. This poses a problem in safety-critical settings. We present a hybrid approach, called iKinQP-RL, that uses an Inverse Kinematics Quadratic Programming (iKinQP) controller to correct actions proposed by an RL policy at runtime. This ensures safe execution in the presence of new obstacles not present during training. Preliminary experiments illustrate our iKinQP-RL framework completely eliminates collisions with new obstacles while maintaining a high task success rate.