Search

The Advanced X-ray Imaging Satellite (AXIS) promises revolutionary science in the X-ray and multi-messenger time domain. AXIS will leverage excellent spatial resolution (80× that of Swift), and a large collecting area (5–10× that of Chandra) across a 24-arcmin diameter field of view at soft X-ray energies (0.3–10.0 keV) to discover and characterize a wide range of X-ray transients from supernova-shock breakouts to tidal disruption events to highly variable supermassive black holes. The observatory’s ability to localize and monitor faint X-ray sources opens up new opportunities to hunt for counterparts to distant binary neutron star mergers, fast radio bursts, and exotic phenomena like fast X-ray transients. AXIS will offer a response time of

We introduce a new capability of the Neil Gehrels Swift Observatory, dubbed “continuous commanding,” that achieves 10 s latency response time on orbit to unscheduled target-of-opportunity requests received on the ground. We show that this will allow Swift to respond to premerger (early-warning) gravitational-wave (GW) detections, rapidly slewing the Burst Alert Telescope (BAT) across the sky to place the GW origin in the BAT field of view at or before merger time. This will dramatically increase the GW/gamma-ray burst (GRB) codetection rate and enable prompt arcminute localization of a neutron star merger. We simulate the full Swift response to a GW early-warning alert, including input sky maps produced at different early-warning times, a complete model of the Swift attitude control system, and a full accounting of the latency between the GW detectors and the spacecraft. 60 s of early warning can double the rate of a prompt GRB detection with arcminute localization, and 140 s guarantees observation anywhere on the unocculted sky, even with localization areas ≫1000 deg ^2 . While 140 s is beyond current GW detector sensitivities, 30–70 s is achievable today. We show that the detection yield is now limited by the latency of LIGO/Virgo cyberinfrastructure and motivate a focus on its reduction. Continuous commanding has been integrated as a general capability of Swift, significantly increasing its versatility in response to the growing demands of time-domain astrophysics. We demonstrate this potential on an externally triggered fast radio burst (FRB), slewing 81° across the sky, and collecting X-ray and UV photons from the source position

We present SN 2023zaw—a subluminous ( M _r = −16.7 mag) and rapidly evolving supernova ( t _1/2, _r = 4.9 days), with the lowest nickel mass (≈0.002 M _⊙ ) measured among all stripped-envelope supernovae discovered to date. The photospheric spectra are dominated by broad He i and Ca near-infrared emission lines with velocities of ∼10,000−12,000 km s ^−1 . The late-time spectra show prominent narrow He i emission lines at ∼1000 km s ^−1 , indicative of interaction with He-rich circumstellar material. SN 2023zaw is located in the spiral arm of a star-forming galaxy. We perform radiation-hydrodynamical and analytical modeling of the lightcurve by fitting with a combination of shock-cooling emission and nickel decay. The progenitor has a best-fit envelope mass of ≈0.2 M _☉ and an envelope radius of ≈50 R _⊙ . The extremely low nickel mass and low ejecta mass (≈0.5 M _⊙ ) suggest an ultrastripped SN, which originates from a mass-losing low-mass He-star (zero-age main-sequence mass < 10 M _⊙ ) in a close binary system. This is a channel to form double neutron star systems, whose merger is detectable with LIGO. SN 2023zaw underscores the existence of a previously undiscovered population of extremely low nickel mass (

We report observations of the optical counterpart of the long gamma-ray burst (GRB) GRB 230812B and its associated supernova (SN) SN 2023pel. The proximity ( z = 0.36) and high energy ( E _γ _,iso ∼ 10 ^53 erg) make it an important event to study as a probe of the connection between massive star core collapse and relativistic jet formation. With a phenomenological power-law model for the optical afterglow, we find a late-time flattening consistent with the presence of an associated SN. SN 2023pel has an absolute peak r -band magnitude of M _r = −19.46 ± 0.18 mag (about as bright as SN 1998bw) and evolves on quicker timescales. Using a radioactive heating model, we derive a nickel mass powering the SN of M _Ni = 0.38 ± 0.01 M _⊙ and a peak bolometric luminosity of L _bol ∼ 1.3 × 10 ^43 erg s ^−1 . We confirm SN 2023pel’s classification as a broad-line Type Ic SN with a spectrum taken 15.5 days after its peak in the r band and derive a photospheric expansion velocity of v _ph = 11,300 ± 1600 km s ^−1 at that phase. Extrapolating this velocity to the time of maximum light, we derive the ejecta mass M _ej = 1.0 ± 0.6 M _⊙ and kinetic energy ${E}_{\mathrm{KE}}={1.3}_{-1.2}^{+3.3}\times {10}^{51}\,\mathrm{erg}$ . We find that GRB 230812B/SN 2023pel has SN properties that are mostly consistent with the overall GRB-SN population. The lack of correlations found in the GRB-SN population between SN brightness and E _γ _,iso for their associated GRBs across a broad range of 7 orders of magnitude provides further evidence that the central engine powering the relativistic ejecta is not coupled to the SN powering mechanism in GRB-SN systems.

We present a systematic analysis of the X-ray emission of a sample of 17 optically selected, X-ray-detected tidal disruption events (TDEs) discovered between 2014 and 2021. The X-ray light curves show a diverse range of temporal behaviors, with most sources not following the expected power-law decline. The X-ray spectra are mostly extremely soft and consistent with thermal emission from the innermost region of an accretion disk, which cools as the accretion rate decreases. Three sources show formation of a hard X-ray corona at late times. The spectral energy distribution shape, probed by the ratio ( L _BB / L _X ) between the UV/optical and X-ray, shows a wide range of L _BB / L _X ∈ (0.5, 3000) at early times and converges to disklike values of L _BB / L _X ∈ (0.5, 10) at late times. We estimate the fraction of optically discovered TDEs with L _X ≥ 10 ^42 erg s ^−1 to be at least 40% and show that X-ray loudness is independent of black hole mass. We argue that distinct disk formation timescales are unlikely to be able to explain the diverse range of X-ray evolution. We combine our sample with X-ray-discovered ones to construct an X-ray luminosity function, best fit by a broken power law, with a break at L _X ≈ 10 ^44 erg s ^−1 . We show that there is no dichotomy between optically and X-ray-selected TDEs; instead, there is a continuum of early-time L _BB / L _X , at least as wide as L _BB / L _X ∈ (0.1, 3000), with optical/X-ray surveys selecting preferentially, but not exclusively, from the higher/lower end of the distribution. Our findings are consistent with unification models for the overall TDE population.

AT2022cmc was recently reported as the first on-axis jetted tidal disruption event (TDE) discovered in the last decade, and the fourth on-axis jetted TDE candidate known so far. In this work, we present NuSTAR hard X-ray (3–30 keV) observations of AT2022cmc, as well as soft X-ray (0.3–6 keV) observations obtained by NICER, Swift, and XMM-Newton. Our analysis reveals that the broadband X-ray spectra can be well described by a broken power law with f _ν ∝ ν ^−0.5 ( f _ν ∝ ν ^−1 ) below (above) the rest-frame break energy of E _bk ∼ 10 keV at the observer frame t _obs = 7.8 and 17.6 days since discovery. At t _obs = 36.2 days, the X-ray spectrum is consistent with either a single power law or a broken power law. By modeling the spectral energy distribution from radio to hard X-ray across the three NuSTAR observing epochs, we find that the submillimeter/radio emission originates from external shocks at large distances ≳10 ^17 cm from the black hole, the UV/optical light comes from a thermal envelope with radius ∼10 ^15 cm, and the X-ray emission is consistent with synchrotron radiation powered by energy dissipation at intermediate radii within the (likely magnetically dominated) jet. We constrain the bulk Lorentz factor of the jet to be of the order 10–100. Our interpretation differs from the model proposed by Pasham et al. where both the radio and X-rays come from the same emitting zone in a matter-dominated jet. Our model for the jet X-ray emission has broad implications on the nature of relativistic jets in other sources such as gamma-ray bursts.

One of the open questions following the discovery of GW170817 is whether neutron star (NS) mergers are the only astrophysical sites capable of producing r -process elements. Simulations have shown that 0.01–0.1 M _⊙ of r -process material could be generated in the outflows originating from the accretion disk surrounding the rapidly rotating black hole that forms as a remnant to both NS mergers and collapsing massive stars associated with long-duration gamma-ray bursts (collapsars). The hallmark signature of r -process nucleosynthesis in the binary NS merger GW170817 was its long-lasting near-infrared (NIR) emission, thus motivating a systematic photometric study of the light curves of broad-lined stripped-envelope (Ic-BL) supernovae (SNe) associated with collapsars. We present the first systematic study of 25 SNe Ic-BL—including 18 observed with the Zwicky Transient Facility and 7 from the literature—in the optical/NIR bands to determine what quantity of r -process material, if any, is synthesized in these explosions. Using semi-analytic models designed to account for r -process production in SNe Ic-BL, we perform light curve fitting to derive constraints on the r -process mass for these SNe. We also perform independent light curve fits to models without the r -process. We find that the r -process-free models are a better fit to the light curves of the objects in our sample. Thus, we find no compelling evidence of r -process enrichment in any of our objects. Further high-cadence infrared photometric studies and nebular spectroscopic analysis would be sensitive to smaller quantities of r -process ejecta mass or indicate whether all collapsars are completely devoid of r -process nucleosynthesis.

We conduct a systematic tidal disruption event (TDE) demographics analysis using the largest sample of optically selected TDEs. A flux-limited, spectroscopically complete sample of 33 TDEs is constructed using the Zwicky Transient Facility over 3 yr (from 2018 October to 2021 September). We infer the black hole (BH) mass ( M _BH ) with host galaxy scaling relations, showing that the sample M _BH ranges from 10 ^5.1 M _⊙ to 10 ^8.2 M _⊙ . We developed a survey efficiency corrected maximum volume method to infer the rates. The rest-frame g -band luminosity function can be well described by a broken power law of $\phi ({L}_{g})\propto {\left[{\left({L}_{g}/{L}_{\mathrm{bk}}\right)}^{0.3}+{\left({L}_{g}/{L}_{\mathrm{bk}}\right)}^{2.6}\right]}^{-1}$ , with L _bk = 10 ^43.1 erg s ^−1 . In the BH mass regime of 10 ^5.3 ≲ ( M _BH / M _⊙ ) ≲ 10 ^7.3 , the TDE mass function follows $\phi ({M}_{\mathrm{BH}})\propto {M}_{\mathrm{BH}}^{-0.25}$ , which favors a flat local BH mass function ( ${{dn}}_{\mathrm{BH}}/d\mathrm{log}{M}_{\mathrm{BH}}\approx \mathrm{constant}$ ). We confirm the significant rate suppression at the high-mass end ( M _BH ≳ 10 ^7.5 M _⊙ ), which is consistent with theoretical predictions considering direct capture of hydrogen-burning stars by the event horizon. At a host galaxy mass of M _gal ∼ 10 ^10 M _⊙ , the average optical TDE rate is ≈3.2 × 10 ^−5 galaxy ^−1 yr ^−1 . We constrain the optical TDE rate to be [3.7, 7.4, and 1.6] × 10 ^−5 galaxy ^−1 yr ^−1 in galaxies with red, green, and blue colors.

We report observations of the optical counterpart of the long gamma-ray burst GRB 221009A. Due to the extreme rarity of being both nearby ( z = 0.151) and highly energetic ( E _γ _,iso ≥ 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 supernova (SN) contributing excess flux to the optical light curve. Though not conclusive, we find moderate evidence ( K _Bayes = 10 ^1.2 ) for the presence of an additional component arising from an associated SN, SN 2022xiw, and find that it must be substantially fainter (

We report the discovery of the unusually bright long-duration gamma-ray burst (GRB), GRB 221009A, as observed by the Neil Gehrels Swift Observatory (Swift), Monitor of All-sky X-ray Image, and Neutron Star Interior Composition Explorer Mission. This energetic GRB was located relatively nearby ( z = 0.151), allowing for sustained observations of the afterglow. The large X-ray luminosity and low Galactic latitude ( b = 4.°3) make GRB 221009A a powerful probe of dust in the Milky Way. Using echo tomography, we map the line-of-sight dust distribution and find evidence for significant column densities at large distances (≳10 kpc). We present analysis of the light curves and spectra at X-ray and UV–optical wavelengths, and find that the X-ray afterglow of GRB 221009A is more than an order of magnitude brighter at T _0 + 4.5 ks than that from any previous GRB observed by Swift. In its rest frame, GRB 221009A is at the high end of the afterglow luminosity distribution, but not uniquely so. In a simulation of randomly generated bursts, only 1 in 10 ^4 long GRBs were as energetic as GRB 221009A; such a large E _γ _,iso implies a narrow jet structure, but the afterglow light curve is inconsistent with simple top-hat jet models. Using the sample of Swift GRBs with redshifts, we estimate that GRBs as energetic and nearby as GRB 221009A occur at a rate of ≲1 per 1000 yr—making this a truly remarkable opportunity unlikely to be repeated in our lifetime.

The host galaxies of tidal disruption events (TDEs) have been shown to possess peculiar properties, including high central light concentrations, unusual star formation histories, and “green” colors. The ubiquity of these large-scale galaxy characteristics among TDE host populations suggests that they may serve to boost the TDE rate in such galaxies by influencing the nuclear stellar dynamics. We present the first population study of integral field spectroscopy for 13 TDE host galaxies across all spectral classes and X-ray brightnesses with the purpose of investigating their large-scale properties. We derive the black hole masses via stellar kinematics (i.e., the M – σ relation) and find masses in the range $5.0\lesssim \mathrm{log}({M}_{\mathrm{BH}}/{M}_{\odot })\lesssim 8.0$ , with a distribution dominated by black holes with M _BH ∼ 10 ^6 M _⊙ . We find one object with M _BH ≳ 10 ^8 M _⊙ , above the “Hills mass”, which if the disrupted star was of solar type, allows a lower limit of a ≳ 0.16 to be placed on its spin, lending further support to the proposed connection between featureless TDEs and jetted TDEs. We also explore the level of rotational support in the TDE hosts, quantified by ( V / σ ) _e , a parameter that has been shown to correlate with the stellar age and may explain the peculiar host-galaxy preferences of TDEs. We find that the TDE hosts exhibit a broad range in ( V / σ ) _e following a similar distribution as E + A galaxies, which have been shown to be overrepresented among TDE host populations.

1ES1927+654 is a nearby active galactic nucleus (AGN) that has shown an enigmatic outburst in optical/UV followed by X-rays, exhibiting strange variability patterns at timescales of months to years. Here we report the unusual X-ray, UV, and radio variability of the source in its postflare state (2022 January–2023 May). First, we detect an increase in the soft X-ray (0.3–2 keV) flux from 2022 May to 2023 May by almost a factor of 5, which we call the bright soft state. The hard X-ray 2–10 keV flux increased by a factor of 2, while the UV flux density did not show any significant changes (≤30%) in the same period. The integrated energy pumped into the soft and hard X-rays during this period of 11 months is ∼3.57 × 10 ^50 erg and 5.9 × 10 ^49 erg, respectively. From the energetics, it is evident that whatever is producing the soft excess (SE) is pumping out more energy than either the UV or hard X-ray source. Since the energy source presumably is ultimately the accretion of matter onto the supermassive black hole, the SE-emitting region must be receiving the majority of this energy. In addition, the source does not follow the typical disk–corona relation found in AGNs, neither in the initial flare (from 2017 to 2019) nor in the current bright soft state (2022–2023). We found that the core (

Despite their discovery about half a century ago, the gamma-ray burst (GRB) prompt emission mechanism is still not well understood. Theoretical modeling of the prompt emission has advanced considerably due to new computational tools and techniques. One such tool is the PLUTO hydrodynamics code, which is used to numerically simulate GRB outflows. PLUTO uses Adaptive Mesh Refinement to focus computational efforts on the portion of the grid that contains the simulated jet. Another tool is the Monte Carlo Radiation Transfer (MCRaT) code, which predicts the electromagnetic signatures of GRBs by conducting photon scatterings within a jet using PLUTO. The effects of the underlying resolution of a PLUTO simulation with respect to MCRaT post-processing radiative transfer results have not yet been quantified. We analyze an analytic spherical outflow and a hydrodynamically simulated GRB jet with MCRaT at varying spatial and temporal resolutions and quantify how decreasing both resolutions affects the resulting mock observations. We find that changing the spatial resolution changes the hydrodynamic properties of the jet, which directly affect the MCRaT mock observable peak energies. We also find that decreasing the temporal resolution artificially decreases the high-energy slope of the mock observed spectrum, which increases both the spectral peak energy and the luminosity. We show that the effects are additive when both spatial and temporal resolutions are modified. Our results allow us to understand how decreased hydrodynamic temporal and spatial resolutions affect the results of post-processing radiative transfer calculations, allowing for the optimization of hydrodynamic simulations for radiative transfer codes.

The dividing line between gamma-ray bursts (GRBs) and ordinary stripped-envelope core-collapse supernovae (SNe) is yet to be fully understood. Observationally mapping the variety of ejecta outcomes (ultrarelativistic, mildly relativistic, or nonrelativistic) in SNe of Type Ic with broad lines (Ic-BL) can provide a key test to stellar explosion models. However, this requires large samples of the rare SN Ic-BL events with follow-up observations in the radio, where fast ejecta can be probed largely free of geometry and viewing angle effects. Here, we present the results of a radio (and X-ray) follow-up campaign of 16 SNe Ic-BL detected by the Zwicky Transient Facility (ZTF). Our radio campaign resulted in four counterpart detections and 12 deep upper limits. None of the events in our sample is as relativistic as SN 1998bw and we constrain the fraction of SN 1998bw-like explosions to

The Swift Burst Alert Telescope (BAT) is a coded-aperture gamma-ray instrument with a large field of view that primarily operates in survey mode when it is not triggering on transient events. The survey data consist of 80-channel detector plane histograms that accumulate photon counts over periods of at least 5 minutes. These histograms are processed on the ground and are used to produce the survey data set between 14 and 195 keV. Survey data comprise >90% of all BAT data by volume and allow for the tracking of long-term light curves and spectral properties of cataloged and uncataloged hard X-ray sources. Until now, the survey data set has not been used to its full potential due to the complexity associated with its analysis and the lack of easily usable pipelines. Here, we introduce the BatAnalysis Python package, a wrapper for HEASoftpy, which provides a modern, open-source pipeline to process and analyze BAT survey data. BatAnalysis allows members of the community to use BAT survey data in more advanced analyses of astrophysical sources, including pulsars, pulsar wind nebula, active galactic nuclei, and other known/unknown transient events that may be detected in the hard X-ray band. We outline the steps taken by the Python code and exemplify its usefulness and accuracy by analyzing survey data of the Crab Nebula, NGC 2992, and a previously uncataloged MAXI transient. The BatAnalysis package allows for ∼18 yr of BAT survey data to be used in a systematic way to study a large variety of astrophysical sources.

Narrow transient emission lines (flash-ionization features) in early supernova (SN) spectra trace the presence of circumstellar material (CSM) around the massive progenitor stars of core-collapse SNe. The lines disappear within days after the SN explosion, suggesting that this material is spatially confined, and originates from enhanced mass loss shortly (months to a few years) prior to the explosion. We performed a systematic survey of H-rich (Type II) SNe discovered within less than 2 days from the explosion during the first phase of the Zwicky Transient Facility survey (2018–2020), finding 30 events for which a first spectrum was obtained within 36% at the 95% confidence level) confirms that elevated mass loss in massive stars prior to SN explosion is common. We find that SNe II showing flash-ionization features are not significantly brighter, nor bluer, nor more slowly rising than those without. This implies that CSM interaction does not contribute significantly to their early continuum emission, and that the CSM is likely optically thin. We measured the persistence duration of flash-ionization emission and find that most SNe show flash features for ≈5 days. Rarer events, with persistence timescales >10 days, are brighter and rise longer, suggesting these may be intermediate between regular SNe II and strongly interacting SNe IIn.

SN 2018aoz is a Type Ia SN with a B -band plateau and excess emission in infant-phase light curves ≲1 day after the first light, evidencing an over-density of surface iron-peak elements as shown in our previous study. Here, we advance the constraints on the nature and origin of SN 2018aoz based on its evolution until the nebular phase. Near-peak spectroscopic features show that the SN is intermediate between two subtypes of normal Type Ia: core normal and broad line. The excess emission may be attributable to the radioactive decay of surface iron-peak elements as well as the interaction of ejecta with either the binary companion or a small torus of circumstellar material. Nebular-phase limits on H α and He i favor a white dwarf companion, consistent with the small companion size constrained by the low early SN luminosity, while the absence of [O i ] and He i disfavors a violent merger of the progenitor. Of the two main explosion mechanisms proposed to explain the distribution of surface iron-peak elements in SN 2018aoz, the asymmetric Chandrasekhar-mass explosion is less consistent with the progenitor constraints and the observed blueshifts of nebular-phase [Fe ii ] and [Ni ii ]. The helium-shell double-detonation explosion is compatible with the observed lack of C spectral features, but current 1D models are incompatible with the infant-phase excess emission, ${B}_{\max }\mbox{--}{V}_{\max }$ color, and weak strength of nebular-phase [Ca ii ]. Although the explosion processes of SN 2018aoz still need to be more precisely understood, the same processes could produce a significant fraction of Type Ia SNe that appear to be normal after ∼1 day.

The detection of gravitational waves from the binary neuron star merger GW170817 and electromagnetic counterparts GRB170817A and AT2017gfo kick-started the field of gravitational-wave multimessenger astronomy. The optically red to near-infrared emission (“red” component) of AT2017gfo was readily explained as produced by the decay of newly created nuclei produced by rapid neutron capture (a kilonova). However, the ultraviolet to optically blue emission (“blue” component) that was dominant at early times (up to 1.5 days) received no consensus regarding its driving physics. Among many explanations, two leading contenders are kilonova radiation from a lanthanide-poor ejecta component and shock interaction (cocoon emission). In this work, we simulate AT2017gfo-like light curves and perform a Bayesian analysis to study whether an ultraviolet satellite capable of rapid gravitational-wave follow-up, could distinguish between physical processes driving the early “blue” component. We find that ultraviolet data starting at 1.2 hr distinguishes the two early radiation models up to 160 Mpc, implying that an ultraviolet mission like Dorado would significantly contribute to insights into the driving emission physics of the postmerger system. While the same ultraviolet data and optical data starting at 12 hr have limited ability to constrain model parameters separately, the combination of the two unlocks tight constraints for all but one parameter of the kilonova model up to 160 Mpc. We further find that a Dorado-like ultraviolet satellite can distinguish the early radiation models up to at least 130 (60) Mpc if data collection starts within 3.2 (5.2) hr for AT2017gfo-like light curves.

The BAT AGN Spectroscopic Survey (BASS) is designed to provide a highly complete census of the key physical parameters of the supermassive black holes (SMBHs) that power local active galactic nuclei (AGNs) ( z ≲ 0.3), including their bolometric luminosity ( L _bol ), black hole (BH) mass ( M _BH ), accretion rates ( L _bol / L _Edd ), line-of-sight gas obscuration ( N _H ), and the distinctive properties of their host galaxies (e.g., star formation rates, masses, and gas fractions). We present an overview of the second data release of BASS (DR2), an unprecedented spectroscopic AGN survey in spectral range, resolution, and sensitivity, including 1449 optical (∼3200 Å–1 μ m) and 233 near-IR (1–2.5 μ m) spectra for the brightest 858 ultrahard X-ray (14–195 keV) selected AGNs across the entire sky and essentially all levels of obscuration. This release provides a highly complete set of key measurements (emission-line measurements and central velocity dispersions), with 99.9% measured redshifts and 98% BH masses estimated (for unbeamed AGNs outside the Galactic plane). The BASS DR2 AGN sample represents a unique census of nearby powerful AGNs, spanning over 5 orders of magnitude in AGN bolometric luminosity ( L _bol ∼ 10 ^40 –10 ^47 erg s ^−1 ), BH mass ( M _BH ∼ 10 ^5 –10 ^10 M _⊙ ), Eddington ratio ( L _bol / L _Edd ≳ 10 ^−5 ), and obscuration ( N _H ∼ 10 ^20 –10 ^25 cm ^−2 ). The public BASS DR2 sample and measurements can thus be used to answer fundamental questions about SMBH growth and its links to host galaxy evolution and feedback in the local universe, as well as open questions concerning SMBH physics. Here we provide a brief overview of the survey strategy, the key BASS DR2 measurements, data sets and catalogs, and scientific highlights from a series of DR2-based works pursued by the BASS team.

Tidal disruption events (TDEs) offer a unique way to study dormant black holes. While the number of observed TDEs has grown thanks to the emergence of wide-field surveys in the past few decades, questions regarding the nature of the observed optical, UV, and X-ray emission remain. We present a uniformly selected sample of 30 spectroscopically classified TDEs from the Zwicky Transient Facility Phase I survey operations with follow-up Swift UV and X-ray observations. Through our investigation into correlations between light-curve properties, we recover a shallow positive correlation between the peak bolometric luminosity and decay timescales. We introduce a new spectroscopic class of TDE, TDE-featureless, which are characterized by featureless optical spectra. The new TDE-featureless class shows larger peak bolometric luminosities, peak blackbody temperatures, and peak blackbody radii. We examine the differences between the X-ray bright and X-ray faint populations of TDEs in this sample, finding that X-ray bright TDEs show higher peak blackbody luminosities than the X-ray faint subsample. This sample of optically selected TDEs is the largest sample of TDEs from a single survey yet, and the systematic discovery, classification, and follow-up of this sample allows for robust characterization of TDE properties, an important stepping stone looking forward toward the Rubin era.