Search

Accurately estimating the parameters of the nanohertz gravitational-wave background is essential for understanding its origin. The background is typically modeled with a power-law spectrum, parametrized with an amplitude $A$, which describes its intensity, and a spectral index $\gamma$, which describes how the background varies with frequency. Different collaborations have produced varied estimates of $\gamma$, some in tension with the value of $\gamma = 13/3$ expected for circular, gravitational-wave-driven binary black holes. However, estimates of $A$ and $\gamma$ can be affected by systematic errors and misspecified noise models. We investigate how systematic errors, which may plausibly be present in pulsar-timing analyses, can shift inferences about $A, \gamma$. We demonstrate that conservatively incorporating noise sources into the model that are not actually present in the data does not produce bias inferences in practice. This addresses concerns that an overly complex noise model might lead to bias from a needlessly conservative prior. Our results highlight the importance of using comprehensive noise models in pulsar timing analyses to ensure accurate and reliable parameter estimation of the gravitational-wave background., Comment: 10 pages, 5 figures, 1 table

Long-period radio transients are a novel class of astronomical objects characterised by prolonged periods ranging from 18 minutes to 54 minutes. They exhibit highly polarised, coherent, beamed radio emission lasting only 10--100 seconds. The intrinsic nature of these objects is subject to speculation, with highly magnetised white dwarfs and neutron stars being the prevailing candidates. Here we present ASKAP J183950.5-075635.0 (hereafter, ASKAP J1839-0756), boasting the longest known period of this class at 6.45 hours. It exhibits emission characteristics of an ordered dipolar magnetic field, with pulsar-like bright main pulses and weaker interpulses offset by about half a period are indicative of an oblique or orthogonal rotator. This phenomenon, observed for the first time in a long-period radio transient, confirms that the radio emission originates from both magnetic poles and that the observed period corresponds to the rotation period. The spectroscopic and polarimetric properties of ASKAP J1839-0756 are consistent with a neutron star origin, and this object is a crucial piece of evidence in our understanding of long-period radio sources and their links to neutron stars., Comment: 44 pages, 14 figures, 2 tables This preprint has not undergone peer review or any post-submission improvements or corrections. The Version of Record of this article is published in Nature Astronomy, and is available online at https://doi.org/10.1038/s41550-024-02452-z

We study an inflation model with nonminimal derivative coupling that features a coupling between the derivative of the inflaton field and the Einstein tensor. This model naturally amplifies curvature perturbations at small scales via gravitationally enhanced friction, a mechanism critical for the formation of primordial black holes and the associated production of potentially detectable scalar-induced gravitational waves. We derive analytical expressions for the primordial power spectrum, enabling efficient exploration of the model parameter space without requiring computationally intensive numerical solutions of the Mukhanov-Sasaki equation. Using the third data release of the Parkes Pulsar Timing Array (PPTA DR3), we constrain the model parameters characterizing the coupling function: $\phi_c = 3.7^{+0.3}_{-0.5} M_\mathrm{P}$, $\log_{10} \omega_L = 7.1^{+0.6}_{-0.3}$, and $\log_{10} \sigma = -8.3^{+0.3}_{-0.6}$ at 90\% confidence level. Our results demonstrate the growing capability of pulsar timing arrays to probe early Universe physics, complementing traditional cosmic microwave background observations by providing unique constraints on inflationary dynamics at small scales., Comment: 14 pages, 5figures

We conduct the first-ever Pulsar Polarization Array (PPA) analysis to detect the ultralight Axion-Like Dark Matter (ALDM) using the polarization data of 22 millisecond pulsars from the third data release of Parkes Pulsar Timing Array. As one of the major dark matter candidates, the ultralight ALDM exhibits a pronounced wave nature on astronomical scales and offers a promising solution to small-scale structure issues within local galaxies. While the linearly polarized pulsar light travels through the ALDM galactic halo, its position angle (PA) can be subject to an oscillation induced by the ALDM Chern-Simons coupling with electromagnetic field. The PPA is thus especially suited for detecting the ultralight ALDM by correlating polarization data across the arrayed pulsars. To accomplish this task, we develop an advanced Bayesian analysis framework that allows us to construct pulsar PA residual time series, model noise contributions properly and search for pulsar cross-correlations. We find that for an ALDM density of $\rho_0=0.4\,\textrm{GeV}/\textrm{cm}^3$, the Parkes PPA offers the best global limits on the ALDM Chern-Simons coupling, namely $\lesssim 10^{-13.5}-10^{-12.2}~{\rm GeV}^{-1}$, for the mass range of $10^{-22} - 10^{-21}~{\rm eV}$. The crucial role of pulsar cross-correlation in recognizing the nature of the derived limits is also highlighted., Comment: 6+15 pages, 10 figures, 2 tables, submitted to the journal

Recently, a class of long-period radio transients (LPTs) has been discovered, exhibiting emission on timescales thousands of times longer than radio pulsars. Several models had been proposed implicating either a strong magnetic field neutron star, isolated white dwarf pulsar, or a white dwarf binary system with a low-mass companion. While several models for LPTs also predict X-ray emission, no LPTs have been detected in X-rays despite extensive searches. Here we report the discovery of an extremely bright LPT (10-20 Jy in radio), ASKAP J1832-0911, which has coincident radio and X-ray emission, both with a 44.2-minute period. The X-ray and radio luminosities are correlated and vary by several orders of magnitude. These properties are unique amongst known Galactic objects and require a new explanation. We consider a $\gtrsim0.5$ Myr old magnetar with a $\gtrsim 10^{13}$ G crustal field, or an extremely magnetised white dwarf in a binary system with a dwarf companion, to be plausible explanations for ASKAP J1832-0911, although both explanations pose significant challenges to formation and emission theories. The X-ray detection also establishes a new class of hour-scale periodic X-ray transients of luminosity $\sim10^{33}$ erg/s associated with exceptionally bright coherent radio emission., Comment: 52 pages, 17 figures, 3 tables

The adiabatic elastocaloric effect relates changes in the strain that a material experiences to resulting changes in its temperature. While elastocaloric materials have been utilized for cooling in room temperature applications, the use of such materials for cryogenic cooling remains relatively unexplored. Here, we use a strain load/unload technique at low temperatures, similar to those employed at room-temperature, to demonstrate a large cooling effect in TmVO$_4$. For strain changes of $1.8 \cdot 10^{-3}$, the inferred cooling reaches approximately 50% of the material's starting temperature at 5 K, justifying the moniker "giant". Beyond establishing the suitability of this class of material for cryogenic elastocaloric cooling, these measurements also provide additional insight to the entropy landscape in the material as a function of strain and temperature, including the behavior proximate to the quadrupolar phase transition.

This paper demonstrates the impact of state-of-the-art instrumental calibration techniques on the precision of arrival times obtained from 9.6 years of observations of millisecond pulsars using the Murriyang 64-m CSIRO Parkes Radio Telescope. Our study focuses on 21-cm observations of 25 high-priority pulsars that are regularly observed as part of the Parkes Pulsar Timing Array (PPTA) project, including those predicted to be the most susceptible to calibration errors. We employ Measurement Equation Template Matching (METM) for instrumental calibration and Matrix Template Matching (MTM) for arrival time estimation, resulting in significantly improved timing residuals with up to a sixfold reduction in white noise compared to arrival times estimated using Scalar Template Matching and conventional calibration based on the Ideal Feed Assumption. The median relative reduction in white noise is 33 percent, and the maximum absolute reduction is 4.5 microseconds. For PSR J0437-4715, METM and MTM reduce the best-fit power-law amplitude (2.7 sigma) and spectral index (1.7 sigma) of the red noise in the arrival time residuals, which can can be tentatively interpreted as mitigation of 1/f noise due to otherwise unmodeled steps in polarimetric response. These findings demonstrate the potential to directly enhance the sensitivity of pulsar timing array experiments through more accurate methods of instrumental calibration and arrival time estimation., Comment: 30 pages, 9 figures, accepted for publication in The Astrophysical Journal

The observation of neutron stars enables the otherwise impossible study of fundamental physical processes. The timing of binary radio pulsars is particularly powerful, as it enables precise characterization of their (three-dimensional) positions and orbits. PSR~J0437$-$4715 is an important millisecond pulsar for timing array experiments and is also a primary target for the Neutron Star Interior Composition Explorer (NICER). The main aim of the NICER mission is to constrain the neutron star equation of state by inferring the compactness ($M_p/R$) of the star. Direct measurements of the mass $M_p$ from pulsar timing therefore substantially improve constraints on the radius $R$ and the equation of state. Here we use observations spanning 26 years from Murriyang, the 64-m Parkes radio telescope, to improve the timing model for this pulsar. Among the new precise measurements are the pulsar mass $M_p=1.418\pm 0.044$ $M_{\odot}$, distance $D=156.96 \pm 0.11$ pc, and orbital inclination angle $i=137.506 \pm 0.016^\circ$, which can be used to inform the X-ray pulse profile models inferred from NICER observations. We demonstrate that these results are consistent between multiple data sets from the Parkes Pulsar Timing Array (PPTA), each modeled with different noise assumptions. Using the longest available PPTA data set, we measure an apparent second derivative of the pulsar spin frequency and discuss how this can be explained either by kinematic effects due to the proper motion and radial velocity of the pulsar or excess low-frequency noise such as a gravitational-wave background., Comment: 13 pages, 3 figures, accepted for publication in Astrophysical Journal Letters

We have searched for anomalous events using 2,520 hours of archival observations from Murriyang, CSIRO's Parkes radio telescope. These observations were originally undertaken to search for pulsars. We used a machine-learning algorithm based on ResNet and Uniform Manifold Approximation and Projection (UMAP) in order to identify parts of the data stream that potentially contain anomalous signals. Many of these anomalous events are radio frequency interference, which were subsequently filtered using multibeam information. We detected 202 anomalous events and provide their positions and event times. Our results show that the UMAP unsupervised machine learning pipeline effectively identifies anomalous signals in high-time-resolution datasets, highlighting its potential for use in future surveys. However, the pipeline is not applicable for standard searches for dispersed single pulses. We classify the detected events and, in particular, we are currently unable to determine the possible origin of events that last multiple seconds. For these we encourage follow-up observations., Comment: 11 pages, 10 figures; accepted by MNRAS

Several sources of repeating coherent bursts of radio emission with periods of many minutes have now been reported in the literature. These "ultra-long period" (ULP) sources have no clear multi-wavelength counterparts and challenge canonical pulsar emission models, leading to debate regarding their nature. In this work we report the discovery of a bright, highly-polarised burst of radio emission at low Galactic latitude as part of a wide-field survey for transient and variable radio sources. ASKAP\,J175534.9$-$252749.1 does not appear to repeat, with only a single intense two-minute $\sim$200-mJy burst detected from 60~hours of observations. The burst morphology and polarisation properties are comparable to those of classical pulsars but the duration is more than one hundred times longer, analogous to ULPs. Combined with the existing ULP population, this suggests that these sources have a strong Galactic latitude dependence and hints at an unexplored population of transient and variable radio sources in the thin disk of the Milky Way. The resemblance of this burst with both ULPs and pulsars calls for a unified coherent emission model for objects with spin periods from milliseconds to tens of minutes. However, whether or not these are all neutron stars or have the same underlying power source remains open for debate.

Fast Radio Bursts (FRBs) are short-duration radio transients that occur at random times in host galaxies distributed all over the sky. Large field of view instruments can play a critical role in the blind search for rare FRBs. We present a concept for an all-sky FRB monitor using a compact all-sky phased array (CASPA), which can efficiently achieve an extremely large field of view of $\sim10^4$ square degrees. Such a system would allow us to conduct a continuous, blind FRB search covering the entire southern sky. Using the measured FRB luminosity function, we investigate the detection rate for this all-sky phased array and compare the result to a number of other proposed large field-of-view instruments. We predict a rate of a few FRB detections per week and determine the dispersion measure and redshift distributions of these detectable FRBs. This instrument is optimal for detecting FRBs in the nearby Universe and for extending the high-end of the FRB luminosity function through finding ultraluminous events. Additionally, this instrument can be used to shadow the new gravitational-wave observing runs, detect high energy events triggered from Galactic magnetars and search for other bright, but currently unknown transient signals., Comment: 10 pages, 5 figures, accepted for publication in PASA

We present the Sydney Radio Star Catalogue, a new catalogue of stars detected at megahertz to gigahertz radio frequencies. It consists of 839 unique stars with 3,405 radio detections, more than doubling the previously known number of radio stars. We have included stars from large area searches for radio stars found using circular polarisation searches, cross-matching, variability searches, and proper motion searches as well as presenting hundreds of newly detected stars from our search of Australian SKA Pathfinder observations. The focus of this first version of the catalogue is on objects detected in surveys using SKA precursor instruments; however we will expand this scope in future versions. The 839 objects in the Sydney Radio Star Catalogue are distributed across the whole sky and range from ultracool dwarfs to Wolf-Rayet stars. We find that the radio luminosities of cool dwarfs are lower than the radio luminosities of more evolved sub-giant and giant stars. We use X-ray detections of 530 radio stars by the eROSITA soft X-ray instrument onboard the SRG spacecraft to show that almost all of the radio stars in the catalogue are over-luminous in the radio, indicating that the majority of stars at these radio frequencies are coherent radio emitters. The Sydney Radio Star Catalogue can be found in Vizier or at https://radiostars.org., Comment: 21 pages, 8 Figures, 7 Tables, Submitted to PASA

A number of pulsar timing arrays have recently reported preliminary evidence for the existence of a nanohertz frequency gravitational-wave background. These analyses rely on detailed noise analyses, which are inherently complex due to the many astrophysical and instrumental factors that contribute to the pulsar noise budget. We investigate whether realistic systematic errors, stemming from misspecified noise models that fail to capture salient features of the pulsar timing noise, could bias the evidence for gravitational waves. We consider two plausible forms of misspecification: small unmodeled jumps and unmodeled chromatic noise. Using simulated data, we calculate the distribution of the commonly used optimal statistic with no signal present and using plausibly misspecified noise models. By comparing the optimal statistic distribution with the distribution created using ``quasi-resampling'' techniques (such as sky scrambles and phase shifts), we endeavor to determine the extent to which plausible misspecification might lead to a false positive. The results are reassuring: we find that quasi-resampling techniques tend to underestimate the significance of pure-noise datasets. We conclude that recent reported evidence for a nanohertz gravitational-wave background is likely robust to the most obvious sources of systematic errors; if anything, the significance of the signal is potentially underestimated., Comment: 8 pages, 10 figures

The suppression of ferroquadrupolar order in TmVO$_4$ in a magnetic field is well-described by the transverse field Ising model, enabling detailed studies of critical dynamics near the quantum phase transition. We describe nuclear magnetic resonance measurements in pure and Y-doped single crystals. The non-Kramers nature of the ground state doublet leads to a unique form of the hyperfine coupling that exclusively probes the transverse field susceptibility. Our results show that this quantity diverges at the critical field, in contrast to the mean-field prediction. Furthermore, we find evidence for quantum critical fluctuations present near Tm-rich regions in Y-doped crystals at levels beyond which long-range order is suppressed, suggesting the presence of quantum Griffiths phases., Comment: 13 pages, 6 figures

Radio astronomy file formats are now required to store wide frequency bandwidths and multiple simultaneous receiver beams and must be able to account for versatile observing modes and numerous calibration strategies. The need to capture and archive high-time and high frequency-resolution data, along with the comprehensive metadata that fully describe the data, implies that a new data format and new processing software are required. This requirement is suited to a well-defined, hierarchically-structured and flexible file format. In this paper we present the Spectral-Domain Hierarchical Data Format (`SDHDF') -- a new file format for radio astronomy data, in particular for single dish or beam-formed data streams. Since 2018, SDHDF has been the primary format for data products from the spectral-line and continuum observing modes at Murriyang, the CSIRO Parkes 64-m radio telescope, and we demonstrate that this data format can also be used to store observations of pulsars and fast radio bursts., Comment: Supplementary material (SDHDF definition): https://ars.els-cdn.com/content/image/1-s2.0-S2213133724000192-mmc1.pdf

We report an investigation of the effect of substitution of Y for Tm in $Tm_{1-x}Y_xVO4$ via low-temperature heat capacity measurements, with the yttrium content $x$ varying from $0$ to $0.997$. Because the Tm ions support a local quadrupolar (nematic) moment, they act as reporters of the local strain state in the material, with the splitting of the ion's non-Kramers crystal field groundstate proportional to the quadrature sum of the in-plane tetragonal symmetry-breaking transverse and longitudinal strains experienced by each ion individually. Analysis of the heat capacity therefore provides detailed insights into the distribution of local strains that arise as a consequence of the chemical substitution. These local strains suppress long-range quadrupole order for $x>0.22$, and result in a broad Schottky-like feature for higher concentrations. Heat capacity data are compared to expectations for a distribution of uncorrelated (random) strains. For dilute Tm concentrations, the heat capacity cannot be accounted for by randomly distributed strains, demonstrating the presence of significant strain correlations between sites. For intermediate Tm concentrations, these correlations must still exist, but the data cannot be distinguished from that which would be obtained from a 2D Gaussian distribution. The cross-over between these limits is discussed in terms of the interplay of key length scales in the substituted material. The central result of this work, that local strains arising from chemical substitution are not uncorrelated, has implications for the range of validity of theoretical models based on random effective fields that are used to describe such chemically substituted materials, particularly when electronic nematic correlations are present.

We study thermalization and thermal transport in single crystals of CsV$_3$Sb$_5$ through the CDW transition by directly measuring thermal diffusivity ($D$), thermal conductivity ($\kappa$), resistivity ($\rho$), and specific heat ($c$). Commensurate with previous reports, we observe a sharp, narrow anomaly in specific heat associated with a first order transition that results in a CDW state below $\sim94$ K. While a corresponding sharp anomaly in thermal diffusivity is also observed, resistivity and thermal conductivity only exhibit small steps at the transition, where the feature is sharp for resistivity and broader for thermal conductivity. Scrutinizing the thermal Einstein relation $\kappa=cD$, we find that this relation is satisfied in the entire temperature range, except in a narrow range around the transition. The Wiedemann-Franz law seems to work outside the critical region as well. Below the transition and persisting below the two-phase regime we find strong resemblance between the resistivity anomaly and the specific heat, which may point to a secondary electronic order parameter that emerges continuously below the transition., Comment: 11 pages, 6 figures

Large widefield surveys make possible the serendipitous discovery of rare sub-classes of pulsars. One such class are "spider"-type pulsar binaries, comprised of a pulsar in a compact orbit with a low-mass (sub)stellar companion. In a search for circularly-polarized radio sources in ASKAP Pilot Survey observations, we discovered highly variable and circularly polarized emission from a radio source within the error region of the $\gamma$-ray source {4FGL}~J1646.5$-$4406. The variability is consistent with the eclipse of a compact, steep-spectrum source behind ablated material from a companion in a $\sim 5.3\,$h binary orbit. Based on the eclipse properties and spatial coincidence with {4FGL} J1646.5$-$4406, we argue that the source is likely a recycled pulsar in a "redback" binary system. Using properties of the eclipses from ASKAP and Murchison Widefield Array observations, we provide broad constraints on the properties of the eclipse medium. We identified a potential optical/infra-red counterpart in archival data consistent with a variable low-mass star. Using the Parkes Radio Telescope "Murriyang" and MeerKAT, we searched extensively for radio pulsations but yielded no viable detections of pulsed emission. We suggest that the non-detection of pulses is due to scattering in the intra-binary material, but scattering from the ISM can also plausibly explain the pulse non-detections if the interstellar dispersion measure exceeds $\sim$600$\,$pc$\,$cm$^{-3}$. Orbital constraints derived from optical observations of the counterpart would be highly valuable for future $\gamma$-ray pulsation searches, which may confirm the source nature as a pulsar., Comment: Submitted to MNRAS. 13 Pages, 10 figures, 3 tables

In a growing number of materials, phonons have been found to generate a thermal Hall effect, but the underlying mechanism remains unclear. Inspired by previous studies that revealed the importance of Tb$^{3+}$ ions in generating the thermal Hall effect in a family of pyrochlores, we investigated the role of Tm$^{3+}$ ions in TmVO$_{4}$, a paramagnetic insulator with a different crystal structure. We observe a negative thermal Hall conductivity in TmVO$_{4}$ with a magnitude such that the Hall angle, $|\kappa_{xy}$/$\kappa_{xx}|$, is approximately 1 x 10$^{-3}$ at $H$ = 15 T and $T$ = 20 K, typical for a phonon-generated thermal Hall effect. In contrast to the negligible $\kappa_{xy}$ found in the nonmagnetic pyrochlore analog (where the Tb$^{3+}$ ions are replaced with Y$^{3+}$), we observe a negative $\kappa_{xy}$ in YVO$_{4}$ with a Hall angle of magnitude comparable to that of TmVO$_{4}$. This shows that the Tm$^{3+}$ ions are not essential for the thermal Hall effect in this family of materials. Interestingly, at an intermediate Y concentration of $x$ = 0.3 in Tm$_{1-x}$Y$_{x}$VO$_{4}$, $\kappa_{xy}$ was found to have a positive sign, pointing to the importance of impurities in the thermal Hall effect of phonons., Comment: Figures modified and added, authors added, appendix added

The Australian, Chinese, European, Indian, and North American pulsar timing array (PTA) collaborations recently reported, at varying levels, evidence for the presence of a nanohertz gravitational wave background (GWB). Given that each PTA made different choices in modeling their data, we perform a comparison of the GWB and individual pulsar noise parameters across the results reported from the PTAs that constitute the International Pulsar Timing Array (IPTA). We show that despite making different modeling choices, there is no significant difference in the GWB parameters that are measured by the different PTAs, agreeing within $1\sigma$. The pulsar noise parameters are also consistent between different PTAs for the majority of the pulsars included in these analyses. We bridge the differences in modeling choices by adopting a standardized noise model for all pulsars and PTAs, finding that under this model there is a reduction in the tension in the pulsar noise parameters. As part of this reanalysis, we "extended" each PTA's data set by adding extra pulsars that were not timed by that PTA. Under these extensions, we find better constraints on the GWB amplitude and a higher signal-to-noise ratio for the Hellings and Downs correlations. These extensions serve as a prelude to the benefits offered by a full combination of data across all pulsars in the IPTA, i.e., the IPTA's Data Release 3, which will involve not just adding in additional pulsars, but also including data from all three PTAs where any given pulsar is timed by more than as single PTA., Comment: 21 pages, 9 figures, submitted to ApJ

Adiabatic decompression of para-quadrupolar materials has significant potential as a cryogenic cooling technology. We focus on TmVO$_4$, an archetypal material that undergoes a continuous phase transition to a ferroquadrupole-ordered state at 2.15 K. Above the phase transition, each Tm ion contributes an entropy of $k_B \ln{2}$ due to the degeneracy of the crystal electric field groundstate. Owing to the large magnetoelastic coupling, which is a prerequisite for a material to undergo a phase transition via the cooperative Jahn-Teller effect, this level splitting, and hence the entropy, can be readily tuned by externally-induced strain. Using a dynamic technique in which the strain is rapidly oscillated, we measure the adiabatic elastocaloric coefficient of single-crystal TmVO$_4$, and thus experimentally obtain the entropy landscape as a function of strain and temperature. The measurement confirms the suitability of this class of materials for cryogenic cooling applications, and provides insight to the dynamic quadrupole strain susceptibility.

The noise in millisecond pulsar (MSP) timing data can include contributions from observing instruments, the interstellar medium, the solar wind, solar system ephemeris errors, and the pulsars themselves. The noise environment must be accurately characterized in order to form the null hypothesis from which signal models can be compared, including the signature induced by nanohertz-frequency gravitational waves (GWs). Here we describe the noise models developed for each of the MSPs in the Parkes Pulsar Timing Array (PPTA) third data release, which have been used as the basis of a search for the isotropic stochastic GW background. We model pulsar spin noise, dispersion measure variations, scattering variations, events in the pulsar magnetospheres, solar wind variability, and instrumental effects. We also search for new timing model parameters and detected Shapiro delays in PSR~J0614$-$3329 and PSR~J1902$-$5105. The noise and timing models are validated by testing the normalized and whitened timing residuals for Gaussianity and residual correlations with time. We demonstrate that the choice of noise models significantly affects the inferred properties of a common-spectrum process. Using our detailed models, the recovered common-spectrum noise in the PPTA is consistent with a power law with a spectral index of $\gamma=13/3$, the value predicted for a stochastic GW background from a population of supermassive black hole binaries driven solely by GW emission., Comment: 18 pages, 10 figures. Accepted for publication in ApJL

We present the third data release from the Parkes Pulsar Timing Array (PPTA) project. The release contains observations of 32 pulsars obtained using the 64-m Parkes "Murriyang" radio telescope. The data span is up to 18 years with a typical cadence of 3 weeks. This data release is formed by combining an updated version of our second data release with $\sim 3$ years of more recent data primarily obtained using an ultra-wide-bandwidth receiver system that operates between 704 and 4032 MHz. We provide calibrated pulse profiles, flux-density dynamic spectra, pulse times of arrival, and initial pulsar timing models. We describe methods for processing such wide-bandwidth observations, and compare this data release with our previous release., Comment: 15 pages, 6 figures. Accepted for publication in PASA

Pulsar timing arrays aim to detect nanohertz-frequency gravitational waves (GWs). A background of GWs modulates pulsar arrival times and manifests as a stochastic process, common to all pulsars, with a signature spatial correlation. Here we describe a search for an isotropic stochastic gravitational-wave background (GWB) using observations of 30 millisecond pulsars from the third data release of the Parkes Pulsar Timing Array (PPTA), which spans 18 years. Using current Bayesian inference techniques we recover and characterize a common-spectrum noise process. Represented as a strain spectrum $h_c = A(f/1 {\rm yr}^{-1})^{\alpha}$, we measure $A=3.1^{+1.3}_{-0.9} \times 10^{-15}$ and $\alpha=-0.45 \pm 0.20$ respectively (median and 68% credible interval). For a spectral index of $\alpha=-2/3$, corresponding to an isotropic background of GWs radiated by inspiraling supermassive black hole binaries, we recover an amplitude of $A=2.04^{+0.25}_{-0.22} \times 10^{-15}$. However, we demonstrate that the apparent signal strength is time-dependent, as the first half of our data set can be used to place an upper limit on $A$ that is in tension with the inferred common-spectrum amplitude using the complete data set. We search for spatial correlations in the observations by hierarchically analyzing individual pulsar pairs, which also allows for significance validation through randomizing pulsar positions on the sky. For a process with $\alpha=-2/3$, we measure spatial correlations consistent with a GWB, with an estimated false-alarm probability of $p \lesssim 0.02$ (approx. $2\sigma$). The long timing baselines of the PPTA and the access to southern pulsars will continue to play an important role in the International Pulsar Timing Array., Comment: 19 pages, 10 figures, Accepted for publication in ApJL

The screening response of vortices in kagome superconductor CsV$_3$Sb$_5$ was measured using the ac mutual inductance technique. Besides confirming the absence of gapless quasiparticles in zero external magnetic field, we observe the peak effect, manifested in enhanced vortex pinning strength and critical current, in a broad intermediate range of magnetic field. The peaks are followed by another crossover from strong to weak pinning, unlike the usual peak effect that diminishes smoothly at $H_{c2}$. Hysteresis in the screening response allows the identification of a vortex glass phase which strongly correlates with the onset of the peaks. A variety of features in the temperature- and field-dependence of the screening response, corroborated by resistance and dc magnetization measurements, have allowed us to extract an $H$-$T$ phase diagram of the vortex states and to infer the irreversibility line $H_\text{irr}(T)$., Comment: Main text: 6 pages, 4 figures; supplementary: 10 pages, 8 figures

Using spin-echo Nuclear Magnetic Resonance in the model Transverse-Field Ising system TmVO$_4$, we show that low frequency quantum fluctuations at the quantum critical point have a very different effect on $^{51}$V nuclear-spins than classical low-frequency noise or fluctuations that arise at a finite temperature critical point. Spin-echos filter out the low frequency classical noise but not the quantum fluctuations. This allows us to directly visualize the quantum critical fan and demonstrate the persistence of quantum fluctuations at the critical coupling strength in TmVO$_4$ to high temperatures in an experiment that remains transparent to finite temperature classical phase transitions. These results show that while dynamical decoupling schemes can be quite effective in eliminating classical noise in a qubit, a quantum critical environment may lead to rapid entanglement and decoherence., Comment: 10 pages, 5 figures

We present results from a radio survey for variable and transient sources on 15-min timescales, using the Australian SKA Pathfinder (ASKAP) pilot surveys. The pilot surveys consist of 505 h of observations conducted at around 1 GHz observing frequency, with a total sky coverage of 1476 deg$^2$. Each observation was tracked for approximately 8-10h, with a typical rms sensitivity of $\sim$30 $\mu$jy/beam and an angular resolution of $\sim$12 arcsec. The variability search was conducted within each 8-10h observation on a 15-min timescale. We detected 38 variable and transient sources. Seven of them are known pulsars, including an eclipsing millisecond pulsar, PSR J2039$-$5617. Another eight sources are stars, only one of which has been previously identified as a radio star. For the remaining 23 objects, 22 are associated with active galactic nuclei or galaxies (including the five intra-hour variables that have been reported previously), and their variations are caused by discrete, local plasma screens. The remaining source has no multi-wavelength counterparts and is therefore yet to be identified. This is the first large-scale radio survey for variables and transient sources on minute timescales at a sub-mJy sensitivity level. We expect to discover $\sim$1 highly variable source per day using the same technique on the full ASKAP surveys., Comment: 20 pages, 11 figures; accepted for publication in MNRAS

The recent observation of a common red-noise process in pulsar timing arrays (PTAs) suggests that the detection of nanohertz gravitational waves might be around the corner. However, in order to confidently attribute this red process to gravitational waves, one must observe the Hellings-Downs curve -- the telltale angular correlation function associated with a gravitational-wave background. This effort is complicated by the complex modelling of pulsar noise. Without proper care, mis-specified noise models can lead to false-positive detections. Background estimation using bootstrap methods such as sky scrambles and phase shifts, which use the data to characterize the noise, are therefore important tools for assessing significance. We investigate the ability of current PTA experiments to estimate their background with "quasi-independent" scrambles -- characterized by a statistical "match" below the fiducial value: $|M|<0.1$. We show that sky scrambling is affected by "saturation" after $O(10)$ quasi-independent realizations; subsequent scrambles are no longer quasi-independent. We show phase scrambling saturates after $O(100)$ quasi-independent realizations. With so few independent scrambles, it is difficult to make reliable statements about the $\gtrsim 5 \sigma$ tail of the null distribution of the detection statistic. We discuss various methods by which one may increase the number of independent scrambles. We also consider an alternative approach wherein one re-frames the background estimation problem so that the significance is calculated using statistically dependent scrambles. The resulting $p$-value is in principle well-defined but may be susceptible to failure if assumptions about the data are incorrect., Comment: Paper accepted for publication on ApJ on the 7th of August 2023

The problem of prediction of behavior of dynamical systems has undergone a paradigm shift in the second half of the 20th century with the discovery of the possibility of chaotic dynamics in simple, physical, dynamical systems for which the laws of evolution do not change in time. The essence of the paradigm is the long term exponential divergence of trajectories. However, that paradigm does not account for another type of unpredictability: the ``Black Swan" event. It also does not account for the fact that short-term prediction is often possible even in systems with exponential divergence. In our framework, the Black Swan type dynamics occurs when an underlying dynamical system suddenly shifts between dynamics of different types. A learning and prediction system should be capable of recognizing the shift in behavior, exemplified by ``confidence loss". In this paradigm, the predictive power is assessed dynamically and confidence level is used to switch between long term prediction and local-in-time prediction. Here we explore the problem of prediction in systems that exhibit such behavior. The mathematical underpinnings of our theory and algorithms are based on an operator-theoretic approach in which the dynamics of the system are embedded into an infinite-dimensional space. We apply the algorithm to a number of case studies including prediction of influenza cases and the COVID-19 pandemic. The results show that the predictive algorithm is robust to perturbations of the available data, induced for example by delays in reporting or sudden increase in cases due to increase in testing capability. This is achieved in an entirely data-driven fashion, with no underlying mathematical model of the disease.

In 2018 an ultra-wide-bandwidth low-frequency (UWL) receiver was installed on the 64-m Parkes Radio Telescope enabling observations with an instantaneous frequency coverage from 704 to 4032 MHz. Here, we present the analysis of a three-year data set of 35 millisecond pulsars observed with the UWL by the Parkes Pulsar Timing Array (PPTA), using wideband timing methods. The two key differences compared to typical narrow-band methods are, firstly, generation of two-dimensional templates accounting for pulse shape evolution with frequency and, secondly, simultaneous measurements of the pulse time-of-arrival (ToA) and dispersion measure (DM). This is the first time that wideband timing has been applied to a uniform data set collected with a single large-fractional bandwidth receiver, for which such techniques were originally developed. As a result of our study, we present a set of profile evolution models and new timing solutions including initial noise analysis. Precision of our ToA and DM measurements is in the range of 0.005 $-$ 2.08 $\mu$s and (0.043$-$14.24)$\times10^{-4}$ cm$^{-3}$ pc, respectively, with 94% of the pulsars achieving a median ToA uncertainty of less than 1 $\mu$s., Comment: 29 pages, 12 figures, Accepted for publication in ApJ

The Deeper, Wider, Faster (DWF) program coordinates observations with telescopes across the electromagnetic spectrum, searching for transients on timescales of milliseconds to days. The tenth DWF observing run was carried out in near real-time during September 2021 and consisted of six consecutive days of observations of the NGC 6744 galaxy group and a field containing the repeating fast radio burst FRB190711 with the Australian Square Kilometre Array Pathfinder, the Dark Energy Camera, the Hard X-ray Modulation Telescope and the Parkes 64m "Murriyang" radio telescope. In this work we present the results of an image-domain search for transient, variable and circularly polarised sources carried out with ASKAP using data from the observing run, along with test observations prior to the run and follow-up observations carried out during and after the run. We identified eight variable radio sources, consisting of one pulsar, six stellar systems (five of which exhibit circularly polarised emission) and one previously uncatalogued source. Of particular interest is the detection of pulses from the ultra-cool dwarf SCR J1845-6357 with a period of $14.2\pm 0.3$ h, in good agreement with the known optical rotation period, making this the slowest rotating radio-loud ultra-cool dwarf discovered.

We present millimeter, optical, and soft X-ray observations of a stellar flare with an energy squarely in the regime of typical X1 solar flares. The flare was observed from Proxima Cen on 2019 May 6 as part of a larger multi-wavelength flare monitoring campaign and was captured by Chandra, LCOGT, du Pont, and ALMA. Millimeter emission appears to be a common occurrence in small stellar flares that had gone undetected until recently, making it difficult to interpret these events within the current multi-wavelength picture of the flaring process. The May 6 event is the smallest stellar millimeter flare detected to date. We compare the relationship between the soft X-ray and millimeter emission to that observed in solar flares. The X-ray and optical flare energies of 10$^{30.3\pm0.2}$ and 10$^{28.9\pm0.1}$ erg, respectively, the coronal temperature of T=11.0$\pm$2.1 MK, and the emission measure of 9.5$\pm$2.2 X 10$^{49}$ cm$^{-3}$ are consistent with M-X class solar flares. We find the soft X-ray and millimeter emission during quiescence are consistent with the Gudel-Benz Relation, but not during the flare. The millimeter luminosity is >100X higher than that of an equivalent X1 solar flare and lasts only seconds instead of minutes as seen for solar flares., Comment: 12 pages, 5 figures. Accepted to ApJ

New classes of astronomical objects are often discovered serendipitously. The enormous data volumes produced by recent high-time resolution, radio-telescope surveys imply that efficient algorithms are required for a discovery. Such algorithms are usually tuned to detect specific, known sources. Existing data sets therefore likely contain unknown astronomical sources, which will remain undetected unless algorithms are developed that can detect a more diverse range of signals. We present the Single-dish PARKES data challenge for finding the uneXpected (SPARKESX), a compilation of real and simulated high-time resolution observations. SPARKESX comprises three mock surveys from the Parkes "Murriyang" radio telescope. A broad selection of simulated and injected expected signals (such as pulsars, fast radio bursts), poorly characterised signals (plausible flare star signatures) and unknown unknowns are generated for each survey. The goal of this challenge is to aid in the development of new algorithms that can detect a wide-range of source types. We show how successful a typical pipeline based on the standard pulsar search software, PRESTO, is at finding the injected signals. The dataset is publicly available at https://doi.org/10.25919/fd4f-0g20., Comment: Accepted for publication in MNRAS. 18 pages, 8 figures, 4 table

We present the results of a radio transient and polarisation survey towards the Galactic Centre, conducted as part of the Australian Square Kilometre Array Pathfinder Variables and Slow Transients pilot survey. The survey region consisted of five fields covering $\sim265\,{\rm deg}^2$ ($350^\circ\lesssim l\lesssim10^\circ$, $\vert b\vert \lesssim 10^\circ$). Each field was observed for 12\,minutes, with between 7 and 9 repeats on cadences of between one day and four months. We detected eight highly variable sources and seven highly circularly-polarised sources (14 unique sources in total). Seven of these sources are known pulsars including the rotating radio transient PSR~J1739--2521 and the eclipsing pulsar PSR~J1723--2837. One of them is a low mass X-ray binary, 4U 1758--25. Three of them are coincident with optical or infrared sources and are likely to be stars. The remaining three may be related to the class of Galactic Centre Radio Transients (including a highly likely one, VAST~J173608.2--321634, that has been reported previously), although this class is not yet understood. In the coming years, we expect to detect $\sim$40 bursts from this kind of source with the proposed four-year VAST survey if the distribution of the source is isotropic over the Galactic fields., Comment: 17 pages, 13 figures, Accepted for publication in MNRAS

Pulsar timing array collaborations have recently reported evidence for a noise process with a common spectrum among the millisecond pulsars in the arrays. The spectral properties of this common-noise process are consistent with expectations for an isotropic gravitational-wave background (GWB) from inspiralling supermassive black-hole binaries. However, recent simulation analyses based on Parkes Pulsar Timing Array data indicate that such a detection may arise spuriously. In this paper, we use simulated pulsar timing array datasets to further test the robustness of the inference methods for spectral and spatial correlations from a GWB. Expanding on our previous results, we find strong support (Bayes factors exceeding $10^5$) for the presence of a common-spectrum noise process in datasets where no common process is present, under a wide range of timing noise prescriptions per pulsar. We show that these results are highly sensitive to the choice of Bayesian priors on timing noise parameters, with priors that more closely match the injected distributions of timing noise parameters resulting in diminished support for a common-spectrum noise process. These results emphasize shortcomings in current methods for inferring the presence of a common-spectrum process, and imply that the detection of a common process is not a reliable precursor to detection of the GWB. Future searches for the nanohertz GWB should remain focussed on detecting spatial correlations, and make use of more tailored specifications for a common-spectrum noise process., Comment: 12 pages, 12 figures. Accepted for publication in MNRAS

Pulsar timing array experiments have recently reported strong evidence for a common-spectrum stochastic process with a strain spectral index consistent with that expected of a nanohertz-frequency gravitational-wave background, but with negligible yet non-zero evidence for spatial correlations required for a definitive detection. However, it was pointed out by the Parkes Pulsar Timing Array (PPTA) collaboration that the same models used in recent analyses resulted in strong evidence for a common-spectrum process in simulations where none is present. In this work, we introduce a methodology to distinguish pulsar power spectra with the same amplitude from noise power spectra of similar but distinct amplitudes. The former is the signature of a spatially uncorrelated pulsar term of a nanohertz gravitational-wave background, whereas the latter could represent ensemble pulsar noise properties. We test the methodology on simulated data sets. We find that the reported common process in PPTA pulsars is indeed consistent with the spectral feature of a pulsar term. We recommend this methodology as one of the validity tests that the real astrophysical and cosmological backgrounds should pass, as well as for inferences about the spatially uncorrelated component of the background., Comment: 9 pages, 2 figures

We describe a new software package for simulating channelised, high-time resolution data streams from radio telescopes. The software simulates data from the telescope and observing system taking into account the observation strategy, receiver system and digitisation. The signatures of pulsars, fast radio bursts and flare stars are modelled, including frequency-dependent effects such as scattering and scintillation. We also simulate more generic signals using spline curves and images. Models of radio frequency interference include signals from satellites, terrestrial transmitters and impulsive, broadband signals. The simulated signals can also be injected into real data sets. Uses of this software include the production of machine learning training data sets, development and testing of new algorithms to search for anomalous patterns and to characterise processing pipelines., Comment: 12 pages, 11 figures, MNRAS in press

We report the discovery of a highly circularly polarized, variable, steep-spectrum pulsar in the Australian Square Kilometre Array Pathfinder (ASKAP) Variables and Slow Transients (VAST) survey. The pulsar is located about $1^\circ$ from the center of the Large Magellanic Cloud, and has a significant fractional circular polarization of $\sim$20%. We discovered pulsations with a period of 322.5 ms, dispersion measure (DM) of 157.5 pc cm$^{-3}$, and rotation measure (RM) of $+456$ rad m$^{-2}$ using observations from the MeerKAT and the Parkes telescopes. This DM firmly places the source, PSR J0523$-$7125, in the Large Magellanic Cloud (LMC). This RM is extreme compared to other pulsars in the LMC (more than twice that of the largest previously reported one). The average flux density of $\sim$1 mJy at 1400 MHz and $\sim$25 mJy at 400 MHz places it among the most luminous radio pulsars known. It likely evaded previous discovery because of its very steep radio spectrum (spectral index $\alpha \approx -3$, where $S_\nu \propto \nu^\alpha$) and broad pulse profile (duty cycle $\gtrsim35$%). We discuss implications for searches for unusual radio sources in continuum images, as well as extragalactic pulsars in the Magellanic Clouds and beyond. Our result highlighted the possibility of identifying pulsars, especially extreme pulsars, from radio continuum images. Future large-scale radio surveys will give us an unprecedented opportunity to discover more pulsars and potentially the most distant pulsars beyond the Magellanic Clouds., Comment: 17 pages, 7 figures; accepted by ApJ

The temperature dependence of the low-energy magnetic excitations in the spin-triplet superconductor UTe$_2$ was measured via inelastic neutron scattering in the normal and superconducting states. The imaginary part of the dynamic susceptibility follows the behavior of interband correlations in a hybridized Kondo lattice with an appropriate characteristic energy. These excitations are a lower-dimensional analogue of phenomena observed in other Kondo lattice materials, such that their presence is not necessarily due to dominance of ferromagnetic or antiferromagnetic correlations. The onset of superconductivity alters the magnetic excitations noticeably on the same energy scales, suggesting that these changes originate from additional electronic structure modification.

The detection of life beyond Earth is an ongoing scientific endeavour, with profound implications. One approach, known as the search for extraterrestrial intelligence (SETI), seeks to find engineered signals (`technosignatures') that indicate the existence technologically-capable life beyond Earth. Here, we report on the detection of a narrowband signal-of-interest at ~982 MHz, recorded during observations toward Proxima Centauri with the Parkes Murriyang radio telescope. This signal, `BLC1', has characteristics broadly consistent with hypothesized technosignatures and is one of the most compelling candidates to date. Analysis of BLC1 -- which we ultimately attribute to being an unusual but locally-generated form of interference -- is provided in a companion paper (Sheikh et al., 2021). Nevertheless, our observations of Proxima Centauri are the most sensitive search for radio technosignatures ever undertaken on a star target., Comment: 14 pages, 4 figures (+3 supplementary figures). Published open-access in Nature Astronomy

The aim of the search for extraterrestrial intelligence (SETI) is to find technologically-capable life beyond Earth through their technosignatures. On 2019 April 29, the Breakthrough Listen SETI project observed Proxima Centauri with the Parkes 'Murriyang' radio telescope. These data contained a narrowband signal with characteristics broadly consistent with a technosignature near 982 MHz ('blc1'). Here we present a procedure for the analysis of potential technosignatures, in the context of the ubiquity of human-generated radio interference, which we apply to blc1. Using this procedure, we find that blc1 is not an extraterrestrial technosignature, but rather an electronically-drifting intermodulation product of local, time-varying interferers aligned with the observing cadence. We find dozens of instances of radio interference with similar morphologies to blc1 at frequencies harmonically related to common clock oscillators. These complex intermodulation products highlight the necessity for detailed follow-up of any signal-of-interest using a procedure such as the one outlined in this work., Comment: 15 pages of main paper followed by 19 pages of supplementary material; 22 figures total; published open-access in Nature Astronomy

Observations of the extragalactic ($z=0.0141$) transient AT2018cow established a new class of energetic explosions shocking a dense medium, which produce luminous emission at millimeter and sub-millimeter wavelengths. Here we present detailed millimeter- through centimeter-wave observations of a similar transient, ZTF20acigmel (AT2020xnd) at $z=0.2433$. Using observations from the NOrthern Extended Millimeter Array and the Very Large Array, we model the unusual millimeter and radio emission from AT2020xnd under several different assumptions, and ultimately favor synchrotron radiation from a thermal electron population (relativistic Maxwellian). The thermal-electron model implies a fast but sub-relativistic ($v\approx0.3c$) shock and a high ambient density ($n_e\approx4\times10^{3}$cm$^{-3}$ at $\Delta t\approx40$ days). The X-ray luminosity of $L_X\approx10^{43}$ erg sec$^{-1}$ exceeds simple predictions from the radio and UVOIR luminosity and likely has a separate physical origin, such as a central engine. Using the fact that month-long luminous ($L_\nu\approx 2\times10^{30}$ erg sec$^{-1}$ Hz$^{-1}$ at 100 GHz) millimeter emission appears to be a generic feature of transients with fast ($t_{1/2}\approx3$ days) and luminous ($M_\mathrm{peak}\approx -21 $mag) optical light curves, we estimate the rate at which transients like AT2018cow and AT2020xnd will be detected by future wide-field millimeter transient surveys such as CMB-S4, and conclude that energetic explosions in dense environments may represent a significant population of extragalactic transients in the 100 GHz sky., Comment: 27 pages, 13 figures. Accepted for publication in ApJ on 23 Jan 2022

We report the discovery of a highly-polarized, highly-variable, steep-spectrum radio source, ASKAP J173608.2-321635, located $\sim$4\,deg from the Galactic center in the Galactic plane. The source was detected six times between 2020 January and 2020 September as part of the Australian Square Kilometre Array Pathfinder Variables and Slow Transients (ASKAP VAST) survey at 888\,MHz. It exhibited a high degree ($\sim 25$\%) of circular polarization when it was visible. We monitored the source with the MeerKAT telescope from 2020 November to 2021 February on a 2--4 week cadence. The source was not detected with MeerKAT before 2021 February 07 when it appeared and reached a peak flux density of 5.6\,mJy. The source was still highly circularly polarized, but also showed up to 80\% linear polarization, and then faded rapidly with a timescale of one day. The rotation measure of the source varied significantly, from $-11.8\pm0.8$\,rad\,m$^{-2}$ to $-64.0\pm1.5$\,rad\,m$^{-2}$, over three days. No X-ray counterpart was found in follow-up \textit{Swift} or \textit{Chandra} observations about a week after the first MeerKAT detection, with upper limits of $\sim 5.0\times10^{31}$\,erg\,s$^{-1}$ (0.3--8\,keV, assuming a distance $\sim10$ kpc). No counterpart is seen in new or archival near-infrared observations down to $J=20.8$\,mag. We discuss possible identifications for ASKAP J173608.2-321635 including a low-mass star/substellar object with extremely low infrared luminosity, a pulsar with scatter-broadened pulses, a transient magnetar, or a Galactic Center Radio Transient: none of these fully explains the observations, which suggests that ASKAP J173608.2-321635 may represent part of a new class of objects being discovered through radio imaging surveys., Comment: 22 pages, 15 figures, Accepted by ApJ

A nanohertz-frequency stochastic gravitational-wave background can potentially be detected through the precise timing of an array of millisecond pulsars. This background produces low-frequency noise in the pulse arrival times that would have a characteristic spectrum common to all pulsars and a well-defined spatial correlation. Recently the North American Nanohertz Observatory for Gravitational Waves collaboration (NANOGrav) found evidence for the common-spectrum component in their 12.5-year data set. Here we report on a search for the background using the second data release of the Parkes Pulsar Timing Array. If we are forced to choose between the two NANOGrav models $\unicode{x2014}$ one with a common-spectrum process and one without $\unicode{x2014}$ we find strong support for the common-spectrum process. However, in this paper, we consider the possibility that the analysis suffers from model misspecification. In particular, we present simulated data sets that contain noise with distinctive spectra but show strong evidence for a common-spectrum process under the standard assumptions. The Parkes data show no significant evidence for, or against, the spatially correlated Hellings-Downs signature of the gravitational-wave background. Assuming we did observe the process underlying the spatially uncorrelated component of the background, we infer its amplitude to be $A = 2.2^{+0.4}_{-0.3} \times 10^{-15}$ in units of gravitational-wave strain at a frequency of $1\, \text{yr}^{-1}$. Extensions and combinations of existing and new data sets will improve the prospects of identifying spatial correlations that are necessary to claim a detection of the gravitational-wave background., Comment: 10 pages, 5 figures

We present the discovery of an extreme flaring event from Proxima Cen by ASKAP, ALMA, HST, TESS, and the du Pont Telescope that occurred on 2019 May 1. In the millimeter and FUV, this flare is the brightest ever detected, brightening by a factor of >1000 and >14000 as seen by ALMA and HST, respectively. The millimeter and FUV continuum emission trace each other closely during the flare, suggesting that millimeter emission could serve as a proxy for FUV emission from stellar flares and become a powerful new tool to constrain the high-energy radiation environment of exoplanets. Surprisingly, optical emission associated with the event peaks at a much lower level with a time delay. The initial burst has an extremely short duration, lasting for <10 sec. Taken together with the growing sample of millimeter M dwarf flares, this event suggests that millimeter emission is actually common during stellar flares and often originates from short burst-like events., Comment: 11 pages, 4 figures, 1 appendix, published in ApJ Letters