Search

We introduce Best-of-N (BoN) Jailbreaking, a simple black-box algorithm that jailbreaks frontier AI systems across modalities. BoN Jailbreaking works by repeatedly sampling variations of a prompt with a combination of augmentations - such as random shuffling or capitalization for textual prompts - until a harmful response is elicited. We find that BoN Jailbreaking achieves high attack success rates (ASRs) on closed-source language models, such as 89% on GPT-4o and 78% on Claude 3.5 Sonnet when sampling 10,000 augmented prompts. Further, it is similarly effective at circumventing state-of-the-art open-source defenses like circuit breakers. BoN also seamlessly extends to other modalities: it jailbreaks vision language models (VLMs) such as GPT-4o and audio language models (ALMs) like Gemini 1.5 Pro, using modality-specific augmentations. BoN reliably improves when we sample more augmented prompts. Across all modalities, ASR, as a function of the number of samples (N), empirically follows power-law-like behavior for many orders of magnitude. BoN Jailbreaking can also be composed with other black-box algorithms for even more effective attacks - combining BoN with an optimized prefix attack achieves up to a 35% increase in ASR. Overall, our work indicates that, despite their capability, language models are sensitive to seemingly innocuous changes to inputs, which attackers can exploit across modalities.

Defending large language models against jailbreaks so that they never engage in a broadly-defined set of forbidden behaviors is an open problem. In this paper, we investigate the difficulty of jailbreak-defense when we only want to forbid a narrowly-defined set of behaviors. As a case study, we focus on preventing an LLM from helping a user make a bomb. We find that popular defenses such as safety training, adversarial training, and input/output classifiers are unable to fully solve this problem. In pursuit of a better solution, we develop a transcript-classifier defense which outperforms the baseline defenses we test. However, our classifier defense still fails in some circumstances, which highlights the difficulty of jailbreak-defense even in a narrow domain., Comment: Accepted to the AdvML-Frontiers and SoLaR workshops at NeurIPS 2024

Humans acquire knowledge by observing the external world, but also by introspection. Introspection gives a person privileged access to their current state of mind (e.g., thoughts and feelings) that is not accessible to external observers. Can LLMs introspect? We define introspection as acquiring knowledge that is not contained in or derived from training data but instead originates from internal states. Such a capability could enhance model interpretability. Instead of painstakingly analyzing a model's internal workings, we could simply ask the model about its beliefs, world models, and goals. More speculatively, an introspective model might self-report on whether it possesses certain internal states such as subjective feelings or desires and this could inform us about the moral status of these states. Such self-reports would not be entirely dictated by the model's training data. We study introspection by finetuning LLMs to predict properties of their own behavior in hypothetical scenarios. For example, "Given the input P, would your output favor the short- or long-term option?" If a model M1 can introspect, it should outperform a different model M2 in predicting M1's behavior even if M2 is trained on M1's ground-truth behavior. The idea is that M1 has privileged access to its own behavioral tendencies, and this enables it to predict itself better than M2 (even if M2 is generally stronger). In experiments with GPT-4, GPT-4o, and Llama-3 models (each finetuned to predict itself), we find that the model M1 outperforms M2 in predicting itself, providing evidence for introspection. Notably, M1 continues to predict its behavior accurately even after we intentionally modify its ground-truth behavior. However, while we successfully elicit introspection on simple tasks, we are unsuccessful on more complex tasks or those requiring out-of-distribution generalization., Comment: 15 pages, 9 figures

In this paper we prove that there exists a compact perturbation of the Ricci flat Taub-Bolt metric that evolves under the Ricci flow into a finite time singularity modelled on the shrinking solition FIK [5]. Moreover, this perturbation can be made arbitrarily L^2-small with respect to the Taub-Bolt metric. The method of proof closely follows the strategy adopted in [14], where the author constructs, via a Wa\.zewski box argument, Ricci flows on compact manifolds which encounter finite singularities modelled on a given asymptotically conical shrinking soliton., Comment: arXiv admin note: text overlap with arXiv:2202.03386 by other authors

The integration of new modalities into frontier AI systems offers exciting capabilities, but also increases the possibility such systems can be adversarially manipulated in undesirable ways. In this work, we focus on a popular class of vision-language models (VLMs) that generate text outputs conditioned on visual and textual inputs. We conducted a large-scale empirical study to assess the transferability of gradient-based universal image ``jailbreaks" using a diverse set of over 40 open-parameter VLMs, including 18 new VLMs that we publicly release. Overall, we find that transferable gradient-based image jailbreaks are extremely difficult to obtain. When an image jailbreak is optimized against a single VLM or against an ensemble of VLMs, the jailbreak successfully jailbreaks the attacked VLM(s), but exhibits little-to-no transfer to any other VLMs; transfer is not affected by whether the attacked and target VLMs possess matching vision backbones or language models, whether the language model underwent instruction-following and/or safety-alignment training, or many other factors. Only two settings display partially successful transfer: between identically-pretrained and identically-initialized VLMs with slightly different VLM training data, and between different training checkpoints of a single VLM. Leveraging these results, we then demonstrate that transfer can be significantly improved against a specific target VLM by attacking larger ensembles of ``highly-similar" VLMs. These results stand in stark contrast to existing evidence of universal and transferable text jailbreaks against language models and transferable adversarial attacks against image classifiers, suggesting that VLMs may be more robust to gradient-based transfer attacks., Comment: NeurIPS 2024 Workshops: RBFM (Best Paper), Frontiers in AdvML (Oral), Red Teaming GenAI (Oral), SoLaR (Spotlight), SATA

The proliferation of generative models, combined with pretraining on web-scale data, raises a timely question: what happens when these models are trained on their own generated outputs? Recent investigations into model-data feedback loops proposed that such loops would lead to a phenomenon termed model collapse, under which performance progressively degrades with each model-data feedback iteration until fitted models become useless. However, those studies largely assumed that new data replace old data over time, where an arguably more realistic assumption is that data accumulate over time. In this paper, we ask: what effect does accumulating data have on model collapse? We empirically study this question by pretraining sequences of language models on text corpora. We confirm that replacing the original real data by each generation's synthetic data does indeed tend towards model collapse, then demonstrate that accumulating the successive generations of synthetic data alongside the original real data avoids model collapse; these results hold across a range of model sizes, architectures, and hyperparameters. We obtain similar results for deep generative models on other types of real data: diffusion models for molecule conformation generation and variational autoencoders for image generation. To understand why accumulating data can avoid model collapse, we use an analytically tractable framework introduced by prior work in which a sequence of linear models are fit to the previous models' outputs. Previous work used this framework to show that if data are replaced, the test error increases with the number of model-fitting iterations; we extend this argument to prove that if data instead accumulate, the test error has a finite upper bound independent of the number of iterations, meaning model collapse no longer occurs.

Simulations of pre-supernova evolution suggest that some intense shell burning can be so active that, in extreme cases, it can merge with the outer shell, changing the initial conditions for the supernova explosion. However, such violent activity in the interior of stars has been difficult to confirm from observations of stars. Here we propose that the elemental composition of O-rich ejecta in supernova remnants can be a tool to test for this kind of intense shell burning activity in the final stages of progenitor evolution. As an example, we discuss the origin of "Mg-rich" ejecta in the supernova remnant N49B. A high Mg/Ne mass ratio $\gtrsim 1$ suggests that the Ne- or O-burning shell has broken into or merged with the outer shell before the collapse. Such Mg-rich (or Ne-poor) ejecta has been identified in some other supernova remnants, supporting the idea that some destratification process, such as a shell merger, does indeed occur in the interiors of some massive stars, although they may not be the majority. Our results suggest that X-ray observations of O-rich ejecta in core-collapse supernova remnants will be a unique tool to probe the shell burning activity during the final stage of a massive star's interior., Comment: Submitted to ApJ and received positive peer review. This is the version in the first submission. 11 pages, 7 figures, 2 tables

The progenitor of the W49B supernova remnant is still under debate. One of the candidates is a jet-driven core-collapse supernova. In such a highly asymmetric explosion, a strong $\alpha$-rich freezeout is expected in local high entropy regions, which should enrich elements synthesized by the capture of $\alpha$-particles such as $^{44}$Ti and $^{48}$Cr (decaying to $^{44}$Ca and $^{48}$Ti, respectively). In the present work, in order to infer the progenitor of the W49B remnant, we constrain the amount of stable Ti ($^{48}$Ti) synthesized, using the {\it Suzaku} observation. We found no firm evidence for the Ti line and set the upper limit of $M_{\rm Ti}/M_{\rm Fe} < 8.2 \times$ 10$^{-4}$ (99\% limit using Xspec) and $M_{\rm Ti}/M_{\rm Fe} < 1.9 \times$ 10$^{-3}$ (99\% limit using SPEX), and thus excluded almost all hypernova/jet-driven supernova models. Our results, as complemented by some previous studies, suggest that a Type Ia supernova from a near-$M_{\rm Ch}$ (Chandrasekhar mass) white dwarf is the most favorable candidate for the origin of W49B. Future observations with X-ray calorimeter missions, such as XRISM, will give us a stronger constraint on the progenitor., Comment: Submitted to ApJ, 7 pages, 3 figures, 1 table

Common methods for aligning large language models (LLMs) with desired behaviour heavily rely on human-labelled data. However, as models grow increasingly sophisticated, they will surpass human expertise, and the role of human evaluation will evolve into non-experts overseeing experts. In anticipation of this, we ask: can weaker models assess the correctness of stronger models? We investigate this question in an analogous setting, where stronger models (experts) possess the necessary information to answer questions and weaker models (non-experts) lack this information. The method we evaluate is debate, where two LLM experts each argue for a different answer, and a non-expert selects the answer. We find that debate consistently helps both non-expert models and humans answer questions, achieving 76% and 88% accuracy respectively (naive baselines obtain 48% and 60%). Furthermore, optimising expert debaters for persuasiveness in an unsupervised manner improves non-expert ability to identify the truth in debates. Our results provide encouraging empirical evidence for the viability of aligning models with debate in the absence of ground truth., Comment: For code please check: https://github.com/ucl-dark/llm_debate

Understanding the behaviour of a system's API can be hard. Giving users access to relevant examples of how an API behaves has been shown to make this easier for them. In addition, such examples can be used to verify expected behaviour or identify unwanted behaviours. Methods for automatically generating examples have existed for a long time. However, state-of-the-art methods rely on either white-box information, such as source code, or on formal specifications of the system behaviour. But what if you do not have access to either? e.g., when interacting with a third-party API. In this paper, we present an approach to automatically generate relevant examples of behaviours of an API, without requiring either source code or a formal specification of behaviour. Evaluation on an industry-grade REST API shows that our method can produce small and relevant examples that can help engineers to understand the system under exploration.

Nebular-phase observations of peculiar Type Ia supernovae (SNe Ia) provide important constraints on progenitor scenarios and explosion dynamics for both these rare SNe and the more common, cosmologically useful SNe Ia. We present observations from an extensive ground-based and space-based follow-up campaign to characterize SN 2022pul, a "super-Chandrasekhar" mass SN Ia (alternatively "03fg-like" SN), from before peak brightness to well into the nebular phase across optical to mid-infrared (MIR) wavelengths. The early rise of the light curve is atypical, exhibiting two distinct components, consistent with SN Ia ejecta interacting with dense carbon-oxygen rich circumstellar material (CSM). In the optical, SN 2022pul is most similar to SN 2012dn, having a low estimated peak luminosity ($M_{B}=-18.9$ mag) and high photospheric velocity relative to other 03fg-like SNe. In the nebular phase, SN 2022pul adds to the increasing diversity of the 03fg-like subclass. From 168 to 336 days after peak $B$-band brightness, SN 2022pul exhibits asymmetric and narrow emission from [O I] $\lambda\lambda 6300,\ 6364$ (${\rm FWHM} \approx 2{,}000$ km s$^{-1}$), strong, broad emission from [Ca II] $\lambda\lambda 7291,\ 7323$ (${\rm FWHM} \approx 7{,}300$ km s$^{-1}$), and a rapid Fe III to Fe II ionization change. Finally, we present the first-ever optical-to-mid-infrared (MIR) nebular spectrum of an 03fg-like SN Ia using data from JWST. In the MIR, strong lines of neon and argon, weak emission from stable nickel, and strong thermal dust emission (with $T \approx 500$ K), combined with prominent [O I] in the optical, suggest that SN 2022pul was produced by a white dwarf merger within carbon/oxygen-rich CSM., Comment: 23 pages, 11 figures, submitted to ApJ

We present an analysis of ground-based and JWST observations of SN~2022pul, a peculiar "03fg-like" (or "super-Chandrasekhar") Type Ia supernova (SN Ia), in the nebular phase at 338d post explosion. Our combined spectrum continuously covers 0.4--14 $\mu$m and includes the first mid-infrared spectrum of an 03fg-like SN Ia. Compared to normal SN Ia 2021aefx, SN 2022pul exhibits a lower mean ionization state, asymmetric emission-line profiles, stronger emission from the intermediate-mass elements (IMEs) argon and calcium, weaker emission from iron-group elements (IGEs), and the first unambiguous detection of neon in a SN Ia. Strong, broad, centrally peaked [Ne II] line at 12.81 $\mu$m was previously predicted as a hallmark of "violent merger'' SN Ia models, where dynamical interaction between two sub-$M_{ch}$ white dwarfs (WDs) causes disruption of the lower mass WD and detonation of the other. The violent merger scenario was already a leading hypothesis for 03fg-like SNe Ia; in SN 2022pul it can explain the large-scale ejecta asymmetries seen between the IMEs and IGEs and the central location of narrow oxygen and broad neon. We modify extant models to add clumping of the ejecta to better reproduce the optical iron emission, and add mass in the innermost region ($< 2000$ km s$^{-1}$) to account for the observed narrow [O I]~$\lambda\lambda6300$, 6364 emission. A violent WD-WD merger explains many of the observations of SN 2022pul, and our results favor this model interpretation for the subclass of 03fg-like SN Ia., Comment: 20 pages, 10 figures, published in ApJ

Neutrino interactions with stellar material are widely believed to be fundamental to the explosion of massive stars. However, this important process has remained difficult to confirm observationally. We propose a new method to verify it using X-ray observations of the supernova remnant Cassiopeia A. The elemental composition in its Fe-rich ejecta that could have been produced at the innermost region of the supernova, where neutrinos are expected to interact, allows us to examine the presence of neutrino interactions. Here we demonstrate that the amount of Mn produced without neutrino nucleosynthesis processes (i.e., the $\nu$- and $\nu$p-process) is too small to explain the Mn/Fe mass ratio we measure (0.14--0.67\%). This result supports the operation of significant neutrino interactions in the Cassiopeia A supernova. If the observed Mn/Fe mass ratio purely reflects the production at the innermost region of the supernova, this would be the first robust confirmation of neutrino-matter interactions in an individual supernova. We further show that the Mn/Fe mass ratio has the potential to constrain supernova neutrino parameters (i.e., total neutrino luminosity, neutrino temperature). Future spatially-resolved, high-resolution X-ray spectroscopy will allow us to investigate the details of neutrino-supernova astrophysics through its signatures in elemental composition not only in Cassiopeia A but also in other remnants., Comment: 17 pages, 10 figures. Accepted for publication in ApJ

Observations of the millimeter sky contain valuable information on a number of signals, including the blackbody cosmic microwave background (CMB), Galactic emissions, and the Compton-$y$ distortion due to the thermal Sunyaev-Zel'dovich (tSZ) effect. Extracting new insight into cosmological and astrophysical questions often requires combining multi-wavelength observations to spectrally isolate one component. In this work, we present a new arcminute-resolution Compton-$y$ map, which traces out the line-of-sight-integrated electron pressure, as well as maps of the CMB in intensity and E-mode polarization, across a third of the sky (around 13,000 sq.~deg.). We produce these through a joint analysis of data from the Atacama Cosmology Telescope (ACT) Data Release 4 and 6 at frequencies of roughly 93, 148, and 225 GHz, together with data from the \textit{Planck} satellite at frequencies between 30 GHz and 545 GHz. We present detailed verification of an internal linear combination pipeline implemented in a needlet frame that allows us to efficiently suppress Galactic contamination and account for spatial variations in the ACT instrument noise. These maps provide a significant advance, in noise levels and resolution, over the existing \textit{Planck} component-separated maps and will enable a host of science goals including studies of cluster and galaxy astrophysics, inferences of the cosmic velocity field, primordial non-Gaussianity searches, and gravitational lensing reconstruction of the CMB., Comment: The Compton-y map and associated products will be made publicly available upon publication of the paper. The CMB T and E mode maps will be made available when the DR6 maps are made public

The community agrees that Type Ia supernovae arise from Carbon/Oxygen white dwarfs undergoing thermonuclear runaway. However, the full progenitor system and the process that prompts the white dwarf to explode remain unknown. Most current models suggest that the white dwarf explodes because of interaction with a binary companion which may survive the process and remain within the resulting remnant of the exploded star. Furthermore, both the pre-supernova interaction process and the explosion of the primary are expected to imprint a significant departure from ordinary stellar radii and temperatures onto the secondary, making the star identifiable against the unrelated stellar population. Identification of a surviving companion inside an SN Ia remnant might confirm a specific corresponding SN Ia progenitor channel based on the identity of the companion. We conducted a surviving companion search of the Type Ia remnant SNR 0509-67.5 based in the Large Magellanic Cloud. The well-constrained distance to and foreground extinction of the Large Magellanic Cloud allow for Bayesian inference of stellar parameters with low correlation and uncertainties. We present a deep catalog of fully characterized stars interior to SNR 0509-67.5 with radii, effective temperatures, and metallicities inferred using combined Hubble Space Telescope photometric observations across multiple visits. We then compile a list of surviving companion models appropriate for the age of the remnant (roughly 400 years after the explosion). We compare these predictions with the inferred stellar parameters and conclude that none of the stars are consistent with the predicted signatures of a surviving companion., Comment: 16 pages, 8 figures, 4 tables

We present weak gravitational lensing measurements of a sample of 157 clusters within the Kilo Degree Survey (KiDS), detected with a $>5\sigma$ thermal Sunyaev-Zel'dovich (SZ) signal by the Atacama Cosmology Telescope (ACT). Using a halo-model approach we constrain the average total cluster mass, $M_{\rm WL}$, accounting for the ACT cluster selection function of the full sample. We find that the SZ cluster mass estimate $M_{\rm SZ}$, which was calibrated using X-ray observations, is biased with $M_{\rm SZ}/M_{\rm WL} = (1-b_{\rm SZ}) = 0.65\pm 0.05$. Separating the sample into six mass bins, we find no evidence of a strong mass-dependency for the mass bias, $(1-b_{\rm SZ})$. Adopting this ACT-KiDS SZ mass-calibration would bring the Planck SZ cluster count into agreement with the counts expected from the {\it Planck} cosmic microwave background $\Lambda$CDM cosmological model, although it should be noted that the cluster sample considered in this work has a lower average mass $M_{\rm SZ, uncor} = 3.64 \times 10^{14} M_{\odot}$ compared to the Planck cluster sample which has an average mass in the range $M_{\rm SZ, uncor} = (5.5-8.5) \times 10^{14} M_{\odot}$, depending on the sub-sample used., Comment: 12 pages, 7 figures

We present cosmological constraints from a gravitational lensing mass map covering 9400 sq. deg. reconstructed from CMB measurements made by the Atacama Cosmology Telescope (ACT) from 2017 to 2021. In combination with BAO measurements (from SDSS and 6dF), we obtain the amplitude of matter fluctuations $\sigma_8 = 0.819 \pm 0.015$ at 1.8% precision, $S_8\equiv\sigma_8({\Omega_{\rm m}}/0.3)^{0.5}=0.840\pm0.028$ and the Hubble constant $H_0= (68.3 \pm 1.1)\, \text{km}\,\text{s}^{-1}\,\text{Mpc}^{-1}$ at 1.6% precision. A joint constraint with CMB lensing measured by the Planck satellite yields even more precise values: $\sigma_8 = 0.812 \pm 0.013$, $S_8\equiv\sigma_8({\Omega_{\rm m}}/0.3)^{0.5}=0.831\pm0.023$ and $H_0= (68.1 \pm 1.0)\, \text{km}\,\text{s}^{-1}\,\text{Mpc}^{-1}$. These measurements agree well with $\Lambda$CDM-model extrapolations from the CMB anisotropies measured by Planck. To compare these constraints to those from the KiDS, DES, and HSC galaxy surveys, we revisit those data sets with a uniform set of assumptions, and find $S_8$ from all three surveys are lower than that from ACT+Planck lensing by varying levels ranging from 1.7-2.1$\sigma$. These results motivate further measurements and comparison, not just between the CMB anisotropies and galaxy lensing, but also between CMB lensing probing $z\sim 0.5-5$ on mostly-linear scales and galaxy lensing at $z\sim 0.5$ on smaller scales. We combine our CMB lensing measurements with CMB anisotropies to constrain extensions of $\Lambda$CDM, limiting the sum of the neutrino masses to $\sum m_{\nu} < 0.13$ eV (95% c.l.), for example. Our results provide independent confirmation that the universe is spatially flat, conforms with general relativity, and is described remarkably well by the $\Lambda$CDM model, while paving a promising path for neutrino physics with gravitational lensing from upcoming ground-based CMB surveys., Comment: 32 pages, 17 figures, replaced with version accepted in ApJ (Feb 2024). Cosmological likelihood data and mass maps are public here: https://lambda.gsfc.nasa.gov/product/act/actadv_prod_table.html ; likelihood software is here: https://github.com/ACTCollaboration/act_dr6_lenslike . Also see companion papers Qu et al and MacCrann et al

We present new measurements of cosmic microwave background (CMB) lensing over $9400$ sq. deg. of the sky. These lensing measurements are derived from the Atacama Cosmology Telescope (ACT) Data Release 6 (DR6) CMB dataset, which consists of five seasons of ACT CMB temperature and polarization observations. We determine the amplitude of the CMB lensing power spectrum at $2.3\%$ precision ($43\sigma$ significance) using a novel pipeline that minimizes sensitivity to foregrounds and to noise properties. To ensure our results are robust, we analyze an extensive set of null tests, consistency tests, and systematic error estimates and employ a blinded analysis framework. The baseline spectrum is well fit by a lensing amplitude of $A_{\mathrm{lens}}=1.013\pm0.023$ relative to the Planck 2018 CMB power spectra best-fit $\Lambda$CDM model and $A_{\mathrm{lens}}=1.005\pm0.023$ relative to the $\text{ACT DR4} + \text{WMAP}$ best-fit model. From our lensing power spectrum measurement, we derive constraints on the parameter combination $S^{\mathrm{CMBL}}_8 \equiv \sigma_8 \left({\Omega_m}/{0.3}\right)^{0.25}$ of $S^{\mathrm{CMBL}}_8= 0.818\pm0.022$ from ACT DR6 CMB lensing alone and $S^{\mathrm{CMBL}}_8= 0.813\pm0.018$ when combining ACT DR6 and Planck NPIPE CMB lensing power spectra. These results are in excellent agreement with $\Lambda$CDM model constraints from Planck or $\text{ACT DR4} + \text{WMAP}$ CMB power spectrum measurements. Our lensing measurements from redshifts $z\sim0.5$--$5$ are thus fully consistent with $\Lambda$CDM structure growth predictions based on CMB anisotropies probing primarily $z\sim1100$. We find no evidence for a suppression of the amplitude of cosmic structure at low redshifts, Comment: 45+22 pages, 50 figures. v2 matches with published version in ApJ. Cosmological likelihood data and lensing maps are here: https://lambda.gsfc.nasa.gov/product/act/actadv_prod_table.html ; likelihood software is here: https://github.com/ACTCollaboration/act_dr6_lenslike . Also see companion papers Madhavacheril et al and MacCrann et al

We conduct a systematic search for transients in three years of data (2017-2019) from the Atacama Cosmology Telescope (ACT). ACT covers 40 percent of the sky at three bands spanning from 77 GHz to 277 GHz. Analysis of 3-day mean-subtracted sky maps, which were match-filtered for point sources, yielded 29 transients detections. Eight of these transients are due to known asteroids, and three others were previously published. Four of these events occur in areas of with poor noise models and thus we cannot be confident they are real transients. We are left with 14 new transient events occurring at 11 unique locations. All of these events are associated with either rotationally variable stars or cool stars. Ten events have flat or falling spectra indicating radiation from synchrotron emission. One event has a rising spectrum indicating a different engine for the flare.

Count data with complex features arise in many disciplines, including ecology, agriculture, criminology, medicine, and public health. Zero inflation, spatial dependence, and non-equidispersion are common features in count data. There are two classes of models that allow for these features -- he mode-parameterized Conway--Maxwell--Poisson (COMP) distribution and the generalized Poisson model. However both require the use of either constraints on the parameter space or a parameterization that leads to challenges in interpretability. We propose a spatial mean-parameterized COMP model that retains the flexibility of these models while resolving the above issues. We use a Bayesian spatial filtering approach in order to efficiently handle high-dimensional spatial data and we use reversible-jump MCMC to automatically choose the basis vectors for spatial filtering. The COMP distribution poses two additional computational challenges -- an intractable normalizing function in the likelihood and no closed-form expression for the mean. We propose a fast computational approach that addresses these challenges by, respectively, introducing an efficient auxiliary variable algorithm and pre-computing key approximations for fast likelihood evaluation. We illustrate the application of our methodology to simulated and real datasets, including Texas HPV-cancer data and US vaccine refusal data.

Parental refusal and delay of childhood vaccination has increased in recent years in the United States. This phenomenon challenges maintenance of herd immunity and increases the risk of outbreaks of vaccine-preventable diseases. We examine US county-level vaccine refusal for patients under five years of age collected during the period 2012--2015 from an administrative healthcare dataset. We model these data with a Bayesian zero-inflated negative binomial regression model to capture social and political processes that are associated with vaccine refusal, as well as factors that affect our measurement of vaccine refusal.Our work highlights fine-scale socio-demographic characteristics associated with vaccine refusal nationally, finds that spatial clustering in refusal can be explained by such factors, and has the potential to aid in the development of targeted public health strategies for optimizing vaccine uptake.

We present JWST near- and mid-infrared spectroscopic observations of the nearby normal Type Ia supernova SN 2021aefx in the nebular phase at $+255$ days past maximum light. Our Near Infrared Spectrograph (NIRSpec) and Mid Infrared Instrument (MIRI) observations, combined with ground-based optical data from the South African Large Telescope (SALT), constitute the first complete optical $+$ NIR $+$ MIR nebular SN Ia spectrum covering 0.3$-$14 $\mu$m. This spectrum unveils the previously unobserved 2.5$-$5 $\mu$m region, revealing strong nebular iron and stable nickel emission, indicative of high-density burning that can constrain the progenitor mass. The data show a significant improvement in sensitivity and resolution compared to previous Spitzer MIR data. We identify numerous NIR and MIR nebular emission lines from iron-group elements and as well as lines from the intermediate-mass element argon. The argon lines extend to higher velocities than the iron-group elements, suggesting stratified ejecta that are a hallmark of delayed-detonation or double-detonation SN Ia models. We present fits to simple geometric line profiles to features beyond 1.2 $\mu$m and find that most lines are consistent with Gaussian or spherical emission distributions, while the [Ar III] 8.99 $\mu$m line has a distinctively flat-topped profile indicating a thick spherical shell of emission. Using our line profile fits, we investigate the emissivity structure of SN 2021aefx and measure kinematic properties. Continued observations of SN 2021aefx and other SNe Ia with JWST will be transformative to the study of SN Ia composition, ionization structure, density, and temperature, and will provide important constraints on SN Ia progenitor and explosion models., Comment: published in ApJ Letters, 17 pages, 12 figures

In this article I recommend a better point estimator for Krippendorff's Alpha agreement coefficient, and develop a jackknife variance estimator that leads to much better interval estimation than does the customary bootstrap procedure or an alternative bootstrap procedure. Having developed the new methodology, I analyze nominal data previously analyzed by Krippendorff, and two experimentally observed datasets: (1) ordinal data from an imaging study of congenital diaphragmatic hernia, and (2) United States Environmental Protection Agency air pollution data for the Philadelphia, Pennsylvania area. The latter two applications are novel. The proposed methodology is now supported in version 2.0 of my open source R package, krippendorffsalpha, which supports common and user-defined distance functions, and can accommodate any number of units, any number of coders, and missingness. Interval computation can be parallelized.

We present our velocity measurements of 59 clumpy, metal-rich ejecta knots in the supernova remnant (SNR) of SN 1572 (Tycho). We use our 450 ks Chandra High Energy Transmission Grating Spectrometer observation to measure the Doppler shift of the He-like Si K$\alpha$ line-center wavelength emitted from these knots to find their line-of-sight (radial) velocities ($v_r$). We find $v_r$ up to $\sim$ 5500 km s$^{-1}$, with roughly consistent speeds between blueshifted and redshifted ejecta knots. We also measure the proper motions (PMs) for our sample based on archival Chandra Advanced CCD Imaging Spectrometer data taken from 2003, 2009, and 2015. We estimate PMs up to 0$"$.35 yr$^{-1}$, which corresponds to a transverse velocity of about 5800 km s$^{-1}$ for the distance of 3.5 kpc to Tycho. Our $v_r$ and transverse velocity measurements imply space velocities of $\sim$ 1900 - 6000 km s$^{-1}$ for the ejecta knots in Tycho. We estimate a new expansion center of R.A.(J2000) = 00$^h$25$^m$18$^s$.725 $\pm$ 1$^s$.157 and decl.(J2000) = +64$^{\circ}$08$'$02$"$.5 $\pm$ 11$"$.2 from our PM measurements, consistent to within $\sim$ 13$"$ of the geometric center. The distribution of space velocities throughout the remnant suggests that the southeast quadrant generally expands faster than the rest of the SNR. We find that blueshifted knots are projected more in the northern shell, while redshifted knots are more in the southern shell. The previously estimated reverse shock position is consistent with most of our estimated ejecta distribution, however some ejecta show deviations from the 1-D picture of the reverse shock., Comment: 16 pages, 7 figures, 3 tables, accepted for publication in ApJ

Using HST narrow-band H{\alpha} images of supernova remnant 0509-67.5 taken {\sim}10 years apart, we measure the forward shock (FS) proper motions (PMs) at 231 rim locations. The average shock radius and velocity are 3.66 {\pm} 0.036 pc and 6315 {\pm} 310 km s{^{-1}}. Hydrodynamic simulations, recast as similarity solutions, provide models for the remnant's expansion into a uniform ambient medium. These are coupled to a Markov chain Monte Carlo (MCMC) analysis to determine explosion parameters, constrained by the FS measurements. For specific global parameters, the MCMC posterior distributions yield an age of 315.5 {\pm} 1.8 yr, a dynamical explosion center at 5h09m31s.16s -67{^\circ}31{^\prime}17.1{^{\prime\prime}} and ambient medium densities at each azimuth ranging from 3.7-8.0 {\times} 10{^{-25}} g cm{^{-3}}. We can detect stellar PMs corresponding to speeds in the LMC of 770 km s{^{-1}} or more using the H${\alpha}$ images. Five stars in the remnant show measurable PMs but none appear to be moving radially from the center, including a prominent red star 4.6{^{\prime\prime}} from the center. Using coronal [Fe XIV] {\lambda}5303 emission as a proxy for the reverse shock location, we constrain the explosion energy (for a compression factor of 4) to a value of E = (1.30 {\pm} 0.41) {\times} 10{^{51}} erg for the first time from shock kinematics alone. Higher compression factors (7 or more) are strongly disfavored based on multiple criteria, arguing for inefficient particle acceleration in the Balmer shocks of 0509-67.5., Comment: 19 pages, 9 figures, 4 tables

We measure the local correlation between radio emission and Compton-$y$ signal across two galaxy clusters, Abell~399 and Abell~401, using maps from the Low-Frequency Array (LOFAR) and the Atacama Cosmology Telescope (ACT) + \Planck. These datasets allow us to make the first measurement of this kind at $\sim$arcminute resolution. We find that the radio brightness scales as $F_{\mathrm{radio}} \propto y^{1.5}$ for Abell~401 and $F_{\mathrm{radio}} \propto y^{2.8}$ for Abell~399. Furthermore, using \XMM data, we derive a sublinear correlation between radio and X-ray brightness for both the clusters ($F_{\mathrm{radio}} \propto F_{\rm X}^{0.7}$). Finally, we correlate the Compton-$y$ and X-ray data, finding that an isothermal model is consistent with the cluster profiles, $y \propto F_{\rm X}^{0.5}$. By adopting an isothermal--$\beta$ model, we are able, for the first time, to jointly use radio, X-ray, and Compton-$y$ data to estimate the scaling index for the magnetic field profile, $B(r) \propto n_{\mathrm{e}}(r)^{\eta}$ in the injection and re-acceleration scenarios. Applying this model, we find that the combined radio and Compton-$y$ signal exhibits a significantly tighter correlation with the X-ray across the clusters than when the datasets are independently correlated. We find $\eta \sim 0.6{-}0.8$. These results are consistent with the upper limit we derive for the scaling index of the magnetic field using rotation measure values for two radio galaxies in Abell~401. We also measure the radio, Compton-$y$, and X-ray correlations in the filament between the clusters but conclude that deeper data are required for a convincing determination of the correlations in the filament., Comment: 16 pages, 8 figures, Accepted for publication in Monthly Notices of the Royal Astronomical Society Main Journal

In the local Universe, OH megamasers (OHMs) are detected almost exclusively in infrared-luminous galaxies, with a prevalence that increases with IR luminosity, suggesting that they trace gas-rich galaxy mergers. Given the proximity of the rest frequencies of OH and the hyperfine transition of neutral atomic hydrogen (HI), radio surveys to probe the cosmic evolution of HI in galaxies also offer exciting prospects for exploiting OHMs to probe the cosmic history of gas-rich mergers. Using observations for the Looking At the Distant Universe with the MeerKAT Array (LADUMA) deep HI survey, we report the first untargeted detection of an OHM at $z > 0.5$, LADUMA J033046.20$-$275518.1 (nicknamed "Nkalakatha"). The host system, WISEA J033046.26$-$275518.3, is an infrared-luminous radio galaxy whose optical redshift $z \approx 0.52$ confirms the MeerKAT emission line detection as OH at a redshift $z_{\rm OH} = 0.5225 \pm 0.0001$ rather than HI at lower redshift. The detected spectral line has 18.4$\sigma$ peak significance, a width of $459 \pm 59\,{\rm km\,s^{-1}}$, and an integrated luminosity of $(6.31 \pm 0.18\,{\rm [statistical]}\,\pm 0.31\,{\rm [systematic]}) \times 10^3\,L_\odot$, placing it among the most luminous OHMs known. The galaxy's far-infrared luminosity $L_{\rm FIR} = (1.576 \pm 0.013) \times 10^{12}\,L_\odot$ marks it as an ultra-luminous infrared galaxy; its ratio of OH and infrared luminosities is similar to those for lower-redshift OHMs. A comparison between optical and OH redshifts offers a slight indication of an OH outflow. This detection represents the first step towards a systematic exploitation of OHMs as a tracer of galaxy growth at high redshifts., Comment: 10 pages, 4 figures. Accepted to ApJ Letters

Instruction plays a critical role in student success. However, most studies of teachers' effects on student learning focus on face-to-face settings. Some aspects of online courses, such as the choice of synchronous or asynchronous instruction and the course structure, might reduce teachers' influence on their students. This study of an online high school biology course offered by Florida Virtual School examined the variation in students' course segment (similar to semester) completion rates, end-of-segment exam scores, and time to segment completion that is attributable to the influence of teachers. Students' scores on the end-of-segment exam varied only slightly across teachers, while completion rates and time to completion varied more. [For the appendixes, see ED609865; for the study snapshot; see ED609866.]

Recent multi-dimensional simulations suggest that high-entropy buoyant plumes help massive stars to explode. Outwardly protruding iron-rich fingers in the galactic supernova remnant Cassiopeia A are uniquely suggestive of this picture. Detecting signatures of specific elements synthesized in the high-entropy nuclear burning regime (i.e., $\alpha$-rich freeze out) would be among the strongest substantiating evidence. Here we report the discovery of such elements, stable Ti and Cr, at a confidence level greater than 5$\sigma$ in the shocked high-velocity iron-rich ejecta of Cassiopeia A. We found the observed Ti/Fe and Cr/Fe mass ratios require $\alpha$-rich freeze out, providing the first observational demonstration for the existence of high-entropy ejecta plumes that boosted the shock wave at explosion. The metal composition of the plumes agrees well with predictions for strongly neutrino-processed proton-rich ejecta. These results support the operation of the convective supernova engine via neutrino heating in the supernova that produced Cassiopeia A., Comment: 34 pages, 13 figures

Models with intractable normalizing functions have numerous applications. Because the normalizing constants are functions of the parameters of interest, standard Markov chain Monte Carlo cannot be used for Bayesian inference for these models. A number of algorithms have been developed for such models. Some have the posterior distribution as their asymptotic distribution. Other ``asymptotically inexact'' algorithms do not possess this property. There is limited guidance for evaluating approximations based on these algorithms. Hence it is very hard to tune them. We propose two new diagnostics that address these problems for intractable normalizing function models. Our first diagnostic, inspired by the second Bartlett identity, is in principle broadly applicable to Monte Carlo approximations beyond the normalizing function problem. We develop an approximate version of this diagnostic that is applicable to intractable normalizing function problems. Our second diagnostic is a Monte Carlo approximation to a kernel Stein discrepancy-based diagnostic introduced by Gorham and Mackey (2017). We provide theoretical justification for our methods and apply them to several algorithms in challenging simulated and real data examples including an Ising model, an exponential random graph model, and a Conway--Maxwell--Poisson regression model, obtaining interesting insights about the algorithms in these contexts.

Compact sources can cause scatter in the scaling relationships between the amplitude of the thermal Sunyaev-Zel'dovich Effect (tSZE) in galaxy clusters and cluster mass. Estimates of the importance of this scatter vary - largely due to limited data on sources in clusters at the frequencies at which tSZE cluster surveys operate. In this paper we present 90 GHz compact source measurements from a sample of 30 clusters observed using the MUSTANG2 instrument on the Green Bank Telescope. We present simulations of how a source's flux density, spectral index, and angular separation from the cluster's center affect the measured tSZE in clusters detected by the Atacama Cosmology Telescope (ACT). By comparing the MUSTANG2 measurements with these simulations we calibrate an empirical relationship between 1.4 GHz flux densities from radio surveys and source contamination in ACT tSZE measurements. We find 3 per cent of the ACT clusters have more than a 20 per cent decrease in Compton-y but another 3 per cent have a 10 per cent increase in the Compton-y due to the matched filters used to find clusters. As sources affect the measured tSZE signal and hence the likelihood that a cluster will be detected, testing the level of source contamination in the tSZE signal using a tSZE selected catalog is inherently biased. We confirm this by comparing the ACT tSZE catalog with optically and X-ray selected cluster catalogs. There is a strong case for a large, high resolution survey of clusters to better characterize their source population., Comment: 13 Pages, 10 figures, 2 tables, with 4 pages of online only figures at end. Published on-line in MNRAS on 22/9/2021

We report a significant detection of the hot intergalactic medium in the filamentary bridge connecting the galaxy clusters Abell 399 and Abell 401. This result is enabled by a low-noise, high-resolution map of the thermal Sunyaev-Zeldovich signal from the Atacama Cosmology Telescope (ACT) and Planck satellite. The ACT data provide the $1.65'$ resolution that allows us to clearly separate the profiles of the clusters, whose centres are separated by $37'$, from the gas associated with the filament. A model that fits for only the two clusters is ruled out compared to one that includes a bridge component at $>5\sigma$. Using a gas temperature determined from Suzaku X-ray data, we infer a total mass of $(3.3\pm0.7)\times10^{14}\,\mathrm{M}_{\odot}$ associated with the filament, comprising about $8\%$ of the entire Abell 399-Abell 401 system. We fit two phenomenological models to the filamentary structure; the favoured model has a width transverse to the axis joining the clusters of ${\sim}1.9\,\mathrm{Mpc}$. When combined with the Suzaku data, we find a gas density of $(0.88\pm0.24)\times10^{-4}\,\mathrm{cm}^{-3}$, considerably lower than previously reported. We show that this can be fully explained by a geometry in which the axis joining Abell 399 and Abell 401 has a large component along the line of sight, such that the distance between the clusters is significantly greater than the $3.2\,\mathrm{Mpc}$ projected separation on the plane of the sky. Finally, we present initial results from higher resolution ($12.7"$ effective) imaging of the bridge with the MUSTANG-2 receiver on the Green Bank Telescope., Comment: 23 pages, 10 figures, 3 tables. This is a pre-copyedited, author-produced PDF of an article accepted for publication in the Monthly Notices of the Royal Astronomical Society following peer review. The version of record is available online at: https://academic.oup.com/mnras/advance-article/doi/10.1093/mnras/stab3391/6442294

We present an improved weak-lensing (WL) study of the high$-z$ $(z=0.87)$ merging galaxy cluster ACT-CL J0102-4915 (El Gordo) based on new wide-field Hubble Space Telescope (HST) imaging data. The new imaging data cover the 3.5$\times$3.5 Mpc region centered on the cluster and enable us to detect WL signals beyond the virial radius, which was not possible in previous studies. We confirm the binary mass structure consisting of the northwestern (NW) and southeastern (SE) subclusters and the 2$\sigma$ dissociation between the SE mass peak and the X-ray cool core. We obtain the mass estimates of the subclusters by simultaneously fitting two Navarro-Frenk-White (NFW) halos without employing mass-concentration relations. The masses are $M_{200c}^{NW} = 9.9^{+2.1}_{-2.2} \times 10^{14} M_{\sun}$ and $M_{200c}^{SE} = 6.5^{+1.9}_{-1.4} \times 10^{14} M_{\sun}$ for the NW and SE subclusters, respectively. The mass ratio is consistent with our previous WL study but significantly different from the previous strong lensing results. This discrepancy is attributed to the use of extrapolation in strong lensing studies because the SE component possesses a higher concentration. By superposing the two best-fit NFW halos, we determine the total mass of El Gordo to be $M_{200c} = 2.13^{+0.25}_{-0.23} \times 10^{15} M_{\sun}$, which is 23% lower than our previous WL result [$M_{200c} =(2.76\pm0.51) \times 10^{15} M_{\sun}$]. Our updated mass is a more direct measurement since we are not extrapolating to $R_{200c}$ as in all previous studies. The new mass is compatible with the current $\Lambda$CDM cosmology., Comment: 27 pages, 18 figures, 3 tables; matches the version accepted for publication in ApJ

We present arcminute-resolution intensity and polarization maps of the Galactic center made with the Atacama Cosmology Telescope (ACT). The maps cover a 32 deg$^2$ field at 98, 150, and 224 GHz with $\vert l\vert\le4^\circ$, $\vert b\vert\le2^\circ$. We combine these data with Planck observations at similar frequencies to create coadded maps with increased sensitivity at large angular scales. With the coadded maps, we are able to resolve many known features of the Central Molecular Zone (CMZ) in both total intensity and polarization. We map the orientation of the plane-of-sky component of the Galactic magnetic field inferred from the polarization angle in the CMZ, finding significant changes in morphology in the three frequency bands as the underlying dominant emission mechanism changes from synchrotron to dust emission. Selected Galactic center sources, including Sgr A*, the Brick molecular cloud (G0.253+0.016), the Mouse pulsar wind nebula (G359.23-0.82), and the Tornado supernova remnant candidate (G357.7-0.1), are examined in detail. These data illustrate the potential for leveraging ground-based Cosmic Microwave Background polarization experiments for Galactic science., Comment: 26 pages, 14 figures, accepted for publication in ApJ

The distribution and kinematics of the circumstellar medium (CSM) around a supernova remnant (SNR) tell us useful information about the explosion of its natal supernova (SN). Kepler's SNR, the remnant of SN1604, is widely regarded to be of Type Ia origin. Its shock is moving through a dense, asymmetric CSM. The presence of this dense gas suggests that its parent progenitor system consisted of a white dwarf and an asymptotic giant branch (AGB) star. In this paper, we analyze a new and long observation with the reflection grating spectrometers (RGS) on board the XMM-Newton satellite, spatially resolving the remnant emission in the cross-dispersion direction. We find that the CSM component is blue-shifted with velocities in the general range 0-500 km/s. We also derive information on the central bar structure and find that the northwest half is blue-shifted, while the southeast half is red-shifted. Our result is consistent with a picture proposed by previous studies, in which a "runaway" AGB star moved to the north-northwest and toward us in the line of sight, although it is acceptable for both single-degenerate and core-degenerate scenarios for the progenitor system., Comment: 19 pages, 5 figures, accepted for publication in ApJ

Galaxy clusters are an important tool for cosmology, and their detection and characterization are key goals for current and future surveys. Using data from the Wide-field Infrared Survey Explorer (WISE), the Massive and Distant Clusters of WISE Survey (MaDCoWS) located 2,839 significant galaxy overdensities at redshifts $0.7\lesssim z\lesssim 1.5$, which included extensive follow-up imaging from the Spitzer Space Telescope to determine cluster richnesses. Concurrently, the Atacama Cosmology Telescope (ACT) has produced large area mm-wave maps in three frequency bands along with a large catalog of Sunyaev-Zeldovich (SZ) selected clusters, as part of its Data Release 5 (DR5). Using the maps and cluster catalog from DR5, we explore the scaling between SZ mass and cluster richness. We use complementary radio survey data from the Very Large Array, submillimeter data from Herschel, and ACT 224~GHz data to assess the impact of contaminating sources on the SZ signals. We then use a hierarchical Bayesian model to fit the mass-richness scaling relation. We find that MaDCoWS clusters have submillimeter contamination which is consistent with a gray-body spectrum, while the ACT clusters are consistent with no submillimeter emission on average. We find the best fit ACT SZ mass vs. MaDCoWS richness scaling relation has a slope of $\kappa = 1.84^{+0.15}_{-0.14}$, where the slope is defined as $M\propto \lambda_{15}^{\kappa}$ where $\lambda_{15}$ is the richness. Additionally, we find that the approximate level of in-fill of the ACT and MaDCoWS cluster SZ signals to be at the percent level, Comment: 25 pages, 17 Figures; accepted for publication in A&A

R package krippendorffsalpha provides tools for measuring agreement using Krippendorff's Alpha coefficient, a well-known nonparametric measure of agreement (also called inter-rater reliability and various other names). This article first develops Krippendorff's Alpha in a natural way, and situates Alpha among statistical procedures. Then the usage of package krippendorffsalpha is illustrated via analyses of two datasets, the latter of which was collected during an imaging study of hip cartilage. The package permits users to apply the Alpha methodology using built-in distance functions for the nominal, ordinal, interval, or ratio levels of measurement. User-defined distance functions are also supported. The fitting function can accommodate any number of units, any number of coders, and missingness. Bootstrap inference is supported, and the bootstrap computation can be carried out in parallel.

The direct Gaussian copula model with discrete marginal distributions is an appealing data-analytic tool but poses difficult computational challenges due to its intractable likelihood. A number of approximations/surrogates for the likelihood have been proposed, including the continuous extension-based approximation (CE) and the distributional transform-based approximation (DT). The continuous extension approach is exact up to Monte Carlo error but does not scale well computationally. The distributional transform approach permits efficient computation but offers no theoretical guarantee that it is exact. In practice, though, the distributional transform-based approximate likelihood is so very nearly exact for some variants of the model as to permit genuine maximum likelihood or Bayesian inference. We demonstrate the exactness of the distributional transform-based objective function for two interesting variants of the model, and propose a quantity that can be used to assess exactness for experimentally observed datasets. Said diagnostic will permit practitioners to determine whether genuine Bayesian inference or ordinary maximum likelihood inference using the DT-based likelihood is possible for a given dataset.

Two recent large data releases for the Atacama Cosmology Telescope (ACT), called DR4 and DR5, are available for public access. These data include temperature and polarization maps that cover nearly half the sky at arcminute resolution in three frequency bands; lensing maps and component-separated maps covering ~ 2,100 deg^2 of sky; derived power spectra and cosmological likelihoods; a catalog of over 4,000 galaxy clusters; and supporting ancillary products including beam functions and masks. The data and products are described in a suite of ACT papers; here we provide a summary. In order to facilitate ease of access to these data we present a set of Jupyter IPython notebooks developed to introduce users to DR4, DR5, and the tools needed to analyze these data. The data products (excluding simulations) and the set of notebooks are publicly available on the NASA Legacy Archive for Microwave Background Data Analysis (LAMBDA); simulation products are available on the National Energy Research Scientific Computing Center (NERSC)., Comment: Accepted to ApJS. 21 pages, 8 figures. Data and notebooks available on LAMBDA https://lambda.gsfc.nasa.gov/product/act/

We present a 5.4$\sigma$ detection of the pairwise kinematic Sunyaev-Zel'dovich (kSZ) effect using Atacama Cosmology Telescope (ACT) and $\it{Planck}$ CMB observations in combination with Luminous Red Galaxy samples from the Sloan Digital Sky Survey (SDSS) DR15 catalog. Results are obtained using three ACT CMB maps: co-added 150 GHz and 98 GHz maps, combining observations from 2008-2018 (ACT DR5), which overlap with SDSS DR15 over 3,700 sq. deg., and a component-separated map using night-time only observations from 2014-2015 (ACT DR4), overlapping with SDSS DR15 over 2,089 sq. deg. Comparisons of the results from these three maps provide consistency checks in relation to potential frequency-dependent foreground contamination. A total of 343,647 galaxies are used as tracers to identify and locate galaxy groups and clusters from which the kSZ signal is extracted using aperture photometry. We consider the impact of various aperture photometry assumptions and covariance estimation methods on the signal extraction. Theoretical predictions of the pairwise velocities are used to obtain best-fit, mass-averaged, optical depth estimates for each of five luminosity-selected tracer samples. A comparison of the kSZ-derived optical depth measurements obtained here to those derived from the thermal SZ effect for the same sample is presented in a companion paper., Comment: 17 pages, 10 figures. Updated to match published version in PRD

We present high signal-to-noise measurements (up to 12$\sigma$) of the average thermal Sunyaev Zel'dovich (tSZ) effect from optically selected galaxy groups and clusters and estimate their baryon content within a 2.1$^\prime$ radius aperture. Sources from the Sloan Digital Sky Survey (SDSS) Baryon Oscillation Spectroscopic Survey (BOSS) DR15 catalog overlap with 3,700 sq. deg. of sky observed by the Atacama Cosmology Telescope (ACT) from 2008 to 2018 at 150 and 98 GHz (ACT DR5), and 2,089 sq. deg. of internal linear combination component-separated maps combining ACT and $\it{Planck}$ data (ACT DR4). The corresponding optical depths, $\bar{\tau}$, which depend on the baryon content of the halos, are estimated using results from cosmological hydrodynamic simulations assuming an AGN feedback radiative cooling model. We estimate the mean mass of the halos in multiple luminosity bins, and compare the tSZ-based $\bar{\tau}$ estimates to theoretical predictions of the baryon content for a Navarro-Frenk-White profile. We do the same for $\bar{\tau}$ estimates extracted from fits to pairwise baryon momentum measurements of the kinematic Sunyaev-Zel'dovich effect (kSZ) for the same data set obtained in a companion paper. We find that the $\bar{\tau}$ estimates from the tSZ measurements in this work and the kSZ measurements in the companion paper agree within $1\sigma$ for two out of the three disjoint luminosity bins studied, while they differ by 2-3$\sigma$ in the highest luminosity bin. The optical depth estimates account for one third to all of the theoretically predicted baryon content in the halos across luminosity bins. Potential systematic uncertainties are discussed. The tSZ and kSZ measurements provide a step towards empirical Compton-$\bar{y}$-$\bar{\tau}$ relationships to provide new tests of cluster formation and evolution models., Comment: 19 pages, 9 figures

We use gravitational lensing of the cosmic microwave background (CMB) to measure the mass of the most distant blindly-selected sample of galaxy clusters on which a lensing measurement has been performed to date. In CMB data from the the Atacama Cosmology Telescope (ACT) and the Planck satellite, we detect the stacked lensing effect from 677 near-infrared-selected galaxy clusters from the Massive and Distant Clusters of WISE Survey (MaDCoWS), which have a mean redshift of $ \langle z \rangle = 1.08$. There are no current optical weak lensing measurements of clusters that match the distance and average mass of this sample. We detect the lensing signal with a significance of $4.2 \sigma$. We model the signal with a halo model framework to find the mean mass of the population from which these clusters are drawn. Assuming that the clusters follow Navarro-Frenk-White density profiles, we infer a mean mass of $\langle M_{500c}\rangle = \left(1.7 \pm 0.4 \right)\times10^{14}\,\mathrm{M}_\odot$. We consider systematic uncertainties from cluster redshift errors, centering errors, and the shape of the NFW profile. These are all smaller than 30% of our reported uncertainty. This work highlights the potential of CMB lensing to enable cosmological constraints from the abundance of distant clusters populating ever larger volumes of the observable Universe, beyond the capabilities of optical weak lensing measurements., Comment: 14 pages, 3 figures, matches version accepted in ApJL, code available at https://github.com/ACTCollaboration/madcows_lensing/

The thermal and kinematic Sunyaev-Zel'dovich effects (tSZ, kSZ) probe the thermodynamic properties of the circumgalactic and intracluster medium (CGM and ICM) of galaxies, groups, and clusters, since they are proportional, respectively, to the integrated electron pressure and momentum along the line-of-sight. We present constraints on the gas thermodynamics of CMASS galaxies in the Baryon Oscillation Spectroscopic Survey (BOSS) using new measurements of the kSZ and tSZ signals obtained in a companion paper. Combining kSZ and tSZ measurements, we measure within our model the amplitude of energy injection $\epsilon M_\star c^2$, where $M_\star$ is the stellar mass, to be $\epsilon=(40\pm9)\times10^{-6}$, and the amplitude of the non-thermal pressure profile to be $\alpha_{\rm Nth}<0.2$ (2$\sigma$), indicating that less than 20% of the total pressure within the virial radius is due to a non-thermal component. We estimate the effects of including baryons in the modeling of weak-lensing galaxy cross-correlation measurements using the best-fit density profile from the kSZ measurement. Our estimate reduces the difference between the original theoretical model and the weak-lensing galaxy cross-correlation measurements in arXiv:1611.08606 by half but does not fully reconcile it. Comparing the tSZ measurements to cosmological simulations, we find that simulations underestimate the CGM pressure at large radii while they fare better in comparison with the kSZ measurements. This suggests that the energy injected via feedback models in the simulations that we compared against does not sufficiently heat the gas at these radii. We do not find significant disagreement at smaller radii. These measurements provide novel tests of current and future simulations. This work demonstrates the power of joint, high signal-to-noise kSZ and tSZ observations, upon which future cross-correlation studies will improve., Comment: Corrected error in the algorithm that calculates the kSZ temperature profile for a given GNFW density model. The value of $\log_{\rm10} \rho_0$ changed by 0.75$\sigma$ in Tab.II and Fig.2, and affected the results in the left panels of Fig. 6. Conclusions are unchanged. Erratum published at https://link.aps.org/doi/10.1103/PhysRevD.107.049903

The scattering of cosmic microwave background (CMB) photons off the free-electron gas in galaxies and clusters leaves detectable imprints on high resolution CMB maps: the thermal and kinematic Sunyaev-Zel'dovich effects (tSZ and kSZ respectively). We use combined microwave maps from the Atacama Cosmology Telescope (ACT) DR5 and Planck in combination with the CMASS and LOWZ galaxy catalogs from the Baryon Oscillation Spectroscopic Survey (BOSS DR10 and DR12), to study the gas associated with these galaxy groups. Using individual reconstructed velocities, we perform a stacking analysis and reject the no-kSZ hypothesis at 6.5$\sigma$, the highest significance to date. This directly translates into a measurement of the electron number density profile, and thus of the gas density profile. Despite the limited signal to noise, the measurement shows at high significance that the gas density profile is more extended than the dark matter density profile, for any reasonable baryon abundance (formally $>90\sigma$ for the cosmic baryon abundance). We simultaneously measure the tSZ signal, i.e. the electron thermal pressure profile of the same CMASS objects, and reject the no-tSZ hypothesis at 10$\sigma$. We combine tSZ and kSZ measurements to estimate the electron temperature to 20% precision in several aperture bins, and find it comparable to the virial temperature. In a companion paper, we analyze these measurements to constrain the gas thermodynamics and the properties of feedback inside galaxy groups. We present the corresponding LOWZ measurements in this paper, ruling out a null kSZ (tSZ) signal at 2.9 (13.9)$\sigma$, and leave their interpretation to future work. Our stacking software ThumbStack is publicly available at https://github.com/EmmanuelSchaan/ThumbStack and directly applicable to future Simons Observatory and CMB-S4 data., Comment: Accepted in Physical Review D, Editors' Suggestion

We measure the projected number density profiles of galaxies and the splashback feature in clusters selected by the Sunyaev--Zeldovich (SZ) effect from the Advanced Atacama Cosmology Telescope (AdvACT) survey using galaxies observed by the Dark Energy Survey (DES). The splashback radius for the complete galaxy sample is consistent with theoretical measurements from CDM-only simulations, and is located at $2.4^{+0.3}_{-0.4}$ Mpc $h^{-1}$. We split the sample based on galaxy color and find significant differences in the profile shapes. Red galaxies and those in the green valley show a splashback-like minimum in their slope profile consistent with theoretical predictions, while the bluest galaxies show a weak feature that appears at a smaller radius. We develop a mapping of galaxies to subhalos in $N$-body simulations by splitting subhalos based on infall time onto the cluster halos. We find that the location of the steepest slope and differences in the shapes of the profiles can be mapped to differences in the average time of infall of galaxies of different colors. The minima of the slope in the galaxy profiles trace a discontinuity in the phase space of dark matter halos. By relating spatial profiles to infall time for galaxies of different colours, we can use splashback as a clock to understand galaxy quenching. We find that red galaxies have on average been in their clusters for over $3.2 ~\rm Gyrs$, green galaxies about $2.2 ~\rm Gyrs$, while blue galaxies have been accreted most recently and have not reached apocenter. Using the information from the complete radial profiles, we fit a simple quenching model and find that the onset of galaxy quenching in clusters occurs after a delay of about a gigayear, and that galaxies quench rapidly thereafter with an exponential timescale of $0.6$ Gyr., Comment: 23 pages, 10 figures, to be submitted to ApJ

We present LCDM cosmological parameter constraints obtained from delensed microwave background power spectra. Lensing maps from a subset of DR4 data from the Atacama Cosmology Telescope (ACT) are used to undo the lensing effect in ACT spectra observed at 150 and 98 GHz. At 150 GHz, we remove the lensing distortion with an effective efficiency of 30% (TT), 30% (EE), 26% (TE) and 20% (BB); this results in detections of the delensing effect at 8.7 sigma (TT), 5.1 sigma (EE), 2.6 sigma (TE), and 2.4 sigma (BB) significance. The combination of 150 and 98 GHz TT, EE, and TE delensed spectra is well fit by a standard LCDM model. We also measure the shift in best-fit parameters when fitting delensed versus lensed spectra; while this shift does not inform our ability to measure cosmological parameters, it does provide a three-way consistency check among the lensing inferred from the best-fit parameters, the lensing in the CMB power spectrum, and the reconstructed lensing map. This shift is predicted to be zero when fitting with the correct model since both lensed and delensed spectra originate from the same region of sky. Fitting with a LCDM model and marginalizing over foregrounds, we find that the shift in cosmological parameters is consistent with zero. Our results show that gravitational lensing of the microwave background is internally consistent within the framework of the standard cosmological model., Comment: 29 pages, 17 figures, version matches that accepted by JCAP

This paper presents a maximum-likelihood algorithm for combining sky maps with disparate sky coverage, angular resolution and spatially varying anisotropic noise into a single map of the sky. We use this to merge hundreds of individual maps covering the 2008-2018 ACT observing seasons, resulting in by far the deepest ACT maps released so far. We also combine the maps with the full Planck maps, resulting in maps that have the best features of both Planck and ACT: Planck's nearly white noise on intermediate and large angular scales and ACT's high-resolution and sensitivity on small angular scales. The maps cover over 18,000 square degrees, nearly half the full sky, at 100, 150 and 220 GHz. They reveal 4,000 optically-confirmed clusters through the Sunyaev Zel'dovich effect (SZ) and 18,500 point source candidates at $> 5\sigma$, the largest single collection of SZ clusters and millimeter wave sources to date. The multi-frequency maps provide millimeter images of nearby galaxies and individual Milky Way nebulae, and even clear detections of several nearby stars. Other anticipated uses of these maps include, for example, thermal SZ and kinematic SZ cluster stacking, CMB cluster lensing and galactic dust science. The method itself has negligible bias. However, due to the preliminary nature of some of the component data sets, we caution that these maps should not be used for precision cosmological analysis. The maps are part of ACT DR5, and are available on LAMBDA at https://lambda.gsfc.nasa.gov/product/act/actpol_prod_table.cfm. There is also a web atlas at https://phy-act1.princeton.edu/public/snaess/actpol/dr5/atlas., Comment: 38 pages, 29 figures, data release on lambda. Published in JCAP

We present new arcminute-resolution maps of the Cosmic Microwave Background temperature and polarization anisotropy from the Atacama Cosmology Telescope, using data taken from 2013-2016 at 98 and 150 GHz. The maps cover more than 17,000 deg$^2$, the deepest 600 deg$^2$ with noise levels below 10 $\mu$K-arcmin. We use the power spectrum derived from almost 6,000 deg$^2$ of these maps to constrain cosmology. The ACT data enable a measurement of the angular scale of features in both the divergence-like polarization and the temperature anisotropy, tracing both the velocity and density at last-scattering. From these one can derive the distance to the last-scattering surface and thus infer the local expansion rate, $H_0$. By combining ACT data with large-scale information from WMAP we measure $H_0 = 67.6 \pm 1.1$ km/s/Mpc, at 68% confidence, in excellent agreement with the independently-measured Planck satellite estimate (from ACT alone we find $H_0 = 67.9 \pm 1.5$ km/s/Mpc). The $\Lambda$CDM model provides a good fit to the ACT data, and we find no evidence for deviations: both the spatial curvature, and the departure from the standard lensing signal in the spectrum, are zero to within 1$\sigma$; the number of relativistic species, the primordial Helium fraction, and the running of the spectral index are consistent with $\Lambda$CDM predictions to within $1.5 - 2.2\sigma$. We compare ACT, WMAP, and Planck at the parameter level and find good consistency; we investigate how the constraints on the correlated spectral index and baryon density parameters readjust when adding CMB large-scale information that ACT does not measure. The DR4 products presented here will be publicly released on the NASA Legacy Archive for Microwave Background Data Analysis., Comment: 33 pages, 24 figures, products available on the NASA LAMBDA website, version accepted for publication in JCAP