-
A More Precise Measurement of the Radius of PSR J0740+6620 Using Updated NICER Data
Authors:
Alexander J. Dittmann,
M. Coleman Miller,
Frederick K. Lamb,
Isiah Holt,
Cecilia Chirenti,
Michael T. Wolff,
Slavko Bogdanov,
Sebastien Guillot,
Wynn C. G. Ho,
Sharon M. Morsink,
Zaven Arzoumanian,
Keith C. Gendreau
Abstract:
PSR J0740+6620 is the neutron star with the highest precisely determined mass, inferred from radio observations to be $2.08\pm0.07\,\rm M_\odot$. Measurements of its radius therefore hold promise to constrain the properties of the cold, catalyzed, high-density matter in neutron star cores. Previously, Miller et al. (2021) and Riley et al. (2021) reported measurements of the radius of PSR J0740+662…
▽ More
PSR J0740+6620 is the neutron star with the highest precisely determined mass, inferred from radio observations to be $2.08\pm0.07\,\rm M_\odot$. Measurements of its radius therefore hold promise to constrain the properties of the cold, catalyzed, high-density matter in neutron star cores. Previously, Miller et al. (2021) and Riley et al. (2021) reported measurements of the radius of PSR J0740+6620 based on Neutron Star Interior Composition Explorer (NICER) observations accumulated through 17 April 2020, and an exploratory analysis utilizing NICER background estimates and a data set accumulated through 28 December 2021 was presented in Salmi et al. (2022). Here we report an updated radius measurement, derived by fitting models of X-ray emission from the neutron star surface to NICER data accumulated through 21 April 2022, totaling $\sim1.1$ Ms additional exposure compared to the data set analyzed in Miller et al. (2021) and Riley et al. (2021), and to data from X-ray Multi-Mirror (XMM-Newton) observations. We find that the equatorial circumferential radius of PSR J0740+6620 is $12.92_{-1.13}^{+2.09}$ km (68% credibility), a fractional uncertainty $\sim83\%$ the width of that reported in Miller et al. (2021), in line with statistical expectations given the additional data. If we were to require the radius to be less than 16 km, as was done in Salmi et al. (2024), then our 68% credible region would become $R=12.76^{+1.49}_{-1.02}$ km, which is close to the headline result of Salmi et al. (2024). Our updated measurements, along with other laboratory and astrophysical constraints, imply a slightly softer equation of state than that inferred from our previous measurements.
△ Less
Submitted 30 June, 2024; v1 submitted 20 June, 2024;
originally announced June 2024.
-
Using graph neural networks to reconstruct charged pion showers in the CMS High Granularity Calorimeter
Authors:
M. Aamir,
B. Acar,
G. Adamov,
T. Adams,
C. Adloff,
S. Afanasiev,
C. Agrawal,
C. Agrawal,
A. Ahmad,
H. A. Ahmed,
S. Akbar,
N. Akchurin,
B. Akgul,
B. Akgun,
R. O. Akpinar,
E. Aktas,
A. AlKadhim,
V. Alexakhin,
J. Alimena,
J. Alison,
A. Alpana,
W. Alshehri,
P. Alvarez Dominguez,
M. Alyari,
C. Amendola
, et al. (550 additional authors not shown)
Abstract:
A novel method to reconstruct the energy of hadronic showers in the CMS High Granularity Calorimeter (HGCAL) is presented. The HGCAL is a sampling calorimeter with very fine transverse and longitudinal granularity. The active media are silicon sensors and scintillator tiles readout by SiPMs and the absorbers are a combination of lead and Cu/CuW in the electromagnetic section, and steel in the hadr…
▽ More
A novel method to reconstruct the energy of hadronic showers in the CMS High Granularity Calorimeter (HGCAL) is presented. The HGCAL is a sampling calorimeter with very fine transverse and longitudinal granularity. The active media are silicon sensors and scintillator tiles readout by SiPMs and the absorbers are a combination of lead and Cu/CuW in the electromagnetic section, and steel in the hadronic section. The shower reconstruction method is based on graph neural networks and it makes use of a dynamic reduction network architecture. It is shown that the algorithm is able to capture and mitigate the main effects that normally hinder the reconstruction of hadronic showers using classical reconstruction methods, by compensating for fluctuations in the multiplicity, energy, and spatial distributions of the shower's constituents. The performance of the algorithm is evaluated using test beam data collected in 2018 prototype of the CMS HGCAL accompanied by a section of the CALICE AHCAL prototype. The capability of the method to mitigate the impact of energy leakage from the calorimeter is also demonstrated.
△ Less
Submitted 30 June, 2024; v1 submitted 17 June, 2024;
originally announced June 2024.
-
Exploring Waveform Variations among Neutron Star Ray-tracing Codes for Complex Emission Geometries
Authors:
Devarshi Choudhury,
Anna L. Watts,
Alexander J. Dittmann,
M. Coleman Miller,
Sharon M. Morsink,
Tuomo Salmi,
Serena Vinciguerra,
Slavko Bogdanov,
Sebastien Guillot,
Michael T. Wolff,
Zaven Arzoumanian
Abstract:
Pulse Profile Modeling (PPM), the technique used to infer mass, radius and geometric parameters for rotation-powered millisecond pulsars using data from the Neutron Star Interior Composition Explorer (NICER), relies on relativistic ray-tracing of thermal X-ray photons from hot spots on the neutron star surface to the observer. To verify our ray-tracing codes we have in the past conducted cross-tes…
▽ More
Pulse Profile Modeling (PPM), the technique used to infer mass, radius and geometric parameters for rotation-powered millisecond pulsars using data from the Neutron Star Interior Composition Explorer (NICER), relies on relativistic ray-tracing of thermal X-ray photons from hot spots on the neutron star surface to the observer. To verify our ray-tracing codes we have in the past conducted cross-tests for simple hot spot geometries, focusing primarily on the implementation of the space-time model. In this paper, we present verification for test problems that explore the more complex hot spot geometries that are now being employed in the NICER PPM analyses. We conclude that the accuracy of our computed waveforms is in general sufficiently high for analyses of current NICER data sets. We have however identified some extreme configurations where extra care may be needed.
△ Less
Submitted 11 June, 2024;
originally announced June 2024.
-
The Effects of Cooling on Boundary Layer Accretion
Authors:
Alexander J. Dittmann
Abstract:
In many cases accretion proceeds from disks onto planets, stars, white dwarfs, and neutron stars via a boundary layer, a region of intense shear where gas transitions from a near-Keplerian speed to that of the surface. These regions are \textit{not} susceptible to the common magnetorotational and Kelvin-Helmholtz instabilities, and instead global modes generated by supersonic shear instabilities a…
▽ More
In many cases accretion proceeds from disks onto planets, stars, white dwarfs, and neutron stars via a boundary layer, a region of intense shear where gas transitions from a near-Keplerian speed to that of the surface. These regions are \textit{not} susceptible to the common magnetorotational and Kelvin-Helmholtz instabilities, and instead global modes generated by supersonic shear instabilities are a leading candidate to govern transport in these regions. This work investigates the dynamics of these systems under a range of thermodynamic conditions, surveying both disk sound speeds and cooling rates. Very fast or very slow cooling has little effect on wave dynamics: in the fast-cooling limit, waves propagate in an effectively isothermal manner, and in the slow limit wave propagation is effectively adiabatic. However, when the cooling timescale is comparable to the wave period, wave damping becomes extreme. In cases with intermediate cooling rates, mass and angular momentum transport can be suppressed by orders of magnitude compared to isothermal and uncooled cases. Cooling in accretion disks leads to a preference for wavenumbers near and below the Mach number of the disk; the corresponding lower frequencies can (in non-isothermal systems) couple to gravity modes within the star, potentially causing low-frequency variability such as dwarf nova and quasi-periodic oscillations in accreting systems.
△ Less
Submitted 30 May, 2024;
originally announced May 2024.
-
Notes on the Practical Application of Nested Sampling: MultiNest, (Non)convergence, and Rectification
Authors:
Alexander J. Dittmann
Abstract:
Nested sampling is a promising tool for Bayesian statistical analysis because it simultaneously performs parameter estimation and facilitates model comparison. MultiNest is one of the most popular nested sampling implementations, and has been applied to a wide variety of problems in the physical sciences. However, MultiNest results are frequently unreliable, and accompanying convergence tests are…
▽ More
Nested sampling is a promising tool for Bayesian statistical analysis because it simultaneously performs parameter estimation and facilitates model comparison. MultiNest is one of the most popular nested sampling implementations, and has been applied to a wide variety of problems in the physical sciences. However, MultiNest results are frequently unreliable, and accompanying convergence tests are a necessary component of any analysis. Using simple, analytically tractable test problems, I illustrate how MultiNest (1) can produce systematically biased estimates of the Bayesian evidence, which are more significantly biased for problems of higher dimensionality; (2) can derive posterior estimates with errors on the order of $\sim100\%$; (3) is more likely to underestimate the width of a credible interval than to overestimate it - to a minor degree for smooth problems, but much more so when sampling noisy likelihoods. Nevertheless, I show how MultiNest can be used to jump-start Markov chain Monte Carlo sampling or more rigorous nested sampling techniques, potentially accelerating more trustworthy measurements of posterior distributions and Bayesian evidences, and overcoming the challenge of Markov chain Monte Carlo initialization.
△ Less
Submitted 25 April, 2024;
originally announced April 2024.
-
Lightweight Multi-System Multivariate Interconnection and Divergence Discovery
Authors:
Mulugeta Weldezgina Asres,
Christian Walter Omlin,
Jay Dittmann,
Pavel Parygin,
Joshua Hiltbrand,
Seth I. Cooper,
Grace Cummings,
David Yu
Abstract:
Identifying outlier behavior among sensors and subsystems is essential for discovering faults and facilitating diagnostics in large systems. At the same time, exploring large systems with numerous multivariate data sets is challenging. This study presents a lightweight interconnection and divergence discovery mechanism (LIDD) to identify abnormal behavior in multi-system environments. The approach…
▽ More
Identifying outlier behavior among sensors and subsystems is essential for discovering faults and facilitating diagnostics in large systems. At the same time, exploring large systems with numerous multivariate data sets is challenging. This study presents a lightweight interconnection and divergence discovery mechanism (LIDD) to identify abnormal behavior in multi-system environments. The approach employs a multivariate analysis technique that first estimates the similarity heatmaps among the sensors for each system and then applies information retrieval algorithms to provide relevant multi-level interconnection and discrepancy details. Our experiment on the readout systems of the Hadron Calorimeter of the Compact Muon Solenoid (CMS) experiment at CERN demonstrates the effectiveness of the proposed method. Our approach clusters readout systems and their sensors consistent with the expected calorimeter interconnection configurations, while capturing unusual behavior in divergent clusters and estimating their root causes.
△ Less
Submitted 12 April, 2024;
originally announced April 2024.
-
Single Transit Detection In Kepler With Machine Learning And Onboard Spacecraft Diagnostics
Authors:
Matthew T. Hansen,
Jason A. Dittmann
Abstract:
Exoplanet discovery at long orbital periods requires reliably detecting individual transits without additional information about the system. Techniques like phase-folding of light curves and periodogram analysis of radial velocity data are more sensitive to planets with shorter orbital periods, leaving a dearth of planet discoveries at long periods. We present a novel technique using an ensemble o…
▽ More
Exoplanet discovery at long orbital periods requires reliably detecting individual transits without additional information about the system. Techniques like phase-folding of light curves and periodogram analysis of radial velocity data are more sensitive to planets with shorter orbital periods, leaving a dearth of planet discoveries at long periods. We present a novel technique using an ensemble of Convolutional Neural Networks incorporating the onboard spacecraft diagnostics of \emph{Kepler} to classify transits within a light curve. We create a pipeline to recover the location of individual transits, and the period of the orbiting planet, which maintains $>80\%$ transit recovery sensitivity out to an 800-day orbital period. Our neural network pipeline has the potential to discover additional planets in the \emph{Kepler} dataset, and crucially, within the $η$-Earth regime. We report our first candidate from this pipeline, KOI 1271.02. KOI 1271.01 is known to exhibit strong Transit Timing Variations (TTVs), and so we jointly model the TTVs and transits of both transiting planets to constrain the orbital configuration and planetary parameters and conclude with a series of potential parameters for KOI 1271.02, as there is not enough data currently to uniquely constrain the system. We conclude that KOI 1271.02 has a radius of 5.32 $\pm$ 0.20 $R_{\oplus}$ and a mass of $28.94^{0.23}_{-0.47}$ $M_{\oplus}$. Future constraints on the nature of KOI 1271.02 require measuring additional TTVs of KOI 1271.01 or observing a second transit of KOI 1271.02.
△ Less
Submitted 5 March, 2024;
originally announced March 2024.
-
The Santa Barbara Binary-Disk Code Comparison
Authors:
Paul C. Duffell,
Alexander J. Dittmann,
Daniel J. D'Orazio,
Alessia Franchini,
Kaitlin M. Kratter,
Anna B. T. Penzlin,
Enrico Ragusa,
Magdalena Siwek,
Christopher Tiede,
Haiyang Wang,
Jonathan Zrake,
Adam M. Dempsey,
Zoltan Haiman,
Alessandro Lupi,
Michal Pirog,
Geoffrey Ryan
Abstract:
We have performed numerical calculations of a binary interacting with a gas disk, using eleven different numerical methods and a standard binary-disk setup. The goal of this study is to determine whether all codes agree on a numerically converged solution, and to determine the necessary resolution for convergence and the number of binary orbits that must be computed to reach an agreed-upon relaxed…
▽ More
We have performed numerical calculations of a binary interacting with a gas disk, using eleven different numerical methods and a standard binary-disk setup. The goal of this study is to determine whether all codes agree on a numerically converged solution, and to determine the necessary resolution for convergence and the number of binary orbits that must be computed to reach an agreed-upon relaxed state of the binary-disk system. We find that all codes can agree on a converged solution (depending on the diagnostic being measured). The zone spacing required for most codes to reach a converged measurement of the torques applied to the binary by the disk is roughly 1% of the binary separation in the vicinity of the binary components. For our disk model to reach a relaxed state, codes must be run for at least 200 binary orbits, corresponding to about a viscous time for our parameters, $0.2 (a^2 Ω_B /ν)$ binary orbits, where $ν$ is the kinematic viscosity. We did not investigate dependence on binary mass ratio, eccentricity, disk temperature, or disk viscosity; therefore, these benchmarks may act as guides towards expanding converged solutions to the wider parameter space but might need to be updated in a future study that investigates dependence on system parameters. We find the most major discrepancies between codes resulted from the dimensionality of the setup (3D vs 2D disks). Beyond this, we find good agreement in the total torque on the binary between codes, although the partition of this torque between the gravitational torque, orbital accretion torque, and spin accretion torque depends sensitively on the sink prescriptions employed. In agreement with previous studies, we find a modest difference in torques and accretion variability between 2D and 3D disk models. We find cavity precession rates to be appreciably faster in 3D than in 2D.
△ Less
Submitted 20 June, 2024; v1 submitted 20 February, 2024;
originally announced February 2024.
-
Quantifying the Effect of Short-timescale Stellar Activity Upon Transit Detection in M Dwarfs
Authors:
Dana Clarice Yaptangco,
Sarah Ballard,
Jason Dittmann
Abstract:
M dwarf stars comprise 70-80% of the galaxy's stars and host most of its rocky planets. They also importantly differ from Sunlike stars in that they are "active" for billions of years or more: rotating quickly, flaring often, and emitting large amounts of UV and X-ray light. The effects of stellar activity upon both photometry and spectroscopy make their exoplanets more difficult to detect, and M…
▽ More
M dwarf stars comprise 70-80% of the galaxy's stars and host most of its rocky planets. They also importantly differ from Sunlike stars in that they are "active" for billions of years or more: rotating quickly, flaring often, and emitting large amounts of UV and X-ray light. The effects of stellar activity upon both photometry and spectroscopy make their exoplanets more difficult to detect, and M dwarfs exhibit this behavior for thousands of times longer than a typical Sunlike star. While activity signals such as flaring and stellar rotation can be more readily modeled or removed from photometry, the contribution of unresolved stellar activity to transit sensitivity is harder to quantify. In this paper, we investigate the difference in the detectability of planetary transits around a sample of M dwarfs observed by NASA's TESS Mission, characterized by a common stellar radius, effective temperature, and TESS magnitude. Our sample is classified as either "active" or "inactive" based upon the presence of H$α$ in emission. After removing the more readily identifiable signatures of activity: stellar rotation and large flares, we perform an injection-and-recovery analysis of transits for each star. We extract detection sensitivity as a function of planetary radius and orbital period for each star in the sample. Then, we produce averaged sensitivity maps for the "active" stars and the "inactive" stars, for the sake of comparison. We quantify the extent to which signal-to-noise is degraded for transit detection, when comparing an active star to an inactive star of the same temperature and apparent brightness. We aim for these sensitivity maps to be useful to the exoplanet community in future M dwarf occurrence rate studies.
△ Less
Submitted 31 January, 2024;
originally announced February 2024.
-
Variations in the Radius Distribution of Single and Compact Multiple Transiting Planets
Authors:
Benjamin T. Liberles,
Jason A. Dittmann,
Stephen M. Elardo,
Sarah Ballard
Abstract:
Previous work has established the enhanced occurrence of compact systems of multiple small exoplanets around metal-poor stars. Understanding the origin of this effect in the planet formation process is a topic of ongoing research. Here we consider the radii of planets residing in systems of multiple transiting planets, compared to those residing in single-transiting systems, with a particular focu…
▽ More
Previous work has established the enhanced occurrence of compact systems of multiple small exoplanets around metal-poor stars. Understanding the origin of this effect in the planet formation process is a topic of ongoing research. Here we consider the radii of planets residing in systems of multiple transiting planets, compared to those residing in single-transiting systems, with a particular focus on late-type host stars. We investigate whether the two radius distributions are consistent with being drawn from the same underlying planetary population. We construct a planetary sample of 290 planets around late-K and M-dwarfs containing 149 planets from single-transiting planetary systems and 141 planets from multi-transiting compact multiple planetary systems (54 compact multiples). We performed a two sample Kolmogorov-Smirnov test, Mann-Whitney U test, and Anderson-Darling K-Sampling test on the radius distributions of our two samples. We find statistical evidence (p < 0.0026) that planets in compact multiple systems are larger, on average, than their single-transiting counterparts. We determine that the offset cannot be explained by detection bias. We investigate whether this effect could be explained via more efficient outgassing of a secondary atmosphere in compact multiple systems due to the stress and strain forces of interplanetary tides on planetary interiors. We find that this effect is insufficient to explain our observations without significant enrichment in H2O compared to Earth-like bulk composition.
△ Less
Submitted 10 December, 2023;
originally announced December 2023.
-
Runaway Eccentricity Growth: A Pathway for Binary Black Hole Mergers in AGN Disks
Authors:
Josh Calcino,
Adam M. Dempsey,
Alexander J. Dittmann,
Hui Li
Abstract:
Binary black holes embedded within the accretion disks that fuel active galactic nuclei (AGN) are promising progenitors for the source of gravitational wave events detected by LIGO/VIRGO. Several recent studies have shown that when these binaries form they should be highly eccentric and retrograde. However, many uncertainties remain concerning the orbital evolution of these binaries as they either…
▽ More
Binary black holes embedded within the accretion disks that fuel active galactic nuclei (AGN) are promising progenitors for the source of gravitational wave events detected by LIGO/VIRGO. Several recent studies have shown that when these binaries form they should be highly eccentric and retrograde. However, many uncertainties remain concerning the orbital evolution of these binaries as they either inspiral towards merger or disassociate. Previous hydrodynamical simulations exploring their orbital evolution have been predominantly two-dimensional, or have been restricted to binaries on nearly circular orbits. We present the first high-resolution, three-dimensional local shearing-box simulations of both prograde and retrograde eccentric binary black holes embedded in AGN disks. We find that retrograde binaries shrink several times faster than their prograde counterparts and exhibit significant orbital eccentricity growth, the rate of which monotonically increases with binary eccentricity. Our results suggest that retrograde binaries may experience runaway orbital eccentricity growth, which may bring them close enough together at pericenter for gravitational wave emission to drive them to coalescence. Although their eccentricity is damped, prograde binaries shrink much faster than their orbital eccentricity decays, suggesting they should remain modestly eccentric as they contract towards merger. Finally, binary precession driven by the AGN disk may dominate over precession induced by the supermassive black hole depending on the binary accretion rate and its location in the AGN disk, which can subdue the evection resonance and von Ziepel-Lidov-Kozai cycles.
△ Less
Submitted 17 December, 2023; v1 submitted 22 November, 2023;
originally announced November 2023.
-
Spatio-Temporal Anomaly Detection with Graph Networks for Data Quality Monitoring of the Hadron Calorimeter
Authors:
Mulugeta Weldezgina Asres,
Christian Walter Omlin,
Long Wang,
David Yu,
Pavel Parygin,
Jay Dittmann,
Georgia Karapostoli,
Markus Seidel,
Rosamaria Venditti,
Luka Lambrecht,
Emanuele Usai,
Muhammad Ahmad,
Javier Fernandez Menendez,
Kaori Maeshima,
the CMS-HCAL Collaboration
Abstract:
The compact muon solenoid (CMS) experiment is a general-purpose detector for high-energy collision at the large hadron collider (LHC) at CERN. It employs an online data quality monitoring (DQM) system to promptly spot and diagnose particle data acquisition problems to avoid data quality loss. In this study, we present semi-supervised spatio-temporal anomaly detection (AD) monitoring for the physic…
▽ More
The compact muon solenoid (CMS) experiment is a general-purpose detector for high-energy collision at the large hadron collider (LHC) at CERN. It employs an online data quality monitoring (DQM) system to promptly spot and diagnose particle data acquisition problems to avoid data quality loss. In this study, we present semi-supervised spatio-temporal anomaly detection (AD) monitoring for the physics particle reading channels of the hadronic calorimeter (HCAL) of the CMS using three-dimensional digi-occupancy map data of the DQM. We propose the GraphSTAD system, which employs convolutional and graph neural networks to learn local spatial characteristics induced by particles traversing the detector, and global behavior owing to shared backend circuit connections and housing boxes of the channels, respectively. Recurrent neural networks capture the temporal evolution of the extracted spatial features. We have validated the accuracy of the proposed AD system in capturing diverse channel fault types using the LHC Run-2 collision data sets. The GraphSTAD system has achieved production-level accuracy and is being integrated into the CMS core production system--for real-time monitoring of the HCAL. We have also provided a quantitative performance comparison with alternative benchmark models to demonstrate the promising leverage of the presented system.
△ Less
Submitted 7 November, 2023;
originally announced November 2023.
-
The Evolution of Accreting Binaries: from Brown Dwarfs to Supermassive Black Holes
Authors:
Alexander J. Dittmann,
Geoffrey Ryan
Abstract:
Circumbinary accretion occurs throughout the universe, from the formation of stars and planets to the aftermath of major galactic mergers. We present an extensive investigation of circumbinary accretion disks, studying circular binaries with mass ratios ($q\equiv M_2/M_1$) from 0.01 to 1 and at each mass ratio probing the effects of disk thickness and viscosity. We study disks with aspect ratios…
▽ More
Circumbinary accretion occurs throughout the universe, from the formation of stars and planets to the aftermath of major galactic mergers. We present an extensive investigation of circumbinary accretion disks, studying circular binaries with mass ratios ($q\equiv M_2/M_1$) from 0.01 to 1 and at each mass ratio probing the effects of disk thickness and viscosity. We study disks with aspect ratios $H/r\in\{0.1, 0.05, 0.033\}$, and vary both the magnitude and spatial dependance of viscosity. Although thin accretion disks have previously been found to promote rapid inspirals of equal-mass binaries, we find that gravitational torques become weaker at lower mass ratios and most binaries with $0.01\leq q\leq0.04$ outspiral, which may delay the coalescence of black hole binaries formed from minor mergers and cause high-mass exoplanets to migrate outwards. However, in a number of cases, the disks accreting onto binaries with mass ratios $\sim 0.07$ fail to develop eccentric modes, leading to extremely rapid inspirals. Variability in black hole accretion correlates with disk eccentricity, and we observe variability above the $\sim10\%$ level even for mass ratios of $0.01$. We demonstrate that the spatial dependence of the viscosity (e.g. $α$ vs constant-$ν$) significantly affects the degree of preferential accretion onto the secondary, resolving discrepancies between previous studies. Colder circumbinary disks remain eccentric even at $q\sim0.01$ and sustain deep, asymmetric cavities.
△ Less
Submitted 11 October, 2023;
originally announced October 2023.
-
The Evolution of Inclined Binary Black Holes in the Disks of Active Galactic Nuclei
Authors:
Alexander J. Dittmann,
Adam M. Dempsey,
Hui Li
Abstract:
The accretion disks that fuel active galactic nuclei (AGN) may house numerous stars and compact objects, formed in situ or captured from nearby star clusters. Embedded neutron stars and black holes may form binaries and eventually merge, emitting gravitational waves detectable by LIGO/VIRGO. AGN disks are a particularly promising environment for the production of high-mass gravitational wave event…
▽ More
The accretion disks that fuel active galactic nuclei (AGN) may house numerous stars and compact objects, formed in situ or captured from nearby star clusters. Embedded neutron stars and black holes may form binaries and eventually merge, emitting gravitational waves detectable by LIGO/VIRGO. AGN disks are a particularly promising environment for the production of high-mass gravitational wave events involving black holes in the pair instability mass gap, and may facilitate electromagnetic counterparts to black hole binary mergers. However, many orders of magnitude separate the typical length scales of binary formation and those on which gravitational waves can drive binary inspirals, making binary mergers inside the disk uncertain. Previous hydrodynamical simulations of binaries have either been restricted to two dimensions entirely, or focused on binaries aligned with the midplane of the disk. Herein we present the first three-dimensional, high-resolution, local shearing-box hydrodynamical simulations of disk-embedded binaries over a range of orbital inclinations. We find that retrograde binaries can shrink up to four times as quickly as prograde binaries, and that all binaries not perfectly aligned (or anti-aligned) with the AGN disk are driven into alignment. An important consequence of this is that initially retrograde binaries will traverse the inclinations where von Zeipel-Lidov-Kozai oscillations can drive binary eccentricities to large values, potentially facilitating mergers. We also find that interactions with the AGN disk may excite eccentricities in retrograde binaries and cause the orbits of embedded binaries to precess.
△ Less
Submitted 29 October, 2023; v1 submitted 5 October, 2023;
originally announced October 2023.
-
The Decoupling of Binaries from Their Circumbinary Disks
Authors:
Alexander J. Dittmann,
Geoffrey Ryan,
M. Coleman Miller
Abstract:
We have investigated, both analytically and numerically, accreting supermassive black hole binaries as they inspiral due to gravitational radiation to elucidate the decoupling of binaries from their disks and inform future multi-messenger observations of these systems. Our numerical studies evolve equal-mass binaries from initial separations of $100 GM/c^2$ until merger, resolving scales as small…
▽ More
We have investigated, both analytically and numerically, accreting supermassive black hole binaries as they inspiral due to gravitational radiation to elucidate the decoupling of binaries from their disks and inform future multi-messenger observations of these systems. Our numerical studies evolve equal-mass binaries from initial separations of $100 GM/c^2$ until merger, resolving scales as small as $\sim0.04 GM/c^2$, where $M$ is the total binary mass. Our simulations accurately capture the point at which the orbital evolution of each binary decouples from that of their circumbinary disk, and precisely resolve the flow of gas throughout the inspiral. We demonstrate analytically and numerically that timescale-based predictions overestimate the binary separations at which decoupling occurs by factors of $\sim3$, and illustrate the utility of a velocity-based decoupling criterion. High-viscosity ($ν\gtrsim0.03 GM/c$) circumbinary systems decouple late ($a_b\lesssim 15 GM/c^2$) and have qualitatively similar morphologies near merger to circumbinary systems with constant binary separations. Lower-viscosity circumbinary disks decouple earlier and exhibit qualitatively different accretion flows, which lead to precipitously decreasing accretion onto the binary. If detected, such a decrease may unambiguously identify the host galaxy of an ongoing event within a LISA error volume. We illustrate how accretion amplitude and variability evolve as binaries gradually decouple from their circumbinary disks, and where decoupling occurs over the course of binary inspirals in the LISA band. We show that, even when dynamically negligible, gas may leave a detectable imprint on the phase of gravitational waves.
△ Less
Submitted 3 May, 2023; v1 submitted 28 March, 2023;
originally announced March 2023.
-
A Sensitive Search for Supernova Emission Associated with the Extremely Energetic and Nearby GRB 221009A
Authors:
Gokul P. Srinivasaragavan,
Brendan O' Connor,
S. Bradley Cenko,
Alexander J. Dittmann,
Sheng Yang,
Jesper Sollerman,
G. C. Anupama,
Sudhanshu Barway,
Varun Bhalerao,
Harsh Kumar,
Vishwajeet Swain,
Erica Hammerstein,
Isiah Holt,
Shreya Anand,
Igor Andreoni,
Michael W. Coughlin,
Simone Dichiara,
Avishay Gal-Yam,
M. Coleman Miller,
Jaime Soon,
Roberto Soria,
Joseph Durbak,
James H. Gillanders,
Sibasish Laha,
Anna M. Moore
, et al. (2 additional authors not shown)
Abstract:
We report observations of the optical counterpart of the long gamma-ray burst (LGRB) GRB 221009A. Due to the extreme rarity of being both nearby ($z = 0.151$) and highly energetic ($E_{γ,\mathrm{iso}} \geq 10^{54}$ erg), GRB 221009A offers a unique opportunity to probe the connection between massive star core collapse and relativistic jet formation across a very broad range of $γ$-ray properties.…
▽ More
We report observations of the optical counterpart of the long gamma-ray burst (LGRB) GRB 221009A. Due to the extreme rarity of being both nearby ($z = 0.151$) and highly energetic ($E_{γ,\mathrm{iso}} \geq 10^{54}$ erg), GRB 221009A offers a unique opportunity to probe the connection between massive star core collapse and relativistic jet formation across a very broad range of $γ$-ray properties. Adopting a phenomenological power-law model for the afterglow and host galaxy estimates from high-resolution Hubble Space Telescope imaging, we use Bayesian model comparison techniques to determine the likelihood of an associated SN contributing excess flux to the optical light curve. Though not conclusive, we find moderate evidence ($K_{\rm{Bayes}}=10^{1.2}$) for the presence of an additional component arising from an associated supernova, SN 2022xiw, and find that it must be substantially fainter ($<$ 67% as bright at the 99% confidence interval) than SN 1998bw. Given the large and uncertain line-of-sight extinction, we attempt to constrain the supernova parameters ($M_{\mathrm{Ni}}$, $M_{\mathrm{ej}}$, and $E_{\mathrm{KE}}$) under several different assumptions with respect to the host galaxy's extinction. We find properties that are broadly consistent with previous GRB-associated SNe: $M_{\rm{Ni}}=0.05$ - $0.25 \, \rm{M_\odot}$, $M_{\rm{ej}}=3.5$ - $11.1 \, \rm{M_\odot}$, and $E_{\rm{KE}} = (1.6$ - $5.2) \times 10^{52} \, \rm{erg}$. We note that these properties are weakly constrained due to the faintness of the supernova with respect to the afterglow and host emission, but we do find a robust upper limit on the $M_{\rm{Ni}}$ of $M_{\rm{Ni}}<0.36\, \rm{M_\odot}$. Given the tremendous range in isotropic gamma-ray energy release exhibited by GRBs (7 orders of magnitude), the SN emission appears to be decoupled from the central engine in these systems.
△ Less
Submitted 28 May, 2023; v1 submitted 22 March, 2023;
originally announced March 2023.
-
Hubble Space Telescope transmission spectroscopy for the temperate sub-Neptune TOI-270d: a possible hydrogen-rich atmosphere containing water vapour
Authors:
Thomas Mikal-Evans,
Nikku Madhusudhan,
Jason Dittmann,
Maximilian N. Guenther,
Luis Welbanks,
Vincent Van Eylen,
Ian J. M. Crossfield,
Tansu Daylan,
Laura Kreidberg
Abstract:
TOI-270d is a temperate sub-Neptune discovered by the Transiting Exoplanet Survey Satellite (TESS) around a bright (J=9.1mag) M3V host star. With an approximate radius of 2RE and equilibrium temperature of 350K, TOI-270d is one of the most promising small exoplanets for atmospheric characterisation using transit spectroscopy. Here we present a primary transit observation of TOI-270d made with the…
▽ More
TOI-270d is a temperate sub-Neptune discovered by the Transiting Exoplanet Survey Satellite (TESS) around a bright (J=9.1mag) M3V host star. With an approximate radius of 2RE and equilibrium temperature of 350K, TOI-270d is one of the most promising small exoplanets for atmospheric characterisation using transit spectroscopy. Here we present a primary transit observation of TOI-270d made with the Hubble Space Telescope (HST) Wide Field Camera 3 (WFC3) spectrograph across the 1.126-1.644 micron wavelength range, and a 95% credible upper limit of $8.2 \times 10^{-14}$ erg s$^{-1}$ cm$^{-2}$ A$^{-1}$ arcsec$^{-2}$ for the stellar Ly-alpha emission obtained using the Space Telescope Imaging Spectrograph (STIS). The transmission spectrum derived from the TESS and WFC3 data provides evidence for molecular absorption by a hydrogen-rich atmosphere at 4-sigma significance relative to a featureless spectrum. The strongest evidence for any individual absorber is obtained for H2O, which is favoured at 3-sigma significance. When retrieving on the WFC3 data alone and allowing for the possibility of a heterogeneous stellar brightness profile, the detection significance of H2O is reduced to 2.8-sigma. Further observations are therefore required to robustly determine the atmospheric composition of TOI-270d and assess the impact of stellar heterogeneity. If confirmed, our findings would make TOI-270d one of the smallest and coolest exoplanets to date with detected atmospheric spectral features.
△ Less
Submitted 28 November, 2022;
originally announced November 2022.
-
Performance of the CMS High Granularity Calorimeter prototype to charged pion beams of 20$-$300 GeV/c
Authors:
B. Acar,
G. Adamov,
C. Adloff,
S. Afanasiev,
N. Akchurin,
B. Akgün,
M. Alhusseini,
J. Alison,
J. P. Figueiredo de sa Sousa de Almeida,
P. G. Dias de Almeida,
A. Alpana,
M. Alyari,
I. Andreev,
U. Aras,
P. Aspell,
I. O. Atakisi,
O. Bach,
A. Baden,
G. Bakas,
A. Bakshi,
S. Banerjee,
P. DeBarbaro,
P. Bargassa,
D. Barney,
F. Beaudette
, et al. (435 additional authors not shown)
Abstract:
The upgrade of the CMS experiment for the high luminosity operation of the LHC comprises the replacement of the current endcap calorimeter by a high granularity sampling calorimeter (HGCAL). The electromagnetic section of the HGCAL is based on silicon sensors interspersed between lead and copper (or copper tungsten) absorbers. The hadronic section uses layers of stainless steel as an absorbing med…
▽ More
The upgrade of the CMS experiment for the high luminosity operation of the LHC comprises the replacement of the current endcap calorimeter by a high granularity sampling calorimeter (HGCAL). The electromagnetic section of the HGCAL is based on silicon sensors interspersed between lead and copper (or copper tungsten) absorbers. The hadronic section uses layers of stainless steel as an absorbing medium and silicon sensors as an active medium in the regions of high radiation exposure, and scintillator tiles directly readout by silicon photomultipliers in the remaining regions. As part of the development of the detector and its readout electronic components, a section of a silicon-based HGCAL prototype detector along with a section of the CALICE AHCAL prototype was exposed to muons, electrons and charged pions in beam test experiments at the H2 beamline at the CERN SPS in October 2018. The AHCAL uses the same technology as foreseen for the HGCAL but with much finer longitudinal segmentation. The performance of the calorimeters in terms of energy response and resolution, longitudinal and transverse shower profiles is studied using negatively charged pions, and is compared to GEANT4 predictions. This is the first report summarizing results of hadronic showers measured by the HGCAL prototype using beam test data.
△ Less
Submitted 27 May, 2023; v1 submitted 9 November, 2022;
originally announced November 2022.
-
The Radius of PSR J0740+6620 from NICER with NICER Background Estimates
Authors:
Tuomo Salmi,
Serena Vinciguerra,
Devarshi Choudhury,
Thomas E. Riley,
Anna L. Watts,
Ronald A. Remillard,
Paul S. Ray,
Slavko Bogdanov,
Sebastien Guillot,
Zaven Arzoumanian,
Cecilia Chirenti,
Alexander J. Dittmann,
Keith C. Gendreau,
Wynn C. G. Ho,
M. Coleman Miller,
Sharon M. Morsink,
Zorawar Wadiasingh,
Michael T. Wolff
Abstract:
We report a revised analysis for the radius, mass, and hot surface regions of the massive millisecond pulsar PSR J0740+6620, studied previously with joint fits to NICER and XMM-Newton data by Riley et al. (2021) and Miller et al. (2021). We perform a similar Bayesian estimation for the pulse-profile model parameters, except that instead of fitting simultaneously the XMM-Newton data, we use the bes…
▽ More
We report a revised analysis for the radius, mass, and hot surface regions of the massive millisecond pulsar PSR J0740+6620, studied previously with joint fits to NICER and XMM-Newton data by Riley et al. (2021) and Miller et al. (2021). We perform a similar Bayesian estimation for the pulse-profile model parameters, except that instead of fitting simultaneously the XMM-Newton data, we use the best available NICER background estimates to constrain the number of photons detected from the source. This approach eliminates any potential issues in the cross-calibration between these two instruments, providing thus an independent check of the robustness of the analysis. The obtained neutron star parameter constraints are compatible with the already published results, with a slight dependence on how conservative the imposed background limits are. A tighter lower limit causes the inferred radius to increase, and a tighter upper limit causes it to decrease. We also extend the study of the inferred emission geometry to examine the degree of deviation from antipodality of the hot regions. We show that there is a significant offset to an antipodal spot configuration, mainly due to the non-half-cycle azimuthal separation of the two emitting spots. The offset angle from the antipode is inferred to be above 25 degrees with 84% probability. This seems to exclude a centered-dipolar magnetic field in PSR J0740+6620.
△ Less
Submitted 13 January, 2023; v1 submitted 26 September, 2022;
originally announced September 2022.
-
The Influence of Disk Composition on the Evolution of Stars in the Disks of Active Galactic Nuclei
Authors:
Alexander J. Dittmann,
Adam S. Jermyn,
Matteo Cantiello
Abstract:
Disks of gas accreting onto supermassive black holes, powering active galactic nuclei (AGN), can capture stars from nuclear star clusters or form stars in situ via gravitational instability. The density and thermal conditions of these disks can result in rapid accretion onto embedded stars, dramatically altering their evolution in comparison to stars in the interstellar medium. Theoretical models…
▽ More
Disks of gas accreting onto supermassive black holes, powering active galactic nuclei (AGN), can capture stars from nuclear star clusters or form stars in situ via gravitational instability. The density and thermal conditions of these disks can result in rapid accretion onto embedded stars, dramatically altering their evolution in comparison to stars in the interstellar medium. Theoretical models predict that, when subjected to sufficiently rapid accretion, fresh gas replenishes hydrogen in the cores of these stars as quickly as it is burned into helium, reaching a quasi-steady state. Such massive, long-lived ("immortal") stars may be capable of dramatically enriching AGN disks with helium, and would increase the helium abundance in AGN broad-line regions relative to that in the corresponding narrow-line regions and hosts. We investigate how the helium abundance of AGN disks alters the evolution of stars embedded therein. We find, in agreement with analytical arguments, that stars at a given mass are more luminous at higher helium mass fractions, and so undergo more radiation-driven mass-loss. We further find that embedded stars tend to be less massive in disks with higher helium mass fractions, and that immortal stars are less common in such disks. Thus, disk composition can alter the rates of electromagnetic and gravitational wave transients as well as further chemical enrichment by embedded stars.
△ Less
Submitted 12 September, 2022;
originally announced September 2022.
-
Effects of an Immortal Stellar Population in AGN Disks
Authors:
Adam S. Jermyn,
Alexander J. Dittmann,
B. McKernan,
K. E. S. Ford,
Matteo Cantiello
Abstract:
Stars are likely embedded in the gas disks of Active Galactic Nuclei (AGN). Theoretical models predict that in the inner regions of the disk these stars accrete rapidly, with fresh gas replenishing hydrogen in their cores faster than it is burned into helium, effectively stalling their evolution at hydrogen burning. We produce order-of-magnitude estimates of the number of such stars in a fiducial…
▽ More
Stars are likely embedded in the gas disks of Active Galactic Nuclei (AGN). Theoretical models predict that in the inner regions of the disk these stars accrete rapidly, with fresh gas replenishing hydrogen in their cores faster than it is burned into helium, effectively stalling their evolution at hydrogen burning. We produce order-of-magnitude estimates of the number of such stars in a fiducial AGN disk. We find numbers of order $10^{2-4}$, confined to the inner $r_{\rm cap} \sim 3000 r_s \sim 0.03\rm pc$. These stars can profoundly alter the chemistry of AGN disks, enriching them in helium and depleting them in hydrogen, both by order-unity amounts. We further consider mergers between these stars and other disk objects, suggesting that star-star mergers result in rapid mass loss from the remnant to restore an equilibrium mass, while star-compact object mergers may result in exotic outcomes and even host binary black hole mergers within themselves. Finally, we examine how these stars react as the disk dissipates towards the end of its life, and find that they may return mass to the disk fast enough to extend its lifetime by a factor of several and/or may drive powerful outflows from the disk. Post-AGN, these stars rapidly lose mass and form a population of stellar mass black holes around $10M_{\odot}$. Due to the complex and uncertain interactions between embedded stars and the disk, their plausible ubiquity, and their order unity impact on disk structure and evolution, they must be included in realistic disk models.
△ Less
Submitted 11 March, 2022;
originally announced March 2022.
-
A Survey of Disc Thickness and Viscosity in Circumbinary Accretion: Binary Evolution, Variability, and Disc Morphology
Authors:
Alexander J. Dittmann,
Geoffrey Ryan
Abstract:
Much of the parameter space relevant to the evolution of astrophysical circumbinary accretion discs remains unexplored. We have carried out a suite of circumbinary disc simulations surveying both disc thickness and kinematic viscosity, using both constant-$ν$ and constant-$α$ prescriptions. We focus primarily on disc aspect ratios between $0.1$ and $0.033$, and on viscosities between $ν=0.0005$ an…
▽ More
Much of the parameter space relevant to the evolution of astrophysical circumbinary accretion discs remains unexplored. We have carried out a suite of circumbinary disc simulations surveying both disc thickness and kinematic viscosity, using both constant-$ν$ and constant-$α$ prescriptions. We focus primarily on disc aspect ratios between $0.1$ and $0.033$, and on viscosities between $ν=0.0005$ and $ν=0.008$ (in units of binary semi-major axis and orbital frequency), and specialise to circular equal-mass binaries. Both factors strongly influence the evolution of the binary semi-major axis: at $ν=0.0005,$ inspirals occur at aspect ratios $\lesssim0.059$, while at $ν=0.004$ inspirals occur only at aspect ratios $\lesssim0.04$. Inspirals occur largely because of the increasingly strong negative torque on the binary by streams of material which lag the binary, with negligible contributions from resonant torques excited in the circumbinary disc. We find that reductions in accretion rate occur when simulations are initialised too far from the eventual quasi-steady state driven by interaction with the binary, rather than being intrinsically linked to the disc aspect ratio. We find not only that the cavity size increases as viscosity is decreased, but that thinner circumbinary discs become more eccentric. Our results suggest that supermassive black hole binaries should be driven, more rapidly than previous estimates, from $\sim$parsec separations to distances where gravitational waves drive their inspiral, potentially reducing the number of binaries observable by pulsar timing arrays.
△ Less
Submitted 6 June, 2022; v1 submitted 19 January, 2022;
originally announced January 2022.
-
An Analytic, Fully Relativistic Framework for Tidal Disruption Event Streams in Schwarzschild Geometry
Authors:
Alexander J. Dittmann
Abstract:
We present an analytic and fully relativistic framework for studying the self-intersection of tidal disruption event (TDE) streams, restricting ourselves to the Schwarzschild spacetime. By taking advantage of the closed-form solution to the geodesic equations in the Schwarzschild metric, we calculate properties of the self-intersection without numerically evaluating the geodesic equations or makin…
▽ More
We present an analytic and fully relativistic framework for studying the self-intersection of tidal disruption event (TDE) streams, restricting ourselves to the Schwarzschild spacetime. By taking advantage of the closed-form solution to the geodesic equations in the Schwarzschild metric, we calculate properties of the self-intersection without numerically evaluating the geodesic equations or making any post-Newtonian approximations. Our analytic treatment also facilitates geometric definitions of the orbital semi-major axis and eccentricity, as opposed to Newtonian formulas which lead to unphysical results for highly-relativistic orbits. Combined with assumptions about energy dissipation during the self-intersection shock, our framework enables the calculation of quantities such as the fraction of material unbound during the self-intersection shock, and the characteristic semi-major axes and eccentricities of the material which remains in orbit after the collision. As an example, we calculate grids of post-intersection properties in stellar and supermassive black hole (SMBH) masses for disruptions of main sequence stars, identifying regions where no material is ejected during self intersection (e.g. SMBH mass $\lesssim 5\times10^6\, {\rm M_\odot}$ for $1\,{\rm M_\odot}$ stars disrupted at the tidal radius), potentially explaining the TDEs observed by SGR/eROSITA which are visible in X-rays but not optical wavelengths. We also identify parameters for which the post-intersection accretion flow has low eccentricity ($e\lesssim0.6$), and find that the luminosity generated by self-intersection shocks only agrees with observed trends in the relationship between light curve decay timescales and peak luminosities over a narrow range of SMBH masses.
△ Less
Submitted 29 December, 2021;
originally announced December 2021.
-
Response of a CMS HGCAL silicon-pad electromagnetic calorimeter prototype to 20-300 GeV positrons
Authors:
B. Acar,
G. Adamov,
C. Adloff,
S. Afanasiev,
N. Akchurin,
B. Akgün,
F. Alam Khan,
M. Alhusseini,
J. Alison,
A. Alpana,
G. Altopp,
M. Alyari,
S. An,
S. Anagul,
I. Andreev,
P. Aspell,
I. O. Atakisi,
O. Bach,
A. Baden,
G. Bakas,
A. Bakshi,
S. Bannerjee,
P. Bargassa,
D. Barney,
F. Beaudette
, et al. (364 additional authors not shown)
Abstract:
The Compact Muon Solenoid Collaboration is designing a new high-granularity endcap calorimeter, HGCAL, to be installed later this decade. As part of this development work, a prototype system was built, with an electromagnetic section consisting of 14 double-sided structures, providing 28 sampling layers. Each sampling layer has an hexagonal module, where a multipad large-area silicon sensor is glu…
▽ More
The Compact Muon Solenoid Collaboration is designing a new high-granularity endcap calorimeter, HGCAL, to be installed later this decade. As part of this development work, a prototype system was built, with an electromagnetic section consisting of 14 double-sided structures, providing 28 sampling layers. Each sampling layer has an hexagonal module, where a multipad large-area silicon sensor is glued between an electronics circuit board and a metal baseplate. The sensor pads of approximately 1 cm$^2$ are wire-bonded to the circuit board and are readout by custom integrated circuits. The prototype was extensively tested with beams at CERN's Super Proton Synchrotron in 2018. Based on the data collected with beams of positrons, with energies ranging from 20 to 300 GeV, measurements of the energy resolution and linearity, the position and angular resolutions, and the shower shapes are presented and compared to a detailed Geant4 simulation.
△ Less
Submitted 31 March, 2022; v1 submitted 12 November, 2021;
originally announced November 2021.
-
A novel measurement of initial-state gluon radiation in hadron collisions using Drell-Yan events
Authors:
CDF Collaboration,
T. Aaltonen,
S. Amerio,
D. Amidei,
A. Anastassov,
A. Annovi,
J. Antos,
G. Apollinari,
J. A. Appel,
T. Arisawa,
A. Artikov,
J. Asaadi,
W. Ashmanskas,
B. Auerbach,
A. Aurisano,
F. Azfar,
W. Badgett,
T. Bae,
A. Barbaro-Galtieri,
V. E. Barnes,
B. A. Barnett,
P. Barria,
P. Bartos,
M. Bauce,
F. Bedeschi
, et al. (375 additional authors not shown)
Abstract:
A study of initial-state gluon radiation (ISR) in hadron collisions is presented using Drell-Yan (DY) events produced in proton-antiproton collisions by the Tevatron collider at a center-of-mass energy of 1.96 TeV. This paper adopts a novel approach which uses the mean value of the Z/$γ^*$ transverse momentum $<p_T^{DY}>$ in DY events as a powerful observable to characterize the effect of ISR. In…
▽ More
A study of initial-state gluon radiation (ISR) in hadron collisions is presented using Drell-Yan (DY) events produced in proton-antiproton collisions by the Tevatron collider at a center-of-mass energy of 1.96 TeV. This paper adopts a novel approach which uses the mean value of the Z/$γ^*$ transverse momentum $<p_T^{DY}>$ in DY events as a powerful observable to characterize the effect of ISR. In a data sample corresponding to an integrated luminosity of 9.4 fb$^{-1}$ collected with the CDF Run II detector, $<p_T^{DY}>$ is measured as a function of the Z/$γ^*$ invariant mass. It is found that these two observables have a dependence, $<p_T^{DY}> = -8 + 2.2 \ln m_{DY}^2$ [GeV/c], where $m_{DY}$ is the value of the Z/$γ^*$ mass measured in units of GeV/$c^2$. This linear dependence is observed for the first time in this analysis. It may be exploited to model the effect of ISR and constrain its impact in other processes.
△ Less
Submitted 28 October, 2021; v1 submitted 28 October, 2021;
originally announced October 2021.
-
On the Terminal Spins of Accreting Stars and Planets: Boundary Layers
Authors:
Alexander J. Dittmann
Abstract:
The origin of the spins of giant planets is an open question in astrophysics. As planets and stars accrete from discs, if the specific angular momentum accreted corresponds to that of a Keplerian orbit at the surface of the object, it is possible for planets and stars to be spun up to near-breakup speeds. However, accretion cannot proceed onto planets and stars in the same way that accretion proce…
▽ More
The origin of the spins of giant planets is an open question in astrophysics. As planets and stars accrete from discs, if the specific angular momentum accreted corresponds to that of a Keplerian orbit at the surface of the object, it is possible for planets and stars to be spun up to near-breakup speeds. However, accretion cannot proceed onto planets and stars in the same way that accretion proceeds through the disk. For example, the magneto-rotational instability cannot operate in the region between the nearly-Keplerian disk and more slowly-rotating surface because of the sign of the angular velocity gradient. Through this boundary layer where the angular velocity sharply changes, mass and angular momentum transport is thought to be driven by acoustic waves generated by global supersonic shear instabilities and vortices. We present the first study of this mechanism for angular momentum transport around rotating stars and planets using 2D vertically-integrated moving-mesh simulations of ideal hydrodynamics. We find that above rotation rates of $\sim0.4-0.6$ times the Keplerian rate at the surface, depending on the gas sound speed, the rate at which angular momentum is transported inwards through the boundary layer by waves decreases by $\sim1-3$ orders of magnitude. We also find that the accretion rate through the boundary layer decreases commensurately and becomes less variable for faster-rotating objects. Our results provide a purely hydrodynamic mechanism for limiting the spins of accreting planets and stars to factors of a few less than the breakup speed.
△ Less
Submitted 17 August, 2021;
originally announced August 2021.
-
Measurement of the charge asymmetry of electrons from the decays of $W$ bosons produced in $p\bar{p}$ collisions at $\sqrt{s}=1.96$ TeV
Authors:
CDF Collaboration,
T. Aaltonen,
S. Amerio,
D. Amidei,
A. Anastassov,
A. Annovi,
J. Antos,
G. Apollinari,
J. A. Appel,
T. Arisawa,
A. Artikov,
J. Asaadi,
W. Ashmanskas,
B. Auerbach,
A. Aurisano,
F. Azfar,
W. Badgett,
T. Bae,
A. Barbaro-Galtieri,
V. E. Barnes,
B. A. Barnett,
P. Barria,
P. Bartos,
M. Bauce,
F. Bedeschi
, et al. (376 additional authors not shown)
Abstract:
At the Fermilab Tevatron proton-antiproton ($p\bar{p}$) collider, high-mass electron-neutrino ($eν$) pairs are produced predominantly in the process $p \bar{p} \rightarrow W(\rightarrow eν) + X$. The asymmetry of the electron and positron yield as a function of their pseudorapidity constrain the slope of the ratio of the $u$- to $d$-quark parton distributions versus the fraction of the proton mome…
▽ More
At the Fermilab Tevatron proton-antiproton ($p\bar{p}$) collider, high-mass electron-neutrino ($eν$) pairs are produced predominantly in the process $p \bar{p} \rightarrow W(\rightarrow eν) + X$. The asymmetry of the electron and positron yield as a function of their pseudorapidity constrain the slope of the ratio of the $u$- to $d$-quark parton distributions versus the fraction of the proton momentum carried by the quarks. This paper reports on the measurement of the electron-charge asymmetry using the full data set recorded by the Collider Detector at Fermilab in 2001--2011 and corresponding to 9.1~fb$^{-1}$ of integrated luminosity. The measurement significantly improves the precision of the Tevatron constraints on the parton-distribution functions of the proton. Numerical tables of the measurement are provided.
△ Less
Submitted 2 November, 2021; v1 submitted 9 July, 2021;
originally announced July 2021.
-
A Revised Taxonomy of Steganography Embedding Patterns
Authors:
Steffen Wendzel,
Luca Caviglione,
Wojciech Mazurczyk,
Aleksandra Mileva,
Jana Dittmann,
Christian Krätzer,
Kevin Lamshöft,
Claus Vielhauer,
Laura Hartmann,
Jörg Keller,
Tom Neubert
Abstract:
Steganography embraces several hiding techniques which spawn across multiple domains. However, the related terminology is not unified among the different domains, such as digital media steganography, text steganography, cyber-physical systems steganography, network steganography (network covert channels), local covert channels, and out-of-band covert channels. To cope with this, a prime attempt ha…
▽ More
Steganography embraces several hiding techniques which spawn across multiple domains. However, the related terminology is not unified among the different domains, such as digital media steganography, text steganography, cyber-physical systems steganography, network steganography (network covert channels), local covert channels, and out-of-band covert channels. To cope with this, a prime attempt has been done in 2015, with the introduction of the so-called hiding patterns, which allow to describe hiding techniques in a more abstract manner. Despite significant enhancements, the main limitation of such a taxonomy is that it only considers the case of network steganography.
Therefore, this paper reviews both the terminology and the taxonomy of hiding patterns as to make them more general. Specifically, hiding patterns are split into those that describe the embedding and the representation of hidden data within the cover object.
As a first research action, we focus on embedding hiding patterns and we show how they can be applied to multiple domains of steganography instead of being limited to the network scenario. Additionally, we exemplify representation patterns using network steganography. Our pattern collection is available under https://patterns.ztt.hs-worms.de.
△ Less
Submitted 16 June, 2021;
originally announced June 2021.
-
The Featureless HST/WFC3 Transmission Spectrum of the Rocky Exoplanet GJ 1132b: No Evidence For A Cloud-Free Primordial Atmosphere and Constraints on Starspot Contamination
Authors:
Jessica E. Libby-Roberts,
Zachory K. Berta-Thompson,
Hannah Diamond-Lowe,
Michael A. Gully-Santiago,
Jonathan M. Irwin,
Eliza M. -R. Kempton,
Benjamin V. Rackham,
David Charbonneau,
Jean-Michel Desert,
Jason A. Dittmann,
Ryan Hofmann,
Caroline V. Morley,
Elisabeth R. Newton
Abstract:
Orbiting a M dwarf 12 pc away, the transiting exoplanet GJ 1132b is a prime target for transmission spectroscopy. With a mass of 1.7 Earth masses and radius of 1.1 Earth radii, GJ 1132b's bulk density indicates that this planet is rocky. Yet with an equilibrium temperature of 580 K, GJ 1132b may still retain some semblance of an atmosphere. Understanding whether this atmosphere exists and its comp…
▽ More
Orbiting a M dwarf 12 pc away, the transiting exoplanet GJ 1132b is a prime target for transmission spectroscopy. With a mass of 1.7 Earth masses and radius of 1.1 Earth radii, GJ 1132b's bulk density indicates that this planet is rocky. Yet with an equilibrium temperature of 580 K, GJ 1132b may still retain some semblance of an atmosphere. Understanding whether this atmosphere exists and its composition will be vital for understanding how the atmospheres of terrestrial planets orbiting M dwarfs evolve. We observe five transits of GJ 1132b with the Wide Field Camera 3 (WFC3) on the Hubble Space Telescope (HST). We find a featureless transmission spectrum from 1.1--1.7 microns, ruling out cloud-free atmospheres with metallicities <300x Solar with >4.8$σ$ confidence. We combine our WFC3 results with transit depths from TESS and archival broadband and spectroscopic observations to find a featureless spectrum from 0.7--4.5 microns. GJ 1132b has a high mean molecular weight atmosphere, possesses a high-altitude aerosol layer, or has effectively no atmosphere. Higher precision observations are required to differentiate between these possibilities. We explore the impact of hot and cold starspots on the observed transmission spectrum GJ 1132b, quantifying the amplitude of spot-induced transit depth features. Using a simple Poisson model we estimate spot temperature contrasts, spot covering fractions, and spot sizes for GJ 1132. These limits, and the modeling framework, may be useful for future observations of GJ 1132b or other planets transiting similarly inactive M dwarfs.
△ Less
Submitted 21 May, 2021;
originally announced May 2021.
-
The Radius of PSR J0740+6620 from NICER and XMM-Newton Data
Authors:
M. C. Miller,
F. K. Lamb,
A. J. Dittmann,
S. Bogdanov,
Z. Arzoumanian,
K. C. Gendreau,
S. Guillot,
W. C. G. Ho,
J. M. Lattimer,
M. Loewenstein,
S. M. Morsink,
P. S. Ray,
M. T. Wolff,
C. L. Baker,
T. Cazeau,
S. Manthripragada,
C. B. Markwardt,
T. Okajima,
S. Pollard,
I. Cognard,
H. T. Cromartie,
E. Fonseca,
L. Guillemot,
M. Kerr,
A. Parthasarathy
, et al. (3 additional authors not shown)
Abstract:
PSR J0740$+$6620 has a gravitational mass of $2.08\pm 0.07~M_\odot$, which is the highest reliably determined mass of any neutron star. As a result, a measurement of its radius will provide unique insight into the properties of neutron star core matter at high densities. Here we report a radius measurement based on fits of rotating hot spot patterns to Neutron Star Interior Composition Explorer (N…
▽ More
PSR J0740$+$6620 has a gravitational mass of $2.08\pm 0.07~M_\odot$, which is the highest reliably determined mass of any neutron star. As a result, a measurement of its radius will provide unique insight into the properties of neutron star core matter at high densities. Here we report a radius measurement based on fits of rotating hot spot patterns to Neutron Star Interior Composition Explorer (NICER) and X-ray Multi-Mirror (XMM-Newton) X-ray observations. We find that the equatorial circumferential radius of PSR J0740$+$6620 is $13.7^{+2.6}_{-1.5}$ km (68%). We apply our measurement, combined with the previous NICER mass and radius measurement of PSR J0030$+$0451, the masses of two other $\sim 2~M_\odot$ pulsars, and the tidal deformability constraints from two gravitational wave events, to three different frameworks for equation of state modeling, and find consistent results at $\sim 1.5-3$ times nuclear saturation density. For a given framework, when all measurements are included the radius of a $1.4~M_\odot$ neutron star is known to $\pm 4$% (68% credibility) and the radius of a $2.08~M_\odot$ neutron star is known to $\pm 5$%. The full radius range that spans the $\pm 1σ$ credible intervals of all the radius estimates in the three frameworks is $12.45\pm 0.65$ km for a $1.4~M_\odot$ neutron star and $12.35\pm 0.75$ km for a $2.08~M_\odot$ neutron star.
△ Less
Submitted 14 May, 2021;
originally announced May 2021.
-
NICER Detection of Thermal X-ray Pulsations from the Massive Millisecond Pulsars PSR J0740+6620 and PSR J1614-2230
Authors:
M. T. Wolff,
S. Guillot,
S. Bogdanov,
P. S. Ray,
M. Kerr,
Z. Arzoumanian,
K. C. Gendreau,
M. C. Miller,
A. J. Dittmann,
W. C. G. Ho,
L. Guillemot,
I. Cognard,
G. Theureau,
K. S. Wood
Abstract:
We report the detection of X-ray pulsations from the rotation-powered millisecond-period pulsars PSR J0740+6620 and PSR J1614-2230, two of the most massive neutron stars known, using observations with the Neutron Star Interior Composition Explorer (NICER). We also analyze X-ray Multi-Mirror Mission (XMM-Newton) data for both pulsars to obtain their time-averaged fluxes and study their respective X…
▽ More
We report the detection of X-ray pulsations from the rotation-powered millisecond-period pulsars PSR J0740+6620 and PSR J1614-2230, two of the most massive neutron stars known, using observations with the Neutron Star Interior Composition Explorer (NICER). We also analyze X-ray Multi-Mirror Mission (XMM-Newton) data for both pulsars to obtain their time-averaged fluxes and study their respective X-ray fields. PSR J0740+6620 exhibits a broad double-peaked profile with a separation of ~0.4 in phase. PSR J1614-2230, on the other hand, has a broad single-peak profile. The broad modulations with soft X-ray spectra of both pulsars are indicative of thermal radiation from one or more small regions of the stellar surface. We show the NICER detections of X-ray pulsations for both pulsars and also discuss the phase relationship to their radio pulsations. In the case of PSR J0740+6620, this paper documents the data reduction performed to obtain the pulsation detection and prepare for pulse profile modeling analysis.
△ Less
Submitted 14 May, 2021;
originally announced May 2021.
-
Constraining the Neutron Star Mass--Radius Relation and Dense Matter Equation of State with NICER. III. Model Description and Verification of Parameter Estimation Codes
Authors:
Slavko Bogdanov,
Alexander J. Dittmann,
Wynn C. G. Ho,
Frederick K. Lamb,
Simin Mahmoodifar,
M. Coleman Miller,
Sharon M. Morsink,
Thomas E. Riley,
Tod E. Strohmayer,
Anna L. Watts,
Devarshi Choudhury,
Sebastien Guillot,
Alice K. Harding,
Paul S. Ray,
Zorawar Wadiasingh,
Michael T. Wolff,
Craig B. Markwardt,
Zaven Arzoumanian,
Keith C. Gendreau
Abstract:
We describe the X-ray pulse profile models we use, and how we use them, to analyze Neutron Star Interior Composition Explorer (NICER) observations of rotation-powered millisecond pulsars to obtain information about the mass-radius relation of neutron stars and the equation of state of the dense matter in their cores. Here we detail our modeling of the observed profile of PSR J0030+0451 that we ana…
▽ More
We describe the X-ray pulse profile models we use, and how we use them, to analyze Neutron Star Interior Composition Explorer (NICER) observations of rotation-powered millisecond pulsars to obtain information about the mass-radius relation of neutron stars and the equation of state of the dense matter in their cores. Here we detail our modeling of the observed profile of PSR J0030+0451 that we analyzed in Miller et al. (2019) and Riley et al. (2019) and describe a cross-verification of computations of the pulse profiles of a star with R/M 3, in case stars this compact need to be considered in future analyses. We also present our early cross-verification efforts of the parameter estimation procedures used by Miller et al. (2019) and Riley et al. (2019) by analyzing two distinct synthetic data sets. Both codes yielded credible regions in the mass-radius plane that are statistically consistent with one another and both gave posterior distributions for model parameter values consistent with the values that were used to generate the data. We also summarize the additional tests of the parameter estimation procedure of Miller et al. (2019) that used synthetic pulse profiles and the NICER pulse profile of PSR J0030+0451. We then illustrate how the precision of mass and radius estimates depends on the pulsar's spin rate and the size of its hot spot by analyzing four different synthetic pulse profiles. Finally, we assess possible sources of systematic error in these estimates made using this technique, some of which may warrant further investigation.
△ Less
Submitted 14 April, 2021;
originally announced April 2021.
-
The TESS Objects of Interest Catalog from the TESS Prime Mission
Authors:
Natalia M. Guerrero,
S. Seager,
Chelsea X. Huang,
Andrew Vanderburg,
Aylin Garcia Soto,
Ismael Mireles,
Katharine Hesse,
William Fong,
Ana Glidden,
Avi Shporer,
David W. Latham,
Karen A. Collins,
Samuel N. Quinn,
Jennifer Burt,
Diana Dragomir,
Ian Crossfield,
Roland Vanderspek,
Michael Fausnaugh,
Christopher J. Burke,
George Ricker,
Tansu Daylan,
Zahra Essack,
Maximilian N. Günther,
Hugh P. Osborn,
Joshua Pepper
, et al. (80 additional authors not shown)
Abstract:
We present 2,241 exoplanet candidates identified with data from the Transiting Exoplanet Survey Satellite (TESS) during its two-year prime mission. We list these candidates in the TESS Objects of Interest (TOI) Catalog, which includes both new planet candidates found by TESS and previously-known planets recovered by TESS observations. We describe the process used to identify TOIs and investigate t…
▽ More
We present 2,241 exoplanet candidates identified with data from the Transiting Exoplanet Survey Satellite (TESS) during its two-year prime mission. We list these candidates in the TESS Objects of Interest (TOI) Catalog, which includes both new planet candidates found by TESS and previously-known planets recovered by TESS observations. We describe the process used to identify TOIs and investigate the characteristics of the new planet candidates, and discuss some notable TESS planet discoveries. The TOI Catalog includes an unprecedented number of small planet candidates around nearby bright stars, which are well-suited for detailed follow-up observations. The TESS data products for the Prime Mission (Sectors 1-26), including the TOI Catalog, light curves, full-frame images, and target pixel files, are publicly available on the Mikulski Archive for Space Telescopes.
△ Less
Submitted 24 March, 2021; v1 submitted 23 March, 2021;
originally announced March 2021.
-
Stellar Evolution in the Disks of Active Galactic Nuclei Produces Rapidly Rotating Massive Stars
Authors:
Adam S. Jermyn,
Alexander J. Dittmann,
Matteo Cantiello,
Rosalba Perna
Abstract:
Stars can either be formed in or captured by the accretion disks in Active Galactic Nuclei (AGN). These AGN stars are irradiated and subject to extreme levels of accretion, which can turn even low-mass stars into very massive ones ($M > 100 {\rm M}_\odot$) whose evolution may result in the formation of massive compact objects ($M > 10 {\rm M}_\odot$). Here we explore the spins of these AGN stars a…
▽ More
Stars can either be formed in or captured by the accretion disks in Active Galactic Nuclei (AGN). These AGN stars are irradiated and subject to extreme levels of accretion, which can turn even low-mass stars into very massive ones ($M > 100 {\rm M}_\odot$) whose evolution may result in the formation of massive compact objects ($M > 10 {\rm M}_\odot$). Here we explore the spins of these AGN stars and the remnants they leave behind. We find that AGN stars rapidly spin up via accretion, eventually reaching near-critical rotation rates. They further maintain near-critical rotation even as they shed their envelopes, become compact, and undergo late stages of burning. This makes them good candidates to produce high-spin massive black holes, such as the ones seen by LIGO-Virgo in GW190521g, as well as long Gamma Ray Bursts (GRBs) and the associated chemical pollution of the AGN disk.
△ Less
Submitted 25 February, 2021;
originally announced February 2021.
-
Accretion onto Stars in the Disks of Active Galactic Nuclei
Authors:
Alexander J. Dittmann,
Matteo Cantiello,
Adam S. Jermyn
Abstract:
Disks of gas accreting onto supermassive black holes are thought to power active galactic nuclei (AGN). Stars may form in gravitationally unstable regions of these disks, or may be captured from nuclear star clusters. Because of the dense gas environment, the evolution of such embedded stars can diverge dramatically from those in the interstellar medium. This work extends previous studies of stell…
▽ More
Disks of gas accreting onto supermassive black holes are thought to power active galactic nuclei (AGN). Stars may form in gravitationally unstable regions of these disks, or may be captured from nuclear star clusters. Because of the dense gas environment, the evolution of such embedded stars can diverge dramatically from those in the interstellar medium. This work extends previous studies of stellar evolution in AGN disks by exploring a variety of ways that accretion onto stars in AGN disks may differ from Bondi accretion. We find that tidal effects from the supermassive black hole significantly alter the evolution of stars in AGN disks, and that our results do not depend critically on assumptions about radiative feedback on the accretion stream. Thus, in addition to depending on $ρ/c_s^3$, the fate of stars in AGN disks depends sensitively on the distance to and mass of the supermassive black hole. This affects where in the disk stellar explosions occur, where compact remnants form and potentially merge to produce gravitational waves, and where different types of chemical enrichment take place.
△ Less
Submitted 2 August, 2021; v1 submitted 24 February, 2021;
originally announced February 2021.
-
Higher Compact Multiple Occurrence Around Metal-Poor M-Dwarfs and Late K-Dwarfs
Authors:
Sophie G. Anderson,
Jason A. Dittmann,
Sarah Ballard,
Megan Bedell
Abstract:
The planet-metallicity correlation serves as a potential link between exoplanet systems as we observe them today and the effects of bulk composition on the planet formation process. Many observers have noted a tendency for Jovian planets to form around stars with higher metallicities; however, there is no consensus on a trend for smaller planets. Here, we investigate the planet-metallicity correla…
▽ More
The planet-metallicity correlation serves as a potential link between exoplanet systems as we observe them today and the effects of bulk composition on the planet formation process. Many observers have noted a tendency for Jovian planets to form around stars with higher metallicities; however, there is no consensus on a trend for smaller planets. Here, we investigate the planet-metallicity correlation for rocky planets in single and multi-planet systems around Kepler M-dwarf and late K-dwarf stars. Due to molecular blanketing and the dim nature of these low mass stars, it is difficult to make direct elemental abundance measurements via spectroscopy. We instead use a combination of accurate and uniformly measured parallaxes and photometry to obtain relative metallicities and validate this method with a subsample of spectroscopically determined metallicities. We use the Kolmogorov-Smirnov (KS) test, Mann-Whitney U test, and Anderson-Darling test to compare the compact multiple planetary systems with single transiting planet systems and systems with no detected transiting planets. We find that the compact multiple planetary systems are derived from a statistically more metal-poor population, with a p-value of 0.015 in the KS test, a p-value of 0.005 in the Mann-Whitney U test, and a value of 2.574 in the Anderson-Darling test statistic, which exceeds the derived threshold for significance by a factor of 25. We conclude that metallicity plays a significant role in determining the architecture of rocky planet systems. Compact multiples either form more readily, or are more likely to survive on Gyr timescales, around metal-poor stars.
△ Less
Submitted 17 February, 2021;
originally announced February 2021.
-
Construction and commissioning of CMS CE prototype silicon modules
Authors:
B. Acar,
G. Adamov,
C. Adloff,
S. Afanasiev,
N. Akchurin,
B. Akgün,
M. Alhusseini,
J. Alison,
G. Altopp,
M. Alyari,
S. An,
S. Anagul,
I. Andreev,
M. Andrews,
P. Aspell,
I. A. Atakisi,
O. Bach,
A. Baden,
G. Bakas,
A. Bakshi,
P. Bargassa,
D. Barney,
E. Becheva,
P. Behera,
A. Belloni
, et al. (307 additional authors not shown)
Abstract:
As part of its HL-LHC upgrade program, the CMS Collaboration is developing a High Granularity Calorimeter (CE) to replace the existing endcap calorimeters. The CE is a sampling calorimeter with unprecedented transverse and longitudinal readout for both electromagnetic (CE-E) and hadronic (CE-H) compartments. The calorimeter will be built with $\sim$30,000 hexagonal silicon modules. Prototype modul…
▽ More
As part of its HL-LHC upgrade program, the CMS Collaboration is developing a High Granularity Calorimeter (CE) to replace the existing endcap calorimeters. The CE is a sampling calorimeter with unprecedented transverse and longitudinal readout for both electromagnetic (CE-E) and hadronic (CE-H) compartments. The calorimeter will be built with $\sim$30,000 hexagonal silicon modules. Prototype modules have been constructed with 6-inch hexagonal silicon sensors with cell areas of 1.1~$cm^2$, and the SKIROC2-CMS readout ASIC. Beam tests of different sampling configurations were conducted with the prototype modules at DESY and CERN in 2017 and 2018. This paper describes the construction and commissioning of the CE calorimeter prototype, the silicon modules used in the construction, their basic performance, and the methods used for their calibration.
△ Less
Submitted 10 December, 2020;
originally announced December 2020.
-
The DAQ system of the 12,000 Channel CMS High Granularity Calorimeter Prototype
Authors:
B. Acar,
G. Adamov,
C. Adloff,
S. Afanasiev,
N. Akchurin,
B. Akgün,
M. Alhusseini,
J. Alison,
G. Altopp,
M. Alyari,
S. An,
S. Anagul,
I. Andreev,
M. Andrews,
P. Aspell,
I. A. Atakisi,
O. Bach,
A. Baden,
G. Bakas,
A. Bakshi,
P. Bargassa,
D. Barney,
E. Becheva,
P. Behera,
A. Belloni
, et al. (307 additional authors not shown)
Abstract:
The CMS experiment at the CERN LHC will be upgraded to accommodate the 5-fold increase in the instantaneous luminosity expected at the High-Luminosity LHC (HL-LHC). Concomitant with this increase will be an increase in the number of interactions in each bunch crossing and a significant increase in the total ionising dose and fluence. One part of this upgrade is the replacement of the current endca…
▽ More
The CMS experiment at the CERN LHC will be upgraded to accommodate the 5-fold increase in the instantaneous luminosity expected at the High-Luminosity LHC (HL-LHC). Concomitant with this increase will be an increase in the number of interactions in each bunch crossing and a significant increase in the total ionising dose and fluence. One part of this upgrade is the replacement of the current endcap calorimeters with a high granularity sampling calorimeter equipped with silicon sensors, designed to manage the high collision rates. As part of the development of this calorimeter, a series of beam tests have been conducted with different sampling configurations using prototype segmented silicon detectors. In the most recent of these tests, conducted in late 2018 at the CERN SPS, the performance of a prototype calorimeter equipped with ${\approx}12,000\rm{~channels}$ of silicon sensors was studied with beams of high-energy electrons, pions and muons. This paper describes the custom-built scalable data acquisition system that was built with readily available FPGA mezzanines and low-cost Raspberry PI computers.
△ Less
Submitted 8 December, 2020; v1 submitted 7 December, 2020;
originally announced December 2020.
-
High-Order Multiderivative IMEX Schemes
Authors:
Alexander J. Dittmann
Abstract:
Recently, a 4th-order asymptotic preserving multiderivative implicit-explicit (IMEX) scheme was developed (Schütz and Seal 2020, arXiv:2001.08268). This scheme is based on a 4th-order Hermite interpolation in time, and uses an approach based on operator splitting that converges to the underlying quadrature if iterated sufficiently. Hermite schemes have been used in astrophysics for decades, partic…
▽ More
Recently, a 4th-order asymptotic preserving multiderivative implicit-explicit (IMEX) scheme was developed (Schütz and Seal 2020, arXiv:2001.08268). This scheme is based on a 4th-order Hermite interpolation in time, and uses an approach based on operator splitting that converges to the underlying quadrature if iterated sufficiently. Hermite schemes have been used in astrophysics for decades, particularly for N-body calculations, but not in a form suitable for solving stiff equations. In this work, we extend the scheme presented in Schütz and Seal 2020 to higher orders. Such high-order schemes offer advantages when one aims to find high-precision solutions to systems of differential equations containing stiff terms, which occur throughout the physical sciences. We begin by deriving Hermite schemes of arbitrary order and discussing the stability of these formulas. Afterwards, we demonstrate how the method of Schütz and Seal 2020 generalises in a straightforward manner to any of these schemes, and prove convergence properties of the resulting IMEX schemes. We then present results for methods ranging from 6th to 12th order and explore a selection of test problems, including both linear and nonlinear ordinary differential equations and Burgers' equation. To our knowledge this is also the first time that Hermite time-stepping methods have been applied to partial differential equations. We then discuss some benefits of these schemes, such as their potential for parallelism and low memory usage, as well as limitations and potential drawbacks.
△ Less
Submitted 11 August, 2020; v1 submitted 7 August, 2020;
originally announced August 2020.
-
Modified Hermite Integrators of Arbitrary Order
Authors:
Alexander J. Dittmann
Abstract:
We present a family of modified Hermite integrators of arbitrary order possessing superior behaviour for the integration of Keplerian and near-Keplerian orbits. After recounting the derivation of Hermite N-body integrators of arbitrary order, we derive a corrector expression that minimises integrated errors in the argument of periapsis for any such integrator. In addition to providing an alternate…
▽ More
We present a family of modified Hermite integrators of arbitrary order possessing superior behaviour for the integration of Keplerian and near-Keplerian orbits. After recounting the derivation of Hermite N-body integrators of arbitrary order, we derive a corrector expression that minimises integrated errors in the argument of periapsis for any such integrator. In addition to providing an alternate derivation of the modified corrector for the 4th-order Hermite integrator, we focus on improved correctors for the 6th- and 8th-order integrators previously featured in the literature. We present a set of numerical examples and find that the higher-order schemes improve performance, even when considering their slightly higher cost in floating point operations. The algorithms presented herein hold promise for systems dominated by central potentials, such as planetary systems and the centres of galaxies. Existing Hermite integrators of any order can be modified to use the expressions presented here with minimal effort. Accordingly the schemes presented herein can be easily implemented on massively parallel architectures.
△ Less
Submitted 21 June, 2020; v1 submitted 3 June, 2020;
originally announced June 2020.
-
Constraining the Neutron Star Mass-Radius Relation and Dense Matter Equation of State with NICER. II. Emission from Hot Spots on a Rapidly Rotating Neutron Star
Authors:
Slavko Bogdanov,
Frederick K. Lamb,
Simin Mahmoodifar,
M. Coleman Miller,
Sharon M. Morsink,
Thomas E. Riley,
Tod E. Strohmayer,
Anna L. Watts,
Alexander J. Dittmann,
Deepto Chakrabarty,
Sebastien Guillot,
Zaven Arzoumanian,
Keith C. Gendreau
Abstract:
We describe the model of surface emission from a rapidly rotating neutron star that is applied to Neutron Star Interior Composition Explorer X-ray data of millisecond pulsars in order to statistically constrain the neutron star mass-radius relation and dense matter equation of state. To ensure that the associated calculations are both accurate and precise, we conduct an extensive suite of verifica…
▽ More
We describe the model of surface emission from a rapidly rotating neutron star that is applied to Neutron Star Interior Composition Explorer X-ray data of millisecond pulsars in order to statistically constrain the neutron star mass-radius relation and dense matter equation of state. To ensure that the associated calculations are both accurate and precise, we conduct an extensive suite of verification tests between our numerical codes for both the Schwarzschild + Doppler and Oblate Schwarzschild approximations, and compare both approximations against exact numerical calculations. We find superb agreement between the code outputs, as well as in comparison against a set of analytical and semi-analytical calculations, which combined with their speed, demonstrates that the codes are well-suited for large-scale statistical sampling applications. A set of verified, high-precision reference synthetic pulse profiles is provided to the community to facilitate testing of other independently developed codes.
△ Less
Submitted 11 December, 2019;
originally announced December 2019.
-
PSR J0030+0451 Mass and Radius from NICER Data and Implications for the Properties of Neutron Star Matter
Authors:
M. C. Miller,
F. K. Lamb,
A. J. Dittmann,
S. Bogdanov,
Z. Arzoumanian,
K. C. Gendreau,
S. Guillot,
A. K. Harding,
W. C. G. Ho,
J. M. Lattimer,
R. M. Ludlam,
S. Mahmoodifar,
S. M. Morsink,
P. S. Ray,
T. E. Strohmayer,
K. S. Wood,
T. Enoto,
R. Foster,
T. Okajima,
G. Prigozhin,
Y. Soong
Abstract:
Neutron stars are not only of astrophysical interest, but are also of great interest to nuclear physicists, because their attributes can be used to determine the properties of the dense matter in their cores. One of the most informative approaches for determining the equation of state of this dense matter is to measure both a star's equatorial circumferential radius $R_e$ and its gravitational mas…
▽ More
Neutron stars are not only of astrophysical interest, but are also of great interest to nuclear physicists, because their attributes can be used to determine the properties of the dense matter in their cores. One of the most informative approaches for determining the equation of state of this dense matter is to measure both a star's equatorial circumferential radius $R_e$ and its gravitational mass $M$. Here we report estimates of the mass and radius of the isolated 205.53 Hz millisecond pulsar PSR J0030+0451 obtained using a Bayesian inference approach to analyze its energy-dependent thermal X-ray waveform, which was observed using the Neutron Star Interior Composition Explorer (NICER). This approach is thought to be less subject to systematic errors than other approaches for estimating neutron star radii. We explored a variety of emission patterns on the stellar surface. Our best-fit model has three oval, uniform-temperature emitting spots and provides an excellent description of the pulse waveform observed using NICER. The radius and mass estimates given by this model are $R_e = 13.02^{+1.24}_{-1.06}$ km and $M = 1.44^{+0.15}_{-0.14}\ M_\odot$ (68%). The independent analysis reported in the companion paper by Riley et al. (2019) explores different emitting spot models, but finds spot shapes and locations and estimates of $R_e$ and $M$ that are consistent with those found in this work. We show that our measurements of $R_e$ and $M$ for PSR J0030$+$0451 improve the astrophysical constraints on the equation of state of cold, catalyzed matter above nuclear saturation density.
△ Less
Submitted 11 December, 2019;
originally announced December 2019.
-
The young planetary system K2-25: constraints on companions and starspots
Authors:
Isabel J. Kain,
Elisabeth R. Newton,
Jason A. Dittmann,
Jonathan M. Irwin,
Andrew W. Mann,
Pa Chia Thao,
David Charbonneau,
Jennifer G. Winters
Abstract:
The abundance of planets with orbital periods of a few to tens of days suggests that exoplanets experience complex dynamical histories. Planets in young stellar clusters or associations have well-constrained ages and therefore provide an opportunity to explore the dynamical evolution of exoplanets. K2-25b is a Neptune-sized planet in an eccentric, 3.48 day orbit around an M4.5 dwarf star in the Hy…
▽ More
The abundance of planets with orbital periods of a few to tens of days suggests that exoplanets experience complex dynamical histories. Planets in young stellar clusters or associations have well-constrained ages and therefore provide an opportunity to explore the dynamical evolution of exoplanets. K2-25b is a Neptune-sized planet in an eccentric, 3.48 day orbit around an M4.5 dwarf star in the Hyades cluster (650 Myr). In order to investigate its non-zero eccentricity and tight orbit, we analyze transit timing variations (TTVs) which could reveal clues to the migration processes that may have acted on the planet. We obtain 12 non-consecutive transits using the MEarth Observatories and long-term photometric monitoring, which we combine with 10 transits from the Spitzer Space Telescope and 20 transits from K2. Tables of MEarth photometry accompany this work. We fit each transit lightcurve independently. We first investigate whether inhomogeneities on the stellar surface (such as spots or plages) are differentially affecting our transit observations. The measured transit depth does not vary significantly between transits, though we see some deviations from the fiducial transit model. We then looked for TTVs as evidence of a non-transiting perturber in the system. We find no evidence for > 1 $M_\oplus$ companions within a 2:1 period ratio, or for > 5 $M_\oplus$ planets within a 7:2 period ratio.
△ Less
Submitted 11 December, 2019;
originally announced December 2019.
-
Star Formation in Accretion Disks and SMBH Growth
Authors:
Alexander J. Dittmann,
M. Coleman Miller
Abstract:
Accretion disks around active galactic nuclei are potentially unstable to star formation at large radii. We note that when the compact objects formed from some of these stars spiral into the central supermassive black hole, there is no radiative feedback and therefore the accretion rate is not limited by radiation forces. Using a set of accretion disk models, we calculate the accretion rate onto t…
▽ More
Accretion disks around active galactic nuclei are potentially unstable to star formation at large radii. We note that when the compact objects formed from some of these stars spiral into the central supermassive black hole, there is no radiative feedback and therefore the accretion rate is not limited by radiation forces. Using a set of accretion disk models, we calculate the accretion rate onto the central supermassive black hole in both gas and compact objects. We find that the timescale for a supermassive black hole to double in mass can decrease by factors ranging from $\sim0.7$ to as low as $\sim0.1$ in extreme cases, compared to gas accretion alone. Our results suggest that the formation of extremely massive black holes at high redshift may occur without prolonged super-Eddington gas accretion or very massive seed black holes. We comment on potential observational signatures as well as implications for other observations of active galactic nuclei.
△ Less
Submitted 19 November, 2019;
originally announced November 2019.
-
Simultaneous Optical Transmission Spectroscopy of a Terrestrial, Habitable-Zone Exoplanet with Two Ground-Based Multi-Object Spectrographs
Authors:
Hannah Diamond-Lowe,
Zachory Berta-Thompson,
David Charbonneau,
Jason Dittmann,
Eliza M. -R. Kempton
Abstract:
Investigating the atmospheres of rocky exoplanets is key to performing comparative planetology between these worlds and the terrestrial planets that reside in the inner solar system. Terrestrial exoplanet atmospheres exhibit weak signals, and attempting to detect them pushes at the boundaries of what is possible for current instrumentation. We focus on the habitable-zone terrestrial exoplanet LHS…
▽ More
Investigating the atmospheres of rocky exoplanets is key to performing comparative planetology between these worlds and the terrestrial planets that reside in the inner solar system. Terrestrial exoplanet atmospheres exhibit weak signals, and attempting to detect them pushes at the boundaries of what is possible for current instrumentation. We focus on the habitable-zone terrestrial exoplanet LHS 1140b. Given its 25-day orbital period and 2 hr transit duration, capturing transits of LHS 1140b is challenging. We observed two transits of this object, approximately 1 yr apart, which yielded four data sets thanks to our simultaneous use of the IMACS and LDSS3C multiobject spectrographs mounted on the twin Magellan telescopes at Las Campanas Observatory. We present a jointly fit white light curve, as well as jointly fit 20 nm wavelength-binned light curves from which we construct a transmission spectrum. Binning the joint white light-curve residuals to 3-minute time bins gives an rms of 145 ppm; binning down to 10-minute time bins gives an rms of 77 ppm. Our median uncertainty in Rp^2/Rs^2 in the 20 nm wavelength bins is 260 ppm, and we achieve an average precision of 1.3x the photon noise when fitting the wavelength-binned light curves with a Gaussian process regression. Our precision on Rp^2/Rs^2 is a factor of four larger than the feature amplitudes of a clear, hydrogen-dominated atmosphere, meaning that we are not able to test realistic models of LHS 1140b's atmosphere. The techniques and caveats presented here are applicable to the growing sample of terrestrial worlds in the Transiting Exoplanet Survey Satellite era, as well as to the upcoming generation of ground-based giant segmented mirror telescopes.
△ Less
Submitted 2 June, 2020; v1 submitted 19 September, 2019;
originally announced September 2019.
-
Absence of a thick atmosphere on the terrestrial exoplanet LHS 3844b
Authors:
Laura Kreidberg,
Daniel D. B. Koll,
Caroline Morley,
Renyu Hu,
Laura Schaefer,
Drake Deming,
Kevin B. Stevenson,
Jason Dittmann,
Andrew Vanderburg,
David Berardo,
Xueying Guo,
Keivan Stassun,
Ian Crossfield,
David Charbonneau,
David W. Latham,
Abraham Loeb,
George Ricker,
Sara Seager,
Roland Vanderspek
Abstract:
Most known terrestrial planets orbit small stars with radii less than 60% that of the Sun. Theoretical models predict that these planets are more vulnerable to atmospheric loss than their counterparts orbiting Sun-like stars. To determine whether a thick atmosphere has survived on a small planet, one approach is to search for signatures of atmospheric heat redistribution in its thermal phase curve…
▽ More
Most known terrestrial planets orbit small stars with radii less than 60% that of the Sun. Theoretical models predict that these planets are more vulnerable to atmospheric loss than their counterparts orbiting Sun-like stars. To determine whether a thick atmosphere has survived on a small planet, one approach is to search for signatures of atmospheric heat redistribution in its thermal phase curve. Previous phase curve observations of the super-Earth 55 Cancri e (1.9 Earth radii) showed that its peak brightness is offset from the substellar point $-$ possibly indicative of atmospheric circulation. Here we report a phase curve measurement for the smaller, cooler planet LHS 3844b, a 1.3 Earth radius world in an 11-hour orbit around a small, nearby star. The observed phase variation is symmetric and has a large amplitude, implying a dayside brightness temperature of $1040\pm40$ kelvin and a nightside temperature consistent with zero kelvin (at one standard deviation). Thick atmospheres with surface pressures above 10 bar are ruled out by the data (at three standard deviations), and less-massive atmospheres are unstable to erosion by stellar wind. The data are well fitted by a bare rock model with a low Bond albedo (lower than 0.2 at two standard deviations). These results support theoretical predictions that hot terrestrial planets orbiting small stars may not retain substantial atmospheres.
△ Less
Submitted 19 August, 2019;
originally announced August 2019.
-
The Random Transiter -- EPIC 249706694/HD 139139
Authors:
S. Rappaport,
A. Vanderburg,
M. H. Kristiansen,
M. R. Omohundro,
H. M. Schwengeler,
I. A. Terentev,
F. Dai,
K. Masuda,
T. L. Jacobs,
D. LaCourse,
D. W. Latham,
A. Bieryla,
C. L. Hedges,
J. Dittmann,
G. Barentsen,
W. Cochran,
M. Endl,
J. M. Jenkins,
A. Mann
Abstract:
We have identified a star, EPIC 249706694 (HD 139139), that was observed during K2 Campaign 15 with the Kepler extended mission that appears to exhibit 28 transit-like events over the course of the 87-day observation. The unusual aspect of these dips, all but two of which have depths of $200 \pm 80$ ppm, is that they exhibit no periodicity, and their arrival times could just as well have been prod…
▽ More
We have identified a star, EPIC 249706694 (HD 139139), that was observed during K2 Campaign 15 with the Kepler extended mission that appears to exhibit 28 transit-like events over the course of the 87-day observation. The unusual aspect of these dips, all but two of which have depths of $200 \pm 80$ ppm, is that they exhibit no periodicity, and their arrival times could just as well have been produced by a random number generator. We show that no more than four of the events can be part of a periodic sequence. We have done a number of data quality tests to ascertain that these dips are of astrophysical origin, and while we cannot be absolutely certain that this is so, they have all the hallmarks of astrophysical variability on one of two possible host stars (a likely bound pair) in the photometric aperture. We explore a number of ideas for the origin of these dips, including actual planet transits due to multiple or dust emitting planets, anomalously large TTVs, S- and P-type transits in binary systems, a collection of dust-emitting asteroids, `dipper-star' activity, and short-lived starspots. All transit scenarios that we have been able to conjure up appear to fail, while the intrinsic stellar variability hypothesis would be novel and untested.
△ Less
Submitted 26 June, 2019;
originally announced June 2019.
-
Lyman-alpha in the GJ 1132 System: Stellar Emission and Planetary Atmospheric Evolution
Authors:
William C. Waalkes,
Zachory Berta-Thompson,
Vincent Bourrier,
Elisabeth Newton,
David Ehrenreich,
Eliza M. -R. Kempton,
David Charbonneau,
Jonathan Irwin,
Jason Dittmann
Abstract:
GJ 1132b, which orbits an M dwarf, is one of the few known Earth-sized planets, and at 12 pc away it is one of the closest known transiting planets. Receiving roughly 19x Earth's insolation, this planet is too hot to be habitable but can inform us about the volatile content of rocky planet atmospheres around cool stars. Using Hubble STIS spectra, we search for a transit in the Lyman-alpha line of…
▽ More
GJ 1132b, which orbits an M dwarf, is one of the few known Earth-sized planets, and at 12 pc away it is one of the closest known transiting planets. Receiving roughly 19x Earth's insolation, this planet is too hot to be habitable but can inform us about the volatile content of rocky planet atmospheres around cool stars. Using Hubble STIS spectra, we search for a transit in the Lyman-alpha line of neutral hydrogen (Ly-alpha). If we were to observe a deep Ly-alpha absorption signature, that would indicate the presence of a neutral hydrogen envelope flowing from GJ 1132b. On the other hand, ruling out deep absorption from neutral hydrogen may indicate that this planet does not have a detectable amount of hydrogen loss, is not losing hydrogen, or lost hydrogen and other volatiles early in the star's life. We do not detect a transit and determine a 2-sigma upper limit on the effective envelope radius of 0.36 R* in the red wing of the Ly-alpha line, which is the only portion of the spectrum we detect after absorption by the ISM. We analyze the Ly-alpha spectrum and stellar variability of GJ1132, which is a slowly-rotating 0.18 solar mass M dwarf with previously uncharacterized UV activity. Our data show stellar variabilities of 5-22%, which is consistent with the M dwarf UV variabilities of up to 41% found by \citet{Loyd2014}. Understanding the role that UV variability plays in planetary atmospheres is crucial to assess atmospheric evolution and the habitability of cooler rocky exoplanets.
△ Less
Submitted 10 June, 2019;
originally announced June 2019.
-
TESS Spots a Compact System of Super-Earths around the Naked-Eye Star HR 858
Authors:
Andrew Vanderburg,
Chelsea X. Huang,
Joseph E. Rodriguez,
Juliette C. Becker,
George R. Ricker,
Roland K. Vanderspek,
David W. Latham,
Sara Seager,
Joshua N. Winn,
Jon M. Jenkins,
Brett Addison,
Allyson Bieryla,
Cesar Briceño,
Brendan P. Bowler,
Timothy M. Brown,
Christopher J. Burke,
Jennifer A. Burt,
Douglas A. Caldwell,
Jake T. Clark,
Ian Crossfield,
Jason A. Dittmann,
Scott Dynes,
Benjamin J. Fulton,
Natalia Guerrero,
Daniel Harbeck
, et al. (26 additional authors not shown)
Abstract:
Transiting Exoplanet Survey Satellite (TESS) observations have revealed a compact multi-planet system around the sixth-magnitude star HR 858 (TIC 178155732, TOI 396), located 32 parsecs away. Three planets, each about twice the size of Earth, transit this slightly-evolved, late F-type star, which is also a member of a visual binary. Two of the planets may be in mean motion resonance. We analyze th…
▽ More
Transiting Exoplanet Survey Satellite (TESS) observations have revealed a compact multi-planet system around the sixth-magnitude star HR 858 (TIC 178155732, TOI 396), located 32 parsecs away. Three planets, each about twice the size of Earth, transit this slightly-evolved, late F-type star, which is also a member of a visual binary. Two of the planets may be in mean motion resonance. We analyze the TESS observations, using novel methods to model and remove instrumental systematic errors, and combine these data with follow-up observations taken from a suite of ground-based telescopes to characterize the planetary system. The HR 858 planets are enticing targets for precise radial velocity observations, secondary eclipse spectroscopy, and measurements of the Rossiter-McLaughlin effect.
△ Less
Submitted 15 July, 2019; v1 submitted 13 May, 2019;
originally announced May 2019.
-
Identifying Exoplanets with Deep Learning III: Automated Triage and Vetting of TESS Candidates
Authors:
Liang Yu,
Andrew Vanderburg,
Chelsea Huang,
Christopher J. Shallue,
Ian J. M. Crossfield,
B. Scott Gaudi,
Tansu Daylan,
Anne Dattilo,
David J. Armstrong,
George R. Ricker,
Roland K. Vanderspek,
David W. Latham,
Sara Seager,
Jason Dittmann,
John P. Doty,
Ana Glidden,
Samuel N. Quinn
Abstract:
NASA's Transiting Exoplanet Survey Satellite (TESS) presents us with an unprecedented volume of space-based photometric observations that must be analyzed in an efficient and unbiased manner. With at least $\sim1,000,000$ new light curves generated every month from full frame images alone, automated planet candidate identification has become an attractive alternative to human vetting. Here we pres…
▽ More
NASA's Transiting Exoplanet Survey Satellite (TESS) presents us with an unprecedented volume of space-based photometric observations that must be analyzed in an efficient and unbiased manner. With at least $\sim1,000,000$ new light curves generated every month from full frame images alone, automated planet candidate identification has become an attractive alternative to human vetting. Here we present a deep learning model capable of performing triage and vetting on TESS candidates. Our model is modified from an existing neural network designed to automatically classify Kepler candidates, and is the first neural network to be trained and tested on real TESS data. In triage mode, our model can distinguish transit-like signals (planet candidates and eclipsing binaries) from stellar variability and instrumental noise with an average precision (the weighted mean of precisions over all classification thresholds) of 97.0% and an accuracy of 97.4%. In vetting mode, the model is trained to identify only planet candidates with the help of newly added scientific domain knowledge, and achieves an average precision of 69.3% and an accuracy of 97.8%. We apply our model on new data from Sector 6, and present 288 new signals that received the highest scores in triage and vetting and were also identified as planet candidates by human vetters. We also provide a homogeneously classified set of TESS candidates suitable for future training.
△ Less
Submitted 14 May, 2019; v1 submitted 4 April, 2019;
originally announced April 2019.