skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Search for: All records

Creators/Authors contains: "Duffell, Paul C"

Note: When clicking on a Digital Object Identifier (DOI) number, you will be taken to an external site maintained by the publisher. Some full text articles may not yet be available without a charge during the embargo (administrative interval).
What is a DOI Number?

Some links on this page may take you to non-federal websites. Their policies may differ from this site.

  1. Abstract Many studies have recently documented the orbital response of eccentric binaries accreting from thin circumbinary disks, characterizing the change in the binary semimajor axis and eccentricity. We extend these calculations to include the precession of the binary’s longitude of periapse induced by the circumbinary disk, and we characterize this precession continuously with binary eccentricityebfor equal mass components. This disk-induced apsidal precession is prograde with a weak dependence on the binary eccentricity wheneb≲ 0.4 and decreases approximately linearly foreb≳ 0.4; yet at allebbinary precession is faster than the rates of change to the semimajor axis and eccentricity by an order of magnitude. We estimate that such precession effects are likely most important for subparsec separated binaries with masses ≲107M, like LISA precursors. We find that accreting, equal-mass LISA binaries withM< 106M(and the most massiveM∼ 107Mbinaries out toz∼ 3) may acquire a detectable phase offset due to the disk-induced precession. Moreover, disk-induced precession can compete with general relativistic precession in a vacuum, making it important for observer-dependent electromagnetic searches for accreting massive binaries—like Doppler boost and binary self-lensing models—after potentially only a few orbital periods. 
    more » « less
  2. Abstract We develop a suite of 3D hydrodynamic models of supernova remnants (SNRs) expanding against the circumstellar medium (CSM). We study the Rayleigh–Taylor instability forming at the expansion interface by calculating an angular power spectrum for each of these models. The power spectra of young SNRs are seen to exhibit a dominant angular mode, which is a diagnostic of their ejecta density profile as found by previous studies. The steep scaling of power at smaller modes and the time evolution of the spectra are indicative of the absence of a turbulent cascade. Instead, as the time evolution of the spectra suggests, they may be governed by an angular mode-dependent net growth rate. We also study the impact of anisotropies in the ejecta and in the CSM on the power spectra of velocity and density. We confirm that perturbations in the density field (whether imposed on the ejecta or the CSM) do not influence the anisotropy of the remnant significantly unless they have a very large amplitude and form large-scale coherent structures. In any case, these clumps can only affect structures on large angular scales. The power spectrum on small angular scales is completely independent of the initial clumpiness and governed only by the growth and saturation of the Rayleigh–Taylor instability. 
    more » « less
  3. Abstract The astrophysical environments capable of triggering heavy-element synthesis via rapid neutron capture (ther-process) remain uncertain. While binary neutron star mergers (NSMs) are known to forger-process elements, certain rare supernovae (SNe) have been theorized to supplement—or even dominate—r-production by NSMs. However, the most direct evidence for such SNe, unusual reddening of the emission caused by the high opacities ofr-process elements, has not been observed. Recent work identified the distribution ofr-process material within the SN ejecta as a key predictor of the ease with which signals associated withr-process enrichment could be discerned. Though this distribution results from hydrodynamic processes at play during the SN explosion, thus far it has been treated only in a parameterized way. We use hydrodynamic simulations to model how disk winds—the alleged locus ofr-production in rare SNe—mix with initiallyr-process-free ejecta. We study mixing as a function of the wind mass, wind duration, and the initial SN explosion energy, and find that it increases with the first two of these and decreases with the third. This suggests that SNe accompanying the longest long-duration gamma-ray bursts are promising places to search for signs ofr-process enrichment. We use semianalytic radiation transport to connect hydrodynamics to electromagnetic observables, allowing us to assess the mixing level at which the presence ofr-process material can be diagnosed from SN light curves. Analytic arguments constructed atop this foundation imply that a wind-drivenr-process-enriched SN model is unlikely to explain standard energetic SNe. 
    more » « less
  4. Abstract Active galactic nuclei (AGN) show a range of morphologies and dynamical properties, which are determined not only by parameters intrinsic to the central engine but also their interaction with the surrounding environment. We investigate the connection of kiloparsec scale AGN jet properties to their intrinsic parameters and surroundings. This is done using a suite of 40 relativistic hydrodynamic simulations spanning a wide range of engine luminosities and opening angles. We explore AGN jet propagation with different ambient density profiles, including r −2 (self-similar solution) and r −1 , which is more relevant for AGN host environments. While confirmation awaits future 3D studies, the Fanaroff–Riley (FR) morphological dichotomy arises naturally in our 2D models. Jets with low energy density compared to the ambient medium produce a center-brightened emissivity distribution, while emissivity from relatively higher energy density jets is dominated by the jet head. We observe recollimation shocks in our simulations that can generate bright spots along the spine of the jet, providing a possible explanation for “knots” observed in AGN jets. We additionally find a scaling relation between the number of knots and the jet-head-to-surroundings energy density ratio. This scaling relation is generally consistent with the observations of the jets in M87 and Cygnus A. Our model also correctly predicts M87 as FRI and Cygnus A as FRII. Our model can be used to relate jet dynamical parameters such as jet head velocity, jet opening angle, and external pressure to jet power, and ambient density estimates. 
    more » « less
  5. ABSTRACT A subset of Type Ic supernovae (SNe Ic), broad-lined SNe Ic (SNe Ic-bl), show unusually high kinetic energies (∼1052 erg) that cannot be explained by the energy supplied by neutrinos alone. Many SNe Ic-bl have been observed in coincidence with long gamma-ray bursts (GRBs) that suggests a connection between SNe and GRBs. A small fraction of core-collapse supernovae form a rapidly rotating and strongly magnetized protoneutron star (PNS), a proto-magnetar. Jets from such magnetars can provide the high kinetic energies observed in SNe Ic-bl and also provide the connection to GRBs. In this work, we use the jetted outflow produced in a 3D general relativistic magnetohydrodynamic CCSN simulation from a consistently formed proto-magnetar as the central engine for full-star explosion simulations. We extract a range of central engine parameters and find that the extracted engine energy is in the range of 6.231 × 1051−1.725 × 1052 erg, the engine time-scale in the range of 0.479−1.159 s and the engine half-opening angle in the range of ∼9°−19°. Using these as central engines, we perform 2D special relativistic (SR) hydrodynamic (HD) and radiation transfer simulations to calculate the corresponding light curves and spectra. We find that these central engine parameters successfully produce SNe Ic-bl that demonstrates that jets from proto-magnetars can be viable engines for SNe Ic-bl. We also find that only the central engines with smaller opening angles (∼11°) form a GRB implying that GRB formation is likely associated with narrower jet outflows and Ic-bl’s without GRBs may be associated with wider outflows. 
    more » « less
  6. Abstract We report multiwavelength observations and characterization of the ultraluminous transient AT 2021lwx (ZTF20abrbeie; aka “Barbie”) identified in the alert stream of the Zwicky Transient Facility (ZTF) using a Recommender Engine For Intelligent Transient Tracking filter on the ANTARES alert broker. From a spectroscopically measured redshift of 0.995, we estimate a peak-observed pseudo-bolometric luminosity of log( L max / [ erg s − 1 ] ) = 45.7 from slowly fading ztf- g and ztf- r light curves spanning over 1000 observer-frame days. The host galaxy is not detected in archival Pan-STARRS observations ( g > 23.3 mag), implying a lower limit to the outburst amplitude of more than 5 mag relative to the quiescent host galaxy. Optical spectra exhibit strong emission lines with narrow cores from the H Balmer series and ultraviolet semi-forbidden lines of Si iii ] λ 1892, C iii ] λ 1909, and  C ii ] λ 2325. Typical nebular lines in Active Galactic Nucleus (AGN) spectra from ions such as [O ii ] and [O iii ] are not detected. These spectral features, along with the smooth light curve that is unlike most AGN flaring activity and the luminosity that exceeds any observed or theorized supernova, lead us to conclude that AT 2021lwx is most likely an extreme tidal disruption event (TDE). Modeling of ZTF photometry with MOSFiT suggests that the TDE was between a ≈14 M ⊙ star and a supermassive black hole of mass M BH ∼ 10 8 M ⊙ . Continued monitoring of the still-evolving light curve along with deep imaging of the field once AT 2021lwx has faded can test this hypothesis and potentially detect the host galaxy. 
    more » « less
  7. null (Ed.)
    ABSTRACT The merger of two neutron stars produces an outflow of radioactive heavy nuclei. Within a second of merger, the central remnant is expected to also launch a relativistic jet, which shock-heats and disrupts a portion of the radioactive ejecta. Within a few hours, emission from the radioactive material gives rise to an ultraviolet, optical, and infrared transient (a kilonova). We use the endstates of a suite of 2D relativistic hydrodynamic simulations of jet–ejecta interaction as initial conditions for multidimensional Monte Carlo radiation transport simulations of the resulting viewing angle-dependent light curves and spectra starting at $$1.5\, \mathrm{h}$$ after merger. We find that on this time-scale, jet shock heating does not affect the kilonova emission for the jet parameters we survey. However, the jet disruption to the density structure of the ejecta does change the light curves. The jet carves a channel into the otherwise spheroidal ejecta, revealing the hot, inner regions. As seen from near (≲30°) the jet axis, the kilonova is brighter by a factor of a few and bluer. The strength of this effect depends on the jet parameters, since the light curves of more heavily disrupted ejecta are more strongly affected. The light curves and spectra are also more heavily modified in the ultraviolet than in the optical. 
    more » « less
  8. null (Ed.)