-
Large dark matter content and steep metallicity profile predicted for Ultra-Diffuse Galaxies formed in high-spin halos
Authors:
José A. Benavides,
Laura V. Sales,
Mario. G. Abadi,
Mark Vogelsberger,
Federico Marinacci,
Lars Hernquist
Abstract:
We study the stellar properties of a sample of simulated ultra-diffuse galaxies (UDGs) with stellar mass $M_\star=10^{7.5}$ - $10^{9} ~ \rm{M_{\odot}}$, selected from the TNG50 simulation, where UDGs form mainly in high-spin dwarf-mass halos. We divide our sample into star-forming and quenched UDGs, finding good agreement with the stellar assembly history measured in observations. Star-forming UDG…
▽ More
We study the stellar properties of a sample of simulated ultra-diffuse galaxies (UDGs) with stellar mass $M_\star=10^{7.5}$ - $10^{9} ~ \rm{M_{\odot}}$, selected from the TNG50 simulation, where UDGs form mainly in high-spin dwarf-mass halos. We divide our sample into star-forming and quenched UDGs, finding good agreement with the stellar assembly history measured in observations. Star-forming UDGs and quenched UDGs with $M_\star \geq 10^8\; \rm M_\odot$ in our sample are particularly inefficient at forming stars, having $2$ - $10$ times less stellar mass than non-UDGs for the same virial mass halo. These results are consistent with recent mass inferences in UDG samples and suggest that the most inefficient UDGs arise from a late assembly of the dark matter mass followed by a stellar growth that is comparatively slower (for star-forming UDGs) or that was interrupted due to environmental removal of the gas (for quenched UDGs). Regardless of efficiency, UDGs are $60\%$ poorer in [Fe/H] than the population of non-UDGs at a fixed stellar mass, with the most extreme objects having metal content consistent with the simulated mass-metallicity relation at $z \sim 2$. Quenched UDGs stop their star formation in shorter timescales than non-UDGs of similar mass and are, as a consequence, alpha-enhanced with respect to non-UDGs. We identify metallicity profiles in UDGs as a potential avenue to distinguish between different formation paths for these galaxies, where gentle formation as a result of high-spin halos would present well-defined declining metallicity radial profiles while powerful-outflows or tidal stripping formation models would lead to flatter or constant metallicity as a function of radius due
△ Less
Submitted 22 July, 2024;
originally announced July 2024.
-
RMS asymmetry: a robust metric of galaxy shapes in images with varied depth and resolution
Authors:
Elizaveta Sazonova,
Cameron Morgan,
Michael Balogh,
Katherine Alatalo,
Jose A. Benavides,
Asa Bluck,
Sarah Brough,
Innocenza Busa,
Ricardo Demarco,
Darko Donevski,
Miguel Figueira,
Garreth Martin,
Vicente Rodriguez-Gomez,
Javier Román,
Kate Rowlands
Abstract:
Structural disturbances, such as galaxy mergers or instabilities, are key candidates for driving galaxy evolution, so it is important to detect and quantify galaxies hosting these disturbances spanning a range of masses, environments, and cosmic times. Traditionally, this is done by quantifying the asymmetry of a galaxy as part of the concentration-asymmetry-smoothness system, $A_{\rm{CAS}}$, and…
▽ More
Structural disturbances, such as galaxy mergers or instabilities, are key candidates for driving galaxy evolution, so it is important to detect and quantify galaxies hosting these disturbances spanning a range of masses, environments, and cosmic times. Traditionally, this is done by quantifying the asymmetry of a galaxy as part of the concentration-asymmetry-smoothness system, $A_{\rm{CAS}}$, and selecting galaxies above a certain threshold as merger candidates. However, in this work, we show that $A_{\rm{CAS}}$, is extremely dependent on imaging properties -- both resolution and depth -- and thus defining a single $A_{\rm{CAS}}$ threshold is impossible. We analyze an alternative root-mean-squared asymmetry, $A_{\rm{RMS}}$, and show that it is independent of noise down to the average SNR per pixel of 1. However, both metrics depend on the resolution. We argue that asymmetry is, by design, always a scale-dependent measurement, and it is essential to define an asymmetry at a given physical resolution, where the limit should be defined by the size of the smallest features one wishes to detect. We measure asymmetry of a set of $z\approx0.1$ galaxies observed with HST, HSC, and SDSS, and show that after matching the resolution of all images to 200 pc, we are able to obtain consistent $A_{\rm{RMS, 200pc}}$ measurements with all three instruments despite the vast differences in the original resolution or depth. We recommend that future studies use $A_{\rm{RMS, x pc}}$ measurement when evaluating asymmetry, where $x$ is defined by the physical size of the features of interest, and is kept consistent across the dataset, especially when the redshift or image properties of galaxies in the dataset vary.
△ Less
Submitted 8 April, 2024;
originally announced April 2024.
-
Size-Mass Relations for Simulated Low-Mass Galaxies: Mock Imaging versus Intrinsic Properties
Authors:
Courtney Klein,
James S. Bullock,
Jorge Moreno,
Francisco J. Mercado,
Philip F. Hopkins,
Rachel K. Cochrane,
Jose A. Benavides
Abstract:
The observationally-inferred size versus stellar-mass relationship (SMR) for low-mass galaxies provides an important test for galaxy formation models. However, the relationship relies on assumptions that relate observed luminosity profiles to underlying stellar mass profiles. Here we use the Feedback in Realistic Environments simulations of low-mass galaxies to explore how the predicted SMR change…
▽ More
The observationally-inferred size versus stellar-mass relationship (SMR) for low-mass galaxies provides an important test for galaxy formation models. However, the relationship relies on assumptions that relate observed luminosity profiles to underlying stellar mass profiles. Here we use the Feedback in Realistic Environments simulations of low-mass galaxies to explore how the predicted SMR changes depending on whether one uses star-particle counts directly or mock observations. We reproduce the SMR found in The Exploration of Local Volume Satellites survey remarkably well only when we infer stellar masses and sizes using mock observations. However, when we use star particles to directly infer stellar masses and half-mass radii, we find that our galaxies are too large and obey a SMR with too little scatter compared to observations. This discrepancy between the "true" galaxy size and mass and those derived in the mock observation approach is twofold. First, our simulated galaxies have higher and more varied MLRs at a fixed colour than those commonly-adopted, which tends to underestimate their stellar masses compared to their true, simulated values. Second, our galaxies have radially increasing MLR gradients therefore using a single MLR tends to under-predict the mass in the outer regions. Similarly, the true half-mass radius is larger than the half-light radius because the light is more concentrated than the mass. If our simulations are accurate representations of the real universe, then the relationship between galaxy size and stellar mass is even tighter for low-mass galaxies than is commonly inferred from observed relations.
△ Less
Submitted 13 June, 2024; v1 submitted 2 April, 2024;
originally announced April 2024.
-
Attenuation proxy hidden in surface brightness-colour diagrams. A new strategy for the LSST era
Authors:
K. Małek,
Junais,
A. Pollo,
M. Boquien,
V. Buat,
S. Salim,
S. Brough,
R. Demarco,
A. W. Graham,
M. Hamed,
J. R. Mullaney,
M. Romano,
C. Sifón,
M. Aravena,
J. A. Benavides,
I. Busà,
D. Donevski,
O. Dorey,
H. M. Hernandez-Toledo,
A. Nanni,
W. J. Pearson,
F. Pistis,
R. Ragusa,
G. Riccio,
J. Román
Abstract:
Large future sky surveys, such as the LSST, will provide optical photometry for billions of objects. This paper aims to construct a proxy for the far ultraviolet attenuation (AFUVp) from the optical data alone, enabling the rapid estimation of the star formation rate (SFR) for galaxies that lack UV or IR data. To mimic LSST observations, we use the deep panchromatic optical coverage of the SDSS Ph…
▽ More
Large future sky surveys, such as the LSST, will provide optical photometry for billions of objects. This paper aims to construct a proxy for the far ultraviolet attenuation (AFUVp) from the optical data alone, enabling the rapid estimation of the star formation rate (SFR) for galaxies that lack UV or IR data. To mimic LSST observations, we use the deep panchromatic optical coverage of the SDSS Photometric Catalogue DR~12, complemented by the estimated physical properties for the SDSS galaxies from the GALEX-SDSS-WISE Legacy Catalog (GSWLC) and inclination information obtained from the SDSS DR7. We restricted our sample to the 0.025-0.1 z-spec range and investigated relations among surface brightness, colours, and dust attenuation in the far UV range for star-forming galaxies obtained from the spectral energy distribution (SED). {Dust attenuation is best correlated with (u-r) colour and the surface brightness in the u band ($\rm μ_{u}$). We provide a dust attenuation proxy for galaxies on the star-forming main sequence, which can be used for the LSST or any other type of broadband optical survey. The mean ratio between the catalogue values of SFR and those estimated using optical-only SDSS data with the AFUVp prior calculated as $Δ$SFR=log(SFR$_{\tiny{\mbox{this work}}}$/SFR$_{\tiny{}\texttt{GSWLC}}$) is found to be less than 0.1~dex, while runs without priors result in an SFR overestimation larger than 0.3~dex. The presence or absence of theAFUVp has a negligible influence on the stellar mass estimation (with $Δ$M$_{star}$ in the range from 0 to $-0.15$ dex). Forthcoming deep optical observations of the LSST Deep Drilling Fields, which also have multi-wavelength data, will enable one to calibrate the obtained relation for higher redshift galaxies and, possibly, extend the study towards other types of galaxies, such as early-type galaxies off the main sequence.
△ Less
Submitted 1 February, 2024; v1 submitted 23 January, 2024;
originally announced January 2024.
-
Imposters among us: globular cluster kinematics and the halo mass of ultra-diffuse galaxies in clusters
Authors:
Jessica E. Doppel,
Laura V. Sales,
José A. Benavides,
Elisa Toloba,
Eric W. Peng,
Dylan Nelson,
Julio F. Navarro
Abstract:
The velocity dispersion of globular clusters (GCs) around ultra-diffuse galaxies (UDGs) in the Virgo cluster spans a wide range, including cases where GC kinematics suggest halos as massive as (or even more massive than) that of the Milky Way around these faint dwarfs. We analyze the catalogs of GCs derived in post-processing from the TNG50 cosmological simulation to study the GC system kinematics…
▽ More
The velocity dispersion of globular clusters (GCs) around ultra-diffuse galaxies (UDGs) in the Virgo cluster spans a wide range, including cases where GC kinematics suggest halos as massive as (or even more massive than) that of the Milky Way around these faint dwarfs. We analyze the catalogs of GCs derived in post-processing from the TNG50 cosmological simulation to study the GC system kinematics and abundance of simulated UDGs in galaxy groups and clusters. UDGs in this simulation reside exclusively in dwarf-mass halos with $M_{200} \sim 10^{11}$ M$_{\odot}$. When considering only GCs gravitationally bound to simulated UDGs, we find GCs properties that overlap well with several observational measurements for UDGs. In particular, no bias towards overly massive halos is inferred from the study of bound GCs, confirming that GCs are good tracers of UDG halo mass. However, we find that contamination by intra-cluster GCs may, in some cases, substantially increase velocity dispersion estimates when performing projected mock observations of our sample. We caution that targets with less than $10$ GC tracers are particularly prone to severe uncertainties.Measuring the stellar kinematics of the host galaxy should help confirm the unusually massive halos suggested by GC kinematics around some UDGs
△ Less
Submitted 6 September, 2023;
originally announced September 2023.
-
DS+: a method for the identification of cluster substructures
Authors:
Jose A. Benavides,
Andrea Biviano,
Mario G. Abadi
Abstract:
The study of cluster substructures is important for the determination of the cluster dynamical status, assembly history, and the evolution of cluster galaxies, and it allows to set of constraints on the nature of dark matter and cosmological parameters. We present and test DS+, a new method for the identification and characterization of group-sized substructures in clusters. Our new method is base…
▽ More
The study of cluster substructures is important for the determination of the cluster dynamical status, assembly history, and the evolution of cluster galaxies, and it allows to set of constraints on the nature of dark matter and cosmological parameters. We present and test DS+, a new method for the identification and characterization of group-sized substructures in clusters. Our new method is based on the projected positions and line-of-sight velocities of cluster galaxies, and it is an improvement and extension of the traditional method of Dressler & Shectman (1988). We test it on cluster-size cosmological halos extracted from the IllustrisTNG simulations, with virial masses $\rm{14 \lesssim \log (M_{200}/M_{\odot}) \lesssim 14.6}$, that contain on average $\sim 190$ galaxies. We also present an application of our method on a real data set, the Bullet cluster. DS+ is able to identify $\sim 80\%$ of real group galaxies as members of substructures, and at least 60\% of the galaxies assigned to substructures belong to real groups. The physical properties of the real groups are significantly correlated with those of the corresponding detected substructures, albeit with significant scatter, and overestimated on average. Application of the DS+ method to the Bullet cluster confirms the presence and main properties of the high-speed collision and identifies other substructures along the main cluster axis. DS+ proves to be a reliable method for the identification of substructures in clusters. The method is made freely available to the community as a Python code.
△ Less
Submitted 30 November, 2022;
originally announced December 2022.
-
Origin and evolution of ultra-diffuse galaxies in different environments
Authors:
Jose A. Benavides,
Laura V. Sales,
Mario. G. Abadi,
Federico Marinacci,
Mark Vogelsberger,
Lars Hernquist
Abstract:
We study the formation of ultra-diffuse galaxies (UDGs) using the cosmological hydrodynamical simulation TNG50 of the Illustris-TNG suite. We define UDGs as dwarf galaxies in the stellar mass range $\rm{7.5 \leq log (M_{\star} / M_{\odot}) \leq 9 }$ that are in the $5\%$ most extended tail of the simulated mass-size relation. This results in a sample of UDGs with half-mass radii…
▽ More
We study the formation of ultra-diffuse galaxies (UDGs) using the cosmological hydrodynamical simulation TNG50 of the Illustris-TNG suite. We define UDGs as dwarf galaxies in the stellar mass range $\rm{7.5 \leq log (M_{\star} / M_{\odot}) \leq 9 }$ that are in the $5\%$ most extended tail of the simulated mass-size relation. This results in a sample of UDGs with half-mass radii $\rm{r_{h \star } \gtrsim 2 \ kpc}$ and surface brightness between $\rm{24.5}$ and $\rm{28 \ mag \ arcsec^{-2}}$, similar to definitions of UDGs in observations. The large cosmological volume in TNG50 allows for a comparison of UDGs properties in different environments, from the field to galaxy clusters with virial mass $\rm{M_{200} \sim 2 \times 10^{14} ~ M_{\odot}}$. All UDGs in our sample have dwarf-mass haloes ($\rm{M_{200}\sim 10^{11} ~ M_{\odot} }$) and show the same environmental trends as normal dwarfs: field UDGs are star-forming and blue while satellite UDGs are typically quiescent and red. The TNG50 simulation predicts UDGs that populate preferentially higher spin haloes and more massive haloes at fixed $\rm{M_{\star}}$ compared to non-UDG dwarfs. This applies also to most satellite UDGs, which are actually ``born" UDGs in the field and infall into groups and clusters without significant changes to their size. We find, however, a small subset of satellite UDGs ($\lesssim 10 \%$) with present-day stellar size a factor $\geq 1.5$ larger than at infall, confirming that tidal effects, particularly in the lower mass dwarfs, are also a viable formation mechanism for some of these dwarfs, although subdominant in this simulation.
△ Less
Submitted 12 April, 2023; v1 submitted 15 September, 2022;
originally announced September 2022.
-
Preparing for low surface brightness science with the Vera C. Rubin Observatory: characterisation of tidal features from mock images
Authors:
G. Martin,
A. E. Bazkiaei,
M. Spavone,
E. Iodice,
J. C. Mihos,
M. Montes,
J. A. Benavides,
S. Brough,
J. L. Carlin,
C. A. Collins,
P. A. Duc,
F. A. Gómez,
G. Galaz,
H. M. Hernández-Toledo,
R. A. Jackson,
S. Kaviraj,
J. H. Knapen,
C. Martínez-Lombilla,
S. McGee,
D. O'Ryan,
D. J. Prole,
R. M. Rich,
J. Román,
E. A. Shah,
T. K. Starkenburg
, et al. (28 additional authors not shown)
Abstract:
Tidal features in the outskirts of galaxies yield unique information about their past interactions and are a key prediction of the hierarchical structure formation paradigm. The Vera C. Rubin Observatory is poised to deliver deep observations for potentially of millions of objects with visible tidal features, but the inference of galaxy interaction histories from such features is not straightforwa…
▽ More
Tidal features in the outskirts of galaxies yield unique information about their past interactions and are a key prediction of the hierarchical structure formation paradigm. The Vera C. Rubin Observatory is poised to deliver deep observations for potentially of millions of objects with visible tidal features, but the inference of galaxy interaction histories from such features is not straightforward. Utilising automated techniques and human visual classification in conjunction with realistic mock images produced using the NEWHORIZON cosmological simulation, we investigate the nature, frequency and visibility of tidal features and debris across a range of environments and stellar masses. In our simulated sample, around 80 per cent of the flux in the tidal features around Milky Way or greater mass galaxies is detected at the 10-year depth of the Legacy Survey of Space and Time (30-31 mag / sq. arcsec), falling to 60 per cent assuming a shallower final depth of 29.5 mag / sq. arcsec. The fraction of total flux found in tidal features increases towards higher masses, rising to 10 per cent for the most massive objects in our sample (M*~10^{11.5} Msun). When observed at sufficient depth, such objects frequently exhibit many distinct tidal features with complex shapes. The interpretation and characterisation of such features varies significantly with image depth and object orientation, introducing significant biases in their classification. Assuming the data reduction pipeline is properly optimised, we expect the Rubin Observatory to be capable of recovering much of the flux found in the outskirts of Milky Way mass galaxies, even at intermediate redshifts (z<0.2).
△ Less
Submitted 7 May, 2022; v1 submitted 15 March, 2022;
originally announced March 2022.
-
SCORPIO, a package for the visualization of galaxy pairs
Authors:
Jose A. Benavides,
Martín Chalela,
Juan B. Cabral,
Bruno O. Sánchez,
Sebastian Gurovich
Abstract:
We present the description of the project \texttt{SCORPIO}, a Python package for retrieving images and associated data of galaxy pairs based on their position, facilitating visual analysis and data collation of multiple archetypal systems. The code ingests information from SDSS, 2MASS, and WISE surveys based on the available bands and is designed for studies of galaxy pairs as natural laboratories…
▽ More
We present the description of the project \texttt{SCORPIO}, a Python package for retrieving images and associated data of galaxy pairs based on their position, facilitating visual analysis and data collation of multiple archetypal systems. The code ingests information from SDSS, 2MASS, and WISE surveys based on the available bands and is designed for studies of galaxy pairs as natural laboratories of multiple astrophysical phenomena such as tidal force deformation of galaxies, pressure gradient induced star formation regions, morphological transformation, to name a few.
△ Less
Submitted 8 November, 2021; v1 submitted 1 November, 2021;
originally announced November 2021.
-
Quiescent Ultra-diffuse galaxies in the field originating from backsplash orbits
Authors:
José A. Benavides,
Laura V. Sales,
Mario. G. Abadi,
Annalisa Pillepich,
Dylan Nelson,
Federico Marinacci,
Michael Cooper,
Ruediger Pakmor,
Paul Torrey,
Mark Vogelsberger,
Lars Hernquist
Abstract:
Ultra-diffuse galaxies (UDGs) are the lowest-surface brightness galaxies known, with typical stellar masses of dwarf galaxies but sizes similar to larger galaxies like the Milky Way. The reason for their extended sizes is debated, with suggested internal processes like angular momentum, feedback or mergers versus external mechanisms or a combination of both. Observationally, we know that UDGs are…
▽ More
Ultra-diffuse galaxies (UDGs) are the lowest-surface brightness galaxies known, with typical stellar masses of dwarf galaxies but sizes similar to larger galaxies like the Milky Way. The reason for their extended sizes is debated, with suggested internal processes like angular momentum, feedback or mergers versus external mechanisms or a combination of both. Observationally, we know that UDGs are red and quiescent in groups and clusters while their counterparts in the field are blue and star-forming. This dichotomy suggests environmental effects as main culprit. However, this scenario is challenged by recent observations of isolated quiescent UDGs in the field. Here we use $Λ$CDM cosmological hydrodynamical simulation to show that isolated quenched UDGs are formed as backsplash galaxies that were once satellites of another galactic, group or cluster halo but are today a few Mpc away from them. These interactions, albeit brief, remove the gas and tidally strip the outskirts of the dark matter haloes of the now quenched seemingly-isolated UDGs, which are born as star-forming field UDGs occupying dwarf-mass dark matter haloes. Quiescent UDGs may therefore be found in non-negligible numbers in filaments and voids, bearing the mark of past interactions as stripped outer haloes devoid of dark matter and gas compared to dwarfs with similar stellar content.
△ Less
Submitted 3 September, 2021;
originally announced September 2021.
-
A unified scenario for the origin of spiral and elliptical galaxy structural scaling laws
Authors:
I. Ferrero,
J. F. Navarro,
M. G. Abadi,
J. A. Benavides,
D. Mast
Abstract:
Elliptical (E) and spiral (S) galaxies follow tight, but different, scaling laws that link their stellar masses, radii, and characteristic velocities. Mass and velocity, for example, scale tightly in spirals with little dependence on galaxy radius (the "Tully-Fisher relation"; TFR). On the other hand, ellipticals appear to trace a 2D surface in size-mass-velocity space (the "Fundamental Plane"; FP…
▽ More
Elliptical (E) and spiral (S) galaxies follow tight, but different, scaling laws that link their stellar masses, radii, and characteristic velocities. Mass and velocity, for example, scale tightly in spirals with little dependence on galaxy radius (the "Tully-Fisher relation"; TFR). On the other hand, ellipticals appear to trace a 2D surface in size-mass-velocity space (the "Fundamental Plane"; FP). Over the years, a number of studies have attempted to understand these empirical relations, usually in terms of variations of the virial theorem for E galaxies and in terms of the scaling relations of dark matter halos for spirals. We use Lambda cold dark matter (LCDM) cosmological hydrodynamical simulations to show that the scaling relations of both ellipticals and spirals arise as the result of (i) a tight galaxy mass-dark halo mass relation and (ii) the self-similar mass profile of cold dark matter halos. In this interpretation, E and S galaxies of a given stellar mass inhabit halos of similar masses, and their different scaling laws result from the varying amounts of dark matter enclosed within their luminous radii. This scenario suggests a new galaxy distance indicator applicable to galaxies of all morphologies and provides simple and intuitive explanations for long-standing puzzles, such as why the TFR is independent of surface brightness, or what causes the "tilt" in the FP. Our results provide strong support for the predictions of LCDM in the strongly non-linear regime, as well as guidance for further improvements to cosmological simulations of galaxy formation.
△ Less
Submitted 3 March, 2021; v1 submitted 8 September, 2020;
originally announced September 2020.
-
Accretion of Galaxy Groups into Galaxy Clusters
Authors:
José. A. Benavides,
Laura. V. Sales,
Mario. G. Abadi
Abstract:
We study the role of group infall in the assembly and dynamics of galaxy clusters in $Λ$CDM. We select $10$ clusters with virial mass $M_{\rm 200} \sim 10^{14} \, M_{\odot}$ from the cosmological hydrodynamical simulation Illustris and follow their galaxies with stellar mass $M_{\star} \geq 1.5 \times 10^8 \, M_{\odot}$. A median of $\sim 38\%$ of surviving galaxies at $z=0$ are accreted as part o…
▽ More
We study the role of group infall in the assembly and dynamics of galaxy clusters in $Λ$CDM. We select $10$ clusters with virial mass $M_{\rm 200} \sim 10^{14} \, M_{\odot}$ from the cosmological hydrodynamical simulation Illustris and follow their galaxies with stellar mass $M_{\star} \geq 1.5 \times 10^8 \, M_{\odot}$. A median of $\sim 38\%$ of surviving galaxies at $z=0$ are accreted as part of groups and did not infall directly from the field, albeit with significant cluster-to-cluster scatter. The evolution of these galaxy associations is quick, with observational signatures of their common origin eroding rapidly in $1$-$3$ Gyr after infall. Substructure plays a dominant role in fostering the conditions for galaxy mergers to happen, even within the cluster environment. Integrated over time, we identify (per cluster) an average of $17 \pm 6$ mergers that occur in infalling galaxy associations, of which $7 \pm 3$ occur well within the virial radius of their cluster hosts. The number of mergers show large dispersion from cluster to cluster, with our most massive system having $42$ mergers above our mass cut-off. These mergers, which are typically gas rich for dwarfs and a combination of gas rich and gas poor for $M_{\star} \sim 10^{11} \, M_{\odot}$, may contribute significantly within $Λ$CDM to the formation of specific morphologies, such as lenticulars (S0) and blue compact dwarfs in groups and clusters.
△ Less
Submitted 2 September, 2020; v1 submitted 11 May, 2020;
originally announced May 2020.
