Showing 1–50 of 1,038 results for author: Anderson, D

Filtrations and recursions for Schubert modules

Authors: David Anderson

Abstract: Revisiting Kraśkiewicz and Pragacz's construction of Schubert modules, we provide a new proof that their characters are equal to Schubert polynomials. The main innovation is a representation-theoretic interpretation of a recurrence relation for Schubert polynomials recently discovered by Nadeau, Spink, and Tewari. Along the way, we review several related constructions, and show that the Nadeau-Spi… ▽ More Revisiting Kraśkiewicz and Pragacz's construction of Schubert modules, we provide a new proof that their characters are equal to Schubert polynomials. The main innovation is a representation-theoretic interpretation of a recurrence relation for Schubert polynomials recently discovered by Nadeau, Spink, and Tewari. Along the way, we review several related constructions, and show that the Nadeau-Spink-Tewari recursion determines the characters of flagged Schur modules coming from a broader class of translucent diagrams. △ Less

Submitted 29 August, 2024; originally announced August 2024.

Comments: 37 pages, many figures

Time-Evolution Images of the Hypergiant RW Cephei During the Re-brightening Phase Following the Great Dimming

Authors: Narsireddy Anugu, Douglas R. Gies, Rachael M. Roettenbacher, John D. Monnier, Miguel Montargés, Antoine Mérand, Fabien Baron, Gail H. Schaefer, Katherine A. Shepard, Stefan Kraus, Matthew D. Anderson, Isabelle Codron, Tyler Gardner, Mayra Gutierrez, Rainer Köhler, Karolina Kubiak, Cyprien Lanthermann, Olli Majoinen, Nicholas J. Scott, Wolfgang Vollmann

Abstract: Stars with initial masses larger than 8 solar masses undergo substantial mass loss through mechanisms that remain elusive. Unraveling the origins of this mass loss is important for comprehending the evolutionary path of these stars, the type of supernova explosion and whether they become neutron stars or black hole remnants. In 2022 December, RW Cep experienced the Great Dimming in its visible bri… ▽ More Stars with initial masses larger than 8 solar masses undergo substantial mass loss through mechanisms that remain elusive. Unraveling the origins of this mass loss is important for comprehending the evolutionary path of these stars, the type of supernova explosion and whether they become neutron stars or black hole remnants. In 2022 December, RW Cep experienced the Great Dimming in its visible brightness, presenting a unique opportunity to understand mass loss mechanisms. Our previous observations of RW Cep from the CHARA Array, taken during the dimming phase, show a compelling asymmetry in the star images, with a darker zone on the west side of the star indicating presence of dust in front of the star in our line of sight. Here, we present multi-epoch observations from CHARA while the star re-brightened in 2023. We created images using three image reconstruction methods and an analytical model fit. Comparisons of images acquired during the dimming and re-brightening phases reveal remarkable differences. Specifically, the west side of RW Cep, initially obscured during the dimming phase, reappeared during the subsequent re-brightening phase and the measured angular diameter became larger by 8%. We also observed image changes from epoch to epoch while the star is brightening indicating the time evolution of dust in front of the star. We suggest that the dimming of RW Cep was a result from a recent surface mass ejection event, generating a dust cloud that partially obstructed the stellar photosphere. △ Less

Submitted 21 August, 2024; originally announced August 2024.

Comments: Accepted by The Astrophysical Journal Letters (ApJL)

Parametric Sensitivity Analysis for Models of Reaction Networks within Interacting Compartments

Authors: David F. Anderson, Aidan S. Howells

Abstract: Models of reaction networks within interacting compartments (RNIC) are a generalization of stochastic reaction networks. It is most natural to think of the interacting compartments as ``cells'' that can appear, degrade, split, and even merge, with each cell containing an evolving copy of the underlying stochastic reaction network. Such models have a number of parameters, including those associated… ▽ More Models of reaction networks within interacting compartments (RNIC) are a generalization of stochastic reaction networks. It is most natural to think of the interacting compartments as ``cells'' that can appear, degrade, split, and even merge, with each cell containing an evolving copy of the underlying stochastic reaction network. Such models have a number of parameters, including those associated with the internal chemical model and those associated with the compartment interactions, and it is natural to want efficient computational methods for the numerical estimation of sensitivities of model statistics with respect to these parameters. Motivated by the extensive work on computational methods for parametric sensitivity analysis in the context of stochastic reaction networks over the past few decades, we provide a number of methods in the RNIC setting. Provided methods include the (unbiased) Girsanov transformation method (also called the Likelihood Ratio method) and a number of coupling methods for the implementation of finite differences. We provide several numerical examples and conclude that the method associated with the ``Split Coupling'' provides the most efficient algorithm. This finding is in line with the conclusions from the work related to sensitivity analysis of standard stochastic reaction networks. We have made all of the Matlab code used to implement the various methods freely available for download. △ Less

Submitted 17 August, 2024; originally announced August 2024.

Comments: 31 pages, a number of images

MSC Class: 92C40; 60J28; 65C05

TOI-2490b- The most eccentric brown dwarf transiting in the brown dwarf desert

Authors: Beth A. Henderson, Sarah L. Casewell, Andrés Jordán, Rafael Brahm, Thomas Henning, Samuel Gill, L. C. Mayorga, Carl Ziegler, Keivan G. Stassun, Michael R. Goad, Jack Acton, Douglas R. Alves, David R. Anderson, Ioannis Apergis, David J. Armstrong, Daniel Bayliss, Matthew R. Burleigh, Diana Dragomir, Edward Gillen, Maximilian N. Günther, Christina Hedges, Katharine M. Hesse, Melissa J. Hobson, James S. Jenkins, Jon M. Jenkins , et al. (18 additional authors not shown)

Abstract: We report the discovery of the most eccentric transiting brown dwarf in the brown dwarf desert, TOI02490b. The brown dwarf desert is the lack of brown dwarfs around main sequence stars within $\sim3$~AU and is thought to be caused by differences in formation mechanisms between a star and planet. To date, only $\sim40$ transiting brown dwarfs have been confirmed. \systemt is a $73.6\pm2.4$ \mjupnos… ▽ More We report the discovery of the most eccentric transiting brown dwarf in the brown dwarf desert, TOI02490b. The brown dwarf desert is the lack of brown dwarfs around main sequence stars within $\sim3$~AU and is thought to be caused by differences in formation mechanisms between a star and planet. To date, only $\sim40$ transiting brown dwarfs have been confirmed. \systemt is a $73.6\pm2.4$ \mjupnospace, $1.00\pm0.02$ \rjup brown dwarf orbiting a $1.004_{-0.022}^{+0.031}$ \msunnospace, $1.105_{-0.012}^{+0.012}$ \rsun sun-like star on a 60.33~d orbit with an eccentricity of $0.77989\pm0.00049$. The discovery was detected within \tess sectors 5 (30 minute cadence) and 32 (2 minute and 20 second cadence). It was then confirmed with 31 radial velocity measurements with \feros by the WINE collaboration and photometric observations with the Next Generation Transit Survey. Stellar modelling of the host star estimates an age of $\sim8$~Gyr, which is supported by estimations from kinematics likely placing the object within the thin disc. However, this is not consistent with model brown dwarf isochrones for the system age suggesting an inflated radius. Only one other transiting brown dwarf with an eccentricity higher than 0.6 is currently known in the brown dwarf desert. Demographic studies of brown dwarfs have suggested such high eccentricity is indicative of stellar formation mechanisms. △ Less

Submitted 8 August, 2024; originally announced August 2024.

Comments: Accepted for publication in MNRAS, 18 pages, 14 figures

Modeling Latent Neural Dynamics with Gaussian Process Switching Linear Dynamical Systems

Authors: Amber Hu, David Zoltowski, Aditya Nair, David Anderson, Lea Duncker, Scott Linderman

Abstract: Understanding how the collective activity of neural populations relates to computation and ultimately behavior is a key goal in neuroscience. To this end, statistical methods which describe high-dimensional neural time series in terms of low-dimensional latent dynamics have played a fundamental role in characterizing neural systems. Yet, what constitutes a successful method involves two opposing c… ▽ More Understanding how the collective activity of neural populations relates to computation and ultimately behavior is a key goal in neuroscience. To this end, statistical methods which describe high-dimensional neural time series in terms of low-dimensional latent dynamics have played a fundamental role in characterizing neural systems. Yet, what constitutes a successful method involves two opposing criteria: (1) methods should be expressive enough to capture complex nonlinear dynamics, and (2) they should maintain a notion of interpretability often only warranted by simpler linear models. In this paper, we develop an approach that balances these two objectives: the Gaussian Process Switching Linear Dynamical System (gpSLDS). Our method builds on previous work modeling the latent state evolution via a stochastic differential equation whose nonlinear dynamics are described by a Gaussian process (GP-SDEs). We propose a novel kernel function which enforces smoothly interpolated locally linear dynamics, and therefore expresses flexible -- yet interpretable -- dynamics akin to those of recurrent switching linear dynamical systems (rSLDS). Our approach resolves key limitations of the rSLDS such as artifactual oscillations in dynamics near discrete state boundaries, while also providing posterior uncertainty estimates of the dynamics. To fit our models, we leverage a modified learning objective which improves the estimation accuracy of kernel hyperparameters compared to previous GP-SDE fitting approaches. We apply our method to synthetic data and data recorded in two neuroscience experiments and demonstrate favorable performance in comparison to the rSLDS. △ Less

Submitted 19 July, 2024; originally announced August 2024.

CHARA Near-Infrared Imaging of the Yellow Hypergiant Star $ρ$ Cassiopeiae: Convection Cells and Circumstellar Envelope

Authors: Narsireddy Anugu, Fabien Baron, John D. Monnier, Douglas R. Gies, Rachael M. Roettenbacher, Gail H. Schaefer, Miguel Montargès, Stefan Kraus, Jean-Baptiste Le Bouquin, Matthew D. Anderson, Theo ten Brummelaar, Isabelle Codron, Christopher D. Farrington, Tyler Gardner, Mayra Gutierrez, Rainer Köhler, Cyprien Lanthermann, Ryan Norris, Nicholas J. Scott, Benjamin R. Setterholm, Norman L. Vargas

Abstract: Massive evolved stars such as red supergiants and hypergiants are potential progenitors of Type II supernovae, and they are known for ejecting substantial amounts of matter, up to half their initial mass, during their final evolutionary phases. The rate and mechanism of this mass loss play a crucial role in determining their ultimate fate and the likelihood of their progression to supernovae. Howe… ▽ More Massive evolved stars such as red supergiants and hypergiants are potential progenitors of Type II supernovae, and they are known for ejecting substantial amounts of matter, up to half their initial mass, during their final evolutionary phases. The rate and mechanism of this mass loss play a crucial role in determining their ultimate fate and the likelihood of their progression to supernovae. However, the exact mechanisms driving this mass ejection have long been a subject of research. Recent observations, such as the Great Dimming of Betelgeuse, have suggested that the activity of large convective cells, combined with pulsation, could be a plausible explanation for such mass loss events. In this context, we conducted interferometric observations of the famous yellow hypergiant, $ρ$ Cassiopeiae using the CHARA Array in H and K-band wavelengths. $ρ$ Cas is well known for its recurrent eruptions, characterized by periods of visual dimming ($\sim$1.5-2 mag) followed by recovery. From our observations, we derived the diameter of the limb-darkened disk and found that this star has a radius of $1.04\pm0.01$ milliarcseconds (mas), or $564 - 700 R_\odot$. We performed image reconstructions with three different image reconstruction software packages, and they unveiled the presence of giant hot and cold spots on the stellar surface. We interpret these prominent hot spots as giant convection cells, suggesting a possible connection to mass ejections from the star's envelope. Furthermore, we detected spectral CO emission lines in the K-band ($λ=2.31-2.38 μ$m), and the image reconstructions in these spectral lines revealed an extended circumstellar envelope with a radius of $1.45\pm0.10$ mas. △ Less

Submitted 7 August, 2024; v1 submitted 5 August, 2024; originally announced August 2024.

Comments: 21 pages, 13 figures, 4 tables, Accepted by ApJ

A Benchmark JWST Near-Infrared Spectrum for the Exoplanet WASP-39b

Authors: A. L. Carter, E. M. May, N. Espinoza, L. Welbanks, E. Ahrer, L. Alderson, R. Brahm, A. D. Feinstein, D. Grant, M. Line, G. Morello, R. O'Steen, M. Radica, Z. Rustamkulov, K. B. Stevenson, J. D. Turner, M. K. Alam, D. R. Anderson, N. M. Batalha, M. P. Battley, D. Bayliss, J. L. Bean, B. Benneke, Z. K. Berta-Thompson, J. Brande , et al. (55 additional authors not shown)

Abstract: Observing exoplanets through transmission spectroscopy supplies detailed information on their atmospheric composition, physics, and chemistry. Prior to JWST, these observations were limited to a narrow wavelength range across the near-ultraviolet to near-infrared, alongside broadband photometry at longer wavelengths. To understand more complex properties of exoplanet atmospheres, improved waveleng… ▽ More Observing exoplanets through transmission spectroscopy supplies detailed information on their atmospheric composition, physics, and chemistry. Prior to JWST, these observations were limited to a narrow wavelength range across the near-ultraviolet to near-infrared, alongside broadband photometry at longer wavelengths. To understand more complex properties of exoplanet atmospheres, improved wavelength coverage and resolution are necessary to robustly quantify the influence of a broader range of absorbing molecular species. Here we present a combined analysis of JWST transmission spectroscopy across four different instrumental modes spanning 0.5-5.2 micron using Early Release Science observations of the Saturn-mass exoplanet WASP-39b. Our uniform analysis constrains the orbital and stellar parameters within sub-percent precision, including matching the precision obtained by the most precise asteroseismology measurements of stellar density to-date, and further confirms the presence of Na, K, H$_2$O, CO, CO$_2$, and SO$_2$ atmospheric absorbers. Through this process, we also improve the agreement between the transmission spectra of all modes, except for the NIRSpec PRISM, which is affected by partial saturation of the detector. This work provides strong evidence that uniform light curve analysis is an important aspect to ensuring reliability when comparing the high-precision transmission spectra provided by JWST. △ Less

Submitted 18 July, 2024; originally announced July 2024.

Comments: 34 pages, 11 figures, 3 tables. Nat Astron (2024)

Nitrogen Abundance Distribution in the inner Milky Way

Authors: Jorge L. Pineda, Shinji Horiuchi, L. D. Anderson, Matteo Luisi, William D. Langer, Paul F. Goldsmith, Thomas B. H. Kuiper, Christian Fischer, Yan Gong, Andreas Brunthaler, Michael Rugel, Karl M. Menten

Abstract: We combine a new Galactic plane survey of Hydrogen Radio Recombination Lines (RRLs) with far-infrared (FIR) surveys of ionized Nitrogen, N+, to determine Nitrogen abundance across Galactic radius. RRLs were observed with NASA DSS-43 70m antenna and the Green Bank Telescope in 108 lines-of-sight spanning -135 degrees < l < 60 degrees, at b=0 degrees. These positions were also observed in [N II] 122… ▽ More We combine a new Galactic plane survey of Hydrogen Radio Recombination Lines (RRLs) with far-infrared (FIR) surveys of ionized Nitrogen, N+, to determine Nitrogen abundance across Galactic radius. RRLs were observed with NASA DSS-43 70m antenna and the Green Bank Telescope in 108 lines-of-sight spanning -135 degrees < l < 60 degrees, at b=0 degrees. These positions were also observed in [N II] 122 um and 205 um lines with the Herschel Space Observatory. Combining RRL and [N II] 122 um and 205 um observations in 41 of 108 samples with high signal-to-noise ratio, we studied ionized Nitrogen abundance distribution across Galactocentric distances of 0-8 kpc. Combined with existing Solar neighborhood and Outer galaxy N/H abundance determinations, we studied this quantity's distribution within the Milky Way's inner 17 kpc for the first time. We found a Nitrogen abundance gradient extending from Galactocentric radii of 4-17 kpc in the Galactic plane, while within 0-4 kpc, the N/H distribution remained flat. The gradient observed at large Galactocentric distances supports inside-out galaxy growth with the additional steepening resulting from variable star formation efficiency and/or radial flows in the Galactic disk, while the inner 4 kpc flattening, coinciding with the Galactic bar's onset, may be linked to radial flows induced by the bar potential. Using SOFIA/FIFI-LS and Herschel/PACS, we observed the [N III] 57 um line to trace doubly ionized gas contribution in a sub-sample of sightlines. We found negligible N++ contributions along these sightlines, suggesting mostly singly ionized Nitrogen originating from low ionization H II region outskirts. △ Less

Submitted 19 July, 2024; v1 submitted 17 July, 2024; originally announced July 2024.

Comments: Accepted for publication at the Astrophysical Journal. 25 pages, 13 figures

In-Situ Infrared Camera Monitoring for Defect and Anomaly Detection in Laser Powder Bed Fusion: Calibration, Data Mapping, and Feature Extraction

Authors: Shawn Hinnebusch, David Anderson, Berkay Bostan, Albert C. To

Abstract: Laser powder bed fusion (LPBF) process can incur defects due to melt pool instabilities, spattering, temperature increase, and powder spread anomalies. Identifying defects through in-situ monitoring typically requires collecting, storing, and analyzing large amounts of data generated. The first goal of this work is to propose a new approach to accurately map in-situ data to a three-dimensional (3D… ▽ More Laser powder bed fusion (LPBF) process can incur defects due to melt pool instabilities, spattering, temperature increase, and powder spread anomalies. Identifying defects through in-situ monitoring typically requires collecting, storing, and analyzing large amounts of data generated. The first goal of this work is to propose a new approach to accurately map in-situ data to a three-dimensional (3D) geometry, aiming to reduce the amount of storage. The second goal of this work is to introduce several new IR features for defect detection or process model calibration, which include laser scan order, local preheat temperature, maximum pre-laser scanning temperature, and number of spatters generated locally and their landing locations. For completeness, processing of other common IR features, such as interpass temperature, heat intensity, cooling rates, and melt pool area, are also presented with the underlying algorithm and Python implementation. A number of different parts are printed, monitored, and characterized to provide evidence of process defects and anomalies that different IR features are capable of detecting. △ Less

Submitted 17 July, 2024; originally announced July 2024.

Comments: 29 Pages, 19 Figures

The Most Sensitive Radio Recombination Line Measurements Ever Made of the Galactic Warm Ionized Medium

Authors: T. M. Bania, Dana S. Balser, Trey V. Wenger, Spencer J. Ireland, L. D. Anderson, Matteo Luisi

Abstract: Diffuse ionized gas pervades the disk of the Milky Way. We detect extremely faint emission from this Galactic Warm Ionized Medium (WIM) using the Green Bank Telescope to make radio recombination line (RRL) observations toward two Milky Way sight lines: G20, $(\ell,{\it b}) = (20^\circ, 0^\circ)$, and G45, $(\ell,{\it b}) = (45^\circ, 0^\circ)$. We stack 18 consecutive Hn$α$ transitions between 4.3… ▽ More Diffuse ionized gas pervades the disk of the Milky Way. We detect extremely faint emission from this Galactic Warm Ionized Medium (WIM) using the Green Bank Telescope to make radio recombination line (RRL) observations toward two Milky Way sight lines: G20, $(\ell,{\it b}) = (20^\circ, 0^\circ)$, and G45, $(\ell,{\it b}) = (45^\circ, 0^\circ)$. We stack 18 consecutive Hn$α$ transitions between 4.3-7.1 GHz to derive ${\rm \langle Hnα\rangle}$ spectra that are sensitive to RRL emission from plasmas with emission measures EM >10 ${\rm \,cm^{-6}\,pc}$. Each sight line has two Gaussian shaped spectral components with emission measures that range between $\sim$100 and $\sim$300 ${\rm \,cm^{-6}\,pc}$. Because there is no detectable RRL emission at negative LSR velocities the emitting plasma must be located interior to the Solar orbit. The G20 and G45 emission measures imply RMS densities of 0.15 and 0.18$\,{\rm cm^{-3}}$, respectively, if these sight lines are filled with homogeneous plasma. The observed ${\rm \langle Hnβ\rangle}$/${\rm \langle Hnα\rangle}$ line ratios are consistent with LTE excitation for the strongest components. The high velocity component of G20 has a narrow line width, 13.5 km s$^{-1}$, that sets an upper limit of <4,000 K for the plasma electron temperature. This is inconsistent with the ansatz of a canonically pervasive, low density, $\sim$ 10,000 K WIM plasma. △ Less

Submitted 8 July, 2024; originally announced July 2024.

Comments: 36 pages, accepted for publication in the Astrophysical Journal

Implications for Galactic Electron Density Structure from Pulsar Sightlines Intersecting HII Regions

Authors: S. K. Ocker, L. D. Anderson, T. J. W. Lazio, J. M. Cordes, V. Ravi

Abstract: Recent radio surveys have revealed pulsars with dispersion and scattering delays induced by ionized gas that are larger than the rest of the observed pulsar population, in some cases with electron column densities (or dispersion measures; DMs) larger than the maximum predictions of Galactic electron density models. By cross-matching the observed pulsar population against HII region catalogs, we sh… ▽ More Recent radio surveys have revealed pulsars with dispersion and scattering delays induced by ionized gas that are larger than the rest of the observed pulsar population, in some cases with electron column densities (or dispersion measures; DMs) larger than the maximum predictions of Galactic electron density models. By cross-matching the observed pulsar population against HII region catalogs, we show that the majority of pulsars with $\rm DM > 600$ pc cm$^{-3}$ and scattering delays $τ(1\ {\rm GHz}) > 10$ ms lie behind HII regions, and that HII region intersections may be relevant to as much as a third of the observed pulsar population. The fraction of the full pulsar population with sightlines intersecting HII regions is likely larger. Accounting for HII regions resolves apparent discrepancies where Galactic electron density models place high-DM pulsars beyond the Galactic disk. By comparing emission measures (EMs) inferred from recombination line observations to pulsar DMs, we show that HII regions can contribute tens to hundreds of pc cm$^{-3}$ in electron column density along a pulsar LOS. We find that nearly all pulsars with significant excess (and deficit) scattering from the mean $τ$-DM relation are spatially coincident with known discrete ionized gas structures, including HII regions. Accounting for HII regions is critical to the interpretation of radio dispersion and scattering measurements as electron density tracers, both in the Milky Way and in other galaxies. △ Less

Submitted 1 August, 2024; v1 submitted 11 June, 2024; originally announced June 2024.

Comments: 23 pages, 10 figures. Accepted to ApJ

Photo-dynamical characterisation of the TOI-178 resonant chain

Authors: A. Leleu, J. -B. Delisle, L. Delrez, E. M. Bryant, A. Brandeker, H. P. Osborn, N. Hara, T. G. Wilson, N. Billot, M. Lendl, D. Ehrenreich, H. Chakraborty, M. N. Günther, M. J. Hooton, Y. Alibert, R. Alonso, D. R. Alves, D. R. Anderson, I. Apergis, D. Armstrong, T. Bárczy, D. Barrado Navascues, S. C. C. Barros, M. P. Battley, W. Baumjohann , et al. (82 additional authors not shown)

Abstract: The TOI-178 system consists of a nearby late K-dwarf transited by six planets in the super-Earth to mini-Neptune regime, with radii ranging from 1.2 to 2.9 earth radius and orbital periods between 1.9 and 20.7 days. All planets but the innermost one form a chain of Laplace resonances. The fine-tuning and fragility of such orbital configurations ensure that no significant scattering or collision ev… ▽ More The TOI-178 system consists of a nearby late K-dwarf transited by six planets in the super-Earth to mini-Neptune regime, with radii ranging from 1.2 to 2.9 earth radius and orbital periods between 1.9 and 20.7 days. All planets but the innermost one form a chain of Laplace resonances. The fine-tuning and fragility of such orbital configurations ensure that no significant scattering or collision event has taken place since the formation and migration of the planets in the protoplanetary disc, hence providing important anchors for planet formation models. We aim to improve the characterisation of the architecture of this key system, and in particular the masses and radii of its planets. In addition, since this system is one of the few resonant chains that can be characterised by both photometry and radial velocities, we aim to use it as a test bench for the robustness of the planetary mass determination with each technique. We perform a global analysis of all available photometry and radial velocity. We also try different sets of priors on the masses and eccentricity, as well as different stellar activity models, to study their effects on the masses estimated by each method. We show how stellar activity is preventing us from obtaining a robust mass estimation for the three outer planets using radial velocity data alone. We also show that our joint photo-dynamical and radial velocity analysis resulted in a robust mass determination for planets c to g, with precision of 12% for the mass of planet c, and better than 10% for planets d to g. The new precisions on the radii range from 2 to 3%. The understanding of this synergy between photometric and radial velocity measurements will be valuable during the PLATO mission. We also show that TOI-178 is indeed currently locked in the resonant configuration, librating around an equilibrium of the chain. △ Less

Submitted 22 May, 2024; originally announced May 2024.

Journal ref: A&A 688, A211 (2024)

Deep Reinforcement Learning for Time-Critical Wilderness Search And Rescue Using Drones

Authors: Jan-Hendrik Ewers, David Anderson, Douglas Thomson

Abstract: Traditional search and rescue methods in wilderness areas can be time-consuming and have limited coverage. Drones offer a faster and more flexible solution, but optimizing their search paths is crucial. This paper explores the use of deep reinforcement learning to create efficient search missions for drones in wilderness environments. Our approach leverages a priori data about the search area and… ▽ More Traditional search and rescue methods in wilderness areas can be time-consuming and have limited coverage. Drones offer a faster and more flexible solution, but optimizing their search paths is crucial. This paper explores the use of deep reinforcement learning to create efficient search missions for drones in wilderness environments. Our approach leverages a priori data about the search area and the missing person in the form of a probability distribution map. This allows the deep reinforcement learning agent to learn optimal flight paths that maximize the probability of finding the missing person quickly. Experimental results show that our method achieves a significant improvement in search times compared to traditional coverage planning and search planning algorithms. In one comparison, deep reinforcement learning is found to outperform other algorithms by over $160\%$, a difference that can mean life or death in real-world search operations. Additionally, unlike previous work, our approach incorporates a continuous action space enabled by cubature, allowing for more nuanced flight patterns. △ Less

Submitted 22 May, 2024; v1 submitted 21 May, 2024; originally announced May 2024.

Comments: 16 pages, 19 figures. Submitted

A Novel Methodology for Autonomous Planetary Exploration Using Multi-Robot Teams

Authors: Sarah Swinton, Jan-Hendrik Ewers, Euan McGookin, David Anderson, Douglas Thomson

Abstract: One of the fundamental limiting factors in planetary exploration is the autonomous capabilities of planetary exploration rovers. This study proposes a novel methodology for trustworthy autonomous multi-robot teams which incorporates data from multiple sources (HiRISE orbiter imaging, probability distribution maps, and on-board rover sensors) to find efficient exploration routes in Jezero crater. A… ▽ More One of the fundamental limiting factors in planetary exploration is the autonomous capabilities of planetary exploration rovers. This study proposes a novel methodology for trustworthy autonomous multi-robot teams which incorporates data from multiple sources (HiRISE orbiter imaging, probability distribution maps, and on-board rover sensors) to find efficient exploration routes in Jezero crater. A map is generated, consisting of a 3D terrain model, traversability analysis, and probability distribution map of points of scientific interest. A three-stage mission planner generates an efficient route, which maximises the accumulated probability of identifying points of interest. A 4D RRT* algorithm is used to determine smooth, flat paths, and prioritised planning is used to coordinate a safe set of paths. The above methodology is shown to coordinate safe and efficient rover paths, which ensure the rovers remain within their nominal pitch and roll limits throughout operation. △ Less

Submitted 21 May, 2024; originally announced May 2024.

Comments: 6 pages. 10 figures. This work has been submitted to the IEEE for possible publication. Copyright may be transferred without notice, after which this version may no longer be accessible

Low-frequency absorption and radio recombination line features of the Galactic Center Lobe

Authors: Natasha Hurley-Walker, L. D. Anderson, M. Luisi, N. M. McClure-Griffiths, Robert A. Benjamin, Michael A. Kuhn, Dylan J. Linville, B. Liu, Catherine Zucker

Abstract: The Galactic center lobe (GCL) is a $\sim 1^\circ$ object located north of the Galactic center. In the mid-infrared (MIR), the GCL appears as two 8.0-micron filaments that roughly define an ellipse. There is strong 24-micron and radio continuum emission in the interior of the ellipse. Due to its morphology and location in the sky, previous authors have argued that the GCL is created by outflows fr… ▽ More The Galactic center lobe (GCL) is a $\sim 1^\circ$ object located north of the Galactic center. In the mid-infrared (MIR), the GCL appears as two 8.0-micron filaments that roughly define an ellipse. There is strong 24-micron and radio continuum emission in the interior of the ellipse. Due to its morphology and location in the sky, previous authors have argued that the GCL is created by outflows from star formation in the central molecular zone or by activity of the central black hole Sgr~A$^*$. We present images of the GCL from the GaLactic and Extragalactic All-sky Murchison Widefield Array survey in radio continuum that show thermal absorption against the Galactic center, incompatible with an interpretation of synchrotron self-absorption. Estimates of the cosmic ray emissivity in this direction allow us to place a distance constraint on the GCL. To be consistent with standard emissivity assumptions, the GCL would be located 2kpc away. At a distance of 8kpc, the synchrotron background emissivity is enhanced by $\sim75$% in the direction of the GCL. We also present radio recombination line data from the Green Bank Telescope that constrains the electron temperature and line widths in this region, which are also more explicable if the GCL lies relatively close. △ Less

Submitted 19 May, 2024; originally announced May 2024.

Comments: 6 figures, accepted to ApJ

The Galactic Center Lobe as an HII Region

Authors: L. D. Anderson, Matteo Luisi, B. Liu, Dylan J. Linville, Robert A. Benjamin, Natasha Hurley-Walker, N. M. McClure-Griffiths, Catherine Zucker

Abstract: The Galactic center lobe (GCL) is an object ~1° across that is located north of the Galactic center. In the mid-infrared (MIR) the GCL appears as two 8.0$μ$m filaments between which is strong 24$μ$m and radio continuum emission. Due to its morphology and location in the sky, previous authors have argued that the GCL is located in the Galactic center region, created by outflows from star formation… ▽ More The Galactic center lobe (GCL) is an object ~1° across that is located north of the Galactic center. In the mid-infrared (MIR) the GCL appears as two 8.0$μ$m filaments between which is strong 24$μ$m and radio continuum emission. Due to its morphology and location in the sky, previous authors have argued that the GCL is located in the Galactic center region, created by outflows from star formation or by activity of the central black hole Sagittarius A*. In an associated paper (Hurley-Walker et al., 2024, in press), low-frequency radio emission indicates that the GCL must instead lie foreground to the Galactic center. If the GCL is foreground to the Galactic center, it is likely to be a type of object common throughout the Galactic disk; we here investigate whether its properties are similar to those of Galactic HII regions. We find that the GCL's MIR morphology, MIR flux densities, dust temperatures, and radio recombination line (RRL) properties as traced by the GBT Diffuse Ionized Gas Survey (GDIGS) are consistent with those of known Galactic HII regions, although the derived electron temperature is low. We search for the ionizing source(s) of the possible HII region and identify a stellar cluster candidate (Camargo #1092/Ryu & Lee #532) and a cluster of young stellar objects (SPICY G359.3+0.3) whose members have Gaia parallaxes distances of 1.7${\pm}$0.4kpc. Taken together, the results of our companion paper and those shown here suggest that the GCL has properties consistent with those of an HII region located ~2kpc from the Sun. △ Less

Submitted 17 May, 2024; originally announced May 2024.

Comments: Accepted to ApJ

Enhancing Reinforcement Learning in Sensor Fusion: A Comparative Analysis of Cubature and Sampling-based Integration Methods for Rover Search Planning

Authors: Jan-Hendrik Ewers, Sarah Swinton, David Anderson, Euan McGookin, Douglas Thomson

Abstract: This study investigates the computational speed and accuracy of two numerical integration methods, cubature and sampling-based, for integrating an integrand over a 2D polygon. Using a group of rovers searching the Martian surface with a limited sensor footprint as a test bed, the relative error and computational time are compared as the area was subdivided to improve accuracy in the sampling-based… ▽ More This study investigates the computational speed and accuracy of two numerical integration methods, cubature and sampling-based, for integrating an integrand over a 2D polygon. Using a group of rovers searching the Martian surface with a limited sensor footprint as a test bed, the relative error and computational time are compared as the area was subdivided to improve accuracy in the sampling-based approach. The results show that the sampling-based approach exhibits a $14.75\%$ deviation in relative error compared to cubature when it matches the computational performance at $100\%$. Furthermore, achieving a relative error below $1\%$ necessitates a $10000\%$ increase in relative time to calculate due to the $\mathcal{O}(N^2)$ complexity of the sampling-based method. It is concluded that for enhancing reinforcement learning capabilities and other high iteration algorithms, the cubature method is preferred over the sampling-based method. △ Less

Submitted 15 August, 2024; v1 submitted 14 May, 2024; originally announced May 2024.

Comments: Submitted to IROS 2024

TOI-2447 b / NGTS-29 b: a 69-day Saturn around a Solar analogue

Authors: Samuel Gill, Daniel Bayliss, Solène Ulmer-Moll, Peter J. Wheatley, Rafael Brahm, David R. Anderson, David Armstrong, Ioannis Apergis, Douglas R. Alves, Matthew R. Burleigh, R. P. Butler, François Bouchy, Matthew P. Battley, Edward M. Bryant, Allyson Bieryla, Jeffrey D. Crane, Karen A. Collins, Sarah L. Casewell, Ilaria Carleo, Alastair B. Claringbold, Paul A. Dalba, Diana Dragomir, Philipp Eigmüller, Jan Eberhardt, Michael Fausnaugh , et al. (41 additional authors not shown)

Abstract: Discovering transiting exoplanets with relatively long orbital periods ($>$10 days) is crucial to facilitate the study of cool exoplanet atmospheres ($T_{\rm eq} < 700 K$) and to understand exoplanet formation and inward migration further out than typical transiting exoplanets. In order to discover these longer period transiting exoplanets, long-term photometric and radial velocity campaigns are r… ▽ More Discovering transiting exoplanets with relatively long orbital periods ($>$10 days) is crucial to facilitate the study of cool exoplanet atmospheres ($T_{\rm eq} < 700 K$) and to understand exoplanet formation and inward migration further out than typical transiting exoplanets. In order to discover these longer period transiting exoplanets, long-term photometric and radial velocity campaigns are required. We report the discovery of TOI-2447 b ($=$ NGTS-29b), a Saturn-mass transiting exoplanet orbiting a bright (T=10.0) Solar-type star (T$_{\rm eff}$=5730 K). TOI-2447 b was identified as a transiting exoplanet candidate from a single transit event of 1.3% depth and 7.29 h duration in $TESS$ Sector 31 and a prior transit event from 2017 in NGTS data. Four further transit events were observed with NGTS photometry which revealed an orbital period of P=69.34 days. The transit events establish a radius for TOI-2447 b of $0.865 \pm 0.010\rm R_{\rm J}$, while radial velocity measurements give a mass of $0.386 \pm 0.025 \rm M_{\rm J}$. The equilibrium temperature of the planet is $414$ K, making it much cooler than the majority of $TESS$ planet discoveries. We also detect a transit signal in NGTS data not caused by TOI-2447 b, along with transit timing variations and evidence for a $\sim$150 day signal in radial velocity measurements. It is likely that the system hosts additional planets, but further photometry and radial velocity campaigns will be needed to determine their parameters with confidence. TOI-2447 b/NGTS-29b joins a small but growing population of cool giants that will provide crucial insights into giant planet composition and formation mechanisms. △ Less

Submitted 12 May, 2024; originally announced May 2024.

Comments: 16 pages, 12 figures. Accepted for publication in MNRAS

Planet Hunters NGTS: New Planet Candidates from a Citizen Science Search of the Next Generation Transit Survey Public Data

Authors: Sean M. O'Brien, Megan E. Schwamb, Samuel Gill, Christopher A. Watson, Matthew R. Burleigh, Alicia Kendall, David R. Anderson, José I. Vines, James S. Jenkins, Douglas R. Alves, Laura Trouille, Solène Ulmer-Moll, Edward M. Bryant, Ioannis Apergis, Matthew P. Battley, Daniel Bayliss, Nora L. Eisner, Edward Gillen, Michael R. Goad, Maximilian N. Günther, Beth A. Henderson, Jeong-Eun Heo, David G. Jackson, Chris Lintott, James McCormac , et al. (13 additional authors not shown)

Abstract: We present the results from the first two years of the Planet Hunters NGTS citizen science project, which searches for transiting planet candidates in data from the Next Generation Transit Survey (NGTS) by enlisting the help of members of the general public. Over 8,000 registered volunteers reviewed 138,198 light curves from the NGTS Public Data Releases 1 and 2. We utilize a user weighting scheme… ▽ More We present the results from the first two years of the Planet Hunters NGTS citizen science project, which searches for transiting planet candidates in data from the Next Generation Transit Survey (NGTS) by enlisting the help of members of the general public. Over 8,000 registered volunteers reviewed 138,198 light curves from the NGTS Public Data Releases 1 and 2. We utilize a user weighting scheme to combine the classifications of multiple users to identify the most promising planet candidates not initially discovered by the NGTS team. We highlight the five most interesting planet candidates detected through this search, which are all candidate short-period giant planets. This includes the TIC-165227846 system that, if confirmed, would be the lowest-mass star to host a close-in giant planet. We assess the detection efficiency of the project by determining the number of confirmed planets from the NASA Exoplanet Archive and TESS Objects of Interest (TOIs) successfully recovered by this search and find that 74% of confirmed planets and 63% of TOIs detected by NGTS are recovered by the Planet Hunters NGTS project. The identification of new planet candidates shows that the citizen science approach can provide a complementary method to the detection of exoplanets with ground-based surveys such as NGTS. △ Less

Submitted 23 April, 2024; originally announced April 2024.

Comments: 42 pages, 20 figures, 17 tables. To be published in AJ

Journal ref: AJ 167 (2024) 238

Chemical mass-action systems as analog computers: implementing arithmetic computations at specified speed

Authors: David F. Anderson, Badal Joshi

Abstract: Recent technological advances allow us to view chemical mass-action systems as analog computers. In this context, the inputs to a computation are encoded as initial values of certain chemical species while the outputs are the limiting values of other chemical species. In this paper, we design chemical systems that carry out the elementary arithmetic computations of: identification, inversion, $m$t… ▽ More Recent technological advances allow us to view chemical mass-action systems as analog computers. In this context, the inputs to a computation are encoded as initial values of certain chemical species while the outputs are the limiting values of other chemical species. In this paper, we design chemical systems that carry out the elementary arithmetic computations of: identification, inversion, $m$th roots (for $m \ge 2$), addition, multiplication, absolute difference, rectified subtraction over non-negative real numbers, and partial real inversion over real numbers. We prove that these ``elementary modules'' have a speed of computation that is independent of the inputs to the computation. Moreover, we prove that finite sequences of such elementary modules, running in parallel, can carry out composite arithmetic over real numbers, also at a rate that is independent of inputs. Furthermore, we show that the speed of a composite computation is precisely the speed of the slowest elementary step. Specifically, the scale of the composite computation, i.e. the number of elementary steps involved in the composite, does not affect the overall asymptotic speed -- a feature of the parallel computing nature of our algorithm. Our proofs require the careful mathematical analysis of certain non-autonomous systems, and we believe this analysis will be useful in different areas of applied mathematics, dynamical systems, and the theory of computation. We close with a discussion on future research directions, including numerous important open theoretical questions pertaining to the field of computation with reaction networks. △ Less

Submitted 5 April, 2024; originally announced April 2024.

MSC Class: 37N25

NGTS-30 b/TOI-4862 b: An 1 Gyr old 98-day transiting warm Jupiter

Authors: M. P. Battley, K. A. Collins, S. Ulmer-Moll, S. N. Quinn, M. Lendl, S. Gill, R. Brahm, M. J. Hobson, H. P. Osborn, A. Deline, J. P. Faria, A. B. Claringbold, H. Chakraborty, K. G. Stassun, C. Hellier, D. R. Alves, C. Ziegler, D. R. Anderson, I. Apergis, D. J. Armstrong, D. Bayliss, Y. Beletsky, A. Bieryla, F. Bouchy, M. R. Burleigh , et al. (41 additional authors not shown)

Abstract: Long-period transiting exoplanets bridge the gap between the bulk of transit- and Doppler-based exoplanet discoveries, providing key insights into the formation and evolution of planetary systems. The wider separation between these planets and their host stars results in the exoplanets typically experiencing less radiation from their host stars; hence, they should maintain more of their original a… ▽ More Long-period transiting exoplanets bridge the gap between the bulk of transit- and Doppler-based exoplanet discoveries, providing key insights into the formation and evolution of planetary systems. The wider separation between these planets and their host stars results in the exoplanets typically experiencing less radiation from their host stars; hence, they should maintain more of their original atmospheres, which can be probed during transit via transmission spectroscopy. Although the known population of long-period transiting exoplanets is relatively sparse, surveys performed by the Transiting Exoplanet Survey Satellite (TESS) and the Next Generation Transit Survey (NGTS) are now discovering new exoplanets to fill in this crucial region of the exoplanetary parameter space. This study presents the detection and characterisation of NGTS-30 b/TOI-4862 b, a new long-period transiting exoplanet detected by following up on a single-transit candidate found in the TESS mission. Through monitoring using a combination of photometric instruments (TESS, NGTS, and EulerCam) and spectroscopic instruments (CORALIE, FEROS, HARPS, and PFS), NGTS-30 b/TOI-4862 b was found to be a long-period (P = 98.29838 day) Jupiter-sized (0.928 RJ; 0.960 MJ) planet transiting a 1.1 Gyr old G-type star. With a moderate eccentricity of 0.294, its equilibrium temperature could be expected to vary from 274 K to 500 K over the course of its orbit. Through interior modelling, NGTS-30 b/TOI-4862 b was found to have a heavy element mass fraction of 0.23 and a heavy element enrichment (Zp/Z_star) of 20, making it metal-enriched compared to its host star. NGTS-30 b/TOI-4862 b is one of the youngest well-characterised long-period exoplanets found to date and will therefore be important in the quest to understanding the formation and evolution of exoplanets across the full range of orbital separations and ages. △ Less

Submitted 3 April, 2024; originally announced April 2024.

Comments: 18 pages, 14 figures, accepted for publication in A&A

Audiosockets: A Python socket package for Real-Time Audio Processing

Authors: Nicolas Shu, David V. Anderson

Abstract: There are many packages in Python which allow one to perform real-time processing on audio data. Unfortunately, due to the synchronous nature of the language, there lacks a framework which allows for distributed parallel processing of the data without requiring a large programming overhead and in which the data acquisition is not blocked by subsequent processing operations. This work improves on p… ▽ More There are many packages in Python which allow one to perform real-time processing on audio data. Unfortunately, due to the synchronous nature of the language, there lacks a framework which allows for distributed parallel processing of the data without requiring a large programming overhead and in which the data acquisition is not blocked by subsequent processing operations. This work improves on packages used for audio data collection with a light-weight backend and a simple interface that allows for distributed processing through a socket-based structure. This is intended for real-time audio machine learning and data processing in Python with a quick deployment of multiple parallel operations on the same data, allowing users to spend less time debugging and more time developing. △ Less

Submitted 14 March, 2024; originally announced March 2024.

Comments: 4 pages, 2 figures

Molecular Gas Tracers in Young and Old Protoplanetary Disks

Authors: Dana E. Anderson, L. Ilsedore Cleeves, Geoffrey A. Blake, Chunhua Qi, Edwin A. Bergin, John M. Carpenter, Kamber R. Schwarz, Claire Thilenius, Ke Zhang

Abstract: Molecular emission is used to investigate both the physical and chemical properties of protoplanetary disks. Therefore, to accurately derive disk properties, we need a thorough understanding of the behavior of the molecular probes we rely on. Here we investigate how the molecular line emission of N$_2$H$^+$, HCO$^+$, HCN, and C$^{18}$O compare to other measured quantities in a set of 20 protoplane… ▽ More Molecular emission is used to investigate both the physical and chemical properties of protoplanetary disks. Therefore, to accurately derive disk properties, we need a thorough understanding of the behavior of the molecular probes we rely on. Here we investigate how the molecular line emission of N$_2$H$^+$, HCO$^+$, HCN, and C$^{18}$O compare to other measured quantities in a set of 20 protoplanetary disks. Overall, we find positive correlations between multiple line fluxes and the disk dust mass and radius. We also generally find strong positive correlations between the line fluxes of different molecular species. However, some disks do show noticeable differences in the relative fluxes of N$_2$H$^+$, HCO$^+$, HCN, and C$^{18}$O. These differences occur even within a single star-forming region. This results in a potentially large range of different disk masses and chemical compositions for systems of similar age and birth environment. While we make preliminary comparisons of molecular fluxes across different star-forming regions, more complete and uniform samples are needed in the future to search for trends with birth environment or age. △ Less

Submitted 7 March, 2024; originally announced March 2024.

Comments: Accepted for publication in The Astrophysical Journal

Mode Consensus Algorithms With Finite Convergence Time

Authors: Chao Huang, Hyungbo Shim, Siliang Yu, Brian D. O. Anderson

Abstract: This paper studies the distributed mode consensus problem in a multi-agent system, in which the agents each possess a certain attribute and they aim to agree upon the mode (the most frequent attribute owned by the agents) via distributed computation. Three algorithms are proposed. The first one directly calculates the frequency of all attributes at every agent, with protocols based on blended dyna… ▽ More This paper studies the distributed mode consensus problem in a multi-agent system, in which the agents each possess a certain attribute and they aim to agree upon the mode (the most frequent attribute owned by the agents) via distributed computation. Three algorithms are proposed. The first one directly calculates the frequency of all attributes at every agent, with protocols based on blended dynamics, and then returns the most frequent attribute as the mode. Assuming knowledge at each agent of a lower bound of the mode frequency as a priori information, the second algorithm is able to reduce the number of frequencies to be computed at every agent if the lower bound is large. The third algorithm further eliminates the need for this information by introducing an adaptive updating mechanism. The algorithms find the mode in finite time, and estimates of convergence time are provided. The proposed first and second algorithms enjoy the plug-and-play property with a dwell time. △ Less

Submitted 29 February, 2024; originally announced March 2024.

Mathematical models of drug delivery via a contact lens during wear

Authors: Daniel M. Anderson, Rayanne A. Luke

Abstract: In this work we develop and investigate mathematical and computational models that describe drug delivery from a contact lens during wear. Our models are designed to predict the dynamics of drug release from the contact lens and subsequent transport into the adjacent pre-lens tear film and post-lens tear film as well as into the ocular tissue (e.g. cornea), into the eyelid, and out of these region… ▽ More In this work we develop and investigate mathematical and computational models that describe drug delivery from a contact lens during wear. Our models are designed to predict the dynamics of drug release from the contact lens and subsequent transport into the adjacent pre-lens tear film and post-lens tear film as well as into the ocular tissue (e.g. cornea), into the eyelid, and out of these regions. These processes are modeled by one dimensional diffusion out of the lens coupled to compartment-type models for drug concentrations in the various accompanying regions. In addition to numerical solutions that are compared with experimental data on drug release in an in vitro eye model, we also identify a large diffusion limit model for which analytical solutions can be written down for all quantities of interest, such as cumulative release of the drug from the contact lens. We use our models to make assessments about possible mechanisms and drug transport pathways through the pre-lens and post-lens tear films and provide interpretation of experimental observations. We discuss successes and limitations of our models as well as their potential to guide further research to help understand the dynamics of ophthalmic drug delivery via drug-eluting contact lenses. △ Less

Submitted 7 February, 2024; originally announced March 2024.

Comments: 44 pages, 20 figures, 4 tables

Square-difference factor absorbing ideals of a commutative ring

Authors: David F. Anderson, Ayman Badawi, Jim Coykendall

Abstract: Let $R$ be a commutative ring with $1 \neq 0$. A proper ideal $I$ of $R$ is a {\it square-difference factor absorbing ideal} (sdf-absorbing ideal) of $R$ if whenever $a^2 - b^2 \in I$ for $0 \neq a, b \in R$, then $a + b \in I$ or $a - b \in I$. In this paper, we introduce and investigate sdf-absorbing ideals. Let $R$ be a commutative ring with $1 \neq 0$. A proper ideal $I$ of $R$ is a {\it square-difference factor absorbing ideal} (sdf-absorbing ideal) of $R$ if whenever $a^2 - b^2 \in I$ for $0 \neq a, b \in R$, then $a + b \in I$ or $a - b \in I$. In this paper, we introduce and investigate sdf-absorbing ideals. △ Less

Submitted 28 February, 2024; originally announced February 2024.

Comments: 18 pages

MSC Class: 13A15; 13F05; 13G05

NGTS-28Ab: A short period transiting brown dwarf

Authors: Beth A. Henderson, Sarah L. Casewell, Michael R. Goad, Jack S. Acton, Maximilian N. Günther, Louise D. Nielsen, Matthew R. Burleigh, Claudia Belardi, Rosanna H. Tilbrook, Oliver Turner, Steve B. Howell, Catherine A. Clark, Colin Littlefield, Khalid Barkaoui, Douglas R. Alves, David R. Anderson, Daniel Bayliss, Francois Bouchy, Edward M. Bryant, George Dransfield, Elsa Ducrot, Philipp Eigmüller, Samuel Gill, Edward Gillen, Michaël Gillon , et al. (21 additional authors not shown)

Abstract: We report the discovery of a brown dwarf orbiting a M1 host star. We first identified the brown dwarf within the Next Generation Transit Survey data, with supporting observations found in TESS sectors 11 and 38. We confirmed the discovery with follow-up photometry from the South African Astronomical Observatory, SPECULOOS-S, and TRAPPIST-S, and radial velocity measurements from HARPS, which allowe… ▽ More We report the discovery of a brown dwarf orbiting a M1 host star. We first identified the brown dwarf within the Next Generation Transit Survey data, with supporting observations found in TESS sectors 11 and 38. We confirmed the discovery with follow-up photometry from the South African Astronomical Observatory, SPECULOOS-S, and TRAPPIST-S, and radial velocity measurements from HARPS, which allowed us to characterise the system. We find an orbital period of ~1.25 d, a mass of 69.0+5.3-4.8 MJ, close to the Hydrogen burning limit, and a radius of 0.95 +- 0.05 RJ. We determine the age to be >0.5 Gyr, using model isochrones, which is found to be in agreement with SED fitting within errors. NGTS-28Ab is one of the shortest period systems found within the brown dwarf desert, as well as one of the highest mass brown dwarfs that transits an M dwarf. This makes NGTS-28Ab another important discovery within this scarcely populated region. △ Less

Submitted 15 February, 2024; originally announced February 2024.

Comments: 20 pages (inc. appendices), 16 figures, accepted for publication in MNRAS

Searching for Giant Exoplanets around M-dwarf Stars (GEMS) I: Survey Motivation

Authors: Shubham Kanodia, Caleb I. Cañas, Suvrath Mahadevan, Eric B. Ford, Ravit Helled, Dana E. Anderson, Alan Boss, William D. Cochran, Megan Delamer, Te Han, Jessica E. Libby-Roberts, Andrea S. J. Lin, Simon Müller, Paul Robertson, Guðmundur Stefánsson, Johanna Teske

Abstract: Recent discoveries of transiting giant exoplanets around M-dwarf stars (GEMS), aided by the all-sky coverage of TESS, are starting to stretch theories of planet formation through the core-accretion scenario. Recent upper limits on their occurrence suggest that they decrease with lower stellar masses, with fewer GEMS around lower-mass stars compared to solar-type. In this paper, we discuss existing… ▽ More Recent discoveries of transiting giant exoplanets around M-dwarf stars (GEMS), aided by the all-sky coverage of TESS, are starting to stretch theories of planet formation through the core-accretion scenario. Recent upper limits on their occurrence suggest that they decrease with lower stellar masses, with fewer GEMS around lower-mass stars compared to solar-type. In this paper, we discuss existing GEMS both through confirmed planets, as well as protoplanetary disk observations, and a combination of tests to reconcile these with theoretical predictions. We then introduce the \textit{Searching for GEMS} survey, where we utilize multi-dimensional nonparameteric statistics to simulate hypothetical survey scenarios to predict the required sample size of transiting GEMS with mass measurements to robustly compare their bulk-density with canonical hot-Jupiters orbiting FGK stars. Our Monte-Carlo simulations predict that a robust comparison requires about 40 transiting GEMS (compared to the existing sample of $\sim$ 15) with 5-$σ$ mass measurements. Furthermore, we discuss the limitations of existing occurrence estimates for GEMS, and provide a brief description of our planned systematic search to improve the occurrence rate estimates for GEMS. △ Less

Submitted 7 February, 2024; originally announced February 2024.

Comments: 16 pages + references, including 7 figures. Accepted in AAS Journals

An analytic approach for understanding mechanisms driving breakthrough infections

Authors: Amanda Brucker, Jillian H Hurst, Emily C O'Brien, Deverick Anderson, Michael E Yarrington, Jay Krishnan, Benjamin A Goldstein

Abstract: Real world data is an increasingly utilized resource for post-market monitoring of vaccines and provides insight into real world effectiveness. However, outside of the setting of a clinical trial, heterogeneous mechanisms may drive observed breakthrough infection rates among vaccinated individuals; for instance, waning vaccine-induced immunity as time passes and the emergence of a new strain again… ▽ More Real world data is an increasingly utilized resource for post-market monitoring of vaccines and provides insight into real world effectiveness. However, outside of the setting of a clinical trial, heterogeneous mechanisms may drive observed breakthrough infection rates among vaccinated individuals; for instance, waning vaccine-induced immunity as time passes and the emergence of a new strain against which the vaccine has reduced protection. Analyses of infection incidence rates are typically predicated on a presumed mechanism in their choice of an "analytic time zero" after which infection rates are modeled. In this work, we propose an explicit test for driving mechanism situated in a standard Cox proportional hazards framework. We explore the test's performance in simulation studies and in an illustrative application to real world data. We additionally introduce subgroup differences in infection incidence and evaluate the impact of time zero misspecification on bias and coverage of model estimates. In this study we observe strong power and controlled type I error of the test to detect the correct infection-driving mechanism under various settings. Similar to previous studies, we find mitigated bias and greater coverage of estimates when the analytic time zero is correctly specified or accounted for. △ Less

Submitted 30 January, 2024; originally announced January 2024.

Discovery of two warm mini-Neptunes with contrasting densities orbiting the young K3V star TOI-815

Authors: Angelica Psaridi, Hugh Osborn, François Bouchy, Monika Lendl, Léna Parc, Nicolas Billot, Christopher Broeg, Sérgio G. Sousa, Vardan Adibekyan, Omar Attia, Andrea Bonfanti, Hritam Chakraborty, Karen A. Collins, Jeanne Davoult, Elisa Delgado-Mena, Nolan Grieves, Tristan Guillot, Alexis Heitzmann, Ravit Helled, Coel Hellier, Jon M. Jenkins, Henrik Knierim, Andreas Krenn, JackJ. Lissauer, Rafael Luque , et al. (108 additional authors not shown)

Abstract: We present the discovery and characterization of two warm mini-Neptunes transiting the K3V star TOI-815 in a K-M binary system. Analysis of the spectra and rotation period reveal it to be a young star with an age of $200^{+400}_{-200}$Myr. TOI-815b has a 11.2-day period and a radius of 2.94$\pm$0.05$\it{R_{\rm\mathrm{\oplus}}}$ with transits observed by TESS, CHEOPS, ASTEP, and LCOGT. The outer pl… ▽ More We present the discovery and characterization of two warm mini-Neptunes transiting the K3V star TOI-815 in a K-M binary system. Analysis of the spectra and rotation period reveal it to be a young star with an age of $200^{+400}_{-200}$Myr. TOI-815b has a 11.2-day period and a radius of 2.94$\pm$0.05$\it{R_{\rm\mathrm{\oplus}}}$ with transits observed by TESS, CHEOPS, ASTEP, and LCOGT. The outer planet, TOI-815c, has a radius of 2.62$\pm$0.10$\it{R_{\rm\mathrm{\oplus}}}$, based on observations of three non-consecutive transits with TESS, while targeted CHEOPS photometry and radial velocity follow-up with ESPRESSO were required to confirm the 35-day period. ESPRESSO confirmed the planetary nature of both planets and measured masses of 7.6$\pm$1.5 $\it{M_{\rm \mathrm{\oplus}}}$ ($ρ_\mathrm{P}$=1.64$^{+0.33}_{-0.31}$gcm$^{-3}$) and 23.5$\pm$2.4$\it{M_{\rm\mathrm{\oplus}}}$ ($ρ_\mathrm{P}$=7.2$^{+1.1}_{-1.0}$gcm$^{-3}$) respectively. Thus, the planets have very different masses, unlike the usual similarity of masses in compact multi-planet systems. Moreover, our statistical analysis of mini-Neptunes orbiting FGK stars suggests that weakly irradiated planets tend to have higher bulk densities compared to those suffering strong irradiation. This could be ascribed to their cooler atmospheres, which are more compressed and denser. Internal structure modeling of TOI-815b suggests it likely has a H-He atmosphere constituting a few percent of the total planet mass, or higher if the planet is assumed to have no water. In contrast, the measured mass and radius of TOI-815c can be explained without invoking any atmosphere, challenging planetary formation theories. Finally, we infer from our measurements that the star is viewed close to pole-on, which implies a spin-orbit misalignment at the 3$σ$ level. △ Less

Submitted 30 January, 2024; v1 submitted 28 January, 2024; originally announced January 2024.

Comments: 24 pages, 27 figures, 6 tables

Principal Component Regression to Study the Impact of Economic Factors on Disadvantaged Communities

Authors: Narmadha M. Mohankumar, Milan Jain, Heng Wan, Sumitrra Ganguli, Kyle D. Wilson, David M. Anderson

Abstract: The Council on Environmental Quality's Climate and Economic Justice Screening Tool defines "disadvantaged communities" (DAC) in the USA, highlighting census tracts where benefits of climate and energy investments are not accruing. We use a principal component generalized linear model, which addresses the intertwined nature of economic factors, income and employment and model their relationship to… ▽ More The Council on Environmental Quality's Climate and Economic Justice Screening Tool defines "disadvantaged communities" (DAC) in the USA, highlighting census tracts where benefits of climate and energy investments are not accruing. We use a principal component generalized linear model, which addresses the intertwined nature of economic factors, income and employment and model their relationship to DAC status. Our study 1) identifies the most significant income groups and employment industries that impact DAC status, 2) provides the probability of DAC status across census tracts and compares the predictive accuracy with widely used machine learning approaches, 3) obtains historical predictions of the probability of DAC status, 4) obtains spatial downscaling of DAC status across block groups. Our study provides valuable insights for policymakers and stakeholders to develop strategies that promote sustainable development and address inequities in climate and energy investments in the USA. △ Less

Submitted 24 January, 2024; originally announced January 2024.

Comments: 13 pages, 9 figures, 2 tables

Covert Quantum Communication Over Optical Channels

Authors: Evan J. D. Anderson, Christopher K. Eyre, Isabel M. Dailey, Boulat A. Bash

Abstract: We explore the problem of covertly communicating qubits over the lossy thermal-noise bosonic channel, which is a quantum-mechanical model of many practical channels, including optical. Covert communication ensures that an adversary is unable to detect the presence of transmissions, which are concealed in channel noise. We investigate an achievable lower bound on quantum covert communication using… ▽ More We explore the problem of covertly communicating qubits over the lossy thermal-noise bosonic channel, which is a quantum-mechanical model of many practical channels, including optical. Covert communication ensures that an adversary is unable to detect the presence of transmissions, which are concealed in channel noise. We investigate an achievable lower bound on quantum covert communication using photonic dual-rail qubits. This encoding has practical significance, as it has been proposed for long-range repeater-based quantum communication over optical channels. △ Less

Submitted 23 January, 2024; v1 submitted 12 January, 2024; originally announced January 2024.

The SARAO MeerKAT 1.3 GHz Galactic Plane Survey

Authors: S. Goedhart, W. D. Cotton, F. Camilo, M. A. Thompson, G. Umana, M. Bietenholz, P. A. Woudt, L. D. Anderson, C. Bordiu, D. A. H. Buckley, C. S. Buemi, F. Bufano, F. Cavallaro, H. Chen, J. O. Chibueze, D. Egbo, B. S. Frank, M. G. Hoare, A. Ingallinera, T. Irabor, R. C. Kraan-Korteweg, S. Kurapati, P. Leto, S. Loru, M. Mutale , et al. (105 additional authors not shown)

Abstract: We present the SARAO MeerKAT Galactic Plane Survey (SMGPS), a 1.3 GHz continuum survey of almost half of the Galactic Plane (251°$\le l \le$ 358°and 2°$\le l \le$ 61°at $|b| \le 1.5°$). SMGPS is the largest, most sensitive and highest angular resolution 1 GHz survey of the Plane yet carried out, with an angular resolution of 8" and a broadband RMS sensitivity of $\sim$10--20 $μ$ Jy/beam. Here we d… ▽ More We present the SARAO MeerKAT Galactic Plane Survey (SMGPS), a 1.3 GHz continuum survey of almost half of the Galactic Plane (251°$\le l \le$ 358°and 2°$\le l \le$ 61°at $|b| \le 1.5°$). SMGPS is the largest, most sensitive and highest angular resolution 1 GHz survey of the Plane yet carried out, with an angular resolution of 8" and a broadband RMS sensitivity of $\sim$10--20 $μ$ Jy/beam. Here we describe the first publicly available data release from SMGPS which comprises data cubes of frequency-resolved images over 908--1656 MHz, power law fits to the images, and broadband zeroth moment integrated intensity images. A thorough assessment of the data quality and guidance for future usage of the data products are given. Finally, we discuss the tremendous potential of SMGPS by showcasing highlights of the Galactic and extragalactic science that it permits. These highlights include the discovery of a new population of non-thermal radio filaments; identification of new candidate supernova remnants, pulsar wind nebulae and planetary nebulae; improved radio/mid-IR classification of rare Luminous Blue Variables and discovery of associated extended radio nebulae; new radio stars identified by Bayesian cross-matching techniques; the realisation that many of the largest radio-quiet WISE HII region candidates are not true HII regions; and a large sample of previously undiscovered background HI galaxies in the Zone of Avoidance. △ Less

Submitted 2 May, 2024; v1 submitted 12 December, 2023; originally announced December 2023.

Comments: Accepted for publication in MNRAS. The data release is live and links can be found in the Data Availability Statement in the paper

Disentangling CO Chemistry in a Protoplanetary Disk Using Explanatory Machine Learning Techniques

Authors: Amina Diop, Ilse Cleeves, Dana Anderson, Jamila Pegues, Adele Plunkett

Abstract: Molecular abundances in protoplanetary disks are highly sensitive to the local physical conditions, including gas temperature, gas density, radiation field, and dust properties. Often multiple factors are intertwined, impacting the abundances of both simple and complex species. We present a new approach to understanding these chemical and physical interdependencies using machine learning. Specific… ▽ More Molecular abundances in protoplanetary disks are highly sensitive to the local physical conditions, including gas temperature, gas density, radiation field, and dust properties. Often multiple factors are intertwined, impacting the abundances of both simple and complex species. We present a new approach to understanding these chemical and physical interdependencies using machine learning. Specifically we explore the case of CO modeled under the conditions of a generic disk and build an explanatory regression model to study the dependence of CO spatial density on the gas density, gas temperature, cosmic ray ionization rate, X-ray ionization rate, and UV flux. Our findings indicate that combinations of parameters play a surprisingly powerful role in regulating CO compared to any singular physical parameter. Moreover, in general, we find the conditions in the disk are destructive toward CO. CO depletion is further enhanced in an increased cosmic ray environment and in disks with higher initial C/O ratios. These dependencies uncovered by our new approach are consistent with previous studies, which are more modeling intensive and computationally expensive. Our work thus shows that machine learning can be a powerful tool not only for creating efficient predictive models, but also for enabling a deeper understanding of complex chemical processes. △ Less

Submitted 8 December, 2023; originally announced December 2023.

Comments: Accepted in ApJ, 17 pages, 13 figures, 5 tables

Neighboring Extremal Optimal Control Theory for Parameter-Dependent Closed-loop Laws

Authors: Ayush Rai, Shaoshuai Mou, Brian D. O. Anderson

Abstract: This study introduces an approach to obtain a neighboring extremal optimal control (NEOC) solution for a closed-loop optimal control problem, applicable to a wide array of nonlinear systems and not necessarily quadratic performance indices. The approach involves investigating the variation incurred in the functional form of a known closed-loop optimal control law due to small, known parameter vari… ▽ More This study introduces an approach to obtain a neighboring extremal optimal control (NEOC) solution for a closed-loop optimal control problem, applicable to a wide array of nonlinear systems and not necessarily quadratic performance indices. The approach involves investigating the variation incurred in the functional form of a known closed-loop optimal control law due to small, known parameter variations in the system equations or the performance index. The NEOC solution can formally be obtained by solving a linear partial differential equation, akin to those encountered in the iterative solution of a nonlinear Hamilton-Jacobi equation. Motivated by numerical procedures for solving these latter equations, we also propose a numerical algorithm based on the Galerkin algorithm, leveraging the use of basis functions to solve the underlying Hamilton-Jacobi equation of the original optimal control problem. The proposed approach simplifies the NEOC problem by reducing it to the solution of a simple set of linear equations, thereby eliminating the need for a full re-solution of the adjusted optimal control problem. Furthermore, the variation to the optimal performance index can be obtained as a function of both the system state and small changes in parameters, allowing the determination of the adjustment to an optimal control law given a small adjustment of parameters in the system or the performance index. Moreover, in order to handle large known parameter perturbations, we propose a homotopic approach that breaks down the single calculation of NEOC into a finite set of multiple steps. Finally, the validity of the claims and theory is supported by theoretical analysis and numerical simulations. △ Less

Submitted 7 December, 2023; originally announced December 2023.

State of Health Estimation for Battery Modules with Parallel-Connected Cells Under Cell-to-Cell Variations

Authors: Qinan Zhou, Dyche Anderson, Jing Sun

Abstract: State of health (SOH) estimation for lithium-ion battery modules with cells connected in parallel is a challenging problem, especially with cell-to-cell variations. Incremental capacity analysis (ICA) and differential voltage analysis (DVA) are effective at the cell level, but a generalizable method to extend them to module-level SOH estimation remains missing, when only module-level measurements… ▽ More State of health (SOH) estimation for lithium-ion battery modules with cells connected in parallel is a challenging problem, especially with cell-to-cell variations. Incremental capacity analysis (ICA) and differential voltage analysis (DVA) are effective at the cell level, but a generalizable method to extend them to module-level SOH estimation remains missing, when only module-level measurements are available. This paper proposes a new method and demonstrates that, with multiple features systematically selected from the module-level ICA and DVA, the module-level SOH can be estimated with high accuracy and confidence in the presence of cell-to-cell variations. First, an information theory-based feature selection algorithm is proposed to find an optimal set of features for module-level SOH estimation. Second, a relevance vector regression (RVR)-based module-level SOH estimation model is proposed to provide both point estimates and three-sigma credible intervals while maintaining model sparsity. With more selected features incorporated, the proposed method achieves better estimation accuracy and higher confidence at the expense of higher model complexity. When applied to a large experimental dataset, the proposed method and the resulting sparse model lead to module-level SOH estimates with a 0.5% root-mean-square error and a 1.5% average three-sigma value. With all the training processes completed offboard, the proposed method has low computational complexity for onboard implementations. △ Less

Submitted 19 May, 2024; v1 submitted 5 December, 2023; originally announced December 2023.

Comments: Addressed reviewer comments: Combined two sections, revised dataset and module-level result sections, corrected a typo in Algorithm 2; Previous Edit Comments: Condensed abstract; Added details in Introduction, Dataset, Module-Level Result Sections; Revised Section I, III & VII, IX; Added the initialization of Phi in Algorithm 2

A resonant sextuplet of sub-Neptunes transiting the bright star HD 110067

Authors: R. Luque, H. P. Osborn, A. Leleu, E. Pallé, A. Bonfanti, O. Barragán, T. G. Wilson, C. Broeg, A. Collier Cameron, M. Lendl, P. F. L. Maxted, Y. Alibert, D. Gandolfi, J. -B. Delisle, M. J. Hooton, J. A. Egger, G. Nowak, M. Lafarga, D. Rapetti, J. D. Twicken, J. C. Morales, I. Carleo, J. Orell-Miquel, V. Adibekyan, R. Alonso , et al. (127 additional authors not shown)

Abstract: Planets with radii between that of the Earth and Neptune (hereafter referred to as sub-Neptunes) are found in close-in orbits around more than half of all Sun-like stars. Yet, their composition, formation, and evolution remain poorly understood. The study of multi-planetary systems offers an opportunity to investigate the outcomes of planet formation and evolution while controlling for initial con… ▽ More Planets with radii between that of the Earth and Neptune (hereafter referred to as sub-Neptunes) are found in close-in orbits around more than half of all Sun-like stars. Yet, their composition, formation, and evolution remain poorly understood. The study of multi-planetary systems offers an opportunity to investigate the outcomes of planet formation and evolution while controlling for initial conditions and environment. Those in resonance (with their orbital periods related by a ratio of small integers) are particularly valuable because they imply a system architecture practically unchanged since its birth. Here, we present the observations of six transiting planets around the bright nearby star HD 110067. We find that the planets follow a chain of resonant orbits. A dynamical study of the innermost planet triplet allowed the prediction and later confirmation of the orbits of the rest of the planets in the system. The six planets are found to be sub-Neptunes with radii ranging from 1.94 to 2.85 Re. Three of the planets have measured masses, yielding low bulk densities that suggest the presence of large hydrogen-dominated atmospheres. △ Less

Submitted 29 November, 2023; originally announced November 2023.

Comments: Published in Nature on November 30, 2023. Supplementary Information can be found in the online version of the paper in the journal

Journal ref: Nature 623, 932-937 (2023)

TESS Duotransit Candidates from the Southern Ecliptic Hemisphere

Authors: Faith Hawthorn, Sam Gill, Daniel Bayliss, Hugh P. Osborn, Ingrid Pelisoli, Toby Rodel, Kaylen Smith Darnbrook, Peter J. Wheatley, David R. Anderson, Ioan nis Apergis, Matthew P. Battley, Matthew R. Burleigh, Sarah L. Casewell, Philipp Eigmüller, Maximilian N. Günther, James S. Jenkins, Monika Lendl, Maximiliano Moyano, Ares Osborn, Gavin Ramsay, Solène Ulmer-Moll, Jose I. Vines, Richard West

Abstract: Discovering transiting exoplanets with long orbital periods allows us to study warm and cool planetary systems with temperatures similar to the planets in our own Solar system. The TESS mission has photometrically surveyed the entire Southern Ecliptic Hemisphere in Cycle 1 (August 2018 - July 2019), Cycle 3 (July 2020 - June 2021) and Cycle 5 (September 2022 - September 2023). We use the observati… ▽ More Discovering transiting exoplanets with long orbital periods allows us to study warm and cool planetary systems with temperatures similar to the planets in our own Solar system. The TESS mission has photometrically surveyed the entire Southern Ecliptic Hemisphere in Cycle 1 (August 2018 - July 2019), Cycle 3 (July 2020 - June 2021) and Cycle 5 (September 2022 - September 2023). We use the observations from Cycle 1 and Cycle 3 to search for exoplanet systems that show a single transit event in each year - which we call duotransits. The periods of these planet candidates are typically in excess of 20 days, with the lower limit determined by the duration of individual TESS observations. We find 85 duotransit candidates, which span a range of host star brightnesses between 8 < $T_{mag}$ < 14, transit depths between 0.1 per cent and 1.8 per cent, and transit durations between 2 and 10 hours with the upper limit determined by our normalisation function. Of these candidates, 25 are already known, and 60 are new. We present these candidates along with the status of photometric and spectroscopic follow-up. △ Less

Submitted 24 January, 2024; v1 submitted 26 October, 2023; originally announced October 2023.

Comments: 25 pages, 16 figures, submitted to Monthly Notices of the Royal Astronomical Society

A hot mini-Neptune and a temperate, highly eccentric sub-Saturn around the bright K-dwarf TOI-2134

Authors: F. Rescigno, G. Hébrard, A. Vanderburg, A. W. Mann, A. Mortier, S. Morrell, L. A. Buchhave, K. A. Collins, C. R. Mann, C. Hellier, R. D. Haywood, R. West, M. Stalport, N. Heidari, D. Anderson, C. X. Huang, M. López-Morales, P. Cortés-Zuleta, H. M. Lewis, X. Dumusque, I. Boisse, P. Rowden, A. Collier Cameron, M. Deleuil, M. Vezie , et al. (42 additional authors not shown)

Abstract: We present the characterisation of an inner mini-Neptune in a 9.2292005$\pm$0.0000063 day orbit and an outer mono-transiting sub-Saturn planet in a 95.50$^{+0.36}_{-0.25}$ day orbit around the moderately active, bright (mv=8.9 mag) K5V star TOI-2134. Based on our analysis of five sectors of TESS data, we determine the radii of TOI-2134b and c to be 2.69$\pm$0.16 R$_{e}$ for the inner planet and 7.… ▽ More We present the characterisation of an inner mini-Neptune in a 9.2292005$\pm$0.0000063 day orbit and an outer mono-transiting sub-Saturn planet in a 95.50$^{+0.36}_{-0.25}$ day orbit around the moderately active, bright (mv=8.9 mag) K5V star TOI-2134. Based on our analysis of five sectors of TESS data, we determine the radii of TOI-2134b and c to be 2.69$\pm$0.16 R$_{e}$ for the inner planet and 7.27$\pm$0.42 R$_{e}$ for the outer one. We acquired 111 radial-velocity spectra with HARPS-N and 108 radial-velocity spectra with SOPHIE. After careful periodogram analysis, we derive masses for both planets via Gaussian Process regression: 9.13$^{+0.78}_{-0.76}$ M$_{e}$ for TOI-2134b and 41.86$^{+7.69}_{-7.83}$ M$_{e}$ for TOI-2134c. We analysed the photometric and radial-velocity data first separately, then jointly. The inner planet is a mini-Neptune with density consistent with either a water-world or a rocky core planet with a low-mass H/He envelope. The outer planet has a bulk density similar to Saturn's. The outer planet is derived to have a significant eccentricity of 0.67$^{+0.05}_{-0.06}$ from a combination of photometry and RVs. We compute the irradiation of TOI-2134c as 1.45$\pm$0.10 times the bolometric flux received by Earth, positioning it for part of its orbit in the habitable sone of its system. We recommend further RV observations to fully constrain the orbit of TOI-2134c. With an expected Rossiter-McLaughlin (RM) effect amplitude of 7.2$\pm$1.3 m/s, we recommend TOI-2134c for follow-up RM analysis to study the spin-orbit architecture of the system. We calculate the Transmission Spectroscopy Metric, and both planets are suitable for bright-mode NIRCam atmospheric characterisation. △ Less

Submitted 20 October, 2023; originally announced October 2023.

Comments: 24 pages, 7 tables, 14 figures

High-angular-momentum Rydberg states in a room-temperature vapor cell for DC electric-field sensing

Authors: Alisher Duspayev, Ryan Cardman, David A. Anderson, Georg Raithel

Abstract: We prepare and analyze Rydberg states with orbital quantum numbers $\ell \le 6$ using three-optical-photon electromagnetically-induced transparency (EIT) and radio-frequency (RF) dressing, and employ the high-$\ell$ states in electric-field sensing. Rubidium-85 atoms in a room-temperature vapor cell are first promoted into the $25F_{5/2}$ state via Rydberg-EIT with three infrared laser beams. Two… ▽ More We prepare and analyze Rydberg states with orbital quantum numbers $\ell \le 6$ using three-optical-photon electromagnetically-induced transparency (EIT) and radio-frequency (RF) dressing, and employ the high-$\ell$ states in electric-field sensing. Rubidium-85 atoms in a room-temperature vapor cell are first promoted into the $25F_{5/2}$ state via Rydberg-EIT with three infrared laser beams. Two RF dressing fields then (near-)resonantly couple $25 \ell$ Rydberg states with high $\ell$. The dependence of the RF-dressed Rydberg-state level structure on RF powers, RF and laser frequencies is characterized using EIT. Furthermore, we discuss the principles of DC-electric-field sensing using high-$\ell$ Rydberg states, and experimentally demonstrate the method using test electric fields of $\lesssim$~50~V/m induced via photo-illumination of the vapor-cell wall. We measure the highly nonlinear dependence of the DC-electric-field strength on the power of the photo-illumination laser. Numerical calculations, which reproduce our experimental observations well, elucidate the underlying physics. Our study is relevant to high-precision spectroscopy of high-$\ell$ Rydberg states, Rydberg-atom-based electric-field sensing, and plasma electric-field diagnostics. △ Less

Submitted 16 October, 2023; originally announced October 2023.

Comments: 12 pages, 5 figures, 1 table

Journal ref: Phys. Rev. Research 6, 023138 (2024)

Methods for Averaging Spectral Line Data

Authors: L. D. Anderson, B. Liu, Dana S. Balser, T. M. Bania, L. M. Haffner, Dylan J. Linville, Matteo Luisi, Trey V. Wenger

Abstract: The ideal spectral averaging method depends on one's science goals and the available information about one's data. Including low-quality data in the average can decrease the signal-to-noise ratio (SNR), which may necessitate an optimization method or a consideration of different weighting schemes. Here, we explore a variety of spectral averaging methods. We investigate the use of three weighting s… ▽ More The ideal spectral averaging method depends on one's science goals and the available information about one's data. Including low-quality data in the average can decrease the signal-to-noise ratio (SNR), which may necessitate an optimization method or a consideration of different weighting schemes. Here, we explore a variety of spectral averaging methods. We investigate the use of three weighting schemes during averaging: weighting by the signal divided by the variance ("intensity-noise weighting"), weighting by the inverse of the variance ("noise weighting"), and uniform weighting. Whereas for intensity-noise weighting the SNR is maximized when all spectra are averaged, for noise and uniform weighting we find that averaging the 35-45% of spectra with the highest SNR results in the highest SNR average spectrum. With this intensity cutoff, the average spectrum with noise or uniform weighting has ~95% of the intensity of the spectrum created from intensity-noise weighting. We apply our spectral averaging methods to GBT Diffuse Ionized Gas (GDIGS) hydrogen radio recombination line (RRL) data to determine the ionic abundance ratio, y+, and discuss future applications of the methodology. △ Less

Submitted 17 October, 2023; v1 submitted 13 October, 2023; originally announced October 2023.

Comments: Accepted by PASP

Cold atomic gas identified by HI self-absorption. Cold atomic clouds toward giant molecular filaments

Authors: J. Syed, H. Beuther, P. F. Goldsmith, Th. Henning, M. Heyer, R. S. Klessen, J. M. Stil, J. D. Soler, L. D. Anderson, J. S. Urquhart, M. R. Rugel, K. G. Johnston, A. Brunthaler

Abstract: Stars form in the dense interiors of molecular clouds. The dynamics and physical properties of the atomic interstellar medium (ISM) set the conditions under which molecular clouds and eventually stars will form. It is, therefore, critical to investigate the relationship between the atomic and molecular gas phase to understand the global star formation process. Using the high angular resolution dat… ▽ More Stars form in the dense interiors of molecular clouds. The dynamics and physical properties of the atomic interstellar medium (ISM) set the conditions under which molecular clouds and eventually stars will form. It is, therefore, critical to investigate the relationship between the atomic and molecular gas phase to understand the global star formation process. Using the high angular resolution data from The HI/OH/Recombination line survey of the Milky Way (THOR), we aim to constrain the kinematic and physical properties of the cold atomic hydrogen gas phase toward the inner Galactic plane. HI self-absorption (HISA) has proven to be a viable method to detect cold atomic hydrogen clouds in the Galactic plane. With the help of a newly developed self-absorption extraction routine (astroSABER), we build upon previous case studies to identify HI self-absorption toward a sample of Giant Molecular Filaments (GMFs). We find the cold atomic gas to be spatially correlated with the molecular gas on a global scale. The column densities of the cold atomic gas traced by HISA are usually of the order of $10^{20}\rm\,cm^{-2}$ whereas those of molecular hydrogen traced by $\rm^{13}CO$ are at least an order of magnitude higher. The HISA column densities are attributed to a cold gas component that accounts for a fraction of $\sim$5% of the total atomic gas budget within the clouds. The HISA column density distributions show pronounced log-normal shapes that are broader than those traced by HI emission. The cold atomic gas is found to be moderately supersonic with Mach numbers of a $\sim$few. In contrast, highly supersonic dynamics drive the molecular gas within most filaments. △ Less

Submitted 3 October, 2023; originally announced October 2023.

Comments: 41 pages, 28 figures, accepted for publication in A&A

Journal ref: A&A 679, A130 (2023)

The GBT Diffuse Ionized Gas Survey (GDIGS): Discrete Sources

Authors: Dylan J. Linville, Matteo Luisi, L. D. Anderson, Bin Liu, T. M. Bania, Dana S. Balser, Trey V. Wenger, L. M. Haffner, J. L. Mascoop

Abstract: The Green Bank Telescope (GBT) Diffuse Ionized Gas Survey (GDIGS) traces ionized gas in the Galactic midplane by observing radio recombination line (RRL) emission from 4-8 GHz. The nominal survey zone is $32.3^{\circ}> {\ell} > -5^{\circ}$, $|b|<0.5^{\circ}$. Here, we analyze GDIGS Hn$α$ ionized gas emission toward discrete sources. Using GDIGS data, we identify the velocity of 35 H II regions tha… ▽ More The Green Bank Telescope (GBT) Diffuse Ionized Gas Survey (GDIGS) traces ionized gas in the Galactic midplane by observing radio recombination line (RRL) emission from 4-8 GHz. The nominal survey zone is $32.3^{\circ}> {\ell} > -5^{\circ}$, $|b|<0.5^{\circ}$. Here, we analyze GDIGS Hn$α$ ionized gas emission toward discrete sources. Using GDIGS data, we identify the velocity of 35 H II regions that have multiple detected RRL velocity components. We identify and characterize RRL emission from 88 H II regions that previously lacked measured ionized gas velocities. We also identify and characterize RRL emission from eight locations that appear to be previously-unidentified H II regions and 30 locations of RRL emission that do not appear to be H II regions based on their lack of mid-infrared emission. This latter group may be a compact component of the Galactic Diffuse Ionized Gas (DIG). There are an additional 10 discrete sources that have anomalously high RRL velocities for their locations in the Galactic plane. We compare these objects' RRL data to 13CO, H I and mid-infrared data, and find that these sources do not have the expected 24 $μ$m emission characteristic of H II regions. Based on this comparison we do not think these objects are H II regions, but we are unable to classify them as a known type of object. △ Less

Submitted 2 October, 2023; originally announced October 2023.

Comments: Accepted for publication by ApJ

Competitive Networked Bivirus SIS spread over Hypergraphs

Authors: Sebin Gracy, Brian D. O. Anderson, Mengbin Ye, Cesar A. Uribe

Abstract: The paper deals with the spread of two competing viruses over a network of population nodes, accounting for pairwise interactions and higher-order interactions (HOI) within and between the population nodes. We study the competitive networked bivirus susceptible-infected-susceptible (SIS) model on a hypergraph introduced in Cui et al. [1]. We show that the system has, in a generic sense, a finite n… ▽ More The paper deals with the spread of two competing viruses over a network of population nodes, accounting for pairwise interactions and higher-order interactions (HOI) within and between the population nodes. We study the competitive networked bivirus susceptible-infected-susceptible (SIS) model on a hypergraph introduced in Cui et al. [1]. We show that the system has, in a generic sense, a finite number of equilibria, and the Jacobian associated with each equilibrium point is nonsingular; the key tool is the Parametric Transversality Theorem of differential topology. Since the system is also monotone, it turns out that the typical behavior of the system is convergence to some equilibrium point. Thereafter, we exhibit a tri-stable domain with three locally exponentially stable equilibria. For different parameter regimes, we establish conditions for the existence of a coexistence equilibrium (both viruses infect separate fractions of each population node). △ Less

Submitted 25 September, 2023; originally announced September 2023.

Longitudinal and transverse spin transfer to $Λ$ and $\overlineΛ$ hyperons in polarized $p$+$p$ collisions at $\sqrt{s} = 200$ GeV

Authors: STAR Collaboration, M. I. Abdulhamid, B. E. Aboona, J. Adam, L. Adamczyk, J. R. Adams, I. Aggarwal, M. M. Aggarwal, Z. Ahammed, D. M. Anderson, E. C. Aschenauer, S. Aslam, J. Atchison, V. Bairathi, W. Baker, J. G. Ball Cap, K. Barish, R. Bellwied, P. Bhagat, A. Bhasin, S. Bhatta, J. Bielcik, J. Bielcikova, J. D. Brandenburg, X. Z. Cai , et al. (357 additional authors not shown)

Abstract: The longitudinal and transverse spin transfers to $Λ$ ($\overlineΛ$) hyperons in polarized proton-proton collisions are expected to be sensitive to the helicity and transversity distributions, respectively, of (anti-)strange quarks in the proton, and to the corresponding polarized fragmentation functions. We report improved measurements of the longitudinal spin transfer coefficient, $D_{LL}$, and… ▽ More The longitudinal and transverse spin transfers to $Λ$ ($\overlineΛ$) hyperons in polarized proton-proton collisions are expected to be sensitive to the helicity and transversity distributions, respectively, of (anti-)strange quarks in the proton, and to the corresponding polarized fragmentation functions. We report improved measurements of the longitudinal spin transfer coefficient, $D_{LL}$, and the transverse spin transfer coefficient, $D_{TT}$, to $Λ$ and $\overlineΛ$ in polarized proton-proton collisions at $\sqrt{s}$ = 200 GeV by the STAR experiment at RHIC. The data set includes longitudinally polarized proton-proton collisions with an integrated luminosity of 52 pb$^{-1}$, and transversely polarized proton-proton collisions with a similar integrated luminosity. Both data sets have about twice the statistics of previous results and cover a kinematic range of $|η_{Λ(\overlineΛ)}|$ $<$ 1.2 and transverse momentum $p_{T,{Λ(\overlineΛ)}}$ up to 8 GeV/$c$. We also report the first measurements of the hyperon spin transfer coefficients $D_{LL}$ and $D_{TT}$ as a function of the fractional jet momentum $z$ carried by the hyperon, which can provide more direct constraints on the polarized fragmentation functions. △ Less

Submitted 7 December, 2023; v1 submitted 25 September, 2023; originally announced September 2023.

The development of HISPEC for Keck and MODHIS for TMT: science cases and predicted sensitivities

Authors: Quinn M. Konopacky, Ashley D. Baker, Dimitri Mawet, Michael P. Fitzgerald, Nemanja Jovanovic, Charles Beichman, Garreth Ruane, Rob Bertz, Hiroshi Terada, Richard Dekany, Larry Lingvay, Marc Kassis, David Anderson, Motohide Tamura, Bjorn Benneke, Thomas Beatty, Tuan Do, Shogo Nishiyama, Peter Plavchan, Jason Wang, Ji Wang, Adam Burgasser, Jean-Baptiste Ruffio, Huihao Zhang, Aaron Brown , et al. (50 additional authors not shown)

Abstract: HISPEC is a new, high-resolution near-infrared spectrograph being designed for the W.M. Keck II telescope. By offering single-shot, R=100,000 between 0.98 - 2.5 um, HISPEC will enable spectroscopy of transiting and non-transiting exoplanets in close orbits, direct high-contrast detection and spectroscopy of spatially separated substellar companions, and exoplanet dynamical mass and orbit measureme… ▽ More HISPEC is a new, high-resolution near-infrared spectrograph being designed for the W.M. Keck II telescope. By offering single-shot, R=100,000 between 0.98 - 2.5 um, HISPEC will enable spectroscopy of transiting and non-transiting exoplanets in close orbits, direct high-contrast detection and spectroscopy of spatially separated substellar companions, and exoplanet dynamical mass and orbit measurements using precision radial velocity monitoring calibrated with a suite of state-of-the-art absolute and relative wavelength references. MODHIS is the counterpart to HISPEC for the Thirty Meter Telescope and is being developed in parallel with similar scientific goals. In this proceeding, we provide a brief overview of the current design of both instruments, and the requirements for the two spectrographs as guided by the scientific goals for each. We then outline the current science case for HISPEC and MODHIS, with focuses on the science enabled for exoplanet discovery and characterization. We also provide updated sensitivity curves for both instruments, in terms of both signal-to-noise ratio and predicted radial velocity precision. △ Less

Submitted 19 September, 2023; originally announced September 2023.

Comments: 25 pages, 9 figures. To appear in the Proceedings of SPIE: Techniques and Instrumentation for Detection of Exoplanets XI, vol. 12680 (2023)

Statistical Signatures of Nanoflare Activity. III. Evidence of Enhanced Nanoflaring Rates in Fully Convective stars as Observed by the NGTS

Authors: S. D. T. Grant, D. B. Jess, C. J. Dillon, M. Mathioudakis, C. A. Watson, J. A. G. Jackman, D. G. Jackson, P. J. Wheatley, M. R. Goad, S. L. Casewell, D. R. Anderson, M. R. Burleigh, R. G. West, J. I. Vines

Abstract: Previous examinations of fully-convective M-dwarf stars have highlighted enhanced rates of nanoflare activity on these distant stellar sources. However, the specific role the convective boundary, which is believed to be present for spectral types earlier than M2.5V, plays on the observed nanoflare rates is not yet known. Here, we utilize a combination of statistical and Fourier techniques to exami… ▽ More Previous examinations of fully-convective M-dwarf stars have highlighted enhanced rates of nanoflare activity on these distant stellar sources. However, the specific role the convective boundary, which is believed to be present for spectral types earlier than M2.5V, plays on the observed nanoflare rates is not yet known. Here, we utilize a combination of statistical and Fourier techniques to examine M-dwarf stellar lightcurves that lie on either side of the convective boundary. We find that fully convective M2.5V (and later sub-types) stars have greatly enhanced nanoflare rates compared with their pre-dynamo mode transition counterparts. Specifically, we derive a flaring power-law index in the region of $3.00 \pm 0.20$, alongside a decay timescale of $200 \pm 100$~s for M2.5V and M3V stars, matching those seen in prior observations of similar stellar sub-types. Interestingly, M4V stars exhibit longer decay timescales of $450 \pm 50$~s, along with an increased power-law index of $3.10 \pm 0.18$, suggesting an interplay between the rate of nanoflare occurrence and the intrinsic plasma parameters, for example, the underlying Lundquist number. In contrast, partially convective (i.e., earlier sub-types from M0V to M2V) M-dwarf stars exhibit very weak nanoflare activity, which is not easily identifiable using statistical or Fourier techniques. This suggests that fully convective stellar atmospheres favor small-scale magnetic reconnection, leading to implications for the flare-energy budgets of these stars. Understanding why small-scale reconnection is enhanced in fully convective atmospheres may help solve questions relating to the dynamo behavior of these stellar sources. △ Less

Submitted 18 September, 2023; originally announced September 2023.

Comments: 16 pages, 5 figures, 5 tables. Accepted to ApJ

A new path method for exponential ergodicity of Markov processes on $\mathbb Z^d$, with applications to stochastic reaction networks

Authors: David F. Anderson, Daniele Cappelletti, Wai-Tong Louis Fan, Jinsu Kim

Abstract: This paper provides a new path method that can be used to determine when an ergodic continuous-time Markov chain on $\mathbb Z^d$ converges exponentially fast to its stationary distribution in $L^2$. Specifically, we provide general conditions that guarantee the positivity of the spectral gap. Importantly, our results do not require the assumption of time-reversibility of the Markov model. We then… ▽ More This paper provides a new path method that can be used to determine when an ergodic continuous-time Markov chain on $\mathbb Z^d$ converges exponentially fast to its stationary distribution in $L^2$. Specifically, we provide general conditions that guarantee the positivity of the spectral gap. Importantly, our results do not require the assumption of time-reversibility of the Markov model. We then apply our new method to the well-studied class of stochastically modeled reaction networks. Notably, we show that each complex-balanced model that is also ``open'' has a positive spectral gap, and is therefore exponentially ergodic. We further illustrate how our results can be applied for models that are not necessarily complex-balanced. Moreover, we provide an example of a detailed-balanced (in the sense of reaction network theory), and hence complex-balanced, stochastic reaction network that is not exponentially ergodic. We believe this to be the first such example in the literature. △ Less

Submitted 29 September, 2023; v1 submitted 13 September, 2023; originally announced September 2023.

Comments: 44 pages

MSC Class: 60J27; 60J28

Protoplanetary Disk Chemistry

Authors: Karin I. Öberg, Stefano Facchini, Dana E. Anderson

Abstract: Planets form in disks of gas and dust around young stars. The disk molecular reservoirs and their chemical evolution affect all aspects of planet formation, from the coagulation of dust grains into pebbles, to the elemental and molecular compositions of the mature planet. Disk chemistry also enables unique probes of disk structures and dynamics, including those directly linked to ongoing planet fo… ▽ More Planets form in disks of gas and dust around young stars. The disk molecular reservoirs and their chemical evolution affect all aspects of planet formation, from the coagulation of dust grains into pebbles, to the elemental and molecular compositions of the mature planet. Disk chemistry also enables unique probes of disk structures and dynamics, including those directly linked to ongoing planet formation. Here we review the protoplanetary disk chemistry of the volatile elements HOCNSP, the associated observational and theoretical methods, and the links between disk and planet chemical compositions. Three takeaways from this review are: (1) The disk chemical composition, including the organic reservoirs, is set by both inheritance and in situ chemistry. (2) Disk gas and solid O/C/N/H elemental ratios often deviate from stellar values due to a combination of condensation of molecular carriers, chemistry, and dynamics. (3) Chemical, physical, and dynamical processes in disks are closely linked, which complicates disk chemistry modeling, but these links also present an opportunity to develop chemical probes of different aspects of disk evolution and planet formation. △ Less

Submitted 11 September, 2023; originally announced September 2023.

Comments: Invited review published in Annual Reviews of Astronomy and Astrophysics; 44 pages, 12 figures

Journal ref: Annual Review of Astronomy and Astrophysics 2023, 61:287-328

Two mini-Neptunes Transiting the Adolescent K-star HIP 113103 Confirmed with TESS and CHEOPS

Authors: Nataliea Lowson, George Zhou, Chelsea X. Huang, Duncan J. Wright, Billy Edwards, Emma Nabbie, Alex Venner, Samuel N. Quinn, Karen A. Collins, Edward Gillen, Matthew Battley, Amaury Triaud, Coel Hellier, Sara Seager, Joshua N. Winn, Jon M. Jenkins, Bill Wohler, Avi Shporer, Richard P. Schwarz, Felipe Murgas, Enric Pallé, David R. Anderson, Richard G. West, Robert A. Wittenmyer, Brendan P. Bowler , et al. (9 additional authors not shown)

Abstract: We report the discovery of two mini-Neptunes in near 2:1 resonance orbits ($P=7.610303$ d for HIP 113103 b and $P=14.245651$ d for HIP 113103 c) around the adolescent K-star HIP 113103 (TIC 121490076). The planet system was first identified from the TESS mission, and was confirmed via additional photometric and spectroscopic observations, including a $\sim$17.5 hour observation for the transits of… ▽ More We report the discovery of two mini-Neptunes in near 2:1 resonance orbits ($P=7.610303$ d for HIP 113103 b and $P=14.245651$ d for HIP 113103 c) around the adolescent K-star HIP 113103 (TIC 121490076). The planet system was first identified from the TESS mission, and was confirmed via additional photometric and spectroscopic observations, including a $\sim$17.5 hour observation for the transits of both planets using ESA CHEOPS. We place $\leq4.5$ min and $\leq2.5$ min limits on the absence of transit timing variations over the three year photometric baseline, allowing further constraints on the orbital eccentricities of the system beyond that available from the photometric transit duration alone. With a planetary radius of $R_{p}=1.829^{+0.096}_{-0.067}\,R_{\oplus}$, HIP 113103 b resides within the radius gap, and this might provide invaluable information on the formation disparities between super-Earths and mini-Neptunes. Given the larger radius $R_{p}=2.40^{+0.10}_{-0.08}\,R_{\oplus}$ for HIP 113103 c, and close proximity of both planets to HIP 113103, it is likely that HIP 113103 b might have lost (or is still losing) its primordial atmosphere. We therefore present simulated atmospheric transmission spectra of both planets using JWST, HST, and Twinkle. It demonstrates a potential metallicity difference (due to differences in their evolution) would be a challenge to detect if the atmospheres are in chemical equilibrium. As one of the brightest multi sub-Neptune planet systems suitable for atmosphere follow up, HIP 113103 b and HIP 113103 c could provide insight on planetary evolution for the sub-Neptune K-star population. △ Less

Submitted 28 January, 2024; v1 submitted 8 September, 2023; originally announced September 2023.

Comments: 18 pages, 12 figures, published in the Monthly Notices of the Royal Astronomical Society