Locally authored papers of the past 5 days

We use the stellar intensity interferometry system implemented with the Very Energetic Radiation Imaging Telescope Array System (VERITAS) at Fred Lawrence Whipple Observatory (FLWO) as a light collector to obtain measurements of the rapid rotator star $\gamma$ Cassiopeiae, at a wavelength of 416 nm. Using data from baselines sampling different position angles, we extract the size, oblateness, and projected orientation of the photosphere. Fitting the data with a uniform ellipse model yields a minor-axis angular diameter of $0.43\pm0.02$ mas, a major-to-minor-radius ratio of $1.28\pm0.04$, and a position angle of $116^\circ\pm5^\circ$ for the axis of rotation. A rapidly-rotating stellar atmosphere model that includes limb and gravity darkening describes the data well with a fitted angular diameter of $0.604^{+0.041}_{-0.034}$ mas corresponding to an equatorial radius of 10.9$^{+0.8}_{-0.6}~R_\odot$, a rotational velocity with a $1~\sigma$ lower limit at $97.7\%$ that of breakup velocity, and a position angle of $114.7^{+6.4}_{-5.7}$ degrees. These parameters are consistent with H$\alpha$ line spectroscopy and infrared-wavelength Michelson interferometric measurements of the star's decretion disk. This is the first measurement of an oblate photosphere using intensity interferometry.

In the era of JWST, with its unprecedented sensitivity and spectral resolution, infrared spectral surveys have revealed a rich inventory of ices, including complex organic molecules (COMs), in young stellar objects (YSOs). However, robust methods to decompose and quantify these absorption features particularly across broad spectral ranges, are still under investigation. We present INDRA (Ice-fitting with NNLS-based Decomposition and Retrieval Algorithm), a fully Python-based tool that performs continuum and silicate removal, global ice fitting using Weighted Non-Negative Least Squares (NNLS), and estimates column densities and statistical significance. We apply INDRA to NGC 1333 IRAS 2A, a target from the JWST Observations of Young protoStars (JOYS+) program previously studied using local fitting. We derive optical depths via polynomial continuum subtraction and remove silicate absorption using a synthetic model, isolating ice features for global MIRI fitting. Our results are consistent with previous local fits, confirming simple species and COMs, and expand the inventory by identifying additional absorption features from CO2 and NH4+. We also propose the presence of organic refractories contributing up to 9.6% in the spectral region of 5-8 microns among the various ice components, whose inclusion significantly improves the global spectral fitting. These broad absorption features, extending across 5.5-11 microns, are likely produced by large, complex molecules containing carbonyl (C=O), hydroxyl (O-H), amine (N-H), and C-H bending modes. Our expanded inventory, now incorporating these organic residues, offers new insights into the chemical evolution of ices in star-forming regions and highlights the importance of global spectral fitting in constraining ice compositions.

We present updated non-adiabatic and inhomogeneous evolution models for Uranus and Neptune, employing an interior composition of methane, ammonia, water, and rocks. Following formation trends of the gas giants, Uranus and Neptune formation models are applied, where both planets begin with layers stable to convection. Both planets are subject to convective mixing throughout their evolution. Uranus undergoes modest convective mixing, preserving much of its primordial internal heat. In contrast, Neptune's interior undergoes extensive mixing, homogenization, and adiabatic cooling of the outer 40\% of its envelope. The subsequent release of internal energy in Neptune, driven by the convective instability of its primordial outer compositional gradient, accounts for its higher luminosity relative to Uranus. Thus, the observed luminosity differences between Uranus and Neptune could be primarily dictated by the convective stability of their outer envelopes. The extensive convective mixing in Neptune leads to a higher metallicity in its outer region compared to Uranus, a feature seen in atmospheric measurements and shown in past interior models of Neptune. Due to Neptune's more pronounced cooling, our models predict favorable conditions for hydrogen-water immiscibility in its envelope.

(abridged) In a series of publications, we describe a comprehensive comparison of Event Horizon Telescope (EHT) data with theoretical models of Sgr A* and M87*. Here, we report on improvements made to our observational data reduction pipeline and present the generation of observables derived from the EHT models. We make use of ray-traced GRMHD simulations that are based on different black hole spacetime metrics and accretion physics parameters. These broad classes of models provide a good representation of the primary targets observed by the EHT. To generate realistic synthetic data from our models, we took the signal path as well as the calibration process, and thereby the aforementioned improvements, into account. We could thus produce synthetic visibilities akin to calibrated EHT data and identify salient features for the discrimination of model parameters. We have produced a library consisting of an unparalleled 962,000 synthetic Sgr A* and M87* datasets. In terms of baseline coverage and noise properties, the library encompasses 2017 EHT measurements as well as future observations with an extended telescope array. We differentiate between robust visibility data products related to model features and data products that are strongly affected by data corruption effects. Parameter inference is mostly limited by intrinsic model variability, which highlights the importance of long-term monitoring observations with the EHT. In later papers in this series, we will show how a Bayesian neural network trained on our synthetic data is capable of dealing with the model variability and extracting physical parameters from EHT observations. With our calibration improvements, our newly reduced EHT datasets have a considerably better quality compared to previously analyzed data.

(abridged) In this second paper in our publication series, we present the open-source Zingularity framework for parameter inference with deep Bayesian artificial neural networks. We carried out out supervised learning with synthetic millimeter very long baseline interferometry observations of the EHT. Our ground-truth models are based on GRMHD simulations of Sgr A* and M87* on horizon scales. We investigated how well Zingularity neural networks are able to infer key model parameters from EHT observations, such as the black hole spin and the magnetic state of the accretion disk, when uncertainties in the data are accurately taken into account. Zingularity makes use of the TensorFlow Probability library and is able to handle large amounts of data with a combination of the efficient TFRecord data format plus the Horovod framework. Our approach is the first analysis of EHT data with Bayesian neural networks, where an unprecedented training data size, under consideration of a closely modeled EHT signal path, and the full information content of the observational data are used. Zingularity infers parameters based on salient features in the data and is containerized. Through parameter surveys and dedicated validation tests, we identified neural network architectures, that are robust against internal stochastic processes and unaffected by noise in the observational and model data. We give examples of how different data properties affect the network training. We show how the Bayesian nature of our networks gives trustworthy uncertainties and uncovers failure modes for uncharacterizable data. It is easy to achieve low validation errors during training on synthetic data with neural networks, particularly when the forward modeling is too simplified. Through careful studies, we demonstrate that our trained networks can generalize well so that reliable results can be obtained from observational data.

(abridged) In the first two papers of this publication series, we present a comprehensive library of synthetic EHT observations and used this library to train and validate Bayesian neural networks for the parameter inference of accreting supermassive black hole systems. The considered models are ray-traced GRMHD simulations of Sgr A* and M87*. In this work, we infer the best-fitting accretion and black hole parameters from 2017 EHT data and predict improvements that will come with future upgrades of the array. Compared to previous EHT analyses, we considered a substantially larger synthetic data library and the most complete set of information from the observational data. We made use of the Bayesian nature of the trained neural networks and apply bootstrapping of known systematics in the observational data to obtain parameter posteriors. Within a wide GRMHD parameter space, we find M87* to be best described by a spin between 0.5 and 0.94 with a retrograde MAD accretion flow and strong synchrotron emission from the jet. Sgr A* has a high spin of $\sim$ 0.8 $-$ 0.9 and a prograde accretion flow beyond the standard MAD/SANE models with a comparatively weak jet emission, seen at a $\sim$ 20$^\circ$ $-$ 40$^\circ$ inclination and $\sim$ 106$^\circ$ $-$ 137$^\circ$ position angle. While previous EHT analyses could rule out specific regions in the model parameter space considered here, we are able to obtain narrow parameter posteriors with our Zingularity framework without being impacted by the unknown foreground Faraday screens and data calibration biases. We further demonstrate that the AMT extension to the EHT will reduce parameter inference errors by a factor of three for non-Kerr models, enabling more robust tests of general relativity. It will be instructive to produce new GRMHD models with the inferred interpolated parameters to study their accretion rate plus jet power.

Supermassive black hole feedback is the currently favoured mechanism to regulate the star formation rate of galaxies and prevent the formation of ultra-massive galaxies ($M_\star>10^{12}M_\odot$). However, the mechanism through which the outflowing energy is transferred to the surrounding medium strongly varies from one galaxy evolution model to another, such that a unified model for AGN feedback does not currently exist. The hot atmospheres of galaxy groups are highly sensitive laboratories of the feedback process, as the injected black hole energy is comparable to the binding energy of halo gas particles. Here we report multi-wavelength observations of the fossil galaxy group SDSSTG 4436. The hot atmosphere of this system exhibits a highly relaxed morphology centred on the giant elliptical galaxy NGC~3298. The X-ray emission from the system features a compact core ($<$10 kpc) and a steep increase in the entropy and cooling time of the gas, with the cooling time reaching the age of the Universe $\sim15$ kpc from the centre of the galaxy. The observed entropy profile implies a total injected energy of $\sim1.5\times10^{61}$ ergs, which given the high level of relaxation could not have been injected by a recent merging event. Star formation in the central galaxy NGC~3298 is strongly quenched and its stellar population is very old ($\sim$10.6 Gyr). The currently detected radio jets have low power and are confined within the central compact core. All the available evidence implies that this system was affected by giant AGN outbursts which excessively heated the neighbouring gas and prevented the formation of a self-regulated feedback cycle. Our findings imply that AGN outbursts can be energetic enough to unbind gas particles and lead to the disruption of cool cores.

Recent observations at low redshift have revealed that some post-starburst galaxies retain significant molecular gas reservoirs despite low ongoing star formation rates, challenging theoretical predictions for galaxy quenching. To test whether this finding holds during the peak epoch of quenching, here we present ALMA CO(2-1) observations of five spectroscopically confirmed post-starburst galaxies at z ~ 1.4 from the HeavyMetal survey. While four galaxies are undetected in CO emission, we detect M_H2 ~ 10^9.7 Msun of molecular gas in one system. The detected system is a close pair of massive (M* = 10^(11.1-11.2) Msun) post-starburst galaxies with no clear tidal features, likely caught in the early stages of a major merger. These results suggest that mergers may be a key factor in retaining molecular gas while simultaneously suppressing star formation in quenched galaxies at high redshift, possibly by driving increased turbulence that decreases star formation efficiency. Unlike previous studies at z < 1, we find no correlation between molecular gas mass and time since quenching. This may be explained by the fact that -- despite having similar UVJ colors -- all galaxies in our sample have post-burst ages older than typical gas-rich quenched systems at low redshift. Our results highlight the importance of major mergers in shaping the cold gas content of quiescent galaxies during the peak epoch of quenching.

Quasi-periodic pulsations (QPPs) at sub-second periods are frequently detected in the time series of X-rays during stellar flares. However, such rapid pulsations are rarely reported in the hard X-ray (HXR) emission of the small solar flare. We explored the QPP patterns with fast-time variations in HXR and radio emissions produced in a small solar flare on 2025 January 19. By applying the Fast Fourier Transform, the fast-variation pulsations at a quasi-period of about 1 s are identified in the HXR channel of 20-80 keV, which were simultaneously measured by the Hard X-ray Imager and the Konus-Wind. The rapid pulsations with a same quasi-period were also detected in the radio emission at a lower frequency range of about 40-100 MHz. The restructured HXR images show that the QPP patterns mainly locate in footpoint areas that connect by hot plasma loops, and they appear in the flare impulsive phase. Our observations suggest that the fast-variation pulsations could be associated with nonthermal electrons that are periodically accelerated by the intermittent magnetic reconnection, and the 1-s period may be modulated by the coalescence instability between current-carrying loops and magnetic islands.

In the standard model of solar eruptive events, coronal mass ejections (CMEs) and flares are associated with each other through magnetic reconnection initiated by erupting flux ropes. Observations also reveal an increasing association ratio between flares and CMEs with flare intensity. However, the fundamental relationship between flares and CMEs, and that between thermal and nonthermal processes, remains unknown. Here we investigate energetic C-class flares (ECFs) -- Geostationary Operational Environmental Satellite (GOES) C-class flares with hard X-ray (HXR) emissions above 30 keV -- using observations from Advanced Space-based Solar Observatory/Hard X-ray Imager (HXI), Solar Dynamic Observatory, and GOES. Among 1331 C-class flares detected by HXI, 127 ECFs (9.5%) were identified for statistical analysis of their properties and associations with CMEs and other flare-related features. Our statistical results reveal that ECFs have relatively shorter durations and harder spectra (the mean electron power-law index is 4.65), with no significant correlation between soft X-ray flux and nonthermal parameters (e.g., HXR peak flux). Among the 127 events, 53 (42%) were associated with type III bursts, 38 (30%) with jets, at least 13 (~11%) with 360 nm brightenings, and only 5 (~4%) with CMEs. Crucially, all five CME events were narrow CMEs associated with jets. The surprising correlation between these ECFs and CMEs suggests that noneruptive or confined magnetic field configurations in these flares may favor electron acceleration, resulting in harder X-ray this http URL discuss the potential formation mechanisms and efficient electron acceleration processes in these atypical flares, providing valuable insights into nonstandard flare behavior.

We present cosmological constraints from Dark Energy Survey Year 3 (DES Y3) weak lensing data using persistent homology, a topological data analysis technique that tracks how features like clusters and voids evolve across density thresholds. For the first time, we apply spherical persistent homology to galaxy survey data through the algorithm TopoS2, which is optimized for curved-sky analyses and HEALPix compatibility. Employing a simulation-based inference framework with the Gower Street simulation suite, specifically designed to mimic DES Y3 data properties, we extract topological summary statistics from convergence maps across multiple smoothing scales and redshift bins. After neural network compression of these statistics, we estimate the likelihood function and validate our analysis against baryonic feedback effects, finding minimal biases (under $0.3\sigma$) in the $\Omega_\mathrm{m}-S_8$ plane. Assuming the $w$CDM model, our combined Betti numbers and second moments analysis yields $S_8 = 0.821 \pm 0.018$ and $\Omega_\mathrm{m} = 0.304\pm0.037$-constraints 70% tighter than those from cosmic shear two-point statistics in the same parameter plane. Our results demonstrate that topological methods provide a powerful and robust framework for extracting cosmological information, with our spherical methodology readily applicable to upcoming Stage IV wide-field galaxy surveys.

Recent observations have revealed slow, coherent temperature fluctuations in AGN disks that propagate both inward and outward at velocities of $\sim 0.01 - 0.1c$, a kind of variability that is distinct from reverberation (mediated by the reprocessing of light) between different regions of the disk. We investigate the origin and nature of these fluctuations using global 3D radiation-magnetohydrodynamic simulations of radiation and magnetic pressure-dominated AGN accretion disks. Disks with a significant turbulent Maxwell stress component exhibit wave-like temperature perturbations, most evident close to the midplane, whose propagation speeds exactly match the local fast magnetosonic speed and are consistent with the speeds inferred in observations. These fluctuations have amplitudes of $2 - 4\%$ in gas temperature, which are also consistent with observational constraints. Disks that are dominated by mean-field Maxwell stresses do not exhibit such waves. While waves may be present in the body of the disk, we do not find them to be present in the photosphere. Although this may in part be due to low numerical resolution in the photosphere region, we discuss the physical challenges that must be overcome for the waves to manifest there. In particular, the fact that such waves are observed implies that the disk photospheres must be magnetically dominated, since radiative damping from photon diffusion smooths out radiation pressure fluctuations. Furthermore, the gas and radiation fluctuations must be out of local thermodynamic equilibrium.

We perform two-dimensional, multi-group radiation hydrodynamic simulations to explore the observational properties of a solar-like star colliding with an accretion disk around a supermassive black hole at separation of $\sim 100$ gravitational radii. We find that the star-disk collision produces ejecta on both sides of the disk. As the ejecta expand and cool, transient flares arise, reaching peak bolometric luminosity of up to $L\gtrsim10^{43}\rm erg~s^{-1}$. We estimate that the typical light curve rises and decays on an hour timescale. The spectral energy distribution (SED) peaks in $20-50$eV. The optical depth in soft X-rays is lower than the frequency-integrated optical depth, yielding $100$eV-$1$KeV luminosity $\nu L_{\nu}\gtrsim10^{42}\rm erg~s^{-1}$. The ejecta aligned with the star's incident direction shows breakout emission, leading to asymmetric SED evolution of the two ejecta. The SED evolution is roughly consistent with those seen in short-period quasi-periodic eruptions (QPEs), which have eruption duration ranging from sub-hour to hours, but the ejecta cooling emission alone may not be sufficient to explain the longer duration flares. Increasing incident velocity generally produces a brighter and harder flare. A larger disk scale height prolongs the breakout emission but leads to a somewhat softer SED. A higher disk surface density can lead to higher ejecta temperature, reducing bound-free opacity and increasing luminosity. When lowering the disk surface density, we find that the ejecta becomes optically thin when the scattering optical depth across disk is at the order of $\tau_{\rm disk}\sim200$, and the ejecta disappear when $\tau_{\rm disk}\sim10$.

We present the Local Group L-Band Survey (LGLBS), a Karl G. Jansky Very Large Array (VLA) survey producing the highest quality 21-cm and 1-2 GHz radio continuum images to date for the six VLA-accessible, star-forming, Local Group galaxies. Leveraging the VLA's spectral multiplexing power, we simultaneously survey the 21-cm line at high 0.4 km/s velocity resolution, the 1-2 GHz polarized continuum, and four OH lines. For the massive spiral M31, the dwarf spiral M33, and the dwarf irregular galaxies NGC6822, IC10, IC1613, and the Wolf-Lundmark-Melotte Galaxy (WLM), we use all four VLA configurations and the Green Bank Telescope to reach angular resolutions of $< 5''$ ($10{-}20$~pc) for the 21-cm line with $<10^{20}$~cm$^{-2}$ column density sensitivity, and even sharper views ($< 2''$; $5{-}10$~pc) of the continuum. Targeting these nearby galaxies ($D\lesssim1$ Mpc) reveals a sharp, resolved view of the atomic gas, including 21-cm absorption, and continuum emission from supernova remnants and HII regions. These datasets can be used to test theories of the abundance and formation of cold clouds, the driving and dissipation of interstellar turbulence, and the impact of feedback from massive stars and supernovae. Here, we describe the survey design and execution, scientific motivation, data processing, and quality assurance. We provide a first look at and publicly release the wide-field 21-cm HI data products for M31, M33, and four dwarf irregular targets in the survey, which represent some of the highest physical resolution 21-cm observations of any external galaxies beyond the LMC and SMC.

Recent JWST mid-infrared (mid-IR) images, tracing polycyclic aromatic hydrocarbons (PAHs) and dust continuum emission, provide detailed views of the interstellar medium (ISM) in nearby galaxies. Leveraging PHANGS-JWST Cycle 1 and PHANGS-MUSE data, we measure the PAH and dust continuum emission lifetimes of gas clouds across 17 nearby star-forming galaxies by analyzing the relative spatial distributions of mid-IR (7.7-11.3$\mu$m) and H$\alpha$ emission at various scales. We find that the mid-IR emitting time-scale of gas clouds in galaxy disks (excluding centers) ranges from 10 to 30Myr. After star formation is detected in H$\alpha$, mid-IR emission persists for 3-7Myr during the stellar feedback phase, covering 70-80% of the H$\alpha$ emission. This significant overlap is due to intense radiation from star-forming regions, illuminating the surrounding PAHs and dust grains. In most galaxies, the mid-IR time-scale closely matches the molecular cloud lifetime measured with CO. Although mid-IR emission is complex as influenced by ISM distribution, radiation, and abundances of dust and PAHs, the similarity between the two time-scales suggests that once gas clouds form with compact mid-IR emission, they quickly provide sufficient shielding for stable CO formation. This is likely due to our focus on molecular gas-rich regions of galaxies with near-solar metallicity. Finally, we find that the mid-IR emitting time-scale is longer in galaxies with well-defined HII regions and less structured backgrounds, allowing photons to more efficiently heat the ambient ISM surrounding the HII regions, rather than contributing to diffuse emission. This suggests that the shape of the ISM also influences mid-IR emission.

Quasi-periodic eruption (QPE) sources in galactic nuclei are often associated with a stellar object orbiting a supermassive black hole with hours-days period, brought in as an extreme mass-ratio inspiral (EMRI). In the presence of an accretion disk, repeated star-disk collisions lead to ablation of a small fraction of the stellar mass during each disk passage. We analytically follow the evolution of the stellar debris as it is tidally stretched outside the EMRI's Hill sphere, forming an elongated, dilute stream, that subsequently collides with the disk, half an orbit after the previous star-disk encounter. At sufficiently long orbital periods ($\gtrsim 12$ hr), the stream is too dilute to penetrate the disk, and is instead strongly shocked and deflected at its surface through a reverse shock. We obtain the resulting emission and explore implications for QPE observations. Due to their low optical depth and prolonged interaction time, radiation from the shocked streams typically dominates over that from shocked disk gas directly impacted by the star or by ejecta confined within its Hill sphere, as was first proposed by Yao et al. 2025. We find that: (1) QPE flare durations reflect the stream-disk collision timescale; (2) Flare luminosities of $10^{42-43}$ erg/s, consistent with observed QPEs, are robustly produced; (3) Soft X-ray flares with temperatures of ${\sim}$100 eV arise when the stream mass is sufficient to sustain a radiation mediated shock at the collision interface. Higher mass streams yield softer flares, typically outshone by the disk, while lower mass streams result in collisionless shocks, which likely produce fainter and harder flares. We discuss observational implications of the temporal evolution of the underlying disk, assuming it is the remnant of a prior tidal disruption event in the same galaxy.

Beyond LISA, proposed space-based gravitational wave (GW) missions aim to explore the sub-millihertz to microhertz frequency band, with one key objective being the detection of massive binary black hole (MBBH) mergers across cosmic distances. In this work, we investigate the detection and localization capabilities of future sub-mHz GW observatories for MBBH coalescences. Including the full galactic foreground noise, we find that signal-to-noise ratios (SNRs) can reach several thousand across a wide range of redshifts. We evaluate three representative orbital configurations--non-precessing and precessing with different inclination angles--and analyze their localization performance for various MBBH populations. In the non-precessing case, a two-hemisphere degeneracy arises when only the dominant (2,2) mode is considered, which is effectively resolved by including higher-order modes. These modes contribute to a more uniform performance across all configurations, thereby mitigating the prior advantage of precessing mission orbits. Sub-mHz missions operating in the [10 $\mu$Hz, 10 mHz] band partially overlap with LISA's range but provide enhanced sensitivity to lower-frequency GWs due to their longer interferometric baselines. This results in significantly improved localization of high-mass MBBHs, enhancing the prospects for multi-messenger astronomy and precision cosmology. Moreover, the high SNRs attainable with sub-mHz detectors could enable stringent tests of general relativity and alternative theories of gravity.

The ALMA survey of Gas Evolution in PROtoplanetary disks (AGE-PRO) Large Program aims to trace the evolution of gas disk mass and size throughout the lifetime of protoplanetary disks. This paper presents Band-6 ALMA observations of 10 embedded (Class I and Flat Spectrum) sources in the Ophiuchus molecular cloud, with spectral types ranging from M3 to K6 stars, which serve as the evolutionary starting point in the AGE-PRO sample. While we find 4 nearly edge on disks (>70 deg.), and 3 highly inclined disks (>60 deg.) in our sample, we show that, as a population, embedded disks in Ophiuchus are not significantly contaminated by more evolved, but highly inclined sources. We derived dust disk masses from the Band 6 continuum and estimated gas disk masses from the C18O and C17O lines. The mass estimates from the C17O line are slightly higher, suggesting C18O emission might be partially optically thick. While the 12CO and 13CO lines are severely contaminated by extended emission and self-absorption, the C18O and C17O lines allowed us to trace the radial extent of the gaseous disks. From these measurements, we found that the C18O and C17O fluxes correlate well with each other and with the continuum fluxes. Furthermore, the C18O and C17O lines present a larger radial extension than disk dust sizes by factors ranging from 1.5 to 2.5, as it is found for Class II disks using the radial extension of the 12CO. In addition, we have detected outflows in three disks from 12CO observations.

Variability of millimetre wavelength continuum emission from Class II protoplanetary disks is extremely rare, and when detected it is usually interpreted as originating from non-thermal emission mechanisms that relate to the host star itself rather than its disk. During observations made as part of the AGE-PRO ALMA Large program, significant variability in the brightness of the 2MASS J16202863-2442087 system was detected between individual executions. We report the observed properties of the variability detected at millimetre wavelengths and investigate potential driving mechanisms. To investigate the nature of the variability we construct a light curve from the continuum observations and analyse imaged constructed from both flaring and quiescent emission. We characterise the dust disk around the star through analysis in the image and visibility plane, and carry out kinematic analysis of the CO(2-1) emission from the gas disk. The continuum flux decays by a factor of 8 in less than an hour, and by a factor of 13 within 8 days. The peak brightness coincides with an expected brightness maximum extrapolated from the periodicity of previously observed optical variability. The flare is most likely the product of synchrotron emission in the close vicinity of the star. The nature of the millimetre flare closely resembles those detected in very close binary systems, and may be due to the interaction of magnetic fields in an as yet undetected binary. Alternatively if the central host is a single-star object, the flare may be due to the interaction of magnetic field loops at the stellar surface or a strong accretion burst.

A handful of stars are known to host both an inner system of multiple transiting planets and an outer giant planet. These systems all feature a prominent gap between the orbits of two of the transiting planets, distinguishing them from typical multiplanet systems with more uniform orbital spacings. The reason for the association between inner gaps and outer giants is unknown. In this paper, we assess whether undiscovered planets might occupy these gaps in systems with outer giants. For each of the four relevant systems - Kepler-48, Kepler-65, Kepler-90, and Kepler-139 - we found that a $\sim 2 - 20 M_\oplus$ planet could reside in the gap without inducing dynamical instability. However, in each case, the gravitational influence of the outer giant planet is insufficient to tilt the orbit of the hypothetical planet by enough to prevent transits, ruling out a proposed theory for the observed gap-giant association. The gaps might instead contain smaller, undetectable planets ($ \lesssim 0.5 - 1\,R_\oplus$), or be entirely devoid of planets.