Abstract visibility
Abstract text size

Papers for Thursday, Jan 06 2022

Papers with local authors

Blake Chellew, Timothy D. Brandt, Brandon S. Hensley, Bruce T. Draine, Eve Matthaey

16 pages, 8 figures. Submitted to ApJ on 12/21/21

1 vote
Paper 17 — arXiv:2201.01378
1 vote
Paper 17 — arXiv:2201.01378

We present a spectrum of the diffuse Galactic light (DGL) between 3700 and 10,000 A, obtained by correlating optical sky intensity with far-infrared dust emission. We use nearly 250,000 blank-sky spectra from BOSS/SDSS-III together with IRIS-reprocessed maps from the IRAS satellite. The larger sample size compared to SDSS-II results in a factor-of-two increase in signal to noise. We combine these data sets with a model for the optical/far-infrared correlation that accounts for self-absorption by dust. The spectral features of the DGL agree remarkably well with features present in stellar spectra. There is evidence for a difference in the DGL continuum between the regions covered by BOSS in the northern and southern Galactic hemisphere. We interpret the difference at red wavelengths as the result of a difference in stellar populations, with mainly old stars in both regions but a higher fraction of young stars in the south. There is also a broad excess in the southern DGL spectrum over the prediction of a simple radiative transfer model, without a clear counterpart in the north. We interpret this excess, centered at ~6500 A, as evidence for luminescence in the form of extended red emission (ERE). The observed strength of the 4000 A break indicates that at most ~7% of the dust-correlated light at 4000 A can be due to blue luminescence. Our DGL spectrum provides constraints on dust scattering and luminescence independent of measurements of extinction.

Francisco Villaescusa-Navarro, Shy Genel, Daniel Anglés-Alcázar, Lucia A. Perez, Pablo Villanueva-Domingo, Digvijay Wadekar, Helen Shao, Faizan G. Mohammad, Sultan Hassan, Emily Moser, Erwin T. Lau, Luis Fernando Machado Poletti Valle, Andrina Nicola, Leander Thiele, Yongseok Jo, Oliver H. E. Philcox, Benjamin D. Oppenheimer, Megan Tillman, ChangHoon Hahn, Neerav Kaushal, Alice Pisani, Matthew Gebhardt, Ana Maria Delgado, Joyce Caliendo, Christina Kreisch, Kaze W.K. Wong, William R. Coulton, Michael Eickenberg, Gabriele Parimbelli, Yueying Ni, Ulrich P. Steinwandel, Valentina La Torre, Romeel Dave, Nicholas Battaglia, Daisuke Nagai, David N. Spergel, Lars Hernquist, Blakesley Burkhart, Desika Narayanan, Benjamin Wandelt, Rachel S. Somerville, Greg L. Bryan, Matteo Viel, Yin Li, et al.

18 pages, 3 figures. More than 350 Tb of data from thousands of simulations publicly available at this https URL

0 votes
Paper 2 — arXiv:2201.01300
0 votes
Paper 2 — arXiv:2201.01300

The Cosmology and Astrophysics with MachinE Learning Simulations (CAMELS) project was developed to combine cosmology with astrophysics through thousands of cosmological hydrodynamic simulations and machine learning. CAMELS contains 4,233 cosmological simulations, 2,049 N-body and 2,184 state-of-the-art hydrodynamic simulations that sample a vast volume in parameter space. In this paper we present the CAMELS public data release, describing the characteristics of the CAMELS simulations and a variety of data products generated from them, including halo, subhalo, galaxy, and void catalogues, power spectra, bispectra, Lyman-$\alpha$ spectra, probability distribution functions, halo radial profiles, and X-rays photon lists. We also release over one thousand catalogues that contain billions of galaxies from CAMELS-SAM: a large collection of N-body simulations that have been combined with the Santa Cruz Semi-Analytic Model. We release all the data, comprising more than 350 terabytes and containing 143,922 snapshots, millions of halos, galaxies and summary statistics. We provide further technical details on how to access, download, read, and process the data at \url{https://camels.readthedocs.io}.

Digvijay Wadekar, Leander Thiele, Francisco Villaescusa-Navarro, J. Colin Hill, David N. Spergel, Miles Cranmer, Nicholas Battaglia, Daniel Anglés-Alcázar, Lars Hernquist, Shirley Ho

11+6 pages, 8+3 figures. The code and data associated with this paper will be uploaded upon the acceptance of this paper

0 votes
Paper 5 — arXiv:2201.01305
0 votes
Paper 5 — arXiv:2201.01305

Complex systems (stars, supernovae, galaxies, and clusters) often exhibit low scatter relations between observable properties (e.g., luminosity, velocity dispersion, oscillation period, temperature). These scaling relations can illuminate the underlying physics and can provide observational tools for estimating masses and distances. Machine learning can provide a systematic way to search for new scaling relations (or for simple extensions to existing relations) in abstract high-dimensional parameter spaces. We use a machine learning tool called symbolic regression (SR), which models the patterns in a given dataset in the form of analytic equations. We focus on the Sunyaev-Zeldovich flux$-$cluster mass relation ($Y_\mathrm{SZ}-M$), the scatter in which affects inference of cosmological parameters from cluster abundance data. Using SR on the data from the IllustrisTNG hydrodynamical simulation, we find a new proxy for cluster mass which combines $Y_\mathrm{SZ}$ and concentration of ionized gas ($c_\mathrm{gas}$): $M \propto Y_\mathrm{conc}^{3/5} \equiv Y_\mathrm{SZ}^{3/5} (1-A\, c_\mathrm{gas})$. $Y_\mathrm{conc}$ reduces the scatter in the predicted $M$ by $\sim 20-30$% for large clusters ($M\gtrsim 10^{14}\, h^{-1} \, M_\odot$) at both high and low redshifts, as compared to using just $Y_\mathrm{SZ}$. We show that the dependence on $c_\mathrm{gas}$ is linked to cores of clusters exhibiting larger scatter than their outskirts. Finally, we test $Y_\mathrm{conc}$ on clusters from simulations of the CAMELS project and show that $Y_\mathrm{conc}$ is robust against variations in cosmology, astrophysics, subgrid physics, and cosmic variance. Our results and methodology can be useful for accurate multiwavelength cluster mass estimation from current and upcoming CMB and X-ray surveys like ACT, SO, SPT, eROSITA and CMB-S4.

A.F.A. Bott, L. Chen, P. Tzeferacos, C.A.J. Palmer, A.R. Bell, R. Bingham, A. Birkel, D.H. Froula, J. Katz, M.W. Kunz, C.-K. Li, H-S. Park, R. Petrasso, J.S. Ross, B. Reville, D. Ryu, F.H. Séguin, T.G. White, A.A. Schekochihin, D.Q. Lamb, G. Gregori

13 pages, 12 figures

0 votes
Paper 6 — arXiv:2201.01705
0 votes
Paper 6 — arXiv:2201.01705

It has recently been demonstrated experimentally that a turbulent plasma created by the collision of two inhomogeneous, asymmetric, weakly magnetised laser-produced plasma jets can generate strong stochastic magnetic fields via the small-scale turbulent dynamo mechanism, provided the magnetic Reynolds number of the plasma is sufficiently large. In this paper, we compare such a plasma with one arising from two pre-magnetised plasma jets whose creation is identical save for the addition of a strong external magnetic field imposed by a pulsed magnetic field generator (`MIFEDS'). We investigate the differences between the two turbulent systems using a Thomson-scattering diagnostic, X-ray self-emission imaging and proton radiography. The Thomson-scattering spectra and X-ray images suggest that the presence of the external magnetic field has a limited effect on the plasma dynamics in the experiment. While the presence of the external magnetic field induces collimation of the flows in the colliding plasma jets and the initial strengths of the magnetic fields arising from the interaction between the colliding jets are significantly larger as a result of the external field, the energy and morphology of the stochastic magnetic fields post-amplification are indistinguishable. We conclude that, for turbulent laser-plasmas with super-critical magnetic Reynolds numbers, the dynamo-amplified magnetic fields are determined by the turbulent dynamics rather than the seed fields and modest changes in the initial flow dynamics of the plasma, a finding consistent with theoretical expectations and simulations of turbulent dynamos.

Shriya Soma, Lingxiao Wang, Shuzhe Shi, Horst Stöcker, Kai Zhou

12 pages, 11 figures

0 votes
Paper 8 — arXiv:2201.01756
0 votes
Paper 8 — arXiv:2201.01756

The equation of state (EoS) of the strongly interacting cold and ultra-dense matter remains a major challenge in the field of nuclear physics. With the advancements in measurements of neutron star masses, radii, and tidal deformabilities from electromagnetic and gravitational wave observations, neutron stars play an important role in constraining the ultra-dense matter EoS. In this work, we present a novel method that exploits deep learning techniques to reconstruct the neutron star EoS from mass-radius (M-R) observations. We employ neural networks (NNs) to represent the EoS in a model-independent way, within the range of 1-7.4 times the nuclear saturation density. In an unsupervised manner, we implement the Automatic Differentiation (AD) framework to optimize the EoS, as to yield through TOV equations an M-R curve that best fits the observations. We demonstrate it in rebuilding the EoS on mock data, i.e., mass-radius pairs derived from a generated set of polytropic EoSs. The results show that it is possible to reconstruct the EoS with reasonable accuracy, using just 11 mock M-R pairs observations, which is close to the current number of observations.

All other papers

Clancy W. James

15 pages, 1 table, 8 figures, submitted to Physical Review D

The nature of the radio-wave radiation generated by particle cascades in both the Earth's atmosphere and dense media such as ice has, historically, been much debated. This situation changed in the early 2010's, with the community converging on the common terminology of "geomagnetic" and "Askaryan" radiation to describe the two emission mechanisms. However, this convergence arose from discussions at various conferences and workshops, and was ultimately reached through agreement between simulation codes and experimental measurements. In this article therefore, I use relatively simple geometrical arguments, and a minimum of calculations based on single particle tracks, to explain the nature of radiation from extensive air showers (EAS) and cascades in dense media such as ice. I identify well-determined frequency regimes where the radiation from the Askaryan effect will be bremsstrahlung-like and Cherenkov-like, being respectively below/above 1 GHz in EAS and 100 MHz in dense media; and where geomagnetic emission will be transverse-current-like and where it will resemble synchrotron radiation, respectively below/above a few GHz in EAS, depending on the height of cascade development. I suggest how these transitions in the nature of the emission may be experimentally observed.

Vijay Varma, Sylvia Biscoveanu, Tousif Islam, Feroz H. Shaik, Carl-Johan Haster, Maximiliano Isi, Will M. Farr, Scott E. Field, Salvatore Vitale

5 pages. Visualizations available at this https URL

The final black hole left behind after a binary black hole merger can attain a recoil velocity, or a "kick", reaching values up to 5000 km/s. This phenomenon has important implications for gravitational wave astronomy, black hole formation scenarios, testing general relativity, and galaxy evolution. We consider the gravitational wave signal from the binary black hole merger GW200129_065458 (henceforth referred to as GW200129), which has been shown to exhibit strong evidence of orbital precession. Using numerical relativity surrogate models, we constrain the kick velocity of GW200129 to $v_f \sim 1542^{+747}_{-1098}$ km/s or $v_f \gtrsim 698$ km/s (one-sided limit), at 90\% credibility. This marks the first identification of a large kick velocity for an individual gravitational wave event. Given the kick velocity of GW200129, we estimate that there is a less than $0.48\%$ ($7.7\%$) probability that the remnant black hole after the merger would be retained by globular (nuclear star) clusters. Finally, we show that kick effects are not expected to cause biases in ringdown tests of general relativity for this event, although this may change in the future with improved detectors.

Christina C. Williams, Stacey Alberts, Justin S. Spilker, Allison G. Noble, Mauro Stefanon, Christopher N. A. Willmer, Rachel Bezanson, Desika Narayanan, Katherine E. Whitaker

Submitted to AAS Journals on December 15, 2021. Comments welcome

We report the serendipitous discovery of an overdensity of CO emitters in an X-ray-identified cluster (Log$_{10}$M$_{\rm halo}/M_{\odot}\sim13.6$ at z=1.3188) using ALMA. We present spectroscopic confirmation of 6 new cluster members exhibiting CO(2-1) emission, adding to 2 existing optical/IR spectroscopic members undetected in CO. This is the lowest mass cluster to date at z>1 with molecular gas measurements, bridging the observational gap between galaxies in the more extreme, well-studied clusters (Log$_{10}$~M$_{\rm halo}/M_{\odot}\gtrsim14$) and those in group or field environments at cosmic noon. The CO sources are concentrated on the sky (within ~1-arcmin diameter) and phase space analysis indicates the gas resides in galaxies already within the cluster environment. We find that CO sources sit in similar phase space as CO-rich galaxies in more massive clusters at similar redshifts (have similar accretion histories) while maintaining field-like molecular gas reservoirs, compared to scaling relations. This work presents the deepest CO survey to date in a galaxy cluster at z>1, uncovering gas reservoirs down to M$_{\rm H_{2}}>1.6\times10^{10}$M$_{\odot}$ (5$\sigma$ at 50% primary beam). Our deep limits rule out the presence of gas content in excess of the field scaling relations; however, combined with literature CO detections, cluster gas fractions in general appear systematically high, on the upper envelope or above the field. This study is the first demonstration that low mass clusters at z~1-2 can host overdensities of CO emitters with surviving gas reservoirs, in line with the prediction that quenching is delayed after first infall while galaxies consume the gas bound to the disk.

R. L. S. Frisbie, M. Donahue, G. M. Voit, K. Lakhchaura, N. Werner, M. Sun

Accepted to ApJ, 16 pages, 9 figures

The Voit et al. (2020) black hole feedback valve model predicts relationships between stellar velocity dispersion and atmospheric structure among massive early-type galaxies. In this work, we test that model using the Chandra archival sample of 49 early-type galaxies from Lakhchaura et al. (2018). We consider relationships between stellar velocity dispersion and entropy profile slope, multiphase gas extent, and the ratio of cooling time to freefall time. We also define subsamples based on data quality and entropy profile properties that clarify those relationships and enable more specific tests of the model predictions. We find that the atmospheric properties of early-type galaxies generally align with the predictions of the Voit et al. (2020) model, in that galaxies with greater stellar velocity dispersion tend to have radial profiles of pressure, gas density, and entropy with steeper slopes and less extended multiphase gas. Quantitative agreement with the model predictions improves when the sample is restricted to have low central entropy and stellar velocity dispersion of between 220 and 300 km/s.

Stacey Alberts, Jéa Adams, Benjamin Gregg, Alexandra Pope, Christina C. Williams, Peter R. M. Eisenhardt

Accepted for publication in ApJ

We present the average gas properties derived from ALMA Band 6 dust continuum imaging of 126 massive (log $M_{\star} / M_{\odot} \gtrsim 10.5$), star-forming cluster galaxies across 11 galaxy clusters at $z=1-1.75$. Using stacking analysis on the ALMA images, combined with UV-far-infrared data, we quantify the average infrared SEDs and gas properties (molecular gas masses, $M_{\rm mol}$; gas depletion timescales, $\tau_{\rm depl}$; and gas fractions, f$_{\rm gas}$) as a function of cluster-centric radius and properties including stellar mass and distance from the Main Sequence. We find a significant dearth in the ALMA fluxes relative to that expected in the field $-$ with correspondingly low $M_{\rm mol}$ and f$_{\rm gas}$ and short $\tau_{\rm depl}$ $-$ with weak or no dependence on cluster-centric radius out to twice the virial radius. The Herschel+ALMA SEDs indicate warmer dust temperatures ($\sim36-38$ K) than coeval field galaxies ($\sim30$ K). We perform a thorough comparison of the cluster galaxy gas properties to field galaxies, finding deficits of 2-3x, 3-4x, and 2-4x in $M_{\rm mol}$, $\tau_{\rm depl}$, and f$_{\rm gas}$ compared to coeval field stacks and larger deficits compared to field scaling relations built primarily on detections. The cluster gas properties derived here are comparable with stacking analyses in (proto-)clusters in the literature and at odds with findings of field-like $\tau_{\rm depl}$ and enhanced f$_{\rm gas}$ reported using CO and dust continuum detections. Our analysis suggests that environment has considerable impact on gas properties out to large radii, in good agreement with cosmological simulations which project gas depletion begins beyond the virial radius and largely completes by first passage of the cluster core.

Nicolas F. Martin, Kim A. Venn, David S. Aguado, Else Starkenburg, Jonay I. González Hernández, Rodrigo A. Ibata, Piercarlo Bonifacio, Elisabetta Caffau, Federico Sestito, Anke Arentsen, Carlos Allende Prieto, Raymond G. Carlberg, Sébastien Fabbro, Morgan Fouesneau, Vanessa Hill, Pascale Jablonka, Georges Kordopatis, Carmela Lardo, Khyati Malhan, Lyudmila I. Mashonkina, Alan W. McConnachie, Julio F. Navarro, Rubén Sánchez Janssen, Guillaume F. Thomas, Zhen Yuan, Alessio Mucciarelli

14 pages, 6 figures, 5 tables. Nature, accepted

Stellar ejecta gradually enrich the gas out of which subsequent stars form, making the least chemically enriched stellar systems direct fossils of structures formed in the early universe. Although a few hundred stars with metal content below one thousandth of the solar iron content are known in the Galaxy, none of them inhabit globular clusters, some of the oldest known stellar structures. These show metal content of at least ~0.2 percent of the solar metallicity ([Fe/H] > -2.7). This metallicity floor appears universal and it has been proposed that proto-galaxies that merge into the galaxies we observe today were simply not massive enough to form clusters that survived to the present day. Here, we report the discovery of a stellar stream, C-19, whose metallicity is less than 0.05 per cent the solar metallicity ([Fe/H]=-3.38 +/- 0.06 (stat.) +/- 0.20 (syst.)). The low metallicity dispersion and the chemical abundances of the C-19 stars show that this stream is the tidal remnant of the most metal-poor globular cluster ever discovered, and significantly below the purported metallicity floor: clusters with significantly lower metallicities than observed today existed in the past and contributed their stars to the Milky Way halo.

Nicolas F. Martin, Rodrigo A. Ibata, Else Starkenburg, Zhen Yuan, Khyati Malhan, Michele Bellazzini, Akshara Viswanathan, David Aguado, Anke Arentsen, Piercarlo Bonifacio, Ray Carlberg, Jonay I. González Hernández, Vanessa Hill, Pascale Jablonka, Georges Kordopatis, Carmela Lardo, Alan W. McConnachie, Julio Navarro, Rubén Sánchez-Janssen, Federico Sestito, Guillaume F. Thomas, Kim A. Venn, Sara Vitali, Karina T. Voggel

25 pages, 29 figures, submitted to MNRAS

We use the photometric metallicities provided by the panoramic Pristine survey to study the veracity and derive the metallicities of the numerous stellar streams found by Ibata et al. (2021) from the application of the STREAMFINDER algorithm to the Gaia EDR3 data. All 26 streams present in Pristine show a clear metallicity distribution function, which provides an independent check of the reality of these structures, supporting the reliability of STREAMFINDER in finding streams and the power of Pristine to measure precise metallicities. We further present 6 candidate structures with coherent phase-space and metallicity signals that are very likely streams. The majority of studied streams are very metal-poor (14 structures with [Fe/H]<-2.0) and include 3 extremely metal-poor systems ([Fe/H]<-3.0; C-11, C-19, and C-20). These streams could be the closest debris of low-luminosity dwarf galaxies or may have originated from globular clusters of significantly lower metallicity than any known current Milky Way globular cluster. Our study shows that the promise of the Gaia data for Galactic Archeology studies can be substantially strengthened by quality photometric metallicities, allowing us to peer back into the earliest epochs of the formation of our Galaxy and its stellar halo constituents.

Gemma González-Torà, Miguel A. Urbaneja, Norbert Przybilla, Stefan Dreizler, Martin M. Roth, Sebastian Kamann, Norberto Castro

8 pages, 9 figures, 2 tables. Accepted for publication in A&A

A quantitative spectral analysis of BA-type supergiants and bright giants in an inner spiral arm region of the nearby spiral galaxy NGC 300 is presented, based on observations with the Multi Unit Spectroscopic Explorer (MUSE) on the European Southern Obsevatory, Very Large Telescope (ESO, VLT). The flux-weighted gravity-luminosity relationship (FGLR), a stellar spectroscopic distance determination method for galaxies, is extended towards stars at lower luminosities. Point spread function fitting 3D spectroscopy was performed with PampelMUSE on the datacube. The 16 stars with the highest signal-to-noise ratios ($S/N s$) are classified with regard to their spectral type and luminosity class using Galactic templates. They were analysed using hybrid non-local thermodynamic equilibrium (non-LTE) model spectra to fit the strongest observed hydrogen, helium, and metal lines in the intermediate-resolution spectra. Supplemented by photometric data, this facilitates fundamental stellar parameters and interstellar reddening which have yet to be determined. Effective temperatures, surface gravities, reddening $E(B-V)$, bolometric magnitudes and luminosities, as well as radii and masses are presented for the sample stars. The majority of the objects follow the FGLR as established from more luminous BA-type supergiants in NGC 300. An increase in the scatter in the flux-weighted gravity-luminosity plane is observed at these lower luminosities, which is in line with predictions from population synthesis models.

Jonathan Holdship, Serena Viti

9 pages, 2 figures. Accepted for publication in A&A. Results table at this https URL

Molecular line emission is a powerful probe of the physical conditions of astrophysical objects but can be complex to model, and it is often unclear which transitions would be the best targets for observers who wish to constrain a given parameter. We therefore produce a list of molecular species for which the gas history can be ignored, removing a major modelling complexity. We then determine the best of these species to observe when attempting to constrain various physical parameters. To achieve this, we use a large set of chemical models with different chemical histories to determine which species have abundances at 1 MYr that are insensitive to the initial conditions. We then use radiative transfer modelling to produce the intensity of every transition of these molecules. We finally compute the mutual information between the physical parameters and all transitions and transition ratios in order to rank their usefulness in determining the value of a given parameter. We find 48 species that are insensitive to the chemical history of the gas, 23 of which have collisional data available. We produce a ranked list of all the transitions and ratios of these species using their mutual information with various gas properties. We show mutual information is an adequate measure of how well a transition can constrain a physical parameter by recovering known probes and demonstrating that random forest regression models become more accurate predictors when high-scoring features are included. Therefore, this list can be used to select target transitions for observations in order to maximize knowledge about those physical parameters.

Allison Youngblood, J. Sebastian Pineda, Thomas Ayres, Kevin France, Jeffrey L. Linsky, Brian E. Wood, Seth Redfield, Joshua E. Schlieder

17 pages, 6 figures, accepted to ApJ

Observations of H I Lyman alpha, the brightest UV emission line of late-type stars, are critical for understanding stellar chromospheres and transition regions, modeling photochemistry in exoplanet atmospheres, and measuring the abundances of neutral hydrogen and deuterium in the interstellar medium. Yet, Lyman alpha observations are notoriously challenging due to severe attenuation from interstellar gas, hindering our understanding of this important emission line's basic morphology. We present high-resolution far- and near-UV spectroscopy of five G, K, and M dwarfs with radial velocities large enough to Doppler shift the stellar Lyman alpha emission line away from much of the interstellar attenuation, allowing the line core to be directly observed. We detect self-reversal in the Lyman alpha emission line core for all targets, and we show that the self-reversal depth decreases with increasing surface gravity. Mg II self-reversed emission line profiles provide some useful information to constrain the Lyman alpha line core, but the differences are significant enough that Mg II cannot be used directly as an intrinsic Lyman alpha template during reconstructions. We show that reconstructions that neglect self-reversal could overestimate intrinsic Lyman alpha fluxes by as much as 60%-100% for G and K dwarfs and 40%-170% for M dwarfs. The five stars of our sample have low magnetic activity and sub-solar metallicity; a larger sample size is needed to determine how sensitive these results are to these factors.

Jaehyun Lee, Taysun Kimm, Jérémy Blaizot, Harley Katz, Wonki Lee, Yun-Kyeong Sheen, Julien Devriendt, Adrianne Slyz

18 pages, 15 figures; Submitted to ApJ

We investigate the formation of jellyfish galaxies using radiation-hydrodynamic simulations of gas-rich dwarf galaxies with a multi-phase interstellar medium (ISM). We find that the ram-pressure-stripped (RPS) ISM is the dominant source of molecular clumps in the near wake within 10 kpc from the galactic plane, while in-situ formation is the major channel for dense gas in the distant tail of the gas-rich galaxy. Only 20% of the molecular clumps in the near wake originate from the intracluster medium (ICM); however, the fraction reaches 50% in the clumps located at 80 kpc from the galactic center since the cooling time of the RPS gas tends to be short due to the ISM-ICM mixing ($\lesssim$ 10 Myr). The tail region exhibits a star formation rate of $0.001-0.01 {\rm M_{\odot}} {\rm yr^{-1}}$, and most of the tail stars are born in the stripped wake within 10 kpc from the galactic plane. These stars induce bright H$\alpha$ blobs in the tail, while H$\alpha$ tails fainter than $6\times10^{38}\,{\rm erg\,s^{-1}\,kpc^{-2}}$ are mostly formed via collisional radiation and heating due to mixing. We also find that the stripped tails have intermediate X-ray to H$\alpha$ surface brightness ratios ($0.5\lesssim F_{\rm X}/F_{\rm H\alpha} \lesssim 10$), compared to the ISM ($\lesssim 0.5$) or pure ICM ($\gg10$). Our results suggest that jellyfish features emerge when the ISM from gas-rich galaxies is stripped by strong ram pressure, mixes with the ICM, and enhances the cooling in the tail.

William T.P. Denman, Samantha K. Trumbo, Michael E. Brown

11 pages, 8 figures, Published in conjunction with two other manuscripts, Planetary Science Journal, in press

Europa's leading hemisphere chaos regions have a spectral feature at 450 nm that has been attributed to absorption by crystal defects in irradiated sodium chloride, known as F-centers. Some discrepancies exist between the laboratory data of irradiated sodium chloride and the observations, including a $\sim$10 nm shift in central wavelength of the F-center band and the lack of the prominent 720 nm absorption on Europa from M-centers, which result from the coalescence of pairs of F-centers. Here, we perform irradiation experiments on sodium chloride in an attempt to understand these discrepancies. We show that careful control of the temperature of the sample at a temperature of 120 K yields F-centers with an absorption wavelength comparable to that of Europa. In addition, we measure the effect of photobleaching -- the destruction of F-centers by photons -- and show that at the energetic particle and photon flux on Europa, an equilibrium will be reached where only a modest F-center absorption develops. The density of F-centers never reaches high enough values for the creation of secondary M-centers. Our experiments predict that F-centers grow during the night on Europa in the absence of photobleaching and then partially decay during the daytime. We show observations from the Hubble Space Telescope consistent with this prediction. All observations of the 450 nm F-center on Europa are now consistent with laboratory measurements of sodium chloride, confirming the presence of this salt on Europa.

Samantha K. Trumbo, Tracy M. Becker, Michael E. Brown, William T. P. Denman, Philippa Molyneux, Amanda Hendrix, Kurt D. Retherford, Lorenz Roth, Juan Alday

6 pages, 2 figures, Accepted to the Planetary Science Journal

Recent visible-wavelength observations of Europa's surface obtained with the Hubble Space Telescope revealed the presence of an absorption feature near 450 nm that appears spatially correlated with leading-hemisphere chaos terrain. This feature was interpreted to reflect the presence of irradiated sodium chloride ultimately sourced from the interior. Here, we use ultraviolet spectra also collected with the Hubble Space Telescope to detect an additional previously unseen absorption near 230 nm, which spatially correlates with the 450 nm feature and with the same leading-hemisphere chaos terrain. We find that the new ultraviolet feature is also well-matched by irradiated sodium chloride at Europa-like conditions. Such confirmation of sodium chloride within geologically young regions has important implications for Europa's subsurface composition.

Michael E. Brown, William T.P. Denman, Samantha K. Trumbo

Planetary Science Journal, in press

Recent observations from the Hubble Space Telescope show a mid-UV absorption feature localized to leading hemisphere chaos regions on Europa. The same regions were previously found to have a visible absorption at 450 nm that was attributed to the presence of irradiated NaCl. The lack of any additional diagnostic absorptions for NaCl in the visible spectrum of these terrains made confirmation of this identification difficult. Here we use laboratory experiments to show that NaCl irradiated at Europa's surface temperatures develops an absorption at $\sim$220 nm consistent with the new detection in Europa's mid-UV spectrum, strongly supporting the NaCl identification. Irradiated NaCl in leading-hemisphere chaos terrain would suggest that sodium and chlorine are important components of Europa's subsurface ocean.

Nearby core-collapse supernovae (CCSNe) are powerful multi-messenger sources for gravitational-wave, neutrino and electromagnetic telescopes as they emit gravitational waves in the ideal frequency band for ground based detectors. Once a CCSN gravitational-wave signal is detected, we will need to determine the parameters of the signal, and understand how those parameters relate to the source's explosion, progenitor and remnant properties. This is a challenge due to the stochastic nature of CCSN explosions, which is imprinted on their time series gravitational waveforms. In this paper, we perform Bayesian parameter estimation of CCSN signals using an asymmetric chirplet signal model to represent the dominant high-frequency mode observed in spectrograms of CCSN gravitational-wave signals. We use design sensitivity Advanced LIGO noise and CCSN waveforms from four different hydrodynamical supernova simulations with a range of different progenitor stars. We determine how well our model can reconstruct time-frequency images of the emission modes, and show how well we can determine parameters of the signal such as the frequency, amplitude, and duration. We show how the parameters of our signal model may allow us to place constraints on the proto-neutron star mass and radius, the turbulent kinetic energy onto the proto-neutron star, and the time of shock revival.

Hsi-An Pan, Eva Schinnerer, Annie Hughes, Adam Leroy, Brent Groves, Ashley Thomas Barnes, Francesco Belfiore, Frank Bigiel, Guillermo A. Blanc, Yixian Cao, Melanie Chevance, Enrico Congiu, Daniel A. Dale, Cosima Eibensteiner, Eric Emsellem, Christopher M. Faesi, Simon C. O. Glover, Kathryn Grasha, Cinthya N. Herrera, I-Ting Ho, Ralf S. Klessen, J. M. Diederik Kruijssen, Philipp Lang, Daizhong Liu, Rebecca McElroy, Sharon E. Meidt, Eric J. Murphy, Jerome Pety, Miguel Querejeta, Alessandro Razza, Erik Rosolowsky, Toshiki Saito, Francesco Santoro, Andreas Schruba, Jiayi Sun, Neven Tomicic, Antonio Usero, Dyas Utomo, Thomas G. Williams

60 pages, 22 figures, 9 tables, accepted for publication in The Astrophysical Journal

The relative distribution of molecular gas and star formation in galaxies gives insight into the physical processes and timescales of the cycle between gas and stars. In this work, we track the relative spatial configuration of CO and H$\alpha$ emission at high resolution in each of our galaxy targets, and use these measurements to quantify the distributions of regions in different evolutionary stages of star formation: from molecular gas without star formation traced by H$\alpha$ to star-forming gas, and to HII regions. The large sample, drawn from the Physics at High Angular resolution in Nearby GalaxieS ALMA and narrowband H$\alpha$ (PHANGS-ALMA and PHANGS-H$\alpha$) surveys, spans a wide range of stellar mass and morphological types, allowing us to investigate the dependencies of the gas-star formation cycle on global galaxy properties. At a resolution of 150 pc, the incidence of regions in different stages shows a dependence on stellar mass and Hubble type of galaxies over the radial range probed. Massive and/or earlier-type galaxies exhibit a significant reservoir of molecular gas without star formation traced by H$\alpha$, while lower-mass galaxies harbor substantial HII regions that may have dispersed their birth clouds or formed from low-mass, more isolated clouds. Galactic structures add a further layer of complexity to relative distribution of CO and H$\alpha$ emission. Trends between galaxy properties and distributions of gas traced by CO and H$\alpha$ are visible only when the observed spatial scale is $\ll$ 500 pc, reflecting the critical resolution requirement to distinguish stages of star formation process.

Chi-Yan Law, Jonathan C. Tan, Prasanta Gorai, Yichen Zhang, Rubén Fedriani, Daniel Tafoya, Kei Tanaka, Giuliana Cosentino, Yao-Lun Yang, Diego Mardones, Maria Teresa Beltrán, Guido Garay

Submitted to ApJ, comments welcome

We report high-resolution 1.3mm continuum and molecular line observations of the massive protostar G28.20-0.05 with ALMA. The continuum image reveals a ring-like structure with 2,000 au radius, similar to morphology seen in archival 1.3cm VLA observations. Based on its spectral index and associated H30$\alpha$ emission, this structure mainly traces ionized gas. However, there is evidence for $\sim$30 M$\odot$ of dusty gas near the main mm continuum peak on one side of the ring, as well as in adjacent regions within 3,000 au. A virial analysis on scales of $\sim$2,000 au from hot core line emission yields a dynamical mass of $\sim$80M$\odot$. A strong velocity gradient in the H30$\alpha$ emission is evidence for a rotating, ionized disk wind, which drives a larger-scale molecular outflow. An infrared SED analysis indicates a current protostellar mass of m$_{star}\sim$24 M$\odot$ forming from a core with initial mass $M_c\sim400\:M_\odot$ in a clump with mass surface density of $\Sigma_{\rm cl}\sim 3\:{\rm g\:cm}^{-2}$. Thus the SED and other properties of the system can be understood in the context of core accretion models. Structure-finding analysis on the larger-scale continuum image indicates G28.20-0.05 is forming in a relatively isolated environment, with no other concentrated sources, i.e., protostellar cores, above $\sim$ 1 M$\odot$ found from $\sim$0.1 to 0.4 pc around the source. This implies that a massive star is able to form in relative isolation and the dearth of other protostellar companions within the $\sim$1 pc environs is a strong constraint on massive star formation theories that predict the presence of a surrounding protocluster.

Alessandro Paduano, Arash Bahramian, James C. A. Miller-Jones, Adela Kawka, Fabian Göttgens, Jay Strader, Laura Chomiuk, Sebastian Kamann, Stefan Dreizler, Craig O. Heinke, Tim-Oliver Husser, Thomas J. Maccarone, Evangelia Tremou, Yue Zhao

16 pages, 4 figures, 7 tables, accepted for publication in MNRAS

The Galactic globular cluster NGC 3201 is the first Galactic globular cluster observed to host dynamically-confirmed stellar-mass black holes, containing two confirmed and one candidate black hole. This result indicates that globular clusters can retain black holes, which has important implications for globular cluster evolution. NGC 3201 has been observed as part of the MAVERIC survey of Galactic globular clusters. We use these data to confirm that there is no radio or X-ray detection of the three black holes, and present the first radio and X-ray limits on these sources. These limits indicate that any accretion present is at an extremely low rate and may be extremely inefficient. In particular, for the system ACS ID #21859, by assuming the system is tidally locked and any accretion is through the capture of the companion's winds, we constrain the radiative efficiency of any accretion to $\leq1.5\times10^{-5}$. We also combine the radio and X-ray source catalogues from the MAVERIC survey with the existing MUSE spectroscopic surveys and the HUGS catalogue of NGC 3201 to provide a catalogue of 42 multiwavelength sources in this cluster. We identify a new red straggler source with X-ray emission, and investigate the multiwavelength properties of the sub-subgiant population in the cluster.

Beibei Guo, Xufen Wu, Guangwen Chen

12 pages, 11 figures, accepted for publication in ApJ

In the local universe, disc galaxies are generally well evolved and Toomre stable. Their collisions with satellite galaxies naturally produce ring structures, which has been observed and extensively studied. In contrast, at high redshifts, disc galaxies are still developing and clumpy. These young galaxies interact with each other more frequently. However, the products of their collisions remain elusive. Here we systematically study the minor collisions between a clumpy galaxy and a satellite on orbits with different initial conditions, and find a new structure that is different from the local collisional ring galaxies. The clumpness of the target galaxy is fine-tuned by the values of Toomre parameter, $Q$. Interestingly, a thick and knotty ring structure is formed without any sign of a central nucleus in the target galaxy. Our results provide a promising explanation of the empty ring galaxy recently observed in R5519 at redshift $z=2.19$. Moreover, we show that the clumpy state of the collided galaxy exists for a much longer timescale, compared to isolated self-evolved clumpy galaxies that have been widely investigated.

Peter J. Rodenkirch, Cornelis P. Dullemond

Accepted for publication in Astronomy & Astrophysics

Magnetically and thermally driven disk winds have gained popularity in the light of the current paradigm of low viscosities in protoplanetary disks that nevertheless present large accretion rates even in the presence of inner cavities. The possibility of dust entrainment in these winds may explain recent scattered light observations and constitutes a way of dust transport towards outer regions of the disk. We aim to study the dust dynamics in these winds and explore the differences between photoevaporation and magnetically driven disk winds in this regard. We quantify maximum entrainable grain sizes, the flow angle, and the general detectability. We used the FARGO3D code to perform global, 2.5D axisymmetric, nonideal MHD simulations including ohmic and ambipolar diffusion. Dust was treated as a pressureless fluid. Synthetic observations were created with the radiative transfer code RADMC-3D. We find a significant difference in the dust entrainment efficiency of warm, ionized winds such as photoevaporation and magnetic winds including X-ray and extreme ultraviolet (XEUV) heating compared to cold magnetic winds. The maximum entrainable grain size varies from $3\,\mu\mathrm{m}$ to $6\,\mu\mathrm{m}$ for ionized winds to $1\,\mu\mathrm{m}$ for cold magnetic winds. Dust grains in cold magnetic winds tend to flow along a shallower angle compared to the warm winds. With increasing distance to the central star, the dust entrainment efficiency decreases. Larger values of the turbulent viscosity increase the maximum grain size radius of possible dust entrainment. Our simulations indicate that diminishing dust content in the outer regions of the wind can be mainly attributed to the dust settling in the disk. In the synthetic images, the dusty wind appears as a faint, conical emission region which is brighter for a cold magnetic wind.

Jiri Stepan, Tanausu del Pino Aleman, Javier Trujillo Bueno

18 pages, 13 figures, accepted for publication in A&A

Inversion of spectropolarimetric observations of the solar upper atmosphere is one of the most challenging goals in solar physics. If we account for all relevant ingredients of the spectral line formation process such as three-dimensional (3D) radiative transfer out of local thermodynamic equilibrium (NLTE), the task becomes extremely computationally expensive. Instead of generalizing 1D methods to 3D, we develop a new approach to the inverse problem. In our meshfree method we do not consider the requirement of 3D NLTE consistency as an obstacle, but as a natural regularization with respect to the traditional pixel-by-pixel methods. This leads to more robust and less ambiguous solutions. We solve the 3D NLTE inverse problem as an unconstrained global minimization problem avoiding repetitive evaluations of the $\Lambda$~operator. Apart from 3D NLTE consistency, the method allows to easily include additional conditions of physical consistency such as zero divergence of the magnetic field. Stochastic ingredients make the method less prone to ending up in local minima of the loss function. Our method is capable of solving the inverse problem by orders of magnitude faster than it would be possible using grid-based methods. The method can provide accurate and physically consistent results if sufficient computing time is available, but also approximate solutions in case of very complex plasma structures or limited computing time.

M. Gordino, F. Auchère, J.-C. Vial, K. Bocchialini, D. M. Hassler, T. Bando, R. Ishikawa, R. Kano, K. Kobayashi, N. Narukage, J. Trujillo Bueno, A. Winebarger

Empirical relations between major UV and extreme UV spectral lines are one of the inputs for models of chromospheric and coronal spectral radiances and irradiances. They are also needed for the interpretation of some of the observations of the Solar Orbiter mission. We aim to determine an empirical relation between the intensities of the HI 121.6 nm and HeII 30.4 nm Lyman $\alpha$ lines. Images at 121.6 nm from the Chromospheric Lyman Alpha Spectro Polarimeter (CLASP) and Multiple XUV Imager (MXUVI) sounding rockets were co-registered with simultaneous images at 30.4 nm from the EIT and AIA orbital telescopes in order to derive a spatially resolved relationship between the intensities. We have obtained a relationship between the HI 121.6 nm and HeII 30.4 nm intensities that is valid for a wide range of solar features, intensities, and activity levels. Additional SUMER data have allowed the derivation of another relation between the HI 102.5 nm (Ly$\beta$) and HeII 30.4 nm lines for quiet-Sun regions. We combined these two relationships to obtain a Ly$\alpha$/Ly$\beta$ intensity ratio that is comparable to the few previously published results. The relationship between the HI 121.6~nm and HeII 30.4 nm lines is consistent with the one previously obtained using irradiance data. We have also observed that this relation is stable in time but that its accuracy depends on the spatial resolution of the observations. The derived Ly$\alpha$/Ly$\beta$ intensity ratio is also compatible with previous results.

A. Soulain, A. Sivaramakrishnan, P. Tuthill, D. Thatte, K. Volk, R. Cooper, L. Albert, É. Artigau, N. Cook, R. Doyon, D. Johnstone, D. Lafrenière, A. Martel

18 pages, 10 figures, invited paper conference

In less than a year, the James Webb Space Telescope (JWST) will inherit the mantle of being the world's pre-eminent infrared observatory. JWST will carry with it an Aperture Masking Interferometer (AMI) as one of the supported operational modes of the Near-InfraRed Imager and Slitless Spectrograph (NIRISS) instrument. Aboard such a powerful platform, the AMI mode will deliver the most advanced and scientifically capable interferometer ever launched into space, exceeding anything that has gone before it by orders of magnitude in sensitivity. Here we present key aspects of the design and commissioning of this facility: data simulations ($\texttt{ami_sim}$), the extraction of interferometeric observables using two different approaches ($\texttt{IMPLANEIA}$ and $\texttt{AMICAL}$), an updated view of AMI's expected performance, and our reference star vetting programs.

G. Ottoni, S. Udry, D. Ségransan, G. Buldgen, C. Lovis, P. Eggenberger, C. Pezzotti, V. Adibekyan, M. Marmier, M. Mayor, N.C. Santos, S.G. Sousa, N. Lagarde, C. Charbonnel

20 pages, 18 figures, 7 tables, to be published in Astronomy & Astrophysics

Context. Following the first discovery of a planet orbiting a giant star in 2002, we started the CORALIE radial-velocity search for companions around evolved stars (CASCADES). We present the observations of three stars conducted at the 1.2 m Leonard Euler Swiss telescope at La Silla Observatory, Chile, using the CORALIE spectrograph. Aims. We aim to detect planetary companions to intermediate-mass G- and K- type evolved stars and perform a statistical analysis of this population. We searched for new planetary systems around the stars HD22532 (TIC 200851704), HD64121 (TIC 264770836), and HD69123 (TIC 146264536). Methods. We have followed a volume-limited sample of 641 red giants since 2006 to obtain high-precision radial-velocity measurements. We used the Data & Analysis Center for Exoplanets (DACE) platform to perform a radial-velocity analysis to search for periodic signals in the line profile and activity indices, to distinguish between planetary-induced radial-velocity variations and stellar photospheric jitter, and to search for significant signals in the radial-velocity time series to fit a corresponding Keplerian model. Results. In this paper, we present the survey in detail, and we report on the discovery of the first three planets of the sample around the giant stars HD22532, HD64121, and HD69123.

G. Buldgen, G. Ottoni, C. Pezzotti, A. Lyttle, P. Eggenberger, S. Udry, D. Ségransan, A. Miglio, M. Mayor, C. Lovis, Y. Elsworth, G.R. Davies, W.H. Ball

Accepted for publication in Astronomy and Astrophysics. Part of the CASCADES series of papers

The advent of asteroseismology as the golden path to precisely characterize single stars naturally led to synergies with the field of exoplanetology. Today, the precise determination of stellar masses, radii and ages for exoplanet-host stars is a driving force in the development of dedicated software and techniques to achieve this goal. However, as various approaches exist, it is clear that they all have advantages and inconveniences and that there is a trade-off between accuracy, efficiency, and robustness of the techniques. We aim to compare and discuss various modelling techniques for exoplanet-host red giant stars for which TESS data are available. The results of the seismic modelling are then used to study the dynamical evolution and atmospheric evaporation of the planetary systems. We study, in detail, the robustness, accuracy and precision of various seismic modelling techniques when applied to four exoplanet-host red giants observed by TESS. We discuss the use of global seismic indexes, the use of individual radial frequencies and that of non-radial oscillations. In each case, we discuss the advantages and inconveniences of the modelling technique. We determine precise and accurate masses of exoplanet-host red giant stars orbited by long-period Jupiter-like planets using various modelling techniques. For each target, we also provide a model-independent estimate of the mass from a mean density inversion combined with radii values from Gaia and spectroscopic data. We show that no engulfment or migration is observed for these targets, even if their evolution is extended beyond their estimated seismic ages up the red giant branch.

Edward J. Parkinson, Christian Knigge, James H. Matthews, Knox S. Long, Nick Higginbottom, Stuart A. Sim, Samuel W. Mangham

Accepted to MNRAS. 19 pages, 11 figures

A significant number of tidal disruption events (TDEs) radiate primarily at optical and ultraviolet (UV) wavelengths, with only weak soft X-ray components. One model for this optical excess proposes that thermal X-ray emission from a compact accretion disc is reprocessed to longer wavelengths by an optically thick envelope. Here, we explore this reprocessing scenario in the context of an optically thick accretion disc wind. Using state-of-the-art Monte Carlo radiative transfer and ionization software, we produce synthetic UV and optical spectra for wind and disc-hosting TDEs. Our models are inspired by observations, spanning a realistic range of accretion rates and wind kinematics. We find that such outflows can efficiently reprocess the disc emission and produce the broad Balmer and helium recombination features commonly seen in TDEs and exhibit asymmetric red wings. Moreover, the characteristic colour temperature of the reprocessed spectral energy distribution (SED) is much lower than that of the accretion disc. We show explicitly how changes in black hole mass, accretion rate and wind properties affect the observed broadband SED and line spectrum. In general, slower, denser winds tend to reprocess more radiation and produce stronger Balmer emission. Most of the outflows we consider are too highly ionized to produce UV absorption features, but this is sensitive to the input SED. For example, truncating the inner disc at just 4 $R_{ISCO}$ lowers the wind ionization state sufficiently to produce UV absorption features for sight lines looking into the wind

Lorenzo Ingoglia, Giovanni Covone, Mauro Sereno, Carlo Giocoli, Sandro Bardelli, Fabio Bellagamba, Gianluca Castignani, Samuel Farrens, Hendrik Hildebrandt, Shahab Joudaki, Eric Jullo, Denise Lanzieri, Giorgio F. Lesci, Federico Marulli, Matteo Maturi, Lauro Moscardini, Lorenza Nanni, Emanuela Puddu, Mario Radovich, Mauro Roncarelli, Feliciana Sapio, Carlo Schimd

Galaxy clusters are biased tracers of the underlying matter density field. At very large radii beyond about 10 Mpc/\textit{h}, the shear profile shows evidence of a second-halo term. This is related to the correlated matter distribution around galaxy clusters and proportional to the so-called halo bias. We present an observational analysis of the halo bias-mass relation based on the AMICO galaxy cluster catalog, comprising around 7000 candidates detected in the third release of the KiDS survey. We split the cluster sample into 14 redshift-richness bins and derive the halo bias and the virial mass in each bin by means of a stacked weak lensing analysis. The observed halo bias-mass relation and the theoretical predictions based on the $\Lambda$CDM standard cosmological model show an agreement within $2\sigma$. The mean measurements of bias and mass over the full catalog give $M_{200c} = (4.9 \pm 0.3) \times 10^{13} M_{\odot}/\textit{h}$ and $b_h \sigma_8^2 = 1.2 \pm 0.1$. With the additional prior of a bias-mass relation from numerical simulations, we constrain the normalization of the power spectrum with a fixed matter density $\Omega_m = 0.3$, finding $\sigma_8 = 0.63 \pm 0.10$.

C. Pezzotti, G. Ottoni, G. Buldgen, A. Lyttle, P. Eggenberger, S. Udry, D. Ségransan, M. Mayor, C. Lovis, M. Marmier, A. Miglio, Y. Elsworth, G.R. Davies, W.H. Ball

Context. Increasing the number of detected exoplanets is far from anecdotal, especially for long-period planets that require a long duration of observation. More detections imply a better understanding of the statistical properties of exoplanet populations, and detailed modelling of their host stars also enables thorough discussions of star-planet interactions and orbital evolution of planetary systems. Aims. In the context of the discovery of a new planetary system, we aim to perform a complete study of HD 29399 and its companion by means of radial-velocity measurements, seismic characterisation of the host-star, and modelling of the orbital evolution of the system. Methods. High-resolution spectra of HD 29399 were acquired with the CORALIE spectrograph mounted on the 1.2- m Swiss telescope located at La Silla Observatory (Chile) as part of the CASCADES survey. We used the moments of the cross-correlation function profile as well as the photometric variability of the star as diagnostics to distinguish between stellar and planetary-induced signals. To model the host star we combined forward modelling with global and local minimisation approaches and inversion techniques. We also studied the orbital history of the system under the effects of both dynamical and equilibrium tides. Results. We present the detection of a long-period giant planet. Combining these measurements with photometric observations by TESS, we are able to thoroughly model the host star and study the orbital evolution of the system. We derive stellar and planetary masses of $1.17 \pm 0.10~ M_{\odot}$ and $1.59 \pm 0.08 ~M_{Jup}$, respectively, and an age for the system of 6.2 Gyr. We show that neither dynamical nor equilibrium tides have been able to affect the orbital evolution of the planet. Moreover, no engulfment is predicted for the future evolution of the system.

Yifei Ge (Peking University), Ke Wang (Peking University)

Accepted for publication in ApJS. 54 pages, 17 figures, 4 tables. See this http URL for figures with full resolution. Code avaliable at this https URL

The interstellar medium has a highly filamentary and hierarchical structure, which may play a significant role in star formation. A systematical study on the large-scale filaments towards their physical parameters, distribution, structures and kinematics will inform us of what kind of filaments have potential to form stars, how the material feed protostars through filaments, and the connection between star formation and Galactic spiral arms. Unlike the traditional "by eyes" searches, we use a customized minimum spanning tree algorithm to identify filaments by linking Galactic clumps from the APEX Telescope Large Area Survey of the Galaxy catalogue. In the inner Galactic plane ($|l| < 60^\circ$), we identify 163 large-scale filaments with physical properties derived, including dense gas mass fraction, and compare them with an updated spiral arm model in position-position-velocity space. Dense gas mass fraction is found not to differ significantly in various Galactic position, neither does it in different spiral arms. We also find that most filaments are inter-arm filaments after adding a distance constraint, and filaments in arm differ a little with those not in. One surprising result is that clumps on and off filaments have no significant distinction in their mass at the same size.

Dynamical features of a massive disk of distant trans-Neptunian objects are considered in the model of the formation of small bodies in the Hill region of a giant gas-dust clump that arose as a result of gravitational instability and fragmentation of the protoplanetary disk. The dynamical evolution of orbits of small bodies under the action of gravitational perturbations from the outer planets and self-gravity of the disk has been studied for a time interval of the order of a billion years. It is shown that the secular effects of the gravitational action of a massive disk of small bodies lead to an increase in the eccentricities of the orbits of individual objects. The result of this dynamical behavior is the creation of a flux of small bodies coming close to the orbit of Neptune. The change in the number of objects surviving in the observable region of distant trans-Neptunian objects (the region of orbits with perihelion distances of 40 < q < 80 AU and semimajor axes 150 < a < 1000 AU), over time depends on the initial mass of the disk. For disks with masses exceeding several Earth masses, there is a tendency to a decrease in the number of distant trans-Neptunian objects surviving in the observable region after evolution for a time interval of the order of the age of the Solar system, with an increase in the initial mass. On the other hand, for most objects, orbital eccentricities decrease under the influence of the self-gravity of the disk. Therefore, the main part of the disk remains in the region of heliocentric distances exceeding 100 AU.

Compact stellar systems such as Ultra-compact dwarfs (UCDs) and Globular Clusters (GCs) around galaxies are known to be the tracers of the merger events that have been forming these galaxies. Therefore, identifying such systems allows to study galaxies mass assembly, formation and evolution. However, in the lack of spectroscopic information detecting UCDs/GCs using imaging data is very uncertain. Here, we aim to train a machine learning model to separate these objects from the foreground stars and background galaxies using the multi-wavelength imaging data of the Fornax galaxy cluster in 6 filters, namely u, g, r, i, J and Ks. The classes of objects are highly imbalanced which is problematic for many automatic classification techniques. Hence, we employ Synthetic Minority Over-sampling to handle the imbalance of the training data. Then, we compare two classifiers, namely Localized Generalized Matrix Learning Vector Quantization (LGMLVQ) and Random Forest (RF). Both methods are able to identify UCDs/GCs with a precision and a recall of >93 percent and provide relevances that reflect the importance of each feature dimension %(colors and angular sizes) for the classification. Both methods detect angular sizes as important markers for this classification problem. While it is astronomical expectation that color indices of u-i and i-Ks are the most important colors, our analysis shows that colors such as g-r are more informative, potentially because of higher signal-to-noise ratio. Besides the excellent performance the LGMLVQ method allows further interpretability by providing the feature importance for each individual class, class-wise representative samples and the possibility for non-linear visualization of the data as demonstrated in this contribution. We conclude that employing machine learning techniques to identify UCDs/GCs can lead to promising results.

Marius Pfeifer, Jessica Agarwal, Matthias Schröter

16 pages, 20 figures

Context. During the post-perihelion phase of the European Space Agency's Rosetta mission to comet 67P, the Optical, Spectroscopic, and Infrared Remote Imaging System on board the spacecraft took numerous image sequences of the near-nucleus coma, with many showing the motion of individual pieces of debris ejected from active surface areas into space. Aims. We aim to track the motion of individual particles in these image sequences and derive their projected velocities and accelerations. This should help us to constrain their point of origin on the surface, understand the forces that influence their dynamics in the inner coma, and predict whether they will fall back to the surface or escape to interplanetary space. Methods. We have developed an algorithm that tracks the motion of particles appearing as point sources in image sequences. Our algorithm employs a point source detection software to locate the particles and then exploits the image sequences' pair-nature to reconstruct the particle tracks and derive the projected velocities and accelerations. We also constrained the particle size from their brightness. Results. Our algorithm identified 2268 tracks in a sample image sequence. Manual inspection not only found that 1187 (~52%) of them are likely genuine, but in combination with runs on simulated data it also revealed a simple criterion related to the completeness of a track to single out a large subset of the genuine tracks without the need for manual intervention. A tentative analysis of a small (n = 89) group of particles exemplifies how our data can be used, and provides first results on the particles' velocity, acceleration, and radius distributions, which agree with previous work.

Hideki Tanimura, Nabila Aghanim, Victor Bonjean, Saleem Zaroubi

8 pages, 6 figures, submitted to A&A. arXiv admin note: text overlap with arXiv:2007.02952

By stacking the latest 217~GHz Planck PR4 map at the positions of 30,431 galaxy clusters from the Wen-Han-Liu (WHL) catalog, we report the detection of the kinetic Sunyaev-Zel'dovich (kSZ) effect in galaxy clusters with a 4.9 sigma significance level and beyond 3xR500. The line of sight velocities of galaxy clusters were estimated with a machine-learning approach, in which the relation between the galaxy distribution around a cluster and its line-of-sight velocity was trained through a convolutional neural network using the simulated galaxies and galaxy clusters in the Magneticum cosmological hydrodynamic simulations. The trained model was applied to the large-scale distribution of the Sloan Digital Sky Survey galaxies to derive the line-of-sight velocities of the WHL galaxy clusters. Assuming a standard beta-model for the intracluster medium, we found that the gas fraction within R500 is fgas,500 = 0.09+-0.02 for the clusters with the mass of M500~1.0x10^14 Msun.

Giovanni Lapenta, David Schriver, Raymond J. Walker, Jean Berchem, Nicole F. Echterling, Mostafa El Alaoui, Pavel Travnicek

The magnetosphere of Mercury is studied using an implicit full particle in cell simulation (PIC). We use a hybrid simulation where ions are full particles and electrons are considered as a fluid to start a full PIC simulation where electrons are also particles and follow their distribution function. This approach allows us to estimate the changes introduced by the electron kinetic physics. We find that the overall macroscopic state of the magnetosphere of Mercury is little affected but several physical processes are significantly modified in the full PIC simulation: the foreshock region is more active with more intense shock reformation, the Kelvin-Helmholtz rippling effects on the nightside magnetopause are sharper, and the magnetotail current sheet becomes thinner than those predicted by the hybrid simulation. The greatest effect of the electron physics, comes from the processes of particle energization. Both species, not just the electrons, are found to gain more energy when kinetic electron processes are taken into account. The region with the most energetic plasma is found on the dusk side of the tail where magnetic flux ropes are formed due to reconnection. We find that the ion and electron energization is associated with the regions of reconnection and the development of kinetic instabilities caused by counter-streaming electron populations. The resulting electron distributions are highly non Maxwellian, a process that neither MHD nor hybrid models can describe.

Leander Thiele, Digvijay Wadekar, J. Colin Hill, Nicholas Battaglia, Jens Chluba, Francisco Villaescusa-Navarro, Lars Hernquist, Mark Vogelsberger, Daniel Anglés-Alcázar, Federico Marinacci

10+5 pages, 7 figures

High-significance measurements of the monopole thermal Sunyaev-Zel'dovich CMB spectral distortions have the potential to tightly constrain poorly understood baryonic feedback processes. The sky-averaged Compton-y distortion and its relativistic correction are measures of the total thermal energy in electrons in the observable universe and their mean temperature. We use the CAMELS suite of hydrodynamic simulations to explore possible constraints on parameters describing the subgrid implementation of feedback from active galactic nuclei and supernovae, assuming a PIXIE-like measurement. The small 25 Mpc/h CAMELS boxes present challenges due to the significant cosmic variance. We utilize machine learning to construct interpolators through the noisy simulation data. Using the halo model, we translate the simulation halo mass functions into correction factors to reduce cosmic variance where required. Our results depend on the subgrid model. In the case of IllustrisTNG, we find that the best-determined parameter combination can be measured to ~2% and corresponds to a product of AGN and SN feedback. In the case of SIMBA, the tightest constraint is ~0.2% on a ratio between AGN and SN feedback. A second orthogonal parameter combination can be measured to ~8%. Our results demonstrate the significant constraining power a measurement of the late-time spectral distortion monopoles would have for baryonic feedback models.

A. Deconto-Machado, R. A. Riffel, G. S. Ilha, S. B. Rembold, T. Storchi-Bergmann, R. Riffel, J. S. Schimoia, D. P. Schneider, D. Bizyaev, S. Feng, D. Wylezalek, L. N. da Costa, J. C. do Nascimento, M. A. G. Maia

13 pages, 11 figures

We analyse the kinematics of 170 AGN host galaxies as compared to those of a matched control sample of non-active galaxies from the MaNGA survey in order to characterise and estimate the extents of the Narrow Line Region (NLR) and of the kinematically disturbed region (KDR) by the AGN. We define the observed NLR radius as the farthest distance from the nucleus within which both [Oiii]/H\beta and [Nii]/H\alpha ratios fall in the AGN region of the BPT diagram and the H\alpha equivalent width is required to be larger than 3.0\AA. The extent of the KDR is defined as the distance from the nucleus within which the AGN hosts galaxies shows a more disturbed gas kinematics than the control galaxies. The kinematics derived from the [Oiii] line profiles reveal that, on average, the most luminous AGN (L[Oiii] > 3.8 * 10^40 erg s^-1) possess higher residual difference between the gaseous and stellar velocities and velocities dispersion than their control galaxies. Spatially resolved NLR's and KDR's were found in 55 and 46 AGN host galaxies, with corrected radii 0.2 < r_KDR,c < 2.3 kpc and 0.4 < r_NLR,c < 10.1 kpc, with a relation between the two given by log r_KDR,c = (0.53\pm0.12) log r_NLR,c + (1.07\pm0.22), respectively. The extension of the KDR corresponds to about 30 per cent of that of the NLR. Assuming that the KDR is due to an AGN outflow, we have estimated ionised gas mass outflow rates that range between 10^-5 and \approx 1 Myr^-1, and kinetic powers that range from 10^34 to 10^40 erg s^-1. Comparing the power of the AGN ionised outflows with the AGN luminosities, they are always below the 0.05 L_AGN model threshold for having an important feedback effect on their respective host galaxies. The mass outflow rates (and power) of our AGN sample correlate with their luminosities, populating the lowest AGN luminosity range of the correlations previously found for more powerful sources.

M.V. del Valle, A. Araudo, F. Suzuki-Vidal

10 pages, 11 figures. Accepted for publication in A&A

The termination regions of non-relativistic jets in protostars and supersonic outflows in classical novae are nonthermal emitters. Given the high densities in these systems, radiative shocks are expected to form. However, in the presence of high velocities, the formation of adiabatic shocks is also possible. A case of interest is when the two types of shocks occur simultaneously. These dense jets/outflows are excellent candidates for laboratory experiments as demonstrated by MHD scaling. We aim at studying the combination of adiabatic and radiative shocks in these systems. We focus on determining the conditions under which this combination is feasible together with its physical implications. We perform an analytical study of the shocks in both types of sources for a set of parameters. The hydrodynamical evolution of a jet colliding with an ambient medium is studied with 2D numerical simulations confirming our initial theoretical estimates. We show that for a wide set of parameters the combination of an adiabatic and a radiative shock is possible at the working surface of the termination region in jets from young stars and novae outflows. We find that instabilities are developed at the contact discontinuity, mixing the shocked materials. Also, we explore the MHD parameter scaling required for studying protostellar jets and novae outflows using laboratory experiments on laser facilities. The coexistence of an adiabatic and a radiative shock is expected at the termination region of protostellar jets and novae outflows. This scenario is very promising for particle acceleration and gamma-ray emission. The parameters for scaled laboratory experiments are very much in line with plasma conditions achievable in currently operating high-power laser facilities. This opens the door to new means for studying novae outflows never considered before.

Samuel Gill, Solene Ulmer-Moll, Peter J. Wheatley, Daniel Bayliss, Matthew R. Burleigh, Jack S. Acton, Sarah L. Casewell, Christopher A. Watson, Monika Lendl, Hannah L. Worters, Ramotholo R. Sefako, David R. Anderson, Douglas R. Alves, François Bouchy, Edward M. Bryant, Philipp Eigmüller, Edward Gillen, Michael R. Goad, Nolan Grieves, Maximilian N. Günther, Beth A. Henderson, James S. Jenkins, Lokesh Mishra, Maximiliano Moyano, Hugh P. Osborn, Rosanna H. Tilbrook, Stéphane Udry, Jose I. Vines, Richard G. West

9 pages, 6 figures, 4 tables. Submitted to MNRAS

We are using precise radial velocities from CORALIE together with precision photometry from the Next Generation Transit Survey (NGTS) to follow up stars with single-transit events detected with the Transiting Exoplanet Survey Satellite (TESS). As part of this survey we identified a single transit on the star TIC-320687387, a bright (T=11.6) G-dwarf observed by TESS in Sector 13 and 27. From subsequent monitoring of TIC-320687387 with CORALIE, NGTS, and Lesedi we determined that the companion, TIC-320687387 B,is a very low-mass star with a mass of $96.2 \pm _{2.0}^{1.9} M_J$ and radius of $1.14 \pm _{0.02}^{0.02} R_J$ placing it close to the hydrogen burning limit ($\sim 80 M_J$). TIC-320687387 B has a wide and eccentric orbit, with a period of 29.77381 days and an eccentricity of $0.366 \pm 0.003$. Eclipsing systems such as TIC-320687387 AB allow us to test stellar evolution models for low-mass stars, which in turn are needed to calculate accurate masses and radii for exoplanets orbiting single low-mass stars. The wide orbit of TIC-320687387 B makes it particularly valuable as its evolution can be assumed to be free from perturbations caused by tidal interactions with its G-type host star.

Gerald Handler, Rahul Jayaraman, Donald W. Kurtz, Jim Fuller, Saul A. Rappaport

6 pages, 6 figures, Proceedings of the 40th meeting of the Polish Astronomical Society, submitted

The tidally tilted pulsators are a new type of oscillating star in close binary systems that have their pulsation axis in the orbital plane because of the tidal distortion caused by their companion. We describe this group of stars on the basis of the first three representatives discovered and illustrate the basic methods used for their analysis. Their value for astrophysical study is rooted in the combination of the strengths of binary star and asteroseismic analyses; pulsational mode identifications can be achieved because the oscillations are visible over nearly 360 degrees of aspect throughout the orbital cycle. An illustrative case of a particularly interesting system is presented.

We study the secular evolution of a particle in deep mean motion resonance (MMR) with a planet in the planar elliptic restricted three body problem. We do not consider any restriction neither in the planet's eccentricity $e_p$ nor in the particle's eccentricity $e$. The methodology used is based on a semi-analytical model that consists on calculating the averaged resonant disturbing function numerically, assuming for this that in the resonant scale of time all the orbital elements of the particle are constant. In order to obtain the secular evolution inside the MMR, we make use of the adiabatic invariance principle, assuming a zero-amplitude resonant libration. We construct two-dimensional surfaces (called $\mathcal{H}$ surfaces) in the three-dimensional space $(\sigma, e, \varpi)$ that allow us to predict the secular evolution of these three variables. The 2:1 MMR is used as example to show some results. We found four apsidal corotation resonance (ACR) families, two symmetric and two asymmetric. One of the symmetric families exists for almost any $e_p$ value. The other one for $e_p>0.3$ and the asymmetric ones for $e_p>0.44$. We corroborate the secular variations in $e$ and $\varpi$ predicted by the model through numerical integrations even when the initial conditions are displaced from those ACR. Some peculiar examples are presented for the 3:1 and 3:2 MMR showing large excursions in eccentricity. As an application, the Planet 9 is investigated as a possible responsible of high eccentric distant TNOs.

S. D. Raducan, M. Jutzi, T. M. Davison, M. E. DeCoster, D. M. Graninger, J. M. Owen, A. M. Stickle, G. S. Collins

The DART spacecraft will impact Didymos's secondary, Dimorphos, at the end of 2022 and cause a change in the orbital period of the secondary. For simplicity, most previous numerical simulations of the impact used a spherical projectile geometry to model the DART spacecraft. To investigate the effects of alternative, simple projectile geometries on the DART impact outcome we used the iSALE shock physics code in two and thee-dimensions to model vertical impacts of projectiles with a mass and speed equivalent to the nominal DART impact, into porous basalt targets. We found that the simple projectile geometries investigated here have minimal effects on the crater morphology and momentum enhancement. Projectile geometries modelled in two-dimensions that have similar surface areas at the point of impact, affect the crater radius and the crater volume by less than 5%. In the case of a more extreme projectile geometry (i.e., a rod, modelled in three-dimensions), the crater was elliptical and 50% shallower compared to the crater produced by a spherical projectile of the same momentum. The momentum enhancement factor in these test cases, commonly referred to as beta, was within 7% for the 2D simulations and within 10% for the 3D simulations, of the value obtained for a uniform spherical projectile. The most prominent effects of projectile geometry are seen in the ejection velocity as a function of launch position and ejection angle of the fast ejecta that resides in the so-called `coupling zone'. These results will inform the LICIACube ejecta cone analysis.

Mireille Louys, Laurent Michel, François Bonnarel, Joann Vetter

4 pages , 2 figures , ADASS XXX1 Conference Proceedings held on 24-28 October 2021

With the success and widespread of the IVOA Table Access Protocol (1) for discovering and querying tabular data in astronomy, more than one hundred of TAP services exposing altogether 22 thousands of tables are accessible from the IVOA Registries at the time of writing. Currently the TAP protocol presents table data and metadata via a {TAP\_SCHEMA} describing the served tables with their columns and possible joins between them. We explore here how to add an information layer, so that values within table columns can be gathered and used to populate instances of objects defined in a selected IVOA data model like Photometry, Coords, Measure, Transform or the proposed MANGO container model. This information layer is provided through annotation tags which tell how the columns' values can be interpreted as attributes of instances of that model. Then when a TAP query is processed, our server add-on interprets the ADQL query string and produces on-the-fly, when possible, the TAP response as an annotated VOTable document. The FIELD elements in the table response are mapped to corresponding model elements templated for this service. This has been prototyped in Java, using the VOLLT package library and a template annotation document representing elements from the MANGO data model. This has been exercised on examples based on Vizier and Chandra catalogs.

Carlos G. Escudero, Arianna Cortesi, Favio R. Faifer, Leandro A. Sesto, Analía V. Smith Castelli, Evelyn J. Johnston, Victoria Reynaldi, Ana L. Chies-Santos, Ricardo Salinas, Karín Menéndez-Delmestre, Thiago S. Gonçalves, Marco Grossi, Claudia Mendes de Oliveira

20 pages, 24 figures. Accepted to be published in MNRAS

NGC 4382 is a merger-remnant galaxy that has been classified as morphological type E2, S0, and even Sa. In this work, we performed a photometric and spectroscopic analysis of the globular cluster (GC) system of this peculiar galaxy in order to provide additional information about its history. We used a combination of photometric data in different filters, and multi-object and long-slit spectroscopic data obtained using the Gemini/GMOS instrument. The photometric analysis of the GC system, using the Gaussian Mixture Model algorithm in the colour plane, reveals a complex colour distribution within $R_\mathrm{gal}<5$ arcmin (26.1 kpc), showing four different groups: the typical blue and red subpopulations, a group with intermediate colours, and the fourth group towards even redder colours. From the spectroscopic analysis of 47 GCs, confirmed members of NGC\,4382 based on radial velocities, we verified 3 of the 4 photometric groups from the analysis of their stellar populations using the ULySS code. NGC 4382 presents the classic blue ($10.4\pm2.8$ Gyr, $\mathrm{[Fe/H]}=-1.48\pm0.18$ dex) and red ($12.1\pm2.3$ Gyr, $\mathrm{[Fe/H]}=-0.64\pm0.26$ dex) GCs formed earlier in the lifetime of the galaxy, and a third group of young GCs ($2.2\pm0.9$ Gyr; $\mathrm{[Fe/H]}=-0.05\pm0.28$ dex). Finally, analysis of long-slit data of the galaxy reveals a luminosity-weighted mean age for the stellar population of $\sim$2.7 Gyr, and an increasing metallicity from [Fe/H]=$-0.1$ to $+0.2$ dex in $R_\mathrm{gal}<10$ arcsec (0.87 kpc). These values, and other morphological signatures in the galaxy, are in good agreement with the younger group of GCs, indicating a common origin as a result of a recent merger.

Edward Nathan, Adam Ingram, Jeroen Homan, Daniela Huppenkothen, Phil Uttley, Michiel van der Klis, Sara Motta, Diego Altamirano, Matthew Middleton

26 pages, 16 figures, accepted to MNRAS

Quasi-periodic oscillations (QPOs) are often present in the X-ray flux from accreting stellar-mass black holes (BHs). If they are due to relativistic (Lense-Thirring) precession of an inner accretion flow which is misaligned with the disc, the iron emission line caused by irradiation of the disc by the inner flow will rock systematically between red and blue shifted during each QPO cycle. Here we conduct phase-resolved spectroscopy of a $\sim2.2$ Hz type-C QPO from the BH X-ray binary GRS 1915+105, observed simultaneously with NICER and NuSTAR. We apply a tomographic model in order to constrain the QPO phase-dependent illumination profile of the disc. We detect the predicted QPO phase-dependent shifts of the iron line centroid energy, with our best fit featuring an asymmetric illumination profile ($>2{\sigma}$ confidence). The observed line energy shifts can alternatively be explained by the spiral density waves of the accretion-ejection instability model. However we additionally measure a significant ($>3{\sigma}$) modulation in reflection fraction, strongly favouring a geometric QPO origin. We infer that the disc is misaligned with previously observed jet ejections, which is consistent with the model of a truncated disc with an inner precessing hot flow. However our inferred disc inner radius is small ($r_\text{in}{\sim} 1.4 GM/c^2$). For this disc inner radius, Lense-Thirring precession cannot reproduce the observed QPO frequency. In fact, this disc inner radius is incompatible with the predictions of all well-studied QPO models in the literature.

Natalia Osorio-Clavijo, Omaira González-Martín, Sebastián F. Sánchez-Sánchez, Donaji Esparza-Arredondo, Josefa Masegosa, César Victoria-Ceballos, Lorena Hernández-García, Yaherlyn Díaz

17 pages, 12 figures, 5 tables. Accepted for publication in MNRAS

According to theory, the torus of active galactic nuclei (AGN) is sustained from a wind coming off the accretion disk, and for low efficient AGN, it has been proposed that such structure disappears. However, the exact conditions for its disappearance remain unclear. This can be studied throughout the reflection component at X-rays, which is associated with distant and neutral material at the inner walls of the torus in obscured AGN. We select a sample of 81 AGNs observed with NuSTAR with a distance limit of D< 200\,Mpc and Eddington rate $\rm{\lambda_{Edd} \equiv L_{bol}/L_{Edd}<10^{-3}}$. We fit the 3-70\,keV spectra using a model accounting for a partial-covering absorber plus a reflection component from neutral material. We find that the existence of the reflection component spans in a wide range of black-hole mass and bolometric luminosities, with only $\sim$13\% of our sample (11 sources) lacking of any reflection signatures. These sources fall in the region in which the torus may be lacking in the L-MBH diagram. For the sources with a detected reflection component, we find that the vast majority of them are highly obscured ($\rm{\log \ N_H > 23}$), with $\rm{\sim 20\%}$ being Compton-thick. We also find an increase on the number of unobscured sources and a tentative increase on the ratio between $\rm{FeK\alpha}$ emission line and Compton-hump luminosities toward $\rm{\lambda_{Edd}=10^{-5}}$, suggesting that the contribution of the $\rm{FeK\alpha}$ line changes with Eddington ratio.

The neutral, flavor singlet scalar uuddss bound state -- the sexaquark, S -- may have a low enough mass to be stable or extremely long-lived. Here we review mass estimates and production expectations and show that laboratory experiments to date do not rule out such a long-lived state. An S with mass below 2054 MeV is either absolutely stable or has a lifetime greater than the age of the Universe. Detection of a stable S in accelerator experiments is very challenging. An examination of the experimental literature shows that such an effectively stable state would have escaped detection. The strongest laboratory constraint on a long-lived S comes from the lower-bound on its formation time in a doubly-strange hypernucleus; this constraint is, however, not stringent enough to exclude a stable S. We develop strategies to discover it. A stable S would be an attractive dark matter candidate. Relevant astrophysical and cosmological observations, which show that sexaquark dark matter is consistent with all current knowledge, are briefly reviewed.

Kadhom F. Fadhel, A. A. Al-Rubaiee

8 pages, 4 figures, Proceedings of the Mustansiriyah International Conference on Applied Physics (MICAP-2021)

In this study, the necessity of the simulation study for exploring the interactions of ultra-high energy particles cosmic rays was examined. Different hadronic interaction models such as (SIBYLL, QGSJET, and EPOS) were simulated by using air showers simulation AIRES system (version 19.04.00). Also, the charged particle density of Extensive Air Showers (EAS) was calculated by estimating the lateral distribution function (LDF). Moreover, the LDF simulation of the two primary particles (proton and iron nuclei) was performed, taking into account their primary energies effect and the zenith angle of charged particles that produced in the EAS, within the energy range (10^17-10^19) eV. At extremely high energies (10^17, 10^18, and 10^19) eV, new parameters as a function of the primary energy were obtained by fitting the lateral distribution curves of EAS using Sigmoidal function (Logistic model). Comparison of the results showed a good agreement between the values obtained from the parameterized LDF using Sigmoidal function with experimental results by AGASA EAS observatory for the primaries proton as well iron nuclei, with the production of (electron positron) pair and the charged muons secondary particles at high energy about 10^19 eV and (theta = 0 degree).

The paper presents the results from a multi-year effort to develop and validate image processing methods for selecting the best physical models based on solar image observations. The approach consists of selecting the physical models based on their agreement with coronal holes extracted from the images. Ultimately, the goal is to use physical models to predict geomagnetic storms. We decompose the problem into three subproblems: (i) coronal hole segmentation based on physical constraints, (ii) matching clusters of coronal holes between different maps, and (iii) physical map classification. For segmenting coronal holes, we develop a multi-modal method that uses segmentation maps from three different methods to initialize a level-set method that evolves the initial coronal hole segmentation to the magnetic boundary. Then, we introduce a new method based on Linear Programming for matching clusters of coronal holes. The final matching is then performed using Random Forests. The methods were carefully validated using consensus maps derived from multiple readers, manual clustering, manual map classification, and method validation for 50 maps. The proposed multi-modal segmentation method significantly outperformed SegNet, U-net, Henney-Harvey, and FCN by providing accurate boundary detection. Overall, the method gave a 95.5% map classification accuracy.

Lawrence M. Krauss (Origins Project Foundation)

3 pages, submitted to Phys. Rev. Lett

We derive two features of axion cosmology that may have cosmological implications, whether or not axions are dark matter: For the full range of allowed axion masses, the evolution of a cosmic axion background allows large CP violation until temperatures as low as $\sim$ 2 GeV, and once the axion field begins to oscillate, the cosmological axion field's relaxation to its ground state can briefly provide a new departure from thermal equilibrium, via time-varying CP violation. During both of these periods, the Strong CP violating parameter $\bar\theta$ can be as large as O(1).

M. G. Boulay, V. Camillo, N. Canci, S. Choudhary, L. Consiglio, A. Flammini, C. Galbiati, C. Ghiano, A. Gola, S. Horikawa, P. Kachru, I. Kochanek, K. Kondo, G. Korga, M. Kuźniak, A. Mazzi, A. Moharana, G. Nieradka, G. Paternoster, A. Razeto, D. Sablone, T.N. Thorpe, C. Türkoğlu, H. Wang, M. Rescigno, S. Sanfilippo

8 pages, 8 figures

SiPM-based readout is becoming the standard for light detection in particle detectors given their superior resolution and simplicity of use respect to vacuum tube photo-multipliers. However, the presence of detection noises such as the dark rate, the cross-talk and after-pulsing may impact significantly the performance. In this work, we document the development of highly reflective single-phase argon chambers capable of displaying light yields up to 32 photo-electrons per keV, out of which only about 12 are primary photo-electrons generated by the argon scintillation and the rest is accounted for by optical cross-talk. Furthermore, the presence of compound processes yields a generalized Fano factor that reaches a value of 9 at higher over-voltages. As a bonus, we have the complete parametrization of the optical cross-talk for the FBK NUV-HD-Cryo SiPMs at 87 K.

Gravitational lensing of a light source by a black hole leads to appearance of higher-order images produced by photons that loop around the black hole before reaching the observer. Higher-order images were widely investigated numerically and analytically, in particular using so-called strong deflection limit of gravitational deflection. After recent observations of the black hole image, attention has been drawn to higher-order rings, which are lensed images of the accreting matter of the black hole environment and can appear near the boundary of the black hole shadow. In this article, we use strong deflection limit technique to investigate higher-order ring images of luminous accretion disc around a Schwarzschild black hole. We consider thin disk given by the inner and outer radii and an observer located far from the black hole on the axis of symmetry. For this configuration, it becomes possible to find the angular radii, thicknesses, and solid angles of higher-order rings in the form of compact analytical expressions. We show that the size of the rings is mainly determined by the position of the inner boundary of the accretion disk, which makes it possible to use them to distinguish between different accretion models. Possible generalizations of our model to non-symmetric images can help to make the estimation of black hole angular momentum. We also present the analytical estimation of fluxes from higher-order images. Our method makes it easy to investigate $n=2$ and $n=3$ higher-order rings, the possible observation of which in future projects is currently being discussed.