We report observations of the ultra-high-energy gamma-ray source LHAASO J2108$+$5157, utilizing VERITAS, HAWC, Fermi-LAT, and XMM-Newton. VERITAS has collected $\sim$ 40 hours of data that we used to set ULs to the emission above 200 GeV. The HAWC data, collected over $\sim 2400$ days, reveal emission between 3 and 146 TeV, with a significance of $7.5~\sigma$, favoring an extended source model. The best-fit spectrum measured by HAWC is characterized by a simple power-law with a spectral index of $2.45\pm0.11_{stat}$. Fermi-LAT analysis finds a point source with a very soft spectrum in the LHAASO J2108+5157 region, consistent with the 4FGL-DR3 catalog results. The XMM-Newton analysis yields a null detection of the source in the 2 - 7 keV band. The broadband spectrum can be interpreted as a pulsar and a pulsar wind nebula system, where the GeV gamma-ray emission originates from an unidentified pulsar, and the X-ray and TeV emission is attributed to synchrotron radiation and inverse Compton scattering of electrons accelerated within a pulsar wind nebula. In this leptonic scenario, our X-ray upper limit provides a stringent constraint on the magnetic field, which is $\lesssim 1.5\ \mu$G.

(Abridged) A simple six-parameter LCDM model provides a successful fit to WMAP data, both when the data are analyzed alone and in combination with other cosmological data. Even so, it is appropriate to search for any hints of deviations from the now standard model of cosmology, which includes inflation, dark energy, dark matter, baryons, and neutrinos. The cosmological community has subjected the WMAP data to extensive and varied analyses. While there is widespread agreement as to the overall success of the six-parameter LCDM model, various "anomalies" have been reported relative to that model. In this paper we examine potential anomalies and present analyses and assessments of their significance. In most cases we find that claimed anomalies depend on posterior selection of some aspect or subset of the data. Compared with sky simulations based on the best fit model, one can select for low probability features of the WMAP data. Low probability features are expected, but it is not usually straightforward to determine whether any particular low probability feature is the result of the a posteriori selection or of non-standard cosmology. We examine in detail the properties of the power spectrum with respect to the LCDM model. We examine several potential or previously claimed anomalies in the sky maps and power spectra, including cold spots, low quadrupole power, quadropole-octupole alignment, hemispherical or dipole power asymmetry, and quadrupole power asymmetry. We conclude that there is no compelling evidence for deviations from the LCDM model, which is generally an acceptable statistical fit to WMAP and other cosmological data.

We investigate the resolved properties of star-forming clumps and their host galaxies at $0.5

Entity is a new-generation, fully open-source particle-in-cell (PIC) code developed to overcome key limitations in astrophysical plasma modeling, particularly the extreme separation of scales and the performance challenges associated with evolving, GPU-centric computing infrastructures. It achieves hardware-agnostic performance portability across various GPU and CPU architectures using the Kokkos library. Crucially, Entity maintains a high standard for usability, clarity, and customizability, offering a robust and easy-to-use framework for developing new algorithms and grid geometries, which allows extensive control without requiring edits to the core source code. This paper details the core general-coordinate special-relativistic module. Entity is the first PIC code designed to solve the Vlasov-Maxwell system in general coordinates, enabling a coordinate-agnostic framework that provides the foundational structure for straightforward extension to arbitrary coordinate geometries. The core methodology achieves numerical stability by solving particle equations of motion in the global orthonormal Cartesian basis, despite using generalized coordinates like Cartesian, axisymmetric spherical, and quasi-spherical grids. Charge conservation is ensured via a specialized current deposition technique using conformal currents. The code exhibits robust scalability and performance portability on major GPU platforms (AMD MI250X, NVIDIA A100, and Intel Max Series), with the 3D particle pusher and the current deposition operating efficiently at about 2 nanoseconds per particle per timestep. Functionality is validated through a comprehensive suite of standard Cartesian plasma tests and the accurate modeling of relativistic magnetospheres in curvilinear axisymmetric geometries.

(Abridged) We present the angular power spectra derived from the 7-year maps and discuss the cosmological conclusions that can be inferred from WMAP data alone. The third acoustic peak in the TT spectrum is now well measured by WMAP. In the context of a flat LambdaCDM model, this improvement allows us to place tighter constraints on the matter density from WMAP data alone, and on the epoch of matter-radiation equality, The temperature-polarization (TE) spectrum is detected in the 7-year data with a significance of 20 sigma, compared to 13 sigma with the 5-year data. The low-l EE spectrum, a measure of the optical depth due to reionization, is detected at 5.5 sigma significance when averaged over l = 2-7. The BB spectrum, an important probe of gravitational waves from inflation, remains consistent with zero. The upper limit on tensor modes from polarization data alone is a factor of 2 lower with the 7-year data than it was using the 5-year data (Komatsu et al. 2010). We test the parameter recovery process for bias and find that the scalar spectral index, ns, is biased high, but only by 0.09 sigma, while the remaining parameters are biased by < 0.15 sigma. The improvement in the third peak measurement leads to tighter lower limits from WMAP on the number of relativistic degrees of freedom (e.g., neutrinos) in the early universe: Neff > 2.7 (95% CL). Also, using WMAP data alone, the primordial helium mass fraction is found to be YHe = 0.28+0.14-0.15, and with data from higher-resolution CMB experiments included, we now establish the existence of pre-stellar helium at > 3 sigma (Komatsu et al. 2010).

The cooling envelope model for tidal disruption events (TDE) postulates that while the stellar debris streams rapidly dissipate their bulk kinetic energy (``circularize"), this does not necessarily imply rapid feeding of the supermassive black hole (SMBH). The bound material instead forms a large pressure-supported envelope which powers optical/UV emission as it undergoes gradual Kelvin-Helmholtz contraction. We present results interpreting a sample of 15 optical TDE within the cooling envelope model in order to constrain the SMBH mass $M_{\rm BH}$, stellar mass $M_{\star}$, and orbital penetration factor $\beta$. The distributions of inferred properties from our sample broadly follow the theoretical expectations of loss-cone analysis assuming a standard stellar initial mass function. However, we find a deficit of events with $M_{\rm BH} \lesssim 5\times 10^{5}M_{\odot}$ and $M_{\star} \lesssim 0.5M_{\odot}$, which could result in part from the reduced detectability of TDEs with these properties. Our model fits also illustrate the predicted long delay between the optical light curve peak and when the SMBH accretion rate reaches its maximum. The latter occurs only once the envelope contracts to the circularization radius on a timescale of months to years, consistent with delayed-rising X-ray and non-thermal radio flares seen in a growing number of TDE.

Supermassive black hole (BH) binaries are thought to produce self-lensing flares (SLF) when the two BHs are aligned with the line-of-sight. If the binary orbit is observed nearly edge-on, we find a distinct feature in the light curve imprinted by the BH shadow from the lensed BH. We study this feature by ray-tracing in a binary model and predict that 1\% of the current binary candidates could show this feature. Our BH tomography method proposed here could make it possible to extract BH shadows that are spatially unresolvable by high-resolution VLBI.

We present observations of the occulted active region AR12222 during the third {\em NuSTAR} solar campaign on 2014 December 11, with concurrent {\em SDO/}AIA and {\em FOXSI-2} sounding rocket observations. The active region produced a medium size solar flare one day before the observations, at $\sim18$UT on 2014 December 10, with the post-flare loops still visible at the time of {\em NuSTAR} observations. The time evolution of the source emission in the {\em SDO/}AIA $335\textrmÅ$ channel reveals the characteristics of an extreme-ultraviolet late phase event, caused by the continuous formation of new post-flare loops that arch higher and higher in the solar corona. The spectral fitting of {\em NuSTAR} observations yields an isothermal source, with temperature $3.8-4.6$ MK, emission measure $0.3-1.8 \times 10^{46}\textrm{ cm}^{-3}$, and density estimated at $2.5-6.0 \times 10^8 \textrm{ cm}^{-3}$. The observed AIA fluxes are consistent with the derived {\em NuSTAR} temperature range, favoring temperature values in the range $4.0-4.3$ MK. By examining the post-flare loops' cooling times and energy content, we estimate that at least 12 sets of post-flare loops were formed and subsequently cooled between the onset of the flare and {\em NuSTAR} observations, with their total thermal energy content an order of magnitude larger than the energy content at flare peak time. This indicates that the standard approach of using only the flare peak time to derive the total thermal energy content of a flare can lead to a large underestimation of its value.

Periodic collisions between a star on an inclined orbit around a supermassive black hole and its accretion disk offers a promising explanation for X-ray "quasi-periodic eruptions" (QPEs). Each passage through the disk shocks and compresses gas ahead of the star, which subsequently re-expands above the disk as a quasi-spherical cloud. We present spherically symmetric Monte Carlo radiation transport simulations which follow the production of photons behind the radiation-mediated shock, Comptonization by hot electrons, and the eventual escape of the radiation through the expanding debris. Such one-dimension calculations are approximately justified for thin disks, through which the star of radius $R_{\star}$ passes faster than the shocked gas can flow around the star. For collision speeds $v_{\rm coll} \gtrsim 0.15 c$ and disk surface densities $\Sigma \sim 10^{3}$ g cm$^{-2}$ characteristic of those encountered by stellar orbits consistent with QPE recurrence times, the predicted transient light curves exhibit peak luminosities $\gtrsim 10^{42}$ erg s$^{-1}$ and Comptonized quasi-thermal (Wien-like) spectra which peak at energies $h\nu \sim 100$ eV, broadly consistent with QPE properties. For these conditions, gas and radiation are out of equilibrium and the emission temperature is harder than the blackbody value due to inefficient photon production behind the shock. Alternatively, for higher disk densities and/or lower shock velocities, QPE emission could instead represent the comparatively brief phase shortly after shock break-out, though in this case the bulk of the radiation is thermalized and occurs in the ultraviolet instead of the X-ray band. In either scenario, reproducing the observed eruption properties (duration, luminosity, temperature) requires a large radius $R_{\star} \gtrsim 10R_{\odot}$, which may point to inflation of the star's atmosphere from repeated collisions.

(Abridged) New full sky temperature and polarization maps based on seven years of data from WMAP are presented. The new results are consistent with previous results, but have improved due to reduced noise from the additional integration time, improved knowledge of the instrument performance, and improved data analysis procedures. The improvements are described in detail. The seven year data set is well fit by a minimal six-parameter flat Lambda-CDM model. The parameters for this model, using the WMAP data in conjunction with baryon acoustic oscillation data from the Sloan Digital Sky Survey and priors on H_0 from Hubble Space Telescope observations, are: Omega_bh^2 = 0.02260 +-0.00053, Omega_ch^2 = 0.1123 +-0.0035, Omega_Lambda = 0.728 +0.015 -0.016, n_s = 0.963 +-0.012, tau = 0.087 +-0.014 and sigma_8 = 0.809 +-0.024 (68 % CL uncertainties). The temperature power spectrum signal-to-noise ratio per multipole is greater that unity for multipoles < 919, allowing a robust measurement of the third acoustic peak. This measurement results in improved constraints on the matter density, Omega_mh^2 = 0.1334 +0.0056 -0.0055, and the epoch of matter- radiation equality, z_eq = 3196 +134 -133, using WMAP data alone. The new WMAP data, when combined with smaller angular scale microwave background anisotropy data, results in a 3 sigma detection of the abundance of primordial Helium, Y_He = 0.326 +-0.075.The power-law index of the primordial power spectrum is now determined to be n_s = 0.963 +-0.012, excluding the Harrison-Zel'dovich-Peebles spectrum by >3 sigma. These new WMAP measurements provide important tests of Big Bang cosmology.

Sub-arcsecond localization of the repeating fast radio burst FRB 121102 revealed its coincidence with a dwarf host galaxy and a steady (`quiescent') non-thermal radio source. We show that the properties of the host galaxy are consistent with those of long-duration gamma-ray bursts (LGRB) and hydrogen-poor superluminous supernovae (SLSNe-I). Both LGRBs and SLSNe-I were previously hypothesized to be powered by the electromagnetic spin-down of newly-formed, strongly-magnetized neutron stars with millisecond birth rotation periods (`millisecond magnetars'). This motivates considering a scenario whereby the repeated bursts from FRB 121102 originate from a young magnetar remnant embedded within a young hydrogen-poor supernova remnant. Requirements on the GHz free-free optical depth through the expanding supernova ejecta (accounting for photo-ionization by the rotationally-powered magnetar nebula), energetic constraints on the bursts, and constraints on the size of the quiescent source all point to an age of less than a few decades to a century. The quiescent radio source can be attributed to synchrotron emission from the shock interaction between the fast outer layer of the supernova ejecta with the surrounding wind of the progenitor star, or from deeper within the magnetar wind nebula. Alternatively, the radio emission could be an orphan afterglow from an initially off-axis LGRB jet, though this might require the source to be too young. The young age of the source can be tested by searching for a time derivative of the dispersion measure and predicted fading of the quiescent radio source. We propose future tests of the SLSNe-I/LGRB/FRB connection, such as searches for FRBs from nearby SLSNe-I/LGRB on timescales of decades after their explosions.

(Abridged) The WMAP 5-year data strongly limit deviations from the minimal LCDM model. We constrain the physics of inflation via Gaussianity, adiabaticity, the power spectrum shape, gravitational waves, and spatial curvature. We also constrain the properties of dark energy, parity-violation, and neutrinos. We detect no convincing deviations from the minimal model. The parameters of the LCDM model, derived from WMAP combined with the distance measurements from the Type Ia supernovae (SN) and the Baryon Acoustic Oscillations (BAO), are: Omega_b=0.0456+-0.0015, Omega_c=0.228+-0.013, Omega_Lambda=0.726+-0.015, H_0=70.5+-1.3 km/s/Mpc, n_s=0.960+-0.013, tau=0.084+-0.016, and sigma_8=0.812+-0.026. With WMAP+BAO+SN, we find the tensor-to-scalar ratio r<0.22 (95% CL), and n_s>1 is disfavored regardless of r. We obtain tight, simultaneous limits on the (constant) equation of state of dark energy and curvature. We provide a set of "WMAP distance priors," to test a variety of dark energy models. We test a time-dependent w with a present value constrained as -0.33<1+w_0<0.21 (95% CL). Temperature and matter fluctuations obey the adiabatic relation to within 8.9% and 2.1% for the axion and curvaton-type dark matter, respectively. The TE and EB spectra constrain cosmic parity-violation. We find the limit on the total mass of neutrinos, sum(m_nu)<0.67 eV (95% CL), which is free from the uncertainty in the normalization of the large-scale structure data. The effective number of neutrino species is constrained as N_{eff} = 4.4+-1.5 (68%), consistent with the standard value of 3.04. Finally, limits on primordial non-Gaussianity are -9

When a star passes close to a supermassive black hole (BH), the BH's tidal forces rip it apart into a thin stream, leading to a tidal disruption event (TDE). In this work, we study the post-disruption phase of TDEs in general relativistic hydrodynamics (GRHD) using our GPU-accelerated code H-AMR. We carry out the first grid-based simulation of a deep-penetration TDE ($\beta$=7) with realistic system parameters: a black-hole-to-star mass ratio of $10^6$, a parabolic stellar trajectory, and a nonzero BH spin. We also carry out a simulation of a tilted TDE whose stellar orbit is inclined relative to the BH midplane. We show that for our aligned TDE, an accretion disk forms due to the dissipation of orbital energy with $\sim$20 percent of the infalling material reaching the BH. The dissipation is initially dominated by violent self-intersections and later by stream-disk interactions near the pericenter. The self-intersections completely disrupt the incoming stream, resulting in five distinct self-intersection events separated by approximately 12 hours and a flaring in the accretion rate. We also find that the disk is eccentric with mean eccentricity e$\approx$0.88. For our tilted TDE, we find only partial self-intersections due to nodal precession near pericenter. Although these partial intersections eject gas out of the orbital plane, an accretion disk still forms with a similar accreted fraction of the material to the aligned case. These results have important implications for disk formation in realistic tidal disruptions. For instance, the periodicity in accretion rate induced by the complete stream disruption may explain the flaring events from Swift J1644+57.

We find line emission from the hydrogen- and/or helium-like ions of Ne, O, N and C in the low and short-on states of Her X-1, using the XMM-Newton Reflection Grating Spectrometer. The emission line velocity broadening is 200 < sigma < 500 km/s. Plasma diagnostics with the Ne IX, O VII and N VI He-alpha lines and the radiative recombination continua of O VII and N VII, indicate the gas is heated by photoionization. We use spectral models to measure the element abundance ratios N/O, C/O, and Ne/O, which quantify CNO processing in HZ Her. Photoexcitation and high-density effects are not differentiated by the measured He-alpha lines. We set limits on the location, temperature and density of the line emission region. The narrow emission lines can be attributed to reprocessing in either an accretion disk atmosphere and corona or on the X-ray illuminated face of HZ Her. In the main-on state, the bright continuum only allows the detection of interstellar absorption, plus O VII He-alpha emission lines with sigma = 3200 +- 700 km/s and complex profiles. Other broad lines may be present. The broad lines may originate in a region near the pulsar magnetosphere. Fe L lines are not detected.

(Abridged) We present WMAP seven-year observations of bright sources which are often used as calibrators at microwave frequencies. Ten objects are studied in five frequency bands (23 - 94 GHz): the outer planets (Mars, Jupiter, Saturn, Uranus and Neptune) and five fixed celestial sources (Cas A, Tau A, Cyg A, 3C274 and 3C58). The seven-year analysis of Jupiter provides temperatures which are within 1-sigma of the previously published WMAP five-year values, with slightly tighter constraints on variability with orbital phase, and limits (but no detections) on linear polarization. Scaling factors are provided which, when multiplied by the Wright Mars thermal model predictions at 350 micron, reproduce WMAP seasonally averaged observations of Mars within ~2%. An empirical model is described which fits brightness variations of Saturn due to geometrical effects and can be used to predict the WMAP observations to within 3%. Seven-year mean temperatures for Uranus and Neptune are also tabulated. Uncertainties in Uranus temperatures are 3%-4% in the 41, 61 and 94 GHz bands; the smallest uncertainty for Neptune is ~8% for the 94 GHz band. Intriguingly, the spectrum of Uranus appears to show a dip at ~30 GHz of unidentified origin, although the feature is not of high statistical significance. Flux densities for the five selected fixed celestial sources are derived from the seven-year WMAP sky maps, and are tabulated for Stokes I, Q and U, along with polarization fraction and position angle. Fractional uncertainties for the Stokes I fluxes are typically 1% to 3%. Source variability over the seven-year baseline is also estimated. Significant secular decrease is seen for Cas A and Tau A: our results are consistent with a frequency independent decrease of about 0.53% per year for Cas A and 0.22% per year for Tau A.

M33 contains a large number of emission nebulae identified as supernova remnants (SNRs) based on the high [S II]:Ha ratios characteristic of shocked gas. Using Chandra data from the ChASeM33 survey with a 0.35-2 keV sensitivity of about 2 x 10**34 ergs/s, we have detected 82 of 137 SNR candidates, yielding confirmation of (or at least strongly support for) their SNR identifications. This provides the largest sample of remnants detected at optical and X-ray wavelengths in any galaxy, including the Milky Way. A spectral analysis of the seven X-ray brightest SNRs reveals that two, G98-31 and G98-35, have spectra that appear to indicate enrichment by ejecta from core-collapse supernova explosions. In general, the X-ray detected SNRs have soft X-ray spectra compared to the vast majority of sources detected along the line of sight to M33. It is unlikely that there are any other undiscovered thermally dominated X-ray SNRs with luminosities in excess of about 4 x 10**35 ergs/s in the portions of M33 covered by the ChASeM33 survey. We have used a combination of new and archival optical and radio observations to attempt to better understand why some objects are detected as X-ray sources and others are not. We have also developed a morphological classification scheme for the optically-identified SNRs, and discuss the efficacy of this scheme as a predictor of X-ray detectability. Finally, we have compared the SNRs found in M33 to those that have been observed in the Galaxy and the Magellanic Clouds. There are no close analogs of Cas A, Kepler's SNR, Tycho's SNR or the Crab Nebula in the regions of M33 surveyed, but we have found an X-ray source with a power law spectrum coincident with a small-diameter radio source that may be the first pulsar-wind nebula recognized in M33.

The first high-energy neutrino source identified by IceCube was a blazar -- an active galactic nucleus driving a relativistic jet towards Earth. Jets driven by accreting black holes are commonly assumed to be needed for high-energy neutrino production. Recently, IceCube discovered neutrinos from Seyfert galaxies, which appears unrelated to jet activity. Here, we show that the observed luminosity ratios of neutrinos and hard X-rays from blazars TXS 0506+056 and GB6 J1542+6129 are consistent with neutrino production in a $\gamma$-obscured region near a central supermassive black hole, with the X-ray flux corresponding to reprocessed $\gamma$-ray emission with flux comparable to that of neutrinos. Similar neutrino - hard X-ray flux ratios are found for four Seyfert galaxies, NGC 1068, NGC 4151, CGCG 420-015 and NGC 3079, raising the possibility of a common neutrino production mechanism that may not involve a strong jet.

The discovery by Advanced LIGO/Virgo of gravitational waves from the binary neutron star (NS) merger GW170817, and subsequently by astronomers of transient counterparts across the electromagnetic (EM) spectrum, has initiated the era of multi-messenger astronomy. Given the slew of papers appearing on this event, I thought it useful to summarize the EM discoveries in the context of theoretical models and present my views on the major take-away lessons from this watershed event. The weak GRB discovered in close time coincidence with GW170817, and potential evidence for a more powerful off-axis relativistic jet (initially beamed away from our line of sight) from the delayed rise of a non-thermal X-ray and radio orphan afterglow, provides the most compelling evidence yet that cosmological short GRBs originate from binary NS mergers. The luminosity and colors of the early optical emission discovered within a day of the merger agrees strikingly well with original predictions (Metzger et al. 2010) for "kilonova" emission powered by the radioactive decay of r-process nuclei, the NS merger origin of which was initially proposed by Lattimer & Schramm 1974. The transition of the spectral energy distribution to NIR wavelengths on timescales of days matches predictions by Barnes & Kasen 2013 and Tanaka & Hotokezaka 2013 if a portion of the ejecta contains heavy r-process nuclei with higher opacities. The "blue" and "red" ejecta components may possess distinct origins (e.g. dynamical ejecta versus accretion disk outflows), with key implications for the merger physics and the properties of neutron stars. I outline the diversity in the EM emission expected from additional mergers-observed with different binary masses and viewing angles-discovered once LIGO/Virgo reach design sensitivity and NS mergers are discovered as frequently as once per week.

In a survey of quasar candidates selected by matching the FIRST and 2MASS catalogs, we have found two extraordinarily red quasars. FIRST J013435.7-093102 is a 1 Jy source at z=2.216 and has B-K > 10, while FIRST J073820.1+275045 is a 2.5 mJy source at z=1.985 with B-K = 8.4. FIRST J073820.1+275045 has strong absorption lines of MgII and CIV in the rest frame of the quasar and is highly polarized in the rest frame ultraviolet, strongly favoring the interpretation that its red spectral energy distribution is caused by dust reddening local to the quasar. FIRST J073820.1+275045 is thus one of the few low radio-luminosity, highly dust-reddened quasars known. The available observational evidence for FIRST J013435.7-093102 leads us to conclude that it too is reddened by dust. We show that FIRST J013435.7-093102 is gravitationally lensed, increasing the number of known lensed, extremely dust-reddened quasars to at least three, including MG0414-0534 and PKS1830-211. We discuss the implications of whether these objects are reddened by dust in the host or lensing galaxies. If reddened by their local environment, then we estimate that between 10 and 20% of the radio-loud quasar population is reddened by dust in the host galaxy. The discovery of FIRST J073820.1+275045 and objects now emerging from X-ray surveys suggests the existence of an analogous radio-quiet red quasar population. Such objects will be entirely missed by standard radio or optical quasar surveys. If dust in the lensing galaxies is primarily responsible for the extreme redness of the lensed quasars, then an untold number of gravitationally lensed quasars are being overlooked.

Einstein@Home aggregates the computer power of hundreds of thousands of volunteers from 193 countries, to search for new neutron stars using data from electromagnetic and gravitational-wave detectors. This paper presents a detailed description of the search for new radio pulsars using Pulsar ALFA survey data from the Arecibo Observatory. The enormous computing power allows this search to cover a new region of parameter space; it can detect pulsars in binary systems with orbital periods as short as 11 minutes. We also describe the first Einstein@Home discovery, the 40.8 Hz isolated pulsar PSR J2007+2722, and provide a full timing model. PSR J2007+2722's pulse profile is remarkably wide with emission over almost the entire spin period. This neutron star is most likely a disrupted recycled pulsar, about as old as its characteristic spin-down age of 404 Myr. However there is a small chance that it was born recently, with a low magnetic field. If so, upper limits on the X-ray flux suggest but can not prove that PSR J2007+2722 is at least ~ 100 kyr old. In the future, we expect that the massive computing power provided by volunteers should enable many additional radio pulsar discoveries.