-
ZTF SN~Ia DR2: Cosmology-independent constraints on Type Ia supernova standardisation from supernova siblings
Authors:
S. Dhawan,
E. Mortsell,
J. Johansson,
A. Goobar,
M. Rigault,
M. Smith,
K. Maguire,
J. Nordin,
G. Dimitriadis,
P. E. Nugent,
L. Galbany,
J. Sollerman,
T. de Jaeger,
J. H. Terwel,
Y. -L. Kim,
Umut Burgaz,
G. Helou,
J. Purdum,
S. L. Groom,
R. Laher,
B. Healy
Abstract:
Understanding Type Ia supernovae (SNe~Ia) and the empirical standardisation relations that make them excellent distance indicators is vital to improving cosmological constraints. SN~Ia ``siblings", i.e. two or more SNe~Ia in the same host or parent galaxy offer a unique way to infer the standardisation relations and their diversity across the population. We analyse a sample of 25 SN~Ia pairs, obse…
▽ More
Understanding Type Ia supernovae (SNe~Ia) and the empirical standardisation relations that make them excellent distance indicators is vital to improving cosmological constraints. SN~Ia ``siblings", i.e. two or more SNe~Ia in the same host or parent galaxy offer a unique way to infer the standardisation relations and their diversity across the population. We analyse a sample of 25 SN~Ia pairs, observed homogeneously by the Zwicky Transient Factory (ZTF) to infer the SNe~Ia light curve width-luminosity and colour-luminosity parameters $α$ and $β$. Using the pairwise constraints from siblings, allowing for a diversity in the standardisation relations, we find $α= 0.218 \pm 0.055 $ and $β= 3.084 \pm 0.312$, respectively, with a dispersion in $α$ and $β$ of $\leq 0.195$ and $\leq 0.923$, respectively, at 95$\%$ C.L. While the median dispersion is large, the values within $\sim 1 σ$ are consistent with no dispersion. Hence, fitting for a single global standardisation relation, we find $α= 0.228 \pm 0.029 $ and $β= 3.160 \pm 0.191$. We find a very small intrinsic scatter of the siblings sample $σ_{\rm int} \leq 0.10$ at 95\% C.L. compared to $σ_{\rm int} = 0.22 \pm 0.04$ when computing the scatter using the Hubble residuals without comparing them as siblings. Splitting the sample based on host galaxy stellar mass, we find that SNe~Ia in both subsamples have consistent $α$ and $β$. The $β$ value is consistent with the value for the cosmological sample. However, we find a higher $α$ by $\sim 2.5 - 3.5 σ$. The high $α$ is driven by low $x_1$ pairs, potentially suggesting that the slow and fast declining SN~Ia have different slopes of the width-luminosity relation. We can confirm or refute this with increased statistics from near future time-domain surveys. (abridged)
△ Less
Submitted 3 June, 2024;
originally announced June 2024.
-
Searching for gravitational wave optical counterparts with the Zwicky Transient Facility: summary of O4a
Authors:
Tomás Ahumada,
Shreya Anand,
Michael W. Coughlin,
Vaidehi Gupta,
Mansi M. Kasliwal,
Viraj R. Karambelkar,
Robert D. Stein,
Gaurav Waratkar,
Vishwajeet Swain,
Theophile Jegou du Laz,
Akash Anumarlapudi,
Igor Andreoni,
Mattia Bulla,
Gokul P. Srinivasaragavan,
Andrew Toivonen,
Avery Wold,
Eric C. Bellm,
S. Bradley Cenko,
David L. Kaplan,
Jesper Sollerman,
Varun Bhalerao,
Daniel Perley,
Anirudh Salgundi,
Aswin Suresh,
K-Ryan Hinds
, et al. (27 additional authors not shown)
Abstract:
During the first half of the fourth observing run (O4a) of the International Gravitational Wave Network (IGWN), the Zwicky Transient Facility (ZTF) conducted a systematic search for kilonova (KN) counterparts to binary neutron star (BNS) and neutron star-black hole (NSBH) merger candidates. Here, we present a comprehensive study of the five high-significance (FAR < 1 per year) BNS and NSBH candida…
▽ More
During the first half of the fourth observing run (O4a) of the International Gravitational Wave Network (IGWN), the Zwicky Transient Facility (ZTF) conducted a systematic search for kilonova (KN) counterparts to binary neutron star (BNS) and neutron star-black hole (NSBH) merger candidates. Here, we present a comprehensive study of the five high-significance (FAR < 1 per year) BNS and NSBH candidates in O4a. Our follow-up campaigns relied on both target-of-opportunity observations (ToO) and re-weighting of the nominal survey schedule to maximize coverage. We describe the toolkit we have been developing, Fritz, an instance of SkyPortal, instrumental in coordinating and managing our telescope scheduling, candidate vetting, and follow-up observations through a user-friendly interface. ZTF covered a total of 2841 deg$^2$ within the skymaps of the high-significance GW events, reaching a median depth of g~20.2 mag. We circulated 15 candidates, but found no viable KN counterpart to any of the GW events. Based on the ZTF non-detections of the high-significance events in O4a, we used a Bayesian approach, nimbus, to quantify the posterior probability of KN model parameters that are consistent with our non-detections. Our analysis favors KNe with initial absolute magnitude fainter than -16 mag. The joint posterior probability of a GW170817-like KN associated with all our O4a follow-ups was 64%. Additionally, we use a survey simulation software, simsurvey, to determine that our combined filtered efficiency to detect a GW170817-like KN is 36%, when considering the 5 confirmed astrophysical events in O3 (1 BNS and 4 NSBH), along with our O4a follow-ups. Following Kasliwal et al. (2020), we derived joint constraints on the underlying KN luminosity function based on our O3 and O4a follow-ups, determining that no more than 76% of KNe fading at 1 mag/day can peak at a magnitude brighter than -17.5 mag.
△ Less
Submitted 20 May, 2024;
originally announced May 2024.
-
An overview of existing and new nuclear and astrophysical constraints on the equation of state of neutron-rich dense matter
Authors:
Hauke Koehn,
Henrik Rose,
Peter T. H. Pang,
Rahul Somasundaram,
Brendan T. Reed,
Ingo Tews,
Adrian Abac,
Oleg Komoltsev,
Nina Kunert,
Aleksi Kurkela,
Michael W. Coughlin,
Brian F. Healy,
Tim Dietrich
Abstract:
Through continuous progress in nuclear theory and experiment and an increasing number of neutron-star observations, a multitude of information about the equation of state (EOS) for matter at extreme densities is available. Here, we apply these different pieces of data individually to a broad set of physics-agnostic candidate EOSs and analyze the resulting constraints. Specifically, we make use of…
▽ More
Through continuous progress in nuclear theory and experiment and an increasing number of neutron-star observations, a multitude of information about the equation of state (EOS) for matter at extreme densities is available. Here, we apply these different pieces of data individually to a broad set of physics-agnostic candidate EOSs and analyze the resulting constraints. Specifically, we make use of information from chiral effective field theory, perturbative quantum chromodynamics, as well as data from heavy-ion collisions and the PREX-II and CREX experiments. We also investigate the impact of current mass and radius measurements of neutron stars, such as radio timing measurements of heavy pulsars, NICER data, and other X-ray observations. We augment these by reanalyses of the gravitational-wave (GW) signal GW170817, its associated kilonova AT2017gfo and gamma-ray burst afterglow, the GW signal GW190425, and the GRB211211A afterglow, where we use improved models for the tidal waveform and kilonova light curves. Additionally, we consider the postmerger fate of GW170817 and its consequences for the EOS. This large and diverse set of constraints is eventually combined in numerous ways to explore limits on quantities such as the typical neutron-star radius, the maximum neutron-star mass, the nuclear symmetry-energy parameters, and the speed of sound. Based on the priors from our EOS candidate set, we find the radius of the canonical 1.4 M$_\odot$ neutron star to be $R_{1.4}= 12.27_{-0.94}^{+0.83}$ km and the TOV mass $M_{\rm TOV}= 2.26_{-0.22}^{+0.45}$ M$_\odot$ at 95% credibility, when including those constraints where systematic uncertainties are deemed small. A less conservative approach, combining all the presented constraints, similarly yields $R_{1.4}= 12.20_{-0.50}^{+0.53}$ km and $M_{\rm TOV}= 2.31_{-0.20}^{+0.08}$ M$_\odot$.
△ Less
Submitted 6 February, 2024;
originally announced February 2024.
-
Spectroscopic observations of progenitor activity 100 days before a Type Ibn supernova
Authors:
S. J. Brennan,
J. Sollerman,
I. Irani,
S. Schulze,
P. Chen,
K. K. Das,
K. De,
C. Fransson,
A. Gal-Yam,
A. Gkini,
K. R. Hinds,
R. Lunnan,
D. Perley,
YJ. Qin,
R. Stein,
J. Wise,
L. Yan,
E. A. Zimmerman,
S. Anand,
R. J. Bruch,
R. Dekany,
A. J. Drake,
C. Fremling,
B. Healy,
V. Karambelkar
, et al. (8 additional authors not shown)
Abstract:
Obtaining spectroscopic observations of the progenitors of core-collapse supernovae is often unfeasible due to an inherent lack of knowledge as to which stars will go supernova and when they will explode. In this letter, we present photometric and spectroscopic observations of the progenitor activity of SN 2023fyq in the preceding 150 days before the He-rich progenitor exploded as a Type Ibn super…
▽ More
Obtaining spectroscopic observations of the progenitors of core-collapse supernovae is often unfeasible due to an inherent lack of knowledge as to which stars will go supernova and when they will explode. In this letter, we present photometric and spectroscopic observations of the progenitor activity of SN 2023fyq in the preceding 150 days before the He-rich progenitor exploded as a Type Ibn supernova. The progenitor of SN 2023fyq shows an exponential rise in flux prior to core-collapse. Complex He I emission line features are observed, with a P-Cygni like profile, as well as an evolving broad base with velocities on the order of 10,000 km/s, possibly due to electron scattering. The luminosity and evolution of SN 2023fyq are consistent with a faint Type Ibn, reaching a peak r-band magnitude of 18.1 mag, although there is some uncertainty in the distance to the host, NGC 4388, located in the Virgo cluster. We present additional evidence of asymmetric He-rich material being present prior to the explosion of SN 2023fyq, as well as after, suggesting this material has survived the ejecta-CSM interaction. Broad [O I] and the Ca II triplet lines are observed at late phases, confirming that SN 2023fyq was a genuine supernova rather than a non-terminal interacting transient. SN 2023fyq provides insight into the final moments of a massive star's life, highlighting that the progenitor is likely highly unstable before core-collapse.
△ Less
Submitted 25 March, 2024; v1 submitted 26 January, 2024;
originally announced January 2024.
-
Quantum-centric Supercomputing for Materials Science: A Perspective on Challenges and Future Directions
Authors:
Yuri Alexeev,
Maximilian Amsler,
Paul Baity,
Marco Antonio Barroca,
Sanzio Bassini,
Torey Battelle,
Daan Camps,
David Casanova,
Young jai Choi,
Frederic T. Chong,
Charles Chung,
Chris Codella,
Antonio D. Corcoles,
James Cruise,
Alberto Di Meglio,
Jonathan Dubois,
Ivan Duran,
Thomas Eckl,
Sophia Economou,
Stephan Eidenbenz,
Bruce Elmegreen,
Clyde Fare,
Ismael Faro,
Cristina Sanz Fernández,
Rodrigo Neumann Barros Ferreira
, et al. (102 additional authors not shown)
Abstract:
Computational models are an essential tool for the design, characterization, and discovery of novel materials. Hard computational tasks in materials science stretch the limits of existing high-performance supercomputing centers, consuming much of their simulation, analysis, and data resources. Quantum computing, on the other hand, is an emerging technology with the potential to accelerate many of…
▽ More
Computational models are an essential tool for the design, characterization, and discovery of novel materials. Hard computational tasks in materials science stretch the limits of existing high-performance supercomputing centers, consuming much of their simulation, analysis, and data resources. Quantum computing, on the other hand, is an emerging technology with the potential to accelerate many of the computational tasks needed for materials science. In order to do that, the quantum technology must interact with conventional high-performance computing in several ways: approximate results validation, identification of hard problems, and synergies in quantum-centric supercomputing. In this paper, we provide a perspective on how quantum-centric supercomputing can help address critical computational problems in materials science, the challenges to face in order to solve representative use cases, and new suggested directions.
△ Less
Submitted 14 December, 2023;
originally announced December 2023.
-
The ZTF Source Classification Project: III. A Catalog of Variable Sources
Authors:
Brian F. Healy,
Michael W. Coughlin,
Ashish A. Mahabal,
Theophile Jegou du Laz,
Andrew Drake,
Matthew J. Graham,
Lynne A. Hillenbrand,
Jan van Roestel,
Paula Szkody,
LeighAnna Zielske,
Mohammed Guiga,
Muhammad Yusuf Hassan,
Jill L. Hughes,
Guy Nir,
Saagar Parikh,
Sungmin Park,
Palak Purohit,
Umaa Rebbapragada,
Draco Reed,
Daniel Warshofsky,
Avery Wold,
Joshua S. Bloom,
Frank J. Masci,
Reed Riddle,
Roger Smith
Abstract:
The classification of variable objects provides insight into a wide variety of astrophysics ranging from stellar interiors to galactic nuclei. The Zwicky Transient Facility (ZTF) provides time series observations that record the variability of more than a billion sources. The scale of these data necessitates automated approaches to make a thorough analysis. Building on previous work, this paper re…
▽ More
The classification of variable objects provides insight into a wide variety of astrophysics ranging from stellar interiors to galactic nuclei. The Zwicky Transient Facility (ZTF) provides time series observations that record the variability of more than a billion sources. The scale of these data necessitates automated approaches to make a thorough analysis. Building on previous work, this paper reports the results of the ZTF Source Classification Project (SCoPe), which trains neural network and XGBoost machine learning (ML) algorithms to perform dichotomous classification of variable ZTF sources using a manually constructed training set containing 170,632 light curves. We find that several classifiers achieve high precision and recall scores, suggesting the reliability of their predictions for 209,991,147 light curves across 77 ZTF fields. We also identify the most important features for XGB classification and compare the performance of the two ML algorithms, finding a pattern of higher precision among XGB classifiers. The resulting classification catalog is available to the public, and the software developed for SCoPe is open-source and adaptable to future time-domain surveys.
△ Less
Submitted 11 April, 2024; v1 submitted 30 November, 2023;
originally announced December 2023.
-
Insight into the Formation of $β$ Pic b through the Composition of Its Parent Protoplanetary Disk as Revealed by the $β$ Pic Moving Group Member HD 181327
Authors:
Henrique Reggiani,
Jhon Yana Galarza,
Kevin C. Schlaufman,
David K. Sing,
Brian F. Healy,
Andrew McWilliam,
Joshua D. Lothringer,
Laurent Pueyo
Abstract:
It has been suggested that $β$ Pic b has a supersolar metallicity and subsolar C/O ratio. Assuming solar carbon and oxygen abundances for the star $β$ Pic and therefore the planet's parent protoplanetary disk, $β$ Pic b's C/O ratio suggests that it formed via core accretion between its parent protoplanteary disk's H$_{2}$O and CO$_{2}$ ice lines. $β$ Pic b's high metallicity is difficult to reconc…
▽ More
It has been suggested that $β$ Pic b has a supersolar metallicity and subsolar C/O ratio. Assuming solar carbon and oxygen abundances for the star $β$ Pic and therefore the planet's parent protoplanetary disk, $β$ Pic b's C/O ratio suggests that it formed via core accretion between its parent protoplanteary disk's H$_{2}$O and CO$_{2}$ ice lines. $β$ Pic b's high metallicity is difficult to reconcile with its mass $M_{\text{p}}~=~11.7~M_{\text{Jup}}$ though. Massive stars can present peculiar photospheric abundances that are unlikely to record the abundances of their former protoplanetary disks. This issue can be overcome for early-type stars in moving groups by infering the elemental abundances of the FGK stars in the same moving group that formed in the same molecular cloud and presumably share the same composition. We infer the photospheric abundances of the F dwarf HD 181327, a $β$ Pic moving group member that is the best available proxy for the composition of $β$ Pic b's parent protoplanetary disk. In parallel, we infer updated atmospheric abundances for $β$ Pic b. As expected for a planet of its mass formed via core-accretion beyond its parent protoplanetary disk's H$_{2}$O ice line, we find that $β$ Pic b's atmosphere is consistent with stellar metallicity and confirm that is has superstellar carbon and oxygen abundances with a substellar C/O ratio. We propose that the elemental abundances of FGK dwarfs in moving groups can be used as proxies for the otherwise difficult-to-infer elemental abundances of early-type and late-type members of the same moving groups.
△ Less
Submitted 20 November, 2023;
originally announced November 2023.
-
Photometric prioritization of neutron star merger candidates
Authors:
E. O. Ofek,
N L. Strotjohann,
I. Arcavi,
A. Gal-Yam,
D. Kushnir,
E. Waxman,
M. M. Kasliwal,
A. Drake,
M. Graham,
J. Purdum,
B. Rusholme,
Y. Sharma,
R. Smith,
A. Wold,
B. F. Healy
Abstract:
Rapid identification of the optical counterparts of Neutron Star (NS) merger events discovered by gravitational wave detectors may require observing a large error region and sifting through a large number of transients to identify the object of interest. Given the expense of spectroscopic observations, a question arises: How can we utilize photometric observations for candidate prioritization, and…
▽ More
Rapid identification of the optical counterparts of Neutron Star (NS) merger events discovered by gravitational wave detectors may require observing a large error region and sifting through a large number of transients to identify the object of interest. Given the expense of spectroscopic observations, a question arises: How can we utilize photometric observations for candidate prioritization, and what kinds of photometric observations are needed to achieve this goal? NS merger kilonova exhibits low ejecta mass (~5x10^-2 solar mass) and a rapidly evolving photospheric radius (with a velocity ~0.2c). As a consequence, these sources display rapid optical-flux evolution. Indeed, selection based on fast flux variations is commonly used for young supernovae and NS mergers. In this study, we leverage the best currently available flux-limited transient survey - the Zwicky Transient Facility Bright Transient Survey - to extend and quantify this approach. We focus on selecting transients detected in a 3-day cadence survey and observed at a one-day cadence. We explore their distribution in the phase space defined by g-r, g-dot, and r-dot. Our analysis demonstrates that for a significant portion of the time during the first week, the kilonova AT 2017gfo stands out in this phase space. It is important to note that this investigation is subject to various biases and challenges; nevertheless, it suggests that certain photometric observations can be leveraged to identify transients with the highest probability of being fast-evolving events. We also find that a large fraction (~0.75) of the transient candidates with |g-dot|>0.7 mag/day, are cataclysmic variables or active galactic nuclei with radio counterparts.
△ Less
Submitted 8 November, 2023;
originally announced November 2023.
-
Feature Aggregation in Joint Sound Classification and Localization Neural Networks
Authors:
Brendan Healy,
Patrick McNamee,
Zahra Nili Ahmadabadi
Abstract:
This study addresses the application of deep learning techniques in joint sound signal classification and localization networks. Current state-of-the-art sound source localization deep learning networks lack feature aggregation within their architecture. Feature aggregation enhances model performance by enabling the consolidation of information from different feature scales, thereby improving feat…
▽ More
This study addresses the application of deep learning techniques in joint sound signal classification and localization networks. Current state-of-the-art sound source localization deep learning networks lack feature aggregation within their architecture. Feature aggregation enhances model performance by enabling the consolidation of information from different feature scales, thereby improving feature robustness and invariance. This is particularly important in SSL networks, which must differentiate direct and indirect acoustic signals. To address this gap, we adapt feature aggregation techniques from computer vision neural networks to signal detection neural networks. Additionally, we propose the Scale Encoding Network (SEN) for feature aggregation to encode features from various scales, compressing the network for more computationally efficient aggregation. To evaluate the efficacy of feature aggregation in SSL networks, we integrated the following computer vision feature aggregation sub-architectures into a SSL control architecture: Path Aggregation Network (PANet), Weighted Bi-directional Feature Pyramid Network (BiFPN), and SEN. These sub-architectures were evaluated using two metrics for signal classification and two metrics for direction-of-arrival regression. PANet and BiFPN are established aggregators in computer vision models, while the proposed SEN is a more compact aggregator. The results suggest that models incorporating feature aggregations outperformed the control model, the Sound Event Localization and Detection network (SELDnet), in both sound signal classification and localization. The feature aggregation techniques enhance the performance of sound detection neural networks, particularly in direction-of-arrival regression.
△ Less
Submitted 27 January, 2024; v1 submitted 29 October, 2023;
originally announced October 2023.
-
Multi-band analyses of the bright GRB 230812B and the associated SN2023pel
Authors:
T. Hussenot-Desenonges,
T. Wouters,
N. Guessoum,
I. Abdi,
A. Abulwfa,
C. Adami,
J. F. Agüí Fernández,
T. Ahumada,
V. Aivazyan,
D. Akl,
S. Anand,
C. M. Andrade,
S. Antier,
S. A. Ata,
P. D'Avanzo,
Y. A. Azzam,
A. Baransky,
S. Basa,
M. Blazek,
P. Bendjoya,
S. Beradze,
P. Boumis,
M. Bremer,
R. Brivio,
V. Buat
, et al. (87 additional authors not shown)
Abstract:
GRB~230812B is a bright and relatively nearby ($z =0.36$) long gamma-ray burst (GRB) that has generated significant interest in the community and has thus been observed over the entire electromagnetic spectrum. We report over 80 observations in X-ray, ultraviolet, optical, infrared, and sub-millimeter bands from the GRANDMA (Global Rapid Advanced Network for Multi-messenger Addicts) network of obs…
▽ More
GRB~230812B is a bright and relatively nearby ($z =0.36$) long gamma-ray burst (GRB) that has generated significant interest in the community and has thus been observed over the entire electromagnetic spectrum. We report over 80 observations in X-ray, ultraviolet, optical, infrared, and sub-millimeter bands from the GRANDMA (Global Rapid Advanced Network for Multi-messenger Addicts) network of observatories and from observational partners. Adding complementary data from the literature, we then derive essential physical parameters associated with the ejecta and external properties (i.e. the geometry and environment) of the GRB and compare with other analyses of this event. We spectroscopically confirm the presence of an associated supernova, SN2023pel, and we derive a photospheric expansion velocity of v $\sim$ 17$\times10^3$ km s$^{-1}$. We analyze the photometric data first using empirical fits of the flux and then with full Bayesian Inference. We again strongly establish the presence of a supernova in the data, with a maximum (pseudo-)bolometric luminosity of $5.75 \times 10^{42}$ erg/s, at $15.76^{+0.81}_{-1.21}$ days (in the observer frame) after the trigger, with a half-max time width of 22.0 days. We compare these values with those of SN1998bw, SN2006aj, and SN2013dx. Our best-fit model favours a very low density environment ($\log_{10}({n_{\rm ISM}/{\rm cm}^{-3}}) = -2.38^{+1.45}_{-1.60}$) and small values for the jet's core angle $θ_{\rm core} = 1.54^{+1.02}_{-0.81} \ \rm{deg}$ and viewing angle $θ_{\rm obs} = 0.76^{+1.29}_{-0.76} \ \rm{deg}$. GRB 230812B is thus one of the best observed afterglows with a distinctive supernova bump.
△ Less
Submitted 17 February, 2024; v1 submitted 22 October, 2023;
originally announced October 2023.
-
An Optically-Discovered Outburst from XTE J1859+226
Authors:
Eric C. Bellm,
Yuankun Wang,
Jan van Roestel,
Rebecca A. Phillipson,
Michael W. Coughlin,
John A. Tomsick,
Steven L. Groom,
Brian Healy,
Josiah Purdum,
Ben Rusholme,
Jesper Sollerman,
Peter Bealo,
Stefano Lora,
Eddy Muyllaert,
Ivo Peretto,
Erik J. Schwendeman
Abstract:
Using the Zwicky Transient Facility, in 2021 February we identified the first known outburst of the Black Hole X-ray Transient XTE J1859+226 since its discovery in 1999. The outburst was visible at X-ray, UV, and optical wavelengths for less than 20 days, substantially shorter than its 320-day full outburst in 1999, and the observed peak luminosity was two orders of magnitude lower. Its peak bolom…
▽ More
Using the Zwicky Transient Facility, in 2021 February we identified the first known outburst of the Black Hole X-ray Transient XTE J1859+226 since its discovery in 1999. The outburst was visible at X-ray, UV, and optical wavelengths for less than 20 days, substantially shorter than its 320-day full outburst in 1999, and the observed peak luminosity was two orders of magnitude lower. Its peak bolometric luminosity was only $2\times 10^{35}$ erg s$^{-1}$, implying an Eddington fraction of about $3\times10^{-4}$. The source remained in the hard spectral state throughout the outburst. From optical spectroscopy measurements we estimate an outer disk radius of 10$^{11}$ cm. The low observed X-ray luminosity is not sufficient to irradiate the entire disk, but we observe a surprising exponential decline in the X-ray lightcurve. These observations highlight the potential of optical and infrared (O/IR) synoptic surveys to discover low-luminosity activity from X-ray transients.
△ Less
Submitted 19 September, 2023;
originally announced September 2023.
-
Chemical Distribution of the Dynamical Ejecta in the Neutron Star Merger GW170817
Authors:
Shreya Anand,
Peter T. H. Pang,
Mattia Bulla,
Michael W. Coughlin,
Tim Dietrich,
Brian Healy,
Thomas Hussenot-Desenonges,
Theophile Jegou du Laz,
Mansi M. Kasliwal,
Nina Kunert,
Ivan Markin,
Kunal Mooley,
Vsevolod Nedora,
Anna Neuweiler
Abstract:
GW170817 and its associated electromagnetic counterpart AT2017gfo continue to be a treasure trove as observations and modeling continue. Recent precision astrometry of AT2017gfo with the Hubble Space Telescope combined with previous constraints from Very Long Baseline Interferometry (VLBI) constraints narrowed down the inclination angle to 19-25 deg (90\% confidence). This paper explores how the i…
▽ More
GW170817 and its associated electromagnetic counterpart AT2017gfo continue to be a treasure trove as observations and modeling continue. Recent precision astrometry of AT2017gfo with the Hubble Space Telescope combined with previous constraints from Very Long Baseline Interferometry (VLBI) constraints narrowed down the inclination angle to 19-25 deg (90\% confidence). This paper explores how the inclusion of precise inclination information can reveal new insights about the ejecta properties, in particular, about the composition of the dynamical ejecta of AT2017gfo. Our analysis relies on updated kilonova modeling, which includes state-of-the-art heating rates, thermalization efficiencies, and opacities and is parameterized by $\bar{Y}_{\rm e,dyn}$, the average electron fraction of the dynamical ejecta component. Using this model, we incorporate the latest inclination angle constraint of AT2017gfo into a light curve fitting framework to derive updated parameter estimates. Our results suggest that the viewing angle of the observer is pointed towards the lanthanide-poor ($Y_{\rm e,dyn}\gtrsim0.25$), squeezed polar dynamical ejecta component, which can explain the early blue emission observed in the light curve of AT2017gfo. In contrast to a recent claim of spherical ejecta powering AT2017gfo, our study indicates that the composition of the dynamical ejecta has a strong angular dependence, with a lanthanide-rich ($Y_{\rm e,dyn}\lesssim0.25$), tidal component distributed around the merger plane with a half-opening angle of $35^\circ$. The inclination angle constraint reduces $\bar{Y}_{\rm e,dyn}$ from $0.24$ to $0.22$, with values $0.17\lesssim Y_{\rm e, dyn} \lesssim0.41$ enabling the robust production of $r$-process elements up to the $3^{\rm rd}$ peak in the tidal dynamical ejecta.
△ Less
Submitted 20 July, 2023;
originally announced July 2023.
-
Encoding a magic state with beyond break-even fidelity
Authors:
Riddhi S. Gupta,
Neereja Sundaresan,
Thomas Alexander,
Christopher J. Wood,
Seth T. Merkel,
Michael B. Healy,
Marius Hillenbrand,
Tomas Jochym-O'Connor,
James R. Wootton,
Theodore J. Yoder,
Andrew W. Cross,
Maika Takita,
Benjamin J. Brown
Abstract:
To run large-scale algorithms on a quantum computer, error-correcting codes must be able to perform a fundamental set of operations, called logic gates, while isolating the encoded information from noise~\cite{Harper2019,Ryan-Anderson2021,Egan2021fault, Chen2022calibrated, Sundaresan2022matching, ryananderson2022implementing, Postler2022demonstration, GoogleAI2023}. We can complete a universal set…
▽ More
To run large-scale algorithms on a quantum computer, error-correcting codes must be able to perform a fundamental set of operations, called logic gates, while isolating the encoded information from noise~\cite{Harper2019,Ryan-Anderson2021,Egan2021fault, Chen2022calibrated, Sundaresan2022matching, ryananderson2022implementing, Postler2022demonstration, GoogleAI2023}. We can complete a universal set of logic gates by producing special resources called magic states~\cite{Bravyi2005universal,Maier2013magic, Chamberland2022building}. It is therefore important to produce high-fidelity magic states to conduct algorithms while introducing a minimal amount of noise to the computation. Here, we propose and implement a scheme to prepare a magic state on a superconducting qubit array using error correction. We find that our scheme produces better magic states than those we can prepare using the individual qubits of the device. This demonstrates a fundamental principle of fault-tolerant quantum computing~\cite{Shor96}, namely, that we can use error correction to improve the quality of logic gates with noisy qubits. Additionally, we show we can increase the yield of magic states using adaptive circuits, where circuit elements are changed depending on the outcome of mid-circuit measurements. This demonstrates an essential capability we will need for many error-correction subroutines. Our prototype will be invaluable in the future as it can reduce the number of physical qubits needed to produce high-fidelity magic states in large-scale quantum-computing architectures.
△ Less
Submitted 13 March, 2024; v1 submitted 22 May, 2023;
originally announced May 2023.
-
A data science platform to enable time-domain astronomy
Authors:
Michael W. Coughlin,
Joshua S. Bloom,
Guy Nir,
Sarah Antier,
Theophile Jegou du Laz,
Stéfan van der Walt,
Arien Crellin-Quick,
Thomas Culino,
Dmitry A. Duev,
Daniel A. Goldstein,
Brian F. Healy,
Viraj Karambelkar,
Jada Lilleboe,
Kyung Min Shin,
Leo P. Singer,
Tomas Ahumada,
Shreya Anand,
Eric C. Bellm,
Richard Dekany,
Matthew J. Graham,
Mansi M. Kasliwal,
Ivona Kostadinova,
R. Weizmann Kiendrebeogo,
Shrinivas R. Kulkarni,
Sydney Jenkins
, et al. (28 additional authors not shown)
Abstract:
SkyPortal is an open-source software package designed to efficiently discover interesting transients, manage follow-up, perform characterization, and visualize the results. By enabling fast access to archival and catalog data, cross-matching heterogeneous data streams, and the triggering and monitoring of on-demand observations for further characterization, a SkyPortal-based platform has been oper…
▽ More
SkyPortal is an open-source software package designed to efficiently discover interesting transients, manage follow-up, perform characterization, and visualize the results. By enabling fast access to archival and catalog data, cross-matching heterogeneous data streams, and the triggering and monitoring of on-demand observations for further characterization, a SkyPortal-based platform has been operating at scale for 2 yr for the Zwicky Transient Facility Phase II community, with hundreds of users, containing tens of millions of time-domain sources, interacting with dozens of telescopes, and enabling community reporting. While SkyPortal emphasizes rich user experiences (UX) across common frontend workflows, recognizing that scientific inquiry is increasingly performed programmatically, SkyPortal also surfaces an extensive and well-documented API system. From backend and frontend software to data science analysis tools and visualization frameworks, the SkyPortal design emphasizes the re-use and leveraging of best-in-class approaches, with a strong extensibility ethos. For instance, SkyPortal now leverages ChatGPT large-language models (LLMs) to automatically generate and surface source-level human-readable summaries. With the imminent re-start of the next-generation of gravitational wave detectors, SkyPortal now also includes dedicated multi-messenger features addressing the requirements of rapid multi-messenger follow-up: multi-telescope management, team/group organizing interfaces, and cross-matching of multi-messenger data streams with time-domain optical surveys, with interfaces sufficiently intuitive for the newcomers to the field. (abridged)
△ Less
Submitted 14 June, 2023; v1 submitted 28 April, 2023;
originally announced May 2023.
-
Zwicky Transient Facility and Globular Clusters: The Period-Luminosity and Period-Wesenheit Relations for SX Phoenicis Variables in the gri-Band
Authors:
Chow-Choong Ngeow,
Anupam Bhardwaj,
Matthew J. Graham,
Brian F. Healy,
Russ R. Laher,
Reed Riddle,
Avery Wold
Abstract:
SX Phoenicis (SXP) variables are short period pulsating stars that exhibit a period-luminosity (PL) relation. We derived the gri-band PL and extinction-free period-Wesenheit (PW) relations, as well as the period-color (PC) and reddening-free period-Q-index (PQ) relations for 47 SXP stars in located in 21 globular clusters using the optical light curves taken from Zwicky Transient Facility (ZTF). T…
▽ More
SX Phoenicis (SXP) variables are short period pulsating stars that exhibit a period-luminosity (PL) relation. We derived the gri-band PL and extinction-free period-Wesenheit (PW) relations, as well as the period-color (PC) and reddening-free period-Q-index (PQ) relations for 47 SXP stars in located in 21 globular clusters using the optical light curves taken from Zwicky Transient Facility (ZTF). These empirically relations were derived for the first time in the gri filters except for the g-band PL relation. We used our gi band PL and PW relations to derive a distance modulus to Crater II dwarf spheroidal which hosts one SXP variable. Assuming that the fundamental and first-overtone pulsation mode for the SXP variable in Crater II, we found distance moduli of $20.03 \pm 0.23$ mag and $20.37 \pm 0.24$ mag, respectively, using the PW relation, where the latter is in excellent agreement with independent RR Lyrae based distance to Crater II dwarf galaxy.
△ Less
Submitted 13 March, 2023;
originally announced March 2023.
-
A Study of Stellar Spins in 15 Open Clusters
Authors:
Brian F. Healy,
P. R. McCullough,
Kevin C. Schlaufman,
Geza Kovacs
Abstract:
We analyze spectroscopic and photometric data to determine the projected inclinations of stars in 11 open clusters, placing constraints on the spin-axis distributions of six clusters. We combine these results with four additional clusters studied by Healy & McCullough (2020) and Healy et al. (2021) to perform an ensemble analysis of their spins. We find that eight out of ten constrained clusters (…
▽ More
We analyze spectroscopic and photometric data to determine the projected inclinations of stars in 11 open clusters, placing constraints on the spin-axis distributions of six clusters. We combine these results with four additional clusters studied by Healy & McCullough (2020) and Healy et al. (2021) to perform an ensemble analysis of their spins. We find that eight out of ten constrained clusters (80%) have spin-axis orientations consistent with isotropy, and we establish a lower limit of four out of ten (40%) isotropic clusters at 75% confidence, assuming no correlation of spins between clusters. We also identify two clusters whose spin-axis distributions can be better described by a model consisting of an aligned fraction of stars combined with an isotropic distribution. However, the inclination values of these stars may be influenced by systematic error, and the small number of stars modeled as aligned in these two clusters precludes the interpretation that their stellar subsets are physically aligned. Overall, no cluster displays an unambiguous signature of spin alignment, and 97% of the stars in our sample are consistent with isotropic orientations in their respective clusters. Our results offer support for the dominance of turbulence over ordered rotation in clumps and do not suggest alignment of rotation axes and magnetic fields in protostars.
△ Less
Submitted 24 January, 2023;
originally announced January 2023.
-
An Optimal Control Strategy for Execution of Large Stock Orders Using LSTMs
Authors:
A. Papanicolaou,
H. Fu,
P. Krishnamurthy,
B. Healy,
F. Khorrami
Abstract:
In this paper, we simulate the execution of a large stock order with real data and general power law in the Almgren and Chriss model. The example that we consider is the liquidation of a large position executed over the course of a single trading day in a limit order book. Transaction costs are incurred because large orders walk the order book, that is, they consume order book liquidity beyond the…
▽ More
In this paper, we simulate the execution of a large stock order with real data and general power law in the Almgren and Chriss model. The example that we consider is the liquidation of a large position executed over the course of a single trading day in a limit order book. Transaction costs are incurred because large orders walk the order book, that is, they consume order book liquidity beyond the best bid/ask. We model the order book with a power law that is proportional to trading volume, and thus transaction costs are inversely proportional to a power of trading volume. We obtain a policy approximation by training a long short term memory (LSTM) neural network to minimize transaction costs accumulated when execution is carried out as a sequence of smaller suborders. Using historical S&P100 price and volume data, we evaluate our LSTM strategy relative to strategies based on time-weighted average price (TWAP) and volume-weighted average price (VWAP). For execution of a single stock, the input to the LSTM is the cross section of data on all 100 stocks, including prices, volumes, TWAPs and VWAPs. By using this data cross section, the LSTM should be able to exploit inter-stock co-dependence in volume and price movements, thereby reducing transaction costs for the day. Our tests on S&P100 data demonstrate that in fact this is so, as our LSTM strategy consistently outperforms TWAP and VWAP-based strategies.
△ Less
Submitted 14 June, 2023; v1 submitted 23 January, 2023;
originally announced January 2023.
-
Bayesian model selection for GRB 211211A through multi-wavelength analyses
Authors:
Nina Kunert,
Sarah Antier,
Vsevolod Nedora,
Mattia Bulla,
Peter T. H. Pang,
Shreya Anand,
Michael Coughlin,
Ingo Tews,
Jennifer Barnes,
Thomas Hussenot-Desenonges,
Brian Healy,
Theophile Jegou du Laz,
Meili Pilloix,
Weizmann Kiendrebeogo,
Tim Dietrich
Abstract:
Although GRB 211211A is one of the closest gamma-ray bursts (GRBs), its classification is challenging because of its partially inconclusive electromagnetic signatures. In this paper, we investigate four different astrophysical scenarios as possible progenitors for GRB~211211A: a binary neutron-star merger, a black-hole--neutron-star merger, a core-collapse supernova, and an r-process enriched core…
▽ More
Although GRB 211211A is one of the closest gamma-ray bursts (GRBs), its classification is challenging because of its partially inconclusive electromagnetic signatures. In this paper, we investigate four different astrophysical scenarios as possible progenitors for GRB~211211A: a binary neutron-star merger, a black-hole--neutron-star merger, a core-collapse supernova, and an r-process enriched core collapse of a rapidly rotating massive star (a collapsar). We perform a large set of Bayesian multi-wavelength analyses based on different models describing these scenarios and priors to investigate which astrophysical scenarios and processes might be related to GRB~211211A. Our analysis supports previous studies in which the presence of an additional component, likely related to $r$-process nucleosynthesis, is required to explain the observed light curves of GRB~211211A, as it can not solely be explained as a GRB afterglow. Fixing the distance to about $350~\rm Mpc$, namely the distance of the possible host galaxy SDSS J140910.47+275320.8, we find a statistical preference for a binary neutron-star merger scenario.
△ Less
Submitted 15 December, 2023; v1 submitted 5 January, 2023;
originally announced January 2023.
-
Evidence that the Hot Jupiter WASP-77 A b Formed Beyond Its Parent Protoplanetary Disk's H2O Ice Line
Authors:
Henrique Reggiani,
Kevin C. Schlaufman,
Brian F. Healy,
Joshua D. Lothringer,
David K. Sing
Abstract:
Idealized protoplanetary disk and giant planet formation models have been interpreted to suggest that a giant planet's atmospheric abundances can be used to infer its formation location in its parent protoplanetary disk. It has recently been reported that the hot Jupiter WASP-77 A b has sub-solar atmospheric carbon and oxygen abundances with a solar C/O abundance ratio. Assuming solar carbon and o…
▽ More
Idealized protoplanetary disk and giant planet formation models have been interpreted to suggest that a giant planet's atmospheric abundances can be used to infer its formation location in its parent protoplanetary disk. It has recently been reported that the hot Jupiter WASP-77 A b has sub-solar atmospheric carbon and oxygen abundances with a solar C/O abundance ratio. Assuming solar carbon and oxygen abundances for its host star WASP-77 A, WASP-77 A b's atmospheric carbon and oxygen abundances possibly indicate that it accreted its envelope interior to its parent protoplanetary disk's H2O ice line from carbon-depleted gas with little subsequent planetesimal accretion or core erosion. We comprehensively model WASP-77 A and use our results to better characterize WASP-77 A b. We show that the photospheric abundances of carbon and oxygen in WASP-77 A are super-solar with a sub-solar C/O abundance ratio, implying that WASP-77 A b's atmosphere has significantly sub-stellar carbon and oxygen abundances with a super-stellar C/O ratio. Our result possibly indicates that WASP-77 A b's envelope was accreted by the planet beyond its parent protoplanetary disk's H2O ice line. While numerous theoretical complications to these idealized models have now been identified, the possibility of non-solar protoplanetary disk abundance ratios confound even the most sophisticated protoplanetary disk and giant planet formation models. We therefore argue that giant planet atmospheric abundance ratios can only be meaningfully interpreted relative to the possibly non-solar mean compositions of their parent protoplanetary disks as recorded in the photospheric abundances of their solar-type dwarf host stars.
△ Less
Submitted 1 February, 2022; v1 submitted 20 January, 2022;
originally announced January 2022.
-
Group boundaries for semidirect products with Z
Authors:
Craig R. Guilbault,
Brendan Burns Healy,
Brian Pietsch
Abstract:
Bestvina's notion of a Z-structure provides a general framework for group boundaries that includes Gromov boundaries of hyperbolic groups and visual boundaries of CAT(0) groups as special cases. A refinement, known as an EZ-structure has proven useful in attacks on the Novikov Conjecture and related problems. Characterizations of groups admitting a Z- or EZ-structure are longstanding open problems…
▽ More
Bestvina's notion of a Z-structure provides a general framework for group boundaries that includes Gromov boundaries of hyperbolic groups and visual boundaries of CAT(0) groups as special cases. A refinement, known as an EZ-structure has proven useful in attacks on the Novikov Conjecture and related problems. Characterizations of groups admitting a Z- or EZ-structure are longstanding open problems. In this paper, we examine semidirect products of a group G with the integers. For example, we show that, if G is torsion-free and admits a Z-structure, then so does every semidirect product of this type. We prove a similar theorem for EZ-structures, under an additional hypothesis.
As applications, we show that all closed 3-manifold groups admit Z-structures, as do all strongly polycyclic groups, and all groups of polynomial growth. In those latter two cases our Z-boundaries are always spheres. This allows one to make strong conclusions about the group cohomology and end invariants of those groups. In another direction, we expand upon the notion of an EZ-structure and discuss new applications to the Novikov Conjecture.
△ Less
Submitted 18 July, 2022; v1 submitted 24 September, 2021;
originally announced September 2021.
-
Stellar Spins in the Pleiades, Praesepe, and M35 Open Clusters
Authors:
Brian F. Healy,
P. R. McCullough,
Kevin C. Schlaufman
Abstract:
We analyze spectroscopic and photometric data to determine the projected inclinations of stars in three open clusters: the Pleiades, Praesepe, and M35. We determine the $\sin i$ values of 42, 35, and 67 stars in each cluster, respectively, and from their distributions we find that isotropic spins and moderate alignment are both consistent with the Pleiades and Praesepe data. While it is difficult…
▽ More
We analyze spectroscopic and photometric data to determine the projected inclinations of stars in three open clusters: the Pleiades, Praesepe, and M35. We determine the $\sin i$ values of 42, 35, and 67 stars in each cluster, respectively, and from their distributions we find that isotropic spins and moderate alignment are both consistent with the Pleiades and Praesepe data. While it is difficult to distinguish between these scenarios for a single cluster, an ensemble of such distributions may facilitate a distinction. The M35 inclination distribution is most consistent with a superposition of isotropic and anisotropic spins, the source of which could be systematic error or a physical grouping of aligned stars. We also study internal cluster kinematics using radial velocities and proper motions. Our kinematics analysis reveals significant plane-of-sky rotation in Praesepe, with a mean velocity of $0.132 \pm 0.022$ km s$^{-1}$ in a clockwise direction.
△ Less
Submitted 5 November, 2021; v1 submitted 17 September, 2021;
originally announced September 2021.
-
Wallpaper Groups and Auxetic Metamaterials
Authors:
Brendan Burns Healy,
Aparna Deshmukh,
Elliott Fairchild,
Caroline J. Merighi,
Konstantin Sobolev
Abstract:
We examine a fundamental material property called Poisson's ratio, which establishes the relationship for the relative deformation of a physical system in orthogonal directions. Architects and engineers have designed advanced systems using repeating patterns that can potentially exhibit auxetic behavior, which is the property of having a negative Poisson's ratio. Because two-dimensional cross sect…
▽ More
We examine a fundamental material property called Poisson's ratio, which establishes the relationship for the relative deformation of a physical system in orthogonal directions. Architects and engineers have designed advanced systems using repeating patterns that can potentially exhibit auxetic behavior, which is the property of having a negative Poisson's ratio. Because two-dimensional cross sections of each of these patterns has an associated wallpaper group, we can look for useful correlations between this geometric information and the two-dimensional response as defined by Poisson's ratio. By analyzing the data, we find two properties of the wallpaper group that correlate with more effective Poisson's ratio required for applications. This paper also contains an introduction to wallpaper groups and orbifold notation and an appendix contains some literature references recreated to use our preferred notation.
△ Less
Submitted 25 February, 2021;
originally announced February 2021.
-
Cusped spaces and quasi-isometries of relatively hyperbolic groups
Authors:
Brendan Burns Healy,
G. Christopher Hruska
Abstract:
A group $Γ$ with a family of subgroups $\mathbb{P}$ is relatively hyperbolic if $Γ$ admits a cusp-uniform action on a proper $δ$--hyperbolic space. We show that any two such spaces for a given group pair are quasi-isometric, provided the spaces have "constant horospherical distortion," a condition satisfied by Groves--Manning's cusped Cayley graph and by all negatively curved symmetric spaces. Con…
▽ More
A group $Γ$ with a family of subgroups $\mathbb{P}$ is relatively hyperbolic if $Γ$ admits a cusp-uniform action on a proper $δ$--hyperbolic space. We show that any two such spaces for a given group pair are quasi-isometric, provided the spaces have "constant horospherical distortion," a condition satisfied by Groves--Manning's cusped Cayley graph and by all negatively curved symmetric spaces. Consequently the Bowditch boundary admits a canonical quasisymmetric structure, which coincides with the "naturally occurring" quasisymmetric structure of the symmetric space when considering lattices in rank one symmetric spaces.
We show that a group $Γ$ is a lattice in a negatively curved symmetric space $X$ if and only if a cusped space for $Γ$ is quasi-isometric to the symmetric space. We also prove an ideal triangle characterization of the $δ$--hyperbolic spaces with uniformly perfect boundary due to Meyer and Bourdon--Kleiner. An appendix concerns the equivalence of several definitions of conical limit point found in the literature.
△ Less
Submitted 8 March, 2021; v1 submitted 19 October, 2020;
originally announced October 2020.
-
Stellar Spins in the Open Cluster NGC 2516
Authors:
Brian F. Healy,
P. R. McCullough
Abstract:
Measuring the distribution of stellar spin axis orientations in a coeval group of stars probes the physical processes underlying the stars' formation. In this paper, we use spectro-photometric observations of the open cluster NGC 2516 to determine the degree of spin alignment among its stars. We combine TESS light curves, ground-based spectroscopy from the Gaia-ESO and GALAH surveys, broad-band st…
▽ More
Measuring the distribution of stellar spin axis orientations in a coeval group of stars probes the physical processes underlying the stars' formation. In this paper, we use spectro-photometric observations of the open cluster NGC 2516 to determine the degree of spin alignment among its stars. We combine TESS light curves, ground-based spectroscopy from the Gaia-ESO and GALAH surveys, broad-band stellar magnitudes from several surveys, and Gaia astrometry to measure 33 stellar inclinations and quantify overall cluster rotation. Our measurements suggest that stellar spins in this cluster are isotropically oriented, while allowing for the possibility that they are moderately aligned. An isotropic distribution of NGC 2516 spins would imply a star-forming environment in which turbulence dominated ordered motion, while a moderately aligned distribution would suggest a more substantial contribution from rotation. We also perform a three-dimensional analysis of the cluster's internal kinematics, finding no significant signatures of overall rotation. Stemming from this analysis, we identify evidence of cluster contraction, suggesting possible ongoing mass segregation in NGC 2516.
△ Less
Submitted 6 November, 2020; v1 submitted 29 September, 2020;
originally announced September 2020.
-
Pinched Curvature in Heintze Groups of Carnot-type
Authors:
Brendan Burns Healy
Abstract:
Rank-one symmetric spaces carry a solvable group model which have a generalization to a larger class of Lie groups that are one-dimensional extensions of nilpotent groups. By examining some metric properties of these symmetric spaces, we motivate and prove the existence of analogous left-invariant, Riemannian metrics on Heintze groups of Carnot-type. These metrics adhere to certain natural curvatu…
▽ More
Rank-one symmetric spaces carry a solvable group model which have a generalization to a larger class of Lie groups that are one-dimensional extensions of nilpotent groups. By examining some metric properties of these symmetric spaces, we motivate and prove the existence of analogous left-invariant, Riemannian metrics on Heintze groups of Carnot-type. These metrics adhere to certain natural curvature pinching properties, and we show in a special case that this pinching is optimal, appealing to a result of Belegradek and Kapovitch.
△ Less
Submitted 23 February, 2021; v1 submitted 11 February, 2020;
originally announced February 2020.
-
PCA for Implied Volatility Surfaces
Authors:
Marco Avellaneda,
Brian Healy,
Andrew Papanicolaou,
George Papanicolaou
Abstract:
Principal component analysis (PCA) is a useful tool when trying to construct factor models from historical asset returns. For the implied volatilities of U.S. equities there is a PCA-based model with a principal eigenportfolio whose return time series lies close to that of an overarching market factor. The authors show that this market factor is the index resulting from the daily compounding of a…
▽ More
Principal component analysis (PCA) is a useful tool when trying to construct factor models from historical asset returns. For the implied volatilities of U.S. equities there is a PCA-based model with a principal eigenportfolio whose return time series lies close to that of an overarching market factor. The authors show that this market factor is the index resulting from the daily compounding of a weighted average of implied-volatility returns, with weights based on the options' open interest (OI) and Vega. The authors also analyze the singular vectors derived from the tensor structure of the implied volatilities of S&P500 constituents, and find evidence indicating that some type of OI and Vega-weighted index should be one of at least two significant factors in this market.
△ Less
Submitted 31 January, 2020;
originally announced February 2020.
-
XO-7 b: A transiting hot Jupiter with a massive companion on a wide orbit
Authors:
Nicolas Crouzet,
Brian F. Healy,
Guillaume Hébrard,
P. R. McCullough,
Doug Long,
Pilar Montañés-Rodríguez,
Ignasi Ribas,
Francesc Vilardell,
Enrique Herrero,
Enrique Garcia-Melendo,
Matthieu Conjat,
Jerry Foote,
Joe Garlitz,
Phillip Vo,
Nuno C. Santos,
Jos de Bruijne,
Hugh P. Osborn,
Shweta Dalal,
Louise D. Nielsen
Abstract:
Transiting planets orbiting bright stars are the most favorable targets for follow-up and characterization. We report the discovery of the transiting hot Jupiter XO-7 b and of a second, massive companion on a wide orbit around a circumpolar, bright, and metal rich G0 dwarf (V = 10.52, $T_{\rm eff} = 6250 \pm 100 \; \rm K$, $\rm[Fe/H] = 0.432 \pm 0.057 \; \rm dex$). We conducted photometric and rad…
▽ More
Transiting planets orbiting bright stars are the most favorable targets for follow-up and characterization. We report the discovery of the transiting hot Jupiter XO-7 b and of a second, massive companion on a wide orbit around a circumpolar, bright, and metal rich G0 dwarf (V = 10.52, $T_{\rm eff} = 6250 \pm 100 \; \rm K$, $\rm[Fe/H] = 0.432 \pm 0.057 \; \rm dex$). We conducted photometric and radial velocity follow-up with a team of amateur and professional astronomers. XO-7 b has a period of $ 2.8641424 \pm 0.0000043$ days, a mass of $0.709 \pm 0.034 \; \rm M_{\rm J}$, a radius of $1.373 \pm 0.026 \; \rm R_{\rm J}$, a density of $0.340 \pm 0.027 \; \rm g \, {cm}^{-3}$, and an equilibrium temperature of $1743 \pm 23 \; \rm K$. Its large atmospheric scale height and the brightness of the host star make it well suited to atmospheric characterization. The wide orbit companion is detected as a linear trend in radial velocities with an amplitude of $\sim100 \; \rm m \, {s}^{-1}$ over two years, yielding a minimum mass of $4 \; \rm M_{\rm J}$; it could be a planet, a brown dwarf, or a low mass star. The hot Jupiter orbital parameters and the presence of the wide orbit companion point towards a high eccentricity migration for the hot Jupiter. Overall, this system will be valuable to understand the atmospheric properties and migration mechanisms of hot Jupiters and will help constrain the formation and evolution models of gas giant exoplanets.
△ Less
Submitted 5 December, 2019;
originally announced December 2019.
-
Touché: Towards Ideal and Efficient Cache Compression By Mitigating Tag Area Overheads
Authors:
Seokin Hong,
Bulent Abali,
Alper Buyuktosunoglu,
Michael B. Healy,
Prashant J. Nair
Abstract:
Compression is seen as a simple technique to increase the effective cache capacity. Unfortunately, compression techniques either incur tag area overheads or restrict data placement to only include neighboring compressed cache blocks to mitigate tag area overheads. Ideally, we should be able to place arbitrary compressed cache blocks without any placement restrictions and tag area overheads.
This…
▽ More
Compression is seen as a simple technique to increase the effective cache capacity. Unfortunately, compression techniques either incur tag area overheads or restrict data placement to only include neighboring compressed cache blocks to mitigate tag area overheads. Ideally, we should be able to place arbitrary compressed cache blocks without any placement restrictions and tag area overheads.
This paper proposes Touché, a framework that enables storing multiple arbitrary compressed cache blocks within a physical cacheline without any tag area overheads. The Touché framework consists of three components. The first component, called the ``Signature'' (SIGN) engine, creates shortened signatures from the tag addresses of compressed blocks. Due to this, the SIGN engine can store multiple signatures in each tag entry. On a cache access, the physical cacheline is accessed only if there is a signature match (which has a negligible probability of false positive). The second component, called the ``Tag Appended Data'' (TADA) mechanism, stores the full tag addresses with data. TADA enables Touché to detect false positive signature matches by ensuring that the actual tag address is available for comparison. The third component, called the ``Superblock Marker'' (SMARK) mechanism, uses a unique marker in the tag entry to indicate the occurrence of compressed cache blocks from neighboring physical addresses in the same cacheline. Touché is completely hardware-based and achieves an average speedup of 12\% (ideal 13\%) when compared to an uncompressed baseline.
△ Less
Submitted 2 September, 2019;
originally announced September 2019.
-
Magnetic Inflation and Stellar Mass III: Revised Parameters for the Component Stars of NSVS 07394765
Authors:
Brian F. Healy,
Eunkyu Han,
Philip S. Muirhead,
Brian Skiff,
Tom Polakis,
Anneliese M. Rilinger,
Jonathan J. Swift
Abstract:
We perform a new analysis of the M dwarf-M dwarf eclipsing binary system NSVS 07394765 in order to investigate the reported hyper-inflated radius of one of the component stars. Our analysis is based on archival photometry from the Wide Angle Search for Planets (WASP), new photometry from the 32 cm {Command Module Observatory (CMO) telescope in Arizona and the 70 cm telescope at Thacher Observatory…
▽ More
We perform a new analysis of the M dwarf-M dwarf eclipsing binary system NSVS 07394765 in order to investigate the reported hyper-inflated radius of one of the component stars. Our analysis is based on archival photometry from the Wide Angle Search for Planets (WASP), new photometry from the 32 cm {Command Module Observatory (CMO) telescope in Arizona and the 70 cm telescope at Thacher Observatory in California}, and new high-resolution infrared spectra obtained with the Immersion Grating Infrared Spectrograph (IGRINS) on the Discovery Channel Telescope. The masses and radii we measure for each component star disagree with previously reported measurements. We show that both stars are early M-type main-sequence stars without evidence for youth or hyper-inflation ($M_1= 0.661\ ^{+0.008}_{-0.036}\ \rm{M_{sun}}$, $M_2= 0.608\ ^{+0.003}_{-0.028}\ \rm{M_{sun}}$, $R_1= 0.599\ ^{+0.032}_{-0.019}\ \rm{R_{sun}}$, $R_2= 0.625\ ^{+0.012}_{-0.027}\ \rm{R_{sun}}$), and we update the orbital period and eclipse ephemerides for the system. We suggest that the likely cause of the initial hyper-inflated result is the use of moderate-resolution spectroscopy for precise radial velocity measurements.
△ Less
Submitted 8 July, 2019;
originally announced July 2019.
-
Data science is science's second chance to get causal inference right: A classification of data science tasks
Authors:
Miguel A. Hernán,
John Hsu,
Brian Healy
Abstract:
Causal inference from observational data is the goal of many data analyses in the health and social sciences. However, academic statistics has often frowned upon data analyses with a causal objective. The introduction of the term "data science" provides a historic opportunity to redefine data analysis in such a way that it naturally accommodates causal inference from observational data. Like other…
▽ More
Causal inference from observational data is the goal of many data analyses in the health and social sciences. However, academic statistics has often frowned upon data analyses with a causal objective. The introduction of the term "data science" provides a historic opportunity to redefine data analysis in such a way that it naturally accommodates causal inference from observational data. Like others before, we organize the scientific contributions of data science into three classes of tasks: Description, prediction, and counterfactual prediction (which includes causal inference). An explicit classification of data science tasks is necessary to discuss the data, assumptions, and analytics required to successfully accomplish each task. We argue that a failure to adequately describe the role of subject-matter expert knowledge in data analysis is a source of widespread misunderstandings about data science. Specifically, causal analyses typically require not only good data and algorithms, but also domain expert knowledge. We discuss the implications for the use of data science to guide decision-making in the real world and to train data scientists.
△ Less
Submitted 6 April, 2019; v1 submitted 28 April, 2018;
originally announced April 2018.
-
Rigidity Properties for Hyperbolic Generalizations
Authors:
Brendan Burns Healy
Abstract:
We make a few observations on the absence of geometric and topological rigidity for acylindrically hyperbolic and relatively hyperbolic groups. In particular, we demonstrate the lack of a well-defined limit set for acylindrical actions on hyperbolic spaces, even under the assumption of universality. We also prove a statement about relatively hyperbolic groups inspired by a remark by Groves, Mannin…
▽ More
We make a few observations on the absence of geometric and topological rigidity for acylindrically hyperbolic and relatively hyperbolic groups. In particular, we demonstrate the lack of a well-defined limit set for acylindrical actions on hyperbolic spaces, even under the assumption of universality. We also prove a statement about relatively hyperbolic groups inspired by a remark by Groves, Manning, and Sisto about the quasi-isometry type of combinatorial cusps. Finally, we summarize these results in a table in order to assert a meta-statement about the decay of metric rigidity as the conditions on actions on hyperbolic spaces are loosened.
△ Less
Submitted 15 July, 2019; v1 submitted 27 March, 2018;
originally announced March 2018.
-
Acylindrical Hyperbolicity of Out($W_n$)
Authors:
Brendan Burns Healy
Abstract:
We prove that the group of outer automorphisms of the free Coxeter group $W_n$ is acylindrically hyperbolic in the sense of Osin. As an application, we observe that any CAT(0) space admitting a geometric action by Out($W_n$) must contain a rank-one geodesic. The theorem proceeds from expanding on a well-known relationship between Out($W_n$) and the outer automorphism group of free groups.
We prove that the group of outer automorphisms of the free Coxeter group $W_n$ is acylindrically hyperbolic in the sense of Osin. As an application, we observe that any CAT(0) space admitting a geometric action by Out($W_n$) must contain a rank-one geodesic. The theorem proceeds from expanding on a well-known relationship between Out($W_n$) and the outer automorphism group of free groups.
△ Less
Submitted 22 October, 2019; v1 submitted 25 October, 2016;
originally announced October 2016.
-
Quantum Monte Carlo calculations of H$_2$ dissociation on Si(001)
Authors:
Claudia Filippi,
Sorcha B. Healy,
P. Kratzer,
E. Pehlke,
M. Scheffler
Abstract:
We present quantum Monte Carlo calculations for various reaction pathways of H$_2$ with Si(001), using large model clusters of the surface. We obtain reaction energies and energy barriers noticeably higher than those from approximate exchange-correlation functionals. In improvement over previous studies, our adsorption barriers closely agree with experimental data. For desorption, the calculatio…
▽ More
We present quantum Monte Carlo calculations for various reaction pathways of H$_2$ with Si(001), using large model clusters of the surface. We obtain reaction energies and energy barriers noticeably higher than those from approximate exchange-correlation functionals. In improvement over previous studies, our adsorption barriers closely agree with experimental data. For desorption, the calculations give barriers for conventional pathways in excess of the presently accepted experimental value, and pinpoint the role of coverage effects and desorption from steps.
△ Less
Submitted 6 May, 2002;
originally announced May 2002.
-
The role of electronic correlation in the Si(100) reconstruction: a quantum Monte Carlo study
Authors:
Sorcha B. Healy,
Claudia Filippi,
P. Kratzer,
E. Penev,
M. Scheffler
Abstract:
Recent low-temperature scanning tunneling experiments have challenged the generally accepted picture of buckled silicon dimers as the ground state reconstruction of the Si(100) surface. Together with the symmetric dimer model of the surface suggested by quantum chemistry calculations on small clusters, these findings question our general understanding of electronic correlations at surfaces and i…
▽ More
Recent low-temperature scanning tunneling experiments have challenged the generally accepted picture of buckled silicon dimers as the ground state reconstruction of the Si(100) surface. Together with the symmetric dimer model of the surface suggested by quantum chemistry calculations on small clusters, these findings question our general understanding of electronic correlations at surfaces and its proper description within density functional theory. We present quantum Monte Carlo calculations on large cluster models of the symmetric and buckled surface, and conclude that buckling remains energetically more favorable even when the present-day best treatment of electronic correlation is employed.
△ Less
Submitted 26 February, 2001;
originally announced February 2001.
-
A non-linear optimal estimation inverse method for radio occultation measurements of temperature, humidity and surface pressure
Authors:
Paul I. Palmer,
J. J. Barnett,
J. R. Eyre,
S. B. Healy
Abstract:
An optimal estimation inverse method is presented which can be used to retrieve simultaneously vertical profiles of temperature and specific humidity, in addition to surface pressure, from satellite-to-satellite radio occultation observations of the Earth's atmosphere. The method is a non-linear, maximum {\it a posteriori} technique which can accommodate most aspects of the real radio occultatio…
▽ More
An optimal estimation inverse method is presented which can be used to retrieve simultaneously vertical profiles of temperature and specific humidity, in addition to surface pressure, from satellite-to-satellite radio occultation observations of the Earth's atmosphere. The method is a non-linear, maximum {\it a posteriori} technique which can accommodate most aspects of the real radio occultation problem and is found to be stable and to converge rapidly in most cases. The optimal estimation inverse method has two distinct advantages over the analytic inverse method in that it accounts for some of the effects of horizontal gradients and is able to retrieve optimally temperature and humidity simultaneously from the observations. It is also able to account for observation noise and other sources of error. Combined, these advantages ensure a realistic retrieval of atmospheric quantities.
A complete error analysis emerges naturally from the optimal estimation theory, allowing a full characterisation of the solution. Using this analysis a quality control scheme is implemented which allows anomalous retrieval conditions to be recognised and removed, thus preventing gross retrieval errors.
The inverse method presented in this paper has been implemented for bending angle measurements derived from GPS/MET radio occultation observations of the Earth. Preliminary results from simulated data suggest that these observations have the potential to improve NWP model analyses significantly throughout their vertical range.
△ Less
Submitted 3 March, 2000;
originally announced March 2000.