Search

The JETSCAPE Collaboration reports a new determination of the jet transport parameter $\hat{q}$ in the Quark-Gluon Plasma (QGP) using Bayesian Inference, incorporating all available inclusive hadron and jet yield suppression data measured in heavy-ion collisions at RHIC and the LHC. This multi-observable analysis extends the previously published JETSCAPE Bayesian Inference determination of $\hat{q}$, which was based solely on a selection of inclusive hadron suppression data. JETSCAPE is a modular framework incorporating detailed dynamical models of QGP formation and evolution, and jet propagation and interaction in the QGP. Virtuality-dependent partonic energy loss in the QGP is modeled as a thermalized weakly-coupled plasma, with parameters determined from Bayesian calibration using soft-sector observables. This Bayesian calibration of $\hat{q}$ utilizes Active Learning, a machine--learning approach, for efficient exploitation of computing resources. The experimental data included in this analysis span a broad range in collision energy and centrality, and in transverse momentum. In order to explore the systematic dependence of the extracted parameter posterior distributions, several different calibrations are reported, based on combined jet and hadron data; on jet or hadron data separately; and on restricted kinematic or centrality ranges of the jet and hadron data. Tension is observed in comparison of these variations, providing new insights into the physics of jet transport in the QGP and its theoretical formulation., Comment: 20 pages, 10 figures, 2 tables, submitted to PRC; updated acknowledgements

We present dynamical modeling of the broad-line region (BLR) for the highly variable AGN SDSS J141041.25+531849.0 ($z = 0.359$) using photometric and spectroscopic monitoring data from the Sloan Digital Sky Survey Reverberation Mapping project and the SDSS-V Black Hole Mapper program, spanning from early 2013 to early 2023. We model the geometry and kinematics of the BLR in the H$\beta$, H$\alpha$, and MgII, emission lines for three different time periods to measure the potential change of structure within the BLR across time and line species. We consistently find a moderately edge-on $(i_{\rm full-state} = 53.29^{\circ} \,{}^{+7.29}_{-6.55})$ thick-disk $(\theta_{\rm opn, \; full-state} = 54.86^{\circ} \,{}^{+5.83}_{-4.74})$ geometry for all BLRs, with a joint estimate for the mass of the supermassive black hole (SMBH) for each of three time periods, yielding $\log_{10}(M_{\rm BH} / M_{\odot}) = 7.66^{+0.12}_{-0.13}$ when using the full dataset. The inferred individual virial factor $f$ $\sim 1$ is significantly smaller than the average factor for a local sample of dynamically modeled AGNs. There is strong evidence for non-virial motion, with over $80\%$ of clouds on inflowing/outflowing orbits. We analyze the change in model parameters across emission lines, finding the radii of BLRs for the emission lines are consistent with the following relative sizes $R_{\rm H\beta} \lesssim R_{\rm MgII } \lesssim R_{\rm H\alpha}$. Comparing results across time, we find $R_{\rm low-state} \lesssim R_{\rm high-state}$, with the change in BLR size for H$\beta$, being more significant than for the other two lines. The data also reveal complex, time-evolving, and potentially transient dynamics of the BLR gas over decade-long timescales, encouraging for future dynamical modeling of fine-scale BLR kinematics., Comment: 32 pages, 16 figures, 4 tables. Submitted to ApJ

A new framework, called x-scape, for the combined study of both hard and soft transverse momentum sectors in high energy proton-proton ($p$-$p$) and proton-nucleus ($p$-$A$) collisions is set up. A dynamical initial state is set up using the 3d-Glauber model with transverse locations of hotspots within each incoming nucleon. A hard scattering that emanates from two colliding hotspots is carried out using the Pythia generator. Initial state radiation from the incoming hard partons is carried out in a new module called I-matter, which includes the longitudinal location of initial splits. The energy-momentum of both the initial hard partons and their associated beam remnants is removed from the hot spots, depleting the energy-momentum available for the formation of the bulk medium. Outgoing showers are simulated using the matter generator, and results are presented for both cases, allowing for and not allowing for energy loss. First comparisons between this hard-soft model and single inclusive hadron and jet data from $p$-$p$ and minimum bias $p$-$Pb$ collisions are presented. Single hadron spectra in $p$-$p$ are used to carry out a limited (in number of parameters) Bayesian calibration of the model. Fair comparisons with data are indicative of the utility of this new framework. Theoretical studies of the correlation between jet $p_T$ and event activity at mid and forward rapidity are carried out., Comment: 18 pages, 15 figures

Biomolecular condensates play a central role in the spatial organization of living matter. Their formation is now well understood as a form of liquid-liquid phase separation that occurs very far from equilibrium. For instance, they can be modeled as active droplets, where the combination of molecular interactions and chemical reactions result in microphase separation. However, so far, models of chemically active droplets are spatially continuous and deterministic. Therefore, the relationship between the microscopic parameters of the models and some crucial properties of active droplets (such as their polydispersity, their shape anisotropy, or their typical lifetime) is yet to be established. In this work, we address this question computationally, using Brownian dynamics simulations of chemically active droplets: the building blocks are represented explicitly as particles that interact with attractive or repulsive interactions, depending on whether they are in a droplet-forming state or not. Thanks to this microscopic and stochastic view of the problem, we reveal how driving the system away from equilibrium in a controlled way determines the fluctuations and dynamics of active emulsions.

The past two decades have yielded several new measurements and reanalysis of older measurements of the neutron lifetime. These have led to a 4.4 standard deviation discrepancy between the most precise measurements of the neutron decay rate producing protons in cold neutron beams and the most precise lifetime measured in neutron storage experiments. Here we publish an analysis of the recently published UCN aimed a searching for an explanation of this difference using the model proposed by Koch and Hummel.

Chemical reaction network theory is a field of applied mathematics concerned with modeling chemical systems, and can be used in other contexts such as in systems biology to study cellular signaling pathways or epidemiology to study the effect of human interaction on the spread of disease. In this paper, we seek to understand a chemical reaction network's equilibrium points through the lens of algebraic geometry by computing the positive part of the steady-state variety defined by polynomial equations arising from the assumption of mass-action kinetics. We provide a systematic classification of positive steady-state varieties produced by 2-species, 2-reaction networks, grounded in combinatorial and algebraic properties. While some (restricted) techniques exist to fully understand the ideal defining the positive steady-state variety, this computation presents a significant challenge in general. Our classification theorems provide a simplification of previous criteria, and aim to provide a foundation for future analysis of larger networks.

We present the results of an investigation of a highly variable CIV broad absorption-line feature in the quasar SBS 1408+544 (z=2.337) that shows a significant shift in velocity over time. This source was observed as a part of the Sloan Digital Sky Survey Reverberation Mapping Project and the SDSS-V Black Hole Mapper Reverberation Mapping Project, and has been included in two previous studies, both of which identified significant variability in a high-velocity CIV broad absorption line (BAL) on timescales of just a few days in the quasar rest frame. Using ~130 spectra acquired over eight years of spectroscopic monitoring with SDSS, we have determined that this BAL is not only varying in strength, but is also systematically shifting to higher velocities. Using cross-correlation methods, we measure the velocity shifts (and corresponding acceleration) of the BAL on a wide range of timescales, measuring an overall velocity shift of delta v = -683 (+89, -84) km s-1 over the 8-year monitoring period. This corresponds to an average rest-frame acceleration of a=1.04 (+0.14, -0.13) cm s-2, though the magnitude of the acceleration on shorter timescales is not constant throughout. We place our measurements in the context of BAL-acceleration models and examine various possible causes of the observed velocity shift., Comment: Published in the Astrophysical Journal

Over 85 million computed tomography (CT) scans are performed annually in the US, of which approximately one quarter focus on the abdomen. Given the current radiologist shortage, there is a large impetus to use artificial intelligence to alleviate the burden of interpreting these complex imaging studies. Prior state-of-the-art approaches for automated medical image interpretation leverage vision language models (VLMs). However, current medical VLMs are generally limited to 2D images and short reports, and do not leverage electronic health record (EHR) data for supervision. We introduce Merlin - a 3D VLM that we train using paired CT scans (6+ million images from 15,331 CTs), EHR diagnosis codes (1.8+ million codes), and radiology reports (6+ million tokens). We evaluate Merlin on 6 task types and 752 individual tasks. The non-adapted (off-the-shelf) tasks include zero-shot findings classification (31 findings), phenotype classification (692 phenotypes), and zero-shot cross-modal retrieval (image to findings and image to impressions), while model adapted tasks include 5-year disease prediction (6 diseases), radiology report generation, and 3D semantic segmentation (20 organs). We perform internal validation on a test set of 5,137 CTs, and external validation on 7,000 clinical CTs and on two public CT datasets (VerSe, TotalSegmentator). Beyond these clinically-relevant evaluations, we assess the efficacy of various network architectures and training strategies to depict that Merlin has favorable performance to existing task-specific baselines. We derive data scaling laws to empirically assess training data needs for requisite downstream task performance. Furthermore, unlike conventional VLMs that require hundreds of GPUs for training, we perform all training on a single GPU., Comment: 18 pages, 7 figures

Numerical solvers of partial differential equations (PDEs) have been widely employed for simulating physical systems. However, the computational cost remains a major bottleneck in various scientific and engineering applications, which has motivated the development of reduced-order models (ROMs). Recently, machine-learning-based ROMs have gained significant popularity and are promising for addressing some limitations of traditional ROM methods, especially for advection dominated systems. In this chapter, we focus on a particular framework known as Latent Space Dynamics Identification (LaSDI), which transforms the high-fidelity data, governed by a PDE, to simpler and low-dimensional latent-space data, governed by ordinary differential equations (ODEs). These ODEs can be learned and subsequently interpolated to make ROM predictions. Each building block of LaSDI can be easily modulated depending on the application, which makes the LaSDI framework highly flexible. In particular, we present strategies to enforce the laws of thermodynamics into LaSDI models (tLaSDI), enhance robustness in the presence of noise through the weak form (WLaSDI), select high-fidelity training data efficiently through active learning (gLaSDI, GPLaSDI), and quantify the ROM prediction uncertainty through Gaussian processes (GPLaSDI). We demonstrate the performance of different LaSDI approaches on Burgers equation, a non-linear heat conduction problem, and a plasma physics problem, showing that LaSDI algorithms can achieve relative errors of less than a few percent and up to thousands of times speed-ups.

The third ML4H symposium was held in person on December 10, 2023, in New Orleans, Louisiana, USA. The symposium included research roundtable sessions to foster discussions between participants and senior researchers on timely and relevant topics for the \ac{ML4H} community. Encouraged by the successful virtual roundtables in the previous year, we organized eleven in-person roundtables and four virtual roundtables at ML4H 2022. The organization of the research roundtables at the conference involved 17 Senior Chairs and 19 Junior Chairs across 11 tables. Each roundtable session included invited senior chairs (with substantial experience in the field), junior chairs (responsible for facilitating the discussion), and attendees from diverse backgrounds with interest in the session's topic. Herein we detail the organization process and compile takeaways from these roundtable discussions, including recent advances, applications, and open challenges for each topic. We conclude with a summary and lessons learned across all roundtables. This document serves as a comprehensive review paper, summarizing the recent advancements in machine learning for healthcare as contributed by foremost researchers in the field., Comment: ML4H 2023, Research Roundtables

Prompt photons are created in the early stages of heavy ion collisions and traverse the QGP medium without any interaction. Therefore, photon-triggered jets can be used to study the jet quenching in the QGP medium. In this work, photon-triggered jets are studied through different jet and jet substructure observables for different collision systems and energies using the JETSCAPE framework. Since the multistage evolution used in the JETSCAPE framework is adequate to describe a wide range of experimental observables simultaneously using the same parameter tune, we use the same parameters tuned for jet and leading hadron studies. The same isolation criteria used in the experimental analysis are used to identify prompt photons for better comparison. For the first time, high-accuracy JETSCAPE results are compared with multi-energy LHC and RHIC measurements to better understand the deviations observed in prior studies. This study highlights the importance of multistage evolution for the simultaneous description of experimental observables through different collision systems and energies using a single parameter tune., Comment: 4 pages, 5 figures. Proceedings of Quark Matter 2023 - XXXth International Conference on Ultra-relativistic Nucleus-Nucleus Collisions, Houston, TX, 3-9 September 2023

The JETSCAPE Collaboration reports the first multi-messenger study of the QGP jet transport parameter $\hat{q}$ using Bayesian inference, incorporating all available hadron and jet inclusive yield and jet substructure data from RHIC and the LHC. The theoretical model utilizes virtuality-dependent in-medium partonic energy loss coupled to a detailed dynamical model of QGP evolution. Tension is observed when constraining $\hat{q}$ for different kinematic cuts of the inclusive hadron data. The addition of substructure data is shown to improve the constraint on $\hat{q}$, without inducing tension with the constraint due to inclusive observables. These studies provide new insight into the mechanisms of jet interactions in matter, and point to next steps in the field for comprehensive understanding of jet quenching as a probe of the QGP., Comment: 4 pages, 2 figures. Proceedings of Quark Matter 2023 - XXXth International Conference on Ultra-relativistic Nucleus-Nucleus Collisions, Houston, TX, 3-9 September 2023

"Changing-look" active galactic nuclei (CL-AGNs) challenge our basic ideas about the physics of accretion flows and circumnuclear gas around supermassive black holes. Using first-year Sloan Digital Sky Survey V (SDSS-V) repeated spectroscopy of nearly 29,000 previously known AGNs, combined with dedicated follow-up spectroscopy, and publicly available optical light curves, we have identified 116 CL-AGNs where (at least) one broad emission line has essentially (dis-)appeared, as well as 88 other extremely variable systems. Our CL-AGN sample, with 107 newly identified cases, is the largest reported to date, and includes $\sim0.4\%$ of the AGNs reobserved in first-year SDSS-V operations. Among our CL-AGNs, 67% exhibit dimming while 33% exhibit brightening. Our sample probes extreme AGN spectral variability on months to decades timescales, including some cases of recurring transitions on surprisingly short timescales ($\lesssim 2$ months in the rest frame). We find that CL events are preferentially found in lower-Eddington-ratio ($f_{Edd}$) systems: Our CL-AGNs have a $f_{Edd}$ distribution that significantly differs from that of a carefully constructed, redshift- and luminosity-matched control sample (Anderson-Darling test yielding $p_{\rm AD}\approx 6\times10^{-5}$; median $f_{Edd}\approx0.025$ vs. $0.043$). This preference for low $f_{Edd}$ strengthens previous findings of higher CL-AGN incidence at lower $f_{Edd}$, found in smaller samples. Finally, we show that the broad MgII emission line in our CL-AGN sample tends to vary significantly less than the broad H$\beta$ emission line. Our large CL-AGN sample demonstrates the advantages and challenges in using multi-epoch spectroscopy from large surveys to study extreme AGN variability and physics., Comment: Submitted to ApJ. Full tables and figure-sets will be published upon acceptance, and can be made available upon request$.$

Considerable information about the early-stage dynamics of heavy-ion collisions is encoded in the rapidity dependence of measurements. To leverage the large amount of experimental data, we perform a systematic analysis using three-dimensional hydrodynamic simulations of multiple collision systems -- large and small, symmetric and asymmetric. Specifically, we perform fully 3D multi-stage hydrodynamic simulations initialized by a parameterized model for rapidity-dependent energy deposition, which we calibrate on the hadron multiplicity and anisotropic flow coefficients. We utilize Bayesian inference to constrain properties of the early- and late- time dynamics of the system, and highlight the impact of enforcing global energy conservation in our 3D model.

By its nature, the so-called social cost of carbon (SCC(t)) will likely not cover the cost induced by climate change (damage cost and abatement cost) if it is used as a CO$_2$-price. It is a marginal price only. We define an implied CO$_2$-price that covers the climate change-induced costs. The price can be interpreted as a \textit{polluter pays principle}. A numerical analysis using a classical DICE model reveals that the cost-implied CO$_2$ price is around 500 \$/tCO$_2$, while the corresponding price associated with the SCC is about 50 \$/tCO$_2$. In addition, we define the internal rate of return of carbon abatement and calculate it for the classical DICE model. This rate is much higher than the model's discount rate, which may suggest the advantage of financing abatement by loans., Comment: 17 pages, 4 figures

Today's decisions on climate change mitigation affect the damage that future generations will bear. Discounting future benefits and costs of climate change mitigation is one of the most critical components of assessing efficient climate mitigation pathways. We extend the DICE model with stochastic discount rates to reflect the uncertain nature of discount rates. Stochastic rates give rise to a stochastic mitigation strategy, resulting in all model quantities becoming stochastic. We show that the classical calibration of the DICE model induces intergenerational inequality: future generations have to bear higher costs relative to GDP. Further, we show that considering stochastic discount rates and stochastic abatement policies, which can be interpreted as successive re-calibration, increases intergenerational inequality (and adds additional risks). Motivated by this, we consider additional financing risks by investigating two modifications of DICE. We find that allowing financing of abatement costs and considering non-linear financing effects for large damages improves intergenerational effort sharing. To conclude our discussion of options to improve intergenerational equity in an IAM, we propose a modified optimization to keep costs below 3% of GDP, resulting in more equal distribution of efforts between generations., Comment: 36 pages, 10 figures, arXiv admin note: text overlap with arXiv:2309.16186

Foundation models hold promise for transforming AI in healthcare by providing modular components that are easily adaptable to downstream healthcare tasks, making AI development more scalable and cost-effective. Structured EHR foundation models, trained on coded medical records from millions of patients, demonstrated benefits including increased performance with fewer training labels, and improved robustness to distribution shifts. However, questions remain on the feasibility of sharing these models across different hospitals and their performance for local task adaptation. This multi-center study examined the adaptability of a recently released structured EHR foundation model ($FM_{SM}$), trained on longitudinal medical record data from 2.57M Stanford Medicine patients. Experiments were conducted using EHR data at The Hospital for Sick Children and MIMIC-IV. We assessed both adaptability via continued pretraining on local data, and task adaptability compared to baselines of training models from scratch at each site, including a local foundation model. We evaluated the performance of these models on 8 clinical prediction tasks. In both datasets, adapting the off-the-shelf $FM_{SM}$ matched the performance of GBM models locally trained on all data while providing a 13% improvement in settings with few task-specific training labels. With continued pretraining on local data, label efficiency substantially improved, such that $FM_{SM}$ required fewer than 1% of training examples to match the fully trained GBM's performance. Continued pretraining was also 60 to 90% more sample-efficient than training local foundation models from scratch. Our findings show that adapting shared EHR foundation models across hospitals provides improved prediction performance at less cost, underscoring the utility of base foundation models as modular components to streamline the development of healthcare AI., Comment: 46 pages, 5 figures, 3 tables, 14 appendices

Synthesizing information from multiple data sources plays a crucial role in the practice of modern medicine. Current applications of artificial intelligence in medicine often focus on single-modality data due to a lack of publicly available, multimodal medical datasets. To address this limitation, we introduce INSPECT, which contains de-identified longitudinal records from a large cohort of patients at risk for pulmonary embolism (PE), along with ground truth labels for multiple outcomes. INSPECT contains data from 19,402 patients, including CT images, radiology report impression sections, and structured electronic health record (EHR) data (i.e. demographics, diagnoses, procedures, vitals, and medications). Using INSPECT, we develop and release a benchmark for evaluating several baseline modeling approaches on a variety of important PE related tasks. We evaluate image-only, EHR-only, and multimodal fusion models. Trained models and the de-identified dataset are made available for non-commercial use under a data use agreement. To the best of our knowledge, INSPECT is the largest multimodal dataset integrating 3D medical imaging and EHR for reproducible methods evaluation and research.

We propose a lean and functional transaction scheme to establish a secure delivery-versus-payment across two blockchains, where a) no intermediary is required and b) the operator of the payment chain/payment system has a small overhead and does not need to store state. The main idea comes with two requirements: First, the payment chain operator hosts a stateless decryption service that allows decrypting messages with his secret key. Second, a "Payment Contract" is deployed on the payment chain that implements a function transferAndDecrypt(uint id, address from, address to, string keyEncryptedSuccess, string keyEncryptedFail) that processes the (trigger-based) payment and emits the decrypted key depending on the success or failure of the transaction. The respective key can then trigger an associated transaction, e.g. claiming delivery by the buyer or re-claiming the locked asset by the seller.

In this talk we review jet production in a large variety of collision systems using the JETSCAPE event generator and Hybrid Hadronization. Hybrid Hadronization combines quark recombination, applicable when distances between partons in phase space are small, and string fragmentation appropriate for dilute parton systems. It can therefore smoothly describe the transition from very dilute parton systems like $e^++e^-$ to full $A+A$ collisions. We test this picture by using JETSCAPE to generate jets in various systems. Comparison to experimental data in $e^++e^-$ and $p+p$ collisions allows for a precise tuning of vacuum baseline parameters in JETSCAPE and Hybrid Hadronization. Proceeding to systems with jets embedded in a medium, we study in-medium hadronization for jet showers. We quantify the effects of an ambient medium, focusing in particular on the dependence on the collective flow and size of the medium. Our results clarify the effects we expect from in-medium hadronization of jets on observables like fragmentation functions, hadron chemistry and jet shape., Comment: 6 pages, 4 figures, Hard Probes 2023 conference, accepted for publication in Proceedings of Science: version 2, references added, typos fixed

Assessing the costs of climate change is essential to finding efficient pathways for the transition to a net-zero emissions economy, which is necessary to stabilise global temperatures at any level. In evaluating the benefits and costs of climate change mitigation, the discount rate converting future damages and costs into net-present values influences the timing of mitigation. Here, we amend the DICE model with a stochastic interest rate model to consider the uncertainty of discount rates in the future. Since abatement reduces future damages, changing interest rates renders abatement investments more or less beneficial. Stochastic interest rates will hence lead to a stochastic abatement strategy. We introduce a simple stochastic abatement model and show that this can increase intergenerational inequality concerning cost and risk. Analysing the sensitivities of the model calibration analytically and numerically exhibits that intergenerational inequality is a consequence of the DICE model calibration (and maybe that of IAMs in general). We then show that introducing funding of abatement costs reduces the variation of future cash-flows, which occur at different times but are off-setting in their net-present value. This effect can be interpreted as improving intergenerational effort sharing, which might be neglected in classical optimisation. This mechanism is amplified, including dependence of the interest rate risk on the amount of debt to be financed, i.e. considering the limited capacity of funding sources. As an alternative policy optimisation method, we propose limiting the total cost of damages and abatement below a fixed level relative to GDP - this modification induces equality between generations compared to their respective economic welfare, inducing early and fast mitigation of climate change to keep the total cost of climate change below 3% of global GDP.

Modeling the couplings between active particles often neglects the possible many-body effects that control the propulsion mechanism. Accounting for such effects requires the explicit modeling of the molecular details at the origin of activity. Here, we take advantage of a recent two-dimensional model of isotropic active particles whose propulsion originates from the interactions between solute particles in the bath. The colloid catalyzes a chemical reaction in its vicinity, which results in a local phase separation of solute particles, and the density fluctuations of solute particles cause the enhanced diffusion of the colloid. In this paper, we investigate an assembly of such active particles, using (i) an explicit model, where the microscopic dynamics of the solute particles is accounted for; and (ii) an implicit model, whose parameters are inferred from the explicit model at infinite dilution. In the explicit solute model, the effective diffusion coefficient of the active colloids strongly decreases with density, an effect which is not captured by the derived implicit model. This suggests that classical models, which usually decouple pair interactions from activity, fail to describe collective dynamics in active colloidal systems driven by solute-solute interactions.

Evaporation of cloud droplets accelerates when turbulence mixes dry air into the cloud, affecting droplet-size distributions in atmospheric clouds, combustion sprays, and jets of exhaled droplets. The challenge is to model local correlations between droplet numbers, sizes, and supersaturation, which determine supersaturation fluctuations along droplet paths (Lagrangian fluctuations). We derived a statistical model that accounts for these correlations. Its predictions are in quantitative agreement with results of direct numerical simulations, and it explains the key mechanisms at play., Comment: 6 pages, 3 figures, supplemental material

This work studies the relationship between accretion-disk size and quasar properties, using a sample of 95 quasars from the SDSS-RM project with measured lags between the $g$ and $i$ photometric bands. Our sample includes disk lags that are both longer and shorter than predicted by the \citet{SS73} model, requiring explanations which satisfy both cases. Although our quasars each have one lag measurement, we explore the wavelength-dependent effects of diffuse broad line region (BLR) contamination through our sample's broad redshift range, $0.1

The proper treatment of hadronic resonances plays an important role in many aspects of heavy ion collisions. This is expected to be the case also for hadronization, due to the large degeneracies of excited states, and the abundant production of hadrons from their decays. We first show how a comprehensive treatment of excited meson states can be incorporated into quark recombination, and in extension, into Hybrid Hadronization. We then discuss the quantum mechanics of forming excited states, utilizing the Wigner distribution functions of angular momentum eigenstates of isotropic 3-D harmonic oscillators. We further describe how resonance decays can be handled, based on a set of minimal assumptions, by creating an extension of hadron decays in PYTHIA 8. Finally, we present first results by simulating $e^+e^-$ collisions using PYTHIA and Hybrid Hadronization with excited mesons up to orbital angular momentum $L=4$ and radial quantum number 2. We find that states up to $L=2$ are produced profusely by quark recombination., Comment: 6 pages, 3 figures; submitted to Proceedings of Science

The ability of large language models (LLMs) to follow natural language instructions with human-level fluency suggests many opportunities in healthcare to reduce administrative burden and improve quality of care. However, evaluating LLMs on realistic text generation tasks for healthcare remains challenging. Existing question answering datasets for electronic health record (EHR) data fail to capture the complexity of information needs and documentation burdens experienced by clinicians. To address these challenges, we introduce MedAlign, a benchmark dataset of 983 natural language instructions for EHR data. MedAlign is curated by 15 clinicians (7 specialities), includes clinician-written reference responses for 303 instructions, and provides 276 longitudinal EHRs for grounding instruction-response pairs. We used MedAlign to evaluate 6 general domain LLMs, having clinicians rank the accuracy and quality of each LLM response. We found high error rates, ranging from 35% (GPT-4) to 68% (MPT-7B-Instruct), and an 8.3% drop in accuracy moving from 32k to 2k context lengths for GPT-4. Finally, we report correlations between clinician rankings and automated natural language generation metrics as a way to rank LLMs without human review. We make MedAlign available under a research data use agreement to enable LLM evaluations on tasks aligned with clinician needs and preferences.

Understanding the modification of jets and high-$p_T$ probes in small systems requires the integration of soft and hard physics. We present recent developments in extending the JETSCAPE framework to build an event generator, which includes correlations between soft and hard partons, to study jet observables in small systems. The multi-scale physics of the collision is separated into different stages. Hard scatterings are first sampled at binary collision positions provided by the Glauber geometry. They are then propagated backward in space-time following an initial-state shower to obtain the initiating partons' energies and momenta before the collision. These energies and momenta are then subtracted from the incoming colliding nucleons for soft-particle production, modeled by the 3D-Glauber + hydrodynamics + hadronic transport framework. This new hybrid approach (X-SCAPE) includes non-trivial correlations between jet and soft particle productions in small systems. We calibrate this framework with the final state hadrons' $p_T$-spectra from low to high $p_T$ in $p$-$p$, and and then compare with the spectra in $p$-$Pb$ collisions from the LHC. We also present results for additional observables such as the distributions of event activity as a function of the hardest jet $p_T$ in forward and mid-rapidity for both $p$-$p$ and $p$-$Pb$ collisions., Comment: 6 pages, 6 figures, Proceedings of Hard Probes 2023, 26-31 March 2023 Aschaffenburg, Germany

We extend the notion of a spectral triple to that of a higher-order relative spectral triple, which accommodates several types of hypoelliptic differential operators on manifolds with boundary. The bounded transform of a higher-order relative spectral triple gives rise to a relative K-homology cycle. In the case of an elliptic differential operator on a compact smooth manifold with boundary, we calculate the K-homology boundary map of the constructed relative K-homology cycle to obtain a generalization of the Baum-Douglas-Taylor index theorem., Comment: 34 pages

We study parton energy-momentum exchange with the quark gluon plasma (QGP) within a multistage approach composed of in-medium DGLAP evolution at high virtuality, and (linearized) Boltzmann Transport formalism at lower virtuality. This multistage simulation is then calibrated in comparison with high $p_T$ charged hadrons, D-mesons, and the inclusive jet nuclear modification factors, using Bayesian model-to-data comparison, to extract the virtuality-dependent transverse momentum broadening transport coefficient $\hat{q}$. To facilitate this undertaking, we develop a quantitative metric for validating the Bayesian workflow, which is used to analyze the sensitivity of various model parameters to individual observables. The usefulness of this new metric in improving Bayesian model emulation is shown to be highly beneficial for future such analyses., Comment: 24 pages, 16 figures

Shower development dynamics for a jet traveling through the quark-gluon plasma (QGP) is a multiscale process, where the heavy flavor mass is an important scale. During the high virtuality portion of the jet evolution in the QGP, emission of gluons from a heavy flavor is modified owing to heavy quark mass. Medium-induced radiation of heavy flavor is sensitive to microscopic processes (e.g. diffusion), whose virtuality dependence is phenomenologically explored in this study. In the lower virtuality part of shower evolution, i.e. when the mass is comparable to the virtuality of the parton, scattering and radiation processes of heavy quarks differ from light quarks. The effects of these mechanisms on shower development in heavy flavor tagged showers in the QGP is explored here. Furthermore, this multiscale study examines dynamical pair production of heavy flavor (via virtual gluon splittings) and their subsequent evolution in the QGP, which is not possible otherwise. A realistic event-by-event simulation is performed using the JETSCAPE framework. Energy-momentum exchange with the medium proceeds using a weak coupling recoil approach. Using leading hadron and open heavy flavor observables, differences in heavy versus light quark energy-loss mechanisms are explored, while the importance of heavy flavor pair production is highlighted along with future directions to study., Comment: 6 pages, 3 figures, contribution to the Hard Probes 2023 proceedings

We utilize event-by-event Monte Carlo simulations within the JETSCAPE framework to examine scale-dependent jet-medium interactions in heavy-ion collisions. The reduction in jet-medium interaction during the early high-virtuality stage, where the medium is resolved at a short distance scale, is emphasized as a key element in explaining multiple jet observables, particularly substructures, simultaneously. By employing the MATTER+LBT setup, which incorporates this explicit reduction of medium effects at high virtuality, we investigate jet substructure observables, such as Soft Drop groomed observables. When contrasted with existing data, our findings spotlight the significant influence of the reduction at the early high-virtuality stages. Furthermore, we study the substructure of gamma-tagged jets, providing predictive insights for future experimental analyses. This broadens our understanding of the various contributing factors involved in modifying jet substructures., Comment: 6 pages, 4 figures, contribution to the Hard Probes 2023 proceedings

While the general machine learning (ML) community has benefited from public datasets, tasks, and models, the progress of ML in healthcare has been hampered by a lack of such shared assets. The success of foundation models creates new challenges for healthcare ML by requiring access to shared pretrained models to validate performance benefits. We help address these challenges through three contributions. First, we publish a new dataset, EHRSHOT, which contains deidentified structured data from the electronic health records (EHRs) of 6,739 patients from Stanford Medicine. Unlike MIMIC-III/IV and other popular EHR datasets, EHRSHOT is longitudinal and not restricted to ICU/ED patients. Second, we publish the weights of CLMBR-T-base, a 141M parameter clinical foundation model pretrained on the structured EHR data of 2.57M patients. We are one of the first to fully release such a model for coded EHR data; in contrast, most prior models released for clinical data (e.g. GatorTron, ClinicalBERT) only work with unstructured text and cannot process the rich, structured data within an EHR. We provide an end-to-end pipeline for the community to validate and build upon its performance. Third, we define 15 few-shot clinical prediction tasks, enabling evaluation of foundation models on benefits such as sample efficiency and task adaptation. Our model and dataset are available via a research data use agreement from our website: https://ehrshot.stanford.edu. Code to reproduce our results are available at our Github repo: https://github.com/som-shahlab/ehrshot-benchmark