Transient Survey Research Articles

Transient electromagnetic data denoising based on cluster analysis and locally weighted linear regression

AbstractIn transient electromagnetic surveys, the collected data inevitably contain noise originating from both natural and cultural sources. This noise has the potential to mask transient electromagnetic responses linked to geological features, thereby posing challenges in accurately interpreting subsurface structures. Hence, the implementation of effective noise reduction techniques is crucial in ensuring the accuracy and reliability of inversion outcomes in transient electromagnetic surveys. This study introduces a novel approach that merges k‐means clustering with locally weighted linear regression to denoise transient electromagnetic data. The results from synthetic examples illustrate that the k‐means locally weighted linear regression method can predict transient electromagnetic data closely resembling true values, similar to the long short‐term memory autoencoder. Occam's inversion results derived from denoised data using both the k‐means locally weighted linear regression and long short‐term memory–autoencoder methods can well reflect the true model. Notably, a key advantage of the k‐means locally weighted linear regression method is its independence from labelled data as the sample set. The k‐means locally weighted linear regression method was applied to field data collected at the Narenbaolige coalfield in Inner Mongolia, China. Occam's inversion models generated from the denoised field data delineate the boundary between the basaltic body and sedimentary rocks, aligning with drilling data. The inversion models derived from the noisy field data also can capture this boundary, but deep section views reveal the presence of numerous intricate high‐resistivity anomalous bodies. These observations highlight the effectiveness of the k‐means locally weighted linear regression method in denoising transient electromagnetic data.

Signatures of circumbinary disc dynamics in multimessenger population studies of massive black hole binaries

ABSTRACT We investigate the effect of the cutting-edge circumbinary disc (CBD) evolution models on massive black hole binary (MBHB) populations and the gravitational wave background (GWB). We show that CBD-driven evolution leaves a tell-tale signature in MBHB populations, by driving binaries towards an equilibrium eccentricity that depends on the binary mass ratio. We find high orbital eccentricities ($e_{\rm b} \sim 0.5$) as MBHBs enter multimessenger observable frequency bands. The CBD-induced eccentricity distribution of MBHB populations in observable bands is independent of the initial eccentricity distribution at binary formation, erasing any memory of eccentricities induced in the large-scale dynamics of merging galaxies. Our results suggest that eccentric MBHBs are the rule rather than the exception in upcoming transient surveys, provided that CBDs regularly form in MBHB systems. We show that the GWB amplitude is sensitive to CBD-driven preferential accretion onto the secondary, resulting in an increase in GWB amplitude $A_{\rm yr^{-1}}$ by over 100 per cent with just 10 per cent Eddington accretion. As we self-consistently allow for binary hardening and softening, we show that CBD-driven orbital expansion does not diminish the GWB amplitude, and instead increases the amplitude by a small amount. We further present detection rates and population statistics of MBHBs with $M_{\rm b} \gtrsim 10^6 \, {\rm M}_{\odot }$ in Laser Interferometer Space Antenna, showing that most binaries have equal mass ratios and can retain residual eccentricities up to $e_{\rm b} \sim 10^{-3}$ due to CBD-driven evolution.

Quantitative suppression of unsteady powerline noise in transient electromagnetic surveys: Adjustment of base-frequency and optimal choice of stacking-times.

Powerline noise is a severe interference source in urban or mine transient electromagnetic (TEM) surveys. TEM systems generally adopt the synchronous detection scheme to suppress the powerline noise. However, when the powerline frequency fluctuates differing from its nominal value (50/60Hz), the considerable powerline noise residue still remains even after synchronous detection. To overcome this problem, this paper proposes the quantitative suppression method for powerline noise taking into account the instability of powerline frequency. The method is based on the adjustment of base-frequency and optimal choice of stacking-times. We represent mathematically the sufficient condition for powerline noise suppression by synchronous detection. It consists of two equations, one with respect to base-frequency and other with respect to stacking-times. The base-frequency is adjusted according to powerline frequency estimate. We first derive the mathematical relationship between frequency estimation accuracy, residual noise amplitude, and stacking-times. Based on it, we develop an efficient algorithm to determine the optimal stacking-times. The algorithm takes as an input the powerline noise estimates and the noise tolerance limit. The tolerance limit is set to a certain value desired by the user, such as 10, 4, and 1 µV. By adjusting base-frequency and determining optimal stacking-times, the powerline noise residue after synchronous detection is reduced to below the desired tolerance limit. We verify the effectiveness of the proposed method for both simulated and actual noise. Experimental results show that the method achieves quantitative suppression of unsteady powerline noise without any damage of effective signal and prevents the measurement time loss due to excessive stacking.

Submillimeter and Mid-Infrared Variability of Young Stellar Objects in the M17SWex Intermediate-mass Star-forming Region

We present a comprehensive analysis of young stellar object (YSO) variability within the M17 Southwest Extension (M17 SWex), using 3.5 yr of monitoring data from the James Clerk Maxwell Telescope (JCMT) Transient Survey at submillimeter and 9 yr from the NEOWISE mission at mid-infrared (mid-IR). Our study encompasses observations of 147 bright submillimeter peaks identified within our deep JCMT coadded map as well as 156 YSOs in NEOWISE W1 and 179 in W2 that were previously identified in Spitzer surveys. We find three robust submillimeter variables: two are candidate YSOs and one is a likely extragalactic source. At mid-IR wavelengths, our analysis reveals secular and stochastic variability in 47 YSOs, with the highest fraction of secular variability occurring at the earliest evolutionary stage. This is similar to what has previously been observed for low-mass YSO variability within the Gould Belt. However, we observe less overall variability in M17 SWex at both the submillimeter and mid-IR. We suspect that this lower fraction is due to the greater distance to M17 SWex. Our findings showcase the utility of multiwavelength observations to better capture the complex variability phenomena inherent to star formation processes and demonstrate the importance of years-long monitoring of a diverse selection of star-forming environments.

The Zwicky Transient Facility Bright Transient Survey. III. BTSbot: Automated Identification and Follow-up of Bright Transients with Deep Learning

The Bright Transient Survey (BTS) aims to obtain a classification spectrum for all bright (m peak ≤ 18.5 mag) extragalactic transients found in the Zwicky Transient Facility (ZTF) public survey. BTS critically relies on visual inspection (“scanning”) to select targets for spectroscopic follow-up, which, while effective, has required a significant time investment over the past ∼5 yr of ZTF operations. We present BTSbot, a multimodal convolutional neural network, which provides a bright transient score to individual ZTF detections using their image data and 25 extracted features. BTSbot is able to eliminate the need for daily human scanning by automatically identifying and requesting spectroscopic follow-up observations of new bright transient candidates. BTSbot recovers all bright transients in our test split and performs on par with scanners in terms of identification speed (on average, ∼1 hr quicker than scanners). We also find that BTSbot is not significantly impacted by any data shift by comparing performance across a concealed test split and a sample of very recent BTS candidates. BTSbot has been integrated into Fritz and Kowalski, ZTF’s first-party marshal and alert broker, and now sends automatic spectroscopic follow-up requests for the new transients it identifies. Between 2023 December and 2024 May, BTSbot selected 609 sources in real time, 96% of which were real extragalactic transients. With BTSbot and other automation tools, the BTS workflow has produced the first fully automatic end-to-end discovery and classification of a transient, representing a significant reduction in the human time needed to scan.

Short-period post-common envelope binaries with Balmer emission from SDSS and LAMOST based on ZTF photometric data

ABSTRACT We present here 55 short-period post-common envelope binaries (PCEBs) containing a hot white dwarf (WD) and a low-mass main sequence (MS). Based on the photometric data from Zwicky Transient Facility survey data Release 19 (ZTF DR19), the light curves are analysed for about 200 WDMS binaries with emission line(s) identified from the Sloan Digital Sky Survey (SDSS) or the Large Sky Area Multi-Object Fibre Spectroscopic Telescope (LAMOST) spectra, in which 55 WDMS binaries are found to exhibit variability in their luminosities with a short period and are thus short-period binaries (i.e. PCEBs). In addition, it is found that the orbital periods of these PCEBs locate in a range from 2.2643 to 81.1526 h. However, only six short-period PCEBs are newly discovered and the orbital periods of 19 PCEBs are improved in this work. Meanwhile, it is found that three objects are newly discovered eclipsing PCEBs, and a object (i.e. SDSS J1541) might be the short-period PCEB with a late M-type star or a brown dwarf companion based on the analysis of its spectral energy distribution. At last, the mechanism(s) being responsible for the emission features in the spectra of these PCEBs are discussed, the emission features arising in their optical spectra might be caused by the stellar activity or an irradiated component owing to a hot WD companion because most of them contain a WD with an effective temperature higher than $\sim$10 000 K.

3D deformation velocity field analysis and TEM method to detect the tectonic influence on the land subsidence in Zamora, Mexico

The city of Zamora is situated in a tectonic basin and, according to previous interferometric studies, has been experiencing sinking up to 13 cm/yr (2007-2011). Although the reported subsidence pattern (WNW-ESE) is similar to the orientation of the regional fault system, there is a lack of detailed studies to establish the connection between this phenomenon, the thickness of sediments, and the basement geometry. Therefore, this paper presents a 3D deformation velocity field analysis for the period of 2014–2021, using 130 images from the Sentinel-2 satellite to calculate the vertical and east-west components of the subsidence. Likewise, a campaign of 28 Transient Electromagnetic soundings was carried out around the city to determine the thickness of sediments and bedrock geometry. The results of the interferometric analysis reveal a WNW-ESE sinking pattern with the maximum vertical velocities observed near the Zamora Museum, amounting to 108 mm/yr (2014–2017) and 102 mm/yr (2018–2021). Meanwhile, the horizontal component has a heterogeneous behavior in both periods with rates reaching up to 1 cm/year. The Transient Electromagnetic survey results show the presence of two grabens, a horst, and a series of half-grabens. The thickness of sediments varies between 94 and 314 m, with the latter corresponding to the main graben located beneath the city. The axis of this structure and the largest sediment deposits are linked to the area of maximum sinking and the spatial pattern of subsidence.

SN 2023zaw: An Ultrastripped, Nickel-poor Supernova from a Low-mass Progenitor

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 (<0.005 M ☉) stripped-envelope supernovae, which can be explored with deep and high-cadence transient surveys.

New ACV variables discovered in the Zwicky Transient Facility survey

The manifestation of surface spots on magnetic chemically peculiar (mCP) stars is most commonly explained by the atomic diffusion theory, which requires a calm stellar atmosphere and only moderate rotation. While very successful and well described, this theory still needs to be revised and fine-tuned to the observations. Our study aims to enlarge the sample of known photometrically variable mCP stars (ACV variables) to pave the way for more robust and significant statistical studies. We derive accurate physical parameters for these objects and discuss our results in the framework of the atomic diffusion theory. We studied 1314 candidate ACV variables that were selected from the Zwicky Transient Factory catalogue of periodic variables based on light curve characteristics. We investigated these objects using photometric criteria, a colour-magnitude diagram, and spectroscopic data from the LAMOST and Gaia missions to confirm their status as ACV variables. We present a sample of 1232 new ACV variables, including information on distance from the Sun, mass, fractional age on the main sequence, fraction of the radius between the zero-age and terminal-age main sequence, and the equatorial velocity and its ratio to the critical velocity. Our results confirm that the employed selection process is highly efficient for detecting ACV variables. We have identified 38 stars with $v_ equ $). This challenges current theories that cannot explain the occurrence of such fast-rotating mCP stars.

Optical Variability of Blazars in the Tomo-e Gozen Northern Sky Transient Survey

We studied the optical variability of 241 BL Lacertae (BL Lacs) and 83 flat-spectrum radio quasars (FSRQs) from the 4LAC catalog using data from the Tomo-e Gozen Northern Sky Transient Survey, with ∼50 epochs per blazar on average. We excluded blazars whose optical variability may be underestimated due to the influence of their host galaxy based on their optical luminosity (L O ). FSRQs with γ-ray photon index greater than 2.6 exhibit very low optical variability, and their distribution of standard deviation of repeated photometry is significantly different from that of the other FSRQs (Kolmogorov–Smirnov test p-value equal to 5 × 10−6). Among a sample of blazars at any particular cosmological epoch, those with lower γ-ray luminosity (L γ ) tend to have lower optical variability, and those FSRQs with a γ-ray photon index greater than 2.6 tend to have low L γ . We also measured the structure function of optical variability and found that the amplitude of the structure function for FSRQs is higher than previously measured and higher than that of BL Lacs at multiple time lags. Additionally, the amplitude of the structure function of FSRQs with high γ-ray photon index is significantly lower than that of FSRQs with low γ-ray photon index. The structure function of FSRQs of high γ-ray photon index shows a characteristic timescale of more than 10 days, which may be the variability timescale of the accretion disk. In summary, we infer that the optical component of FSRQs with high γ-ray photon index may be dominated by the accretion disk.

NEural Engine for Discovering Luminous Events (NEEDLE): identifying rare transient candidates in real time from host galaxy images

ABSTRACT Known for their efficiency in analysing large data sets, machine learning-based classifiers have been widely used in wide-field sky survey pipelines. The upcoming Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST) will generate millions of real-time alerts every night, enabling the discovery of large samples of rare events. Identifying such objects soon after explosion will be essential to study their evolution. Using ∼5400 transients from the Zwicky Transient Facility (ZTF) Bright Transient Survey as training and test data, we develop NEEDLE (NEural Engine for Discovering Luminous Events), a novel hybrid (convolutional neural network + dense neural network) classifier to select for two rare classes with strong environmental preferences: superluminous supernovae (SLSNe) preferring dwarf galaxies, and tidal disruption events (TDEs) occurring in the centres of nucleated galaxies. The input data includes (i) cutouts of the detection and reference images, (ii) photometric information contained directly in the alert packets, and (iii) host galaxy magnitudes from Pan-STARRS (Panoramic Survey Telescope and Rapid Response System). Despite having only a few tens of examples of the rare classes, our average (best) completeness on an unseen test set reaches 73 per cent (86 per cent) for SLSNe and 80 per cent (87 per cent) for TDEs. While very encouraging for completeness, this may still result in relatively low purity for the rare transients, given the large class imbalance in real surveys. However, the goal of NEEDLE is to find good candidates for spectroscopic classification, rather than to select pure photometric samples. Our system will be deployed as an annotator on the UK alert broker, Lasair, to provide predictions of real-time alerts from ZTF and LSST to the community.

The JCMT Transient Survey: Six Year Summary of 450/850 μm Protostellar Variability and Calibration Pipeline Version 2.0

The James Clerk Maxwell Telescope (JCMT) Transient Survey has been monitoring eight Gould Belt low-mass star-forming regions since 2015 December and six somewhat more distant intermediate-mass star-forming regions since 2020 February with the Submillimeter Common User Bolometer Array 2 on board JCMT at 450 and 850 μm and with an approximately monthly cadence. We introduce our pipeline v2 relative calibration procedures for image alignment and flux calibration across epochs, improving on our previous pipeline v1 by decreasing measurement uncertainties and providing additional robustness. These new techniques work at both 850 and 450 μm, where version 1 only allowed investigation of the 850 μm data. Pipeline v2 achieves better than 0.″5 relative image alignment, less than a tenth of the submillimeter beam widths. The version 2 relative flux calibration is found to be 1% at 850 μm and <5% at 450 μm. The improvement in the calibration is demonstrated by comparing the two pipelines over the first 4 yr of the survey and recovering additional robust variables with version 2. Using the full 6 yr of the Gould Belt survey, the number of robust variables increases by 50%, and at 450 μm we identify four robust variables, all of which are also robust at 850 μm. The multiwavelength light curves for these sources are investigated and found to be consistent with the variability being due to dust heating within the envelope in response to accretion luminosity changes from the central source.

Dramatic Rebrightening of the Type-changing Stripped-envelope Supernova SN 2023aew

Multipeaked supernovae with precursors, dramatic light-curve rebrightenings, and spectral transformation are rare, but are being discovered in increasing numbers by modern night-sky transient surveys like the Zwicky Transient Facility. Here, we present the observations and analysis of SN 2023aew, which showed a dramatic increase in brightness following an initial luminous (−17.4 mag) and long (∼100 days) unusual first peak (possibly precursor). SN 2023aew was classified as a Type IIb supernova during the first peak but changed its type to resemble a stripped-envelope supernova (SESN) after the marked rebrightening. We present comparisons of SN 2023aew’s spectral evolution with SESN subtypes and argue that it is similar to SNe Ibc during its main peak. P-Cygni Balmer lines are present during the first peak, but vanish during the second peak’s photospheric phase, before Hα resurfaces again during the nebular phase. The nebular lines ([O i], [Ca ii], Mg i], Hα) exhibit a double-peaked structure that hints toward a clumpy or nonspherical ejecta. We analyze the second peak in the light curve of SN 2023aew and find it to be broader than that of normal SESNe as well as requiring a very high 56Ni mass to power the peak luminosity. We discuss the possible origins of SN 2023aew including an eruption scenario where a part of the envelope is ejected during the first peak and also powers the second peak of the light curve through interaction of the SN with the circumstellar medium.

The Next Generation Virgo Cluster Survey. XXXVII. Distant RR Lyrae Stars and the Milky Way Stellar Halo Out to 300 kpc

RR Lyrae stars are standard candles with characteristic photometric variability and serve as powerful tracers of Galactic structure, substructure, accretion history, and dark matter content. Here we report the discovery of distant RR Lyrae stars, including some of the most distant stars known in the Milky Way halo, with Galactocentric distances of ∼300 kpc. We use time-series u*g′i′z′ Canada–France–Hawaii Telescope/MegaCam photometry from the Next Generation Virgo Cluster Survey (NGVS). We use a template light-curve fitting method based on empirical Sloan Digital Sky Survey Stripe 82 RR Lyrae data to identify RR Lyrae candidates in the NGVS data set. We eliminate several hundred suspected quasars and identify 180 RR Lyrae candidates with heliocentric distances of ∼20–300 kpc. The halo stellar density distribution is consistent with an r −4.09±0.10 power-law radial profile over most of this distance range with no signs of a break. The distribution of ab-type RR Lyrae in a period–amplitude plot (Bailey diagram) suggests that the mean metallicity of the halo decreases outward. Compared to other recent RR Lyrae surveys, like Pan-STARRS1, the High Cadence Transient Survey, and the Dark Energy Survey, our NGVS study has better single-epoch photometric precision and a comparable number of epochs but smaller sky coverage. At large distances, our RR Lyrae sample appears to be relatively pure and complete, with well-measured periods and amplitudes. These newly discovered distant RR Lyrae stars are important additions to the few secure stellar tracers beyond 150 kpc in the Milky Way halo.

Taking the pulse of the outer Milky Way with the Halo Outskirts With Variable Stars (HOWVAST) survey: an RR Lyrae density profile out to &amp;gt;200 kpc

ABSTRACT In order to constrain the evolutionary history of the Milky Way, we hunt for faint RR Lyrae stars (RRLs) using Dark Energy Camera data from the High cadence Transient Survey and the Halo Outskirts With Variable Stars survey. We report the detection of $\sim$500 RRLs, including previously identified stars and $\sim$90 RRLs not yet reported. We identify nine new RRLs beyond 100 kpc from the Sun, most of which are classified as fundamental-mode pulsators. The periods and amplitudes of the distant RRLs do not place them in either one of the two classical Oosterhoff groups, but in the Oosterhoff intermediate region. We detect two groups of clumped distant RRLs with similar distances and equatorial coordinates, which we interpret as an indication of their association with undiscovered bound or unbound satellites. We study the halo density profile using spheroidal and ellipsoidal ($q=0.7$) models, following a Markov chain Monte Carlo methodology. For a spheroidal halo, our derived radial profile is consistent with a broken power law with a break at $18.1^{+2.1}_{-1.1}$ kpc separating the inner and the outer halo, and an outer slope of $-4.47^{+0.11}_{-0.18}$. For an ellipsoidal halo, the break is located at $24.3^{+2.6}_{-3.2}$ kpc and the outer slope is $-4.57^{+0.17}_{-0.25}$. The break in the density profile is a feature visible in different directions of the halo. The similarity of these radial distributions with previous values reported in the literature seems to depend on the regions of the sky surveyed (direction and total area) and halo tracer used. Our findings are compatible with simulations and observations that predict that the outer regions of Milky Way-like galaxies are mainly composed of accreted material.

NuSTAR Observations of Candidate Subparsec Binary Supermassive Black Holes

We present an analysis of NuSTAR X-ray observations of three active galactic nuclei (AGN) that were identified as candidate subparsec binary supermassive black hole (SMBH) systems in the Catalina Real-Time Transient Survey based on apparent periodicity in their optical light curves. Simulations predict that close-separation accreting SMBH binaries will have different X-ray spectra than single accreting SMBHs. We previously observed these AGN with Chandra and found no differences between their low-energy X-ray properties and the larger AGN population. However, some models predict differences to be more prominent at energies higher than probed by Chandra. We find that even at the higher energies probed by NuSTAR, the spectra of these AGN are indistinguishable from the larger AGN population. This could rule out models predicting large differences in the X-ray spectra in the NuSTAR bands. Alternatively, it might mean that these three AGN are not binary SMBHs.

Searching for precursor activity of Type IIn Supernovae

We conducted a search for luminous outbursts prior to the explosion of Type IIn supernovae (SNe IIn). We built a sample of 27 objects spectroscopically classified as SNe IIn and all located at $z&lt;0.015$. Using deep archival SN fields images ---taken up to nearly 20 years prior to the SN explosions themselves--- from transient surveys (PTF, ZTF, DES and CHASE) and major astronomical observatories (European Southern Observatory, ESO and National Optical Astronomy Observatory, NOAO), we found at least one outburst years to months before the explosion of seven SNe IIn, with the earliest precursor being 10 years prior to the explosion of SN 2019bxq. The maximum absolute magnitudes of the outbursts range between $-$11.5 mag and $-$15 mag, and the eruptive phases last for a few weeks to a few years. The $g-r$ colour measured for three objects during their outburst is relatively red, with $g-r$ ranging between 0.5 and 1.0 mag. This is similar to the colour expected during the eruptions of luminous blue variables. We note that the light curves of SNe with pre-SN outbursts have faster decline rates than those of the SNe that do not show pre-SN outbursts. SN 2011fh is remarkable, as it is still visible 12 years after the luminous SN-like event, indicating that the progenitor possibly survived, or that the interaction is still ongoing. We detect precursor activity in 29&lt;!PCT!&gt; of bona fide SNe IIn in our sample. However, a quantitative assessment of the observational biases affecting the sample suggests that this fraction is an underestimation of the intrinsic precursor occurrence rate.

The distortion and misinterpretation of TEM responses caused by the IP effect

Abstract Transient electromagnetic surveys are commonly conducted to map the distribution of resistivity, a key physical property in mineral exploration and other geological prospecting problems. However, the responses obtained in regions associated with chargeable minerals are always distorted by the induced polarization effects. In this study, the distorted responses are initially simulated in the frequency domain employing the Cole–Cole complex resistivity model and subsequently converted into the time domain through a time-frequency transformation method. A uniform half-space model is employed to validate the algorithm and illustrate the distortion characteristics of the responses in polarizable formations. A three-layer model is designed to estimate the misinterpretation of slightly complicated models. An actual misinterpretation is demonstrated by field responses containing induced polarization effects collected in the Wulong Gold Mine. The results show that the distortions under different geoelectrical conditions are consistent, enhancing the responses in the early stage and counteracting the responses in the late stage. The strong induced polarizable effects distort the responses by causing explicit sign reversals, whereas the weak induced polarizable effects only distort the decay rate of the responses. These distortions are prone to causing misinterpretations and resulting in excessively intricate geological structures.

A Pan-STARRS1 Search for Planet Nine

We present a search for Planet Nine using the second data release of the Pan-STARRS1 survey. We rule out the existence of a Planet Nine with the characteristics of that predicted in Brown & Batygin to a 50% completion depth of V = 21.5. This survey, along with previous analyses of the Zwicky Transient Facility and Dark Energy Survey data, rules out 78% of the Brown & Batygin parameter space. Much of the remaining parameter space is at V > 21 in regions near and in the area where the northern galactic plane crosses the ecliptic.

Unraveling the dust activity of naked-eye comet C/2022 E3 (ZTF)

A morphological and photometric analysis of the naked-eye long-period comet C/2022 E3 (ZTF) before perihelion is presented in this study. The observation images taken by the Zwicky Transient Facility survey telescope from July 2022 to October 2022 show a gradually brightening dust coma and a tail with a clear structure. The morphology of the dust coma reveals nonsteady-state emission with an ejection velocity lower than 14 m s−1 for particles larger than 100 µm. According to the syndyne-synchrone analysis, dust particles larger than about 10 µm contribute significantly to the observed tail. The model simulations of the 10 October 2022 image suggest that the radii of large particles lingering near the nucleus range from 0.1 to 1 mm. Assuming that the nucleus of comet E3 is a homogeneous sphere with an albedo of 0.1, the photometry analysis sets the lower and upper limits of the nucleus radius to be 0.81 ± 0.07 km and 2.79 ± 0.01 km, respectively. The dust production rates increased continuously from 241 ± 3 kg s−1 in July to 476 ± 9 kg s−1 in October. The dependence of the ejection velocity v⊥ perpendicular to the orbital plane of comet E3 on the particle size a can be simplified as v⊥ ∝ a−1/2, which indicates that the dust emission is likely driven by gas. The water-production rate is inferred as ~368 ± 72 kg s−1 in October 2022, which is sustained by an equilibrium-sublimating area of 8.2 × 106 m2 at least. The comparative analysis of the characteristics of comet E3 with those of comets belonging to different types shows that the activity profile of long-period comet E3 surprisingly aligns more closely with those of short-period comets within a heliocentric distance range of about [1.7, 3.4] AU, where the images of comet E3 that we used in this study were taken.