Search

Accurate and robust modelling of large deformation three dimensional contact interaction is an important area of engineering, but it is also challenging from a computational mechanics perspective. This is particularly the case when there is significant interpenetration and evolution of the contact surfaces, such as the case of a relatively rigid body interacting with a highly deformable body. The numerical challenges come from several non-linear sources: large deformation mechanics, history dependent material behaviour and slip/stick frictional contact. In this paper the Material Point Method (MPM) is adopted to represent the deformable material, combined with a discretised rigid body which provides an accurate representation of the contact surface. The three dimensional interaction between the bodies is detected though the use of domains associated with each material point. This provides a general and consistent representation of the extent of the deformable body without introducing boundary representation in the material point method. The dynamic governing equations allows the trajectory of the rigid body to evolve based on the interaction with the deformable body and the governing equations are solved within an efficient implicit framework. The performance of the method is demonstrated on a number of benchmark problems with analytical solutions. The method is also applied to the specific case of soil-structure interaction, using geotechnical centrifuge experimental data that confirms the veracity of the proposed approach.

The color code has been invaluable for the development of the theory of fault-tolerant logic gates using transversal rotations. Three-dimensional examples of the color code have shown us how its structure, specifically the intersection of the supports of logical operators, can give rise to non-Clifford $T$ and $CCZ$ gates. Here we present a color code with a logical control-$S$ gate that is accomplished with transversal $T$ and $T^\dagger$ rotations on its physical qubits., Comment: 7 pages, 2 figures, comments welcome

We report the discovery of three ultracompact binary white dwarf systems hosting accretion disks, with orbital periods of 7.95, 8.68, and 13.15 minutes. This significantly augments the population of mass-transferring binaries at the shortest periods, and provides the first evidence that accretors in ultracompacts can be dense enough to host accretion disks even below 10 minutes (where previously only direct-impact accretors were known). In the two shortest-period systems, we measured changes in the orbital periods driven by the combined effect of gravitational wave emission and mass transfer; we find $\dot{P}$ is negative in one case, and positive in the other. This is only the second system measured with a positive $\dot{P}$, and it the most compact binary known that has survived a period minimum. Using these systems as examples, we show how the measurement of $\dot{P}$ is a powerful tool in constraining the physical properties of binaries, e.g. the mass and mass-radius relation of the donor stars. We find that the chirp masses of ultracompact binaries at these periods seem to cluster around $\mathcal{M}_c \sim 0.3 M_\odot$, perhaps suggesting a common origin for these systems or a selection bias in electromagnetic discoveries. Our new systems are among the highest-amplitude known gravitational wave sources in the millihertz regime, providing exquisite opportunity for multi-messenger study with future space-based observatories such as \textit{LISA} and TianQin; we discuss how such systems provide fascinating laboratories to study the unique regime where the accretion process is mediated by gravitational waves., Comment: Accepted for publication in ApJ

Different choices of quantum error-correcting codes can reduce the demands on the physical hardware needed to build a quantum computer. To achieve the full potential of a code, we must develop practical decoding algorithms that can correct errors that have occurred with high likelihood. Matching decoders are very good at correcting local errors while also demonstrating fast run times that can keep pace with physical quantum devices. We implement variations of a matching decoder for a three-dimensional, non-CSS, low-density parity check code known as the Chamon code, which has a non-trivial structure that does not lend itself readily to this type of decoding. The non-trivial structure of the syndrome of this code means that we can supplement the decoder with additional steps to improve the threshold error rate, below which the logical failure rate decreases with increasing code distance. We find that a generalized matching decoder that is augmented by a belief-propagation step prior to matching gives a threshold of 10.5% for depolarising noise., Comment: 9 pages, 7 figures

In this work we present source-tailored WISE mid-infrared photometry (at 3.4$\mu$m, 4.6$\mu$m, 12$\mu$m, and 23$\mu$m) of 2812 galaxies in the extended Spitzer Survey of Stellar Structure in Galaxies (S$^4$G) sample, and characterise the mid-infrared colors and dust properties of this legacy nearby galaxy data set. Informed by the relative emission between W3 (12$\mu$ m) and W4 (23$\mu$ m), we re-derive star formation rate (SFR) scaling relations calibrated to L$_{\rm TIR}$, which results in improved agreement between the two tracers. By inverse-variance weighting the W3 and W4-derived SFRs, we generate a combined mid-infrared SFR that is a broadly robust measure of star formation activity in dusty, star-forming galaxies in the nearby Universe. In addition, we investigate the use of a W3-derived dust density metric, $\Sigma_{\rm 12\mu m}$ (L$_\odot$/kpc$^2$), to estimate the SFR deficit of low mass, low dust galaxies. This is achieved by combining WISE with existing GALEX ultraviolet (UV) photometry, which we further use to explore the relationship between dust and UV emission as a function of morphology. Finally, we use our derived SFR prescriptions to examine the location of galaxies in the log SFR - log M$_\textrm{stellar}$ plane, as a function of morphological type, which underscores the complexity of dust-derived properties seen in galaxies of progressively earlier type., Comment: Accepted to ApJ

This paper introduces the Bayesian Confidence Estimator (BACON) for deep neural networks. Current practice of interpreting Softmax values in the output layer as probabilities of outcomes is prone to extreme predictions of class probability. In this work we extend Waagen's method of representing the terminal layers with a geometric model, where the probability associated with an output vector is estimated with Bayes' Rule using validation data to provide likelihood and normalization values. This estimator provides superior ECE and ACE calibration error compared to Softmax for ResNet-18 at 85% network accuracy, and EfficientNet-B0 at 95% network accuracy, on the CIFAR-10 dataset with an imbalanced test set, except for very high accuracy edge cases. In addition, when using the ACE metric, BACON demonstrated improved calibration error when estimating probabilities for the imbalanced test set when using actual class distribution fractions., Comment: 14 pages, 15 figures (10 of which include sub-figures)

The diffuse stellar component of galaxy clusters made up of intergalactic stars is termed the intracluster light (ICL). Though there is a developing understanding of the mechanisms by which the ICL is formed, no strong consensus has yet been reached on which objects the stars of the ICL are primarily sourced from. We investigate the assembly of the ICL starting approximately $10$ Gyr before $z=0$ in 11 galaxy clusters (halo masses between $\sim1\times 10^{14}$ M$_{\odot}$ and $\sim7\times 10^{14}$ M$_{\odot}$ at $z\approx0$) in the Horizon-AGN simulation. By tracking the stars of galaxies that fall into these clusters past cluster infall, we are able to link almost all of the $z\approx0$ ICL back to progenitor objects. Satellite stripping, mergers, and pre-processing are all found to make significant contributions to the ICL, but any contribution from in-situ star-formation directly into the ICL appears negligible. Even after compensating for resolution effects, we find that approximately $90$ per cent of the stacked ICL of the 11 clusters that is not pre-processed should come from galaxies infalling with stellar masses above $10^{9}$ M$_{\odot}$, with roughly half coming from infalling galaxies with stellar masses within half a dex of $10^{11}$ M$_{\odot}$. The fact that the ICL appears largely sourced from such massive objects suggests that the ICL assembly of any individual cluster may be principally stochastic., Comment: 11 pages, 8 figures, accepted for publication in MNRAS

Fault-tolerant implementation of non-Clifford gates is a major challenge for achieving universal fault-tolerant quantum computing with quantum error-correcting codes. Magic state distillation is the most well-studied method for this but requires significant resources. Hence, it is crucial to tailor and optimize magic state distillation for specific codes from both logical- and physical-level perspectives. In this work, we perform such optimization for two-dimensional color codes, which are promising due to their higher encoding rates compared to surface codes, transversal implementation of Clifford gates, and efficient lattice surgery. We propose two distillation schemes based on the 15-to-1 distillation circuit and lattice surgery, which differ in their methods for handling faulty rotations. Our first scheme uses faulty T-measurement, offering resource efficiency when the target infidelity is above a certain threshold ($\sim 35p^3$ for physical error rate $p$). To achieve lower infidelities while maintaining resource efficiency, our second scheme exploits a distillation-free fault-tolerant magic state preparation protocol, achieving significantly lower infidelities (e.g., $\sim 10^{-19}$ for $p = 10^{-4}$) than the first scheme. Notably, our schemes outperform the best existing magic state distillation methods for color codes by up to about two orders of magnitude in resource costs for a given achievable target infidelity., Comment: 42 pages (22 pages for main text), 21 figures, 3 tables; v2 - updated combined MSD scheme (without autocorrection qubits) thanks to Sam Roberts's suggestion & additional comparison with a previous color code MSD scheme in Fig. 14

In this work, we apply state-of-the-art self-supervised learning techniques on a large dataset of seafloor imagery, \textit{BenthicNet}, and study their performance for a complex hierarchical multi-label (HML) classification downstream task. In particular, we demonstrate the capacity to conduct HML training in scenarios where there exist multiple levels of missing annotation information, an important scenario for handling heterogeneous real-world data collected by multiple research groups with differing data collection protocols. We find that, when using smaller one-hot image label datasets typical of local or regional scale benthic science projects, models pre-trained with self-supervision on a larger collection of in-domain benthic data outperform models pre-trained on ImageNet. In the HML setting, we find the model can attain a deeper and more precise classification if it is pre-trained with self-supervision on in-domain data. We hope this work can establish a benchmark for future models in the field of automated underwater image annotation tasks and can guide work in other domains with hierarchical annotations of mixed resolution.

The key difficulty faced by 2D models for planet-disc interaction is in appropriately accounting for the impact of the disc's vertical structure on the dynamics. 3D effects are often mimicked via softening of the planet's potential; however, the planet-induced flow and torques often depend strongly on the choice of softening length. We show that for a linear adiabatic flow perturbing a vertically isothermal disc, there is a particular vertical average of the 3D equations of motion which exactly reproduces 2D fluid equations for arbitrary adiabatic index. There is a strong connection here with the Lubow-Pringle 2D mode of the disc. Correspondingly, we find a simple, general prescription for the consistent treatment of planetary potentials embedded within '2D' discs. The flow induced by a low-mass planet involves large-scale excited spiral density waves which transport angular momentum radially away from the planet, and 'horseshoe streamlines' within the co-orbital region. We derive simple linear equations governing the flow which locally capture both effects faithfully simultaneously. We present an accurate co-orbital flow solution allowing for inexpensive future study of corotation torques, and predict the vertical structure of the co-orbital flow and horseshoe region width for different values of adiabatic index, as well as the vertical dependence of the initial shock location. We find strong agreement with the flow computed in 3D numerical simulations, and with 3D one-sided Lindblad torque estimates, which are a factor of 2 to 3 times lower than values from previous 2D simulations., Comment: Accepted MNRAS, 18 pages, 6 figures

UW Coronae Borealis (UW CrB) is a low mass X-ray binary that shows both Type 1 X-ray and optical bursts, which typically last for 20 s. The system has a binary period of close to 2 hours and is thought to have a relatively high inclination due to the presence of an eclipse in the optical light curve. There is also evidence that an asymmetric disc is present in the system, which precesses every 5.5 days based on changes in the depth of the eclipse. In this paper, we present optical photometry and spectroscopy of UW CrB taken over 2 years. We update the orbital ephemeris using observed optical eclipses and refine the orbital period to 110.97680(1) min. A total of 17 new optical bursts are presented, with 10 of these bursts being resolved temporally. The average $e$-folding time of $19\pm3$s for the bursts is consistent with the previously found value. Optical bursts are observed during a previously identified gap in orbital phase centred on $\phi=0.967$, meaning the reprocessing site is not eclipsed as previously thought. Finally, we find that the apparent P-Cygni profiles present in some of the atomic lines in the optical spectra are due to transient absorption., Comment: 14 pages, 17 figures, submitted to OJAp

We propose a technique for frequency locking a laser to the Zeeman sublevel transitions between the 5P$_{3/2}$ intermediate and 32D$_{5/2}$ Rydberg states in $^{87}$Rb. This method allows for continuous frequency tuning over 0.6 GHz by varying an applied external magnetic field. In the presence of the applied field, the electromagnetically induced transparency (EIT) spectrum of an atomic vapor splits via the Zeeman effect according to the strength of the magnetic field and the polarization of the pump and probe lasers. We show that the 480 nm pump laser, responsible for transitions between the Zeeman sublevels of the intermediate state and the Rydberg state, can be locked to the Zeeman-split EIT peaks. The short-term frequency stability of the laser lock is 0.15 MHz and the long-term stability is within 0.5 MHz. The linewidth of the laser lock is ~0.8 MHz and ~1.8 MHz in the presence and in the absence of the external magnetic field, respectively. In addition, we show that in the absence of an applied magnetic field and adequate shielding, the frequency shift of the lock point has a peak-to-peak variation of 1.6 MHz depending on the polarization of the pump field, while when locked to Zeeman sublevels this variation is reduced to 0.6 MHz. The proposed technique is useful for research involving Rydberg atoms, where large continuous tuning of the laser frequency with stable locking is required.

PLATO (PLAnetary Transits and Oscillations of stars) is ESA's M3 mission designed to detect and characterise extrasolar planets and perform asteroseismic monitoring of a large number of stars. PLATO will detect small planets (down to <2 R_(Earth)) around bright stars (<11 mag), including terrestrial planets in the habitable zone of solar-like stars. With the complement of radial velocity observations from the ground, planets will be characterised for their radius, mass, and age with high accuracy (5 %, 10 %, 10 % for an Earth-Sun combination respectively). PLATO will provide us with a large-scale catalogue of well-characterised small planets up to intermediate orbital periods, relevant for a meaningful comparison to planet formation theories and to better understand planet evolution. It will make possible comparative exoplanetology to place our Solar System planets in a broader context. In parallel, PLATO will study (host) stars using asteroseismology, allowing us to determine the stellar properties with high accuracy, substantially enhancing our knowledge of stellar structure and evolution. The payload instrument consists of 26 cameras with 12cm aperture each. For at least four years, the mission will perform high-precision photometric measurements. Here we review the science objectives, present PLATO's target samples and fields, provide an overview of expected core science performance as well as a description of the instrument and the mission profile at the beginning of the serial production of the flight cameras. PLATO is scheduled for a launch date end 2026. This overview therefore provides a summary of the mission to the community in preparation of the upcoming operational phases.

We present high-cadence photometric and spectroscopic observations of supernova (SN) 2024ggi, a Type II SN with flash spectroscopy features which exploded in the nearby galaxy NGC 3621 at $\sim$7 Mpc. The light-curve evolution over the first 30 hours can be fit by two power law indices with a break after 22 hours, rising from $M_V \approx -12.95$ mag at +0.66 days to $M_V \approx -17.91$ mag after 7 days. In addition, the densely sampled color curve shows a strong blueward evolution over the first few days and then behaves as a normal SN II with a redward evolution as the ejecta cool. Such deviations could be due to interaction with circumstellar material (CSM). Early high- and low-resolution spectra clearly show high-ionization flash features from the first spectrum to +3.42 days after the explosion. From the high-resolution spectra, we calculate the CSM velocity to be 37 $\pm~4~\mathrm{km\,s^{-1}} $. We also see the line strength evolve rapidly from 1.22 to 1.49 days in the earliest high-resolution spectra. Comparison of the low-resolution spectra with CMFGEN models suggests that the pre-explosion mass-loss rate of SN 2024ggi falls in a range of $10^{-3}$ to $10^{-2}$ M$_{\odot}$ yr$^{-1}$, which is similar to that derived for SN 2023ixf. However, the rapid temporal evolution of the narrow lines in the spectra of SN 2024ggi ($R_\mathrm{CSM} \sim 2.7 \times 10^{14} \mathrm{cm}$) could indicate a smaller spatial extent of the CSM than in SN 2023ixf ($R_\mathrm{CSM} \sim 5.4 \times 10^{14} \mathrm{cm}$) which in turn implies lower total CSM mass for SN 2024ggi., Comment: 23 pages, 15 figures, 4 tables, accepted for publication in ApJL

Understanding the nature of the luminous 1991T-like supernovae is of great importance to supernova cosmology as they are likely to have been more common in the early universe. In this paper we explore the observational properties of 1991T-like supernovae to study their relationship to other luminous, slow-declining Type~Ia supernovae (SNe Ia). From the spectroscopic and photometric criteria defined in Phillips et al. (1992), we identify 17 1991T-like supernovae from the literature. Combining these objects with ten 1991T-like supernovae from the Carnegie Supernova Project-II, the spectra, light curves, and colors of these events, along with their host galaxy properties, are examined in detail. We conclude that 1991T-like supernovae are closely related in essentially all of their UV, optical, and near-infrared properties -- as well as their host galaxy parameters -- to the slow-declining subset of Branch core-normal supernovae and to the intermediate 1999aa-like events, forming a continuum of luminous SNe Ia. The overriding difference between these three subgroups appears to be the extent to which $^{56}$Ni mixes into the ejecta, producing the pre-maximum spectra dominated by Fe III absorption, the broader UV light curves, and the higher luminosities that characterize the 1991T-like events. Nevertheless, the association of 1991T-like SNe with the rare Type Ia CSM supernovae would seem to run counter to this hypothesis, in which case 1991T-like events may form a separate subclass of SNe Ia, possibly arising from single-degenerate progenitor systems., Comment: Accepted for publication in ApJS

We present new frontiers in the modelling of the spectral energy distributions (SED) of active galaxies by introducing the radio-to-X-ray fitting capabilities of the publicly available Bayesian code AGNfitter. The new code release, called AGNfitter-rx, models the broad-band photometry covering the radio, infrared (IR), optical, ultraviolet (UV) and X-ray bands consistently, using a combination of theoretical and semi-empirical models of the AGN and host galaxy emission. This framework enables the detailed characterization of four physical components of the active nuclei: the accretion disk, the hot dusty torus, the relativistic jets/core radio emission, and the hot corona; alongside modeling three components within the host galaxy: stellar populations, cold dust, and the radio emission from the star-forming regions. Applying AGNfitter-rx to a diverse sample of 36 AGN SEDs at z<0.7 from the AGN SED ATLAS, we investigate and compare the performance of state-of-the-art torus and accretion disk emission models on fit quality and inferred physical parameters. We find that clumpy torus models that include polar winds and semi-empirical accretion disk templates including emission line features significantly increase the fit quality in 67% of the sources, by effectively reducing by $2\sigma$ fit residuals in the $1.5-5 \mu \rm m$ and $0.7 \mu \rm m$ regimes.We demonstrate that, by applying AGNfitter-rx on photometric data, we are able to estimate inclination and opening angles of the torus, consistent with spectroscopic classifications within the AGN unified model, as well as black hole mass estimates in agreement with virial estimates based on H$\alpha$. The wavelength coverage and the flexibility for the inclusion of state-of-the-art theoretical models make AGNfitter-rx a unique tool for the further development of SED modelling for AGNs in present and future radio-to-X-ray galaxy surveys., Comment: 20 pages, 10 figures, Accepted for publication by A&A

Monolithic large language models (LLMs) like GPT-4 have paved the way for modern generative AI applications. Training, serving, and maintaining monolithic LLMs at scale, however, remains prohibitively expensive and challenging. The disproportionate increase in compute-to-memory ratio of modern AI accelerators have created a memory wall, necessitating new methods to deploy AI. Composition of Experts (CoE) is an alternative modular approach that lowers the cost and complexity of training and serving. However, this approach presents two key challenges when using conventional hardware: (1) without fused operations, smaller models have lower operational intensity, which makes high utilization more challenging to achieve; and (2) hosting a large number of models can be either prohibitively expensive or slow when dynamically switching between them. In this paper, we describe how combining CoE, streaming dataflow, and a three-tier memory system scales the AI memory wall. We describe Samba-CoE, a CoE system with 150 experts and a trillion total parameters. We deploy Samba-CoE on the SambaNova SN40L Reconfigurable Dataflow Unit (RDU) - a commercial dataflow accelerator architecture that has been co-designed for enterprise inference and training applications. The chip introduces a new three-tier memory system with on-chip distributed SRAM, on-package HBM, and off-package DDR DRAM. A dedicated inter-RDU network enables scaling up and out over multiple sockets. We demonstrate speedups ranging from 2$\times$ to 13$\times$ on various benchmarks running on eight RDU sockets compared with an unfused baseline. We show that for CoE inference deployments, the 8-socket RDU Node reduces machine footprint by up to 19$\times$, speeds up model switching time by 15$\times$ to 31$\times$, and achieves an overall speedup of 3.7$\times$ over a DGX H100 and 6.6$\times$ over a DGX A100., Comment: 2024 57th IEEE/ACM International Symposium on Microarchitecture (MICRO)

Advances in underwater imaging enable the collection of extensive seafloor image datasets that are necessary for monitoring important benthic ecosystems. The ability to collect seafloor imagery has outpaced our capacity to analyze it, hindering expedient mobilization of this crucial environmental information. Recent machine learning approaches provide opportunities to increase the efficiency with which seafloor image datasets are analyzed, yet large and consistent datasets necessary to support development of such approaches are scarce. Here we present BenthicNet: a global compilation of seafloor imagery designed to support the training and evaluation of large-scale image recognition models. An initial set of over 11.4 million images was collected and curated to represent a diversity of seafloor environments using a representative subset of 1.3 million images. These are accompanied by 2.6 million annotations translated to the CATAMI scheme, which span 190,000 of the images. A large deep learning model was trained on this compilation and preliminary results suggest it has utility for automating large and small-scale image analysis tasks. The compilation and model are made openly available for use by the scientific community at https://doi.org/10.20383/103.0614.

Quantum error correcting codes protect quantum information, allowing for large quantum computations provided that physical error rates are sufficiently low. We combine post-selection with surface code error correction through the use of a parameterized family of exclusive decoders, which are able to abort on decoding instances that are deemed too difficult. We develop new numerical sampling methods to quantify logical failure rates with exclusive decoders as well as the trade-off in terms of the amount of post-selection required. For the most discriminating of exclusive decoders, we demonstrate a threshold of 50\% under depolarizing noise for the surface code (or $32(1)\%$ for the fault-tolerant case with phenomenological measurement errors), and up to a quadratic improvement in logical failure rates below threshold. Furthermore, surprisingly, with a modest exclusion criterion, we identify a regime at low error rates where the exclusion rate decays with code distance, providing a pathway for scalable and time-efficient quantum computing with post-selection. We apply our exclusive decoder to the 15-to-1 magic state distillation protocol, and report a $75\%$ reduction in the number of physical qubits required, and a $60\%$ reduction in the total spacetime volume required, including accounting for repetitions required for post-selection. We also consider other applications, as an error mitigation technique, and in concatenated schemes. Our work highlights the importance of post-selection as a powerful tool in quantum error correction., Comment: 20 pages, 17 figures, comments welcome

We present a JWST MIRI/MRS spectrum (5-27 $\mathrm{\mu}$m) of the Type Ia supernova (SN Ia), SN 2021aefx at $+415$ days past $B$-band maximum. The spectrum, which was obtained during the iron-dominated nebular phase, has been analyzed in combination with previous JWST observations of SN 2021aefx, to provide the first JWST time series analysis of an SN Ia. We find the temporal evolution of the [Co III] 11.888 $\mathrm{\mu}$m feature directly traces the decay of $^{56}$Co. The spectra, line profiles, and their evolution are analyzed with off-center delayed-detonation models. Best fits were obtained with White Dwarf (WD) central densities of $\rho_c=0.9-1.1\times 10^9$g cm$^{-3}$, a WD mass of M$_{\mathrm{WD}}$=1.33-1.35M$_\odot$, a WD magnetic field of $\approx10^6$G, and an off-center deflagration-to-detonation transition at $\approx$ 0.5 $M_\odot$ seen opposite to the line of sight of the observer (-30). The inner electron capture core is dominated by energy deposition from $\gamma$-rays whereas a broader region is dominated by positron deposition, placing SN 2021aefx at +415 d in the transitional phase of the evolution to the positron-dominated regime. The formerly `flat-tilted' profile at 9 $\mathrm{\mu}$m now has significant contribution from [Ni IV], [Fe II], and [Fe III] and less from [Ar III], which alters the shape of the feature as positrons excite mostly the low-velocity Ar. Overall, the strength of the stable Ni features in the spectrum is dominated by positron transport rather than the Ni mass. Based on multi-dimensional models, our analysis is consistent with a single-spot, close-to-central ignition with an indication for a pre-existing turbulent velocity field, and excludes a multiple-spot, off-center ignition., Comment: Accepted for publication in ApJ

We present a method of extrapolating the spectroscopic behavior of Type Ia supernovae (SNe Ia) in the near-infrared (NIR) wavelength regime up to 2.30 $\mu$m using optical spectroscopy. Such a process is useful for accurately estimating K-corrections and other photometric quantities of SNe Ia in the NIR. Principal component analysis is performed on data consisting of Carnegie Supernova Project I & II optical and near-infrared FIRE spectra to produce models capable of making these extrapolations. This method differs from previous spectral template methods by not parameterizing models strictly by photometric light-curve properties of SNe Ia, allowing for more flexibility of the resulting extrapolated NIR flux. A difference of around -3.1% to -2.7% in the total integrated NIR flux between these extrapolations and the observations is seen here for most test cases including Branch core-normal and shallow-silicon subtypes. However, larger deviations from the observation are found for other tests, likely due to the limited high-velocity and broad-line SNe Ia in the training sample. Maximum-light principal components are shown to allow for spectroscopic predictions of the color-stretch light-curve parameter, $s_{BV}$, within approximately $\pm$0.1 units of the value measured with photometry. We also show these results compare well with NIR templates, although in most cases the templates are marginally more fitting to observations, illustrating a need for more concurrent optical+NIR spectroscopic observations to truly understand the diversity of SNe Ia in the NIR., Comment: 25 pages, 16 figures, ApJ, in press

We propose a new method for identifying active galactic nuclei (AGN) in low mass ($\rm M_*\leq10^{10}M_\odot$) galaxies. This method relies on spectral energy distribution (SED) fitting to identify galaxies whose radio flux density has an excess over that expected from star formation alone. Combining data in the Galaxy and Mass Assembly (GAMA) G23 region from GAMA, Evolutionary Map of the Universe (EMU) early science observations, and Wide-field Infrared Survey Explorer (WISE), we compare this technique with a selection of different AGN diagnostics to explore the similarities and differences in AGN classification. We find that diagnostics based on optical and near-infrared criteria (the standard BPT diagram, the WISE colour criterion, and the mass-excitation, or MEx diagram) tend to favour detection of AGN in high mass, high luminosity systems, while the ``ProSpect'' SED fitting tool can identify AGN efficiently in low mass systems. We investigate an explanation for this result in the context of proportionally lower mass black holes in lower mass galaxies compared to higher mass galaxies and differing proportions of emission from AGN and star formation dominating the light at optical and infrared wavelengths as a function of galaxy stellar mass. We conclude that SED-derived AGN classification is an efficient approach to identify low mass hosts with low radio luminosity AGN., Comment: 16 pages, 12 figures, 5 tables. Accepted for publication in PASA

Radio emission has been detected from tens of white dwarfs, in particular in accreting systems. Additionally, radio emission has been predicted as a possible outcome of a planetary system around a white dwarf. We searched for 3 GHz radio continuum emission in 846,000 candidate white dwarfs previously identified in Gaia using the Very Large Array Sky Survey (VLASS) Epoch 1 Quick Look Catalogue. We identified 13 candidate white dwarfs with a counterpart in VLASS within 2". Five of those were found not to be white dwarfs in follow-up or archival spectroscopy, whereas seven others were found to be chance alignments with a background source in higher-resolution optical or radio images. The remaining source, WDJ204259.71+152108.06, is found to be a white dwarf and M-dwarf binary with an orbital period of 4.1 days and long-term stochastic optical variability, as well as luminous radio and X-ray emission. For this binary, we find no direct evidence of a background contaminant, and a chance alignment probability of only ~2 per cent. However, other evidence points to the possibility of an unfortunate chance alignment with a background radio and X-ray emitting quasar, including an unusually poor Gaia DR3 astrometric solution for this source. With at most one possible radio emitting white dwarf found, we conclude that strong (> 1-3 mJy) radio emission from white dwarfs in the 3 GHz band is virtually nonexistent outside of interacting binaries., Comment: 18 pages, 17 figures. Updated to match version accepted for publication in MNRAS

Approximately 150 low-mass white dwarfs, with masses below 0.4Msun, have been discovered. The majority of these low-mass WDs are observed in binary systems as they cannot be formed through single-star evolution within the Hubble time. In this study, we present a comprehensive analysis of the double low-mass WD eclipsing binary system J2102-4145. Our investigation involved an extensive observational campaign, resulting in the acquisition of approximately 28 hours of high-speed photometric data across multiple nights using NTT/ULTRACAM, SOAR/Goodman, and SMARTS-1m telescopes. These observations have provided critical insights into the orbital characteristics of this system, including parameters such as inclination and orbital period. To disentangle the binary components of J2102-4145, we employed the XT GRID spectral fitting method with GMOS/Gemini-South and X-Shooter data. Additionally, we used the PHOEBE package for light curve analysis on NTT/ULTRACAM high-speed time-series photometry data to constrain the binary star properties. Our analysis reveals remarkable similarities between the two components of this binary system. For the primary star, we determined Teff1 = 13688 +- 65 K, log g1 = 7.36 +- 0.01, R1 = 0.0211 +- 0.0002 Rsun, and M1 = 0.375 +- 0.003 Msun, while the secondary star is characterized by Teff2 = 12952 +- 53 K, log g2 = 7.32 +- 0.01, R2 = 0.0203 +- 0.0002 Rsun, and M2 = 0.31 +- 0.003 Msun. Furthermore, we observe a notable discrepancy between Teff and R of the less massive WD compared to evolutionary sequences for WDs from the literature, which has significant implications for our understanding of WD evolution. We discuss a potential formation scenario for this system that might explain this discrepancy and explore its future evolution. We predict that this system will merge in about 800 Myr, evolving into a helium-rich hot subdwarf star and later into a hybrid He/CO WD.

The feasibility of extrapolation of completely monotone functions can be quantified by examining the worst case scenario, whereby a pair of completely monotone functions agree on a given interval to a given relative precision, but differ as much as it is theoretically possible at a given point. We show that extrapolation is impossible to the left of the interval, while the maximal discrepancy to the right exhibits a power law typical for extrapolation of similar classes of complex analytic functions. The power law exponent is derived explicitly, and shows a precipitous drop immediately beyond the right end-point, with a subsequent decay to zero inversely proportional to the distance from the interval. The local extrapolation problem, where the worst discrepancy from a given completely monotone function is sought, is also analyzed. In this case explicit and easily verifiable optimality conditions are derived, enabling us to solve the problem exactly for a single decaying exponential. In the general case, our approach leads to a natural algorithm for computing solutions to the local extrapolation problem numerically. The methods developed in this paper can easily be adapted to other classes of analytic functions represented as integral transforms of positive measures with analytic kernels., Comment: 40 pages, 1 figure

We report the All-Sky Automated Survey for SuperNovae discovery of the tidal disruption event (TDE) ASASSN-23bd (AT 2023clx) in NGC 3799, a LINER galaxy with no evidence of strong AGN activity over the past decade. With a redshift of $z = 0.01107$ and a peak UV/optical luminosity of $(5.4\pm0.4)\times10^{42}$ erg s$^{-1}$, ASASSN-23bd is the lowest-redshift and least-luminous TDE discovered to date. Spectroscopically, ASASSN-23bd shows H$\alpha$ and He I emission throughout its spectral time series, and the UV spectrum shows nitrogen lines without the strong carbon and magnesium lines typically seen for AGN. Fits to the rising ASAS-SN light curve show that ASASSN-23bd started to brighten on MJD 59988$^{+1}_{-1}$, $\sim$9 days before discovery, with a nearly linear rise in flux, peaking in the $g$ band on MJD $60000^{+3}_{-3}$. Scaling relations and TDE light curve modelling find a black hole mass of $\sim$10$^6$ $M_\odot$, which is on the lower end of supermassive black hole masses. ASASSN-23bd is a dim X-ray source, with an upper limit of $L_{0.3-10\,\mathrm{keV}} < 1.0\times10^{40}$ erg s$^{-1}$ from stacking all \emph{Swift} observations prior to MJD 60061, but with soft ($\sim 0.1$ keV) thermal emission with a luminosity of $L_{0.3-2 \,\mathrm{keV}}\sim4\times10^{39}$ erg s$^{-1}$ in \emph{XMM-Newton} observations on MJD 60095. The rapid $(t < 15$ days) light curve rise, low UV/optical luminosity, and a luminosity decline over 40 days of $\Delta L_{40}\approx-0.7$ make ASASSN-23bd one of the dimmest TDEs to date and a member of the growing ``Low Luminosity and Fast'' class of TDEs., Comment: 17 pages, 13 figures, submitted to MNRAS

We present comprehensive optical observations of SN~2021gmj, a Type II supernova (SN~II) discovered within a day of explosion by the Distance Less Than 40~Mpc (DLT40) survey. Follow-up observations show that SN~2021gmj is a low-luminosity SN~II (LL~SN~II), with a peak magnitude $M_V = -15.45$ and Fe~II velocity of $\sim 1800 \ \mathrm{km} \ \mathrm{s}^{-1}$ at 50 days past explosion. Using the expanding photosphere method, we derive a distance of $17.8^{+0.6}_{-0.4}$~Mpc. From the tail of the light curve we obtain a radioactive nickel mass of $0.014 \pm 0.001$ M$_{\odot}$. The presence of circumstellar material (CSM) is suggested by the early-time light curve, early spectra, and high-velocity H$\alpha$ in absorption. Analytical shock-cooling models of the light curve cannot reproduce the fast rise, supporting the idea that the early-time emission is partially powered by the interaction of the SN ejecta and CSM. The inferred low CSM mass of 0.025 M$_{\odot}$ in our hydrodynamic-modeling light curve analysis is also consistent with our spectroscopy. We observe a broad feature near 4600 \AA, which may be high-ionization lines of C, N, or/and He~II. This feature is reproduced by radiation-hydrodynamic simulations of red supergiants with extended atmospheres. Several LL~SNe~II show similar spectral features, implying that high-density material around the progenitor may be common among them., Comment: Accepted version at ApJ

Local ionospheric density anomalies have been reported in the days prior to major earthquakes. This global study statistically investigates whether consistent ionospheric anomalies occur in the 24 hours prior to earthquakes across different regions, magnitudes, temporal and spatial scales. We match earthquake data to Total Electron Content (TEC) data from 2000-2020 at a higher resolution and cadence than previous assessed. Globally, no significant, consistent anomaly is found. Regionally, statistically significant ionospheric anomalies arise in the 12 hours prior to earthquakes with $p \leq 0.01$ following Wilcoxon tests. For the Japanese region we find a median negative ionospheric anomaly of around 0.5 TECU between 3 and 8 hours before earthquakes. For the South American region, the median TEC is enhanced by up to ~ 2 TECU, between 7 and 10 hours before an event. We show that the results are robust to different definitions of the ''local'' region and earthquake magnitude. This demonstrates the promise of monitoring the ionosphere as part of a multimodal earthquake forecasting system., Comment: 12 pages, 4 figures. Presented at AGU fall meeting 2022 (https://agu.confex.com/agu/fm22/meetingapp.cgi/Paper/1142329)

We demonstrate the importance of radio selection in probing heavily obscured galaxy populations. We combine Evolutionary Map of the Universe (EMU) Early Science data in the Galaxy and Mass Assembly (GAMA) G23 field with the GAMA data, providing optical photometry and spectral line measurements, together with Wide-field Infrared Survey Explorer (WISE) infrared (IR) photometry, providing IR luminosities and colours. We investigate the degree of obscuration in star forming galaxies, based on the Balmer decrement (BD), and explore how this trend varies, over a redshift range of 0

All very massive early-type galaxies contain supermassive blackholes but are these blackholes all sufficiently active to produce detectable radio continuum sources? We have used the 887.5~MHz Rapid ASKAP Continuum Survey DR1 to measure the radio emission from morphological early-type galaxies brighter than $K_S=9.5$ selected from the 2MASS Redshift Survey, HyperLEDA and RC3. In line with previous studies, we find median radio power increases with infrared luminosity, with $P_{1.4} \propto L_K^{2.2}$, although the scatter about this relation spans several orders of magnitude. All 40 of the $M_K<-25.7$ early-type galaxies in our sample have measured radio flux densities that are more than $2\sigma$ above the background noise, with $1.4~{\rm GHz}$ radio powers spanning $\sim 3 \times 10^{20}$ to $\sim 3\times 10^{25}~{\rm W~Hz^{-1}}$. Cross matching our sample with integral field spectroscopy of early-type galaxies reveals that the most powerful radio sources preferentially reside in galaxies with relatively low angular momentum (i.e. slow rotators). While the infrared colours of most galaxies in our early-type sample are consistent with passive galaxies with negligible star formation and the radio emission produced by active galactic nuclei or AGN remnants, very low levels of star formation could power the weakest radio sources with little effect on many other star formation rate tracers., Comment: Accepted for publication in the Publications of the Astronomical Society of Australia. 9 pages, 7 figures, 1 table

AM CVn-type systems are ultracompact, helium-accreting binary systems which are evolutionarily linked to the progenitors of thermonuclear supernovae and are expected to be strong Galactic sources of gravitational waves detectable to upcoming space-based interferometers. AM CVn binaries with orbital periods $\lesssim$ 20--23 min exist in a constant high state with a permanently ionised accretion disc. We present the discovery of TIC 378898110, a bright ($G=14.3$ mag), nearby ($309.3 \pm 1.8$ pc), high-state AM CVn binary discovered in TESS two-minute-cadence photometry. At optical wavelengths this is the third-brightest AM CVn binary known. The photometry of the system shows a 23.07172(6) min periodicity, which is likely to be the `superhump' period and implies an orbital period in the range 22--23 min. There is no detectable spectroscopic variability. The system underwent an unusual, year-long brightening event during which the dominant photometric period changed to a shorter period (constrained to $20.5 \pm 2.0$ min), which we suggest may be evidence for the onset of disc-edge eclipses. The estimated mass transfer rate, $\log (\dot{M} / \mathrm{M_\odot} \mathrm{yr}^{-1}) = -6.8 \pm 1.0$, is unusually high and may suggest a high-mass or thermally inflated donor. The binary is detected as an X-ray source, with a flux of $9.2 ^{+4.2}_{-1.8} \times 10^{-13}$ erg cm$^{-2}$ s$^{-1}$ in the 0.3--10 keV range. TIC 378898110 is the shortest-period binary system discovered with TESS, and its large predicted gravitational-wave amplitude makes it a compelling verification binary for future space-based gravitational wave detectors., Comment: 15 pages, 14 figures. Accepted to MNRAS

We study parallel fault-tolerant quantum computing for families of homological quantum low-density parity-check (LDPC) codes defined on 3-manifolds with constant or almost-constant encoding rate. We derive generic formula for a transversal $T$ gate of color codes on general 3-manifolds, which acts as collective non-Clifford logical CCZ gates on any triplet of logical qubits with their logical-$X$ membranes having a $\mathbb{Z}_2$ triple intersection at a single point. The triple intersection number is a topological invariant, which also arises in the path integral of the emergent higher symmetry operator in a topological quantum field theory: the $\mathbb{Z}_2^3$ gauge theory. Moreover, the transversal $S$ gate of the color code corresponds to a higher-form symmetry supported on a codimension-1 submanifold, giving rise to exponentially many addressable and parallelizable logical CZ gates. We have developed a generic formalism to compute the triple intersection invariants for 3-manifolds and also study the scaling of the Betti number and systoles with volume for various 3-manifolds, which translates to the encoding rate and distance. We further develop three types of LDPC codes supporting such logical gates: (1) A quasi-hyperbolic code from the product of 2D hyperbolic surface and a circle, with almost-constant rate $k/n=O(1/\log(n))$ and $O(\log(n))$ distance; (2) A homological fibre bundle code with $O(1/\log^{\frac{1}{2}}(n))$ rate and $O(\log^{\frac{1}{2}}(n))$ distance; (3) A specific family of 3D hyperbolic codes: the Torelli mapping torus code, constructed from mapping tori of a pseudo-Anosov element in the Torelli subgroup, which has constant rate while the distance scaling is currently unknown. We then show a generic constant-overhead scheme for applying a parallelizable universal gate set with the aid of logical-$X$ measurements., Comment: 40 pages, 32 figures. In the updated version v2, we have simplified the TQFT derivation of the logical gates via an operator-valued cochain formalism in Sec. III, which also gives rise to the explicit construction of constant-depth circuits corresponding to logical CCZ and CZ gates in three copies of identical toric codes defined on arbitrary 3-manifolds